Differential expression of microRNAs in mouse liver under aberrant energy metabolic status.

Despite years of effort, exact pathogenesis of nonalcoholic fatty liver disease (NAFLD) remains obscure. To gain an insight into the regulatory roles of microRNAs (miRNAs) in aberrant energy metabolic status and pathogenesis of NAFLD, we analyzed the expression of miRNAs in livers of ob/ob mice, streptozotocin (STZ)-induced type 1 diabetic mice, and normal C57BL/6 mice by miRNA microarray. Compared with normal C57BL/6 mice, ob/ob mice showed upregulation of eight miRNAs and downregulation of four miRNAs in fatty livers. Upregulation of miR-34a and downregulation of miR-122 was found in livers of STZ-induced diabetic mice. These results demonstrate that distinct miRNAs are strongly dysregulated in NAFLD and hyperglycemia. Comparison between miRNA expressions in livers of ob/ob mice and STZ-administered mice further revealed upregulation of four miRNAs and downregulation of two miRNAs in livers of ob/ob mice, indicating that these miRNAs may represent a molecular signature of NAFLD. A distinctive miRNA expression pattern was identified in ob/ob mouse liver, and hierarchical clustering of this pattern could clearly discriminate ob/ob mice from either normal C57BL/6 mice or STZ-administered mice. These findings suggest an important role of miRNAs in hepatic energy metabolism and implicate the participation of miRNAs in the pathophysiological processes of NAFLD.

Nonalcoholic fatty liver disease (NAFLD) is the most common form of chronic liver disease worldwide and is becoming a major public health concern in modern society (1)(2)(3). NAFLD is a clinicopathological syndrome characterized by excess fat accumulation in liver, which ranges from simple steatosis to steatohepatitis and cirrhosis in the absence of heavy alcohol consumption ( 4 ). Disorders related to metabolic syndrome, such as obesity, type 2 diabetes mellitus, and dyslipidemia, were identifi ed as the main risk factors for the development of NAFLD (5)(6)(7). Although a widely accepted two-hit hypothesis ( 8 ) may partially explain the progressive liver damage by nonalcoholic steatosis and steatohepatitis, much of the pathogenesis of NAFLD remains undiscovered and requires further study.
Recently, a new class of RNA regulatory genes known as microRNAs (miRNAs) has been found to introduce a whole new layer of gene regulation in eukaryotes. miRNAs are endogenous noncoding RNAs of 19-24 nucleotides in length that play an important role in the negative regulation of gene expression by base-pairing to complementary sites on the target mRNAs, causing a block of translation or degradation of the target mRNA ( 9 ). In addition to the fundamental roles in diverse biological and pathological processes, including developmental timing, apoptosis, proliferation, differentiation, organ development, carcinogenesis, and immune response (10)(11)(12), miRNAs are also reported to play important roles in energy metabolism, both in invertebrates and vertebrate animals. The role of miRNAs in energy metabolism was fi rst indicated by a study in the fruit fl y Drosophlia melanogaster , suggesting an important role of miR-14 in energy metabolism on the Abstract Despite years of effort, exact pathogenesis of nonalcoholic fatty liver disease (NAFLD) remains obscure. To gain an insight into the regulatory roles of microRNAs (miRNAs) in aberrant energy metabolic status and pathogenesis of NAFLD, we analyzed the expression of miRNAs in livers of ob/ob mice, streptozotocin (STZ)-induced type 1 diabetic mice, and normal C57BL/6 mice by miRNA microarray. Compared with normal C57BL/6 mice, ob/ob mice showed upregulation of eight miRNAs and downregulation of four miRNAs in fatty livers. Upregulation of miR-34a and downregulation of miR-122 was found in livers of STZ-induced diabetic mice. These results demonstrate that distinct miRNAs are strongly dysregulated in NAFLD and hyperglycemia. Comparison between miRNA expressions in livers of ob/ob mice and STZ-administered mice further revealed upregulation of four miRNAs and downregulation of two miRNAs in livers of ob/ob mice, indicating that these miRNAs may represent a molecular signature of NAFLD. A distinctive miRNA expression pattern was identifi ed in ob/ ob mouse liver, and hierarchical clustering of this pattern could clearly discriminate ob/ob mice from either normal C57BL/6 mice or STZ-administered mice. These fi ndings suggest an important role of miRNAs in hepatic energy metabolism and implicate the participation of miRNAs in the pathophysiological processes of NAFLD. -Li, S., X. Chen ( 23,24 ). Briefl y, 50 g of total RNA was purifi ed using the mirVANA miRNA isolation kit (Ambion, Austin, TX) to enrich small RNA fraction. Purifi ed RNA was labeled with fl uorescein, and hybridization was carried out on the CapitalBio Mammalian miRNA Array V 3.0 (CapitalBio, Beijing, China) containing 509 probes in triplicate, corresponding to 435 human and mouse mature miRNA genes. Liver RNA sample from each mouse was analyzed on individual chip. Finally, hybridization signals were detected, and scanner images were quantifi ed. These quantifi ed signal intensity values of microarray were normalized to per-chip mean values. The miRNA microarray data have been deposited into Gene Expression Omnibus (GEO accession number GSE13840) and can be prepared in a fully MIAME-compliant manner.

Real-time quantitative RT-PCR of mature miRNAs
Stem-loop quantitative RT-PCR (qRT-PCR) assays to quantify the mature miRNAs were performed as previously described ( 25,26 ) using fl uorescent nucleic acid dye. The RT primers and realtime PCR primers were designed as described ( 25 ). Briefl y, 1 g of total RNA was reverse transcribed under the following conditions: 16°C for 15 min, 42°C for 60 min, and 85°C for 5 min. The 20 l PCR included 1 l RT product and 1 l EvaGreen dye (Biotium, Hayward, CA). The conditions for the PCR reaction were as follows: 95°C for 5 min followed by 40 cycles of 95°C for 15 s and 60°C for 1 min using an ABI PRISM 7300 thermal cycler. All reactions were run in triplicate. The threshold cycle (C T ) is defi ned as the fractional cycle number at which the fl uorescence passes the fi xed threshold. The miRNA expression levels were normalized to U6 RNA. The relative expression was calculated using the comparative ⌬ ⌬ C T method, and the values were expressed as 2 Ϫ ⌬ ⌬ CT ( 27 ).

Statistical analysis
Differentially expressed miRNAs were identifi ed using the t -test procedure within signifi cance analysis of microarrays (SAM) ( 28 ) with a 5% false discovery rate (FDR) threshold. Cluster and Java TreeView were used to build the unsupervised tree. The genes and arrays were mean centered, and hierarchical trees were built using correlation metrics. Real-time qRT-PCR assays were performed in triplicate. Data were expressed as means ± SE, and the differences with P < 0.05 were considered statistically signifi cant using a two-sided unpaired Student's t -test.

Distinctive patterns of miRNA expression in livers of ob/ob mice
To investigate the potential involvement of miRNAs in NAFLD and hyperglycemia, we performed a microarraybased analysis of miRNAs in livers of ten ob/ob mice and eight C57BL/6 mice (served as controls), respectively. As shown in Table 1 , the body weight, blood glucose levels, serum insulin levels, and liver triglyceride concentration of the ob/ob mice were all signifi cantly elevated compared with those of the C57BL/6 mice, suggesting the whole-animal level ( 13 ). Another study showed an involvement of miR-278 in energy homeostasis of Drosophila ( 14 ). For vertebrates, miRNAs were found to regulate adipocyte differentiation ( 15,16 ), and glucose or potassium stimulated insulin secretion from pancreatic islet cells ( 17,18 ). Esau et al. ( 19 ) applied an miR-122 antagonist intraperitoneally into mice and found that this procedure could reduce plasma cholesterol level, increase hepatic fatty acid oxidation, and decrease hepatic fatty acid and cholesterol synthesis rate. Jin et al. ( 20 ) also reported that miR-122 level was over 2-fold upregulated in NAFLD rat. Together, these fi ndings suggested a strong connection between miRNAs and energy metabolism. However, whether miRNAs play a role in the pathophysiological processes of NAFLD remains to be elucidated.
To address this question, we employed miRNA microarray to identify differentially expressed miRNAs in livers of ob/ob mice [a well-established model of NAFLD with hyperglycemia ( 21,22 )] and streptozotocin (STZ)-induced type 1 diabetic mice (a model simply under high blood glucose without fatty liver) to characterize the potential roles of miRNAs in livers under NAFLD and hyperglycemia, respectively.

Animal models
Adult (aged 11-12 weeks) obese Lep ob/ob C57BL/6 mice (six male and four female) and normal control C57BL/6 mice (aged 8-12 weeks; four male and four female) were purchased from Model Animal Research Center of Nanjing University. For STZ-induced type 1 diabetic model, 8-week-old wild-type C57BL/6 mice (eight male) received a single intraperitoneal injection of 150 mg/kg STZ dissolved in citrate buffer at pH 4.5. Mice were fasted overnight before glucose measurements. Blood samples were taken from the tail vein 3, 7, and 14 days after STZ injection, respectively. Animals were euthanized 14 days after STZ injection. The mice with fasting blood glucose у 12 mmol/L were considered diabetic. All mice were housed in a temperature-controlled environment with a 12 h light/dark cycle and free access to water and a standard chow diet containing 60% carbohydrate, 13% fat, and 27% protein on a caloric basis. All animals were weighed and taken fasting blood glucose measurements and serum samples before euthanized. The liver tissues were dissected and immersed in liquid nitrogen immediately and stored at Ϫ 80°C until used for subsequent analysis. Serum insulin levels of these animal models were quantifi ed by the Rat/Mouse Insulin ELISA Kit (Linco Research, St. Charles, MO). Liver triglyceride concentration was determined using the Triglyceride Kit (Jiancheng Bioengineering Institute, Nanjing, China). All animal care and handling procedures were carried out in accordance with the Institute for Laboratory Animal Research Guide for Care and Use of Laboratory Animals and approved by the Institutional Review Board of Nanjing University.

RNA isolation and microarray experiments
Frozen liver tissues were dissected and homogenized, and total RNA were extracted using Trizol Reagent (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. RNA labeling and hybridization on miRNA microarray chips were conducted and stayed high in the following week, while the body weight of these mice decreased gradually (data not shown). On the 14th day after STZ injection, the mean fasting blood glucose level of STZ-administered mice was almost 4-fold higher than that of control C57BL/6 mice. By contrast, the average body weight of STZ-administered mice was signifi cantly decreased compared with that of C57BL/6 mice ( Table 1 ). At this stage, mice were killed and liver tissues were dissected for RNA extraction and miRNA microarray detection subsequently. SAM analysis of miRNA microarray data (with a 2-fold change threshold and 5% FDR) showed that only miR-34a was upregulated in liver samples of STZ-induced diabetic mice compared with samples of C57BL/6 mice ( Table 2 ). Unsupervised hierarchical clustering of miRNA expression from liver samples of STZ-induced diabetic mice and C57BL/6 mice is depicted in Fig. 1B . Sample dendrogram in columns showed no clear separation between STZ samples and C57BL/6 samples, indicating that general miRNA expression in mice livers under 2-week hyperglycemia was not signifi cantly altered, except for miR-34a.

Differentially expressed miRNAs in livers of ob/ob mice versus STZ-induced diabetic mice
To fi nd out the miRNAs that may be only involved in the pathophysiological processes of NAFLD, we next compared the miRNA expression pattern in livers of ob/ob mice with that of STZ-induced diabetic mice. Here ob/ob mice served as an animal model under both NAFLD and hyperglycemia conditions, while STZ-induced diabetic mice served as a model under a simple status of hyperglycemia. As shown in Fig. 2 , unsupervised hierarchical clustering of liver miRNA expression of ob/ob mice and STZ-induced diabetic mice generated a sample dendrogram with two major branches in columns, clearly discriminating these two groups of samples. This result suggests a unique miRNA expression pattern in nonalcoholic fatty livers, excluding the infl uence of hyperglycemia on the hepatic miRNA expression. Under 5% FDR and 2-fold change threshold, SAM algorithm identifi ed four upregulated miRNAs (miR-103, miR-31, miR-107, and miR-126-3p) and two downregulated miRNAs (miR-100 and miR-29c) in ob/ob samples compared with STZ samples ( Table 3 ), indicating that these miRNAs may represent a molecular signature of NAFLD.

Dysregulated liver miRNAs under aberrant metabolic status segregated by gender
In this study, both male and female mice (four male and four female normal C57BL/6 mice; six male and four female ob/ob mice; and eight male STZ-induced diabetic mice) were used. To further analyze differentially expressed liver miRNAs in these animal models of different gender, we compared the liver miRNA expression in different animal models of male and female, respectively. As shown in supplementary Tables I and II, SAM analysis of miRNA microarray data (with a 2-fold change threshold and 5% FDR) listed 18 dysregulated miRNAs in livers of male ob/ob mice versus male normal aberrant energy metabolic status of ob/ob mice. Next, total RNAs were extracted from liver tissues of ob/ob mice and C57BL/6 mice and subsequently used in miRNA microarray analysis. The discrimination of miRNA expression in livers of ob/ob mice and C57BL/6 mice was clearly revealed by a two-way (genes against samples) unsupervised hierarchical clustering analysis of the miRNA microarray data. As displayed in Fig. 1A , the sample dendrogram generated by hierarchical cluster analysis showed two major branches in columns (ob/ob mice vs. C57BL/6 mice), indicating a specifi c miRNA expression pattern in the fatty livers of ob/ob mice versus the normal livers of C57BL/6 mice.
To further identify differentially expressed miRNAs, a statistical method termed SAM ( 28 ) was employed. SAM is an algorithm that calculates a score for each gene and therefore identifi es genes that are signifi cantly associated with an outcome variable, such as the type of analyzed tissue (NAFLD vs. normal). Employing two-class unpaired analysis within SAM, the miRNA expression levels in livers of ob/ob mice were compared with those of C57BL/6 mice. In this test, an FDR ( q -value) <5% was selected. The expressions of miRNAs were considered signifi cantly altered only when they fulfi lled two criteria: 1 ) mean fold change > 2 or < 0.5; 2 ) q -value < 5%. Based on these principles, SAM analysis generated a list of 11 miRNAs that were differentially expressed in livers of ob/ob mice ( Table 2 ). Among these miRNAs, eight miRNAs (miR-34a, miR-31, miR-103, miR-107, miR-194, miR-335-5p, miR-221, and miR-200a) were at least 2-fold upregulated, and three miRNAs (miR-29c, miR-451, and miR-21) were at least 0.5fold downregulated in ob/ob mice liver samples. Collectively, these data demonstrate that distinct miRNAs are signifi cantly regulated in NAFLD with hyperglycemia, implicating the potential roles of miRNAs as functional modulators in the pathophysiological processes of these diseases.

miRNA expression in livers of STZ-induced type 1 diabetic mice
By causing rapid and irreversible necrosis of pancreatic ␤ -cells ( 32 ), STZ-induced type 1 diabetic mice provide us an animal model of high blood glucose but without fatty liver. This diabetic mouse model was employed here to study the relationship between hepatic miRNA expression and hyperglycemia. Within 7 days of STZ injection, fasting blood glucose levels of the mice were all above 12 mmol/L a Fasting blood glucose levels and fasting serum insulin levels measured before the mice were euthanized; ** P < 0.01. induced diabetic mice. By unsupervised hierarchical clustering of miRNA expression in different liver samples segregated by gender, clear separation of ob/ob mouse samples from normal C57BL/6 mouse samples and of ob/ob mouse samples from STZ-induced diabetic mouse C57BL/6 mice, 23 dysregulated miRNAs in livers of female ob/ob mice versus female C57BL/6 mice, four dysregulated miRNAs in livers of male STZ-induced diabetic mice versus male C57BL/6 mice, and 11 dysregulated miRNAs in livers of male ob/ob mice versus male STZ-

Validation of differentially expressed miRNAs by real-time qRT-PCR
To validate the accuracy of the miRNA microarray data, we carried out stem-loop qRT-PCR assay ( 25,26 ) with the same RNA preparations used in microarray analysis. As shown in Fig. 3 , differentially expressed miRNAs identifi ed by SAM algorithm based on the microarray results were analyzed using qRT-PCR assay. Results showed that the directions of the changes in miRNA expression were concordant between the two platforms.
Since miR-122, a liver-specifi c miRNA, makes up 70% of all miRNAs in liver ( 33,34 ) and is reported to regulate hepatic lipid metabolism ( 19,35,36 ), we analyzed miR-122 variation using SAM algorithm with 2-fold change as a threshold and also quantifi ed it using real-time qPCR. As shown in Fig. 3 , miR-122 was signifi cantly reduced in livers of both ob/ob mice and STZ-induced diabetic mice compared with C57BL/6 mice, suggesting the involvement of miR-122 in hepatic energy metabolism.

Predicted gene targets of differentially expressed miRNAs
Lack of reliable and specifi c methods for biological target validation hampers the full understanding of the mechanisms by which miRNAs execute their functions. Only a few miRNAs have so far been assigned target mRNAs, and the Ϫ 0.75 0.85 0.00 a q -value: false discovery rate, the expected percentage of genes identifi ed by chance.
b These miRNAs were not considered as signifi cantly changed by SAM due to insuffi cient fold change but proved to be signifi cantly changed by real-time qRT-PCR. conventional methodologies are still labor intensive. Three computer-aided algorithms, including TargetScan ( 29 ), miRanda ( 31 ), and PicTar ( 30 ), however, have partially overcome this limitation. With these bioinformatic approaches, we obtained a list of predicted targets that may be potentially involved in metabolic diseases ( Table 4 ). As illustrated by underlines in Table 4 , several such putative targets related to energy metabolism or NAFLD were derived by using various algorithms. To further elucidate the correlation between the miRNAs and the predicted targets, these prediction data were compared with previously published liver mRNA expression profi les of obese Lep ob/ob C57BL/6 mice ( 37,38 ). Some of the putative targets were shown dysregulated in livers of ob/ob mice ( Table 4 ). Further identifi cation of these target genes and their characterization will be greatly helpful for better understanding the potential roles of these differentially expressed miRNAs in hepatic energy metabolism.

DISCUSSION
Although progress has been made in the understanding of causative mechanisms of nonalcoholic steatosis, steatohepatitis, and cirrhosis, the knowledge about underlying regulatory networks of these pathophysiological conditions is still insuffi cient. Compelling evidence has demonstrated the substantial regulatory role of miRNAs in energy metabolism and liver functions (17)(18)(19)39 ). Because of the spatial and temporal expression pattern of miRNAs and their conservation across species, miRNAs are believed to play similar roles in all animal cells (40)(41)(42). In light of these fi ndings, we have identifi ed mouse miRNAs that are possibly involved in the pathophysiological processes of NAFLD. Most of the varied miRNAs are highly conserved in mammals and are expected to play similar roles in the pathology of mouse and human NAFLD. Our results strongly suggest a regulatory role of miRNAs in hepatic energy metabolism.

A distinctive miRNA expression pattern in fatty livers of ob/ob mice
Due to inherited defi ciency of the appetite-suppressing hormone, leptin ( 43 ), ob/ob mice spontaneously develop obesity, insulin resistance ( 44 ), and fatty livers ( 45 ). Therefore, ob/ob mice represent a naturally occurring model of NAFLD with excessive fat accumulation in liver ( 21,22 ).
In this study, we found that 11 miRNAs in liver samples of ob/ob mice were signifi cantly dysregulated (>2-fold change) compared with those of wild-type C57BL/6 mice. Among those altered miRNAs, eight miRNAs (miR-34a, miR-31, miR-103, miR-107, miR-194, miR-335-5p, miR-221, and miR-200a) were upregulated and three miRNAs (miR-29c, miR-451, and miR-21) were downregulated. A distinctive miRNA expression pattern was identifi ed in ob/ob mouse liver, and hierarchical clustering of this pattern could clearly discriminate ob/ob mice from normal C57BL/6 mice and STZ-administered mice, indicating the participation of miRNAs in the pathophysiological processes of NA-FLD with hyperglycemia. db/db mice, another animal model of type 2 diabetes mellitus with leptin receptor mutations, also develop NAFLD. Further analysis of miRNA expression in livers of db/db mice are required in the future to elucidate the primary effect of leptin defi ciency on liver miRNA expression.

miRNA expression in livers of STZ-induced type 1 diabetic mice
Since ob/ob mice can spontaneously developed obesity and insulin resistance, they are also considered as animal models of type 2 diabetes mellitus with high blood glucose levels. Glucose has been demonstrated to regulate transcription of genes encoding not only lipogenic and glycolytic enzymes but also proteins involved in global cell functions ( 46 ). However, it still remains unknown whether high glucose exerts an impact on the miRNA expression. STZ, an antibiotic extracted from Streptomyces acromogenes ( 47 ), is widely used to make experimental animal models of type 1 diabetes mellitus via inducing rapid and irreversible necrosis of ␤ -cells ( 32 ). STZ-induced type 1 diabetic mice develop hyperglycemia with decreased body weight but no fatty liver. In this study, we employed the STZadministered mice as a model of hyperglycemia without fatty liver. To illustrate the potential effect of high blood glucose on the liver miRNA expression, we compared the miRNA expressions in livers of STZ-induced type 1 diabetic mice with those of control C57BL/6 mice. Surprisingly, only one miRNA, miR-34a, showed an altered expression above 2-fold change threshold. Clustering analysis also failed to distinguish the STZ-induced diabetic mice from the normal C57BL/6 mice, indicating that short-term (2 weeks) high blood glucose does not affect liver miRNA expression globally but only affect certain miRNA expression in liver. Comparing the varied miRNAs under these two pathologic conditions, miR-34a is upregulated and miR-122 is downregulated in both ob/ob mice and STZ-induced diabetic mice, suggesting that miR-34a and miR-122 may be mainly related to the regulation of glucose metabolism. To discover the differentially expressed miRNAs in NAFLD without the infl uence of hyperglycemia, we further compared liver miRNA expression of ob/ob mice with that of STZ-induced diabetic mice. The result reveals six miRNAs with over 2-fold changed expression, implicating that these miRNAs are more likely to be involved in lipid metabolism and contribute to the pathology of NAFLD. duced plasma cholesterol levels, increased hepatic fatty acid oxidation, and decreased hepatic fatty acid and cholesterol synthesis rate ( 19,35,36,39 ). Considering the abundance of miR-122 in liver and the potential role in energy metabolism, we tested its expression levels in realtime qRT-PCR assay, although the alteration of miR-122 was not >2-fold change threshold in SAM analysis. We found that miR-122 was signifi cantly downregulated in livers of both ob/ob mice and STZ-induced diabetic mice. This outcome seems to be controversial with the previously reported regulatory function of miR-122 in energy metabolism, and we postulate that the reduction of miR-122 in NAFLD and diabetic livers may be part of a self-protection mechanism of liver cells against lipid or glucose overload.

Liver miRNAs associated with aberrant metabolic status
A number of researchers have reported the abnormal miR-34a expressions in different types of cancer ( 48,49 ), suggesting miR-34a as a miRNA component of the p53 network ( 48,50 ). miR-21, another extensively studied miRNA, has been observed to be upregulated in tumors of various origins (51)(52)(53). In our study, miR-34a was the most elevated miRNA in livers of both ob/ob mice and STZ-induced diabetic mice, while miR-21 was signifi cantly downregulated in livers of ob/ob mice, raising intriguing questions about the functional role of these two miRNAs in energy metabolism of liver.
Previous studies also demonstrated that suppressed liverspecifi c miR-122 in livers of mice would resulted in re- Fig. 3. Validation of microarray data by real-time qRT-PCR. A, B: Real-time qRT-PCR analysis of differentially expressed miRNAs in liver tissue RNA samples previously profi led by microarray. Triplicate assays were carried out for each RNA sample, and the relative amount of each miRNA was normalized to U6 snRNA. Data are shown as fold changes of miRNA levels in livers of ob/ob mice (A) and STZ-induced diabetic mice (B) compared with C57BL/6 mice, which is set as 1 (mean ± SE, n = 8) (* P < 0.05; ** P < 0.01). The expression fold changes of the validated miRNAs by SAM algorithm in microarray analysis from ob/ob mice (C) and STZ-induced diabetic mice (D) compared with that from C57BL/6 mice, which is set as 1. FDR ( q -value) <5%.

Target and function prediction of differentially expressed miRNAs related to NAFLD and energy metabolism
Currently, the major diffi culty for understanding the biological functions of miRNAs is determining their specifi c target genes at the posttranscriptional level. Bioinformatics has been a powerful tool in target and function prediction of miRNAs. miR-103 and miR-107, two closely related miRNAs that differ at a single nucleotide ( 54 ), were previously predicted by bioinformatics to affect multiple mRNA targets in pathways that involve cellular acetyl-CoA and lipid levels ( 55 ). In accordance with this prediction, we found that miR-103 and miR-107 were signifi cantly upregulated in ob/ob mice compared with normal C57BL/6 mice and STZ-administered mice. These results strongly suggest the potential role of miR-103 and miR-107 in the pathogenesis of NAFLD and confi rm the importance of bioinformatics in exploring miRNA functions. Employing three major algorithms for miRNA target prediction, TargetScan, PicTar, and miRanda, we analyzed the predicted target genes of the differentially expressed miRNAs identifi ed in livers of ob/ob mice and STZ-induced diabetic mice. The predicted target genes of the dysregulated miRNAs include some genes that are closely related to glucose or lipid metabolism in liver and may play important roles in NAFLD. The interaction of miRNAs and their target genes may lead to translational repression and, in many cases, mRNA degradation ( 56 ). Thus, we compared the predicted targets with previously published liver mRNA expression profi les of obese Lep ob/ob C57BL/6 mice ( 37,38 ). mRNA expression of some putative targets was found dysregulated in ob/ob mice, and the variation directions of some genes were in accordance with the changes of miRNAs, suggesting these genes to be Genes with double underlines are predicted by all the three algorithms, genes with single underlines are predicted by two algorithms, and genes without underline are predicted by only one algorithm. mRNA expression of genes in bold were found dysregulated ( ↑ , upregulated; ↓ , downregulated) in livers of ob/ob mice compared with those of C57BL/6 mice by previous studies ( 37,38 ). possible targets that may be repressed by miRNAs via mRNA degradation. Other predicted target genes that were not shown dysregulated in ob/ob mice on mRNA levels may possibly be regulated by miRNAs on protein levels. Citrate synthase, a Krebs tricarboxylic acid cycle enzyme that catalyzes the synthesis of citrate from oxaloacetate and acetyl-CoA, was predicted by all the three major algorithms to be a target gene of miR-122. As miR-122 was shown to be downregulated in livers of ob/ob mice in our study, we speculate citrate synthase to be upregulated in ob/ob mice. In accordance with this speculation, a previous study reported that citrate synthase activity was significantly elevated in livers of ob/ob mice compared with that of the wild-type control ( 57 ), indicating the targeting correlation between miR-122 and citrate synthase. Future studies are required to validate these potential target genes of the differentially expressed miRNAs and discover the downstream regulatory functions of miRNAs in NAFLD and hyperglycemia.
In conclusion, this study presents for the fi rst time an extensive analysis of miRNA expression in livers of ob/ob mice, STZ-induced diabetic mice, and normal C57BL/6 mice. The study has identifi ed a distinctive miRNA expression pattern in fatty livers of ob/ob mice, which could clearly distinguish ob/ob mice from either normal C57BL/6 mice or STZ-induced diabetic mice. The fi ndings of this study demonstrate a novel miRNA-based mechanism in regulating pathophysiological processes of NAFLD and provide potential targets for future clinical applications.
We thank Zhi-Xiang Fang for support in statistical analysis and Ke-Hui Wang, Chengfu Xu, Zheng Bian, and Jing Li for their technical help.