Abstract
BACKGROUND:
Non-alcoholic fatty liver disease (NAFLD) comorbidity with adipose tissue dysfunction is not new and studies have focused on how adipose tissue influences NAFLD pathophysiology.
OBJECTIVE:
Quantification of nature and magnitude of the association between serum resistin and also insulin resistance, by calculating pooled Homeostatic Model Assessment for Insulin Resistance (HOMA-IR) score, with NAFLD pathophysiology was the objective of the current study.
METHODS:
Using systematic review and meta-analysis and standardized mean difference (SMD) as the effect size, the levels of resistin and HOMA-IR scores have been investigated in NAFLD subjects in comparison with controls in the case-control studies using random-effects models.
RESULTS:
This meta-analysis retrieved a total number of 665 and 522 cases and 671 and 555 control subjects until May 2020 for serum levels of resistin and HOMA-IR score until May 2020. The final analyses demonstrated that pooled SMD of resistin and HOMA-IR score was 0.687 (95% confidence interval, 0.070–1.304) and 1.368 (95% confidence interval, 1.080–1.655); respectively. Moreover, the p-value for the test of significance for each pooled SMD was examined by the z-test and calculated as 0.029 and 0.000 for resistin and HOMA-IR score (clearly considered as statistically significant).
CONCLUSION:
Based on the findings, the HOMA-IR score and the serum levels of resistin in NAFLD subjects are associated with disease pathogenesis.
Background
Non-alcoholic fatty liver disease (NAFLD) is highly prevalent worldwide and manifests in pathological stages comprising a wide progressive spectrum of diseases ranging from simple steatosis (i.e., more than 5% of the liver weight is composed of fat) to non-alcoholic steatohepatitis (NASH), in which steatosis is combined with inflammation and fibrosis with a significant risk at later stages for the development of cirrhosis [1]. NAFLD is associated with significant adverse outcomes both through liver-specific complications which promotes both morbidity and mortality [2]. It is believed that NAFLD is closely related to obesity and insulin resistance (type 2 diabetes), and these conditions augment the progressiveness of the disease-associated complications toward more advanced stages [3]. NAFLD pathogenesis is considered multifactorial; however, the mechanisms that play role in hepatic steatosis and liver predisposition to inflammation, fibrosis and consequently cirrhosis are still not fully understood but reflect complex interactions among metabolic target tissues including hepatocytes, adipose, skeletal muscle, and immune system cells [4]. It is believed that NAFLD is a hepatic manifestation of metabolic syndrome. Factors like genetic predisposition, increased gut microbial-derived components, dietary intakes (e.g., high fat and fructose intake), inflammation, insulin resistance and lipotoxicity are presumed to contribute to the NAFLD development and progression [5–7]. The adipose tissue, which is not only a lipid-storage organ but also an endocrine organ, closely interacts with liver by secreting polypeptides known as adipokines [8, 9]. A growing body of literature demonstrates that adipokines are involved in various processes, such as immunity, inflammation, insulin sensitivity, hepatic steatosis, and NASH [10]. Resistin, leptin and adiponectin are known as main adipokines playing role in NAFLD development [10]. In the current study, the serum levels of resistin, and the reported HOMA-IR scores in the included records were investigated in NAFLD pathogenesis using systematic review and meta-analysis.
Methods
Protocol of the systematic review and meta-analysis
The PRISMA Checklist 2009 was used to conduct this systematic review and meta-analysis [11]. This investigation was approved by the Ethics Committee of Ilam University of Medical Sciences (ethical code: IR.MEDILAM.REC.1398.183).
Information sources and search strategies
The literature sources of information for the systematic review were Pubmed, Pubmed Central, Medline, Google Scholar, Embase and SCOPUS specialized biomedical databases and independently carried out by two authors (A. M., A. A.) without the beginning date restriction until May 2020. The search key terms were framed along with three major categories: NAFLD pathological stages; resistin blood levels and study type. Publications using the MeSH and non-MeSH terms “non-alcoholic fatty liver disease”, “NAFLD”, “simple steatosis”, “steatohepatitis”, “cirrhosis” in combination with “resistin” with additional keywords such as “peripheral blood samples”, “circulatory” and “case-control study” were identified. Then, the search results were restricted to the English language and the review possesses were limited to case-control studies by two authors independently, which fulfilled our inclusion and exclusion criteria. Apart from the databases aforementioned, the harvested references from retrieved studies were also subjected to a manual search of references.
Eligibility criteria
Only original articles that had been published in peer-reviewed journals were pre-included. On the other hand, studies applying NAFLD reviews and theories and models, to discuss phases of NAFLD development and management, conference papers, non-English language articles were excluded. Studies were pre-included if they used standardized methods such as enzyme-linked immunosorbent assay (ELISA) and chemiluminescent methods for the detection of serum resistin levels. No limitation was applied to the subtypes of NAFLD, the severity of NAFLD, sex or race of the study participants. However, studies were excluded if they enrolled populations other than NAFLD as well as reports describing NAFLD interventional therapies. Moreover, only records reporting serum resistin levels and their corresponding HOMA-IR score were considered for effect size analyses in each study in the case and control groups.
Study selection
First, records excluded with no NAFLD data, records with no full text available (conference proceedings), studies that did not report original data (editorials or reviews), records that did not report the number of participants, and case reports. Publications were assessed based on their titles and then based on abstracts; to exclude studies that did not meet the inclusion and exclusion criteria. The full texts of the remaining papers were archived and read, and at this stage, studies that not fulfilling the inclusion criteria were rejected. Only the case-control studies providing enough information about their results in such a way to compute an estimate of the effect size (such as mean and standard deviation or standardized error of the mean, etc.) for serum resistin concentrations were included. The data must be reported in an exact numerical format and not extracted from the figures of the record by estimation and if the data not reported in such a way, the record was excluded. Due to the different protocol definition used for assessment of serum resistin levels, only studies which determined the levels of resistin in NAFLD in comparison with non-NAFLD subjects were considered to be included in the current meta-analysis.
Data collection process
The first author of the selected records, date of publication, the serum levels of resistin in NAFLD patients and non-NAFLD controls, evaluated criteria for the NAFLD and also the total number of cases and controls and other related information were extracted from the publications which have been provided from the systematic review for the meta-analytical process. Moreover, the HOMA-IR was also meta-analyzed for the selected studies, if the finally included records reported such information.
Summary measures and synthesis of results
Stata version 14.0 (Stata Corporation, College Station, TX, USA) was employed for data analyses of the systematic review and meta-analysis. Between-study heterogeneity was assessed using the χ2 based Q-test and also I2. The Q test and I2 statistics were applied to evaluate the inconsistencies and heterogeneities among the records. A significant Q suggests the existence of heterogeneities and I2 estimates the magnitude of the inconsistencies among the studies [12]. For analysis of the estimated total effect size, the random-effects model and SMD were used. Data were shown as the estimated SMD with a 95% confidence interval (CI) for each study. The significance of the total SMD was examined by the z-test and p < 0.05 was considered as statistically significant.
Risk of bias across studies
For the risk of bias across records, the records were scrutinized for method validation and data processing. For effect size estimation for each selected study, the funnel plot was developed. For interpretation of any publication bias among records, visual inspections of the generated funnel plot were employed to evaluate the plot symmetry. In this plot, the Y and X axes represent the standard error and standard difference in means, respectively. Besides, for the investigation of the HOMA-IR score and serum resistin relationship, the meta-regression (under random-effects model) was performed in such a way that the calculated SMD for HOMA-IR defined as a moderator variable and it was plotted against the SMD for resistin serum levels for the finally included studies.
Results
Study selection
A flowchart, representing the systematic review process for study selection, has been shown in Fig. 1. The initial search for serum resistin concentrations in NAFLD retrieved a total number of 5180 potentially eligible records and, then 543 studies removed being as duplicates. Then, 4600 studies were excluded after reading titles or abstracts as being obviously irrelevant to the goal of this meta-analysis. Moreover, because of insufficient data presented by the authors of such publications for calculation of the SMD and 95% CI and because of poor quality, 3 papers were excluded additionally for serum resistin levels in NAFLD. Finally, the records which fulfilled exclusion and inclusion criteria (14 studies, 18 including within study subgroups, which have been detailed in Table 1A and Table 1B cited in the references part of the current study as [13–26]); were included in the meta-analytical process. The data for HOMA-IR calculation was extracted from the selected studies which have reported serum resistin levels as well in this systematic review and then subjected to meta-analysis (10 studies, 13 including within-study subgroups, which have been detailed in Table 2 and cited in the references part of the current study as [13–15, 26]).

Systematic review process. The flowchart represents the systematic review process for study selection according to the inclusion and exclusion criteria. The results have been shown for serum resistin concentrations and HOMA-IR scores in non-alcoholic fatty liver disease patients vs controls; respectively.
Demographic data for finally included studies reported serum resistin concentrations
Data extracted from the finally included studies to calculate pooled standardized mean difference (SMD) for serum resistin concentrations
* -: NAFLD subtype not reported by authors, NASH: Nonalcoholic steatohepatitis, SS: simple steatosis.
Data extracted from the finally included studies to calculate pooled standardized mean difference (SMD) for HOMA-IR score as reported by finally included records for serum resistin levels
* -: NAFLD subtype not reported by authors, NASH: Nonalcoholic steatohepatitis, SS: simple steatosis.
For each record, the demographic and disease characteristics data which have been reported by the authors were extracted. The number of NAFLD subjects, NAFLD pathological state, number of non-NAFLD controls, NAFLD subtype (if reported by authors), body mass index (BMI), mean age in control subjects and NAFLD patients have been presented in Table 1A. Actually, this meta-analysis retrieved a total number of 665 and 522 cases and 671 and 555 control subjects for serum levels of resistin and HOMA-IR score until May 2020.
Risk of bias within studies
The results showed that all of the selected studies were not homogenous and actually they were inconsistent in their corresponding results. The Q test was calculated as 417.428 (24.713) and 49.975 (14.860) for serum resistin and HOMA-IR score under the fixed-effect model (random-effects model); respectively. It is well accepted that the Q statistic test could be only applied for testing the heterogeneity existence among studies, but it is not suitable for calculation of the heterogeneity extent. Moreover, the I2% tests were calculated as 95.927 (31.211) and 75.988 (18.417) for resistin and HOMA-IR under the fixed-effect model (random-effects model); respectively. As the I2 index and the between-studies variance, τ2, are directly related: the higher the τ2, the higher the I2 index [12]. However, it has been proposed a tentative classification of I2 values with the purpose to interpret heterogeneity magnitude. Thus, the percentages of nearly 25% (I2 = 25), 50% (I2 = 50), and 75% (I2 = 75) would mean low, medium, and high heterogeneity, respectively [12]. Therefore, the random-effects model of the meta-analysis was applied for the presentation of the forest plots of the finally included studies for serum resistin levels as well as HOMA-IR score effect size calculation in NAFLD subjects as compared to controls.
Synthesis of results
The forest plots showed the weighted SMD and standard deviation including the corresponding confidence interval for each study and the overall SMD as well (Fig. 2). The final analyses demonstrated that mean pooled SMD for serum resistin level was 0.687 (95% confidence interval, 0.070–1.304). Moreover, the p-value for the test of significance for the pooled SMD was examined by the z-test and calculated as 0.029, demonstrating a statistically significant effect. The weighted pooled SMD for the association between serum resistin concentration with NAFLD pathogenesis was calculated for 14 studies (18 including within study subgroups, as presented in Fig. 2A). Moreover, the forest plot for the HOMA-IR score which has been reported by the included studies measuring serum resistin levels was developed. The analyses demonstrated that the weighted pooled SMD for the association between the HOMA-IR score and NAFLD pathogenesis was 1.368 (95% confidence interval, 1.080–1.655) applying the random-effects model for 10 studies (13 including within study subgroups) as presented in Fig. 2B. The p-value of the HOMA-IR score for the significance of the combined SMD calculated 0.000, as examined by the corresponding z-test.

Forest plots. For serum resistin (2A), HOMA-IR score (2B), the final analyses demonstrated that standardized mean differences (SMDs) of serum resistin concentration and HOMA-IR score was 0.687 (95% confidence interval, 0.070–1.304) and 1.368 (95% confidence interval, 1.080–1.655); respectively. Moreover, the p-value for the test of significance for pooled SMD was examined by the z-test and calculated as 0.029 and 0.000 for resistin and HOMA-IR score, respectively (clearly considered as statistically significant).
Shapes of the funnel plots were considered to be moderately asymmetrical confirming that there was some publication bias in the related records in reporting serum resistin concentrations and HOMA-IR scores (Fig. 3A and Fig. 3B). For serum resistin (Fig. 3A) and HOMA-IR score (Fig. 3B), this bias mainly refers to the right part of each plot that has been mainly occupied by the publications reporting the higher SMD in NAFLD subjects in comparison with corresponding non-NAFLD controls. The meta-regression analysis showed that the SMD for serum resistin levels was directly proportional to the calculated HOMA-IR score (the moderator) using the random-effects model (Fig. 4).

Funnel plots. For serum resistin (3A), HOMA-IR score (3B), the publication bias mainly refers to the right part of each plot that has been mainly occupied by the publications demonstrating the higher standardized mean difference (SMD) in NAFLD subjects in comparison with corresponding controls.

Meta-regression for calculated HOMA-IR score and standardized mean difference (SMD) for serum resistin levels. As it is evident, the meta-regression analysis showed that SMD for serum resistin levels is directly proportional to the calculated HOMA-IR (the moderator) using the random-effects model and the meta-regression curve roughly showed a linear relationship i.e. higher SMD for serum resistin levels corresponds to higher HOMA-IR and thus more severe insulin resistance.
NAFLD comorbidity with adipose tissue dysfunction is not new and studies have focused on how adipokines might be involved in the NAFLD pathophysiology. Quantification of nature and magnitude of the association between adipokines and also insulin resistance by calculating pooled HOMA-IR score and NAFLD pathophysiology was the objective of the current study. In this study, the serum levels of resistin and the reported HOMA-IR scores in the included records were investigated for their associated roles in NAFLD pathogenesis using a systematic review and meta-analysis. The records were collected by systematic review and then filtered according to inclusion and exclusion criteria. Actually, the number of filtered records for resistin and HOMA-IR was sufficient. According to recently published studies that have been included in the current investigation, the results demonstrated that serum levels of resistin and the calculated HOMA-IR score were clearly considered to be statistically significant. Interestingly, the highest mean weighted pooled SMD between serum level of resistin and the calculated HOMA-IR factors for the association between HOMA-IR and NAFLD pathogenesis was HOMA-IR score (mean pooled SMD = 1.368, 95% confidence interval, 1.080–1.655). Besides, the adipose-derived serum resistin levels were higher in NAFLD than healthy controls implying the adverse effects of this adipokine in NAFLD development. The results show the important role of resistin, even though, the role of other adipokines cannot be ruled out. Another factor we considered was HOMA-IR and the highest mean weighted pooled SMD for the association between HOMA-IR and NAFLD pathogenesis was HOMA-IR score (mean pooled SMD = 1.368). HOMA-IR score is calculated as follows: HOMA-IR score = (Fasting insulin)*(Fasting glucose) / 405. The HOMA-IR is used to measure the severity of insulin resistance, though normal insulin resistance varies depending on the population. Common reference levels for HOMA-IR insulin resistance range from 0.7–2 [27]. As it is evident, the meta-regression analysis showed that SMD of serum resistin levels are directly proportional to the calculated HOMA-IRs (the moderator) using the random-effects model and is actually evident that the curve roughly showed a linear relationship i.e. higher SMD for serum resistin levels causes higher HOMA-IR score and thus more severe insulin resistance. Recent study showed that HOMA-IR is an independent predictor of advanced liver fibrosis [28]. Our data are in accordance with the narrative reviews discussing the relationships among NAFLD, resistin hormone levels, HOMA-IR score and insulin resistance development [29–34]. To best our knowledge, this is the first meta-analysis that reports the serum levels of resistin and their corresponding HOMA-IR score in this regard. Although, the precise mechanisms by which resistin affects fatty liver pathogenesis are unclear, however, investigations are ongoing and some of them shed light on its mechanisms of action. Resent study shows that inflammation appears to play a pivotal role in determination of resistin circulating levels and it proposed resistin as a molecular link between body fat deposition and insulin resistance [35]. Moreover, intracellular signaling for resistin induces Nuclear factor kappa B (NF-κB) and mitogen-activated protein kinase (MAPK or MAP kinase) activations and also it causes increased oxidative stress [35]. Recently published data suggests that acute elevated resistin levels exacerbate mitochondrial damage and aggravates liver steatosis through 5’ AMP-activated protein kinase (AMPK)/Peroxisome proliferator-activated receptor gamma coactivator 1-alpha (PGC-1α) signaling pathway in male NAFLD mice [36]. Moreover, it has been proposed that MiR-34a regulates mitochondrial content and fat ectopic deposition induced by resistin through the AMPK/Peroxisome proliferator-activated receptor alpha (PPAR-α) pathway in HepG2 cells [37]. Future studies may focus on the mechanism of action of resistin and also other adipokines and their relationship with NAFLD in depth. In conclusion, the HOMA-IR score and the serum levels of resistin in NAFLD subjects are associated with disease pathogenesis and one of the implications of this study is that resistin may be considered as a useful diagnostic and prognostic tool in NAFLD.
Footnotes
Acknowledgments
None declared.
Funding
Not funded.
Conflict of interest
The authors do not have any disclosures to report.
Author contributions
A. N. performed the study, acquisition and analysis of data and writing the manuscript. A. M., A. A. and A. S-N performed literature database searching.
