Introduction
Non-alcoholic fatty liver disease (NAFLD) is the clinicopathological syndrome, characterized by excessive intracellular fat deposition that occurs with the exclusion of alcohol and other specific liver damaging factors, such as viruses, drugs, autoimmune factors, and genetic factors.1 According to the clinical statistic results, coronary artery disease (CAD) causes the highest mortality in the patients with NAFLD, and the incidence of CAD in patients with NAFLD has increased significantly in recent years.2,3 A previous study showed that NAFLD can act as one of the independent risk factors for CAD, after excluding hypertension, obesity, diabetes, and other factors.4 Both NAFLD and CAD are influenced by the interaction of genetic and environmental factors, and they share some pathological mechanisms, such as obesity, insulin resistance, dyslipidemia, inflammation, and oxidative stress.3,5,6 A large number of studies have investigates the role of NAFLD in the development of CAD. NAFLD was reported to be related to the formation of coronary artery plaques and impaired coronary blood flow reserve.4,7,8 In addition, genetic pathogenesis in NAFLD and CAD is drawing more attention, and the gene polymorphisms of ADIPOQ, LEPR, APOC3, PPAR, SREBP, TM6SF2, and MTTP have been reported as risk factors for the development of NAFLD and CAD.5
TCF7L2 is located on chromosome 10q25 and encodes the transcription factor-4 (TCF-4) (also known as lymphocyte factor-4), which is an important member of the transcription factor family.9–11 The TCF-4 family can participate in a variety of physiological pathways in different types of cells.12,13 In 2009, Musso et al.14 first reported that TCF7L2 rs7903146 T allele frequency was higher in NAFLD patients than in healthy controls. In addition, they found that TCF7L2 rs7903146 T allele carriers possessed higher levels of low-density lipoprotein (LDL), very low-density lipoprotein, triglyceride (TG), and cytokeratin 18 than non-carriers in the NASH subgroup. Subsequently, Giovanni et al.15 reported that TCF7L2 rs7903146 T allele carriers had the prolonged elevation of post-meal glucose-dependent insulinotropic polypeptide (GIP). The elevated GIP was demonstrated in other studies to exacerbate liver steatosis and increase the levels of postprandial resistin and free fatty acids.16,17 Besides, some evidence has suggested that TCF7L2 rs7903146 is associated with serum lipid indexes, including TG, total cholesterol (TC), LDL, and apolipoprotein B, all of which are associated with risk of CAD development.18,19 Corella et al.20 demonstrated that the relationship between TCF7L2 rs7903146 and serum lipid profiles was influenced by adherence to the MedDiet. When adherence was low, the levels of TC and LDL were higher in TT carriers than CT/CC carriers. But when adherence was high, no difference was observed between TT carriers and CT/CC carriers. Up to now, the association of TCF7L2 rs7903146 with the risk of developing NAFLD and CAD was investigated only in a limited amount of studies, and no related report was available in Chinese populations. In addition, the effect of TCF7L2 rs7903146 polymorphism on serum lipid profiles is still unclear.21–23 In consideration of the increased incidence of NAFLD and CAD in China, illuminating the association of TCF7L2 rs7903146 with the risk of NAFLD and CAD in Chinese may provide some new ideas for the treatment of NAFLD and CAD.
The aim of this study was to investigate the relationship of TCF7L2 rs7903146 with the risk of developing NAFLD and CAD in a Chinese Han population, and explore the effect of TCF7L2 rs7903146 on the serum lipid profiles.
Methods
Subjects
This study was conducted according to the principles of the declaration of Helsinki and its appendices.24 Subjects were composed of 143 patients with NAFLD, 159 patients with CAD, 131 patients with NAFLD + CAD, and 212 healthy controls, who were enrolled at the Qingdao Municipal Hospital (China) from December 2017 to December 2018. The Guidelines for the Diagnosis and Treatment of Non-alcoholic Fatty Liver Disease (2010 Revision) and the Guidelines for the Diagnosis and Treatment of Acute ST-Segment Elevation Myocardial Infarction (2015 Version) were used as exclusion and inclusion criteria for NAFLD and CAD. Patients with NAFLD were included from the Department of Gastroenterology, in a random manner, and diagnosed by B-type ultrasonography. Combining the medical history and laboratory testing, subjects with alcoholic fatty liver, viral hepatitis, autoimmune liver disease, drug-induced liver disease and hereditary metabolic disease were excluded. Patients with CAD were included from the Department of Cardiology, in a random manner, and diagnosed by coronary angiography. We divided the CAD patients into with or without NAFLD groups, respectively. The healthy controls were included from the Physical Examination Center of the Qingdao Municipal Hospital, in a random manner. All the subjects signed informed consent forms after participating in this study.
Biochemical analyses
The demographic data (gender, age, height and weight) were obtained by questionnaire. The body mass index (BMI) was calculated by weight (kg)/height (m)2. Fasting venous blood samples of each subject were collected after a 12-h overnight fasting and kept in ethylene diamine tetraacetic acid-containing tubes. The biochemical data, such as TG, TC, fasting plasma glucose (FPG), high-density lipoprotein (HDL), LDL, alkaline phosphatase (ALP), gamma-glutamyltransferase (GGT), and total bilirubin (TBIL) were measured through standard laboratory method.
Genomic DNA extraction and genotyping
Genomic DNA was extracted from each blood sample and stored at -80oC until use. The primers for PCR amplification were designed and synthesized by Bomiao Biotechnology (Beijing, China) as: 5’-ACGTTGGATGAACTAAGGGTGCCTCATACG-3’ and 5’-ACGTTGGATGGCCTCAAAACCTAGCACAGC-3’. The detailed process of PCR amplification consists of the following steps: an initial denaturation at 94oC for 5 m, followed by 45 cycles of denaturing at 94oC for 20 s, annealing at 56oC for 30 s, and extending at 72oC for 1 m. After that, a final extension at 72 oC was conducted for 5 m. The genotype TCF7L2 rs7903146 was detected by DNA sequence using the MALDI-TOF-MS (MassARRAY System; Agena Bioscience, Shanghai, China) and raw data was acquired by TYPER4.0.
Statistical analysis
The measurement data were expressed as mean ± standard deviation (SD) for those conforming to the normal distribution, otherwise as median quartile. Homogeneity of variance was tested between each of two groups. The independent samples t-test was used for data conforming to normal distribution with homogeneity variances, while the rest were tested by rank sum test. Chi-square test was used to analyze the gender distribution differences. Hardy-Weinberg equilibrium was performed by chi-square test. The distribution of genotypes and allele frequency between two groups were also tested by chi-square. The correlation between genotype and risk of diseases were estimated by binary logistic regression modeling. The demographic and biochemical data of different genotypes were tested by independent sample t-test and rank sum test. A p-value of <0.05 was considered to indicate statistical significance. All the statistical analyses were conducted by SPSS 23.0 software (IBM Corp., Armonk, NY, USA).
Results
Demographic and biochemical characteristics of subjects
A total of 457 subjects were included in this study. As the results show in Table 1, differences of gender distribution were observed between the NAFLD group and control group, and between the NAFLD + CAD group and the control group. The NAFLD group had higher levels of BMI, TC, TG, LDL, alanine aminotransferase, GGT and FPG than the control group, while age and HDL level were lower in the NAFLD group compared to the control group (all p < 0.05). No significant differences were observed for aspartate aminotransferase, ALP and TBIL between the NAFLD group and the control group. Age, BMI and levels of TG, alanine aminotransferase, ALP, GGT, and FPG in the CAD group were high compared to those in the control group, and the levels of TC, HDL, and LDL were low in the CAD group compared to those in the control group (all p < 0.05). There were higher age, BMI, and levels of TG, alanine aminotransferase, ALP, GGT and FPG, and lower levels of TC, HDL and LDL in the NAFLD + CAD group compared to the control group (all p < 0.05). The NAFLD + CAD group had higher age and levels of ALP and FPG, and lower BMI and levels of HDL, LDL, alanine aminotransferase and TC than the NAFLD group (all p < 0.05).
Table 1.Demographic and biochemical characteristics of subjects
| Control group, n = 212 | NAFLD group, n = 143 | CAD group, n = 159 | NAFLD + CAD group, n = 131 | p1 | p2 | p3 | p4 |
Male (%) | 121 (57.06) | 100 (69.93) | 96 (60.38) | 90 (68.70) | 0.014 | 0.523 | 0.032 | 0.826 |
Age (yr) | 47 (41-57.75) | 42 (37-45) | 66 (61-74) | 62 (55-67) | <0.001 | <0.001 | <0.001 | <0.001 |
BMI (kg/m2) | 23.90 (21.82-25.91) | 26.61 ± 5.50 | 25.28 ± 3.25 | 25.19 ± 2.54 | <0.001 | <0.001 | <0.001 | <0.001 |
TC (mmol/L) | 5.11 (4.59-5.74) | 5.48 ± 0.82 | 4.64 ± 1.28 | 4.31 (3.76-5.54) | <0.001 | <0.001 | <0.001 | <0.001 |
TG (mmol/L) | 1.12 (0.87-1.69) | 1.54 (1.16-2.20) | 1.43 (1.03-1.85) | 1.46 (1.01-2.18) | <0.001 | 0.002 | <0.001 | 0.328 |
HDL (mmol/L) | 1.30 ± 0.29 | 1.21 (1.08-1.35) | 1.01 (0.86-1.16) | 1.02 (0.85-1.17) | 0.031 | <0.001 | <0.001 | <0.001 |
LDL (mmol/L) | 3.10 (2.67-3.58) | 3.25 ± 0.61 | 2.64 (2.09-3.41) | 2.64 (2.12-3.43) | 0.047 | <0.001 | <0.001 | <0.001 |
ALT (U/L) | 18.28 (12.87-25.87) | 23.80 (18.42-43.65) | 20.73 (14.72-31.36) | 22.67 (15.63-33.44) | <0.001 | 0.015 | 0.001 | 0.021 |
AST (U/L) | 20.79 (18.15-24.68) | 21.64 (18.33-26.78) | 21.70 (16.99-31.08) | 22.17 (16.80-32.72) | 0.139 | 0.335 | 0.107 | 0.659 |
ALP (U/L) | 71.17 ± 20.35 | 68.63 (58.99-83.39) | 79.78 ± 22.87 | 82.50 (70.52-98.63) | 0.785 | <0.001 | <0.001 | <0.001 |
GGT (U/L) | 22.12 (16.19-30.21) | 30.63 (20.77-47.73) | 26.35 (18.16-36.84) | 27.70 (18.54-39.98) | <0.001 | 0.001 | <0.001 | 0.084 |
TBIL (μmol/l) | 13.20 (9.43-16.88) | 11.9 (9.9-14.9) | 13.50 (10.20-18.28) | 13.10 (10.48-16.30) | 0.398 | 0.169 | 0.448 | 0.111 |
Glu (mmol/L) | 4.48 (4.04-4.96) | 4.79 (4.50-5.18) | 5.20 (4.55-6.43) | 5.42 (4.75-6.40) | <0.001 | <0.001 | <0.001 | <0.001 |
TCF7L2 rs7903146 genotype and allele frequency distribution
The genotype frequency distribution of TCF7L2 rs7903146 was in accordance with Hardy-Weinberg equilibrium in each group (all p > 0.05) (Table 2). As the results show in Table 3, there was no significant difference in the genotype distribution of TCF7L2 rs7903146 (CC vs. CT + TT) between the NAFLD group and the control group, the CAD group and the control group, the NAFLD + CAD group and the control group, and the NAFLD + CAD group and the NAFLD group (all p > 0.05). In addition, the difference of TCF7L2 rs7903146 allele frequency distribution in each group was also not significant (all p > 0.05). When we analyzed the association of TCF7L2 rs7903146 genotype and the risk of NAFLD, CAD, and NAFLD + CAD, the results suggested that the CT + TT genotype did not increase the risk of developing NAFLD, CAD, and NAFLD + CAD in the general population, and the CT + TT genotype also did not increase the risk of developing CAD in the patients with NAFLD (Table 4).
Table 2.Result of Hardy-Weinberg Equilibrium analysis
Gene locus | Group | X2 | p |
rs7903146 | NAFLD group | 0.39 | 0.53 |
| Control group | 0.63 | 0.43 |
| CAD group | 0.03 | 0.86 |
| NALFD + CAD group | 0.07 | 0.81 |
Table 3.Distribution of genotype and allele frequency of TCF7L2 rs7903146 in each group*
| Genotype frequency, n (%) | p1 | p2 | p3 | p4 | Allele frequency, n (%) | p1 | p2 | p3 | p4 |
CC | CT + TT | C | T |
NAFLD | 125 (89.93) | 14 (10.07) | 0.93 | 0.28 | 0.06 | 0.08 | 263 (94.60) | 15 (5.40) | 0.90 | 0.22 | 0.06 | 0.06 |
Control | 190 (89.62) | 22 (10.38) | 402 (94.81) | 22 (5.19) |
CAD | 134 (85.90) | 22 (14.10) | 289 (92.63) | 23 (7.37) |
NAFLD + CAD | 110 (95.65) | 5 (4.35) | 225 (97.83) | 5 (2.17) |
Table 4.Association between genotype and the risk of NAFLD, CAD, and CAD & NAFLD*
| OR (95% CI) | p1 | OR (95% CI) | p2 | OR (95% CI) | p3 | OR (95% CI) | p4 |
CC | 1 | 0.93 | 1 | 0.29 | 1 | 0.07 | 1 | 0.09 |
CT + TT | 0.97 (0.48-1.96) | 1.42 (0.75-2.67) | 0.39 (0.15-1.07) | 0.41 (0.14-1.16) |
Association of TCF7L2 rs7903146 genotype with the risk of NAFLD in the non-obese subjects
The difference of clinical parameters between the T allele carriers and non-carriers in each group was analyzed, and no significant difference was observed in the NAFLD group, the CAD group, and the NAFLD + CAD group (data not shown). However, the BMI of the T allele carriers was lower than in non-carriers in the control group (p = 0.027). To investigate the potential relationship of TCF7FL2 rs7903146 genotype with the BMI value, we divided the whole subject population into two groups: obese subjects (BMI ≥25) and non-obese subjects (BMI <25). In the obese group, no significant difference was observed for the genotype distribution of TCF7L2 rs7903146 (CC vs. CT + TT) between the NAFLD patients and non-NAFLD patients. In the non-obese group, distribution of the CT + TT genotype and CC genotype was significantly different (p = 0.035) (Table 5). Logistic regression model analysis suggested that the CT + TT genotype was a protective factor against the development of NAFLD in the non-obese subjects (odds ratio=0.359, 95% confidence interval: 0.134-0.961, p = 0.041). After adjustment for gender and age, and serum lipids, the CT + TT genotype was still a protective factor against the development of NAFLD in the non-obese subjects (odds ratio = 0.245, 95% confidence interval: 0.072-0.837, p = 0.025) (Table 6).
Table 5.Distribution of genotype frequency and allele frequency of TCF7L2 rs7903146 in obese and non-obese subjects*
| | NAFLD group, n (%) | Non-NAFLD group, n (%) | X2 | p |
Obese | | | | |
| Genotype | | | | |
| CC | 151 (91.5) | 149 (90.9) | 0.045 | 0.832 |
| CT + TT | 14 (8.5) | 15 (9.1) |
| Allele | | | | |
| C | 316 (95.76) | 313 (95.43) | 0.043 | 0.836 |
| T | 14 (4.24) | 15 (4.57) |
Non-obese | | | | |
| Genotype | | | | |
| CC | 84 (94.38) | 175 (85.78) | 4.466 | 0.035 |
| CT + TT | 5 (5.63) | 29 (14.22) |
| Allele | | | | |
| C | 171 (96.07) | 378 (92.65) | 2.451 | 0.117 |
| T | 7 (3.93) | 30 (7.35) |
Table 6.Association of TCF7L2 rs7903146 genotype and NAFLD in non-obese subjects*
| Unadjusted | | Adjusted* |
| OR (95% CI) | P | | OR (95% CI) | P |
CC | 1 | 0.041 | CC | 1 | 0.025 |
CT + TT | 0.359 (0.134-0.961) | CT + TT | 0.245 (0.072-0.837) |
Discussion
Genetic factors significantly participate in the development and progression of NAFLD and CAD.25,26 In this study, we investigated the relationship of TCF7L2 rs7903146 with the risk of developing NAFLD and CAD in a Chinese Han population. As the results showed, there were no differences in the TCF7L2 rs7903146 genotype and allele distribution between the NAFLD group and the control group, the CAD group and the control group, the NAFLD + CAD group and the control group, and the NAFLD + CAD group and the NAFLD group. In addition, the CT + TT genotype did not increase the risk of developing NAFLD, CAD, and NALFD + CAD in the general population, and did not increase the risk of developing CAD in the patients with NAFLD. However, the TCF7L2 rs7903146 CT + TT genotype was a protective factor for NAFLD in the non-obese subjects.
As a component of the Wnt/beta-catenin signaling pathway, TCF7L2 participates in a variety of physiological pathways, such as cell proliferation, differentiation, apoptosis, and oxidative stress, all the effect of TCF7L2 which can be exerted throughout growth, zoning, heterogenesis, and other intrinsic metabolism processes of liver.27–29 The relationship of TCF7L2 with metabolism was first found by Grant et al.,30 who found that the microsatellite of TCF7L2 intron 3, DG10S478, was associated with type 2 diabetes in an Icelandic population. The gene’s expression product, a high-mobility family box, contains transcription factors that influence blood glucose homeostasis. Subsequent studies revealed several potential functions of TCF7L2, including decreasing insulin secretion, impairing incretin effect, and increasing insulin resistance.31–33 In addition, TCF7L2 rs7903146 has been reported to be relevant to lipid disorders, NAFLD, obesity and hypertension.14,18,34 Corella et al.20 reported the correlation between TCF7L2 rs7903146 and multiple lipid indexes, but the potential molecular mechanism by which TCF7L2 rs7903146 induces abnormal hepatic metabolism remains unknown.
Some researchers reported that TCF7L2 rs7930146 can regulate open chromatin and genetic transcription in pancreas, thus increasing the risk of diabetes.31,35 In addition, Dorota et al.36 found that the expression of TCF7L2 in liver may be regulated by weight loss. However, the evidence of how TCF7L2 rs7903146 participates in hepatic lipid metabolism remains unclear. Musso et al.14 found that the TCF7L2 rs7903146 CT/TT genotype was a risk factor for NAFLD. In that study, the results indicated that the TCF7L2 rs7903146 polymorphism did not associate with the risk of developing NAFLD, CAD, or NAFLD + CAD in the general population, and did not associate with the risk of developing NAFLD + CAD. Except for BMI in the control group, the differences of clinical parameters between the TCF7L2 rs7903146 T allele carriers and non-carriers in each group were also not significant. Further analysis found that the distribution of the CT + TT genotype and the CC genotype was significantly different between the NAFLD patients and non-NAFLD patients in the non-obese group. After adjustment for gender and age, and serum lipids, the CT + TT genotype was still a protective factor for the development of NAFLD in the non-obese subjects. Previous study determined that the minor allele frequency of the rs7903146 T allele varied according to the difference of region (the minor allele frequencies in Vietnam, Estonia, Sweden, Caucasia, India, East Asia were 4.63, 19.82, 22.5, 29, 28.1, and 3.2, respectively).37 In this study, the minor allele frequency of rs7903146 was 4.63, which was similar to the minor allele frequency of East Asia. Development of NAFLD was affected by ethnic, genetic, dietary, and environmental factors; therefore, the difference of the effect of TCF7L2 rs7903146 polymorphism on the risk of developing NAFLD in Chinese and other countries may be affected by ethnic and genetic factors.
There were several limitations in this study. First, all the NAFLD patients were diagnosed by ultrasound rather than the liver biopsy; therefore, the diagnosis of NAFLD may not be very accurate. Second, the sample size was relatively small, which may affect the conclusion. Third, our conclusion cannot reflect the association between TCF7L2 rs7903146 polymorphism and the risk of developing NAFLD in other countries, due to all the subjects in this study being of the Chinese Han population.
Conclusions
This study investigated the relationship of TCF7L2 rs7903146 polymorphism with the risk of NAFLD, CAD, and NAFLD + CAD in a Chinese Han population, for the first time. The results suggest that TCF7L2 rs7903146 was not associated with the risk of developing NAFLD, CAD, and NAFLD + CAD in our Chinese Han population, and TCF7L2 rs7903146 did not affect the serum lipid metabolism. In the non-obese subpopulation, the TCF7L2 rs7903146 CT + TT genotype was a protective factor for the development of NAFLD. More subjects of other ethnicity should be included to further investigate the association of TCF7L2 rs7903146 with the risk of NAFLD, and the potential mechanism underlying the TCF7L2 rs7903146 variant’s affect the lipid metabolism should be illustrated.
Abbreviations
- ALP:
alkaline phosphatase
- BMI:
body mass index
- CAD:
coronary artery disease
- FPG:
fasting plasma glucose
- GGT:
gamma-glutamyltransferase
- GIP:
glucose-dependent insulinotropic polypeptide
- HDL:
high-density lipoprotein
- LDL:
low-density lipoprotein
- NAFLD:
non-alcoholic fatty liver disease
- SD:
standard deviation
- TBIL:
total bilirubin
- TC:
total cholesterol
- TCF-4:
transcription factor-4
- TG:
triglyceride
Declarations
Funding
This study was supported by a grant from the National Natural Science Foundation of China (31770837).
Conflict of interest
The authors have no conflict of interests related to this publication.
Authors’ contributions
Study concept and design (YX and SL), acquisition of the data (XY, WJ, JZ, and MW), analysis and interpretation of the data (YX and WJ), drafting of the manuscript: (XY), critical revision of the manuscript for important intellectual content (YX and LS), supervision (YX). All the authors read and approved the final manuscript.