ABSTRACT
Background Interleukin-6 (IL-6) plays important and dynamic roles in inflammation associated with fatty liver disease over all stages, from simple steatosis to steatohepatitis, cirrhosis and cancer. IL-6 signals locally, but also circulates with multiple co-factors that control paracrine and endocrine signaling. As inflammation is a main driver of liver fibrosis, we investigated relationships between circulating components of the interleukin-6 signaling pathway (IL-6, sIL-6R and sgp130) and liver pathology in subjects with metabolically associated fatty liver disease (MAFLD) or steatohepatitis (MASH).
Methods Predictive performances of plasma IL-6, sIL-6R and sgp130 were investigated in two independent cohorts: 1) patients with biopsy-confirmed MASH (n=49), where magnetic resonance spectroscopy (MRS), imaging (MRI) and elastography (MRE) assessed liver fat, volume and stiffness; and 2) patients with morbid obesity (n=245) undergoing bariatric surgery where histological staging of steatosis, activity, and fibrosis determined MASH severity. Correlations were evaluated between IL-6, sIL-6R and sgp130 and anthropomorphic characteristics, plasma markers of metabolic disease or liver pathology.
Results In patients with MASH, plasma IL-6 and sgp130 strongly correlated with liver stiffness, which for sgp130 was independent of age, sex, BMI, diabetes, hyperlipidemia, hypertension or history of HCC. Plasma sgp130 was the strongest predictor of liver stiffness compared to common predictors and risk scores. Plasma sIL-6R correlated with liver volume independent of age, sex, and BMI. In patients with morbid obesity, circulating sgp130 correlated with advanced liver fibrosis.
Conclusion Levels of circulating sgp130 can predict progressing MASH and may be used alone or in combination with other predictors as a non-invasive measure of liver disease severity.
INTRODUCTION
Nonalcoholic fatty liver disease (NAFLD) is rapidly increasing worldwide, with global prevalence of disease reaching approximately 25% (1–3). NAFLD is a multistep, progressive disorder beginning with simple steatosis that can evolve to non-alcoholic steatohepatitis (NASH), characterized by hepatocellular ballooning, lobular inflammation, and a range of fibrosis. Simple steatosis is relatively benign, while NASH can progress to cirrhosis and/or hepatocellular carcinoma (HCC) (4, 5). Steatohepatitis is now considered the second most common indication for liver transplantation and HCC in Canada (6).
Fatty liver disease stemming from metabolic disorders has garnered the unique sub-designation of MAFLD or MASH, due to its distinct etiology and progression. The transition from simple steatosis to steatohepatitis is an important stage of fatty liver disease correlating strongly with poor prognosis (7). Liver fibrosis is a primary predictor of morbidity and mortality (8). Each advancing stage adversely affects the health and survival of patients, yet we do not yet know or understand all the factors that promote these steps. Our ability to identify these factors is limited by the asymptomatic nature of the simple steatosis to MASH transition, preventing early diagnosis and study of the disease.
Initial diagnoses often use liver risk scores calculated from age, body mass index (BMI), diabetes status, and plasma liver enzymes (alanine and aspartate aminotransferases, ALT and AST) to predict NAFLD progression and severity. However, recent evidence suggests that BMI poorly correlates with NAFLD severity (9) and as high as 80% of subjects with NAFLD and 19% of subjects with biopsy-proven NASH have normal liver enzyme levels (10–13). Thus, relying on obesity and increased liver enzymes to screen for NAFLD can lead to underdiagnosis and/or underestimation of liver damage. Liver inflammation plays a key role in the transition to steatohepatitis, yet currently, there are no non-invasive, sensitive indicators of the early inflammatory stages of NASH.
Interleukin-6 (IL-6) is an inflammatory cytokine closely associated with metabolic disease (14–18). It has both pro- and anti-inflammatory properties and may play a role in the progression of NAFLD to NASH (19, 20), cirrhosis (21) and liver cancer (22). IL-6 signaling involves either activation of a membrane-bound receptor (classical) or formation of a signaling complex with a soluble IL-6 receptor (sIL-6R) found in circulation (trans-signaling). Both types of signaling require the ligand/receptor complex to bind a co-receptor, glycoprotein 130 (gp130), on the surface of cells. Circulating sIL-6R, shed from receptor-expressing cells, can dock to gp130 on distant cells, allowing IL-6 (trans)signaling in tissues that do not express the IL-6 receptor. A soluble, secreted form of the co-receptor (sgp130) also circulates and can inhibit trans-signaling by sequestering the IL-6/sIL-6R complex (23). Based on expression patterns, hepatocytes and hepatic stellate cells may be a rich source of sIL-6R and sgp130, respectively, suggesting the liver as a potential major player in IL-6 trans-signaling (24, 25). There is evidence linking increased IL-6 trans-signaling to other metabolic diseases (17, 26, 27), as well as to alcohol- and infection-induced chronic liver disease (28).
Given the involvement of IL-6 signaling at all stages of NAFLD and its link to metabolic disease, we hypothesized that circulating mediators of the IL-6 trans-signaling pathway may predict fatty liver disease severity in populations with MASH and/or morbid obesity. In this study, we investigated relationships between circulating IL-6, sIL-6R and sgp130, and liver pathology associated with metabolic syndrome in two human cohorts.
RESULTS
Characteristics of patients with MASH
Anthropometric, metabolic, and clinical characteristics of our first cohort (28 women and 21 men) are presented in Table 1. Subjects chosen met the current criteria of MAFLD/MASH, which requires either obesity, diagnosis of type 2 diabetes mellitus, or evidence of metabolic dysregulation (9) and elimination of other possible NAFLD causes. Among the population: 44.9% of subjects had diabetes; 49.0% had obesity; 28.6% had hyperlipidemia; 44.9% had hypertension; and, 40.8% consumed alcohol in moderation. Hepatic histological evaluation of patients revealed that 75.0% had stage 2 steatosis, 70.0% had an activity score of 2, and 87.5% had stage 1-2 fibrosis at the time of the biopsy (Table 1), suggesting that most were in the early stages of MASH. MRI/MRE measures of liver fat, volume and stiffness were affected by comorbidities in patients with MASH. Liver fat fraction was lower in patients with a previous history of HCC (p = 0.004); liver volume was significantly higher in patients with diabetes (p = 0.004) and/or obesity (p = 0.004); and liver stiffness was higher in patients with diabetes (p = 0.005), obesity (p = 0.003), hypertension (p = 0.022) and/or a previous history of HCC (p = 0.003) (Supplementary Figure 2A-R). As expected in MASH, plasma globulin, ALT, GGT, glycemia and HbA1c levels were higher than normal ranges in the Canadian population (29) (Table 1). In line with expected sex differences, blood platelet count was lower (p = 0.036), and plasma AST (p = 0.034) and GGT were higher (p = 0.017), in men compared to women (29).
In healthy subjects, normal levels of plasma IL-6 are <3 pg/ml (29, 30) and average plasma concentrations are 35 ng/ml for sIL-6R and 217 ng/ml for sgp130 (30, 31). In this cohort, average plasma IL-6 and spg130 were above those reported for healthy subjects, while average levels of sIL-6R were close to normal (Table 1). Plasma IL-6 was significantly higher in MASH patients with diabetes (p = 0.029), hypertension (p = 0.005) or previous history of HCC (p = 0.025) compared to MASH patients without these co-morbidities (Figure 1A, 1J, and 1P). Plasma IL-6 also correlated positively with BMI (r = 0.34, p = 0.018, Supplementary Table 1) as previously reported in literature (15, 16). Plasma sIL-6R was higher in patients with diabetes (p = 0.007) (Figure 1B) while plasma sgp130 was higher in patients with hypertension (p = 0.034) or history of HCC (p < 0.0001) (Figure 1L and 1R).
Plasma sgp130 predicts the severity of liver stiffness independent of age, sex, adiposity, and comorbidities in patients with MASH
Liver volume and liver stiffness (resistance of the tissue to deformation) are both accurately measured by MRI/MRE over the entire liver (32) and measurements of liver stiffness by MRI/MRE have good prognostic value to predict liver disease severity (33–35). Increased liver volume is influenced by several factors including fat content, inflammation, and/or edema, while increased liver stiffness is more indicative of increased inflammation and/or fibrosis (36). Of the three components of the IL-6 pathway, plasma sgp130 negatively correlated with liver fat fraction (r = −0.31, p = 0.031) (Figure 2C), while plasma sIL-6R positively correlated with liver volume (r = 0.36, p = 0.011) (Figure 2E). Both plasma IL-6 (r = 0.43, p = 0.002) and plasma sgp130 positively correlated with liver stiffness (Figure 2G and 2I), with plasma sgp130 showing the strongest association (r = 0.77, p<0.0001).
Inclusion of subjects with a previous history of HCC allowed broader representation of liver disease severity within our cohort. While these subjects had no detectable cancer at the time of blood sampling, they had expectedly higher levels of liver stiffness (5416.5 ± 2517.9 Pa) and lower liver fat (4.6 ± 2.8 %). Concerned that this variable could influence the associations, we performed similar analysis while excluding subjects with history of HCC. In this subset (n = 39), the correlation between plasma sIL-6R and liver volume was strengthened (r = 0.457, p = 0.004) and the association between sgp130 and liver stiffness remained robust (r = 0.692, p<0.0001). The relationship between plasma IL-6 and liver stiffness also persisted (r = 0.340, p = 0.0342); however, the association between sgp130 and fat content disappeared (r = −0.003, p = 0.985) and a positive correlation between plasma sIL-6R and fat content (r = 0.358, p = 0.0223) was revealed. While this showed that general relationships between IL-6 signaling components and liver volume and stiffness in MASH were not significantly influenced by a history of HCC, this analysis promoted us to evaluate whether associations were dependent on other NAFLD comorbidities.
Stepwise forward regression analysis was performed to explore the influence of age, sex, BMI, and the other diseases associated with metabolic syndrome (diabetes, hyperlipidemia, hypertension, or previous HCC) on associations between IL-6 components and MASH severity (liver fat, volume and stiffness) (Table 2). Adjustment for any of these factors eliminated the associations of plasma IL-6 with all MRE/MRI measures of liver disease. On the other hand, plasma sIL-6R predicted liver volume independent of age, sex, BMI, and hyperlipidemia, hypertension or previous HCC history, but the association was influenced by diabetes status. Plasma sgp130 predicted liver stiffness independent of age, sex, BMI and any comorbidity measured, including diabetes (Table 2).
Histological evaluation of liver biopsy (steatosis, active liver inflammation and liver fibrosis) is the gold standard to diagnose MASH and define the range of disease severity. Histological scoring data (steatosis grade, activity score or fibrosis stage) were available from biopsies performed 6-12 months prior to imaging and blood collection. Plasma IL-6, sIL-6R or sgp130 were not significantly different between subjects with varying MASH severity by histological grading (Supplementary Figure 3). While this appears inconsistent with our findings above using imaging, there may have been changes in liver health during the time between liver biopsy and blood collection. Consistent with this possibility, there were also no correlations between histological scores and plasma ALT (Supplementary Figure 4A, 4B, 4C) or liver stiffness (Supplementary Figure 4D, 4E, 4F).
In contrast, imaging and blood collection were performed on the same day, and multiple serum predictors of liver disease showed strong correlations with MRI/MRE measured liver stiffness, including INR-PT (r = 0.61, p < 0.0001), platelet count (r = −0.59, p < 0.0001), ALP (r = 0.55, p < 0.0001), GGT (r = 0.67, p < 0.0001), albumin (r = −0.52, p < 0.0001) and CK-18 (r = 0.38, p = 0.039) (Table 3). In addition, liver stiffness correlated strongly with all fibrosis risk scores. For these reasons, we focused on imaging-based measurements of liver disease severity for this cohort. We noted that the association of liver stiffness with plasma sgp130 was stronger than all other metabolic risk factors, liver disease scores or individual predictors measured (Table 3). These data suggest that the level of sgp130 in plasma is tightly associated with the extent of liver damage and stage of liver disease.
Plasma sgp130 predicts the severity of liver stiffness better than currently used plasma predictors of liver disease or calculated liver risk scores in patients with MASH
We next evaluated the strength of relationships between plasma components of the IL-6 pathway and: 1) metabolic risk factors (age, BMI, plasma glucose, glycated hemoglobin and lipids); 2) common risk scores used to estimate severity of the liver disease (NAFLD fibrosis score, FIB-4 index and APRI score); and, 3) individual plasma parameters used to estimate NASH risk/severity and calculate risk scores (Supplementary Table 1). Plasma IL-6 correlated with common plasma predictors of liver disease including ALP (r = 0.47, p = 0.001), GGT (r = 0.30, p = 0.044), and albumin (r = −0.56, p < 0.0001), as well as scores used to assess MAFLD severity (NAFLD fibrosis score and FIB-4 index, r = 0.50, p = 0.001; r = 0.37, p = 0.011 respectively). In contrast, plasma sIL-6R correlated more readily with metabolic risk factors including glycemia (r = 0.55, p < 0.0001) and triglycerides (r = 0.37, p = 0.011) compared to other IL-6 components, while also correlating with the liver damage marker GGT (r = 0.50, p < 0.0001) and fibrosis risk scores. Plasma sgp130 correlated very strongly with multiple markers of liver damage including platelet count, INR-PT (r = 0.65, p < 0.0001), ALP (r = 0.62, p < 0.0001), GGT (r = 0.55, p < 0.0001), albumin (r = −0.66, p < 0.0001) and all fibrosis risk scores, but not with glycemia or lipid parameters. These data are in line with IL-6 signaling being a major player in metabolic disease (15, 17, 18, 26, 37–39) and suggest that different components of the IL-6 pathway (i.e. ligand and soluble receptors) may be differentially involved in the metabolic versus inflammatory aspects of fatty liver disease.
Given the strength of the associations of the three IL-6 signaling mediators with image-based measurements of liver disease and multiple markers of metabolic liver disease, we used stepwise linear regression analysis to identify the best predictors of liver fat, volume and stiffness, without or with adjustment for sex (Table 4). Sex was considered as a variable since many metabolic diseases including NAFLD and NASH are influenced by sex, and the liver is one of the most sexually dimorphic tissues (40–42). Independent parameters entered in the model were: the 3 components of the IL-6 pathway (IL-6, sIL-6R and sgp130); plasma parameters related to liver disease (globulin, INR-PT, ALP, total bilirubin and GGT); and liver disease scores (NAFLD fibrosis score, FIB-4 index and APRI score). Age and BMI were not adjusted for as they are used to calculate NAFLD fibrosis score and FIB-4 index.
Plasma ALP and GGT were the primary predictors of liver fat, predicting 26% of the intersubject variability, followed by plasma spg130, predicting an additional 9% (Table 4). All other parameters were excluded and adjustment for sex had little effect. Plasma sgp130 was retained in the regression model to predict liver volume after plasma GGT, FIB-4 index and NAFLD fibrosis score; however, adjustment for sex excluded sgp130 from this model. On the other hand, plasma sgp130 was the primary predictor of liver stiffness, alone explaining 64% of intersubject variability, while NAFLD fibrosis score and plasma globulin together explained an additional 9%. All other independent parameters were excluded and adjustment for sex had little effect to predict liver stiffness.
Plasma sgp130 correlates with advanced liver disease in patients with morbid obesity
Anthropometric, metabolic, and clinical characteristics of patients with morbid obesity are presented in Table 5. In line with published correlations between obesity and NAFLD (1), 98% of these patients had hepatic steatosis. Among the population, 40.4% had diabetes, 81.9% had MAFLD, and 16.5% had MASH. However, despite extreme obesity (average BMI = 48), most were classified as MAFLD (not MASH) based on Bedossa scoring (43), with a high percentage of simple steatosis (70.6%, G0-G1), low activity score (81.9%, A0-A1), and low fibrosis stage (59.6%, F0-F1) assessed from the liver biopsies. Histological evaluation showed generally low scores for immune cell infiltration and hepatocyte ballooning, with average NAFLD activity scores (NAS) falling between 2-3 (Table 5). Taken together, the overall degree of liver disease severity was low in this cohort of morbid obesity, suggesting that this cohort had less active MASH compared to our first cohort at the time of liver biopsy.
Average plasma IL-6, sIL-6R and spg130 were all above levels reported in healthy subjects (30, 31) (Table 5). Similar to the MASH cohort, plasma IL-6 (p = 0.029), sIL-6R (p = 0.043) and sgp130 (p < 0.0001) were higher in subjects with diabetes (Supplementary Figure 5A-C). Consistent with being a risk factor for advanced disease, subjects with diabetes had higher steatosis grade, activity score and fibrosis stage determined using liver biopsy samples (Supplementary Figure 6A-C); however, BMI did not vary with any of these histological parameters (Supplementary Figure 6D-F). After adjusting for age, sex, BMI and diabetes, neither IL-6, sIL-6R nor sgp130 corresponded with steatosis grade or activity score. However, plasma sgp130 was higher in subjects with an F4 fibrosis stage compared to subjects at all other stages (F0-F3) (p = 0.0001) (Figure 3). In contrast, higher plasma ALT corresponded with G3 steatosis grade (p = 0.019), A1 activity score (p = 0.0249), and F2 fibrosis (p = 0.0149), but levels plateaued quickly at early stages (Supplementary Figure 6G, H, I). Other serum parameters of liver damage were not available from this biobank. Thus, in this cohort, our analysis suggests that increased circulating sgp130 corresponded to and was specific for advanced liver fibrosis.
DISCUSSION
In this study we explored relationships between three components of the IL-6 signaling pathway (IL-6, sIL-6R and sgp130) and metabolic fatty liver disease evaluated by liver biopsy, MRE/MRI, plasma predictors and/or liver risk scores in two separate cohorts with diagnosed MASH or morbid obesity. We provide novel data showing that: 1) plasma concentrations of IL-6 transsignaling mediators are differentially increased in obesity, diabetes, hypertension and/or previous history of HCC in patients with MASH; 2) plasma sgp130 strongly predicts MRE/MRI measured liver stiffness in MASH, independent of age, sex, BMI and any comorbidity (diabetes, hypertension, hyperlipidemia or history of HCC), performing superior to many other plasma predictors and liver risk scores; 3) plasma sgp130 is increased in advanced fibrosis independent of age, sex, BMI and diabetes in patients with morbid obesity; and finally, 4) plasma sIL-6R predicts MRE/MRI measures of liver volume, independent of age, sex and BMI.
Circulating inflammatory factors that associate tightly with MAFLD are not well known. We found higher levels of circulating IL-6 and sgp130 in patients with MAFLD/MASH, and that all three circulating IL-6 transignaling components are further increased by concomitant diabetes. This is in line with previous associations of these factors with metabolic syndrome (17–20, 39). However, our new data reveal a strong link between plasma levels of these factors and fatty liver disease, a pathology often occurring concurrently with diabetes and other metabolic diseases. Our data support liver as a possible major source for circulating components of IL-6 transsignaling, and that transignaling may contribute to the development or progression of MAFLD and other metabolic diseases.
The role of IL-6 transsignaling on liver function is controversial. Activation of IL-6 transsignaling promotes tumor formation in a mouse model of HCC (44) and selective blockade of IL-6 transsignaling by sgp130 decreases liver regeneration following partial hepatectomy (45). These studies suggest that IL-6 transsignaling may promote hepatocyte proliferation. Ablation of total hepatic IL-6 signaling in mice causes steatosis and fibrosis (46), but the contribution of classical versus trans-IL-signaling is not clear. In line with their human data (46), we also show no significant association between plasma sIL-6R and liver fibrosis stages in our cohorts; however, we find circulating sIL-6R significantly associated with multiple metabolic aspects of MAFLD including glycemia and lipidemia, hepatic fat fraction, as well as total liver volume. Interestingly, larger liver volume was recently associated with a 3-fold increase in risk of all-cause mortality in NAFLD (47). While we observe strong associations of sgp130 with liver stiffness and advanced fibrosis in humans with NASH, sgp130 treatment does not impact NASH development in mice fed a western diet (48). More work is needed to determine whether local or systemic secretion of IL-6-transsignaling factors has a direct effect on the progression of MAFLD/MASH.
Using stepwise linear regression analysis, we show for the first time that plasma sgp130 alone is a strong predictor of liver stiffness in MASH with a wide range of disease severity, better than commonly used serum predictors and risk scores. When combined with NAFLD fibrosis score and plasma globulin, the predicative power of the model increased. Thus, we propose that plasma sgp130 (alone or in combination with NAFLD fibrosis score) could be an effective, non-invasive method to predict the severity of liver disease in patients with MASH. Identification of non-invasive, reliable predictors of hepatic inflammation and fibrosis, especially at early stages (F1 and F2), will facilitate diagnosis and possibly increase success rates of emerging interventions and treatment strategies. We present this as a hypothesis-generating study and recognize that a larger sample size is required to confirm IL-6 transsignaling proteins as biomarkers of MAFLD.
One limitation of our study is the use of different methods to quantify liver disease severity across cohorts. For the first, we relied mainly on MRI and MRE analysis to assess severity of liver disease, while for the second cohort we used scoring data from liver biopsies. Liver biopsy data was available for the MASH cohort; however, the sample was taken up to 12 months prior to the blood draw and disease severity may have changed over that time. Liver biopsy remains gold standard, yet accumulating studies show that MRI/MRE also has high diagnostic accuracy for liver fibrosis, that is similar to biopsy in subjects with NAFLD (33–35). Regardless, the two methods of assessing liver disease severity limited our ability to directly compare liver disease parameters (steatosis grade, activity scores, and fibrosis stage) between the cohorts.
Liver fibrosis and liver stiffness are also not interchangeable. Liver fibrosis in biopsies is scored based on collagen staining, while MRE-determined liver stiffness is influenced by influenced to varying degress by fibrosis, inflammation, and steatosis (33–35). This may explain discrepancies between sgp130 sensitivity to predict stiffness versus fibrosis across the two cohorts. Interestingly, if we consider this differences, our data could also suggest that spg130 is a sensitive marker of active liver inflammation and damage, but not a predictor of collagen deposition per se. Sgp130 is a component of an inflammatory signaling pathway, which could explain stronger correlation with liver stiffness versus collagen staining (a late consequence of damaging stimuli). In line with this theory, sgp130 was significantly increased only in late fibrosis stages determined by histology, while the strong linear correlation between plasma sgp130 and liver stiffness measured by MRI/MRE spanned across all stages of disease severity.
In conclusion, our data support that circulating components of the IL-6 signaling system may play a role in MAFLD/NASH pathogenesis and have potential to serve as sensitive predictors of liver damage associated with metabolic disease. Importantly, relationships between plasma sIL-6R and particularly spg130 with metabolic liver disease severity are independent of sex, age and BMI. Their individual associations with either metabolic or inflammatory aspects of MAFLD suggest interesting mechanistic roles in disease progression. Their strong associations with liver fat, volume and stiffness may also be useful for non-invasive monitoring of the early MAFLD to MASH transition.
METHODS
Patients with MASH
MASH had been previously confirmed by liver biopsies performed 6-12 months prior to MRI/MRE, according to the clinical standard of care using 16-G or 18-G core needles. Hematoxylin and eosin (H&E) slides were scored by liver pathologists, with fibrosis stage, inflammation grade, and steatosis grade assessed according to the NASH Clinical Research Network (NASH CRN) histological scoring system(49).
Subjects were aged 18 years and older, diagnosed with MASH or with HCC on a MASH background, able to undergo magnetic resonance imaging (MRI) without administration of a contrast agent and to understand instructions in either French or English. Subjects were excluded if they: had high alcohol consumption (>10 drinks/week for women and >15 for men); had liver disease other than MASH; were taking medications associated with steatosis (e.g. amiodarone, valproate, tamoxifen, methotrexate or corticosteroids); physically unable to fit in the MRI machine; had contraindications to MRI; or were pregnant or wished to be pregnant during the study-year.
Of the 493 patients, 133 patients were eligible for this study, of whom 89 were diagnosed with MASH and 44 with MASH and a previous history of HCC. Between May 2018 and June 2019, which represents six to twelve months after the liver biopsies and diagnosis of MASH, the 133 patients were invited back for an MRI/MRE scan to assess liver fat fraction, volume and liver stiffness. Among these, 49 subjects with MASH alone and 34 with HCC and MASH did not participate in the study for the following reasons: refusal (n=40), unreachable (n=25), cancellation (n=8), distance to the hospital (n=8), language issues (n=1) or other reasons (comorbidities including amputations) (n=1). Thus, this analysis includes 40 patients with NASH (28 women, 22 men) and 10 patients with MASH and a history of HCC (4 women and 6 men) from the original registry of 493 patients. All patients (n=50) signed a consent form to be included in this study, approved by the Centre hospitalier de l’Université de Montréal human ethics committee (IRB #17.031). A flowchart of inclusion/exclusion criteria used is shown in Supplementary Figure 1.
MRI/MRE examinations were performed using a 3.0 T clinical scanner (Skyra; Siemens Healthineers, Mountain View, California). Proton density fat fraction (PDFF), liver volume (voxels, cm3), and liver stiffness (Pa) were measured as quantitative predictors of liver fat, volume and fibrosis, respectively. Average PDFF values for the entire liver volume were obtained using the LiverLab package (Magnetom Aera, Software version VE11C, Siemens Healthcare GmbH, Erlangen, Germany). Liver stiffness measurements by MRE were performed according to previously described methods (50). Fasting plasma samples were collected on the day of the MRI/MRE and stored at −80°C until measurement of plasma IL-6 parameters and other predictors for calculation of the liver disease scores.
Calculation of fasting NAFLD fibrosis score, FIB-4 index and APRI score
The following equations were used to calculate liver scores in the NASH population:
NAFLD fibrosis score = −1.675 + [0.037 x age (years)] + [0.094 x BMI (kg/m2)] + [1.13 × impaired fasting glucose/diabetes (yes=1, no=0)] + [0.99 × AST/ALT ratio] - [0.013 × platelet count (× 109/L)] – [0.66 × albumin (g/dl)] (51).
Fibrosis-4 (FIB-4) index = [Age (years)] × [AST (U/L)] / [platelet count (109/L) × sqrt (ALT) (U/L)] (52).
AST to platelet count ratio index (APRI): [AST (IU/L)/40 IU/L] / [platelet count (× 109/L)]] × 100 (53).
Patients with morbid obesity
Adult men and women (18 years and older) were selected from a registry of 4,781 patients undergoing bariatric surgery. Exclusion criteria were: having high alcohol consumption (>10 drinks/week for women and >15 for men); or having liver disease other than NASH (e.g. autoimmune hepatitis, Wilson disease, hemochromatosis or HBV, HCV or human immunodeficiency viruses). Thus, this retrospective analysis included a subpopulation of 245 subjects (123 men, 122 women) selected based on availability of serum and histological scoring data.
Random blood samples were collected on the night before bariatric surgery and stored immediately at −80°C until time of analysis. Sampling procedure and position were standardized among all surgeons. Liver samples were obtained by incisional biopsy of the left lobe and were not cauterized. Grading and staging of histological liver sections were performed using the protocol of Brunt et al (54) by pathologists blinded to the study objectives. Bedossa algorithm (43) was used to diagnose NASH, using liver biopsy histological scores for hepatocellular ballooning stage (0-2), lobular inflammation (0-2), steatosis grade (G0-G3), activity score (A0-A4) and fibrosis stage (F0-F4). These were also used to calculate NAFLD activity (NAS) and steatosis, activity and fibrosis (SAF) scores (43).
Measurements of plasma IL-6, sIL-6R, sgp130, Cytokeratin-18 and Alanine aminotransferase (ALT)
Commercial ELISA kits were used to measure plasma concentrations of IL-6, sIL-6R and sgp130 (R&D systems, Human Quantikine ELISA kits, D6050, DR600, and DGP00 respectively) and Cytokeratin-18 (Peviva, M30 Apoptosense ELISA, 10011). For sgp130 and sIL-6R, samples were diluted 1:100, while for IL-6 and Cytokeratin-18, undiluted samples were used. ALT was measured in patients with morbid obesity using a commercial kit (SGPT liquid ALT reagent set, Pointe scientific A7526). All assays performed according to manufacturer’s instructions and quality controls/standard provided with kits were included. Sample analysis was blinded using subject identification code.
Statistics
Data in Table 1 and Table 5 are presented as mean±SD for continuous variables and as the number of subjects (n) and percentage (%) within the subpopulation for categorial variables. Normality was evaluated using a Kolmogrov-Smirnov test. When normality failed, data was log transformed (log10). Outlier and influencer points were identified using SPSS. One subject in the NASH cohort was a strong influencer for plasma IL-6 and sIL-6R in all analysis. Thus, this subject was excluded from analyses. Unpaired t-test, one-way ANOVA and Pearson correlation were used to analyze parametric or log10 transformed data. Given the variability of some continuous data in tables 1 and 5, sensitivity analysis was performed using nonparametric Mann-Whitney U test for intergroup differences to validate significant findings. For categorical variables, chi-square test was used for count >5 in each cell, otherwise Fisher’s exact test was used.
Clinical endpoints used were liver fat fraction, volume and stiffness for the NASH cohort, and liver steatosis grade, activity score and fibrosis stage for the morbid obesity cohort. To assess relationships between circulating IL-6, sIL-6R and sgp130 and these endpoints, stepwise forward regression analysis was used to predict measurements of liver disease severity using log10 [plasma IL-6], log10 [plasma sIL-6R] and log10 [plasma sgp130] as independent variables with adjustment for age, sex, and BMI and comorbidities (diabetes, hyperlipidemia, hypertension or previous HCC history) for the NASH cohort. For the morbid obesity cohort, log10 [plasma IL-6], sIL-6R and sgp130 were used as independent variables with adjustment for age, sex, BMI and diabetes. For categorical variables (NASH cohort: liver steatosis stage, activity score and fibrosis stage. Morbid obesity cohort: liver steatosis stage, activity score and fibrosis stage, Bedossa NASH diagnosis, NAS and SAF score), univariate analysis was used with adjustment for age, sex, BMI and diabetes. Data was analyzed using IBM SPSS (Version 27) and GraphPad Prism (Version 8) and significance was set at P < 0.05.
Study approval
Patients with MASH
Participants were selected from a registry of 493 patients with available liver biopsies collected under ethical approval (IRB: #15.147) at Centre hospitalier de l’Université de Montréal that was initiated in February 2016. All selected participants (n = 50) in the current study (IRB # 17.031) provided written, informed consent allowing preservation and subsequent use of their data. All participants provided written, informed consent allowing preservation and subsequent use of their plasma samples and data.
Patients with morbid obesity
Plasma samples and matching liver biopsies were obtained from the Biobank of the Institut universitaire de cardiologie et de pneumologie de Québec – Université Laval in compliance with Institutional Review Board-approved management policies initiated in 2002 and still ongoing. Liver biopsy samples were collected at the time of bariatric surgery under ethical approval (IRB: #1142).
Data Availability
All data produced in the present study are available upon reasonable request to the authors
Author Contributions
JLE and AT designed studies. CH and AB managed the NASH clinical study. L Biertho and AL designed and carried out the bariatric surgery study. AT, L Bilodeau, AG and C Baldwin acquired, analyzed and interpreted data. MF assisted with data analysis. JLE, L Bilodeau, AT, MF, ML, EG, MB, C Bémeur contributed intellectual content. CYC scored liver histological samples. AG, MF and JLE wrote the manuscript. All authors reviewed the manuscript.
SUPPLEMENTARY MATERIAL
Acknowledgements
Authors acknowledge the invaluable collaboration of the surgery team, bariatric surgeons, and biobank staff of the IUCPQ. We would also like to thank Paule Bodson Clermont, Hannah Zhang, Mélissa Léveillé and Stewart Jeromson for their assistance.
Footnotes
Abbreviations: NAFLD, non-alcoholic fatty liver disease; NASH, non-alcoholic steatohepatitis; MAFLD, metabolically-associated fatty liver disease; MASH, metabolically-associated steatohepatitis, MASH; hepatocellular carcinoma, (HCC); interleukin-6, IL-6; soluble interleukin-6 receptor (sIL-6R); glycoprotein 130, gp130; magnetic resonance spectroscopy, MRS; imaging, MRI; elastography, MRE; alanine and aspartate aminotransferases, ALT and AST; gamma glutamyl transferase, GGT; NAFLD activity score NAS; NAFLD activity fibrosis score, SAF; FIB-4 index, fibrosis-4 index; APRI score, AST to platelet platelet ratio.
Conflict of interest statement: Laurent Biertho receives funding from Johnson & Johnson Medical Companies, Medtronic, Bodynov and GI Windows for studies on bariatric surgery. An Tang and Jennifer Estall received a speaking honorarium from Eli Lilly. A Gunes, L Bilodeau, C Huet, A Belblidia, C Baldwin, JM Giard, A Lafortune, CY Couture, BN Nguyen, E Galun, C Bémeur, M Bilodeau, M Laplante, and May Faraj have nothing to disclose.
Data Transparency: The data, analytic methods and study materials is available for other researchers upon request.
Funding: This work was supported by unrestricted operating grants from the Merck, Sharpe and Dohme Corporation/University of Montreal and the International Development Research Centre (108591-001). JLE and AT are supported by Chercheurs-boursiers: Senior awards from the Fonds de recherche du Québec Santé. The ULaval Biobank is supported by the IUCPQ Foundation and Research Center.
Figures have been reformatted. Introduction and discussion have been modified to include mechanistic implications of our work.