- Open Access
Circulating microRNAs are deregulated in overweight/obese children: preliminary results of the I.Family study
© The Author(s) 2016
Received: 1 December 2015
Accepted: 8 March 2016
Published: 21 March 2016
MicroRNAs (miRNAs) are small non-coding RNAs involved in the modulation of gene expression and in the control of numerous cell functions. Alterations of miRNA patterns frequently occur in cancer and metabolic disorders, including obesity. Recent studies showed remarkable stability of miRNAs in both plasma and serum making them suitable as potential circulating biomarkers for a variety of diseases and conditions.
The aim of this study was to assess the profile of circulating miRNAs expressed in plasma samples of overweight or obese (OW/Ob) and normal weight (NW) prepubertal children from a European cohort (www.ifamilystudy.eu). The project, aimed to assess the determinants of eating behavior in children and adolescents of eight European countries, is built on the IDEFICS cohort (www.ideficsstudy.eu), established in 2006. Among the participants of the I.Family Italian Cohort, ten OW/Ob (age 10.7 ± 1.5 years, BMI 31.6 ± 4.3 kg/m2) and ten NW (age 10.5 ± 2.7 years, BMI 16.4 ± 1.7 kg/m2) children were selected for the study. Gene arrays were employed to differentially screen the expression of 372 miRNAs in pooled plasma samples. Deregulated miRNAs (p < 0.05) were further validated in the individual samples using a real-time PCR (RT-qPCR) approach.
Using a significance threshold of p < 0.05 and a fold-change threshold of ± 4.0, we preliminarily identified in the pooled samples eight miRNAs that differed between the OW/Ob and NW groups. The validation by RT-qPCR in the individual plasma samples showed a twofold upregulation of miR-31-5p, a threefold upregulation of miR-2355-5p, and a 0.5-fold downregulation of miR-206 in OW/Ob as compared with NW. The molecular functions of these differentially expressed plasma miRNAs as well as their expected mRNA targets were predicted by bioinformatics tools.
This pilot study shows that three circulating miRNAs are differentially regulated in OW/Ob as compared with NW children. Although causal pathways cannot be firmly inferred by these results, that deserve confirmation in larger samples, it is conceivable that circulating miRNAs may be novel biomarkers of obesity and related metabolic disturbances.
The obesity epidemic represents a major health challenge worldwide since it is associated with severe adverse consequences for human health [1–3]. Excess energy intake and lack of physical activity are the main factors that drive the development of obesity. Individual traits, such as the genetic profile, likewise contribute to obesity. A quantity of studies has found a number of genes associated with obesity and obesity-related phenotypes [4–6].
Adipose tissue, the storage site of triglycerides, acts as an endocrine organ contributing to regulate energy homeostasis [7, 8]. Anomalous fat accumulation in obesity increases the risk of severe diseases such as metabolic syndrome, type 2 diabetes, atherosclerosis, and cancer . Evidence suggests that in obesity, adipose tissues stay in a state of subclinical chronic inflammation  mainly prompted by the massive recruitment of macrophages into adipose tissues . Adipogenesis is tightly controlled by a mixture regulatory signals including endocellular transcription factors and circulating hormones . Additional regulators of adipogenesis also include microRNAs (miRNAs) [13–17]. These small non-coding RNAs, which have a length of only 20–24 nucleotides, are involved in the modulation of gene expression and, consequently, in the control of numerous cell functions [18–20]. At present, more than 2.000 different miRNAs have been described in humans and their number is still increasing (miRBase) (http://www.mirbase.org). miRNAs are important elements of the cell epigenetic machinery that, opposing to the other known players of epigenetic regulation (DNA methylation and post-translational histone modifications), can affect gene expression by binding to the 3′ untranslated sequence of a target messenger RNA (mRNA) . Each miRNA can target many transcripts, and transcription levels of a gene may be regulated by multiple miRNAs [21, 22]. Computational and experimental analyses indicate that endogenous miRNAs regulate the expression of up to 30 % of mouse and human genes . miRNAs have been shown to be involved in fundamental cellular processes including proliferation, differentiation, DNA repair, apoptosis, and metabolism [24, 25].
As a consequence, remarkable evidence indicates that dysregulation of miRNAs is causative and/or indicator of a number of disease processes including cancer . Altered circulating miRNA profiles have been linked to numerous cardiometabolic diseases, including hypertension (let-7e), type 2 diabetes (miR-126), hepatic injury (miR-122), and atherosclerosis (miR-223) .
The presence of extracellular circulating miRNAs in blood and other biological fluids is well established [27, 28], but the reason for the nuclease resistance of miRNAs outside the cell remained unclear for a long period . The discovery that cell-free miRNAs are detectable in serum and plasma and that their expression varies as a result of disease offers great potential for the discovery of novel health/disease biomarkers and candidate therapeutic targets .
Various miRNAs have been recognized to play a role in adipogenesis . Preliminary evidence showed that some circulating miRNAs are associated with obesity and related metabolic disturbances in adults as well as in children (see [16, 32] for review).
The primary aim of our study was to identify circulating miRNAs potentially associated to early obesity via an integrated study comprising miRNA signatures and bioinformatic analyses in children participating to the I.Family survey. Here, we present the results of a pilot study on a relatively small sample of children belonging to the Italian cohort of the I.Family survey.
Results and discussion
Distinct miRNA expression in plasma of OW/Ob and NW children
Anthropometric and metabolic variables of the individuals included in this study
10.50 ± 2.76
10.70 ± 1.57
16.45 ± 1.71
31.68 ± 4.32
53.80 ± 22.08
84.20 ± 30.01
146.80 ± 24.10
163.40 ± 33.87
62.30 ± 12.33
52.10 ± 12.37
77.60 ± 22.51
100.70 ± 30.57
91 ± 6.16
94.6 ± 3.53
The geometric mean Ct value of multiple selected housekeeping genes was 26.13 for NW and 26.38 for OW/Ob. Reverse transcription control average C t (RTC) was 20.68 for NW and 20.75 for OW/Ob. Absent calls were 1.82 % for NW and 2.34 % for OW/Ob.
Differential miRNA expression profile determined by gene-arrays screening
Differential miRNA profile confirmed by RT-qPCR
miRNAs target prediction, KEGG pathway analysis
The specialized algorithm miRPath v2.0  was used to understand the molecular functions of deregulated miRNAs and to predict their putative mRNA targets. To recognize biologically relevant molecular networks and pathways, we further explored miRNA targets by analyzing molecular interactions on the comprehensive knowledge base including the Kyoto Encyclopedia of Genes and Genomes (KEGG)  and Ingenuity Pathways Analysis. Both IPA and KEGG technologies provided hints for the functions of miRNAs.
Arrhythmogenic right ventricular cardiomyopathy (ARVC)
Notch signaling pathway
Regulation of actin cytoskeleton
Dorso-ventral axis formation
Bacterial invasion of epithelial cells
Endocrine and other factor-regulated calcium reabsorption
Glycosaminoglycan biosynthesis—heparan sulfate/heparin
Pathways in cancer
miRPath analysis of Hsa-mir-31-5p interactions
(A) KEGG-predicted pathway derived from DIANA-microT-CDS database
Metabolism of xenobiotics by cytochrome P450
Inositol phosphate metabolism
Fatty acid metabolism
(B) KEGG-predicted targets provided by the DIANA-TarBase v6.0 database
Regulation of actin cytoskeleton
Bacterial invasion of epithelial cells
Pathogenic Escherichia coli infection
Chemokine signaling pathway
Hypertrophic cardiomyopathy (HCM)
T cell receptor signaling pathway
Leukocyte transendothelial migration
Arrhythmogenic right ventricular cardiomyopathy (ARVC)
Dorso-ventral axis formation
Computational prediction of miR-31 identified the following targets in the “fatty acid degradation” pathway: (a) the enoyl-CoA hydratase, a key enzyme in the breakdown of fatty acids essential to metabolizing fatty acids to produce both acetyl CoA and energy, and (b) the alcohol dehydrogenase 1A (Additional file 1: Figure S1).
Among the top predicted targets of miR-31-5p, we also identify the CCAAT-enhancer-binding protein-α (C/EBPα), a main modulator of the expression of genes involved in cell cycle regulation as well as adipocyte functions during adipogenesis (Fig. 2).
As relevant miR-206 targets, we recognized the “insulin signaling,” the “pentose phosphate” pathways, and the liver X receptor-α (LXRα), a ligand-dependent transcription factor playing a relevant role in the metabolism and homeostasis of lipids, cholesterol, bile acids, and steroid hormones with a preeminent expression in the liver.
Further, we identified the oleoyl-ACP hydrolase, a medium-chain acyl-[acyl-carrier-protein] hydrolase in the “fatty acid biosynthesis” pathway as a target of miR-2355 (p = 0.005).
Finally, the target prediction of concurrently deregulated miRNAs gave evidence for the “Notch” signaling (p = 0.0008) (Additional file 2: Figure S2) and the “Arrhythmogenic right ventricular cardiomyopathy” pathways (p = 2.905898 × 10−05) (Additional file 3: Figure S3).
We identified three circulating miRNAs differentially expressed in plasma samples of OW/Ob children. Bioinformatic exploratory analysis predicted their involvement in critical pathways including lipid metabolism and adipocyte differentiation.
Adipose tissue is involved in pathophysiological processes because it influences metabolism. In adipogenesis, pluripotent adipose-derived mesenchymal stem cells differentiate into adipocytes. This process is finely regulated by both hormones and transcription factors in a complex signaling network. A number of miRNAs and their targets were recently identified to play regulatory roles in adipocyte differentiation [35–37]. To date, up to 221 miRNAs have been found to be expressed or dysregulated in adipose tissue or adipocytes in mammals .
miR-31 directly targets C/EBPα transcript whose levels are in vivo downregulated by miR-31 . Notably, miR-31 expression is inversely correlated with adipogenic differentiation status , and consequently, it has a decisive anti-adipogenic function.
We found that miR-31 was significantly upregulated in plasma of OW/Ob children identifying C/EBPα as a privileged target of miR-31. miR-31 upregulation may be considered as an apparent contradiction, but it probably reflects the well-established role of circulating miRNAs as “endocrine signals.” The concept that miRNAs may be transferred from donor to recipient cells is now well established with the relevant implication that miRNAs could affect gene expression distally, so acting as paracrine or endocrine signaling molecules. Accordingly, a panel of “endocrine miRNAs” released from fat cells that may be used as a marker of disturbed adipose tissue function has been characterized . The idea that miR-31 acts as an obesity controlling signal seems to be supported by the literature. Actually, in analyzing the expression of miRNAs in the livers of leptin-deficient obese (ob/ob) mice compared to normal C57BL/6 mice, Li et al. found that ob/ob mice exhibited an upregulation of miR-31 , despite its well-established anti-adipogenic role.
Computational prediction of miR-31 recognized several targets in the “fatty acid degradation” pathway (p = 0.027) (Additional file 1: Figure S1). Whether miR-31 deregulation plays a causative role or is merely an epiphenomenon in the etiology of obesity still remains to be determined.
In this study, we also observed that miR-206 was significantly downregulated in plasma of OW/Ob children. miR-206 belongs to the cluster of the so-called myomiRs and represents one of the most studied tissue-specific miRNAs involved in skeletal muscle differentiation. Several studies indicate that miR-206 plays a key role in the growth and development of skeletal muscle, promoting myogenic differentiation . Moreover, it has been well established that it is related to the pathogenesis of numerous diseases including heart failure, chronic obstructive pulmonary disease, and Alzheimer’s disease and in numerous kinds of cancers. Interestingly, in most of these diseases, miR-206 levels are downregulated, which leads to the recognition of miRNA as a “disease-avoiding” molecule. Moreover, miR-206 expression is abundant in brown adipocytes but missing in white adipocytes . Remarkably, it has been shown that miR-206 suppresses liver LXRα, a gene target of PPAR, consequently inhibiting lipogenesis and controlling lipid metabolism . The insulin signaling pathway (p = 0.01, 11 target genes) and the pentose phosphate pathway (p = 0.01, 4 target genes) were also identified as candidate targets of miR-206. Therefore, the reduced level of circulating miR-206 in OW/Ob children may be either representative of a functional decreased synthesis  or an epiphenomenon reflecting the reduced relative skeletal muscle mass in the OW/Ob group.
Among the top pathways that were predicted to be targeted by the three differentially expressed miRNAs, “Notch” signaling appears relevant (p = 0.0008, 5 target genes, Additional file 2: Figure S2). Notch signaling represents a basic biological pathway primarily related to cell communication and development with extensive integration and crosstalk with other signaling pathways. Perturbations in this pathway have been linked to many genetic disorders and cancers. Accordingly, the direct involvement of Notch signaling has recently been demonstrated in obesity and type 2 diabetes : through transcriptional binding partners FoxO1 and Rbpj, Notch activation enhances pathologic glucose output and fatty liver in obese, in combination with a reduced Akt activation by insulin. Activation of the Notch pathway in obesity results in heightened diabetes and fatty liver production, as a consequence of an augmented lipogenesis that is not blocked despite the reduced insulin signaling. Notably, blocking Notch signaling in mice caused white fat cells to transform into beige fat consequently reducing the risk of obesity and related health problems .
While we recognize the limitations and the relative reliability of our target prediction analyses, this hypothesis-generating exercise suggested the involvement of the detected miRNAs in metabolic relevant pathways, to be obviously confirmed by larger and ad hoc designed studies.
Our results confirm those of recent papers indicating that circulating miRNAs are deregulated in obese children  and are associated with obesity and lipid abnormalities in adolescents . However, differentially detected miRNA signals in refs. 46 and 47 are not consistent with those identified in our study. Of note, all these miRNAs are present in the 372 miRNA arrays we used in our study but did not pass the threshold we set to consider them as differentially expressed. A solid explanation of these discrepancies is not at hand. Considering the rapidly changing scenario in this field, the interpretation and discussion of the inconsistent findings regarding the association of circulating miRNAs with obesity in children [45, 46] is far beyond the objective of this study. However, recent reviews extensively summarized potentiality, but also uncertainties, of the role of extracellular miRNAs as potential regulators and/or biomarkers of obesity and metabolic diseases [16, 32]. In particular, the inconsistencies of the findings may be ascribed to inter-population differences, but even, and perhaps most likely, to methodological differences in sampling, RNA isolation, detection, and normalization . The need for consistent analysis strategies and quality control has been recently discussed .
This pilot study showed that three circulating miRNAs are differentially regulated in OW/Ob as compared with NW children. Although causal pathways cannot be firmly inferred—but only suggested—by these preliminary results, that deserve confirmation in other populations, it is conceivable that circulating miRNAs may be novel biomarkers of obesity and related metabolic disturbances.
The planned completion of the miRNA studies in the much larger sample of children from the eight European countries participating to the I.Family project will certainly contribute to the ongoing debate whether circulating miRNAs play a role in disease pathogenesis, are (early) indicators of a metabolic dysfunction, or both.
The I.Family project aimed to assess the determinants of eating behavior in children and adolescents of eight European countries and related health outcomes was built on the IDEFICS cohort (http://www.ideficsstudy.eu), established in 2006. We selected our sample among OW/Ob and NW participants of the Italian cohort of the project. The weight status of children was defined according to age- and sex-specific BMI categories . Ten OW/Ob (four females and six males, age 10.7 ± 1.5 years, BMI 31.7 ± 4.3 kg/m2) and ten NW (five females and five males, age 10.5 ± 2.7 years, BMI 16.4 ± 1.7 kg/m2) children were included in the study.
The study protocol was approved by the local Ethics Committee of the local Health Authority (ASL Avellino) and informed written parental consent was obtained for each participant.
Whole blood was collected according to the standard operating procedures in BD Vacutainer® blood collection tubes. The optimal approach for miRNA purification from whole blood necessitates collecting fresh blood and processes the sample as quickly as possible. Accordingly, plasma was isolated by centrifugation at 1900×g for 10 min at 4 °C in an Eppendorf benchtop centrifuge (Eppendorf, Germany), aliquoted into 1.5 mL Eppendorf tubes, and promptly stored at −80 °C, in the absence of freeze-thaw cycles, until processing.
miRNA extraction, reverse transcription, and pre-screening
Prior to miRNA extraction, spectrophotometry was carried out on plasma samples to test for hemolysis by measuring the absorbance of free hemoglobin at 414 nm; samples with OD414 greater than 0.2 were excluded from the study. RNA was extracted from both plasma pooled samples and from individual plasma samples using the commercial column-based system miRNeasy serum/plasma Kit (Qiagen, Germany) according to the manufacturer’s instructions with minor modifications. Plasma was thawed on ice and spun for 10 min at 13,000×g at 4 °C in an Eppendorf benchtop centrifuge (Eppendorf, Germany) to pellet any debris. An aliquot of 200 μL of plasma was transferred to a new microcentrifuge tube and 800 μL of a Qiazol (Qiagen) mixture, containing 1.25 μg/mL of MS2 bacteriophage RNA (Roche Applied Science), was added to the sample. Caenorhabditis elegans miR-39 (C.el-miR-39) was spiked into the sample (5.6 × 108 molecules) before the extraction process to assess RNA recovery. A rinse step (500 μL Qiagen RPE buffer) was repeated two times. Total RNA was eluted by adding 14 μL of RNase-free water to the membrane of the column and incubating for 1 min before centrifugation at 13,000×g for 1 min at room temperature. The obtained RNA was stored at −80 °C.
Complementary DNA (cDNA) was generated by means of the miScript RTII kit (Qiagen), in HiSpec Buffer with miRNA-specific stem-looped RT primers. Nine microliters of isolated RNA were reverse transcribed in 20-μL reactions according to the manufacturer’s recommendations. The obtained cDNA was diluted 1:5.5 in RNase-free water. This procedure was repeated three times. To identify hemolyzed samples, the levels of miR-451, highly abundant in RBCs, were preliminarily assessed in cDNA samples by qPCR. Hemolyzed samples were excluded from the analysis.
Circulating miRNA profiling
Selected samples from ten OW/Ob and ten NW children (Table 1) were pooled, and miRNAs were extracted as described above. miRNA profiling was performed by using the Serum and Plasma 384HC miScript miRNA PCR Arrays (SABiosciences, Qiagen) to assess the expression of 372 miRNAs typically detectable in serum and plasma. The SYBR green-based qPCR assay was performed with cDNA (dil. 1:225) in 10-uL reaction volume using a 7900HT fast Real-Time PCR System instrument (Applied Biosystems). Reaction conditions were as follows: 15 min at 95 °C and 40 cycles of 15 s at 94 °C, 30 s at 60 °C, and 30 s at 72 °C with a set of technical controls on each plate which included the spike-in control and small non-coding RNAs (SNORDs) as reference genes. References for data normalization included the spike-in C.el-miR-39-3p, SNORD61, SNORD68, SNORD72, SNORD95, and SNORD96A, which have been verified to hold relatively stable expression levels across tissues and cell types. Normalization of expression was done using the geometric mean of the endogenous controls. All assays were inspected for distinct melting curves, and the Tm was checked. The amplification efficiency was also calculated. Ct values >35 were considered as negative amplification. Data that did not pass these criteria were omitted from further analysis. The efficiency of reverse transcription was assessed with miRTC.
Statistically significant deregulated miRNAs with a fold-change ≥4.0 were selected and further validated by individual microRNA assays (SABiosciences, Qiagen) performed on the individual plasma samples of ten OW/Ob and ten NW children, using the miScript SYBR green PCR kit (Qiagen) according to the manufacturer’s instructions. RT-qPCR was performed in triplicates using the following conditions: 95 °C for 15 min, followed by 40 cycles of 94 °C for 15 s, 55 °C for 30 s, and 70 °C for 30 s. Relative miRNA levels were determined by the ΔΔCt method using the C.el-miR-39 as endogenous normalizer.
On pooled samples, RT-qPCR quantifications were performed by using SDS 2.3 (Life Technologies, USA) to generate Ct values.
Gene Arrays data were organized in Excel (Microsoft), and the integrated web-based software package for miScript Arrays analysis was exploited. This software performs quantification using the ΔΔCt method and interpretation of the control assays. Tools are available at the SABiosciences data analysis web portal.
Data on the individual samples are presented as the average of at least three independent experiments ± s.e.. Data Assist 3.1 software packages (Life Technologies, USA) were used to generate relative “fold-change” values. Student’s t test, performed with the GraphPad Prism 6 software, was used to determine the significance of any difference in the levels of miRNA expression between NW and OW/Ob subjects with alpha <0.05.
Obesity-associated miRNAs were tested using the well-established target prediction tool miRPath v2.0 . miRPath achieves advanced analysis pipelines, such as hierarchical clustering of miRNAs and pathways based on the levels of their interactions. miRNA targets (in CDS or 3′-UTR regions) were predicted by the DIANA-microT-CDS algorithm or even experimentally validated miRNA interactions derived from DIANA-TarBase v6.0. Predicted and/or validated interactions were subsequently combined with merging and meta-analysis algorithms.
All the predicted targets were further analyzed through the use of Kyoto Encyclopedia of Genes and Genomes (KEGG)  and QIAGEN’s Ingenuity® Pathway Analysis (IPA®, QIAGEN Redwood City, www.qiagen.com/ingenuity). IPA identified lists of genes that satisfy specific biological criteria scoring these genes towards pathways in the Ingenuity Knowledge Base which embraces a large database of biological and chemical relationships extracted from scientific literature.
The work was done as part of the IDEFICS study (http://www.idefics.eu) and the I.Family Study (http://www.ifamilystudy.eu/). We gratefully acknowledge the financial support of the European Community within the Sixth RTD Framework Programme Contract No. 016181 (FOOD) for the IDEFICS study and within the Seventh RTD Framework Programme Contract No. 266044 for the I.Family study.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Zimmet P, Alberti KG, Shaw J. Global and societal implications of the diabetes epidemic. Nature. 2001;414:782–7.View ArticlePubMedGoogle Scholar
- Kahn SE, Hull RL, Utzschneider KM. Mechanisms linking obesity to insulin resistance and type 2 diabetes. Nature. 2006;444:840–6.View ArticlePubMedGoogle Scholar
- Despres J-P, Lemieux I. Abdominal obesity and metabolic syndrome. Nature. 2006;444:881–7.View ArticlePubMedGoogle Scholar
- Berndt SI, Gustafsson S, Mägi R, Ganna A, Wheeler E, Feitosa MF, et al. Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture. Nat Genet. 2013;45:501–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Kunej T, Jevsinek Skok D, Zorc M, Ogrinc A, Michal JJ, Kovac M, et al. Obesity gene atlas in mammals. J Genomics. 2013;1:45–55.View ArticlePubMedPubMed CentralGoogle Scholar
- Loos RJF, Yeo GSH. The bigger picture of FTO: the first GWAS-identified obesity gene. Nat Rev Endocrinol. 2014;10:51–61.View ArticlePubMedPubMed CentralGoogle Scholar
- Rosen ED, Spiegelman BM. Adipocytes as regulators of energy balance and glucose homeostasis. Nature. 2006;444:847–53.View ArticlePubMedPubMed CentralGoogle Scholar
- Schaffler A, Muller-Ladner U, Scholmerich J, Buchler C. Role of adipose tissue as an inflammatory organ in human diseases. Endocr Rev. 2006;27:449–67.View ArticlePubMedGoogle Scholar
- Kopelman PG. Obesity as a medical problem. Nature. 2000;404:635–43.PubMedGoogle Scholar
- Kammoun HL, Kraakman MJ, Febbraio MA. Adipose tissue inflammation in glucose metabolism. Rev Endocr Metab Disord. 2014;15:31–44.View ArticlePubMedGoogle Scholar
- Kraakman MJ, Murphy AJ, Jandeleit-Dahm K, Kammoun HL. Macrophage polarization in obesity and type 2 diabetes: weighing down our understanding of macrophage function? Front Immunol. 2014;5:470.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang QQ, Lane MD. Adipogenesis: from stem cell to adipocyte. Annu Rev Biochem. 2012;81:715–36.View ArticlePubMedGoogle Scholar
- Xie H, Lim B, Lodish HF. MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity. Diabetes. 2009;58:1050–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Rottiers V, Naar AM. MicroRNAs in metabolism and metabolic disorders. Nat Rev Mol Cell Biol. 2012;13:239–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim HJ, Cho H, Alexander R, Patterson HC, Gu M, Lo KA, et al. MicroRNAs are required for the feature maintenance and differentiation of brown adipocytes. Diabetes. 2014;63:4045–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Arner P, Kulyte A. MicroRNA regulatory networks in human adipose tissue and obesity. Nat Rev Endocrinol. 2015;11:276–88.View ArticlePubMedGoogle Scholar
- Li J, Zhou C, Li J, Su Z, Sang H, Jia E, et al. Global correlation analysis for microRNA and gene expression profiles in human obesity. Pathol Res Pract. 2015;211:361–8.View ArticlePubMedGoogle Scholar
- Bartel DP. MicroRNAs: target recognition and regulatory functions. Cell. 2009;136:215–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghildiyal M, Zamore PD. Small silencing RNAs: an expanding universe. Nat Rev Genet. 2009;10:94–108.View ArticlePubMedPubMed CentralGoogle Scholar
- Ha M, Kim VN. Regulation of microRNA biogenesis. Nat Rev Mol Cell Biol. 2014;15:509–24.View ArticlePubMedGoogle Scholar
- Ebert MS, Sharp PA. Roles for microRNAs in conferring robustness to biological processes. Cell. 2012;149:515–24.View ArticlePubMedPubMed CentralGoogle Scholar
- Krek A, Grün D, Poy MN, Wolf R, Rosenberg L, Epstein EJ, et al. Combinatorial microRNA target predictions. Nat Genet. 2005;37:495–500.View ArticlePubMedGoogle Scholar
- Lewis BP, Burge CB, Bartel DP. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell. 2005;120:15–20.View ArticlePubMedGoogle Scholar
- Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.View ArticlePubMedGoogle Scholar
- van Rooij E. The art of microRNA research. Circ Res. 2011;108:219–34.View ArticlePubMedGoogle Scholar
- Mendell JT, Olson EN. MicroRNAs in stress signaling and human disease. Cell. 2012;148:1172–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Park NJ, Zhou H, Elashoff D, Henson BS, Kastratovic DA, Abemayor E, et al. Salivary microRNA: discovery, characterization, and clinical utility for oral cancer detection. Clin Cancer Res. 2009;15:5473–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Kosaka N, Iguchi H, Ochiya T. Circulating microRNA in body fluid: a new potential biomarker for cancer diagnosis and prognosis. Cancer Sci. 2010;101:2087–92.View ArticlePubMedGoogle Scholar
- Valadi H, Ekstrom K, Bossios A, Sjostrand M, Lee JJ, Lotvall JO. Exosome-mediated transfer of mRNAs and microRNAs is a novel mechanism of genetic exchange between cells. Nat Cell Biol. 2007;9:654–9.View ArticlePubMedGoogle Scholar
- Turchinovich A, Samatov TR, Tonevitsky AG, Burwinkel B. Circulating miRNAs: cell-cell communication function? Front Genet. 2013;4:119.View ArticlePubMedPubMed CentralGoogle Scholar
- McGregor RA, Choi MS. microRNAs in the regulation of adipogenesis and obesity. Curr Mol Med. 2011;11:304–16.View ArticlePubMedPubMed CentralGoogle Scholar
- Deiuliis JA. MicroRNAs as regulators of metabolic disease: pathophysiologic significance and emerging role as biomarkers and therapeutics. Int J Obes (Lond). 2016;40:88–101.View ArticleGoogle Scholar
- Vlachos IS, Kostoulas N, Vergoulis T, Georgakilas G, Reczko M, Maragkakis M, et al. DIANA miRPath v.2.0: investigating the combinatorial effect of microRNAs in pathways. Nucleic Acids Res. 2012;40(Web Server issue):W498–504.View ArticlePubMedPubMed CentralGoogle Scholar
- Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42:D199–205.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen L, Song J, Cui J, Hou J, Zheng X, Li C, et al. microRNAs regulate adipocyte differentiation. Cell Biol Intern. 2013;37:533–46.View ArticleGoogle Scholar
- Zhang R, Wang D, Xia Z, Chen C, Cheng P, Xie H, et al. The role of microRNAs in adipocyte differentiation. Front Med. 2013;7:223–30.View ArticlePubMedGoogle Scholar
- Son YH, Ka S, Kim AY, Kim JB. Regulation of adipocyte differentiation via microRNAs. Endocrinol Metab (Seoul). 2014;29:122–35.View ArticleGoogle Scholar
- Linhart HG, Ishimura-Oka K, DeMayo F, Kibe T, Repka D, Poindexter B, et al. C/EBPalpha is required for differentiation of white, but not brown, adipose tissue. Proc Natl Acad Sci U S A. 2001;98:12532–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang YF, Zhang Y, Li XY, Li C, Tian W, Liu L. Expression of miR-31, miR-125b-5p, and miR-326 in the adipogenic differentiation process of adipose-derived stem cells. OMICS. 2009;13:331–6.View ArticlePubMedGoogle Scholar
- Li S, Chen X, Zhang H, Liang X, Xiang Y, Yu C, et al. Differential expression of microRNAs in mouse liver under aberrant energy metabolic status. J Lipid Res. 2009;50:1756–65.View ArticlePubMedPubMed CentralGoogle Scholar
- Novak J, Kruzliak P, Bienertova-Vasku J, Slaby O, Novak M. MicroRNA-206: a promising theranostic marker. Theranostics. 2014;4:119–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Walden TB, Timmons JA, Keller P, Nedergaard J, Cannon B. Distinct expression of muscle-specific microRNAs (myomirs) in brown adipocytes. J Cell Physiol. 2009;218:444–9.View ArticlePubMedGoogle Scholar
- Zhong D, Huang G, Zhang Y, Zeng Y, Xu Z, Zhao Y, et al. MicroRNA-1 and microRNA-206 suppress LXRα-induced lipogenesis in hepatocytes. Cell Signal. 2013;25:1429–37.View ArticlePubMedGoogle Scholar
- Bi P, Shan T, Liu W, Yue F, Yang X, Liang XR, et al. Inhibition of Notch signaling promotes browning of white adipose tissue and ameliorates obesity. Nat Med. 2014;20:911–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Prats-Puig A, Ortega FJ, Mercader JM, Moreno-Navarrete JM, Moreno M, Bonet N, et al. Changes in circulating microRNAs are associated with childhood obesity. J Clin Endocrinol Metab. 2013;98:E1655–60.View ArticlePubMedGoogle Scholar
- Can U, Buyukinan M, Yerlikaya FH. The investigation of circulating microRNAs associated with lipid metabolism in childhood obesity. Pediatr Obes. 2015. doi:10.1111/ijpo.12050.
- Mestdagh P, Hartmann N, Baeriswyl L, Andreasen D, Bernard N, Chen C, et al. Evaluation of quantitative miRNA expression platforms in the microRNA quality control (miRQC) study. Nat Methods. 2014;11:809–15.View ArticlePubMedGoogle Scholar
- Cole TJ, Lobstein T. Extended international (IOTF) body mass index cut-offs for thinness, overweight and obesity. Pediatr Obes. 2012;7:284–94.View ArticlePubMedGoogle Scholar