Figures
Abstract
Obesity is a major health concern worldwide which is associated with increased risk of chronic diseases such as metabolic syndrome, cardiovascular disease and cancer. The elucidation of the molecular mechanisms involved in adipogenesis and obesogenesis is of essential importance as it could lead to the identification of novel biomarkers and therapeutic targets for the development of anti-obesity drugs. MicroRNAs (miRNAs) have been shown to play regulatory roles in several biological processes. They have become a growing research field and consist of promising pharmaceutical targets in various fields such as cancer, metabolism, etc. The present study investigated the possible implication of miRNAs in adipose tissue during the development of obesity using as a model the C57BLJ6 mice fed a high-fat diet.
C57BLJ6 wild type male mice were fed either a standard (SD) or a high-fat diet (HFD) for 5 months. Total RNA was prepared from white adipose tissue and was used for microRNA profiling and qPCR.
Twenty-two of the most differentially expressed miRNAs, as identified by the microRNA profiling were validated using qPCR. The results of the present study confirmed previous results. The up-regulation of mmu-miR-222 and the down-regulation of mmu-miR-200b, mmu-miR-200c, mmu-miR-204, mmu-miR-30a*, mmu-miR-193, mmu-miR-378 and mmu-miR-30e* after HFD feeding has also been previously reported. On the other hand, we show for the first time the up-regulation of mmu-miR-342-3p, mmu-miR-142-3p, mmu-miR-142-5p, mmu-miR-21, mmu-miR-146a, mmu-miR-146b, mmu-miR-379 and the down-regulation of mmu-miR-122, mmu-miR-133b, mmu-miR-1, mmu-miR-30a*, mmu-miR-192 and mmu-miR-203 during the development of obesity. However, future studies are warranted in order to understand the exact role that miRNAs play in adipogenesis and obesity.
Citation: Chartoumpekis DV, Zaravinos A, Ziros PG, Iskrenova RP, Psyrogiannis AI, Kyriazopoulou VE, et al. (2012) Differential Expression of MicroRNAs in Adipose Tissue after Long-Term High-Fat Diet-Induced Obesity in Mice. PLoS ONE 7(4): e34872. https://doi.org/10.1371/journal.pone.0034872
Editor: Zhongjun Zhou, The University of Hong Kong, Hong Kong
Received: October 11, 2011; Accepted: March 10, 2012; Published: April 4, 2012
Copyright: © 2012 Chartoumpekis et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: The authors have no support or funding to report.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Obesity has emerged as an epidemic in modern societies [1] and has detrimental effects on the quality of life and the life expectancy of people as it increases the risk of cardiovascular disease, metabolic syndrome and cancer [2], [3].
Obesity is characterized by an excess of white adipose tissue mass. Two types of adipose tissue are found in mammals, white and brown adipose tissue (WAT and BAT). WAT is found below the skin (subcutaneous WAT) and inside the abdomen (abdominal WAT) while BAT is found in the interscapular region and axillae in rodents and in upper chest and neck in humans. The WAT has been characterized as an endocrine organ itself [4]. It produces endocrine acting peptides such as leptin and is metabolically important as WAT excess is associated with metabolic syndrome [5], [6], [7]. BAT is specialized in energy expenditure and its activation can counteract obesity [8], [9], [10]. Despite the research that has been made in the field of obesity and adipocyte differentiation, the exact mechanisms involved in these processes remain to be elucidated.
The non-coding microRNAs (miRNAs) which were discovered in the last decade [11], [12] and are currently known to post-transcriptionally regulate genes that are involved in various biological processes (cell growth, apoptosis, metabolism) [13], [14], have come in the spotlight of research in obesity and metabolic syndrome [15], [16], [17]. The existing studies on miRNAs and obesity indicate that these molecules may play an important role in the development of obesity and can possibly be promising targets for future therapeutic interventions.
A popular method to identify miRNAs that may play a role in adipogenesis and obesogenesis is the performance of miRNA arrays in animal (mouse) models of obesity (genetically or diet-induced) and in preadipocyte cell-culture models that are induced to differentiate into adipocytes [18]. In this way, miRNAs that change significantly during these processes are identified and then they can be checked with reference to their exact role in obesity or adipocyte differentiation. The majority of the existing studies with miRNA arrays in obesity and adipogenesis models [17], [19], [20] use the C57BLJ6 mouse after HFD feeding for 3 months as a model.
In the present study we investigated the changes that occur in miRNA expression in white adipose tissue after long-term (5 months) high fat diet (HFD) feeding in C57BLJ6 mice, using miRNA profiling and qPCR. We verified previously reported changes in miRNA expression. Moreover, we identified for the first time some new miRNAs that change significantly during long-term diet-induced obesity. These miRNAs could be proved to be valuable pharmaceutical targets for anti-obesity treatment.
Materials and Methods
Ethics statement
All animal procedures were approved by the institutional review board of the University of Patras Medical School (approval ID 15/2008) and were in accordance with EC Directive 86/609/EEC.
Mice
The C57BL6J wild type male mice (8 weeks old) were used. Eight mice were fed a standard diet (SD) (10 kcal% fat) and other eight were fed an HFD (60 kcal% fat; Research Diets, New Brunswick, NJ) for 5 months. Mice were housed in the animal facility of the University of Patras Medical School in temperature-, light- and humidity-controlled rooms with a 12-h light/dark cycle. Weight measurements were performed on a weekly basis to confirm that the HFD-fed mice gain more weight than the SD-fed ones. Evaluation of blood glucose levels and of insulin sensitivity was performed as described previously [21] at the beginning and in the end of the experiment, in both groups, so as to confirm that the HFD-fed mice developed an obese and insulin resistant phenotype.
White adipose tissue total RNA extraction
White adipose tissue (WAT) was excised from mice from epididymal fat pads and was submerged immediately in RNA later solution (Ambion, Foster City, CA). Total RNA was isolated from WAT from SD- or HFD-fed mice using the TRIzol reagent (Invitrogen, Carlsbad, CA), following the manufacturer's instructions. The individual RNA samples from SD- or HFD-fed mice were pooled for microarray experiments. The Reverse-Transcription PCR was performed in individual samples. The quality of the total RNA for microRNA profiling was assessed using an Agilent 2100 Bioanalyzer profile (Agilent Technologies, Santa Clara, CA). Total RNA concentration was measured with a Nanodrop instrument, calculating the ratio A260/A280.
microRNA profiling
One µg total RNA from sample and reference were labeled with Hy3™ and Hy5™ fluorescent label, respectively, using the miRCURY™ LNA Array power labeling kit (Exiqon,Denmark) following the procedure described by the manufacturer. The Hy3™-labeled samples and the Hy5™-labeled sample were mixed pair-wise and hybridized to the miRCURY™ LNA array version 10.0 (Exiqon,Denmark ), which contains capture probes targeting all miRNAs for all species registered in the miRBASE version 11.0 at the Sanger Institute. The hybridization was performed according to the miRCURY™ LNA array manual using a Tecan HS4800 hybridization station (Tecan, Austria). After hybridization, the microarray slides were scanned and stored in an ozone free environment (ozone level below 2.0 ppb) in order to prevent potential bleaching of the fluorescent dyes. The miRCURY™ LNA array microarray slides were scanned using the Agilent G2565BA Microarray Scanner System (Agilent Technologies, Inc., USA) and the image analysis was carried out using the ImaGene 8.0 software (BioDiscovery, Inc., USA).
Microarray Data Filtering, Background Correction and normalization
Filtering was performed based on the signal intensity. Spots with no signal above the background (flags 1 and 2) were detected and removed. Background correction was performed to remove non-biological contributions (“background”) to the measured signal, by subtracting the median global background from the median local background from the signal intensity. A threshold of 2 was set as cut-off, meaning that spot intensity for at least one channel should be twice as much as that of the background. Normalization was performed using the Lowess regression algorithm, in order to remove certain systematic biases from microarray data, such as dye effects or intensity dependence (Fig. S1 and Fig. S2). The positive effect from normalization is illustrated on each slide sheet with an M-A plot before and after normalization (Fig. S3). After normalization the spots appeared symmetrically scattered around the horizontal line (M = 0). Normalized data were extracted, pre-processed and sorted with Microsoft Excel ®.
Slide quality check using spike-in controls
Spike-in CV values were calculated between 32 replicates for each of the spike-in controls on each slide. CVs did not exceed 30%, thus reflecting lack of spatial effects on the array (smear across the slide) or misplacement of the grid for annotating the spots.
Clustering
We performed clustering of the miRNAs using four methods: Hierarchical clustering (HCL), k-means clustering, Principle Components Analysis (PCA) and TIGR clustering (terrain classification) using 3-dimensional representation of the clustering results.
GEO accession numbers
Our miRNA profiling data were MIAME compliant and deposited at the Gene Expression Omnibus (GEO) repository, with accession numbers from GSM802367 through GSM802369 (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=btghlsqqmoessvc&acc=GSE32419).
Reverse-Transcription (RT) and qPCR validation
The microRNA in the total RNA sample was converted to cDNA by reverse transcription using the miRCURY™ LNA microRNA PCR system First-strand cDNA synthesis kit and a microRNA-specific primer. The cDNA was amplified by real-time PCR using SYBR® Green master mix and LNA™ microRNA-specific primers on a Step One Plus Instrument (Applied Biosystems, Foster City, CA) using the designated PCR optical tubes. The following 22 murine microRNAs were selected for qPCR validation of their expression: mmu-miR-1, mmu-miR-21, mmu-miR-30a*, mmu-miR-30e*, mmu-miR-122, mmu-miR-130a, mmu-miR-133b, mmu-miR-141, mmu-miR-142-3p, mmu-miR-142-5p, mmu-miR-146a, mmu-miR-146b, mmu-miR-192, mmu-miR-193a-3p, mmu-miR-200b, mmu-miR-200c, mmu-miR-203, mmu-miR-204, mmu-miR-222, mmu-miR-342-3p, mmu-miR-378 and mmu-miR-379. All runs included no template and RT-minus controls. A 5-fold serial dilution of pooled cDNA samples was generated for each assay, in order to calculate the PCR efficiency. All reactions were performed in triplicates. Before proceeding with data normalization, a technical quality assessment was performed based on results of the melting curve, serial dilution curve and no-template-controls. The stability and ranking of the endogenous controls was also calculated with the SLqPCR algorithim. For every endogenous control gene, the pair-wise variation with all other endogenous controls was determined as a gene stability measurement M. An M value below 1.5 was recommended and genes with expression stability above 1.5 were considered unstable across the samples and unsuitable for endogenous controls in this experiment. U6 snRNA, mmu-191, mmu-423-5p, mmu-361 and mmu-103 were considered acceptable and used for normalization [22]. The gene 5S rRNA was unstable across the samples and therefore excluded from the normalization (Fig. S4). Relative expression was performed using the ΔΔCt method [23], [24]. Data were standardized by log2 transformation.
Gene Ontology (GO) and enrichment analysis
Gene Ontology (GO) analysis is essential in the deduction of conclusions from microarray data. GO is a database with curated annotations for known genes i.e. gene biological processes, molecular functions and cellular components. Firstly, using the miRWalk software [25] the validated target genes of the confirmed differentially expressed miRNAs were identified. Then, GO analysis was performed for the validated miRNA target genes, using the Genesis 1.7.2 software [26] and the WebGestalt web-tool (http://bioinfo.vanderbilt.edu/webgestalt) [27]). All gene definitions and functions were based on the National Institute of Health databases (http://www.ncbi.nlm.nih.gov/sites/entrez/). Further enrichment analysis for the validated miRNA targets was also investigated, using the WebGestalt web-tool (http://bioinfo.vanderbilt.edu/webgestalt) [27]. The hypergeometric test, with Bonferroni correction was used for enrichment evaluation analysis. The R function adjP was used in order to adjust the nominal p values of the large number of categories at the same time. The significance level for the adjusted p-value was set at 0.01 and the number of minimum genes for a category was set at 2.
Statistical analysis
Normality of the data distribution was checked by Kolmogorov-Smirnov test. Differences in the expression levels between SD- and HFD-fed mice in WAT were explored using the Mann-Whitney U test. Numerical values are expressed as the mean±standard deviation. Statistical significance was set at the 95% level (p<0.05). In weight and blood glucose measurements in mice, as well as in qPCR, the Student t-test or the one-way ANOVA followed by Tukey's test was used for statistical analysis. These data were expressed as the mean±SEM or mean±standard deviation, as described in the corresponding figure legends. Statistical significance was set at the 95% level (p<0.05). The statistical package GraphPad Prism was used (GraphPad Software, La Jolla, CA).
Results
In the present study, we examined the expression levels of 530 mmu-miRNAs in adipose tissue received from SD-fed mice and compared them with the corresponding expression in adipose tissue from HFD-fed mice for a 5-month period.
Evaluation of the obese and insulin resistant phenotype in mice after 5 months on HFD
At the beginning of the experiment, mice were age- and weight-matched, and thus had no statistically significant differences in their weights (Fig. S5). After 5 months on HFD or SD, the HFD-fed mice were significantly heavier (p<0.0001) than the SD-fed mice, as expected (Fig. S5). Moreover, although the two mouse groups did not have any differences in glucose tolerance and insulin resistance at the beginning of the experiment, the HFD-fed mice developed a significantly less glucose tolerant and more insulin resistant phenotype than the SD-fed mice after 5 months on the diet, as revealed by the IPGTT (Intraperitoneal Glucose Tolerance Test) and by the IPITT (Intraperitoneal Insulin Tolerance Test), correspondingly (Fig. S6).
MicroRNA profiling
Based on our miRCURY™ LNA array microarray data, we identified 26 differentially expressed miRNAs in adipose tissue between SD and HFD mice. The following miRNAs were found to be up-regulated in WAT after HFD feeding: mmu-miR-342-3p, mmu-miR-222, mmu-miR-221, mmu-miR-142-3p, mmu-miR-142-5p, mmu-miR-21, mmu-miR-335-5p, mmu-miR-146a, mmu-miR-146b, mmu-miR-647* and mmu-miR-379. On the contrary, the following miRNAs were down-regulated in WAT after HFD feeding: mmu-miR-141, mmu-miR-200a, mmu-miR-200b, mmu-miR-200c, mmu-miR-122, mmu-miR-204, mmu-miR-133b, mmu-miR-1, mmu-miR-30a*, mmu-miR-130a, mmu-miR-192, mmu-miR-193a-3p, mmu-miR-203, mmu-miR-378 and mmu-miR-30e*. The log2 fold change between SD and HFD mice is depicted in Figure 1.
The SD #1 and #2 indicate technical replicates. The diagram shows the top ranking 26 differentially expressed miRNAs based on the comparison samples versus control. Each bar represents the fold change between the sample mice versus the control mouse. Since slide no. 3 is a dye swap (control vs. samples) the values have been reversed to fit slides no. 1 and 2. A fold change >1 indicates up-regulation after HFD feeding, whereas a fold change <1 indicates down-regulation after HFD feeding. SD; standard diet, HFD; high-fat diet.
Clustering
We performed Hierarchical clustering using Average Linkage Clustering and Euclidian Distance. The condensed tree and cluster image are depicted in Figure 2. Two big branches could be noticed, and each one was composed of smaller ones. The first branch (5 miRNA groups) contained mmu-miRNAs which were under-expressed in the SD vs. the HFD mice. The second branch was composed of 2 big groups, partitioned into 4 sub-clusters. On the contrary, the first group of the second big branch contained miRNAs which were over-expressed in the SD vs. the HFD mice. Finally, the second group of the second big branch contained miRNAs which were over-expressed both in SD and HFD mice. Each group is distinct by a different coloring.
Color saturation is directly proportional to measured expression ratio magnitude. Rows represent individual probe set. Columns represent experimental sample. Red bars indicate high expression. Green bars indicate low expression. Slides 1 and 2: control (SD) mice, slide3: sample (HFD) mice. SD; standard diet. HFD; high-fat diet.
Similar clustering was performed on the principal components (PCs) outputted by Principal Component Analysis (PCA) on the same original data (where the gene expression measurements are the variables and the samples are the observations). In the PCA method, the observed samples were represented as a linear combination of the PCs with associated gene scores. PCA showed an even larger number of miRNA groups. Thirteen clusters could be observed in a three-dimensional clustering (Fig. 3). The clustering results obtained by HCL and PCA are depicted side-by-side in Figure 2 for an immediate comparison between them. K-means clustering identified 10 clusters (Fig. 4). The number of genes in each cluster, the (%) share of genes in each cluster, the average distance from the cluster mean, the next neighbor variance, as well as the within cluster variance, are depicted in detail, in Table 1.
The observed samples were represented as a linear combination of the principal components with associated gene scores. The three-dimensional PCA clustering revealed 13 miRNA clusters.
This non-hierarchial method initially takes the number of components of the population equal to the final required number of clusters. In this step itself the final required number of clusters is chosen such that the points are mutually farthest apart. Next, it examines each component in the population and assigns it to one of the clusters depending on the minimum distance. The centroid's position is recalculated everytime a component is added to the cluster and this continues until all the components are grouped into the final required number of clusters. A. K-means clustering of profiles and B. Centroids.
A three-dimensional TIGR clustering analysis was also performed. The terrain (map) analysis, depicting the gene clusters, as well as the links among the clustered genes, is depicted in Figure 5. The predictive power of the K-means algorithm was estimated using a figure of merit (FOM). A FOM value vs. number of clusters was computed by removing each sample in turn from the data set, clustering genes based on the remaining data, and calculating the fit of the withheld sample to the clustering pattern obtained from the other samples (Fig. 6).
A. Each miRNA cluster is depicted with a different color. B. The correlations among the various miRNAs are depicted with lines.
A figure of merit is an estimate of the predictive power of a clustering algorithm. The lower the adjusted FOM value is, the higher the predictive power of the algorithm. The value of the adjusted FOM for the k-means run decreases steeply until the number of clusters reaches 7, after which it levels out. This suggests that, for this data set, k-means performs optimally for 7 clusters and that any additional clusters produced will not add to the predictive value of the algorithm.
qPCR validation of expression
Twenty-two of the most differentially expressed miRNAs after HFD feeding, as revealed by our microRNA profiling, were also tested by qPCR. After normalizing the expression levels with the corresponding geometric mean value of the reference genes U6 snRNA, mmu-miR-191, mmu-miR-423-5p, mmu-miR-361 and mmu-miR-103, the samples were checked for outliers to be excluded. The qPCR verified the expression profile of the 22 DE miRNAs. The fold changes were calculated using the ΔΔCt comparative quantification method. Data were corrected for assay specific PCR efficiency (microRNA or endogenous control) and normalized with a normalization factor based on the five validated endogenous controls. Mmu-miR-141 was not detected in the SD mice (control), hence no fold change was calculated. The relative miRNA expression levels for each miRNA are depicted in Figure 7. Moreover, we performed HCL and PCA analysis for the qPCR validation data of 18 miRNAs (Fig. 8). The clustering results were similar to those acquired by the microarray data.
The miRNA levels were measured by quantitative RT-PCR in white adipose tissue from mice fed a standard or a high-fat diet for 5 months (n = 8 for each diet type). The RT-PCR was performed in triplicate wells for each individual sample. Bars show means±standard deviation. * p<0.0001, † p<0.001. SD; standard diet, HFD; high-fat diet.
HCL (A) and PCA analysis (B) for 18 mmu-miRs validated by qPCR. The clustering results were similar to those acquired by the miRNA profiling.
Correlation between microrray and qPCR results
We identified an agreement between the microarray and the qPCR results, showing the relative expression of array data compared to the relative expression of qPCR data (Fig. 9A). The relative expressions were based on three technical replicates. The Pearson correlation value between microarray and qPCR data was 0.888 (Fig. 9B).
A. The fold changes were calculated using the ΔΔCt comparative quantification method. Mmu-miR-141 was not detected in “control mice” (SD), hence no fold change was calculated. We identified an agreement between the microarrays and the qPCR results. B. Scatterplot and Correlation between qPCR and microarrays (Pearson's correlation = 0.888).
Gene Ontology (GO)
GO analysis was performed in order to approach the functionality of the target genes of the studied miRNAs. Firstly, the validated target genes of the differentially expressed miRNAs were identified using the miRWalk algorithms shown in Table S1, which is divided into the up- and the down-regulated miRNAs. GO analysis was initially performed for the target genes of the up-regulated miRNAs and then for the down-regulated ones. The genes were initially separated into three groups, based on their molecular function, biological process and cellular component. The enriched processes of the up-regulated miRNA targets, along with their statistical significance, are depicted in Figure 10. Cell differentiation (adjP = 1.41e-30) and metabolic processes (regulation of metabolic process, adjP = 1.85e-25) are among the most significant biological processes related to the development of obesity. Table S2 includes the genes that participate in all the depicted enriched processes. Figure 11 exhibits the enriched processes for the target genes of the down-regulated miRNAs. Again the enrichment of processes related to obesity, such as regulation of metabolic process (adjP = 3.21e-25) or cell differentiation (adjP = 4.89e-35), is high. Table S3 depicts the genes that participate in all the depicted enriched processes.
The most significant functions of the enriched miRNAs are highlighted in red color.
The most significant functions of the enriched miRNAs are highlighted in red color.
Discussion
The current study investigated the changes that occur in miRNA expression in the adipose tissue, in mice after a long-term (5-month) HFD-induced obesity. For this purpose, we performed miRNA profiling and qPCR validation of our results. The expression of some miRNAs was found either increased or decreased in HFD-fed mice, compared to their SD-fed counterparts. Some of the results are in accordance with previous studies, such as the up-regulation of mmu-miR-221 and mmu-miR222 cluster and the down-regulation of the mmu-miR-200 family, as well as of mmu-miR-204, mmu-miR-30a*, mmu-miR-193, mmu-miR-378 and mmu-miR-30e*. These findings have also been confirmed in other independent studies. Their importance lies in that they suggest new miRNA players in the field of obesity. However, further in vivo and in vitro experimental investigation is warranted in order to describe them as potential pharmaceutical targets.
The rest of the differentially expressed miRNAs, after HFD feeding, has not been previously reported, regarding adipogenesis and obesogenesis. These include on the one hand the up-regulated miRNAs: mmu-miR-342-3p, mmu-miR-142-3p, mmu-miR-142-5p, mmu-miR-21, mmu-miR-335-5p, mmu-miR-146a, mmu-miR-146b, mmu-miR-674 and mmu-miR-379; and on the other hand the down-regulated ones after HFD-induced obesity: mmu-miR-122, mmu-miR-133p, mmu-miR-1, mmu-miR-30a, mmu-miR-192 and mmu-miR-203.
To get a deeper insight into the potential functional significance of the differentially expressed miRNAs in white adipose tissue, their validated target genes were identified using the miRWalk algorithm. Then, Gene Ontology analysis classified the potential enriched functions of these genes in groups. A single miRNA is likely to regulate many genes and each target gene may be regulated by more than one miRNAs, as revealed by the miRWalk algorithm. Therefore, the potential network of miRNAs and genes that are involved in adipose tissue during HFD-induced obesity seems to be highly complicated. Nevertheless, special focus was given to certain highly enriched functions related to obesity and adipogenesis.
Firstly, a significant function of the differentially expressed miRNA targets is cell differentiation and the regulation of cell differentiation (Figs.10, 11). This is linked to the fact that obesity is characterized of both hypertrophy (enlargement of the existing adipocytes) and hyperplasia (adipogenesis or adipocyte differentiation) [28]. It is also important to point out that miRNAs that are induced during adipogenesis, are usually repressed during obesity development [19].
The mmu-miR-142-5p-regulated genes, Bmp-4 and Fgf10, are involved in adipogenesis regulation [29]. Furthermore, Akt1, a target of mmu-miR-142-3p, is also involved in this regulation [30]. Both mmu-miR-142-3p and mmu-miR-142-5p derive from the stem loop precursor of mmu-miR-142. This is the first study to connect mmu-miR-142 with obesity; however, further experiments are warranted to elucidate its role in this process.
Mmu-miR-21 and mmu-miR-146a are also potentially involved in adipocyte differentiation by targeting C/EBP beta [31] and ApoE, respectively [32]. Mmu-miR-146a, which shares 91% homology with mmu-miR-146b in a length of 22 nt [33], is mainly implicated in inflammatory response by being induced by NFκB [34] and by targeting interferon gamma (as revealed by miRWalk in this study). In contrast with mmu-miR-21 for which there are studies describing it to increase during adipocyte differentiation [35], [36], no reference is available regarding mmu-miR-146a, mmu-miR-146b and obesity. Taking into account the up-regulation of mmu-miR-21, mmu-miR-146a and mmu-miR-146b after HFD-induced obesity in the present study, these miRNAs are suggested as a research subject for future studies in obesity.
Mmu-miR-222 also belongs to the miRNA group being highly enriched during adipogenesis. It has already been shown to be up-regulated in obesity and down-regulated in adipogenesis [19], [37]. Hif-1α and TNF-α are two potential mmu-miR-222 targets participating in the regulation of cell differentiation. Both genes have been described to be indirectly implicated in adipogenesis [38], [39].
Mmu-miR-342-3p and mmu-miR-379 were both found to be up-regulated after HFD-induced obesity in the present study. These miRNAs have not been described before with relevance to adipogenesis and obesity. A possible link between mmu-miR-342-3p and obesity is the glucagon that is a potential target of this miRNA and is implicated in obesogenesis [40]. Further studies are needed to clarify any potential role of these miRNAs in these processes.
Next, focus was given on the enriched process “regulation of metabolic process”. Mmu-miR-200b and mmu-miR-200c(members of the mmu-miR-200 family), mmu-miR-203 and mmu-miR-192 target Zeb1 and Zeb2 that regulate epithelial to mesenchymal transition [41] and have recently been implicated in adipogenesis and obesity [42]. The down-regulation of mmu-miR-200b and mmu-miR-200c after HFD-induced obesity, is in accordance with a previous study which showed that the mmu-miR-200 family promotes adipogenesis [43]. The down-regulation of mmu-miR-203 and mmu-miR-192 is described for the first time and warrants further elucidation.
Of the down-regulated miRNAs after HFD-induced obesity, the putative targets of mmu-miR-1, mmu-miR-204 and mmu-miR-133b, are implicated in cell differentiation. Mmu-miR-1 has mainly been characterized as a muscle-specific miRNA [44] and it has been found to be expressed in brown adipocytes [45]. This is in accordance with the theory that brown adipocytes and muscle share common precursor cells. Mmu-miR-1 targets UCP-1, which uncouples oxidative phosphorylation from the production of ATP in brown adipocytes and energy is dissipated as heat [46]. The expression of mmu-miR-1 that we detected in white adipose tissue, can be due to the presence of brown adipocytes in white adipose tissue [47] and may be important in the metabolic adaptation that occurs in the mouse during HFD-feeding. Mmu-miR-204 and mmu-miR-133b target Runx2 which promotes osteogenesis and represses adipogenesis [48]. Mmu-miR-204 has already been described to promote adipogenesis and repress osteogenesis through Runx2 repression [49].
The down-regulation of mmu-miR-30a*, mmu-miR-30e*, mmu-miR-193 and mmu-miR-378 during HFD-induced obesity is consistent with previous studies [19], [37], [50]. This indicates that these miRNAs should be further studied, so as to check if they can successfully be targets for obesity treatment. On the other hand, the down-regulation of mmu-miR-122 and mmu-miR-130a during HFD-induced obesity that we report has not been described before. Mmu-miR-122 is mainly expressed in the liver and makes up for the 70% of all the liver miRNAs [51]. It is known to be implicated in cholesterol biosynthetic pathway and fatty acid metabolism [52]. Its expression in adipose tissue has been found to be 200-fold less than in liver [53] and further characterization of its potential role in obesogenesis is warranted. Similarly, mmu-miR-130a is known to be a pro-angiogenic miRNA and its role in adipose tissue has not been described before [54]. Taking these data into account, it seems that miRNAs play a regulatory role in the development of obesity. The present study confirms the results of previous reports with reference to the differential expression of miRNAs, after HFD feeding and during adipogenesis. Moreover, we present novel miRNAs which are involved in the aforementioned processes. Their exact role remains to be further investigated, so as to come to a conclusion regarding their usefulness as therapeutical targets for the prevention and treatment of obesity, diabetes type 2 and metabolic syndrome.
Supporting Information
Figure S1.
Densities of individual-channel intensities for two-color microarray data for all slides. The density plot shows the effect of the normalization as the signal distribution of the two channels is more similar after normalization than before. (Hy3; green, Hy5; red).
https://doi.org/10.1371/journal.pone.0034872.s001
(TIF)
Figure S2.
Box plots illustrating the distribution of signals measured on the capture probes before and after normalization. The box-plots show that the majority of capture probes have similar signals in both channels suggesting that the miRNAs have similar expression levels in the paired samples on each slide (Capture probes with a log2 median ratio of “0” on the Y-axis correspond to miRNAs that are equally expressed in the 2 different samples on this slide). The lower boundary of the box indicates the 25th percentile, the line within the box shows the median, and the upper edge of the box marks the 75th percentile. Whiskers above and below each box indicate the 95th and 5th percentiles. All data points that lie outside the 5th and 95th percentiles are shown as symbols.
https://doi.org/10.1371/journal.pone.0034872.s002
(TIF)
Figure S3.
The positive effect from normalization is illustrated on each slide with an M-A plot before and after normalization. After normalization the spots appear symmetrically scattered around the horizontal line M = 0. The difference between the two channels (M) is now independent of the average intensity level of the two channels. (slide 1-control 1, A; control 2-sample 2, B; slide 3-sample, C)
https://doi.org/10.1371/journal.pone.0034872.s003
(TIF)
Figure S4.
Gene stability measurement (M) for endogenous control genes. For every endogenous control gene, the pair-wise variation with all other endogenous controls was determined as a gene stability measurement M. An M value below 1.5 was recommended and genes with expression stability above 1.5 were considered unstable across the samples and unsuitable for endogenous controls in this experiment. U6 snRNA, mmu-191, mmu-423-5p, mmu-361 and mmu-103 were considered acceptable and used for normalization. The gene 5S rRNA was unstable across the samples and therefore excluded from the normalization.
https://doi.org/10.1371/journal.pone.0034872.s004
(TIF)
Figure S5.
Weight measurements of mice before and after HFD feeding (5 months). Wild type mice were fed a SD or a HFD diet for 5 months. At the beginning of the experiment the two groups of mice had no difference in their weights. After 5 months on SD or HFD, the mice of the HFD group were significantly heavier than those of the SD group, as expected. Data show mean±SEM values. N = 8 for each diet type. *p<0.0001. SD; standard diet, HFD; high-fat diet.
https://doi.org/10.1371/journal.pone.0034872.s005
(TIF)
Figure S6.
Wild type mice fed a HFD for 5 months become more insulin resistant and less glucose tolerant than the SD-fed ones. Results from glucose and insulin tolerance tests are shown at the beginning of the experiment (day 0) and after 5 months on SD or HFD feeding. Data show means±SEM. N = 8 for each diet type. *p<0.01, **p<0.001. IPGTT: intraperitoneal glucose tolerance test, IPITT: intraperitoneal insulin tolerance test. SD; standard diet, HFD; high-fat diet.
https://doi.org/10.1371/journal.pone.0034872.s006
(TIF)
Table S1.
Upregulated and downregulated miRNAs with their corresponding validated target genes based on miRWalk algorithms
https://doi.org/10.1371/journal.pone.0034872.s007
(XLS)
Table S2.
Target genes of the upregulated miRNAs and the enriched processes in which they participate.
https://doi.org/10.1371/journal.pone.0034872.s008
(XLS)
Table S3.
Target genes of the downregulated miRNAs and the enriched processes in which they participate.
https://doi.org/10.1371/journal.pone.0034872.s009
(XLS)
Author Contributions
Conceived and designed the experiments: DVC PGZ IGH. Performed the experiments: DVC PGZ RPI IGH. Analyzed the data: AZ DVC IGH. Contributed reagents/materials/analysis tools: AZ AIP VEK. Wrote the paper: DVC AZ IGH AIP VEK.
References
- 1. James PT, Leach R, Kalamara E, Shayeghi M (2001) The worldwide obesity epidemic. Obes Res 9: Suppl 4228S–233S.
- 2. Calle EE, Kaaks R (2004) Overweight, obesity and cancer: epidemiological evidence and proposed mechanisms. Nat Rev Cancer 4: 579–591.
- 3. Bergman RN, Kim SP, Hsu IR, Catalano KJ, Chiu JD, et al. (2007) Abdominal obesity: role in the pathophysiology of metabolic disease and cardiovascular risk. Am J Med 120: S3–8; discussion S29–32.
- 4. Kershaw EE, Flier JS (2004) Adipose tissue as an endocrine organ. J Clin Endocrinol Metab 89: 2548–2556.
- 5. Despres JP (1993) Abdominal obesity as important component of insulin-resistance syndrome. Nutrition 9: 452–459.
- 6. Despres JP, Lemieux I (2006) Abdominal obesity and metabolic syndrome. Nature 444: 881–887.
- 7. Wajchenberg BL (2000) Subcutaneous and visceral adipose tissue: their relation to the metabolic syndrome. Endocr Rev 21: 697–738.
- 8. Yoshida T, Sakane N, Wakabayashi Y, Umekawa T, Kondo M (1994) Anti-obesity and anti-diabetic effects of CL 316,243, a highly specific beta 3-adrenoceptor agonist, in yellow KK mice. Life Sci 54: 491–498.
- 9. Umekawa T, Yoshida T, Sakane N, Saito M, Kumamoto K, et al. (1997) Anti-obesity and anti-diabetic effects of CL316,243, a highly specific beta 3-adrenoceptor agonist, in Otsuka Long-Evans Tokushima Fatty rats: induction of uncoupling protein and activation of glucose transporter 4 in white fat. Eur J Endocrinol 136: 429–437.
- 10. Clarke IJ, Henry BA (2010) Targeting energy expenditure in muscle as a means of combating obesity. Clin Exp Pharmacol Physiol 37: 121–124.
- 11. Couzin J (2008) MicroRNAs make big impression in disease after disease. Science 319: 1782–1784.
- 12. Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T (2001) Identification of novel genes coding for small expressed RNAs. Science 294: 853–858.
- 13. Osada H, Takahashi T (2007) MicroRNAs in biological processes and carcinogenesis. Carcinogenesis 28: 2–12.
- 14. He L, Hannon GJ (2004) MicroRNAs: small RNAs with a big role in gene regulation. Nat Rev Genet 5: 522–531.
- 15. McGregor RA, Choi MS (2011) microRNAs in the regulation of adipogenesis and obesity. Curr Mol Med 11: 304–316.
- 16. Kolfschoten IG, Roggli E, Nesca V, Regazzi R (2009) Role and therapeutic potential of microRNAs in diabetes. Diabetes Obes Metab 11: Suppl 4118–129.
- 17. Heneghan HM, Miller N, Kerin MJ (2010) Role of microRNAs in obesity and the metabolic syndrome. Obes Rev 11: 354–361.
- 18. Alexander R, Lodish H, Sun L (2011) MicroRNAs in adipogenesis and as therapeutic targets for obesity. Expert Opin Ther Targets 15: 623–636.
- 19. Xie H, Lim B, Lodish HF (2009) MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity. Diabetes 58: 1050–1057.
- 20. Bengestrate L, Virtue S, Campbell M, Vidal-Puig A, Hadaschik D, et al. (2011) Genome-wide profiling of microRNAs in adipose mesenchymal stem cell differentiation and mouse models of obesity. PLoS One 6: e21305.
- 21. Chartoumpekis DV, Ziros PG, Psyrogiannis AI, Papavassiliou AG, Kyriazopoulou VE, et al. (2011) Nrf2 Represses FGF21 During Long-Term High-Fat Diet-Induced Obesity in Mice. Diabetes 60: 2465–2473.
- 22. Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, et al. (2002) Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol 3: RESEARCH0034.
- 23. Livak KJ, Schmittgen TD (2001) Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods 25: 402–408.
- 24. Radojicic J, Zaravinos A, Vrekoussis T, Kafousi M, Spandidos DA, et al. (2011) MicroRNA expression analysis in triple-negative (ER, PR and Her2/neu) breast cancer. Cell Cycle 10: 507–517.
- 25. Dweep H, Sticht C, Pandey P, Gretz N (2011) miRWalk–database: prediction of possible miRNA binding sites by “walking” the genes of three genomes. J Biomed Inform 44: 839–847.
- 26. Sturn A, Quackenbush J, Trajanoski Z (2002) Genesis: cluster analysis of microarray data. Bioinformatics 18: 207–208.
- 27. Zhang B, Kirov S, Snoddy J (2005) WebGestalt: an integrated system for exploring gene sets in various biological contexts. Nucleic Acids Res 33: W741–748.
- 28. Nishimura S, Manabe I, Nagasaki M, Hosoya Y, Yamashita H, et al. (2007) Adipogenesis in obesity requires close interplay between differentiating adipocytes, stromal cells, and blood vessels. Diabetes 56: 1517–1526.
- 29. Gimble JM, Morgan C, Kelly K, Wu X, Dandapani V, et al. (1995) Bone morphogenetic proteins inhibit adipocyte differentiation by bone marrow stromal cells. J Cell Biochem 58: 393–402.
- 30. Yun SJ, Kim EK, Tucker DF, Kim CD, Birnbaum MJ, et al. (2008) Isoform-specific regulation of adipocyte differentiation by Akt/protein kinase Balpha. Biochem Biophys Res Commun 371: 138–143.
- 31. Tanaka T, Yoshida N, Kishimoto T, Akira S (1997) Defective adipocyte differentiation in mice lacking the C/EBPbeta and/or C/EBPdelta gene. EMBO J 16: 7432–7443.
- 32. Huang ZH, Gu D, Mazzone T (2009) Role of adipocyte-derived apoE in modulating adipocyte size, lipid metabolism, and gene expression in vivo. Am J Physiol Endocrinol Metab 296: E1110–1119.
- 33. Lagos-Quintana M, Rauhut R, Yalcin A, Meyer J, Lendeckel W, et al. (2002) Identification of tissue-specific microRNAs from mouse. Curr Biol 12: 735–739.
- 34. Taganov KD, Boldin MP, Chang KJ, Baltimore D (2006) NF-kappaB-dependent induction of microRNA miR-146, an inhibitor targeted to signaling proteins of innate immune responses. Proc Natl Acad Sci U S A 103: 12481–12486.
- 35. Kim YJ, Hwang SJ, Bae YC, Jung JS (2009) MiR-21 regulates adipogenic differentiation through the modulation of TGF-beta signaling in mesenchymal stem cells derived from human adipose tissue. Stem Cells 27: 3093–3102.
- 36. Keller P, Gburcik V, Petrovic N, Gallagher IJ, Nedergaard J, et al. (2011) Gene-chip studies of adipogenesis-regulated microRNAs in mouse primary adipocytes and human obesity. BMC Endocr Disord 11: 7.
- 37. Ortega FJ, Moreno-Navarrete JM, Pardo G, Sabater M, Hummel M, et al. (2010) MiRNA expression profile of human subcutaneous adipose and during adipocyte differentiation. PLoS One 5: e9022.
- 38. Xu H, Sethi JK, Hotamisligil GS (1999) Transmembrane tumor necrosis factor (TNF)-alpha inhibits adipocyte differentiation by selectively activating TNF receptor 1. J Biol Chem 274: 26287–26295.
- 39. Yun Z, Maecker HL, Johnson RS, Giaccia AJ (2002) Inhibition of PPAR gamma 2 gene expression by the HIF-1-regulated gene DEC1/Stra13: a mechanism for regulation of adipogenesis by hypoxia. Dev Cell 2: 331–341.
- 40. Day JW, Ottaway N, Patterson JT, Gelfanov V, Smiley D, et al. (2009) A new glucagon and GLP-1 co-agonist eliminates obesity in rodents. Nat Chem Biol 5: 749–757.
- 41. Gregory PA, Bert AG, Paterson EL, Barry SC, Tsykin A, et al. (2008) The miR-200 family and miR-205 regulate epithelial to mesenchymal transition by targeting ZEB1 and SIP1. Nat Cell Biol 10: 593–601.
- 42. Saykally JN, Dogan S, Cleary MP, Sanders MM (2009) The ZEB1 transcription factor is a novel repressor of adiposity in female mice. PLoS One 4: e8460.
- 43. Kennell JA, Gerin I, MacDougald OA, Cadigan KM (2008) The microRNA miR-8 is a conserved negative regulator of Wnt signaling. Proc Natl Acad Sci U S A 105: 15417–15422.
- 44. Yang B, Lin H, Xiao J, Lu Y, Luo X, et al. (2007) The muscle-specific microRNA miR-1 regulates cardiac arrhythmogenic potential by targeting GJA1 and KCNJ2. Nat Med 13: 486–491.
- 45. Walden TB, Timmons JA, Keller P, Nedergaard J, Cannon B (2009) Distinct expression of muscle-specific microRNAs (myomirs) in brown adipocytes. J Cell Physiol 218: 444–449.
- 46. Feldmann HM, Golozoubova V, Cannon B, Nedergaard J (2009) UCP1 ablation induces obesity and abolishes diet-induced thermogenesis in mice exempt from thermal stress by living at thermoneutrality. Cell Metab 9: 203–209.
- 47. Cousin B, Cinti S, Morroni M, Raimbault S, Ricquier D, et al. (1992) Occurrence of brown adipocytes in rat white adipose tissue: molecular and morphological characterization. J Cell Sci 103(Pt 4): 931–942.
- 48. Enomoto H, Furuichi T, Zanma A, Yamana K, Yoshida C, et al. (2004) Runx2 deficiency in chondrocytes causes adipogenic changes in vitro. J Cell Sci 117: 417–425.
- 49. Huang J, Zhao L, Xing L, Chen D (2010) MicroRNA-204 regulates Runx2 protein expression and mesenchymal progenitor cell differentiation. Stem Cells 28: 357–364.
- 50. Gerin I, Bommer GT, McCoin CS, Sousa KM, Krishnan V, et al. (2010) Roles for miRNA-378/378* in adipocyte gene expression and lipogenesis. Am J Physiol Endocrinol Metab 299: E198–206.
- 51. Chang J, Nicolas E, Marks D, Sander C, Lerro A, et al. (2004) miR-122, a mammalian liver-specific microRNA, is processed from hcr mRNA and may downregulate the high affinity cationic amino acid transporter CAT-1. RNA Biol 1: 106–113.
- 52. Esau C, Davis S, Murray SF, Yu XX, Pandey SK, et al. (2006) miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab 3: 87–98.
- 53. Zhao E, Keller MP, Rabaglia ME, Oler AT, Stapleton DS, et al. (2009) Obesity and genetics regulate microRNAs in islets, liver, and adipose of diabetic mice. Mamm Genome 20: 476–485.
- 54. Chen Y, Gorski DH (2008) Regulation of angiogenesis through a microRNA (miR-130a) that down-regulates antiangiogenic homeobox genes GAX and HOXA5. Blood 111: 1217–1226.