Skip to main content
  • Research article
  • Open access
  • Published:

A circulating microRNA signature as noninvasive diagnostic and prognostic biomarkers for nonalcoholic steatohepatitis

Abstract

Background

Noninvasive biomarkers are urgently needed for patients with nonalcoholic steatohepatitis (NASH) to assist in diagnosis, monitoring disease progression and assessing treatment response. Recently several exploratory studies showed that circulating level of microRNA is associated with NASH and correlated with disease severity. Although these data were encouraging, the application of circulating microRNA as biomarkers for patient screening and stratification need to be further assessed under well-controlled conditions.

Results

The expression of circulating microRNAs were profiled in diet-induced NASH progression and regression models to assess the diagnostic and prognostic values and the translatability between preclinical mouse model and men. Since these mice had same genetic background and were housed in the same conditions, there were minimal confounding factors. Histopathological lesions were analyzed at distinct disease progression stages along with microRNA measurement which allows longitudinal assessment of microRNA as NASH biomarkers. Next, differentially expressed microRNAs were identified and validated in an independent cohorts of animals. Thirdly, these microRNAs were examined in a NASH regression model to assess whether they would respond to NASH treatment.

MicroRNA profiling in two independent cohorts of animals validated the up-regulation of 6 microRNAs (miR-122, miR-192, miR-21, miR-29a, miR-34a and miR-505) in NASH mice, which was designated as the circulating microRNA signature for NASH. The microRNA signature could accurately distinguish NASH mice from lean mice, and it responded to chow diet treatment in a NASH regression model. To further improve the performance of microRNA-based biomarker, a new composite biomarker was proposed, which consists of miR-192, miR-21, miR-505 and ALT. The new composite biomarker outperformed the microRNA signature in predicting NASH mice which had NAS > 3, and deserves further validations in large scale studies.

Conclusion

The present study supported the translation of circulating microRNAs between preclinical models and humans in NASH pathogenesis and progression. The microRNA-based composite biomarker may be used for non-invasive diagnosis, clinical monitoring and assessing treatment response for NASH.

Background

Nonalcoholic steatohepatitis (NASH) is histologically characterized by hepatic steatosis, ballooning, inflammation, and varying amount of pericellular fibrosis [1]. Patients are frequently either asymptomatic or have nonspecific symptoms, which impose a great obstacle for clinical diagnosis. Moreover, NASH patients are often associated with higher risk of developing both cardiovascular disease and type 2 diabetes, therefore early diagnosis and intervention would greatly benefit the patient by preventing the progression of major hepatic and extrahepatic manifestations [2, 3]. Currently, liver biopsy is still required in clinical diagnosis and assessment of the treatment response of NASH, despite its invasive nature, variability in sampling and potential complications [4]. Thus identification and development of circulating biomarkers for clinical applications are urgently needed.

MicroRNAs are small non-coding RNA molecules which are evolutionally conserved epigenetic regulators. A single microRNA could post-transcriptionally regulate hundreds of mRNAs through binding to the 3′-untranslated region and leading to its translational repression [5]. MicroRNAs are also present in extracellular spaces, and serve as intercellular messengers to affect many physiological processes [6, 7]. Circulating microRNA have emerged as attractive diagnostic tools for the assessment of the pathological state of their origin tissue [8,9,10]. Studies suggest that the circulating level of microRNA is associated with liver steatosis, inflammation, fibrosis and carcinogenesis [11, 12]. The performance of circulating microRNA alone or in combination with other serum biomarkers in NASH diagnosis has been explored in a number of clinical studies [13,14,15]. However, the significance of these studies is limited by the heterogeneity in disease status, patients’ genetic variations, and different analysis platforms [16]. For example, hepatic miR-34a is differentially expressed in NASH and physiologically associated with lipid metabolism, inflammation, and hormone signaling in NASH. However, the circulating level of miR-34a in NASH patients is not conclusive [17,18,19]. Given the important role of microRNA in NASH pathophysiology, the value and feasibility of circulating microRNA as non-invasive biomarkers need to be further investigated.

The present study aims at evaluating the diagnostic and prognostic value of circulating microRNA in a well-studied diet-induced NASH mouse model, in which potential genetic and environmental confounding effects are minimized. We started by reviewing the published data on circulating microRNA profile in NAFLD or NASH patients, and shortlisted 25 microRNAs which were reported to be differentially expressed [12, 15, 20]. Then we quantitatively measured the expression of these 25 microRNAs by qPCR in a time-dependent NASH progression model, to understand the association between circulating microRNA and disease evolution. A panel of differentially expressed microRNAs were identified and confirmed in an independent study, and was designated as circulating microRNA signature for NASH. The microRNA signature responded to NASH treatment and correlated well with disease regression. Based on these intriguing findings, the diagnostic power of circulating microRNA signature was assessed and a microRNA-based composite biomarker was proposed. Collectively, we established the translation of circulating microRNA between mice and men in NASH disease progression and regression. The new composite microRNA-based biomarker holds great promise in NASH diagnosis and prognosis, and warrants further evaluation in large scale translational studies.

Methods

Animal model generation

Male C57BL/6 mice were purchased from Nanjng Biomedical Research Institute. Mice were housed in a 12 h light-12 h dark cycle in a 21-23 °C facility. To induce NASH, mice at 4 weeks of age were fed ad libitum a high fat, high fructose and high cholesterol diet (3H diet, Research Diets Inc. D09100301). Control mice were fed a standard chow diet (LabDiet, 5C02C). Control and 3H diet mice were anesthetized in 5% isoflurane at varying time points following initiation of dietary intervention to assess liver histological lesions as described before [21]. Blood was collected by cardiac puncture, and plasma was separated by centrifugation at 10000 rpm for 10 min at 4 °C. Plasma samples were aliquoted and stored in − 80 °C freezer for biochemical assays and microRNA analysis. All procedures were performed according to protocols approved by the Institutional Animal Care and Use Committee (IACUC) of Lilly animal care and use ethical committees. A detailed description of diet treatment, histopathological analysis and biochemical assays were provided in an independent manuscript under preparation.

Circulating microRNA detection

Total RNA was isolated from 50 μL of mouse plasma using miRNeasy serum/plasma kit (QIAGEN, 217184) following manufacturer’s instructions. MicroRNAs were measured by RT-qPCR method using customized IDEAL microRNA assay kits provided by MiRXES Singapore. Briefly, 5 μL of spike-in control microRNA was added into each RNA isolation reaction. The RNA quality and concentration were checked by NanoDrop 8000, and 135 ng total RNA were converted to cDNA using microRNA-specific primers in 40 μL reaction. The cDNA was further diluted 11 times, and 5 μL diluted cDNA were used in subsequent qPCR assays for each microRNA. The PCR reaction was run using ABI7900HT thermal cycler with the following settings: 95 °C-10 min; 40 °C-5 min; 95 °C-10 s, 60 °C-30 s, repeated for 40 cycles. The qPCR data were expressed as minus delta Ct using spike-in control as reference gene.

Bioinformatics and statistical analysis

The qPCR data were expressed as minus delta Ct value using spike-in control as reference gene. The qPCR data were analyzed by hierarchical clustering, principle component analysis, random forest algorithms using MetaboAnalyst 3.0 web based data analysis suite without sample normalization, data transformation or data scaling [22]. Univariate/multivariate receiver operating characteristics (ROC) curve analysis were performed using MetaboAnalyst 3.0A as well, in which ALT concentrations were log transformed before analysis. All association studies were conducted using R-studio software.

Results

Circulating microRNA profiling revealed differential expression pattern between lean and NASH mice

The present study assessed the change in circulating microRNA expression in C57BL/6 mice fed on high fat, high cholesterol and high fructose diet (3H diet) in a time-course study. Mice on chow or 3H diet were terminated at 2 months (3H_2m, n = 5), 4 months (3H_4m, n = 5) or 7 months (3H_7m, n = 10; Lean, n = 10) for disease model characterization. Biochemical and histological studies showed that these mice developed obesity, dyslipidemia, hyperinsulinemia and liver injury in a time-dependent manner. Progressive hepatic steatosis and perisinusoidal fibrosis with concomitant elevations on ALT and AST were observed in 3H mice, which recapitulated the histological changes seen in human with NASH (The manuscript on detailed animal model characterization is in preparation.). Total 25 microRNAs were selected for profiling based on literature review of human studies (Complete list of microRNA was included in Additional file 1: Table S1). The microRNA was measured using ultrasensitive RT-qPCR assays which include three mature microRNA-specific primers in the detection system, and data were expressed as minus delta Ct with reference to spike-in control microRNA. Hierarchical clustering analysis revealed that there was distinct expression pattern of microRNAs between lean and 3H mice (Fig. 1a). Lean mice had low expression of microRNAs and were clustered together. Majority of 3H mice had high expression of microRNAs and were clustered together. In addition, a higher out of bag error rate was observed in mice on 3H diet for 2 and 4 months by random forest analysis, indicating that there is high degree of heterogeneity in NASH disease progression (Fig. 1b). It was also noted that the misclassified animal 3H_7m_9 had similar histopathological scores compared to other animals in 3H_7m group but with low expression of all microRNAs, which may imply for technical errors during microRNA isolation or detection. These data reveals distinct expression of circulating microRNAs between 3H diet treated and lean mice.

Fig. 1
figure 1

Circulating microRNA expression pattern in mice fed on 3H diet for varying time points compared with chow diet control mice. a Hierarchical clustering analysis revealed differential expression pattern between chow and 3H diet fed mice in circulating microRNA profile. The distance between features/samples were calculated by Euclidean method, and clustered using the Ward algorithm. b Random forest analysis on the microRNA expression data demonstrated a greater out-of-bag error in mice fed 3H diet for 2 and 4 months, suggesting heterogeneity in disease progression. Total number of trees were 500, and number of predictors were 7. Abbreviations: 3H, high fat, high fructose and high cholesterol; 3H_2m, 3H_4m, 3H_7m, mice on 3H diet for 2 months, 4 months, and 7 months respectively; lean, mice on standard chow diet; class. Error, classification error

A microRNA signature was identified and it accurately distinguished NASH mice from lean mice

A second circulating microRNA profiling study was performed in an independent cohort to validate the data obtained from study one. Lean mice and NASH mice (on 3H diet for 7 months) were clearly separated on the heat map by hierarchical clustering analysis (Additional file 2: Figure S1). To further understand which microRNA(s) contribute most to the classification of NASH mice, the relative importance of individual microRNA was ranked by their contributions to classification accuracy in random forest analysis for both studies (Fig. 2a). The top six common microRNAs including miR-21, miR-122, miR-192, miR-29a, miR-34a and miR-505, were designated as the circulating microRNA signature for NASH (Fig. 2b). These microRNAs were significantly up-regulated in NASH mice in both studies, and their expression was correlated with the duration of 3H diet treatment (Additional file 3: Figure S2). Dimensional reduction by principle component analysis revealed that lean mice and NASH mice were distantly related based on their microRNA signature. As shown in Fig. 2c, the two groups of animals were clearly separated on the 3D-PCA plot, in which 96.2% of variance could be explained between NASH and lean mice. This set of data suggests that microRNA signature has excellent prediction power in distinguishing NASH mice, although the effect of sample homogeneity needs to be considered. Since these microRNAs were validated in human studies, our data established the translation of circulating microRNA between mice and men in NASH disease evolution and progression.

Fig. 2
figure 2

A panel of 6 microRNAs was designated as circulating microRNA signature for NASH. a Top seven microRNAs were listed and ranked by their classification accuracy in two independent studies by random forest analysis. b The expression level of the microRNA signature in lean and NASH mice in study 2. Data were expressed as -ΔCt (minus delta Ct) with reference to spike-in control microRNA. c Principle component analysis (PCA) of miR-192, miR-122, miR-21, miR-29a, miR-34a, and miR-505 in study 2. Red dots represented NASH mice (mice on 3H diet for 7 months), and green dots represented lean mice. Numbers in parenthesis mean the percentage of variance explained by each principle component

The microRNA signature responded to chow diet treatment in NASH regression model

We then assessed whether the circulating microRNA signature could be reversed in a disease regression model. To test our hypothesis, the expression of microRNA signature was examined in NASH regression model, in which established NASH mice (with 7 months of 3H diet treatment) either continued on 3H diet or switched to chow diet for 42 and 84 days before termination for analysis. Three dimensional PCA plot revealed clear separation of lean mice from 3H mice groups (3H_3H_42 and 3H_3H_84). More importantly, 94.4% of mice switched to chow diet (3H_chow_42/84) were clustered with lean mice (Fig. 3a), indicating that the microRNA signature in diet-switching groups was normalized to an extent which is indistinguishable from that of lean mice. To further understand whether decrease in microRNA signature expression predicts disease regression, scatter plots on miR-192 expression against histological scores were generated (Fig. 3b). It was clear that diet-switching group animals had lower level of circulating miR-192 accompanied with improved liver histological lesions. Similar results were obtained for other signature microRNAs (data not shown). Altogether, these data independently validated the value of microRNA signature as diagnostic biomarkers for NASH, and further suggests the potential application in assessing treatment effects.

Fig. 3
figure 3

The microRNA signature responded to chow diet treatment in NASH mice. Mice were on 3H diet for 7 months for NASH development. Half of mice were switched to chow diet and were sacrificed on day 42 (3H_chow_42) or day 84 (3H_chow_84); half of mice were continued on 3H diet and were sacrificed on day 42 (3H_3H_42) or day 84 (3H_3H_84) as control. Lean mice were on standard chow diet and sacrificed together with the 3H_3H_84 group. a Principle component analysis of the microRNA signature in separating diet-switching animals (3H_chow_42, 3H_chow_84) from 3H mice (3H_3H_42, 3H_3H_84). b Scatter plots of miR-192 expression (reference to spike-in control) against histopathological scores: inflammation, macrovesicular vacuolation, perisinusoidal fibrosis and portal fibrosis

Inclusion of ALT improved the diagnostic power of microRNA signature

Clinical diagnosis of NASH requires the histological examination of liver biopsies [23]. NAFLD activity score was developed by the pathology subcommittee of the NASH Clinical Research Network, and was validated by a group of expert hepatopathologists. In clinical practice, patients with NAS score of ≥5 strongly correlated with a diagnosis of “definite NASH”, whereas NAS ≤3 correlated with a diagnosis of “not NASH” [24]. We then assessed the performance of microRNA signature in discriminating mice which had NAS above 3 in study 1 because this cohort of mice had heterogeneous disease pathology. Univariate ROC curve analysis showed that miR-192 and miR-505 achieved the greatest AUROC of 0.923 and 0.919 respectively in discriminating mice had NAS > 3, while miR-122, miR-29a, miR-34a and miR-21 had an AUROC of 0.88, 0.84, 0.80 and 0.79 (Fig. 4a and b). To avoid potential overfitting of the data, multivariate ROC analysis was performed, and incorporation of all six microRNAs achieved the highest prediction power (AUROC = 0.897, 95% confidence interval: 0.75- 1) (Fig. 4b).

Fig. 4
figure 4

The diagnostic performance of microRNA signature in discriminating mice having NAS > 3 was improved by incorporating ALT. a, b Univariate ROC analysis of miR-192 and miR-505, which was created by plotting the true positive rate (sensitivity) against false positive rate (1-specificity). The right panel depicted the expression level of miR-192 and miR-505 in mice having NAS > 3 compared with mice having NAS ≤ 3. c ROC curves of the microRNA signature in predicting mice having NAS > 3 (AUC = 0.897, 95% confidence interval: 0.75 - 1) in study 1. The confusion matrix was depicted in the inset (numbers in rows are actual classification, numbers in columns are predicted classification). d ROC curves of the microRNA signature plus ALT in predicting mice having NAS > 3 (AUC = 0.913, 95% confidence interval: 0.75 - 1) in study 1. The confusion matrix was depicted in the inset. Abbreviations: NAS, NAFLD activity score; ROC, receiver operating characteristics; AUC, area under curve

Liver transaminases (ALT and AST) have been widely used in clinical settings to predict the presence of liver diseases including NASH [25]. Elevations in circulating ALT level are indicative of hepatocyte damage, and ALT is an important component of several composite NASH scoring systems [26, 27]. As ALT and microRNAs are independent measurements of liver functions, we then sought to incorporate ALT into the microRNA signature to see if the performance could be further improved. As shown in Fig. 4c, the AUROC reached 0.913, and there was a reduction of 2 misclassified animals compared with microRNA signature model. These data demonstrated the feasibility of developing reliable circulating microRNA and ALT-based biomarkers to identify patients “at risk” for NASH thereby decreasing the need of liver biopsy during patient screening.

Development of composite biomarkers for non-invasive diagnosis and prognosis of NASH

A robust non-invasive biomarker for NASH should be sensitive in disease discrimination, predictive of disease severity, and inexpensive for detection. The microRNA signature consists of 6 microRNAs, which may reduce the robustness and pose hurdles in future development. Thus, we sought to reduce the number of microRNAs by studying the associations among the microRNAs, and between microRNAs, and disease pathologies (Fig. 5a). MiR-192 and miR-122 fell into the same category, and both were implicated in liver injury and hepatocyte death [28, 29]. Since miR-192 had the highest association with other microRNAs and disease pathologies (NAS, r = 0.78), it was chosen to build the new composite biomarker. MiR-505 was selected because of its high correlations with disease pathologies (NAS, r = 0.77), but relatively low associations with other microRNAs. MiR-21 and miR-29a represented another cluster of microRNA, and both were functionally related to hepatic stellate cell activation and liver fibrosis [30,31,32,33]. We included miR-21 in the new biomarker because of its stronger human validation [13]. MiR-34a was weakly associated with disease pathologies or other microRNAs, thus it was omitted for further analysis. We then tested the performance of the new composite biomarker consisting of miR-192, miR-505, miR-21 and ALT in discriminating NASH animals (NAS > 3) from healthy mice (NAS ≤ 3) in study 1 and 3. As shown in Fig. 5b, the new biomarker outperformed the microRNA signature plus ALT model, and achieved AUROC of 0.958 with further reduction of one misclassified animal in study one. The new composite biomarker also accomplished excellent prediction accuracy in study three (AUROC of 0.932). Collectively, our data strongly validated the new composite biomarker in NASH diagnosis, and the performance of the new biomarker warrants further evaluation in large scale studies.

Fig. 5
figure 5

Development of non-invasive biomarker for NASH. a Associations between microRNAs and ALT was calculated. ALT concentration was log transformed before analysis. b, c ROC curves of the miR-192, miR-21, miR-505 and ALT in predicting mice having NAS > 3 in study 1 (AUC = 0.948, 95% confidence interval: 0.827 - 1) and study 3 (AUC = 0.931, 95% confidence interval: 0.845 - 1). The confusion matrix of study 1 was depicted in the inset

Discussion

The present study evaluated the diagnostic and prognostic value of circulating microRNA in a well-studied diet-induced NASH mouse model. Expression of 25 selected microRNA was measured using a unique and ultrasensitive RT-qPCR based platform. A panel of 6 microRNA was validated to be differentially expressed in NASH mice and was correlated with disease severity. Based on these findings, a microRNA-based composite biomarker consisting of miR-192, miR-21, miR-505 and ALT was proposed, which demonstrated great performance in distinguishing between lean and NASH mice (NAS > 3). The composite biomarker mirrors the histological and molecular events in the liver, and potentially have both preclinical and clinical applications. In preclinical drug discovery, the composite biomarker could be used to evaluate disease progression and to assess treatment effects in mouse models. In clinical settings, the composite biomarker could be used for screening patients at risk for NASH, and/or used for monitoring treatment response, thus reducing the use of liver biopsies.

The microRNA based biomarkers demonstrated robust sensitivity in NASH diagnosis. In view of the biological function of each feature, the new composite biomarker could potentially reflect the status of the liver from the following perspectives: miR-192 and ALT serve as independent indicator of hepatocyte function, miR-21 suggests of stellate cell activation and liver fibrosis, miR-505 implies for pathological manifestations. MiR-192 regulates lipid synthesis in hepatocytes, and it was independently validated as biomarker for NASH by several groups using in vitro assays, animal models and clinical samples [12, 13, 34,35,36]. MiR-21 as biomarkers for hepatic fibrosis has been validated in both animal models and clinical samples [13, 37]. Studies have also shown that miR-21 regulates the peroxisome proliferator-activated receptor α (PPARα) activity and associates with stellate cell activation [31, 38]. Considering the association of miR-505 with various types of cancers, the biological function of miR-505 in NASH and liver carcinogenesis deserves further investigation [39]. Since these microRNAs originated from liver, their expression in NASH liver should also be considered. Previous studies have shown that liver miR-122 preferentially expressed in lipid-laden hepatocytes and was down-regulated in NASH patients [12]. Similarly, we observed slightly decrease in the liver miR-122 and miR-192 expression (not statistically significant) in NASH mice compared to lean mice. As miR-122 was very abundant in hepatocytes, a small number of hepatocyte degeneration may lead to significant release of miR-122 into circulation. Notably, we detected significant up-regulation of liver miR-21, miR-34a and miR-505 in NASH mice compared to lean mice (data was not shown). Since these microRNAs have been implicated in hepatic lipid homeostasis, inflammation, stellate cell activation and carcinogenesis [9, 37, 40, 41]; our results strongly support that the circulating level of microRNA could reflect the status of the liver, and hence could be employed as non-invasive diagnostic tools for NASH.

MicroRNAs are epigenetic regulators, which post-transcriptionally modulate hundreds of target genes. The experimental validated target genes of the microRNA signature were identified by miRTarBase [42]. Ingenuity Pathway Analysis of the target genes reveal that the top network was endocrine disorders, gastrointestinal disease, and metabolic disorders. It also suggests that EGFR signaling pathway may be a potential target for liver regeneration following chronic liver damage in NASH. Indeed, Erlotinib treatment in carbon tetrachloride-induced mouse fibrosis model attenuated liver fibrosis and the development of hepatocellular carcinoma [43]. These findings suggest that functional microRNA – mRNA network analysis may represent a novel approach in the discovery of new drug target for the treatment of NASH.

Considering future directions, we think that development of reliable internal normalization strategy is critical. We used spike-in microRNA control as internal reference, however, it might not be ideal for clinical use as the relative microRNA expression might be influenced by batch-to-batch variance resulting from RNA isolation. Therefore we propose that identification of house-keeping microRNAs or generation of statistical internal normalization algorithms are urgently needed for the absolute quantification of circulating microRNAs in NASH patients for accurate diagnosis and staging of the disease.

MicroRNAs in serum or plasma are relatively stable as they are associated with membrane-bound vesicles or proteins. The ease of sample handling and the availability of inexpensive qPCR-assays make them ideal for use in clinic settings for patient screening and stratification. Although our data are promising, the main limitation was that the current study focused on a selected panel of microRNAs. Future effort should be devoted to global screening studies to identify microRNAs with the best diagnostic and prognostic performance specific for NASH.

Conclusion

The present circulating microRNA profiling studies assessed the expression of 25 selected microRNAs in a diet-induced NASH mouse model. Our data validated the differential expression of a panel of microRNAs in NASH models, and supported the translation between preclinical and clinical studies. The circulating microRNA based biomarker proposed here holds great promise for clinical application, and deserves further evaluation in clinical studies.

Abbreviations

ALT:

Alanine transaminase

AST:

Asparte transaminase

AUROC:

Area under ROC

MicroRNA:

microRNA

NAFLD:

Nonalcoholic fatty liver disease

NAS:

NAFLD Activity score

NASH:

Nonalcoholic steatohepatitis

PCA:

Principle component analysis

qPCR:

Quantitative polymerase chain reaction

ROC:

Receiver operating characteristic

References

  1. Neuschwander-Tetri BA, Clark JM, Bass NM, Van Natta ML, Unalp-Arida A, Tonascia J, Zein CO, Brunt EM, Kleiner DE, McCullough AJ, et al. Clinical, laboratory and histological associations in adults with nonalcoholic fatty liver disease. Hepatology. 2010;52(3):913–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Anstee QM, Targher G, Day CP. Progression of NAFLD to diabetes mellitus, cardiovascular disease or cirrhosis. Nat Rev Gastroenterol Hepatol. 2013;10(6):330–44.

    Article  CAS  PubMed  Google Scholar 

  3. Sanyal AJ, Friedman SL, McCullough AJ, Dimick-Santos L, American Association for the Study of Liver D, United States F, Drug A. Challenges and opportunities in drug and biomarker development for nonalcoholic steatohepatitis: findings and recommendations from an American Association for the Study of Liver Diseases-U.S. Food and Drug Administration joint workshop. Hepatology. 2015;61(4):1392–405.

    Article  PubMed  PubMed Central  Google Scholar 

  4. Ratziu V, Charlotte F, Heurtier A, Gombert S, Giral P, Bruckert E, Grimaldi A, Capron F, Poynard T, Group LS. Sampling variability of liver biopsy in nonalcoholic fatty liver disease. Gastroenterology. 2005;128(7):1898–906.

    Article  PubMed  Google Scholar 

  5. Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet. 2007;8(2):93–103.

    Article  CAS  PubMed  Google Scholar 

  6. Momen-Heravi F, Saha B, Kodys K, Catalano D, Satishchandran A, Szabo G. Increased number of circulating exosomes and their microRNA cargos are potential novel biomarkers in alcoholic hepatitis. J Transl Med. 2015;13:261.

    Article  PubMed  PubMed Central  Google Scholar 

  7. Baldassarre A, Felli C, Prantera G, Masotti A. Circulating microRNAs and bioinformatics tools to discover novel diagnostic biomarkers of pediatric diseases. Genes (Basel). 2017;8(9):234.

    Article  Google Scholar 

  8. Esteller M. Non-coding RNAs in human disease. Nat Rev Genet. 2011;12(12):861–74.

    Article  CAS  PubMed  Google Scholar 

  9. Rottiers V, Naar AM. MicroRNAs in metabolism and metabolic disorders. Nat Rev Mol Cell Biol. 2012;13(4):239–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  10. Deiuliis JA. MicroRNAs as regulators of metabolic disease: pathophysiologic significance and emerging role as biomarkers and therapeutics. Int J Obes. 2016;40(1):88–101.

    Article  CAS  Google Scholar 

  11. Roderburg C, Luedde T. Circulating microRNAs as markers of liver inflammation, fibrosis and cancer. J Hepatol. 2014;61(6):1434–7.

    Article  PubMed  Google Scholar 

  12. Pirola CJ, Fernandez Gianotti T, Castano GO, Mallardi P, San Martino J, Mora Gonzalez Lopez Ledesma M, Flichman D, Mirshahi F, Sanyal AJ, Sookoian S. Circulating microRNA signature in non-alcoholic fatty liver disease: from serum non-coding RNAs to liver histology and disease pathogenesis. Gut. 2015;64(5):800–12.

    Article  CAS  PubMed  Google Scholar 

  13. Becker PP, Rau M, Schmitt J, Malsch C, Hammer C, Bantel H, Mullhaupt B, Geier A. Performance of serum microRNAs −122, −192 and −21 as biomarkers in patients with non-alcoholic steatohepatitis. PLoS One. 2015;10(11):e0142661.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Chen YJ, Zhu JM, Wu H, Fan J, Zhou J, Hu J, Yu Q, Liu TT, Yang L, Wu CL, et al. Circulating microRNAs as a fingerprint for liver cirrhosis. PLoS One. 2013;8(6):e66577.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Celikbilek M, Baskol M, Taheri S, Deniz K, Dogan S, Zararsiz G, Gursoy S, Guven K, Ozbakir O, Dundar M, et al. Circulating microRNAs in patients with non-alcoholic fatty liver disease. World J Hepatol. 2014;6(8):613–20.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Pogribny IP. MicroRNAs as biomarkers for clinical studies. Exp Biol Med (Maywood). 2018;243(3):283–90.

    Article  CAS  Google Scholar 

  17. Sarkar N. Circulating microRNAs as markers for liver histological lesions in patients with nonalcoholic steatohepatitis. AASLD LiverLearning. 2016;143957. Abstract ID:1063.

  18. Francque S. LBP-535 - next-generation sequencing (NGS) of two independent cohorts identifies eleven circulating miRNAs for diagnosis of NASH and fibrosis. J Hepatol. 2017;66(1):S110–1.

    Article  Google Scholar 

  19. DiStefano JK, Gerhard GS. Circulating microRNAs in nonalcoholic fatty liver disease. Expert Rev Gastroenterol Hepatol. 2016;10(2):161–3.

    Article  CAS  PubMed  Google Scholar 

  20. Yamada H, Suzuki K, Ichino N, Ando Y, Sawada A, Osakabe K, Sugimoto K, Ohashi K, Teradaira R, Inoue T, et al. Associations between circulating microRNAs (miR-21, miR-34a, miR-122 and miR-451) and non-alcoholic fatty liver. Clin Chim Acta. 2013;424:99–103.

    Article  CAS  PubMed  Google Scholar 

  21. Clapper JR, Hendricks MD, Gu G, Wittmer C, Dolman CS, Herich J, Athanacio J, Villescaz C, Ghosh SS, Heilig JS, et al. Diet-induced mouse model of fatty liver disease and nonalcoholic steatohepatitis reflecting clinical disease progression and methods of assessment. Am J Physiol Gastrointest Liver Physiol. 2013;305(7):G483–95.

    Article  CAS  PubMed  Google Scholar 

  22. Xia J, Wishart DS. Using MetaboAnalyst 3.0 for comprehensive metabolomics data analysis. Curr Protoc Bioinformatics. 2016;55:14 10 11-14 10 91.

    Google Scholar 

  23. Ludwig J, Viggiano TR, McGill DB, Oh BJ. Nonalcoholic steatohepatitis: Mayo Clinic experiences with a hitherto unnamed disease. Mayo Clin Proc. 1980;55(7):434–8.

    CAS  PubMed  Google Scholar 

  24. Kleiner DE, Brunt EM, Van Natta M, Behling C, Contos MJ, Cummings OW, Ferrell LD, Liu YC, Torbenson MS, Unalp-Arida A, et al. Design and validation of a histological scoring system for nonalcoholic fatty liver disease. Hepatology. 2005;41(6):1313–21.

    Article  PubMed  Google Scholar 

  25. Xia MF, Lin HD, Yan HM, Bian H, Chang XX, Zhang LS, He WY, Gao X. The association of liver fat content and serum alanine aminotransferase with bone mineral density in middle-aged and elderly Chinese men and postmenopausal women. J Transl Med. 2016;14:11.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Fallatah HI, Akbar HO, Fallatah AM. Fibroscan compared to FIB-4, APRI, and AST/ALT ratio for assessment of liver fibrosis in Saudi patients with nonalcoholic fatty liver disease. Hepat Mon. 2016;16(7):e38346.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Rath MM, Panigrahi MK, Pattnaik K, Bhuyan P, Kar SK, Misra B, Misra D, Meher C, Agrawal O, Rath J, et al. Histological evaluation of non-alcoholic fatty liver disease and its correlation with different noninvasive scoring systems with special reference to fibrosis: a single center experience. J Clin Exp Hepatol. 2016;6(4):291–6.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Roy S, Benz F, Alder J, Bantel H, Janssen J, Vucur M, Gautheron J, Schneider A, Schuller F, Loosen S, et al. Down-regulation of miR-192-5p protects from oxidative stress-induced acute liver injury. Clin Sci. 2016;130(14):1197–207.

    Article  CAS  PubMed  Google Scholar 

  29. Roderburg C, Benz F, Vargas Cardenas D, Koch A, Janssen J, Vucur M, Gautheron J, Schneider AT, Koppe C, Kreggenwinkel K, et al. Elevated miR-122 serum levels are an independent marker of liver injury in inflammatory diseases. Liver Int. 2015;35(4):1172–84.

    Article  CAS  PubMed  Google Scholar 

  30. Wu K, Ye C, Lin L, Chu Y, Ji M, Dai W, Zeng X, Lin Y. Inhibiting miR-21 attenuates experimental hepatic fibrosis by suppressing both the ERK1 pathway in HSC and hepatocyte EMT. Clin Sci. 2016;130(16):1469–80.

    Article  CAS  PubMed  Google Scholar 

  31. Wei J, Feng L, Li Z, Xu G, Fan X. MicroRNA-21 activates hepatic stellate cells via PTEN/Akt signaling. Biomed Pharmacother. 2013;67(5):387–92.

    Article  CAS  PubMed  Google Scholar 

  32. Knabel MK, Ramachandran K, Karhadkar S, Hwang HW, Creamer TJ, Chivukula RR, Sheikh F, Clark KR, Torbenson M, Montgomery RA, et al. Systemic delivery of scAAV8-encoded MiR-29a ameliorates hepatic fibrosis in carbon tetrachloride-treated mice. PLoS One. 2015;10(4):e0124411.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Matsumoto Y, Itami S, Kuroda M, Yoshizato K, Kawada N, Murakami Y. MiR-29a assists in preventing the activation of human stellate cells and promotes recovery from liver fibrosis in mice. Mol Ther. 2016;24(10):1848–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Di Mauro S, Ragusa M, Urbano F, Filippello A, Di Pino A, Scamporrino A, Pulvirenti A, Ferro A, Rabuazzo AM, Purrello M, et al. Intracellular and extracellular miRNome deregulation in cellular models of NAFLD or NASH: clinical implications. Nutr Metab Cardiovasc Dis. 2016;26(12):1129–39.

    Article  CAS  PubMed  Google Scholar 

  35. Povero D, Eguchi A, Li H, Johnson CD, Papouchado BG, Wree A, Messer K, Feldstein AE. Circulating extracellular vesicles with specific proteome and liver microRNAs are potential biomarkers for liver injury in experimental fatty liver disease. PLoS One. 2014;9(12):e113651.

    Article  PubMed  PubMed Central  Google Scholar 

  36. Liu XL, Cao HX, Wang BC, Xin FZ, Zhang RN, Zhou D, Yang RX, Zhao ZH, Pan Q, Fan JG. miR-192-5p regulates lipid synthesis in non-alcoholic fatty liver disease through SCD-1. World J Gastroenterol. 2017;23(46):8140–51.

    Article  PubMed  PubMed Central  Google Scholar 

  37. Takeuchi-Yorimoto A, Yamaura Y, Kanki M, Ide T, Nakata A, Noto T, Matsumoto M. MicroRNA-21 is associated with fibrosis in a rat model of nonalcoholic steatohepatitis and serves as a plasma biomarker for fibrotic liver disease. Toxicol Lett. 2016;258:159–67.

    Article  CAS  PubMed  Google Scholar 

  38. Rodrigues PM, Afonso MB, Simao AL, Carvalho CC, Trindade A, Duarte A, Borralho PM, Machado MV, Cortez-Pinto H, Rodrigues CM, et al. miR-21 ablation and obeticholic acid ameliorate nonalcoholic steatohepatitis in mice. Cell Death Dis. 2017;8(4):e2748.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Chen S, Sun KX, Liu BL, Zong ZH, Zhao Y. MicroRNA-505 functions as a tumor suppressor in endometrial cancer by targeting TGF-alpha. Mol Cancer. 2016;15:11.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Loyer X, Paradis V, Henique C, Vion AC, Colnot N, Guerin CL, Devue C, On S, Scetbun J, Romain M, et al. Liver microRNA-21 is overexpressed in non-alcoholic steatohepatitis and contributes to the disease in experimental models by inhibiting PPARalpha expression. Gut. 2016;65(11):1882–94.

    Article  CAS  PubMed  Google Scholar 

  41. Lu L, Qiu C, Li D, Bai G, Liang J, Yang Q. MicroRNA-505 suppresses proliferation and invasion in hepatoma cells by directly targeting high-mobility group box 1. Life Sci. 2016;157:12–8.

    Article  CAS  PubMed  Google Scholar 

  42. Hsu SD, Lin FM, Wu WY, Liang C, Huang WC, Chan WL, Tsai WT, Chen GZ, Lee CJ, Chiu CM, et al. miRTarBase: a database curates experimentally validated microRNA-target interactions. Nucleic Acids Res. 2011;39(Database issue):D163–9.

    Article  CAS  PubMed  Google Scholar 

  43. Fuchs BC, Hoshida Y, Fujii T, Wei L, Yamada S, Lauwers GY, McGinn CM, DePeralta DK, Chen X, Kuroda T, et al. Epidermal growth factor receptor inhibition attenuates liver fibrosis and development of hepatocellular carcinoma. Hepatology. 2014;59(4):1577–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

Not applicable.

Availability of data and materials

The datasets for the current study are available from the corresponding author on reasonable request.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

JL, HH, ZD, MK, and JS conceived the idea and designed the study; JL, YX, and SY conducted the experiments; JL, XW, LJ, ZF performed the bioinformatic and statistical analysis;JL, HH, ZD, MK, and JS analyzed the data; all authors participated in the manuscript writing; JL, HH and JS revised the manuscript. All authors reviewed and approved the manuscript for publication.

Corresponding authors

Correspondence to Jie Liu or Jianyong Shou.

Ethics declarations

Ethics approval

All protocols related to animal studies were approved by the Institutional Animal Care and Use Committee (IACUC) of Lilly animal care and use ethical committees.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

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

Additional files

Additional file 1:

Table S1. Complete list of microRNAs analyzed. (DOCX 17 kb)

Additional file 2:

Figure S1. Hierarchical clustering analysis revealed differential expression of circulating microRNAs between lean and 3H mice in an independent study (study 2). (DOCX 347 kb)

Additional file 3:

Figure S2. Time-dependent expression of miR-122, miR-192, miR-21, miR-29a, miR-34a, and miR-505 in study 1. Data were expressed as minus delta Ct with reference to spike-in control miRNA. (DOCX 149 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Liu, J., Xiao, Y., Wu, X. et al. A circulating microRNA signature as noninvasive diagnostic and prognostic biomarkers for nonalcoholic steatohepatitis. BMC Genomics 19, 188 (2018). https://doi.org/10.1186/s12864-018-4575-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-018-4575-3

Keywords