Home
JournalsCollections
For Authors For Reviewers For Editorial Board Members
Article Processing Charges Open Access
Ethics Advertising Policy
Editorial Policy Resource Center
Company Information Contact Us
OPEN ACCESS

Association of GCKR Gene Polymorphisms with the Risk of Nonalcoholic Fatty Liver Disease and Coronary Artery Disease in a Chinese Northern Han Population

  • Hui Gao#,1,
  • Shousheng Liu#,2,4,
  • Zhenzhen Zhao2,4,
  • Xinjuan Yu2,
  • Qun Liu1,
  • Yongning Xin*,1,3,4 and
  • Shiying Xuan*,1,3,4
Journal of Clinical and Translational Hepatology   2019;7(4):297-303

doi: 10.14218/JCTH.2019.00030

Received:

Revised:

Accepted:

Published online:

 Author information

Citation: Gao H, Liu S, Zhao Z, Yu X, Liu Q, Xin Y, et al. Association of GCKR Gene Polymorphisms with the Risk of Nonalcoholic Fatty Liver Disease and Coronary Artery Disease in a Chinese Northern Han Population. J Clin Transl Hepatol. 2019;7(4):297-303. doi: 10.14218/JCTH.2019.00030.

Abstract

Background and Aims: Accumulated studies have evaluated the effects of glucokinase regulatory protein (GCKR) gene polymorphisms on the risk of nonalcoholic fatty liver disease (NAFLD) and coronary artery disease (CAD), but the association of GCKR polymorphisms with the risk of NAFLD and CAD in the Chinese Han population have remained unclear. The aim of this study was to investigate the association between GCKR gene polymorphisms (rs780094 and rs1260326) and the risk of NAFLD and CAD in NAFLD patients in a Chinese Northern Han population.

Methods:GCKR rs780094 and rs1260326 gene polymorphisms were genotyped by polymerase chain reaction sequencing for B-type ultrasonography-proven NAFLD patients with (n = 82) or without (n = 142) CAD, and in healthy controls (n = 152). Serum lipid profiles’ levels were determined using biochemical methods. Statistical analyses were conducted using SPSS 22.0 statistical software.

Results: As the results showed, significant differences in the serum lipid profiles existed between each group. No significant differences were observed in the distributions of genotypes and alleles of GCKR rs780094 and rs1260326 in each group. The GCKR rs780094 T and rs1260326 T allele carriers possessed decreased body mass index value, and serum fasting plasma glucose and TG levels in the overall subjects, respectively. In addition, the GCKR rs780094 T allele carriers possessed decreased serum fasting plasma glucose level in the controls and NAFLD + CAD patients.

Conclusions:GCKR rs780094 and rs1260326 polymorphisms were found to be not associated with the risk of NAFLD nor of CAD in NAFLD patients in this Chinese Northern Han population. GCKR rs780094 T and rs1260326 T alleles could affect the body mass index value and serum fasting plasma glucose and triglyceride levels.

Keywords

Non-alcoholic fatty liver disease, Glucokinase regulatory protein gene, Polymorphism, Coronary artery disease

Introduction

Nonalcoholic fatty liver disease (NAFLD) is the most prevalent chronic liver disease in the world and is regarded as a severe public health concern.1 The overall prevalence of NAFLD is approximately 25% worldwide and is 23%, 31% and 27% in Europe, South America and Asia, respectively.2–4 Many factors contribute to the development of NAFLD, such as aging, hyperlipidemia, insulin resistance, type 2 diabetes, diet, genetics and so on.5–8 Accumulated lines of evidence have suggested that coronary artery disease (CAD) plays an important role in the progression of NAFLD and is associated with the severity of NAFLD.

CAD has become the most important cause of mortality among NAFLD patients.9 Although some defects of liver biopsy exist, such as invasive, sample errors and operator dependence, it remains the gold standard for diagnosis of NAFLD.10,11 In recent years, some noninvasive diagnostic methods have been developed and their effects have been tested in different countries and populations. Genome-wide association studies (GWAS) have identified several important single nucleotide polymorphism (SNP) sites which are tightly associated with the risk of development of NAFLD; these include PNPLA3 rs738409, TM6SF2 rs58542926, and LYPLAL1 rs12137855.12 In consideration of the differences of genetic background of NAFLD patients, it is necessary to identify the relationship of gene polymorphisms with the risk of NAFLD in different countries and races.

Glucokinase (GCK) is a phosphorylating enzyme which can regulate hepatic glucose metabolism and activate hepatic lipogenesis.13 The GCKR gene encodes the glucokinase regulatory protein (GCKRP), which can bind to GCK allosterically and regulate the activity of GCK.14 In 2011, Speliotes et al.15 identified the GCKR rs780094 by GWAS as a significant SNP site associated with the risk of NAFLD. Subsequent studies in different districts reported controversial conclusions of GCKR rs780094 in NAFLD. The GCKR rs780094 T allele was found to be tightly associated with susceptibility to NAFLD in Americans and Malaysians, acting as a risk factor for NAFLD.15,16 However, no significant association between GCKR rs780094 and the risk of NAFLD was observed in the community-based study in Suzhou and Shanghai, and in obese adolescents in Beijing, China.17–20 On the contrary, Lin et al.21 reported that GCKR rs780094 was associated with an increased risk of NAFLD in obese children in Taiwan.22 Therefore, large-scale population-based studies in different districts are needed to investigate the effect of GCKR rs780094 on the development of NAFLD in China.

In addition, GWAS analyses conducted in Finnish, Swedish and Danish populations found a strong linkage disequilibrium of GCKR rs780094 with GCKR rs1260326.23,24 In Japan, Kawaguchi et al.25 found that GCKR rs1260326 was a significant risk factor for the development of NAFLD. In France, Petit et al.26 reported that GCKR rs1260326 influences the liver fat content in patients with type 2 diabetes. Di et al.27 reported that GCKR rs1260326 could increase the liver fat content, although the detailed role of GCKR rs1260326 in NAFLD remains controversial. Besides, Simons et al.28 conducted a meta-analysis to investigate the relationship of GCKR polymorphisms with the risk of CAD, and a tight association of GCKR variants with the risk of CAD was observed.

In consideration of the unclear effects of GCKR rs780094 and rs1260326 on NAFLD in Chinese, and the tight association of CAD and NAFLD, it is of interest to explore the effects of GCKR rs780094 and rs1260326 on the development of NAFLD and the development of CAD in NAFLD patients from this population.

The aim of this study was to explore the relationship between GCKR gene polymorphisms (rs780094 and rs1260326) with the risk of NAFLD, the risk of CAD in NAFLD patients in a Chinese Northern Han population, and the effects of GCKR rs780094 and rs1260326 on the levels of serum lipid profiles.

Methods

Subjects

This case-control study was conducted according to the principles of the Declaration of Helsinki and its appendices.29 This study was approved by the ethics committee of Qingdao Municipal Hospital (Qingdao, China). All the subjects included in this study were of the Chinese Northern Han population, and each signed a written informed consent form.

From June 2018 to December 2018, 225 unrelated adult Chinese Northern Han patients with NAFLD, of both genders, that were diagnosed by B-type ultrasonography were recruited. Among these NAFLD patients, 82 suffered from CAD (41 males, 41 females; mean age of 55.62 ± 7.187 years) and 142 did not have CAD (75 males, 67 females; mean age of 50.14 ± 11.842 years). A total of 152 healthy controls, matched for sex and age (75 males, 77 females; mean age of 51.38 ± 13.276 years) were included. All the patients were recruited from the Departments of Gastroenterology and Cardiology and the Health Examination Center of Qingdao Municipal Hospital.

Standard clinical evaluations of the NAFLD patients were performed according to the criteria of American Association for the Study of Liver Diseases.21 Subjects with the following symptoms were excluded from the NAFLD cohort: 1) excessive alcohol intake (males >210 g/w, females >140 g/w); 2) viral hepatitis, such as hepatitis B or hepatitis C; 3) drug-induced hepatitis; 4) other endocrine and metabolic disorders, such as renal disease. CAD was diagnosed using the findings from percutaneous coronary angiogram conducted by two experienced interventional cardiologists, and defined as the presence of at least 50% stenosis in at least one of the coronary arteries. The healthy controls were also subjected to the same diagnosis routines for NAFLD and CAD, as well as the laboratory and general examinations.

Biochemical analyses

The basic clinical pathological information (gender, age, height, and weight) was obtained by a standard study questionnaire. The body mass index (BMI) was calculated as mass weight (kg)/height2 (m2). All the subjects underwent a 12-h overnight fast, and then the blood samples were collected into ethylene diaminetetraacetic acid-containing tubes. Serum fasting plasma glucose (FPG), total cholesterol (TC), triglycerides (TG), high-density lipoprotein (HDL), low-density lipoprotein (LDL), alanine aminotransferase (ALT), aspartate aminotransferase (AST), glutamyltransferase (r-GT), alkaline phosphatase (ALP), and bilirubin (BIL) were measured using standard clinical laboratory techniques.30

Genomic DNA extraction and genotyping

Genomic DNA was extracted and stored as previously described.31 The primers for PCR amplification containing rs780094 were 5′-ACGTTGGATGAAGTGGTGGGATTACAGGAG-3′ and 5′-ACGTTGGATGAGGGCCCCAGTTTTTTAGAC-3′, rs1260326 5′-ACGTTGGATGACCTGGGTCCCTTTGTCAC-3′ and 5′-ACGTTGGATGAGAAGACCAACCACATCCAG-3′ and were synthesized by Beijing BoMiao Biotech Company (China). The PCR amplification profile was as follows: predenaturation at 94°C for 5 m, followed by 40 cycles of denaturation at 94°C for 20 s, annealing at 56°C for 30 s, extending at 72°C for 60 s, followed by the final extension at 72°C for 3 m. The target amplified fragments were detected by 2% gel electrophoresis, with the predicted fragment length. The genotypes of rs780094 and rs1260326 were detected by direct DNA sequencing using the Applied Biosystems Inc. (USA) Veriti-384 Prism Sequence Detection System, and the raw data were analyzed using MassARRAY TYPER4.0 software (Agena Bioscience, USA). Genotyping was performed in a blinded fashion and the success rates were >95%.

Statistical analysis

Statistical analysis was conducted using the SPSS 22.0 statistical software (SPSS Inc., USA). Student’s t-test or χ2 test were used to analyzed the differences in characteristics of different groups. First, continuous variables were tested by normal distribution and F test; if yes, the baseline characteristics were expressed as the mean ± standard deviation (SD), and if no, the baseline characteristics were expressed as the median ± quartile. Hardy-Weinberg equilibrium was between expected and observed genotype distributions and was estimated using the χ2 test. Genotype and allele frequencies were assessed by counting the DNA sequencing data, and the distributions in NAFLD versus controls were analyzed by Pearson’s χ2 test or Fisher’s exact test, where appropriate. Association between polymorphisms and NAFLD/CAD were analyzed by logistic regression and odds ratios (ORs) with 95% confidence intervals (CIs). A p-value <0.05 was considered statistically significant.

Results

Characteristics of the study participants

The baseline characteristics of the three groups (NAFLD, NAFLD with CAD, and healthy controls) according to the experimental requirements are shown in Table 1. There were no significant differences in age or gender among the three groups (all p > 0.05). The increased serum levels of BMI, TC, TG, LDL, ALT, AST, and r-GT were observed in the NAFLD patients compared to the healthy controls (all p < 0.05). In the NAFLD + CAD group, the BMI value and serum levels of FPG, TC, TG, LDL, ALT, AST, r-GT, and ALP were higher and the serum HDL level was lower than in healthy controls (all p < 0.05). Compared to the NAFLD patients, the patients with NAFLD and CAD possessed higher serum levels of FPG, AST and ALP, and lower BMI value and serum levels of TC, HDL, and LDL (all p < 0.05).

Table 1.

Baseline characteristics of the study participantsa

CharacteristicControlsNAFLDNAFLD + CADpp1p2p3
Age (y)51.38 ± 13.27650.14 ± 11.84255.62 ± 7.1870.156---
Gender (F/M)77/7567/7541/410.174---
BMI (kg/m2)23.45 ± 3.37626.74 ± 2.80525.24 ± 2.606<0.001<0.001<0.0010.003
FPG (mmol/L)4.80 ± 1.4795.00 ± 1.3435.90 ± 2.180<0.0010.010<0.0010.001
TC (mmol/L)4.99 ± 1.1925.48 ± 0.8224.56 ± 1.274<0.0010.010<0.001<0.001
TG (mmol/L)1.45 ± 1.0091.90 ± 1.9032.23 ± 4.396<0.0010.003<0.0011.000
HDL (mmol/L)1.28 ± 0.3871.24 ± 0.1981.05 ± 0.294<0.0011.000<0.001<0.001
LDL (mmol/L)3.10 ± 0.7083.27 ± 0.5552.75 ± 0.972<0.0010.0340.006<0.001
ALT (U/L)19.90 ± 11.15534.40 ± 26.16528.90 ± 25.044<0.001<0.0010.0060.019
AST (U/L)21.93 ± 8.87724.16 ± 8.99444.67 ± 69.6490.067---
r-GT (mmol/L)26.30 ± 24.36739.63 ± 26.48147.89 ± 85.986<0.001<0.0010.0100.019
ALP (mmol/L)75.14 ± 24.66171.62 ± 18.27285.79 ± 21.840<0.0010.852<0.001<0.001
BIL (mmol/L)13.20 ± 5.35813.30 ± 4.95714.44 ± 7.0200.475---

GCKR rs780094 and rs1260326 genotype and allele distribution

The genotype distributions of GCKR rs780094 and rs1260326 were in accordance with the Hardy-Weinberg equilibrium in each group (all p > 0.05) (Table 2). As the results show in Table 3, there were no significant differences in the genotypic and allelic distributions of GCKR rs780094 and rs1260326 between the NAFLD patients and healthy controls, NAFLD + CAD patients and healthy controls, nor NAFLD + CAD patients and NAFLD patients, respectively (all p > 0.05). These results suggested that GCKR rs780094 and rs1260326 polymorphisms did not associate with the risk of NAFLD in this Northern Chinese Han population, and GCKR rs780094 and rs1260326 polymorphisms did not increase the risk of CAD in these NAFLD patients.

Table 2.

Results of the Hardy-Weinberg equilibrium testa

Gene locus groupsCCCTTTχ2p value
rs780094
 Controls2681451.0400.308
 NAFLD2671450.0470.828
 NAFLD + CAD1839250.0140.708
rs1260326
 Controls2382461.8900.169
 NAFLD2473450.3720.542
 NAFLD + CAD1541260.0280.868
Table 3.

Distribution of the GCKR rs780094 and rs1260326 polymorphisms in the study groupsa,b

GenotypeControlsNAFLDNAFLD + CADOR (95%CI)χ2p1OR (95%CI)χ2p2OR(95%CI)χ2p3
rs780094
 CC26 (17.1%)26 (18.3%)18 (22.0%)
 CT + TT126 (82.9%)116 (81.7%)64 (78.0%)0.837 (0.455–1.538)0.0730.7870.734 (0.375–1.473)0.8190.3651.255 (0.640–2.642)0.4370.509
Alleles
 C133 (43.8%)123 (43.3%)75 (45.7%)
 T171 (56.2%)161 (56.7%)89 (54.3%)0.982 (0.709–1.361)0.0120.9141.083 (0.740–1.587)0.1690.6811.103 (0.749–1.624)0.2470.619
rs1260326
 CC23 (15.1%)24 (16.9%)15 (18.3%)
 CT + TT129 (84.9%)118 (83.1%)67 (81.7%)0.877 (0.470–1.636)0.1710.6791.256 (0.615–2.565)0.3510.5320.908 (0.446–1.850)0.070.791
Alleles
 C128 (42.1%)121 (42.6%)71 (43.3%)
 T176 (57.9%)163 (57.4%)93 (56.7%)1.021 (0.736–1.416)0.0150.9021.050 (0.715–1.541)0.0610.8041.028 (0.698–1.516)0.020.887

Association of GCKR polymorphisms with clinical parameters in each group

The effects of GCKR rs780094 T or rs1260326 T allele on the clinical characteristics of NAFLD patients, NAFLD + CAD patients, and overall series were analyzed. As the results show in Table 4 and Table 5, the BMI value and serum FPG and TG levels were significantly lower in the rs780094 T and rs1260326 T allele carriers than in the non-carriers in the overall series (all p < 0.05), respectively. In addition, the serum level of FPG in rs780094 T carriers was higher than in non-carriers in the NAFLD + CAD group (p < 0.05). No marked difference of other clinical parameters was observed between the carriers and non-carriers of the rs780094 T or rs1260326 T allele in each group (p > 0.05). Multivariate logistic regression analysis showed that GCKR rs780094 is the dependent factor associated with FPG (p < 0.05).

Table 4.

Clinical characteristics of GCKR rs780094 T carriers and non-carriers in the study populationa

CharacteristicControlsNAFLDNAFLD + CAD
Carriers, n = 126Non-carriers,n = 26p valueCarriers, n = 117Non-carriers,n = 26p valueCarriers,n = 64Non-carriers,n = 18p value
BMI (kg/m2)22.91 ± 4.3624.44 ± 3.400.00126.35 ± 3.5027.05 ± 5.870.659#25.29 ± 2.42225.05 ± 3.2500.73#
FPG (mmol/L)4.45 ± 1.004.74 ± 1.000.0434.74 ± 0.654.86 ± 0.600.4055.18 ± 1.386.04 ± 3.160.001
TC (mmol/L)5.04 ± 1.005.09 ± 2.000.8605.49 ± 1.035.48 ± 1.280.9474.52 ± 1.236#4.69 ± 1.4300.616#
TG (mmol/L)1.08 ± 1.001.64 ± 2.000.0021.53 ± 1.011.75 ± 1.190.2521.49 ± 1.201.65 ± 1.330.751
HDL (mmol/L)1.29 ± 0.001.15 ± 0.000.2971.23 ± 0.271.23 ± 0.210.3251.03 ± 0.331.06 ± 0.280.396
LDL (mmol/L)3.07 ± 0.6973.23 ± 0.7540.295#3.25 ± 0.5473.36 ± 0.5860.349#2.74 ± 0.9482.78 ± 1.0830.866#
ALT (U/L)16.83 ± 11.0019.29 ± 11.000.99624.81 ± 23.6621.77 ± 40.270.95022.78 ± 18.5920.41 ± 13.030.320
AST (U/L)22.73 ± 6.0020.33 ± 6.000.24722.20 ± 8.5020.38 ± 13.110.52721.30 ± 18.6120.34 ± 15.070.251
r-GT (mmol/L)18.89 ± 13.0022.02 ± 14.000.14329.40 ± 26.7637.35 ± 23.790.29425.62 ± 16.7323.19 ± 16.950.406
ALP (mmol/L)69.92 ± 26.0077.06 ± 36.000.41766.87 ± 24.0774.80 ± 22.730.05679.85 ± 23.1084.77 ± 28.110.809
BIL (mmol/L)12.70 ± 6.0011.50 ± 4.000.43211.90 ± 5.1012.30 ± 3.600.58713.05 ± 6.1512.65 ± 5.200.397
Table 5.

Clinical characteristics of GCKR rs1260326 T carriers and non-carriers in the study populationa

CharacteristicControlsNAFLDNAFLD + CAD
Carriers, n = 129Non-carriers, n = 23p valueCarriers, n = 119Non-carriers, n = 24p valueCarriers, n = 67Non-carriers, n = 15p value
BMI (kg/m2)22.96 ± 4.3524.46 ± 3.810.00326.35 ± 3.4827.05 ± 6.060.665#25.27 ± 2.61625.11 ± 2.6480.836#
FPG (mmol/L)4.49 ± 1.004.76 ± 1.000.0544.75 ± 0.674.86 ± 0.560.5865.22 ± 1.436.02 ± 3.150.051
TC (mmol/L)5.03 ± 1.005.13 ± 2.000.8965.44 ± 1.005.89 ± 4.330.7134.51 ± 1.2244.76 ± 1.5060.496#
TG (mmol/L)1.08 ± 1.001.61 ± 2.000.0061.54 ± 1.001.75 ± 1.170.3661.46 ± 1.231.88 ± 1.210.710
HDL (mmol/L)1.29 ± 0.001.12 ± 0.000.2601.23 ± 0.281.23 ± 0.210.3691.03 ± 0.331.06 ± 0.300.564
LDL (mmol/L#)3.08 ± 0.7083.20 ± 0.7090.446#3.25 ± 0.5423.37 ± 0.610.324#2.72 ± 0.9352.85 ± 1.150.635#
ALT (U/L)16.81 ± 11.0019.70 ± 10.000.72924.81 ± 22.5221.77 ± 42.580.91822.70 ± 18.1422.63 ± 15.310.488
AST (U/L)20.30 ± 6.0022.91 ± 6.000.08522.20 ± 8.2720.38 ± 14.610.64721.29 ± 18.3920.36 ± 20.790.358
r-GT (mmol/L)18.79 ± 13.0022.30 ± 13.000.13529.40 ± 25.5937.35 ± 27.630.29625.79 ± 16.9222.91 ± 14.910.503
ALP (mmol/L)69.71 ± 26.0080.25 ± 32.000.25766.87 ± 22.9675.67 ± 23.550.04080.20 ± 23.3680.02 ± 30.260.802
BIL (mmol/L)12.50 ± 6.0011.70 ± 5.000.59311.90 ± 6.0012.30 ± 3.800.63613.10 ± 5.7011.10 ± 5.200.400

Discussion

NAFLD can be caused by many risk factors, such as diet, genetics and so on. Genetic susceptibility genes are a type of risk factor that can contribute to the development of NAFLD. In this study, we investigated the relationship between GCKR rs780094 and rs1260326 gene polymorphisms and the risk of NAFLD in a general Chinese Northern Han population, and the risk of CAD in NAFLD patients in this population for the first time. Our results showed that the genotype distributions of GCKR rs780094 and rs1260326 were not associated with the risk of NAFLD, and not associated with the risk of CAD in NAFLD patients in this Chinese Northern Han population. In addition to the BMI value and serum FPG and TG levels, the T allele of GCKR rs780094 and rs1260326 also did not significantly affect the clinical parameters in NAFLD patients and NAFLD + CAD patients. These data suggest that GCKR rs780094 and rs1260326 polymorphisms are not associated with the risk of NAFLD in the general Chinese Northern Han population, and GCKR rs780094 and rs1260326 polymorphisms do not increase the risk of CAD in NAFLD patients in the Chinese Northern Han population.

Accumulated studies have explored the effects of GCKR polymorphisms on the risk of NAFLD but the results have been inconsistent.16–18,32 Some previous studies and our study suggested that GCKR rs780094 and rs1260326 polymorphisms are not associated with the risk of NAFLD, the probable reason for this finding may be due to the ethnic differences among the NAFLD patients studied. The relationship between GCKR polymorphism and NAFLD or nonalcoholic steatohepatitis may be modified by the change and confluence of race, therefore the effects of a GCKR polymorphism may become more dominant and significant in European populations than they are in Asian populations. Therefore, the GCKR rs780094 T allele might increase the risk of NAFLD in European populations but not in Asian populations.15,17,18,33 A higher prevalence of NAFLD has been observed in some Indians that possess European ancestry, as opposed to those that possess Asian ancestry, which might explain the above hypothesis.34

The GCKR gene encodes GCKRP, an enzyme that plays a regulatory role in hepatic GCK activity.35,36 The GCKR rs1260326 functional variant P446L increases GCK activity by down-regulating the level of fructose 6-phosphate.37 Increased GCK activity is tightly accompanied by elevated hepatic glycolytic flux, de novo lipogenesis, and hepatic TG level.13,38 In this study, our results showed that the GCKR rs780094 and rs1260326 T allele were significantly associated with decreasing serum FPG level among the overall patients and the NAFLD + CAD patients. These results suggest that GCKR rs780094 and rs1260326 may repress the activity of GCKRP to decrease the serum FPG level.

Previous studies also investigated the association of GCKR polymorphisms with the risk of CAD. In a dietary intervention study, Shen et al.39 found that the GCKR rs1260326 T allele may increase the risk of atherosclerosis. Lian et al.40 demonstrated the positive association of GCKR rs780093 with the risk of CAD in an older Chinese Han population, through a case-control study. On the contrary, Járomi et al.41 did not detect any association of GCKR rs1260326 with the susceptibility of stroke. Bi et al.42 also found that the GCKR rs780094 SNP had no significant association with the incidence of CAD or stroke. SNPs of GCKR rs780094 showed strong linkage disequilibrium with rs1260326.24,39,40 Our results showed that GCKR rs780094 and rs1260326 genotype distributions were not associated with the occurrence of CAD in NAFLD patients but that the rs780094 T allele could decrease the serum FPG levels in the NAFLD + CAD group. These results indicated that GCKR rs780094 and rs1260326 were not the risk factor of CAD in NAFLD patients but that the rs780094 T allele could affect serum FPG level in the NAFLD + CAD patients of our Chinese Northern Han population.

Some limitations of our study should be acknowledged. First, all the subjects in this study were of the Chinese Northern Han population, so that the applicability of our conclusion to other ethnic populations requires further study for confirmation. Second, ultrasonography was used to diagnose NAFLD because of the difficulty in conducting liver biopsy. Third, our study did not grade the severity of NAFLD patients. Finally, a greater number of subjects should be included in the further studies to verify this conclusion in the Chinese Han population and other ethnic populations.

Conclusions

In summary, we investigated the relationship of GCKR rs780094 and rs1260326 gene polymorphisms with the risk of NAFLD as well as the risk of CAD in NAFLD patients in the Chinese Northern Han population. We concluded that GCKR rs780094 and rs1260326 gene polymorphisms did not associate with the risk of NAFLD nor with the risk of CAD in NAFLD patients in the Chinese Northern Han population. In addition, the GCKR rs780094 T and rs1260326 T alleles did affect the BMI value and serum FPG and TG levels in the overall subjects and NAFLD + CAD patients.

Abbreviations

ALP: 

alkaline phosphatase

ALT: 

alanine aminotransferase

AST: 

aspartate aminotransferase

BMI: 

body mass index

BIL: 

bilirubin

CAD: 

coronary artery disease

CI: 

confidence intervals

FPG: 

fasting plasma glucose

GCK: 

glucokinase

GCKR: 

glucokinase regulatory protein

GWAS: 

genome-wide association studies

HDL: 

high-density lipoprotein

LDL: 

low-density lipoprotein

NAFLD: 

nonalcoholic fatty liver disease

OR: 

odds ratio

PCR: 

polymerase chain reaction

r-GT: 

glutamyltransferase

SD: 

standard deviation

SNP: 

single nucleotide polymorphism

TC: 

total cholesterol

TG: 

triglyceride

Declarations

Acknowledgement

This study was supported by grants from the National Natural Science Foundation of China (31770837) and the Qingdao People’s Livelihood Science and Technology Plan (18-6-1-68-nsh).

Conflict of interest

The authors have no conflict of interests related to this publication.

Authors’ contributions

Study concept and design (YX and SX), acquisition of the data (HG, SL, ZZ, XY, and QL), analysis and interpretation of the data (HG, SL, and ZZ), drafting of the manuscript (HG and SL), critical revision of the manuscript for important intellectual content (YX and SX) supervision (YX and SX). All the authors read and approved the final manuscript.

References

  1. Danford CJ, Yao ZM, Jiang ZG. Non-alcoholic fatty liver disease: a narrative review of genetics. J Biomed Res 2018;32:389-400 View Article PubMed/NCBI
  2. Younossi Z, Anstee QM, Marietti M, Hardy T, Henry L, Eslam M. Global burden of NAFLD and NASH: trends, predictions, risk factors and prevention. Nat Rev Gastroenterol Hepatol 2018;15:11-20 View Article PubMed/NCBI
  3. Ibrahim SH, Hirsova P, Gores GJ. Non-alcoholic steatohepatitis pathogenesis: sublethal hepatocyte injury as a driver of liver inflammation. Gut 2018;67:963-972 View Article PubMed/NCBI
  4. Fiorucci S, Biagioli M, Distrutti E. Future trends in the treatment of non-alcoholic steatohepatitis. Pharmacol Res 2018;134:289-298 View Article PubMed/NCBI
  5. Younossi ZM, Stepanova M, Negro F, Hallaji S, Younossi Y, Lam B. Nonalcoholic fatty liver disease in lean individuals in the United States. Medicine (Baltimore) 2012;91:319-327 View Article PubMed/NCBI
  6. Liangpunsakul S, Chalasani N. Lipid mediators of liver injury in nonalcoholic fatty liver disease. Am J Physiol Gastrointest Liver Physiol 2019;316:G75-G81 View Article PubMed/NCBI
  7. Feng S, Dai Z, Liu AB, Huang J, Narsipur N, Guo G. Intake of stigmasterol and β-sitosterol alters lipid metabolism and alleviates NAFLD in mice fed a high-fat western-style diet. Biochim Biophys Acta Mol Cell Biol Lipids 2018;1863:1274-1284 View Article PubMed/NCBI
  8. Anstee QM, Day CP. The genetics of NAFLD. Nat Rev Gastroenterol Hepatol 2013;10:645-655 View Article PubMed/NCBI
  9. Vilar CP, Cotrim HP, Florentino GS, Bragagnoli G, Schwingel PA, Barreto CP. Nonalcoholic fatty liver disease in patients with coronary disease from a Brazil northeast area. Arq Gastroenterol 2015;52:111-116 View Article PubMed/NCBI
  10. Krishan S, Jain D, Bathina Y, Kale A, Saraf N, Saigal S. Non-invasive quantification of hepatic steatosis in living, related liver donors using dual-echo Dixon imaging and single-voxel proton spectroscopy. Clin Radiol 2016;71:58-63 View Article PubMed/NCBI
  11. Zhang Q, Zhang HM, Qi WQ, Zhang YG, Zhao P, Jiao J. 3.0T 1H magnetic resonance spectroscopy for assessment of steatosis in patients with chronic hepatitis C. World J Gastroenterol 2015;21:6736-6744 View Article PubMed/NCBI
  12. León-Mimila P, Vega-Badillo J, Gutiérrez-Vidal R, Villamil-Ramírez H, Villareal-Molina T, Larrieta-Carrasco E. A genetic risk score is associated with hepatic triglyceride content and non-alcoholic steatohepatitis in Mexicans with morbid obesity. Exp Mol Pathol 2015;98:178-183 View Article PubMed/NCBI
  13. Peter A, Stefan N, Cegan A, Walenta M, Wagner S, Königsrainer A. Hepatic glucokinase expression is associated with lipogenesis and fatty liver in humans. J Clin Endocrinol Metab 2011;96:E1126-E1130 View Article PubMed/NCBI
  14. Iynedjian PB. Molecular physiology of mammalian glucokinase. Cell Mol Life Sci 2009;66:27-42 View Article PubMed/NCBI
  15. Speliotes EK, Yerges-Armstrong LM, Wu J, Hernaez R, Kim LJ, Palmer CD. Genome-wide association analysis identifies variants associated with nonalcoholic fatty liver disease that have distinct effects on metabolic traits. PLoS Genet 2011;7:e1001324 View Article PubMed/NCBI
  16. Tan HL, Zain SM, Mohamed R, Rampal S, Chin KF, Basu RC. Association of glucokinase regulatory gene polymorphisms with risk and severity of non-alcoholic fatty liver disease: an interaction study with adiponutrin gene. J Gastroenterol 2014;49:1056-1064 View Article PubMed/NCBI
  17. Yang H, Chen G, Song C, Li D, Ma Q, Chen G. A novel index including SNPs for the screening of nonalcoholic fatty liver disease among elder Chinese: A population-based study. Medicine (Baltimore) 2018;97:e0272 View Article PubMed/NCBI
  18. Wang X, Liu Z, Wang K, Wang Z, Sun X, Zhong L. Additive effects of the risk alleles of PNPLA3 and TM6SF2 on non-alcoholic fatty liver disease (NAFLD) in a Chinese population. Front Genet 2016;7:140 View Article PubMed/NCBI
  19. Yang Z, Wen J, Tao X, Lu B, Du Y, Wang M. Genetic variation in the GCKR gene is associated with non-alcoholic fatty liver disease in Chinese people. Mol Biol Rep 2011;38:1145-1150 View Article PubMed/NCBI
  20. Shang XR, Song JY, Liu FH, Ma J, Wang HJ. GWAS-identified common variants with nonalcoholic fatty liver disease in Chinese children. J Pediatr Gastroenterol Nutr 2015;60:669-674 View Article PubMed/NCBI
  21. Lin YC, Chang PF, Chang MH, Ni YH. Genetic variants in GCKR and PNPLA3 confer susceptibility to nonalcoholic fatty liver disease in obese individuals. Am J Clin Nutr 2014;99:869-874 View Article PubMed/NCBI
  22. Kumar A, Shalimar , Walia GK, Gupta V, Sachdeva MP. Genetics of nonalcoholic fatty liver disease in Asian populations. J Genet 2019;98:29 View Article PubMed/NCBI
  23. Saxena R, Voight BF, Lyssenko V, Burtt NP, de Bakker PI, Chen H. Genome-wide association analysis identifies loci for type 2 diabetes and triglyceride levels. Science 2007;316:1331-1336 View Article PubMed/NCBI
  24. Sparsø T, Andersen G, Nielsen T, Burgdorf KS, Gjesing AP, Nielsen AL. The GCKR rs780094 polymorphism is associated with elevated fasting serum triacylglycerol, reduced fasting and OGTT-related insulinaemia, and reduced risk of type 2 diabetes. Diabetologia 2008;51:70-75 View Article PubMed/NCBI
  25. Kawaguchi T, Shima T, Mizuno M, Mitsumoto Y, Umemura A, Kanbara Y. Risk estimation model for nonalcoholic fatty liver disease in the Japanese using multiple genetic markers. PLoS One 2018;13:e0185490 View Article PubMed/NCBI
  26. Petit JM, Masson D, Guiu B, Rollot F, Duvillard L, Bouillet B. GCKR polymorphism influences liver fat content in patients with type 2 diabetes. Acta Diabetol 2016;53:237-242 View Article PubMed/NCBI
  27. Di Costanzo A, Pacifico L, D’Erasmo L, Polito L, Martino MD, Perla FM. Nonalcoholic fatty liver disease (NAFLD), but not its susceptibility gene variants, influences the decrease of kidney function in overweight/obese children. Int J Mol Sci 2019;20:E4444 View Article PubMed/NCBI
  28. Simons PIHG, Simons N, Stehouwer CDA, Schalkwijk CG, Schaper NC, Brouwers MCGJ. Association of common gene variants in glucokinase regulatory protein with cardiorenal disease: A systematic review and meta-analysis. PLoS One 2018;13:e0206174 View Article PubMed/NCBI
  29. Rickham PP. Human experimentation. Code of ethics of the world medical association. declaration of helsinki. Br Med J 1964;2:177 View Article PubMed/NCBI
  30. Kotronen A, Westerbacka J, Bergholm R, Pietiläinen KH, Yki-Järvinen H. Liver fat in the metabolic syndrome. J Clin Endocrinol Metab 2007;92:3490-3497 View Article PubMed/NCBI
  31. Xu Y, Lu LL, Liu SS, Du SX, Zhu HL, Dong QJ. Apolipoprotein A5 gene polymorphisms are associated with non-alcoholic fatty liver disease. Hepatobiliary Pancreat Dis Int 2018;17:214-219 View Article PubMed/NCBI
  32. Flores YN, Velázquez-Cruz R, Ramírez P, Bañuelos M, Zhang ZF, Yee HF. Association between PNPLA3 (rs738409), LYPLAL1 (rs12137855), PPP1R3B (rs4240624), GCKR (rs780094), and elevated transaminase levels in overweight/obese Mexican adults. Mol Biol Rep 2016;43:1359-1369 View Article PubMed/NCBI
  33. Kitamoto A, Kitamoto T, Nakamura T, Ogawa Y, Yoneda M, Hyogo H. Association of polymorphisms in GCKR and TRIB1 with nonalcoholic fatty liver disease and metabolic syndrome traits. Endocr J 2014;61:683-689 View Article PubMed/NCBI
  34. Bamshad M, Kivisild T, Watkins WS, Dixon ME, Ricker CE, Rao BB. Genetic evidence on the origins of Indian caste populations. Genome Res 2001;11:994-1004 View Article PubMed/NCBI
  35. Hayward BE, Dunlop N, Intody S, Leek JP, Markham AF, Warner JP. Organization of the human glucokinase regulator gene GCKR. Genomics 1998;49:137-142 View Article PubMed/NCBI
  36. Grimsby J, Coffey JW, Dvorozniak MT, Magram J, Li G, Matschinsky FM. Characterization of glucokinase regulatory protein-deficient mice. J Biol Chem 2000;275:7826-7831 View Article PubMed/NCBI
  37. Beer NL, Tribble ND, McCulloch LJ, Roos C, Johnson PR, Orho-Melander M. The P446L variant in GCKR associated with fasting plasma glucose and triglyceride levels exerts its effect through increased glucokinase activity in liver. Hum Mol Genet 2009;18:4081-4088 View Article PubMed/NCBI
  38. Rees MG, Wincovitch S, Schultz J, Waterstradt R, Beer NL, Baltrusch S. Cellular characterisation of the GCKR P446L variant associated with type 2 diabetes risk. Diabetologia 2012;55:114-122 View Article PubMed/NCBI
  39. Shen H, Pollin TI, Damcott CM, McLenithan JC, Mitchell BD, Shuldiner AR. Glucokinase regulatory protein gene polymorphism affects postprandial lipemic response in a dietary intervention study. Hum Genet 2009;126:567-574 View Article PubMed/NCBI
  40. Lian J, Guo J, Chen Z, Jiang Q, Ye H, Huang X. Positive association between GCKR rs780093 polymorphism and coronary heart disease in the aged Han Chinese. Dis Markers 2013;35:863-868 View Article PubMed/NCBI
  41. Járomi L, Csöngei V, Polgár N, Szolnoki Z, Maász A, Horvatovich K. Functional variants of glucokinase regulatory protein and apolipoprotein A5 genes in ischemic stroke. J Mol Neurosci 2010;41:121-128 View Article PubMed/NCBI
  42. Bi M, Kao WH, Boerwinkle E, Hoogeveen RC, Rasmussen-Torvik LJ, Astor BC. Association of rs780094 in GCKR with metabolic traits and incident diabetes and cardiovascular disease: the ARIC Study. PLoS One 2010;5:e11690 View Article PubMed/NCBI