Next Article in Journal
Single Nucleotide Polymorphisms in the Leptin-a Gene and Associations with Growth Traits in the Orange-Spotted Grouper (Epinephelus coioides)
Next Article in Special Issue
The Role of microRNA in Gastric Malignancy
Previous Article in Journal
A One-Step Homogeneous Sandwich Immunosensor for Salmonella Detection Based on Magnetic Nanoparticles (MNPs) and Quantum Dots (QDs)
Previous Article in Special Issue
Advances in the Techniques for the Prediction of microRNA Targets
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 14
Issue 4
10.3390/ijms14048611
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

MicroRNA Transcriptomes Relate Intermuscular Adipose Tissue to Metabolic Risk

by
Jideng Ma
1,†,
Shuzhen Yu
1,†,
Fengjiao Wang
1,
Lin Bai
1,
Jian Xiao
1,
Yanzhi Jiang
2,
Lei Chen
3,
Jinyong Wang
3,
Anan Jiang
1,
Mingzhou Li
1,* and
Xuewei Li
1,*
1
Institute of Animal Genetics & Breeding, College of Animal Science & Technology, Sichuan Agricultural University, Ya'an 625014, China
2
College of Life and Basic Sciences, Sichuan Agricultural University, Ya'an 625014, China
3
Chongqing Academy of Animal Science, Chongqing 402460, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this work.
Int. J. Mol. Sci. 2013, 14(4), 8611-8624; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms14048611
Submission received: 20 March 2013 / Revised: 15 April 2013 / Accepted: 17 April 2013 / Published: 22 April 2013
(This article belongs to the Special Issue Regulation by non-coding RNAs 2013)
Browse Figures
Versions Notes

Abstract

:
Intermuscular adipose tissue is located between the muscle fiber bundles in skeletal muscles, and has similar metabolic features to visceral adipose tissue, which has been found to be related to a number of obesity-related diseases. Although various miRNAs are known to play crucial roles in adipose deposition and adipogenesis, the microRNA transcriptome of intermuscular adipose tissue has not, until now, been studied. Here, we sequenced the miRNA transcriptomes of porcine intermuscular adipose tissue by small RNA-sequencing and compared it to a representative subcutaneous adipose tissue. We found that the inflammation- and diabetes-related miRNAs were significantly enriched in the intermuscular rather than in the subcutaneous adipose tissue. A functional enrichment analysis of the genes predicted to be targeted by the enriched miRNAs also indicated that intermuscular adipose tissue was associated mainly with immune and inflammation responses. Our results suggest that the intermuscular adipose tissue should be recognized as a potential metabolic risk factor of obesity.

1. Introduction

Adipose tissues (ATs) play a vital role in energy homeostasis and process the largest energy reserve in the body of animals. The rapidly expanding adipokine family is secreted by ATs [1], and, as a result, AT has been identified as an endocrine organ that influences a variety of physiological and pathological processes (such as immunity and inflammation) [2,3] that are involved in the development of metabolic diseases such as cardiovascular disease and type 2 diabetes mellitus [4–6]. Functional and metabolic differences between the visceral and subcutaneous ATs have been well documented. Subcutaneous AT mainly affects metabolic processes, while visceral AT has been identified as a metabolic risk factor for obesity. Recent studies have revealed that the intermuscular adipose tissue (IMAT), which is located between the muscle fiber bundles in skeletal muscles, has similar functional and metabolic features as the visceral ATs [7,8]. Indeed, IMAT was found in greater amounts than visceral AT in acromegaly patients despite their increased muscle mass, suggesting that increased amounts of AT in muscles might be associated with growth hormone-induced insulin resistance [9].
MicroRNAs (miRNAs) are endogenous small non-coding RNAs that modulate gene expression at a post-transcriptional level by binding to the 3′ untranslated region (3′-UTR) of the target mRNAs [10]. During the past decade, various miRNAs that play crucial roles in adipose deposition and adipogenesis have been identified. Typically, miR-143 was identified as a pro-adipogenic modulator during pre-adipocyte differentiation [11,12]. MiR-103 [13] and the miR-17-92 cluster [14] were reported to accelerate adipocyte differentiation. MiR-27a [15], miR-27b [16], miR-448 [17] and miR-15a [18] were demonstrated to suppress adipogenic differentiation. MiR-519d [19], miR-335 [20] and miR-377 [21] were associated with lipid metabolism disorders. However, features of the miRNA transcriptome of IMAT have yet to be investigated.
Sus scrofa (pig) is emerging as an ideal biomedical model for obesity and metabolic disorders in human because of the similarity in metabolic features and proportional organ sizes in these two species [22]. To decipher the unique metabolic and functional features of IMAT, we sequenced the miRNA transcriptomes of porcine IMAT by small RNA-sequencing and compared it with a representative subcutaneous adipose tissue, superficial abdominal subcutaneous adipose tissue (sASAT). We identified various known, conserved, and putative novel porcine miRNAs in these two tissues. Notably, the sASAT-enriched miRNAs were related mainly to lipid metabolic homeostasis, while the IMAT-enriched miRNAs were related mainly to inflammation and diabetes, and the target genes of the IMAT-enriched miRNAs were primarily associated with inflammatory and diabetes processes. Together, these findings indicated the metabolic risk of IMAT. Our results will contribute to studies into the role of IMAT in obesity-related metabolic disease.

2. Results and Discussion

2.1. Transcriptome Sequencing Data

We used a small RNA-sequencing approach to sequence the miRNA transcriptomes of porcine IMAT and sASAT and obtained 17.76 million (M) and 18.50 M raw reads, respectively. More than 80% of the raw reads passed the quality filters (see Methods) and were termed the high-quality reads (IMAT: 14.93 M, 84.09%; sASAT: 14.92 M, 81.59%) (Supplementary Table S1). The high-quality reads in both transcriptomes exhibited the canonical size range distribution that is common to mammalian miRNAs (Figure 1a). The vast majority of the reads were 21–23 nucleotides (nt) in length. The 22-nt reads accounted for 61.50% of all the high-quality reads, followed by the 21-nt (14.00%) and 23-nt (13.99%) reads. This result indicated the reliability of using the small RNA-sequencing approach to generate miRNA reads as candidates for further analysis.

2.2. MiRNA Profiling of sASAT and IMAT

A total of 597 mature miRNAs corresponding to 453 miRNA precursors (pre-miRNAs) were identified in the two libraries by mapping them to the pig genome. In agreement with previous reports [23,24], we found that all the miRNA classes consisted of multiple mature variants (the isomiRs). The most abundant isomiR in each class was picked as the reference sequence for that class [23] based on the evidence that there was a significant positive correlation between the counts of most abundant isomiR and the total counts of all isomiRs in the same class (IMAT: Spearman’s r = 0.98, p < 10−5; sASAT: Spearman’s r = 0.97, p < 10−4).
The identified mature miRNAs and their precursors were divided into three subgroups according to alignment criteria (Supplementary Table S2) as: (1) Porcine known miRNAs: 297 miRNAs mapped to 176 known porcine pre-miRNAs; specifically, 210 were in miRBase 18.0 [25] and 87 were novel miRNA*s; (2) Porcine conserved miRNAs: 107 miRNAs mapped to 71 other known mammalian pre-miRNAs in miRBase 18.0 and these pre-miRNAs mapped to the pig genome. These miRNAs were labeled with the names of the corresponding conserved miRNAs; (3) Porcine putative new miRNAs: 230 miRNAs (longer than 18 nt and unmapped to any known mammalian pre-miRNAs in miRBase 18.0) encompassing 206 candidate pre-miRNAs that were predicted RNA hairpins derived from the pig genome, and were labeled PPN (Porcine putative new). Notably, there are the distinct pre-miRNAs coding the identical mature miRNAs, which resulting in 617 miRNAs (i.e., reference sequence) corresponding to 597 unique miRNA sequences (Supplementary Table S3).
The identified miRNAs exhibited a large dynamic range of read counts ranging from 3 to millions. The vast majority of miRNAs (IMAT: 61.11%; sASAT: 61.03%) were in low abundance (3 to 100 read counts) and belonged mainly to the porcine conserved and putative new miRNA groups. Only a few miRNAs (IMAT: 2.11%; sASAT: 1.17%) were in high abundance (>100,000 read counts) and they belonged mainly to the porcine known miRNA group (Figure 1b,c). This result suggests that the low-abundance conserved and putative new miRNAs may have escaped from previous detection efforts.
We found that the top ten miRNAs with the highest abundance contributed 67.46% and 82.94% of the total counts in the IMAT and sASAT libraries, respectively, and eight miRNAs were shared by two libraries in the top 10 positions (Figure 1d). The high abundance of these miRNAs implies that they may have housekeeping cellular roles and may be the main regulatory miRNAs in adipogenesis [11,26,27] and cellular basal metabolism [28,29]. For example, let-7a-5p [12], miR-148a-3p [26], miR-21-5p [27], miR-143-3p [30] and miR-101-3p [13] have been reported to be up-regulated during 3T3-L1 pre-adipocyte differentiation, whereas miR-27b-3p was found to be down-regulated during adipogenesis of human multipotent adipose-derived stem cells [31] and miR-199a-5p was up-regulated in subcutaneous AT in obese versus non-obese individuals [13].

2.3. Inflammation- and Diabetes-Related miRNAs Enriched in IMAT

More than half of the unique miRNAs (351 of 597, 58.79%) were co-expressed in IMAT and sASAT. Only 171 (28.64%) and 75 (12.56%) of the unique miRNAs were expressed specifically in IMAT and sASAT, respectively (Figure 2a and Supplementary Table S3). It was well-known that miRNAs function in a dose-dependent manner [32], thus the less abundant miRNAs (<1000 read counts in both libraries) were considered to be less important and were filtered out. Of the 110 more abundant unique miRNAs (>1000 read counts in either library), 53 (48.18%) were determined to be differentially expressed (DE) between IMAT and sASAT using the IDEG6 program [33] (Figure 2a and Supplementary Table S4). The changes in expression patterns of the top 14 DE miRNAs with the highest read counts showed significant positive correlations between the q-PCR results and the small RNA-sequencing data (Person’s r = 0.894, p < 10−4), again highlighting the reliability of the small RNA-sequencing approach (Figure 2b). Moreover, in the process of q-PCR validation, we also found that all expression levels of selected miRNAs obtained by q-PCR within the biological replicates were highly correlated and with very low deviation, which not only indicated the high repeatability and reliability of the q-PCR approach but also reflected the high purity of our experimental samples (Supplementary Table S5).
Notably, many of the DE miRNAs (19 out of 53, 35.84%) were associated with inflammation and diabetes based on the annotations assigned using the Pathway Central database (SA Biosciences, Frederick, MD, USA) (Figure 2c). Eight inflammation-related and 9 diabetes-related miRNAs were present in higher abundance in the IMAT transcriptome compared with the sASAT transcriptome (Figure 2c). MiR-21 was found to be over-expressed at the inflammation site [34] and it has been suggested that miR-21 could act as a biomarker for inflammation in the aging process and cardiovascular disease [35]. MiR-101 was reported to be related to inflammation and chondrocyte extracellular matrix degradation [36]. Circulating miR-30a was up-regulated in diabetes patients and has been associated with insulin resistance [20]. The ectopic high expression of miR-103 could induce impaired glucose homeostasis or, conversely, the silencing of miR-103 could improve glucose homeostasis and insulin sensitivity [37].
In addition, various IMAT-enriched miRNAs related to pathological responses were found (Supplementary Table S4). For example, miR-182 and miR-183 (members of the miR-183-96-182 cluster) are well-characterized oncomiRs that can promote the clonal expansion of activated helper T lymphocytes [38,39]. MiR-200b, miR-200c and miR-141 (members of the miR-200 family) were reported to be significantly altered in bladder [40] and breast cancers [41]. MiR-191 was suggested as a biomarker for the diagnosis and prognosis of acute myeloid leukemia [42]. These results suggest that IMAT is associated mainly with inflammation- and diabetes-related responses, and should be deemed as a potential metabolic risk factor of obesity.
In contrast, the sASAT-enriched miRNAs (Supplementary Table S4) were mainly related to adipogenesis and lipid metabolism. For example, miR-378 was up-regulated in adipogenesis of human AT-derived stromal cells [43] and the over-expression of miR-378 in ST2 cell line was reported to promote lipid accumulation by enhancing de novo lipogenesis [44]. MiR-365 was revealed as a central regulator of brown fat differentiation and adipogenesis [45]. MiR-146a regulated mainly lipid accumulation induced by oxidized low-density lipoprotein [46].

2.4. Functional Enrichment Analyses of miRNA Target Genes

To further highlight the distinct functional features of IMAT and sASAT, the target genes of the top eight DE miRNAs enriched in sASAT (967 mRNA genes) and IMAT (1707 mRNA genes) were predicted using PicTar [47], TargetScan human 6.2 [48] and MicroCosm Targets (version 5.0) [49] (Supplementary Table S6), and analyzed using DAVID [50] to determine whether or not they were enriched for specific functional categories and pathways. Similar to the finding for the DE miRNAs, the target genes of the IMAT-enriched miRNAs were primarily associated with inflammatory and diabetes-related processes, such as “inflammatory response” (77 genes, p = 1.21 × 10−15), “cellular response to insulin stimulus” (20 genes, p = 4.14 × 10−6), “lymphocyte differentiation” (19 genes, p = 3.75 × 10−3), “regulation of interleukin-6 production” (9 genes, p = 1.17 × 10−2), “macrophage activation during immune response” (4 genes, p = 2.74 × 10−2), “chemokine and toll-like signaling pathways” (42 genes, p = 1.61 × 10−7) and “insulin signaling pathway” (22 genes, p = 1.51 × 10−2). In contrast, the target genes of the sASAT-enriched miRNAs were mainly associated with lipid and energy metabolism, such as “glycerophospholipid metabolic process” (46 genes, p = 5.10 × 10−18), “lipid biosynthetic process” (38 genes, p = 7.36 × 10−6), “glucose metabolic process” (16 genes, p = 1.53 × 10−2) and “Wnt signaling pathway” (20 genes, p = 5.27 × 10−4) (Figure 3). These results further suggested that while sASAT is mainly involved in metabolic homeostasis, IMAT is susceptive to inflammation and should be regarded as a potential metabolic risk factor.

3. Experimental Section

3.1. Animals and Sample Collection

Three 210-day-old female Landrace pigs with normal weight (111.67 ± 1.15 kg) were used. The piglets were weaned simultaneously at 28 ± 1 day of age. A starter diet provided 3.40 Mcal·kg−1 metabolisable energy (ME), 20.00% crude protein and 1.15% lysine from the thirtieth to sixtieth day after weaning. From the 61st to the 120th day, the diet contained 3.40 Mcal·kg−1 ME, 17.90% crude protein and 0.83% lysine. From the 121st to 210th day, the diet contained 3.40 Mcal kg−1 ME, 15.00% crude protein and 1.15% lysine. The animals were allowed access to feed and water ad libitum and lived under the same normal conditions.
The macroscopic IMAT were directly separated from the regions that were beneath the biceps femoris muscle fascia of porcine hind leg. Since IMAT preparation could be easily contaminated by its surrounding tissues, we paid maximum attention to eliminate the others especially, such as connective tissue and muscle tissue, and all samples were resected from central part of tissue block. The sASAT were from the subcutaneous tissue of central abdomen near the last rib. All samples were immediately frozen in liquid nitrogen and stored at −80 °C before total RNA extraction.

3.2. Small RNA Libraries Construction and High-throughput Sequencing

Total RNA was extracted using the mirVana™ miRNA isolation kit (Ambion, Austin, USA) following the manufacturer’s protocol. The integrity of total RNA was also tested via analysis by Bioanalyzer 2100 and RNA 6000 Nano LabChip Kit (Agilent, Palo Alto, CA, USA) with RIN number >6.0.
For a certain adipose tissue, equal amounts (5 μg) of total RNA isolated from three pigs were mixed. Approximately 15 μg of small RNA-enriched total RNA was prepared for Illumina sequencing. In general, the processing by Illumina consisted of the following successive steps: the small RNA ranged from 14 to 40 nt were purified by polyacrylamide gel electrophoresis (PAGE) and ligated specific adapters followed by polyacrylamide gel purification. Then the modified small RNA was reverse transcripted and amplified by RT-PCR. Finally, the enriched cDNA was sequenced on Genome Analyzer Instrument (GAI, Illumina, San Diego, CA, USA). The small RNA-sequencing data discussed in this publication have been deposited in NCBI’s Gene Expression Omnibus and are accessible through GEO Series accession number GSE30334.

3.3. Analysis of Small RNA-Sequencing Data

The raw reads were processed using Illumina’s Genome Analyzer Pipeline software and subsequently handled as described by Li et al. with some improvement [23]. After trimming off the sequencing adapters, the resulting reads was successively filtered by read length (only the read with the of 14 to 27 nt were retained), sequence component (containing <80% A, C, G or T; containing no more than two N (undetermined bases)) and copy numbers (the low-abundance reads (only the read with >3 counts were retained). Then the retained reads were searched against the NCBI [51], Rfam [52] and Repbase database [53] to remove porcine known classes of RNAs (i.e., mRNA, rRNA, tRNA, snRNA, snoRNA and repeats). The sequencing reads survived from above strict filter rules were deemed as “high-quality reads”.
The high-quality reads were mapped to the pig genome (Sscrofa9) using NCBI Local BLAST following five steps in order: (1) map the high-quality reads to the 228 known porcine pre-miRNAs (encoding 257 miRNAs) and then to 6716 known pre-miRNAs (encoding 7952 miRNAs) from 24 other mammals in miRBase 18.0; (2) map the mapped high-quality reads to pig genome to obtain their genomic locations and annotations in Ensembl release 59 (Sscrofa 9, April 2009); (3) cluster the unmapped sequences in step 1 that mapped to the pig genome as putative novel miRNAs; and (4) predict hairpin RNA structures of the high-quality reads in step 3 from the adjacent 60 nt sequences in either direction from the pig genome using UNAFold [54]. To avoid ambiguous reads that have been assigned to multiple positions in pig genome, only reads longer than 18 nt in length were included in step 4.

3.4. miRNA Differential Expression Analysis

Program IDEG6 [33] was employed for detecting the DE miRNAs between two libraries. A unique miRNA is considered to be differentially expressed when it simultaneously obtains p < 0.001 under three statistical tests (a Audic-Claverie test, a Fisher exact test and a Chi-squared 2 × 2 test) with the Bonferroni correction.

3.5. Prediction and Functional Annotation of miRNA Target Genes

The potential targets of a certain miRNA were predicted by PicTar [47], TargetScan human 6.2 [48] and, MicroCosm Targets Version 5.0 [49], and the pairwise overlaps of results from three programs composed the final predicted targets. The predictions were according to the interactions of human mRNA-miRNA due to the absence of porcine miRNAs in current version of above-mentioned algorithm. The gene ontology biological process (GO-BP) terms and KEGG pathway terms enriched in predicted target genes were determined using a DAVID bioinformatics resources [50].

3.6. Q-PCR Validation

The expression changes of 14 selected miRNAs were validated by an EvaGreen-based High-Specificity miRNA qRT-PCR Detection Kit (Stratagene, La Jolla, CA, USA) on the CFX96™ Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). The q-PCR validation were carried out on three biological replicates. The primer pairs were available in Supplementary Table S7. Three endogenous control genes (U6 snRNA, 18S rRNA and Met-tRNA) [23] were used in this assay. The ΔΔCt method was used to determine the expression level differences between surveyed samples. Normalized factors (NF) of three endogenous control genes and relative quantities of objective miRNAs were analyzed using the qBase software [55].

4. Conclusions

We have generated reliable miRNA transcriptomes of porcine sASAT and IMAT, and identified many known and novel miRNAs using small RNA-sequencing approach. We found that inflammationand diabetes-related miRNAs were enriched in IMAT compared with sASAT, which indicated the metabolic risk of the IMAT. A functional enrichment analysis of genes targeted by the enriched miRNAs also indicated that IMAT was mainly associated with the immune and inflammation response and may be a potential metabolic risk factor of obesity. The current study provides data that can be used in future studies to investigate the metabolic role of IMAT in obesity-related metabolic dysfunction. Our findings will also help promote the further development of the pig model for human metabolic research. It is also worth noting that further detailed comparision of IMAT between obese and non-obese individuals will be necessary and beneficial to decipher the role of miRNAs in adipogenesis and IMAT-related metabolic diseases.

Acknowledgments

This work was supported by grants from the National High Technology Research and Development Program of China (863 Program) (2013AA102502), the Specialized Research Fund of Ministry of Agriculture of China (NYCYTX-009), the Project of Provincial Twelfth Five Years’ Animal Breeding of Sichuan Province (2011YZGG15), and the National Special Foundation for Transgenic Species of China (2011ZX08006-003) to X.L. and M.L., the Chongqing Fund for Distinguished Young Scientists (CSTC2010BA1007).

Conflict of Interest

The authors declare no conflict of interest.

References

  1. Zhang, Y.; Proenca, R.; Maffei, M.; Barone, M.; Leopold, L.; Friedman, J.M. Positional cloning of the mouse obese gene and its human homologue. Nature 1994, 372, 425–432. [Google Scholar]
  2. Kadowaki, T.; Yamauchi, T.; Kubota, N.; Hara, K.; Ueki, K.; Tobe, K. Adiponectin and adiponectin receptors in insulin resistance, diabetes, and the metabolic syndrome. J. Clin. Invest 2006, 116, 1784–1792. [Google Scholar]
  3. Lago, F.; Dieguez, C.; Gómez-Reino, J.; Gualillo, O. The emerging role of adipokines as mediators of inflammation and immune responses. Cytokine Growth Factor Rev 2007, 18, 313–325. [Google Scholar]
  4. Arner, P. Insulin resistance in type 2 diabetes-role of the adipokines. Curr. Mol. Med 2005, 5, 333–339. [Google Scholar]
  5. Dogru, T.; Sonmez, A.; Tasci, I.; Bozoglu, E.; Yilmaz, M.I.; Genc, H.; Erdem, G.; Gok, M.; Bingol, N.; Kilic, S. Plasma visfatin levels in patients with newly diagnosed and untreated type 2 diabetes mellitus and impaired glucose tolerance. Diabetes Res. Clin. Pract. Suppl 2007, 76, 24–29. [Google Scholar]
  6. Berg, A.H.; Scherer, P.E. Adipose tissue, inflammation, and cardiovascular disease. Circ. Res 2005, 96, 939–949. [Google Scholar]
  7. Gallagher, D.; Kuznia, P.; Heshka, S.; Albu, J.; Heymsfield, S.B.; Goodpaster, B.; Visser, M.; Harris, T.B. Adipose tissue in muscle: A novel depot similar in size to visceral adipose tissue. Am. J. Clin. Nutr 2005, 81, 903–910. [Google Scholar]
  8. Boettcher, M.; Machann, J.; Stefan, N.; Thamer, C.; Häring, H.U.; Claussen, C.D.; Fritsche, A.; Schick, F. Intermuscular adipose tissue (IMAT): Association with other adipose tissue compartments and insulin sensitivity. J. Magn. Reson. Imaging 2009, 29, 1340–1345. [Google Scholar]
  9. Freda, P.U.; Shen, W.; Heymsfield, S.B.; Reyes-Vidal, C.M.; Geer, E.B.; Bruce, J.N.; Gallagher, D. Lower visceral and subcutaneous but higher intermuscular adipose tissue depots in patients with growth hormone and insulin-like growth factor I excess due to acromegaly. J. Clin. Endocrinol. Metab 2008, 93, 2334–2343. [Google Scholar]
  10. Nelson, P.; Kiriakidou, M.; Sharma, A.; Maniataki, E.; Mourelatos, Z. The microRNA world: Small is mighty. Trends Biochem. Sci 2003, 28, 534–540. [Google Scholar]
  11. Xie, H.; Lim, B.; Lodish, H.F. MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity. Diabetes 2009, 58, 1050–1057. [Google Scholar]
  12. Kajimoto, K.; Naraba, H.; Iwai, N. MicroRNA and 3T3-L1 pre-adipocyte differentiation. RNA 2006, 12, 1626–1632. [Google Scholar]
  13. Ortega, F.J.; Moreno-Navarrete, J.M.; Pardo, G.; Sabater, M.; Hummel, M.; Ferrer, A.; Rodriguez-Hermosa, J.I.; Ruiz, B.; Ricart, W.; Peral, B. MiRNA expression profile of human subcutaneous adipose and during adipocyte differentiation. PLoS One 2010, 5, e9022. [Google Scholar]
  14. Wang, Q.; Li, Y.C.; Wang, J.; Kong, J.; Qi, Y.; Quigg, R.J.; Li, X. miR-17-92 cluster accelerates adipocyte differentiation by negatively regulating tumor-suppressor Rb2/p130. Proc. Natl. Acad. Sci. USA 2008, 105, 2889–2894. [Google Scholar]
  15. Lin, Q.; Gao, Z.; Alarcon, R.M.; Ye, J.; Yun, Z. A role of miR-27 in the regulation of adipogenesis. FEBS J 2009, 276, 2348–2358. [Google Scholar]
  16. Kim, S.Y.; Kim, A.Y.; Lee, H.W.; Son, Y.H.; Lee, G.Y.; Lee, J.-W.; Lee, Y.S.; Kim, J.B. miR-27a is a negative regulator of adipocyte differentiation via suppressing PPARγ expression. Biochem. Biophys. Res. Commun 2010, 392, 323–328. [Google Scholar]
  17. Kinoshita, M.; Ono, K.; Horie, T.; Nagao, K.; Nishi, H.; Kuwabara, Y.; Takanabe-Mori, R.; Hasegawa, K.; Kita, T.; Kimura, T. Regulation of adipocyte differentiation by activation of serotonin (5-HT) receptors 5-HT2AR and 5-HT2CR and involvement of microRNA-448-mediated repression of KLF5. Mol. Endocrinol 2010, 24, 1978–1987. [Google Scholar]
  18. Andersen, D.C.; Jensen, C.H.; Schneider, M.; Nossent, A.Y.; Eskildsen, T.; Hansen, J.L.; Teisner, B.; Sheikh, S.P. MicroRNA-15a fine-tunes the level of Delta-like 1 homolog (DLK1) in proliferating 3T3-L1 preadipocytes. Exp. Cell Res 2010, 316, 1681–1691. [Google Scholar]
  19. Martinelli, R.; Nardelli, C.; Pilone, V.; Buonomo, T.; Liguori, R.; Castanò, I.; Buono, P.; Masone, S.; Persico, G.; Forestieri, P. miR-519d overexpression is associated with human obesity. Obesity 2012, 18, 2170–2176. [Google Scholar]
  20. Ferland-McCollough, D.; Ozanne, S.; Siddle, K.; Willis, A.; Bushell, M. The involvement of microRNAs in Type 2 diabetes. Biochem. Soc. Trans 2010, 38, 1565. [Google Scholar]
  21. Wang, Q.; Wang, Y.; Minto, A.W.; Wang, J.; Shi, Q.; Li, X.; Quigg, R.J. MicroRNA-377 is up-regulated and can lead to increased fibronectin production in diabetic nephropathy. FASEB J 2008, 22, 4126–4135. [Google Scholar]
  22. Spurlock, M.E.; Gabler, N.K. The development of porcine models of obesity and the metabolic syndrome. J. Nutr 2008, 138, 397–402. [Google Scholar]
  23. Li, M.; Xia, Y.; Gu, Y.; Zhang, K.; Lang, Q.; Chen, L.; Guan, J.; Luo, Z.; Chen, H.; Li, Y. MicroRNAome of porcine pre-and postnatal development. PLoS One 2010, 5, e11541. [Google Scholar]
  24. Xie, S.-S.; Li, X.-Y.; Liu, T.; Cao, J.-H.; Zhong, Q.; Zhao, S.-H. Discovery of porcine microRNAs in multiple tissues by a Solexa deep sequencing approach. PLoS One 2011, 6, e16235. [Google Scholar]
  25. Kozomara, A.; Griffiths-Jones, S. miRBase: Integrating microRNA annotation and deep-sequencing data. Nucleic Acids Res 2011, 39, D152–D157. [Google Scholar]
  26. Qin, L.; Chen, Y.; Niu, Y.; Chen, W.; Wang, Q.; Xiao, S.; Li, A.; Xie, Y.; Li, J.; Zhao, X. A deep investigation into the adipogenesis mechanism: Profile of microRNAs regulating adipogenesis by modulating the canonical Wnt/β-catenin signaling pathway. BMC Genomics 2010, 11, 320. [Google Scholar]
  27. Kim, Y.J.; Hwang, S.J.; Bae, Y.C.; Jung, J.S. MiR-21 regulates adipogenic differentiation through the modulation of TGF-β signaling in mesenchymal stem cells derived from human adipose tissue. Stem Cells 2009, 27, 3093–3102. [Google Scholar]
  28. Hoekstra, M.; van der Lans, C.A.; Halvorsen, B.; Gullestad, L.; Kuiper, J.; Aukrust, P.; van Berkel, T.J.; Biessen, E.A. The peripheral blood mononuclear cell microRNA signature of coronary artery disease. Biochem. Biophys. Res. Commun 2010, 394, 792–797. [Google Scholar]
  29. Ro, S.; Park, C.; Young, D.; Sanders, K.M.; Yan, W. Tissue-dependent paired expression of miRNAs. Nucleic Acids Res 2007, 35, 5944–5953. [Google Scholar]
  30. Esau, C.; Kang, X.; Peralta, E.; Hanson, E.; Marcusson, E.G.; Ravichandran, L.V.; Sun, Y.; Koo, S.; Perera, R.J.; Jain, R. MicroRNA-143 regulates adipocyte differentiation. J. Biol. Chem 2004, 279, 52361–52365. [Google Scholar]
  31. Karbiener, M.; Fischer, C.; Nowitsch, S.; Opriessnig, P.; Papak, C.; Ailhaud, G.; Dani, C.; Amri, E.-Z.; Scheideler, M. microRNA miR-27b impairs human adipocyte differentiation and targets PPARγ. Biochem. Biophys. Res. Commun 2009, 390, 247–251. [Google Scholar]
  32. Carlsbecker, A.; Lee, J.-Y.; Roberts, C.J.; Dettmer, J.; Lehesranta, S.; Zhou, J.; Lindgren, O.; Moreno-Risueno, M.A.; Vatén, A.; Thitamadee, S. Cell signalling by microRNA165/6 directs gene dose-dependent root cell fate. Nature 2010, 465, 316–321. [Google Scholar]
  33. Romualdi, C.; Bortoluzzi, S.; d’Alessi, F.; Danieli, G.A. IDEG6: A web tool for detection of differentially expressed genes in multiple tag sampling experiments. Physiol. Genomics 2003, 12, 159–162. [Google Scholar]
  34. Okayama, H.; Schetter, A.; Harris, C. MicroRNAs and inflammation in the pathogenesis and progression of colon cancer. Dig. Dis 2012, 30, 9–15. [Google Scholar]
  35. Olivieri, F.; Spazzafumo, L.; Santini, G.; Lazzarini, R.; Albertini, M.C.; Rippo, M.R.; Galeazzi, R.; Abbatecola, A.M.; Marcheselli, F.; Monti, D. Age-related differences in the expression of circulating microRNAs: miR-21 as a new circulating marker of inflammaging. Mech. Ageing Dev 2012, 133, 675–685. [Google Scholar]
  36. Dai, L.; Zhang, X.; Hu, X.; Zhou, C.; Ao, Y. Silencing of microRNA-101 prevents IL-1b-induced extracellular matrix degradation in chondrocytes. Arthritis Res. Ther 2012, 14, R268. [Google Scholar]
  37. Trajkovski, M.; Hausser, J.; Soutschek, J.; Bhat, B.; Akin, A.; Zavolan, M.; Heim, M.H.; Stoffel, M. MicroRNAs 103 and 107 regulate insulin sensitivity. Nature 2011, 474, 649–653. [Google Scholar]
  38. Stittrich, A.-B.; Haftmann, C.; Sgouroudis, E.; Kühl, A.A.; Hegazy, A.N.; Panse, I.; Riedel, R.; Flossdorf, M.; Dong, J.; Fuhrmann, F. The microRNA miR-182 is induced by IL-2 and promotes clonal expansion of activated helper T lymphocytes. Nat. Immunol 2010, 11, 1057–1062. [Google Scholar]
  39. Sarver, A.L.; Li, L.; Subramanian, S. MicroRNA miR-183 functions as an oncogene by targeting the transcription factor EGR1 and promoting tumor cell migration. Cancer Res 2010, 70, 9570–9580. [Google Scholar]
  40. Wiklund, E.D.; Bramsen, J.B.; Hulf, T.; Dyrskjøt, L.; Ramanathan, R.; Hansen, T.B.; Villadsen, S.B.; Gao, S.; Ostenfeld, M.S.; Borre, M. Coordinated epigenetic repression of the miR-200 family and miR-205 in invasive bladder cancer. Int. J. Cancer 2011, 128, 1327–1334. [Google Scholar]
  41. Isacke, C. MicroRNA-200 family modulation in distinct breast cancer phenotypes. PLoS One 2012, 7, e47709. [Google Scholar]
  42. Garzon, R.; Volinia, S.; Liu, C.-G.; Fernandez-Cymering, C.; Palumbo, T.; Pichiorri, F.; Fabbri, M.; Coombes, K.; Alder, H.; Nakamura, T. MicroRNA signatures associated with cytogenetics and prognosis in acute myeloid leukemia. Blood 2008, 111, 3183–3189. [Google Scholar]
  43. Zaragosi, L.-E.; Wdziekonski, B.; Brigand, K.L.; Villageois, P.; Mari, B.; Waldmann, R.; Dani, C.; Barbry, P. Small RNA sequencing reveals miR-642a-3p as a novel adipocyte-specific microRNA and miR-30 as a key regulator of human adipogenesis. Genome Biol 2011, 12, R64. [Google Scholar]
  44. Gerin, I.; Bommer, G.T.; McCoin, C.S.; Sousa, K.M.; Krishnan, V.; MacDougald, O.A. Roles for miRNA-378/378* in adipocyte gene expression and lipogenesis. Am. J. Physiol. Endocrinol. Metab 2010, 299, E198–E206. [Google Scholar]
  45. Sun, L.; Xie, H.; Mori, M.A.; Alexander, R.; Yuan, B.; Hattangadi, S.M.; Liu, Q.; Kahn, C.R.; Lodish, H.F. Mir193b-365 is essential for brown fat differentiation. Nat. Cell Biol 2011, 13, 958–965. [Google Scholar]
  46. Yang, K.; He, Y.S.; Wang, X.Q.; Lu, L.; Chen, Q.J.; Liu, J.; Sun, Z.; Shen, W.F. MiR-146a inhibits oxidized low-density lipoprotein-induced lipid accumulation and inflammatory response via targeting toll-like receptor 4. FEBS Lett 2011, 585, 854–860. [Google Scholar]
  47. Krek, A.; Grün, D.; Poy, M.N.; Wolf, R.; Rosenberg, L.; Epstein, E.J.; MacMenamin, P.; da Piedade, I.; Gunsalus, K.C.; Stoffel, M. Combinatorial microRNA target predictions. Nat. Genet 2005, 37, 495–500. [Google Scholar]
  48. Lewis, B.P.; Burge, C.B.; Bartel, D.P. Conserved seed pairing, often flanked by adenosines, indicates that thousands of human genes are microRNA targets. Cell 2005, 120, 15–20. [Google Scholar]
  49. Griffiths-Jones, S.; Saini, H.K.; van Dongen, S.; Enright, A.J. miRBase: Tools for microRNA genomics. Nucleic Acids Res 2008, 36, D154–D158. [Google Scholar]
  50. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc 2008, 4, 44–57. [Google Scholar]
  51. Pruitt, K.D.; Tatusova, T.; Klimke, W.; Maglott, D.R. NCBI Reference Sequences: Current status, policy and new initiatives. Nucleic Acids Res 2009, 37, D32–D36. [Google Scholar]
  52. Gardner, P.P.; Daub, J.; Tate, J.G.; Nawrocki, E.P.; Kolbe, D.L.; Lindgreen, S.; Wilkinson, A.C.; Finn, R.D.; Griffiths-Jones, S.; Eddy, S.R.; et al. Rfam: Updates to the RNA families database. Nucleic Acids Res 2009, 37, D136–D140. [Google Scholar]
  53. Kohany, O.; Gentles, A.J.; Hankus, L.; Jurka, J. Annotation, submission and screening of repetitive elements in Repbase: RepbaseSubmitter and Censor. BMC Bioinforma 2006, 7, 474. [Google Scholar]
  54. Markham, N.R.; Zuker, M. UNAFold: Software for nucleic acid folding and hybridization. Methods Mol. Biol 2008, 453, 3–31. [Google Scholar]
  55. Hellemans, J.; Mortier, G.; de Paepe, A.; Speleman, F.; Vandesompele, J. qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biol 2007, 8, R19. [Google Scholar]
Ijms 14 08611f1 1024
Figure 1. Description of miRNAs in two adipose tissues. (a) Length distribution of the high-quality reads; (b) Distribution of read counts of the identified miRNAs; (c) Distribution of read counts in the three defined miRNA groups; (d) Copy numbers of the top 10 miRNAs with highest read counts. IMAT, intermuscular adipose tissue; sASAT, superficial abdominal subcutaneous adipose tissue.
Figure 1. Description of miRNAs in two adipose tissues. (a) Length distribution of the high-quality reads; (b) Distribution of read counts of the identified miRNAs; (c) Distribution of read counts in the three defined miRNA groups; (d) Copy numbers of the top 10 miRNAs with highest read counts. IMAT, intermuscular adipose tissue; sASAT, superficial abdominal subcutaneous adipose tissue.
Ijms 14 08611f1
Ijms 14 08611f2a 1024Ijms 14 08611f2b 1024
Figure 2. Characteristics of the differentially expressed (DE) miRNAs between porcine sASAT and IMAT. (a) Distribution of 597 unique miRNAs between sASAT (blue) and IMAT (yellow). The red circle represents the 110 miRNAs with read counts >1000 in either of the two libraries. The dashed circles indicate the 45 IMAT-enriched (left) and eight sASAT-enriched (right) miRNAs (p < 0.001); (b) q-PCR validation for the top 14 DE miRNAs with highest read counts between IMAT and sASAT. Pearson’s correlation was used to determine the relationship between the q-PCR and small RNA-seq results for miRNA expression levels. IMAT-NE and sASAT-NE represent normalized expression levels for the miRNAs in the IMAT and sASAT libraries, respectively; (c) The differential expression of 19 inflammation- and diabetes-related miRNAs between IMAT and sASAT.
Figure 2. Characteristics of the differentially expressed (DE) miRNAs between porcine sASAT and IMAT. (a) Distribution of 597 unique miRNAs between sASAT (blue) and IMAT (yellow). The red circle represents the 110 miRNAs with read counts >1000 in either of the two libraries. The dashed circles indicate the 45 IMAT-enriched (left) and eight sASAT-enriched (right) miRNAs (p < 0.001); (b) q-PCR validation for the top 14 DE miRNAs with highest read counts between IMAT and sASAT. Pearson’s correlation was used to determine the relationship between the q-PCR and small RNA-seq results for miRNA expression levels. IMAT-NE and sASAT-NE represent normalized expression levels for the miRNAs in the IMAT and sASAT libraries, respectively; (c) The differential expression of 19 inflammation- and diabetes-related miRNAs between IMAT and sASAT.
Ijms 14 08611f2aIjms 14 08611f2b
Ijms 14 08611f3 1024
Figure 3. KEGG pathways and gene ontology biological process (GO-BP) categories enriched in the target genes of the top eight sASAT- and IMAT-enriched miRNAs. GO-BP is the GO terms under the biological process ontology.
Figure 3. KEGG pathways and gene ontology biological process (GO-BP) categories enriched in the target genes of the top eight sASAT- and IMAT-enriched miRNAs. GO-BP is the GO terms under the biological process ontology.
Ijms 14 08611f3

Share and Cite

MDPI and ACS Style

Ma, J.; Yu, S.; Wang, F.; Bai, L.; Xiao, J.; Jiang, Y.; Chen, L.; Wang, J.; Jiang, A.; Li, M.; et al. MicroRNA Transcriptomes Relate Intermuscular Adipose Tissue to Metabolic Risk. Int. J. Mol. Sci. 2013, 14, 8611-8624. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms14048611

AMA Style

Ma J, Yu S, Wang F, Bai L, Xiao J, Jiang Y, Chen L, Wang J, Jiang A, Li M, et al. MicroRNA Transcriptomes Relate Intermuscular Adipose Tissue to Metabolic Risk. International Journal of Molecular Sciences. 2013; 14(4):8611-8624. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms14048611

Chicago/Turabian Style

Ma, Jideng, Shuzhen Yu, Fengjiao Wang, Lin Bai, Jian Xiao, Yanzhi Jiang, Lei Chen, Jinyong Wang, Anan Jiang, Mingzhou Li, and et al. 2013. "MicroRNA Transcriptomes Relate Intermuscular Adipose Tissue to Metabolic Risk" International Journal of Molecular Sciences 14, no. 4: 8611-8624. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms14048611

Article Metrics

Back to TopTop