Introduction
Nonalcoholic fatty liver disease (NAFLD) is defined as the pathological accumulation of lipids in hepatocytes without significant alcohol consumption, use of steatogenic medication or hereditary disorders [1]. The pathologic classification of NAFLD encompasses isolated hepatic steatosis, nonalcoholic steatohepatitis (NASH) and cryptogenic fibrosis/cirrhosis [2]. As the global obesity epidemic grows, the incidence of NAFLD is increasing; in morbidly obese adults, the prevalence of NAFLD is more than 90% [3], and the prevalence is thought to reach up to 50% in obese children [4]. Unfortunately, patients with NAFLD show increased rates of overall mortality compared to matched control populations [5,6]. Although much progress has been made towards gaining a better understanding of the pathogenesis of NAFLD, the exact molecular mechanisms remain unclear.
MicroRNAs (miRNAs), which are 19-25-nucleotide (nt) non-coding small RNAs, modulate gene expression by binding to the 3' untranslated region (3'UTR) of target messenger RNAs (mRNAs). MiRNAs are thought to be involved in various biological and pathological processes, including the growth and development of organisms, cancer and other diseases [7,8,9,10]. Recent research has demonstrated that miRNAs are able to regulate metabolism-related genes that may be involved in the pathogenesis of NAFLD. For example, miR-122 is the most highly expressed miRNA in the liver and can regulate cholesterol and fatty acid metabolism [11,12], and miR-10b has been shown to regulate the level of cellular steatosis through peroxisome proliferator-activated receptor-alpha [13]. In addition, several other miRNAs have been shown to regulate insulin resistance [14], mitochondrial dysfunction [15], oxidative stress [16], and inflammation [17], which are all generally regarded as critical components of NAFLD.
Therefore, we hypothesised that certain miRNAs may regulate steatosis in NAFLD. In the present study, we employed the Illumina HiSeq deep sequencing platform to screen the miRNA expression profiles in high-fat diet (HFD)-induced fatty liver tissue and control tissue from animals fed a standard diet. From our sequencing results, we identified differentially expressed miRNAs between the normal control group and the model group. We then investigated the expression of certain miRNAs in HepG2 cells and human hepatocytes cultured with a high concentration of free fatty acids (FFAs) with or without proinflammatory factors. We found that miR-200a, miR-200b, miR-200c, miR-146a, miR-146b and miR-152 were up-regulated both in fatty liver tissues and in cells cultured with FFAs and proinflammatory factors.
Materials and Methods
Animal experiments
This study was approved by the University Committee for the Ethical Care and Use of Animals in Experiments of ZheJiang University (Permit Number: 20100022). Forty male Sprague-Dawley rats were obtained from the Experimental Animal Center of Nanjing Medical University (Nanjing, China). The animals were divided randomly into the following four groups, with each group matched for body weight and age: (1) normal control group A (n=10) and (2) normal control group B (n=10), which were fed a standard diet for 4 and 12 weeks, respectively; and (3) model group C (n=10) and (4) model group D (n=10), which were fed a HFD (standard diet + 10% lard + 2% cholesterol) for 4 and 12 weeks, respectively. All of the animals were given ad libitum access to food and purified water and were housed in sterilised cages with a 12-h light/dark cycle. At the end of week 4 or 12, after one night of fasting, the rats were weighed, and blood samples were collected through the inferior vena cava under general anaesthesia using 10% chloral hydrate inhalation. Intact livers were removed from the abdominal cavity and washed with ice-cold 0.9% saline. The weights of the livers were measured, and the hepatic index was calculated (wet weight of liver/body weight). The right lobes of the livers were dissected, and pieces were frozen immediately in liquid nitrogen and stored at -80°C for further analysis. The left lobes were soaked in 4% paraformaldehyde for histological observation by haematoxylin and eosin (H&E) staining.
The diagnosis of NASH and classification of severity were made according to the 2008 criteria of the Chinese Medical Association Committee of Nonalcoholic Fatty Liver Disease [18].
Serum biochemical analysis
The levels of serum alanine aminotransferase (ALT), aspartate aminotransferase (AST), total cholesterol (TC) and triglycerides (TG) were determined using a fully automatic biochemical analyser (Roche Inc., USA). Hepatic TG levels were assayed using a TG assay kit (Applygen Technologies Inc., Beijing, China).
Deep sequencing
Total RNA from the livers of normal control group B and model group D rats was obtained using the Trizol (Invitrogen) reagent. Subsequently, small RNAs ranging from 18 to 30 nt were isolated using denaturing polyacrylamide gel electrophoresis. We then ligated synthetic oligonucleotide adapters of known sequences to the 3' and 5' ends of the small RNA pool using T4 RNA ligases, which can introduce primer-binding sites for reverse transcription and PCR amplification. Then, the modified small RNAs were reversed transcribed into cDNA and amplified by PCR. Thereafter, purified samples were obtained for sequencing-by-synthesis (Solexa). Mireap software (http://sourceforge.net/projects/mireap/) was used to identify both known and novel microRNAs from small RNA libraries that were deep-sequenced using Solexa technology. Finally, differentially expressed microRNAs with statistical significance were identified through volcano plot filtering. The threshold used to screen up-regulated and down-regulated microRNAs was a 1.5-fold increase and a 0.5-fold decrease, respectively (P<0.05).
Cell culture
Human hepatocellular carcinoma cells (HepG2) were purchased from ATCC (Rockville, MD, USA), and human hepatocytes were obtained from Sciencell. HepG2 cells were cultured in Dulbecco's modified Eagle's medium (DMEM) (Gibco) supplemented with 10% heat-inactivated foetal bovine serum (FBS), 100 U/ml penicillin, and 100 µg/ml streptomycin at 37°C and 5% CO2. Human hepatocytes were cultured in hepatocyte medium (Sciencell) containing 5% FBS, 10 µg/ml bovine serum albumin (BSA), 10 µg/ml apotransferrin, 5 µg/ml insulin, 2 ng/ml epidermal growth factor (EGF), 2 ng/ml fibroblast growth factor-2 (FGF-2), 2 ng/ml insulin-like growth factor-1 (IGF-1), 0.7 nM retinoic acid, and 1 µg/ml hydrocortisone. The cultures were used at 60-80% confluency.
Lipid accumulation and Oil Red O stain
To induce fat overloading of cells, human hepatocytes and HepG2 cells at 80% confluency were cultured for 24 h in complete medium containing 1 mM long-chain mixed FFAs (oleate acid and palmitate (Sigma-Aldrich) mixed at a final ratio of 2:1 and bound to 1% defatted BSA (Equitech-Bio) [19]. To detect lipid droplet accumulation, cells were stained with Oil Red O. Briefly, cells were rinsed three times in PBS, fixed in 4% paraformaldehyde, rinsed three times in PBS, stained with 0.3% Oil Red O for 30 min at 37°C, rinsed with PBS and photographed. Then integrated option density (IOD) determined by Image Pro-Plus according to the staining intensity by visual assessment was utilized to compare the lipid droplet accumulation. To quantify lipid accumulation, cultured cells with or without 1mM FFAs were scraped in PBS and disrupted by cell lysis buffer. The intracellular triglycerides were assayed using a commercial kit (GPO-POD; Applygen Technologies Inc., Beijing, China) according to the manufacturer's recommended protocol.
Proinflammatory factor assay
Inflammation is an important component of NAFLD. Specifically, tumour necrosis factor-alpha (TNF-α) and interleukin-6 (IL-6) are known to be significantly increased in NAFLD patients compared to healthy controls [3,20]. To assess the effect of proinflammatory cytokines on the expression of miRNAs, we used recombinant human IL-6 and TNF-α (R&D Systems, Minneapolis, MN, USA) to treat HepG2 cells. Specifically, cells were seeded in normal growth medium at 60,000 cells/ml, with 2.0 ml/well in 6-well culture plates, and grown for 2 days. On day 3, the cells were washed with serum-free DMEM and then incubated in serum-free DMEM for an additional 24 h. Subsequently, the cells were washed and then incubated with DMEM containing 10 ng/ml TNF-α [21], 25 ng/ml IL-6 [22], or both of these compounds for 0, 4, 8 or 24 h. All of the cell experiments were performed in triplicate.
Quantitative miRNA expression analysis
Total RNA was extracted from 10 mg of liver tissue using a miRNeasy mini kit (Qiagen, Valencia, CA, USA), according to the manufacturer's instructions. Briefly, the tissue was dispersed using an electric homogeniser after adding Qiazol lysis reagent, and the cells were directly lysed using Qiazol lysis reagent. Then, chloroform was added, and the RNA phase was transferred to a column. After washing with the provided buffer, the RNA was dissolved in nuclease-free water. The RNA was quantified using a NanoDrop spectrophotometer (Thermo Scientific, USA), and the purity was evaluated by measuring the A260/280 ratio.
The reverse transcriptase (RT) primer, PCR primers, and TaqMan probe for miR-200a, miR-200b, miR-200c, miR-146a, miR-146b and miR-152 were purchased from Applied Biosystems (Foster City, CA, USA). Total RNA was reverse-transcribed using the TaqMan MicroRNA reverse transcription kit (Applied Biosystems). U6 snRNA served as an endogenous control for miRNA expression. Reverse transcription was performed in a thermal cycler using the following reaction protocol: 30 min at 16°C (annealing), 30 min at 42°C (cDNA synthesis), and 5 min at 85°C (enzyme denaturation). For real-time PCR, 10 µl of TaqMan 2× Universal PCR Master Mix II, 1 µl of cDNA, 1 µl of the TaqMan miRNA assay mix (Applied Biosystems), and 8 µl of water were mixed for 10 min at 95°C, followed by 40 cycles of 15 sec at 95°C and 1 min at 60°C. An ABI 7500 real-time PCR system (Applied Biosystems) was used for the real-time PCR experiments. Data analysis was carried out using the 2-ΔCt Method. Difference in the Ct between miRNA and U6 was called ΔCt. Data are 2-ΔCt levels calculated relative to the U6 snRNA expression.
Bioinformatic analysis
Predicted gene targets of all differentially expressed miRNAs were obtained using two computer-aided algorithms, including TargetScan (http://www.targetscan.org/) and miRanda (http://www.microrna.org/). To further understand the function, classification and metabolic regulatory networks of the miRNA targets, functional analyses were performed based on the Gene Ontology (GO) database (http://www.geneontology. org/) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway database (http://www.genome. jp/kegg/). The GO database is extensively used to analyse all types of high-throughput experiments, such as large-scale sequencing efforts [23], microarrays [24] protein-protein interactions [25], protein-DNA interactions [26], and knockouts [27]. The KEGG pathway is a collection of manually drawn pathway maps representing molecular interactions and reaction networks and has been widely used for biological interpretation of higher-level systemic functions [28].
Statistical analysis
All of the data are presented as the mean ± standard deviation. The statistical significance between two groups was compared using the Student's t-test, and one-way analysis of variance (ANOVA) was used for comparisons of more than three groups. A P value <0.05 was considered statistically significant.
Results
Analysis of the animal model
Twenty male Sprague-Dawley rats were fed a HFD (standard diet + 10% lard + 2% cholesterol) for 4 or 12 weeks. All of the rats developed NASH at the end of 12 weeks according to the criteria for NASH diagnosis. Histological examinations showed normal liver tissue or mild microvesicular steatosis without inflammatory infiltration or fibrosis in the normal group (Fig. 1A and 1C). Microvesicular steatosis, consisting predominantly of swollen, ballooned hepatocytes, with occasional foci of inflammatory cells in the hepatic lobules was observed in the 4-week model group (group C) (Fig. 1B). In the 12-week model group (group D), we observed mild to severe macro- and micro-vesicular steatosis as well as infiltration of mild lobular mixed polymorphonuclear and mononuclear cells. The 12-week group also exhibited focal necrosis, perisinusoidal and portal fibrosis, as well as severe hepatocellular ballooning (Fig. 1D). Furthermore, at the end of 4 or 12 weeks, the liver TG level and liver index (liver weight/body weight) were significantly higher than those in controls. Body weight, serum transaminase levels (ALT, AST) and lipids (TG, TC) were also remarkably elevated in the model group D compared with the control and model C groups (Table 1).
Table 1
Serum and liver biochemistry. All data are expressed as the mean ± standard deviation; Control: the normal control animals fed a standard diet; Model: the model animals fed a HFD; 4W: 4 weeks; 12W: 12 weeks; * P<0.05 vs. control; ** P<0.01 vs. control
Fig. 1
High-fat feeding induces hepatic steatosis and steatohepatitis. (A): Rats fed a standard diet for 4 weeks. Normal liver without vacuolar degeneration and no distinct inflammatory cell infiltration as detected by H&E staining (image magnification: ×200). (B): Rats fed a HFD for 4 weeks. Pronounced microvesicular steatosis and predominantly swollen, ballooned hepatocytes with occasional foci of inflammatory cells in the hepatic lobules, as detected by H&E staining (image magnification: ×200). (C): Rats fed a standard diet for 12 weeks. Mild microvesicular steatosis without inflammatory infiltration or fibrosis, as detected by H&E staining (image magnification: ×200). (D): Rats fed a HFD for 12 weeks. Pronounced steatohepatitis with some hepatocytes showing severe ballooning, degeneration and foci of inflammatory cell infiltration in the hepatic lobules, as detected by H&E staining (image magnification: ×200).
Deep sequencing revealed unique liver miRNA expression profiles
Deep sequencing has been used to profile miRNA profiles in several species, including plants, animals and humans [29,30,31]. Although the expression profile of miRNAs in rats and humans with NAFLD has been previously reported [32,33], we applied Illumina HiSeq deep sequencing to study microRNA expression levels. Briefly, this technique provides a quantitative measure of the presence and prevalence of transcripts from known and previously unknown genes. In total, the two libraries contained 49,478,704 reads, representing 1,665,677 unique sequences (Fig. 2). We focused our analysis on the 33,870,904 reads that matched the mouse genome with 100% identity over their entire length. The analysis showed that 57.1% (11,689,671) of all reads from the control library and 50% (14,476,560) of all reads from the model library matched the annotated miRNA hairpins. Meanwhile, 53 unannotated reads from the control and 64 unannotated reads from the model library were previously uncharacterised RNAs, referred to as “novel” sequences (data not shown).
Fig. 2
Quality control of small RNA deep sequencing. (A): Small RNA annotation of the total reads. The above pie chart displays the composition and proportion of small RNAs, with a total of 20,478,816 reads in the control group. The lower chart contains the 28,999,888 reads for NAFLD model group. (B): Length distribution of deep sequencing results. This figure shows the percent frequency of different length clean reads in each group.
When matching the miRBase 16.0, 239 mature miRNAs were identified in both the control and model groups (Fig. 3). Among these known miRNAs, 44 miRNAs were up-regulated and 12 were down-regulated in the model group compared with the control group (Table 2). Although our sequencing results showed that the expression levels of miR-146b and miR-152 in the model group were not 2-fold higher than those in the control group, a recent study indicated that miR-146a and miR-146b regulate the inflammatory response [34]. In addition, miR-152 was found to be up-regulated in many metabolic disorders, such as type 1 diabetes [35]. Thus, we selected miR-200a, miR-200b, miR-200c, miR-146a, miR-146b and miR-152 for further analysis.
Table 2
The dysregulated miRNAs in NAFLD, as determined by deep sequencing. The miRNA fold change presented here represents the mean value of sequenced samples from the two groups. “Control” represents the normal 12-week liver group, and “Model” represents the model 12-week liver group. Bold text indicates a change >1.5; Italicised text indicates a change <0.5
Fig. 3
Scatter plot of all miRNAs by deep sequencing. Comparison of known miRNA expression between the two groups to identify differentially expressed miRNAs. The X-axis and Y-axis show the expression levels of miRNAs in the control and NAFLD samples, respectively. Each point in the figure represents one miRNA.
Differential miRNA expression in HFD-induced fatty liver
Stem-loop real-time PCR was used to confirm the deep sequencing results in six randomly selected samples from normal group B and model group D. All of the six chosen miRNAs showed the same expression patterns as observed in the sequencing experiments, further supporting the effectiveness of deep sequencing for screening dysregulated miRNAs (Fig. 4A). To detect miRNA expression in different stages of NAFLD, we also detected the expression of the six miRNAs in 4-week liver samples (Fig. 4B). Four of the selected miRNAs (miR-200a, miR-200b, miR-146b and miR-152) were overexpressed in the model group (P<0.05), whereas two were not (miR-200c, miR-146b). To address the role of miRNAs in NAFLD, we further examined miRNA expression in HepG2 cells and human hepatocytes cultured in complete medium with or without 1 mM of FFA mixture (oleate/palmitate) for 24 h. HepG2 cells (Fig. 5A and 5D) and human hepatocytes (Fig. 5B and 5E) exposed to the FFA mixture developed lipid droplet accumulation. Consistent with this finding, the hepatic TG content was significantly higher in cells treated with the FFA mixture (Fig. 5C). The expression of certain miRNAs was also increased in fat-overloaded HepG2 cells and human hepatocytes (Fig. 5F and 5G).
Fig. 4
Validation of the differentially expressed miRNAs in NAFLD. miR-200a, miR-200b, miR-200c, miR-146a, miR-146b and miR-152 were chosen to evaluate the NAFLD models using real-time PCR. (A): The relative expression levels of the six chosen miRNAs in rats fed a HFD for 12 weeks were strikingly greater than those in the control group fed the standard diet for 12 weeks (n=10). (B): Except for miR-200c and miR-146a, the expression levels of the miRNAs in rats fed a HFD for 4 weeks were significantly increased over the control group (n=10). The data are expressed as the mean ± standard deviation, * P<0.05.
Fig. 5
Lipid accumulation test using HepG2 cells and human hepatocytes. (A, B) Representative Oil Red O staining. (C) The TG content in HepG2 cells and human hepatocytes. (D, E) Integrated option density (IOD) determined by Image Pro-Plus according to the staining intensity by visual assessment. (F, G) The expression levels of the six miRNAs in HepG2 cells cultured in medium with the FFA mixture (1 mM) (C) and human hepatocytes (D).* P<0.05.
The effect of proinflammatory factors on miRNA expression
As both TNF-α and IL-6 play important roles in NAFLD [20,36,37,38], we assessed whether these proinflammatory factors could alter miRNA expression. HepG2 cells were cultured for 24 h in complete medium containing FFAs and then cultured for an additional 24 h with the proinflammatory factors TNF-α (10 ng/ml) [21] and IL-6 (25 ng/ml) [22]. The cells were collected after 0, 4, 8, and 24 h to determine the expression levels of the six selected miRNAs. Increased expression levels of the six miRNA were observed in HepG2 cells after treatment with TNF-α and IL-6 (Fig. 6).
Fig. 6
The miRNA expression in HepG2 cells cultured with the FFA mixture and TNF-α or IL-6 for 24 h. The data are expressed as the mean ± standard deviation, * P<0.05 vs. 0 h. (Experiments performed 3 times) (A): 10 ng/ml TNF-α treated HepG2 cells for 0, 4, 8, and 24 h. (B): 25 ng/ml IL-6 treatment led to increased expression of the six miRNAs at 8 h and 24 h, compared with 0 h.
GO and KEGG pathway analysis
Predicted gene targets of all differentially expressed miRNAs were obtained using two databases: TargetScan and miRanda. A total of 2,821 target genes of the up-regulated miRNAs and 507 target genes of the down-regulated miRNAs were identified. All of the predicted targets were categorised into different functional classes, which were identified using gene ontology. In our survey of existing data, the target gene function of up-regulated miRNAs mainly involved the following processes: (1) GO:0045944 - positive regulation of transcription from RNA polymerase II promoter; (2) GO:0006470 - protein amino acid dephosphorylation; (3) GO:0007156 - homophilic cell adhesion; (4) GO:0006813 - potassium ion transport; (5) GO:0030335 - positive regulation of cell migration; and (6) GO:0042493 - response to drug. Meanwhile, the target gene function of down-regulated miRNAs mainly involved the following processes: (1) GO:0016192 - vesicle-mediated transport; (2) GO:0006886 -intracellular protein transport; (3) GO:0009401 - phosphoenolpyruvate-dependent sugar phosphotransferase system; (4) GO:0060021 - palate development; (5) GO:0006813 -potassium ion transport; and (6) GO:0006606 - protein import into the nucleus (Fig. 7). GO enrichment analysis showed that the mRNA clusters were significantly enriched in basal biological processes, such as material transport, differentiation and development of tissues and organs, and transcriptional regulation.
Fig. 7
Significant miRNA targets by GO analysis. The GOs targeted by up-regulated miRNAs. The GOs targeted by down-regulated miRNAs. All of these GOs show increased enrichment. The vertical axis is the GO category, and the horizontal axis is the enrichment of GO.
We also utilised KEGG pathway analysis to address higher-level functions associated with organism behaviours and cellular processes from our sequencing data. In our survey of existing data, the target gene function of up-regulated miRNAs mainly involved the following functions: (1) MAPK signalling, (2) focal adhesion signalling, (3) insulin signalling, (4) leukocyte transendothelial migration, (5) TGF-beta signalling, (6) cytokine-cytokine receptor interactions, (7) adherens junction adipocytokine signalling, and (8) Wnt signalling. Meanwhile, the target gene function of down-regulated miRNAs mainly involved the following functions: (1) MAPK signalling, (2) long-term depression, (3) insulin signalling, (4) SNARE interactions in vesicular transport, (5) gonadotropin-releasing hormone (GnRH) signalling, (6) long-term potentiation, (7) TGF-beta signalling, and (8) regulation of the actin cytoskeleton (Fig. 8). The significant enrichment of pathways associated with insulin signalling and TGF-β signalling indicated that genes with functions potentially relevant to the pathogenesis of NAFLD are targeted by the detected miRNAs [39]. Furthermore, target genes of the up-regulated miRNAs were found to be involved in cytokine-cytokine receptor interactions and the adipocytokine signalling pathway (Fig. 8). Taken together, this information demonstrates that the up-regulated and down-regulated miRNAs co-regulate genes involved in the pathogenesis of NAFLD.
Fig. 8
Pathway analysis based on miRNA target genes. Significant pathways targeted by up-regulated (A) and down-regulated (B) miRNAs. The vertical axis is the pathway category, and the horizontal axis is the pathway enrichment.
Discussion
NAFLD is a multifactorial disease that is predominantly regulated by genetic predisposition coupled with environmental factors [40]. Insulin resistance and abnormal lipid metabolism have been accepted as the main driving force behind this disease process. Additionally, inflammation, apoptosis and oxidative stress all play critical roles in NAFLD. The growing application of high-throughput technology has brought miRNAs into focus as potentially important mediators of NAFLD.
There are several animal models of NAFLD, each with advantages and disadvantages [41]. Here, we adopted a commonly used nutritional model, which utilises a high-fat and high-cholesterol diet to induce fatty liver disease. Using this model, our animals developed visceral obesity and hyperlipidaemia after 4 weeks with high-fat feeding, in addition to hepatic steatosis as detected by liver histology. At the end of 12 weeks, the rats showed an abnormal hepatic enzyme spectrum, and histological analysis revealed the presence of hepatitis. All of the above data indicate that we were able to use a HFD to create a NASH model that reflects the major features of human NAFLD.
To identify the miRNA expression profile of NAFLD, we employed miRNA deep sequencing and validated selected, aberrantly expressed miRNAs by stem-loop RT-PCR. Compared to microarray platforms, the data obtained by deep sequencing is more precise and sensitive [42]. Furthermore, deep sequencing allows the prediction of novel miRNAs using bioinformatic algorithms; in this regard, we predicted more than 100 novel miRNAs that require further experimental verification (data not shown). Overall, there were 44 up-regulated miRNAs and 12 down-regulated miRNAs. Within the up-regulated group, the miR-200 family was the most notable. Interestingly, each member in this family, including miR-200a, miR-200a*, miR-200b, miR-200b*, miR-200c, miR-141 and miR-429, was significantly up-regulated in the NAFLD group. This finding is partially consistent with the characteristics of NAFLD previously observed in diet-induced Sprague-Dawley rats [43], ob/ob mice [43], methyl-deficient mice [44], and NASH patients [33]. In the report by Jin et al., miR-200b expression was increased in the early stages of NAFLD but was decreased in the inflammatory stage of disease [32]. In addition, miR-146 miRNAs (miR146a and miR-146b), first identified as immune system regulators, were also up-regulated [45]. In agreement with this result, miR-146b was reported to be overexpressed in humans with NASH [33]. However, three up-regulated miRNAs in our study (miR-880, miR-881 and miR-741-3p) were not previously reported. For the down-regulated group, miR-182 and miR-183 have been reported to exhibit decreased expression in rodents with alcoholic steatohepatitis and increased expression in mice fed a methionine choline-deficient diet [46]; these differences may be attributed to the different models of NASH. Another down-regulated miRNA, miR-451, was also shown to be decreased in diet-induced NAFLD [43], as well as in ob/ob mice [44]. Consequently, miR-200a/b/c, miR-146a/b, and miR-152 were chosen for further validation, and our validation experiment using the animal model of NASH revealed that the expression levels of the six miRNAs were increased during disease development. Our data also demonstrated that these selected miRNAs regulate the gene expression profile involved in NAFLD.
TNF-α and IL-6 are produced by various cell types, including monocytes/macrophages, T cells, fibroblasts, and endothelial cells. TNF-α and IL-6 are also secreted by adipocytes and are therefore referred to as adipokines [47]. In particular, TNF-α is a key factor in the development of NAFLD because it can induce inflammation and insulin resistance by activating proinflammatory pathways, such as the NF-κB and c-jun N-terminal kinase pathways [36]. An elevated TNF-α level has also been demonstrated in several studies of fatty liver disease [20]. The role of IL-6 in liver pathology is complex, although it has been associated with protection from steatosis via gp130 signalling and activation of the Jak/Stat pathway in animal models, similar to leptin [37]. Nevertheless, IL-6 is considered a predictive marker of insulin resistance via the induction of SOCS3 expression [48], and serum IL-6 levels are elevated in patients with NAFLD [38]. As such, both TNF-α and IL-6 may be associated with hepatic lipogenesis, obesity and insulin resistance, supporting a central role for these inflammatory factors in NAFLD. To determine the effects of these representative proinflammatory factors on miRNA expression, we used them to treat HepG2 cells. Similar to the results obtained from the sequencing results, both TNF-α and IL-6 up-regulated the expression of the six chosen miRNAs. We interpret this result as indicating that proinflammatory factors are key mediators of the changes in miRNA expression in NAFLD.
As TNF-α and IL-6 were found to mediate the aberrant expression of miRNAs associated with NAFLD, we also analysed the gene mapping and transcription factor binding sites around the differentially expressed miRNAs (Fig. 9). The gene mapping analyses showed that miR-200a, miR-200b and miR-200c have an adjacent genomic location and may be transcribed as a cluster. Alternatively, miR-146a and miR-146b have their own transcriptional regulatory region. We also found that miR-152 is parasitic in the COPZ2 gene, and both function as tumour suppressor genes. Furthermore, miR-30a, miR-30b, and miR-30d share a transcribed region, and miR-30c and miR-30e are both located in the NFYC gene, which encodes MHC I and MHC II molecules. Therefore, these miRNAs are likely to play an important role in antigen presentation. The down-regulated miRNAs, including miR-182, miR-183 and miR-96, are all located in the NRF1 gene, a key gene that mediates oxidative stress. Additionally, miR-1 and miR-133 are both parasitic in the MIB1 gene, which regulates apoptosis [49]. Furthermore, NF-κB can regulate the expression of miR-200b, miR-30c and miR-152 [50,51,52], whereas miR-182, miR-183, miR-96 and miR-1 can be regulated by FOXO3 [53,54].
Fig. 9
Potential mechanism for miRNAs associated with NAFLD. The red frames represent up-regulated genes, and the green frames represent down-regulated genes. The yellow boxes represent inhibited biological processes, and the purple boxes represent stimulatory biological process. Two pathways associated with NAFLD were proposed; please refer to the Discussion for details. TNF-α: tumour necrosis factor-α; NF-κB: nuclear factor kappa-light-chain-enhancer of activated B cells; PPAR: peroxisome proliferator-activated receptors; ATG: autophagy; PI3K: phosphatidyl inositol 3-kinase; AKT: also known as protein kinase B (PKB); FOXO3: Forkhead box O3.
When combining the effect of these two inflammatory factors, two pathways associated with NAFLD were proposed (Fig. 9). In the first pathway, TNF-α stimulates NF-κB activation, which may up-regulate mir-200/mir-30/mir-152 expression, causing abnormal glucose and lipid metabolism/autophagy through targeting PPAR/ATG genes. In the second pathway, IL-6 may activate the PI3K/AKT pathway, thereby down-regulating the expression of the downstream transcription factor FOXO3. This down-regulation would ultimately lead to a lower expression of mir-1/mir-183, causing Smad4 gene expression de-repression and subsequent fibrosis. However, these hypotheses require further experimental evidence.
In summary, this study has shown that NAFLD is associated with changes in the expression of hepatic miRNAs and that the target genes of aberrantly expressed miRNAs are related to the pathogenesis of NAFLD. We also demonstrated that inflammatory factors are crucial mediators of the aberrant expression of miRNAs and propose two pathways leading to the described changes. Although this study was based on a small sample size and included only rodents, we believe that our findings offer novel insights into the pathogenesis of NAFLD in humans.
Acknowledgments
This work was supported by grants from the National Natural Science Foundation of China (81270938), Zhejiang Provincial Key Medical Disciplines (Innovation Discipline, 11-CX24), Zhejiang Province key scientific and technological innovation team (2010R50050), and the National Key Technology R&D Program of China (2012BAI02B03).
Disclosure Statement
The authors declare no conflicts of interest in relation to the content of this article.


Get Permission
