Kuan-Chun Hsueh1,2
Shun-Fa Yang
Yung-Po Liaw- 1Department of Surgery, Tungs’ Taichung Metro Harbor Hospital, Taichung, Taiwan
- 2Department of Post-Baccalaureate Medicine, College of Medicine, National Chung Hsing University, Taichung, Taiwan
- 3Department of Public Health and Institute of Public Health, Chung Shan Medical University, Taichung, Taiwan
- 4Institute of Medicine, Chung Shan Medical University, Taichung, Taiwan
- 5Department of Medical Research, Chung Shan Medical University Hospital, Taichung, Taiwan
- 6Department of Medical Imaging, Chung Shan Medical University Hospital, Taichung City, Taiwan
Type 2 diabetes (T2D) and liver cirrhosis remain significant public health threats in Taiwan. These conditions are reported to be associated with the rs738409 polymorphism of the patatin-like phospholipase domain-containing protein three gene (PNPLA3) in European populations. We assessed the effect of T2D and PNPLA3 rs738409 polymorphism on liver cirrhosis among Taiwan Biobank (TWB) participants. In total, 17,985 participants in TWB had their health records linked to the National Health Insurance Research Database (NHIRD). Participants included those who visited the assessment centers between 2008 and 2015, with an age range between 30 and 70 years of age. We performed logistic regression analysis to investigate the odds ratios (OR) for liver cirrhosis among participants based on the T2D status and rs738409 genotypes. Genotyping was performed using the Axiom Genome-Wide TWB Array Plate. In our analysis, 150 of the 17,619 eligible participants were identified as cirrhosis cases. Based on the univariate analysis, liver cirrhosis was positively associated with T2D (OR, 1.83; 95% CI 1.23–2.70) whereas, the variant rs738409 was not (regardless of the genetic model). The variant and T2D, however, showed significant interactions in the additive, genotype, and dominant models (p values of 0.0302, 0.0395, and 0.0455, respectively). We observed a statistically significant association between T2D and liver cirrhosis and variant rs738409 with an OR of 1.71 (95% CI, 1.03–2.84) for individuals carrying a G allele compared to those with a C allele and 2.92 (95% CI 1.07–7.99) for GG compared to CC individuals. According to our study, Taiwanese adults with T2D and the rs738409 GG genotype are more likely to develop liver cirrhosis.
Introduction
Globally, liver diseases are one of the most important public health issues (Asrani et al., 2019). They remain a leading cause of morbidity and mortality. Cirrhosis is a potentially fatal complication of chronic liver disease, which may manifest 20–40 years after infection (Smith et al., 2019). It is the 11th leading cause of death globally (Asrani et al., 2019) and has a considerable economic impact (Smith et al., 2019). It is commonly caused by hepatitis, alcoholism, and non-alcoholic liver disease, but other causes remain unknown in approximately 5%–30% of those with this condition. (Nalbantoglu and Jain, 2019; Patel et al., 2019). As a result, exploring other possible variables may increase our understanding of the disease pathology.
Cirrhosis is strongly linked to viral hepatitis in several countries, including Taiwan (Mokdad et al., 2014; Sepanlou et al., 2020; Yang et al., 2020). In spite of this, the underlying causes and disease burden vary across demographic groups and geographic regions (Sepanlou et al., 2020). Diabetes is among the chronic conditions that have previously been linked to liver cirrhosis and mortality in various ethnic groups (Holstein et al., 2002; Quintana et al., 2016; Goh et al., 2017), even though the underlying pathogenic mechanisms are not entirely clear (Tolman et al., 2007). Mortality risk was found to be higher in patients with compensated liver cirrhosis who had T2D (Quintana et al., 2016).
Aside from the common causes, genetic factors have also been associated with liver cirrhosis (Stickel et al., 2020). However, the association between liver cirrhosis and a number of polymorphic variants has been inconsistent (Basyte-Bacevice et al., 2019; Chen et al., 2020). Among the previously identified polymorphisms, rs738409 variant of the PNPLA3 gene has been strongly associated with cirrhosis and indices of liver damage, especially in Caucasians (Romeo et al., 2010). This variant has also been linked to liver-related events and death among Italians with nonalcoholic fatty liver disease (NAFLD) (Grimaudo et al., 2020), which is becoming more common in the eastern regions of Asia (Zhang et al., 2015). Despite its role in aggravating liver steatosis, the rs738409 variant has also been associated with a decreased risk of T2D in metabolically unhealthy middle-aged and elderly Chinese populations (Xia et al., 2019).
On the contrary, an earlier study found no association between diabetes and the PNPLA3 variant (Speliotes et al., 2010). Of note, diabetes and cirrhosis remain public health concerns in Taiwan. These two conditions have been associated with rs738409 polymorphism mainly in other populations, as cited above. Considering that the rs738409 variant has not been studied extensively among Asians, we examined the effect of this variant in combination with T2D on liver cirrhosis in TWB participants.
Materials and Methods
Data Sources
Through the data management system of the Health and Welfare Data Science Center (HWDC), we linked participants’ genotype data in TWB to their medical records in the NHIRD using personal identification numbers. Participants included those who visited the assessment centers between 2008 and 2015, with an age range between 30 and 70 years of age. Written informed consent was obtained from each participant before they were assessed. Disease information in the NHIRD was available from 1998 to 2015. Approval was obtained from the Institutional Review Board of Chung Shan Medical University. All authors had access to the study data and had reviewed and approved the final manuscript.
Definition of the Variables
Data were available for 17,957 subjects in TWB. We excluded those with an incomplete questionnaire (n = 271) and genotyping information (n = 27), as well as cases of cirrhosis identified before the diagnosis of diabetes (n = 40). Finally, we identified 150 participants as cirrhosis cases and 17,469 as control individuals. Liver cirrhosis was diagnosed based on the International Classification of Diseases, Ninth Revision, Clinical Modification (ICD-9-CM), codes 571.2, 571.5, 571.6, 572.2, 572.3, 572.4, 572.8, 573.0) and either one-time hospitalization or two outpatient visits. Detailed information on the diagnoses of T2D has been described elsewhere (Nfor et al., 2018).
Genotyping
The PNPLA3 rs738409 was genotyped using the custom Taiwan Biobank chip (Axiom™ Genome-Wide Array Plate System (Affymetrix, Santa Clara, CA, United States). The SNP call rate was >95% and was in Hardy-Weinberg equilibrium (p-value >1.0 × 10–3. The minor allele frequency (MAF) was >0.05.
Statistical Analysis
Our data analyses were conducted with the software package PLINK 1.09 and the SAS software (version 9.4). We compared baseline parameters using the x2 test. We investigated the effects of T2D and PNPLA3 rs738409 genetic variant on liver cirrhosis based on univariate and multivariate logistic regression analyses. The additive, dominant, and recessive models were used to determine the allele and genotype-specific ORs. Adjusted variables included sex, age, educational level, smoking, alcohol intake, body mass index (BMI), physical activity, hepatitis C virus (HCV), hepatitis B virus (HBV), and hyperlipidemia. In order to reduce potential bias, patients with cirrhosis before the diagnosis of diabetes were excluded from the final logistic regression models. We tested for interaction by including the term rs738409*T2D in the regression models. The ORs were estimated along with their 95% CI.
Results
Table 1 describes the general characteristics of study participants. Overall, 150 of the 17,469 eligible participants were identified as cirrhosis cases. Among the cases of cirrhosis, 60 were of the rs738409 CC genotype, 62 were of the CG genotype, and 28 were of the GG genotype. About 21.33% (n = 32) of participants with liver cirrhosis had T2D. According to the univariate analysis, liver cirrhosis was positively associated with T2D (OR, 1.83; 95% CI 1.23–2.70) but the variant rs738409 was not despite the genetic model used (Table 2). Multivariate analyses, however, found no significant correlation between cirrhosis and rs738409 or T2D regardless of the model used. There was, however, a significant interaction between the variant and T2D in the additive, genotype, and dominant models (with p values of 0.0302, 0.0395, and 0.0455, respectively).
TABLE 1. Basic data of study participants.
TABLE 2. Association of liver cirrhosis with PNPLA3 rs738409 and T2D.
According to the univariate analyses, HCV and HBV infections both contributed significantly to cirrhosis risk, with ORs of 6.89 (95% CI, 4.3–11.05) and 4.19 (95% CI, 2.98–5.88), respectively. In spite of the analysis model, both variables remained significantly associated with liver cirrhosis as shown in Table 2. For subjects aged 40–49, 50–59, and 60–70, the ORs and 95% CIs were respectively 2.07 (1.15–3.72), 3.54 (2.06–6.09), and 3.86 (2.19–6.79) when compared with those aged 30–39. Our analysis of cases with and without T2D based on rs738409 genotypes (Table 3) indicated that diabetic individuals with the G allele had a higher cirrhosis risk compared to those with the C allele (OR, 1.71; 95% CI, 1.03–2.84), as did diabetic individuals with the GG genotype compared to those carrying the CC genotype (OR, 2.92; 95% CI, 1.07–7.99).
TABLE 3. Association of liver cirrhosis with rs738409 genotypes in individuals with and without T2D.
Discussion
Using the data of TWB participants collected from 2008 to 2015, we found that the PNPLA3 gene variant rs738409 was not associated with liver cirrhosis while type 2 diabetes was. Despite these, we observed a significant interaction between the variant and T2D as shown in the additive, genotype, and dominant models. Furthermore, our analyses showed that rs738409 GG individuals who had T2D were more prone to liver cirrhosis than their CC counterparts.
Taiwan still faces a substantial number of liver diseases, including cirrhosis. It should be noted that cirrhosis is more prevalent among people with viral hepatitis (Mokdad et al., 2014). In contrast, 12.3–57% of patients with cirrhosis were also diagnosed with T2D (Trombetta et al., 2005). In a recent review, Lee and colleagues claimed that there is a higher prevalence of diabetes in patients with liver cirrhosis than previously believed (Lee et al., 2019). Numerous studies have already suggested a link between liver cirrhosis and T2D. However, there is little data about genetic links between these two health conditions. To gain more insight into the relationship, we replicated the rs738409 variant and showed that it may interact with T2D to increase cirrhosis risk. This variant showed associations with liver disease, T2D, and other disorders in diverse populations as mentioned above. Additionally, it has been linked to non-viral hepatitis-associated hepatocellular carcinoma among T2D patients (Ueyama et al., 2016).
The relationship between viral hepatitis and liver cirrhosis is well known. Our univariate and multivariate models also showed that participants with HBV and HCV infection had an increased risk for cirrhosis. Cirrhosis is believed to develop in most viral hepatitis patients around the age of 65 or so (Pradat et al., 2007). Our study found that cirrhosis prevalence was higher among individuals aged 50 to 59 and 60 to 70.
Our study has the advantage of a large sample size. By using personal identification numbers, we were able to link genetic, clinical, and lifestyle information from the TWB database to medical records from the NHIRD; these data enabled more detailed and accurate analyses at the individual level. Our study was limited, however, by the lack of information regarding the severity of liver cirrhosis. There is a need for further research into the role of T2D and PNPLA3 rs738409 in the development of liver cirrhosis.
This study indicates that Taiwanese adults with T2D who have the PNPLA3 rs738409GG genotype might have a higher risk for liver cirrhosis than those with the CC genotype. Findings like these contribute to a better understanding of the genetic bases of liver cirrhosis, which may be useful in personalized medicine.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding authors.
Ethics Statement
The study protocol followed the guidelines for human studies and was conducted ethically in accordance with the World Medical Association Declaration of Helsinki. Ethical approval was obtained from the Institutional Review Board of Chung Shan Medical University. Participants provided written informed consent before being assessed.
Author Contributions
K-CH, ON, S-YH, S-FY, and Y-PL conceived and designed the study. S-YH, ON, and Y-PL analyzed and interpreted data. K-CH drafted the manuscript with support from S-FY. ON, S-YF, and Y-PL critically reviewed the manuscript for important intellectual content. All authors approved the final version of the manuscript.
Funding
This work was supported by Chung Shan Medical University and Tungs’ Taichung Metro Habor Hospital (grant number, CSMU-TTM-107-01) and the Ministry of Science and Technology of Taiwan (MOST 109-2121-M-040-002 and MOST 110-2121-M-040-002).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
Asrani, S. K., Devarbhavi, H., Eaton, J., and Kamath, P. S. (2019). Burden of Liver Diseases in the World. J. Hepatol. 70 (1), 151–171. doi:10.1016/j.jhep.2018.09.014
Basyte-Bacevice, V., Skieceviciene, J., Valantiene, I., Sumskiene, J., Petrenkiene, V., Kondrackiene, J., et al. (2019). TM6SF2 and MBOAT7 Gene Variants in Liver Fibrosis and Cirrhosis. Ijms 20 (6), 1277. doi:10.3390/ijms20061277
Chen, V. L., Chen, Y., Du, X., Handelman, S. K., and Speliotes, E. K. (2020). Genetic Variants that Associate with Cirrhosis Have Pleiotropic Effects on Human Traits. Liver Int. 40 (2), 405–415. doi:10.1111/liv.14321
Goh, G. B.-B., Pan, A., Chow, W.-C., Yuan, J.-M., and Koh, W.-P. (2017). Association between Diabetes Mellitus and Cirrhosis Mortality: the Singapore Chinese Health Study. Liver Int. 37 (2), 251–258. doi:10.1111/liv.13241
Grimaudo, S., Pipitone, R. M., Pennisi, G., Celsa, C., Cammà, C., Di Marco, V., et al. (2020). Association between PNPLA3 Rs738409 C>G Variant and Liver-Related Outcomes in Patients with Nonalcoholic Fatty Liver Disease. Clin. Gastroenterol. Hepatol. 18 (4), 935–944. doi:10.1016/j.cgh.2019.08.011
Holstein, A., Hinze, S., Thiessen, E., Plaschke, A., and Egberts, E.-H. (2002). Clinical Implications of Hepatogenous Diabetes in Liver Cirrhosis. J. Gastroenterol. Hepatol. 17 (6), 677–681. doi:10.1046/j.1440-1746.2002.02755.x
Lee, W. G., Wells, C. I., McCall, J. L., Murphy, R., and Plank, L. D. (2019). Prevalence of Diabetes in Liver Cirrhosis: A Systematic Review and Meta-Analysis. Diabetes Metab. Res. Rev. 35 (6), e3157. doi:10.1002/dmrr.3157
Mokdad, A. A., Lopez, A. D., Shahraz, S., Lozano, R., Mokdad, A. H., Stanaway, J., et al. (2014). Liver Cirrhosis Mortality in 187 Countries between 1980 and 2010: a Systematic Analysis. BMC Med. 12 (1), 145. doi:10.1186/s12916-014-0145-y
Nalbantoglu, I., and Jain, D. (2019). “Cryptogenic Cirrhosis: Old and New Perspectives in the Era of Molecular and Genomic Medicine,” in Seminars in Diagnostic Pathology (St Louis: Elsevier) 36 (6), 389–394.
Nfor, O. N., Wu, M.-F., Lee, C.-T., Wang, L., Liu, W.-H., Tantoh, D. M., et al. (2018). Body Mass index Modulates the Association between CDKAL1 Rs10946398 Variant and Type 2 Diabetes Among Taiwanese Women. Sci. Rep. 8 (1), 13235. doi:10.1038/s41598-018-31415-4
Patel, N., Sharma, B., and Samant, H. (2019). “Cryptogenic Cirrhosis,” in: StatPearls. Treasure Island (FL): StatPearls Publishing.
Pradat, P., Voirin, N., Tillmann, H. L., Chevallier, M., and Trépo, C. (2007). Progression to Cirrhosis in Hepatitis C Patients: an Age-dependent Process. Liver Int. 27 (3), 335–339. doi:10.1111/j.1478-3231.2006.01430.x
Quintana, J. O., García-Compean, D., González, J. A., Pérez, J. Z., González, F. J., Espinosa, L. E., et al. (2016). The Impact of Diabetes Mellitus in Mortality of Patients with Compensated Liver Cirrhosis-A Prospective Study. Ann. Hepatol. 10 (1), 56–62. doi:10.1016/S1665-2681(19)31588-1
Romeo, S., Sentinelli, F., Dash, S., Yeo, G. S. H., Savage, D. B., Leonetti, F., et al. (2010). Morbid Obesity Exposes the Association between PNPLA3 I148M (Rs738409) and Indices of Hepatic Injury in Individuals of European Descent. Int. J. Obes. 34 (1), 190–194. doi:10.1038/ijo.2009.216
Sepanlou, S. G., Safiri, S., Bisignano, C., Ikuta, K. S., Merat, S., Saberifiroozi, M., et al. (2020). The Global, Regional, and National burden of Cirrhosis by Cause in 195 Countries and Territories, 1990-2017: a Systematic Analysis for the Global Burden of Disease Study 2017. Lancet Gastroenterol. Hepatol. 5 (3), 245–266. doi:10.1016/S2468-1253(19)30349-8
Smith, A., Baumgartner, K., and Bositis, C. (2019). Cirrhosis: Diagnosis and Management. Am. Fam. Physician 100 (12), 759–770.
Speliotes, E. K., Butler, J. L., Palmer, C. D., Voight, B. F., Hirschhorn, J. N., Crn, N., et al. (2010). PNPLA3variants Specifically Confer Increased Risk for Histologic Nonalcoholic Fatty Liver Disease but Not Metabolic Disease. Hepatology 52 (3), 904–912. doi:10.1002/hep.23768
Stickel, F., Lutz, P., Buch, S., Nischalke, H. D., Silva, I., Rausch, V., et al. (2020). Genetic Variation in HSD17B13 Reduces the Risk of Developing Cirrhosis and Hepatocellular Carcinoma in Alcohol Misusers. Hepatology 72 (1), 88–102. doi:10.1002/hep.30996
Tolman, K. G., Fonseca, V., Dalpiaz, A., and Tan, M. H. (2007). Spectrum of Liver Disease in Type 2 Diabetes and Management of Patients with Diabetes and Liver Disease. Diabetes care 30 (3), 734–743. doi:10.2337/dc06-1539
Trombetta, M., Spiazzi, G., Zoppini, G., and Muggeo, M. (2005). Review Article: Type 2 Diabetes and Chronic Liver Disease in the Verona Diabetes Study. Aliment. Pharmacol. Ther. 22, 24–27. doi:10.1111/j.1365-2036.2005.02590.x
Ueyama, M., Nishida, N., Korenaga, M., Korenaga, K., Kumagai, E., Yanai, H., et al. (2016). The Impact of PNPLA3 and JAZF1 on Hepatocellular Carcinoma in Non-viral Hepatitis Patients with Type 2 Diabetes Mellitus. J. Gastroenterol. 51 (4), 370–379. doi:10.1007/s00535-015-1116-6
Xia, M.-F., Lin, H.-D., Chen, L.-Y., Wu, L., Ma, H., Li, Q., et al. (2019). The PNPLA3 Rs738409 C>G Variant Interacts with Changes in Body Weight over Time to Aggravate Liver Steatosis, but Reduces the Risk of Incident Type 2 Diabetes. Diabetologia 62 (4), 644–654. doi:10.1007/s00125-018-4805-x
Yang, T.-W., Wang, C.-C., Tsai, M.-C., Wang, Y.-T., Tseng, M.-H., and Lin, C.-C. (2020). Comorbidities and Outcome of Alcoholic and Non-alcoholic Liver Cirrhosis in Taiwan: A Population-Based Study. Ijerph 17 (8), 2825. doi:10.3390/ijerph17082825
Keywords: liver, observational study, type 2 diabetes, polymorphism, chronic disease BUN (mg/dl) 0.996 0.954–1.040
Citation: Hsueh K-C, Nfor ON, Hsu S-Y, Yang S-F and Liaw Y-P (2022) Type 2 Diabetes, PNPLA3 rs738409 Polymorphism, and the Risk of Liver Cirrhosis: Analysis of Taiwan Biobank. Front. Genet. 13:822700. doi: 10.3389/fgene.2022.822700
Received: 30 November 2021; Accepted: 22 February 2022;
Published: 07 March 2022.
Edited by:
Robert Klein, Icahn School of Medicine at Mount Sinai, United StatesReviewed by:
Umberto Vespasiani-Gentilucci, Campus Bio-Medico University, ItalyLanlan Chen, Jilin University, China
Copyright © 2022 Hsueh, Nfor, Hsu, Yang and Liaw. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shun-Fa Yang, ysf@csmu.edu.tw; Yung-Po Liaw, Liawyp@csmu.edu.tw