Skip to main content

PNPLA3 rs738409 risk genotype decouples TyG index from HOMA2-IR and intrahepatic lipid content

Abstract

Background

Recent reports suggested a different predictive value for TyG index compared to HOMA-IR in coronary artery calcification (CAC) and other atherosclerotic outcomes, despite that both indices are proposed as surrogate markers of insulin resistance. We hypothesized a key role for liver pathology as an explanation and therefore assessed the relationship among the two indices and the intrahepatic lipid content stratified by PNPLA3 rs738409 genotypes as a known non-alcoholic fatty liver disease (NAFLD) genetic risk.

Methods

Thirty-nine women from a prior GDM-genetic study were recalled with PNPLA3 rs738409 CC and GG genotypes for metabolic phenotyping and to assess hepatic triglyceride content (HTGC). 75 g OGTT was performed, fasting lipid, glucose, insulin levels and calculated insulin resistance indices (TyG and HOMA2-IR) were used. HTGC was measured by MR based methods. Mann–Whitney-U, χ2 and for the correlation analysis Spearman rank order tests were applied.

Results

The PNPLA3 rs738409 genotype had a significant effect on the direct correlation between the HOMA2-IR and TyG index: the correlation (R = 0.52, p = 0.0054) found in the CC group was completely abolished in those with the GG (NAFLD) risk genotype. In addition, the HOMA2-IR correlated with HTGC in the entire study population (R = 0.69, p < 0.0001) and also separately in both genotypes (CC R = 0.62, p = 0.0006, GG: R = 0.74, p = 0.0058). In contrast, the correlation between TyG index and HTGC was only significant in rs738409 CC genotype group (R = 0.42, p = 0.0284) but not in GG group. A similar pattern was observed in the correlation between TG and HTGC (CC: R = 0.41, p = 0.0335), when the components of the TyG index were separately assessed.

Conclusions

PNPLA3 rs738409 risk genotype completely decoupled the direct correlation between two surrogate markers of insulin resistance: TyG and HOMA2-IR confirming our hypothesis. The liver lipid content increased in parallel with the HOMA2-IR independent of genotype, in contrast to the TyG index where the risk genotype abolished the correlation. This phenomenon seems to be related to the nature of hepatic fat accumulation and to the different concepts establishing the two insulin resistance markers.

Introduction

Insulin resistance (IR) refers to a decreased sensitivity of peripheral tissues to insulin and is a hallmark pathophysiologic feature of type 2 diabetes mellitus (T2DM) that clinically develops when the β-cell compensatory mechanisms might no longer overcome the increased insulin need due to the peripheral IR [1,2,3]. IR is closely related to the sedentary lifestyle, obesity and non-alcoholic fatty liver disease (NAFLD) that are highly prevalent [4, 5].

The most established method to assess the IR is the euglycemic–hyperinsulinemic clamp [6] developed as early as 1979, however its cost and complexity limited its everyday clinical use. There are many IR indices that are easier to assess in the clinical practice out of which HOMA-IR index appears to be the most widely used since its first description in 1985 till date [7, 8]. HOMA-IR was developed as a breakthrough model of the glucose-insulin feedback system in the homeostatic (overnight fasting) state and hepatic insulin resistance under these conditions is the major determinant of HOMA-IR [9, 10].

Nevertheless there are other surrogate markers of IR such as the triglyceride-glucose (TyG) index [11, 12] that are increasingly used due to the clinicians need to meaningfully assess IR without measuring the insulin levels. TyG index was produced in a logarithmized manner from the fasting triglycerides and glucose values and was reported to be a predictor of adverse cardiovascular outcomes in acute coronary syndrome patients with, but also without T2DM [13, 14].

Unsurprisingly both the HOMA-IR [15, 16] and the TyG index [17] are closely related to NAFLD. NAFLD is not only the most common chronic liver disease currently, but also its prevalence is doubled in patients with T2DM [4] and due to the increased IR, decreased insulin clearance [18,19,20] and the key metabolic functions of the liver the presence of NAFLD is a crucial factor in the speed-up of the metabolic deterioration in T2DM development and worsening of atherogenic dyslipidaemia [5, 21, 22].

Although both the HOMA-IR and TyG indices are proposed as surrogate markers of IR, recently the TyG index was found to be superior compared to the HOMA-IR in predicting coronary artery calcification (CAC) [23].

To the best of our current knowledge, the PNPLA3 rs738409 variant shows the strongest association with NAFLD development [24]. The rs738409 gene variant and particularly the GG genotype is also associated with the progression of NAFLD throughout its entire spectrum [15, 25, 26], including complications occurring in patients with more advanced liver disease [27, 28].

The risk G allele occurs with high enough allele frequency in many of the populations worldwide, including the European (RAF: 23% [29]) and Hungarian subjects (RAF: 22%—[30]) to have a potential influence at population level.

Therefore, we have assessed the relationship between hepatic fat content, the rs738409 gene variant and the potential genotype effect on the correlation between the HOMA-IR and TyG indeces in a T2DM prone middle aged Hungarian female population.

Patients and methods

Participants and study design

We have access to data of over 600 Hungarian female subjects genotyped for 77 gene variants including PNPLA3 rs738409 in a prior study [30]. A genetic based recall study (GBR) design was applied due to that GBR studies allow the assessment of genotype–phenotype associations in studies with substantially smaller samples on the basis of pre-existing genetic data resulting in a substantially higher power compared to same-sized conventional studies [31,32,33]. Thirty-nine women were recruited with known PNPLA3 rs738409 genotype data: only individuals with CC or GG genotype were eligible for this study-part. A patients flow-chart is indicated on Fig. 1.

Fig. 1
figure 1

Flow chart of the genotype-based recall (GBR) study. *reference number [30]

The study was conducted according to the Declaration of Helsinki after receiving approval from the relevant institutional bodies (Semmelweis University, Regional and Institutional Committee of Science, the Medical Research Council, Scientific and Research Committee of Hungary, 14486-6/2017/EKU). Participants had previously given informed consent for the whole project [30].

Exclusion criteria

Exclusion criteria from the original case-control study [30] were extended to returning volunteers to the GBR study and included: prediabetes/T2DM diagnosed since delivery, claustrophobia and/or MR-incompatible metal implants, weight over 200 kg, waist/hip circumference ≥ 135 cm/148 cm, taking medication with a known impact on glycemic and lipid traits or HTGC, significant alcohol consumption (> 20 g/day), malignancy, or any other causes of fatty liver as detailed in the new MAFLD definition consensus statement [34]. In addition, patients with other major chronic/acute diseases or ongoing pregnancy/breastfeeding were also excluded.

Phenotyping

Anthropometric data (age, weight, height, BMI) was recorded. 75 g OGTT was performed. Insulin and glucose levels were measured in addition to fasting „routine” laboratory parameters, including HbA1c, lipid profile and liver enzymes. From fasting trigliceride (TG), glucose and insulin levels insulin resistance indices: Trigliceride-Glucose index: TyG = ln[38,67*Tg(mg/dl)*18*glucose(mg/dl)/2], and homeostasis model assessement 2 insulin resistance index (HOMA2-IR) values were calculated [35, 36]. Diabetes mellitus (DM) and prediabetes (impaired fasting glucose (IFG)/impaired glucose tolerance (IGT)) were diagnosed according to the WHO, 2016 American Diabetes Association (ADA), and Hungarian guidelines [37,38,39] based on the 75 g OGTT and/or HbA1c results.

HTGC measurement

All MR imaging (MRI) sessions were acquired on a clinical 3 T MRI system (Prisma, Siemens Healthineers, Erlangen, Germany) with the subject in a supine position. For MRI protocol, a standard body-array of 18 channel flexible coils was positioned on the liver region and combined with a spine array coil located below the subject. Imaging proton density fat-fraction (PDFF), unenhanced axial images were obtained by using a low–flip-angle, six-echo two-dimensional spoiled gradient-recalled-echo sequence with all array coil elements (TE = 2, 4.1, 6.2, 8.2, 10.2, and 12.3 ms). The repetition time and flip angle were chosen to avoid T1 weighting: TR: 15 ms, Flip-angle: 11°, FOV = 240 × 400 mm in plane, matrix = 240 × 130, slice thickness = 3.5 mm, space between slices 4.3 mm. Phase and magnitude images were systematically saved. For MR imaging using custom-made MATLAB routines, multi-section liver PDFF maps were generated offline from the source images via joint estimation of water and fat images and field maps [40]. This method uses a complex signal model similar to advanced multipoint DIXON/IDEAL algorithms (incorporating a multipeak/multifrequency fat model and T2* decay) analyzing six-echo FLASH complex images. Voxels representing liver tissue were delineated with a freehand ROI defining tool by an expert radiologist on multi-section PDFF maps excluding vessels, liver and edges, and artifacts. In each subject the mean value of the selected voxels was calculated, representing the average liver fat fraction. Fatty liver was diagnosed if the liver fat fraction exceeded 5.5% [15].

Statistics

Due to the limited sample size and the non-normal distribution of the majority of variables we used non-parametric tests. Data are expressed in median and 25th-75th percentiles. For evaluation of the association between two parameters Spearman’s correlation test and for comparing two (genotype) groups Mann–Whitney U/χ2 tests were applied. TIBCO Statistica (version 13.4.0.14, TIBCO Software Inc.) software was used.

Results

Study population clinical characteristics

Age and BMI values of study participants were 37.0 (34–40) years and 26.2 (22.8–32.5) kg/m2, respectively. We diagnosed thirteen individuals with prediabetes and one with overt T2DM. No patient had established cardiovascular disease at enrolment. No significant difference of anthropometric data, glycaemic or lipid values was detected between genotype groups, except of the HTGC values (and NAFLD prevalence). The characteristic of the study population is indicated in Table 1.

Table 1 Population characteristic

Correlation between TyG and HOMA2-IR indices stratified by PNPLA3 rs738409 genotype

Between the two insulin resistance indices we found an association (R = 0.42, p = 0.0078), however a significant correlation was only detectable in individuals with the rs738409 CC genotype (R = 0.52, p = 0.0054) but not with GG (r = − 0.16, p = 0.62). The scatterplot is indicated in Fig. 2.

Fig. 2
figure 2

Correlation between HOMA2-IR and TyG indices stratified by PNPLA3 rs738409 genotypes. Spearman rank order test results indicated in the figure

Genotype effect on correlations among insulin resistance indices, their components and HTGC

HTGC values significantly correlated with both markers of insulin resistance: (TyG R = 0.36, p = 0.0252, HOMA2-IR: R = 0.69, p < 0.0001), and their components (fasting glucose: R = 0.51, p = 0.0008, fasting insulin: R = 0.66, p < 0.0001, fasting TG: R = 0.32, p = 0.0445).

After genotype stratification In contrast to HOMA2-IR and its components, where the correlations were nearly parallel in the two genotype groups after genotype stratification, the correlations of TyG index and fasting TG were only significant in the rs738409 CC genotype group indicating different associations between the parameters according to the genotypes to HOMA2-IR and its components, where the correlations were nearly parallel in the two genotype groups, the correlations of TyG index and fasting TG were only significant in the rs738409 CC genotype group indicating different associations between the parameters according to the genotypes. Similarly, the correlation between HOMA2-IR and fasting TG was disrupted in individuals with GG genotype. Correlations between HTGC and TyG, HTGC and HOMA2-IR, HTGC and TG, HOMA2-IR and TG are indicated in Fig. 3A, B, C, D, respectively.

Fig. 3
figure 3

Correlations among HTGC, insulin resistance indices and TG values stratified by PNPLA3 rs738409 genotypes (CC vs GG): A: between HTGC and TyG index, B: between HTGC and HOMA2-IR, C: between HTGC and TG, D: between HOMA2-IR and TG. Spearman rank order test results indicated in the figure. HTGC: hepatic triacylglycerol content, TG: fasting serum triglyceride

Discussion

We assessed intrahepatic lipid content measured by MR method and two indices proposed as surrogate markers of insulin resistance indices (TyG and HOMA2-IR) in 39 middle-aged women prone to T2DM development. Genotype based recall study design was applied on the basis of PNPLA3 NAFLD risk genotype to decrease the number of patients needed to be enrolled [31, 33]. Therefore, only patients with rs738409 homozygous genotypes (GG vs CC) were included. Despite the two IR indices were suggested to have a linear relationship [17, 41], we found compelling difference in the relationship between TyG and HOMA2-IR when the data were stratified by the PNPLA3 genotypes. The profound genetic effect on the correlations between TyG, TG and HTGC and also between TG and HOMA2-IR was observed.

Insulin resistance and hyperinsulinemia are classified as major factors contributing to pathology of T2DM development and its cardiovascular complications. Many concepts and methods were developed to measure the “insulin resistance” (IR), however no single parameter could meet the need of a precise pathophysiological assessment combined with an easy use for the everyday clinics. Out of the many IR indices HOMA-IR index appears to be the most widely used in the clinical practice, including the recent recommendation to establish the diagnosis of metabolic (dysfunction) associated fatty liver disease (MAFLD) [34]. The HOMA-IR and the further developed HOMA2-IR are proposed to be mainly associated with hepatic insulin resistance under fasting steady-state conditions [8, 36, 42]. Despite HOMA-IR is easily calculated from only two parameters there is a need from clinicians to find alternative surrogate markers of IR which do not require the measurement of fasting insulin levels such as the triglyceride-glucose (TyG) index [12].

Two markers of insulin resistance, the HOMA-IR and TyG were compared in a recent study on the prognostic value of T2DM development and atherosclerosis and the authors concluded that TyG index performed better in both endpoints [43].

Nevertheless there is still no clear-cut opinion in the literature regarding the predictive powers of the TyG and HOMA-IR index in CVD, and studies suggested only a moderate coronary artery disease (CAD) predictive power for TyG, although its inclusion into predictive models of MACEs could improve the predictive accuracy in patients with ACS [44].

Due to that hepatic insulin resistance is significantly increased in NAFLD that is the most prevalent liver disease and strongly related to obesity, T2DM, dyslipidaemia and metabolic syndrome, the pathophysiology of the liver lipid accumulation could be a key factor in (hepatic) IR associated to the common metabolic diseases. Therefore we assessed the effect of PNPLA3 rs738409 gene variant on the correlation between TyG and HOMA2-IR. This gene variant was identified in fatty liver genome-wide association studies (GWAS) [24] with pathogenic role throughout the entire NAFLD spectrum [25, 27, 28].

We astonishingly found that TyG and HOMA2-IR indices were completely dissociated in individuals with the rs738409 GG risk genotype (M148M) in contrast to those with the CC (I148I) genotype where the strong positive correlation was conserved.

This finding should mean that despite both the TyG and the HOMA2-IR are considered as surrogate markers of insulin resistance they display fundamentally different characteristics according to the genetic background of the liver lipid accumulation confirming the hypothesis of the PNPLA3 rs738409 genotype effect.

Although prior studies with human liver tissue samples reported that the hepatic content of diacylglycerol (DAG) species (FA18:1, implicated in diminished insulin signaling) remained unaltered in rs738409 homozygous and heterozygous G allele carriers (M148M and I148M) and it was hypothetically suggested that the G risk allele carriers might be protected from insulin resistance [45] we report here that the HOMA2-IR values remained correlated with HTGC. In addition, our observation is consistent is consistent with prior results where HOMA-IR/HOMA2-IR and HTGC values were directly measured and reported to be correlated in subjects with PNPLA3 risk genotypes [16]. In contrast to the HOMA2-IR only the TyG index displayed the theoretically expected characteristic in correlation with HTGC in individuals homozygous for the PNPLA3 NAFLD risk genotype based on the reported hepatic DAG species content critical in disturbing the insulin signal [45].

We observed that the correlation between TyG and HTGC was diminished in those with the rs738409 GG risk genotype. The explanation could be attributed to that the input parameters and their PNPLA3 risk genotype associated changes are largely different. The insulin concentrations were increasing in parallel with HTGC, in contrast to the glucose levels that were not different. Consistently, the fasting insulin levels and the HOMA-IR indices were significantly higher and the glucose levels were unchanged in a severely obese large European cohort in the rs738409 GG genotype group [46] that was presumably presented with higher HTGC values as well. Opposingly, the serum TG levels in the morbid obese patients with the GG genotype was found to be lower compared to CC genotype group [46] and confirmed subsequently [47, 48]. The lower circulating TG levels in the PNPLA3 risk genotype group could be due to the lower hepatic TG efflux in the fasting state associated to the genetic effect on intrahepatocellular lipolysis [46, 48, 49]. Based on that the hepatic VLDL output occurs in parallel with the TG efflux in the fasting state [50] these observations might well explain that the TyG index could be an useful prognostic marker for CAD, especially in high CVD risk patients with T2DM [43, 44, 51,52,53]. Furthermore, in a recent study with over 4000 participants undergoing cardiac CT the TyG index was found to be superior compared to the HOMA-IR in predicting CAC [23]. It could be outlined that the rs738409 G allele frequencies in the East Asian populations are substantially higher compared to populations with European or African origin [29] which is consistent with the results of the latter East Asian study [23] and could contribute to the difference found between TyG index and HOMA-IR in predicting CAC.

The current state of the art is that HOMA-IR primarily reflects hepatic insulin resistance [10]. The increase of HOMA2-IR occurred in parallel with the increase of HTGC in those with the PNPLA3 rs738409 GG risk genotype and partially behaves like a biomarker of intrahepatic steatosis that is dissociated from the other surrogate marker of IR, the TyG index composed of only metabolic parameters but does not directly include serum insulin level.

As a consequence, it may also be raised, that HOMA-IR could be significantly dissociated from the TyG index in those populations where the rs738409 G NAFLD risk allele is the major PNPLA3 allele [29, 54].

The authors suggest that the genetic background of the pathologic liver lipid accumulation and its relation to HOMA-IR and TyG index are important insights towards a better understanding, more precise interpretation and perhaps towards a better clinical use of these IR estimates.

Limitations

Our results were limited by the sample size and the lack of supplementary methods to study the underlying pathology. Further studies are needed to replicate our findings and to confirm our results in other populations.

Availability of data and materials

The datasets used and analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

BMI:

Body mass index

FIB-4:

Fibrosis 4 score

GBR:

Genetic-based recall (study)

GDM:

Gestational diabetes mellitus

HOMA:

Homeostasis model assessment

IFG:

Impaired fasting glucose

HTGC:

Hepatic triglyceride content

IGT:

Impaired glucose tolerance

IR:

Insulin resistance

MR:

Magnetic resonance

NAFLD:

Non-alcoholic fatty liver disease

NFS:

NAFLD fibrosis score

OGTT:

Oral glucose tolerance test

PDFF:

Proton density fat-fraction

PNPLA3 :

Patatin-like phospholipase domain-containing protein 3

TyG index:

Triglyceride-glucose index

T2DM:

Type 2 diabetes mellitus

TG:

Triacylglycerol

W/H:

Waist-to-hip ratio

References

  1. Reaven GM. Role of insulin resistance in human disease. Diabetes. 1988;37:1595–607.

    Article  CAS  PubMed  Google Scholar 

  2. Petersen MC, Shulman GI. Mechanisms of insulin action and insulin resistance. Physiol Rev Am Physiol Soc. 2018;98:2133–223.

    CAS  Google Scholar 

  3. Taylor R. Insulin resistance and type 2 diabetes. Diabetes. 2012;61:778–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Younossi ZM, Golabi P, de Avila L, Paik JM, Srishord M, Fukui N, et al. The global epidemiology of NAFLD and NASH in patients with type 2 diabetes: a systematic review and meta-analysis. J Hepatol. 2019;71:793–801.

    Article  PubMed  Google Scholar 

  5. Firneisz G. Non-alcoholic fatty liver disease and type 2 diabetes mellitus: the liver disease of our age? World J Gastroenterol. 2014;20:9072–89.

    PubMed  PubMed Central  Google Scholar 

  6. DeFronzo RA, Tobin JD, Andres R. Glucose clamp technique: a method for quantifying insulin secretion and resistance. American journal of physiology-endocrinology and metabolism. Am Physiol Soc. 1979;237:214.

    Google Scholar 

  7. Park SY, Gautier J-F, Chon S. Assessment of insulin secretion and insulin resistance in human. Diabetes Metab J Korean Diabetes Assoc. 2021;45:641–54.

    Article  Google Scholar 

  8. Matthews DR, Hosker JP, Rudenski AS, Naylor BA, Treacher DF, Turner RC. Homeostasis model assessment: insulin resistance and β-cell function from fasting plasma glucose and insulin concentrations in man. Diabetologia. 1985;28:412–9.

    Article  CAS  PubMed  Google Scholar 

  9. Tripathy D, Almgren P, Tuomi T, Groop L. Contribution of insulin-stimulated glucose uptake and basal hepatic insulin sensitivity to surrogate measures of insulin sensitivity. Diabetes Care. 2004;27:2204–10.

    Article  CAS  PubMed  Google Scholar 

  10. Retnakaran R, Qi Y, Ye C, Hanley AJG, Connelly PW, Sermer M, et al. Hepatic insulin resistance is an early determinant of declining β-cell function in the first year postpartum after glucose intolerance in pregnancy. Diabetes Care. 2011;34:2431–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Simental-Mendía LE, Rodríguez-Morán M, Guerrero-Romero F. The product of fasting glucose and triglycerides as surrogate for identifying insulin resistance in apparently healthy subjects. Metab Syndr Relat Disord. 2008;6:299–304.

    Article  PubMed  Google Scholar 

  12. Guerrero-Romero F, Simental-Mendía LE, González-Ortiz M, Martínez-Abundis E, Ramos-Zavala MG, Hernández-González SO, et al. The product of triglycerides and glucose, a simple measure of insulin sensitivity comparison with the euglycemic-hyperinsulinemic clamp. J Clin Endocrinol Metab. 2010;95:3347–51.

    Article  CAS  PubMed  Google Scholar 

  13. Ma X, Dong L, Shao Q, Cheng Y, Lv S, Sun Y, et al. Triglyceride glucose index for predicting cardiovascular outcomes after percutaneous coronary intervention in patients with type 2 diabetes mellitus and acute coronary syndrome. Cardiovasc Diabetol. 2020;19:31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Zhao Q, Zhang T-Y, Cheng Y-J, Ma Y, Xu Y-K, Yang J-Q, et al. Triglyceride-glucose index as a surrogate marker of insulin resistance for predicting cardiovascular outcomes in nondiabetic patients with non-st-segment elevation acute coronary syndrome undergoing percutaneous coronary intervention. J Atheroscler Thromb. 2021;28:1175–94.

    Article  CAS  PubMed  Google Scholar 

  15. European Association for the Study of the Liver (EASL) European Association for the Study of Diabetes (EASD) European Association for the Study of Obesity (EASO). EASL–EASD–EASO clinical practice guidelines for the management of non-alcoholic fatty liver disease. J Hepatol. 2016;64:1388–402.

    Article  Google Scholar 

  16. Isokuortti E, Zhou Y, Peltonen M, Bugianesi E, Clement K, Bonnefont-Rousselot D, et al. Use of HOMA-IR to diagnose non-alcoholic fatty liver disease: a population-based and inter-laboratory study. Diabetologia. 2017;60:1873–82.

    Article  PubMed  Google Scholar 

  17. Lee SB, Kim MK, Kang S, Park K, Kim JH, Baik SJ, et al. Triglyceride glucose index is superior to the homeostasis model assessment of insulin resistance for predicting Nonalcoholic fatty liver disease in Korean adults. Endocrinol Metab Korean Endocrine Soc. 2019;34:179–86.

    Article  Google Scholar 

  18. Kotronen A, Vehkavaara S, Seppälä-Lindroos A, Bergholm R, Yki-Järvinen H. Effect of liver fat on insulin clearance. Am J Physiol Endocrinol Metab. 2007;293:E1709–15.

    Article  CAS  PubMed  Google Scholar 

  19. Kotronen A, Juurinen L, Tiikkainen M, Vehkavaara S, Yki-Järvinen H. Increased liver fat, impaired insulin clearance, and hepatic and adipose tissue insulin resistance in type 2 diabetes. Gastroenterology. 2008;135:122–30.

    Article  CAS  PubMed  Google Scholar 

  20. Trico D, Galderisi A, Galuppo B, Pierpont B, Samuels S, Santoro N, et al. Fatty liver, irrespective of ethnicity, is associated with reduced insulin clearance and insulin resistance in obese youths. Endocrine Abstracts. Bioscientifica; 2020. https://www.endocrine-abstracts.org/ea/0070/ea0070aep256. Accessed 8 Jun 2022

  21. Targher G, Corey KE, Byrne CD, Roden M. The complex link between NAFLD and type 2 diabetes mellitus—mechanisms and treatments nat rev gastroenterol hepatol. Nature Publ Group. 2021;18:599–612.

    Google Scholar 

  22. Watt MJ, Miotto PM, De Nardo W, Montgomery MK. The liver as an endocrine organ—linking NAFLD and insulin resistance. Endocr Rev. 2019;40:1367–93.

    Article  PubMed  Google Scholar 

  23. Kim MK, Ahn CW, Kang S, Nam JS, Kim KR, Park JS. Relationship between the triglyceride glucose index and coronary artery calcification in Korean adults. Cardiovasc Diabetol. 2017;16:108.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Romeo S, Kozlitina J, Xing C, Pertsemlidis A, Cox D, Pennacchio LA, et al. Genetic variation in PNPLA3 confers susceptibility to nonalcoholic fatty liver disease. Nat Genet. 2008;40:1461–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. Sookoian S, Pirola CJ. Meta-analysis of the influence of I148M variant of patatin-like phospholipase domain containing 3 gene (PNPLA3) on the susceptibility and histological severity of nonalcoholic fatty liver disease. Hepatology. 2011;53:1883–94.

    Article  CAS  PubMed  Google Scholar 

  26. Kotronen A, Johansson LE, Johansson LM, Roos C, Westerbacka J, Hamsten A, et al. A common variant in PNPLA3, which encodes adiponutrin, is associated with liver fat content in humans. Diabetologia. 2009;52:1056–60.

    Article  CAS  PubMed  Google Scholar 

  27. Mandorfer M, Scheiner B, Stättermayer AF, Schwabl P, Paternostro R, Bauer D, et al. Impact of patatin-like phospholipase domain containing 3 rs738409 G/G genotype on hepatic decompensation and mortality in patients with portal hypertension. Aliment Pharmacol Ther. 2018;48:451–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Hassan MM, Kaseb A, Etzel CJ, El-Serag H, Spitz MR, Chang P, et al. Genetic variation in the PNPLA3 gene and hepatocellular carcinoma in USA: risk and prognosis prediction. Mol Carcinog. 2013;52:139–47.

    Article  Google Scholar 

  29. Cunningham F, Allen JE, Allen J, Alvarez-Jarreta J, Amode MR, Armean IM, et al. Ensembl 2022. Nucleic Acids Res. 2022;50:D988–95.

    Article  CAS  PubMed  Google Scholar 

  30. Rosta K, Al-Aissa Z, Hadarits O, Harreiter J, Nádasdi Á, Kelemen F, et al. Association study with 77 SNPs confirms the robust role for the rs10830963/G of MTNR1B variant and identifies two novel associations in gestational diabetes mellitus development. PLoS ONE. 2017;12:e0169781.

    Article  PubMed  PubMed Central  Google Scholar 

  31. Franks PW, Timpson NJ. Genotype-based recall studies in complex cardiometabolic traits. Circ Genom Precis Med. 2018. https://doi.org/10.1161/CIRCGEN.118.001947.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Corbin LJ, Tan VY, Hughes DA, Wade KH, Paul DS, Tansey KE, et al. Formalising recall by genotype as an efficient approach to detailed phenotyping and causal inference. Nat Commun. 2018;9:711.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Atabaki-Pasdar N, Ohlsson M, Shungin D, Kurbasic A, Ingelsson E, Pearson ER, et al. Statistical power considerations in genotype-based recall randomized controlled trials. Sci Rep. 2016;6:37307.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Eslam M, Newsome PN, Sarin SK, Anstee QM, Targher G, Romero-Gomez M, et al. A new definition for metabolic dysfunction-associated fatty liver disease: an international expert consensus statement. J Hepatol. 2020;73:202–9.

    Article  PubMed  Google Scholar 

  35. Alizargar J, Hsieh N-C, Wu S-FV. The correct formula to calculate triglyceride-glucose index (TyG). J Pediatr Endocrinol Metab. 2020;33:945–6.

    Article  PubMed  Google Scholar 

  36. Levy JC, Matthews DR, Hermans MP. Correct homeostasis model assessment (HOMA) evaluation uses the computer program. Diabetes Care. 1998;21:2191–2.

    Article  CAS  PubMed  Google Scholar 

  37. World Health Organization & International Diabetes Federation (2006). Definition and diagnosis of diabetes mellitus and intermediate hyperglycaemia : report of a WHO/IDF consultation. World Health Organization. https://apps.who.int/iris/handle/10665/43588.

  38. American Diabetes Association. 1. Improving care and promoting health in populations: standards of medical care in diabetes—2018. Diabetes Care. 2018;41:S7-12.

    Article  Google Scholar 

  39. Gaál Z, Gerő L, Hidvégi T, Jermendy G, Kempler P, Winkler G, et al. Clinical practice guideline—diagnosis, antihyperglycaemic treatment and care of patients with diabetes in adulthood. Jermendy G (edr). Diabetologia Hungarica. 2017; 25:3–77.

  40. Hernando D, Kellman P, Haldar JP, Liang Z-P. Robust water/fat separation in the presence of large field inhomogeneities using a graph cut algorithm. Magn Reson Med. 2010;63:79–90.

    PubMed  PubMed Central  Google Scholar 

  41. Guerrero-Romero F, Villalobos-Molina R, Jiménez-Flores JR, Simental-Mendia LE, Méndez-Cruz R, Murguía-Romero M, et al. Fasting triglycerides and glucose index as a diagnostic test for insulin resistance in young adults. Arch Med Res. 2016;47:382–7.

    Article  CAS  PubMed  Google Scholar 

  42. Reaven GM. What do we learn from measurements of HOMA-IR? Diabetologia. 2013;56:1867–8.

    Article  CAS  PubMed  Google Scholar 

  43. Son D-H, Lee HS, Lee Y-J, Lee J-H, Han J-H. Comparison of triglyceride-glucose index and HOMA-IR for predicting prevalence and incidence of metabolic syndrome. Nutr Metab Cardiovasc Dis. 2022;32:596–604.

    Article  CAS  PubMed  Google Scholar 

  44. Tao L-C, Xu J, Wang T, Hua F, Li J-J. Triglyceride-glucose index as a marker in cardiovascular diseases: landscape and limitations. Cardiovasc Diabetol. 2022;21:68.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Franko A, Merkel D, Kovarova M, Hoene M, Jaghutriz BA, Heni M, et al. Dissociation of fatty liver and insulin resistance in I148M PNPLA3 carriers: differences in diacylglycerol (DAG) FA18:1 lipid species as a possible explanation. Nutrients. 2018;10:1314.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Palmer CNA, Maglio C, Pirazzi C, Burza MA, Adiels M, Burch L, et al. Paradoxical lower serum triglyceride levels and higher type 2 diabetes mellitus susceptibility in obese individuals with the PNPLA3 148M variant. PLoS ONE. 2012;7:e39362.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Hyysalo J, Gopalacharyulu P, Bian H, Hyötyläinen T, Leivonen M, Jaser N, et al. Circulating triacylglycerol signatures in nonalcoholic fatty liver disease associated with the I148M variant in PNPLA3 and with obesity. Diabetes. 2014;63:312–22.

    Article  CAS  PubMed  Google Scholar 

  48. Liu DJ, Peloso GM, Yu H, Butterworth AS, Wang X, Mahajan A, et al. Exome-wide association study of plasma lipids in >300,000 individuals. Nat Genet Nat Publ Group. 2017;49:1758–66.

    Article  CAS  Google Scholar 

  49. Romeo S, Savage DB. Lipase tug of war: PNPLA3 sequesters ABHD5 from ATGL. Nat Metab. 2019;1:505.

    Article  CAS  PubMed  Google Scholar 

  50. Gotto AM, Pownall HJ, Havel RJ. [1] Introduction to the plasma lipoproteins. Methods in Enzymology. Academic Press; 1986. p. 3–41. https://www.sciencedirect.com/science/article/pii/0076687986280611. Accessed 19 Nov 2022

  51. Won K-B, Park EJ, Han D, Lee JH, Choi S-Y, Chun EJ, et al. Triglyceride glucose index is an independent predictor for the progression of coronary artery calcification in the absence of heavy coronary artery calcification at baseline. Cardiovasc Diabetol. 2020;19:34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Lee EY, Yang HK, Lee J, Kang B, Yang Y, Lee S-H, et al. Triglyceride glucose index, a marker of insulin resistance, is associated with coronary artery stenosis in asymptomatic subjects with type 2 diabetes. Lipids Health Dis. 2016;15:155.

    Article  PubMed  PubMed Central  Google Scholar 

  53. Song S, Choi S, Park HE, Han H, Park SH, Sung J, et al. Incremental prognostic value of triglyceride glucose index additional to coronary artery calcium score in asymptomatic low-risk population. Cardiovasc Diabetol. 2022;21:193.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  54. Cunningham F, Allen JE, Allen J, Alvarez-Jarreta J, Amode MR, Armean IM, et al. Ensembl 2022. Nucleic Acids Res. 2022; 50(1):D988-D995

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

Open access funding provided by Semmelweis University. EFSD New Horizons Grant, Wörwag Pharma Doctoral Fellowship Award, Hungarian Diabetes Association Grant to Support the Clinical Diabetes Research Conducted in Doctoral Schools—2022.

Author information

Authors and Affiliations

Authors

Contributions

Conceptualization, GF; formal analysis, ÁN and GF; funding acquisition, GF; investigation, ÁN, VG, and GF; project administration, AS and GF; resources, TM; supervision, AS and GF; validation, GF; writing—original draft, ÁN and GF; writing—review & editing, GF. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Gábor Firneisz.

Ethics declarations

Ethics approval and consent to participate

The study was conducted according to the Declaration of Helsinki after receiving approval from the relevant institutional bodies (Semmelweis University, Regional and Institutional Committee of Science, the Medical Research Council, Scientific and Research Committee of Hungary, 14486-6/2017/EKU). Verbal and written consent was obtained from all participants.

Consent for publication

All named authors meet the International Committee of Medical Journal Editors (ICMJE) criteria for authorship for this article, take responsibility for the integrity of the work as a whole, and had no objection for this version of the manuscript to be submitted for publication.

Competing interests

The authors declare no conflicts of interest regarding this manuscript.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nádasdi, Á., Gál, V., Masszi, T. et al. PNPLA3 rs738409 risk genotype decouples TyG index from HOMA2-IR and intrahepatic lipid content. Cardiovasc Diabetol 22, 64 (2023). https://doi.org/10.1186/s12933-023-01792-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12933-023-01792-w

Keywords