Integrating miRNA and mRNA Expression Profiling Uncovers miRNAs Underlying Fat Deposition in Sheep.
Adipocytes are cells that form and store fat globules in the body. Adipose tissue deposits are located in different regions of the body, primarily in the form of subcutaneous fat or intramuscular fat in domestic animals. Though the volume of adipocytes varies from breed to breed and is influenced by diet, genetic factors are considered major determinants for the formation of adipocytes. Fat-tailed sheep, which exhibit distinctive large tails and hindquarters, comprise approximately 25% of the world's sheep population  and are raised commercially for meat, milk, fat, or wool production. Fat-tail is regarded as an adaptive response to the harsh challenges of desert life and the fat deposits provide a valuable energy reservation during drought seasons. The Kazakhstan sheep (KS, fat-tail breed) and Tibetan sheep (TS, thin-tail breed) are two native breeds raised in the extremely arid regions of western China, exhibiting distinct phenotypes of tails as previously described .
MicroRNAs (miRNAs) are a set of small noncoding RNAs (~22 nt) that bind to partially complementary sequences on target mRNAs resulting in posttranscriptional regulation of gene expression . MicroRNAs play key roles in regulating numerous biological processes in developmental, cell differentiation, and disease processes. These include adipogenesis and obesity in humans [4-7], as well as lipid metabolism and fat deposition in livestock [8, 9]. Additionally, evidence suggests that miRNAs regulate the formation of adipose tissue. For instance, miR-378 was strongly associated with the thickness of back fat in beef cattle . In recent years, next generation sequencing has been widely used to explore molecular mechanism of adipogenesis in sheep. However, miRNA-mRNA correlation analysis relating to adipogenesis in sheep has not been examined as in other livestock species [11,12]. This research provides an understanding of the process in sheep and contributes to the understanding of obesity and related diseases.
We have previously described specific changes of transcriptomic differences with respect to adipose tissues from the representative KS and TS breeds and identified a list of candidate genes underlying the phenotypic differences of fat/ thin tails in sheep . In this study, we conducted miRNA sequencing using the same tissues, in attempt to elucidate the roles of miRNAs in the presentation of these two extreme phenotypes. We further employed a systems biology approach to examine the correlation between miRNA and mRNA expression data to identify potential miRNA-associated target genes. Here, we have implicated interactions from particular cellular processes such as the MAPK signaling pathway, FoxO and Wnt signaling pathway, and focal adhesion, gaining insights into the potential roles of miRNAs in the regulation of these pivotal pathways.
2. Materials and Methods
2.1. Ethics Statement. The experiments were conducted following the guidelines of the Animal Ethics Committee at Northwest A&F University under the document 2011-31101684. The sampling procedures complied with the "Guidelines on Ethical Treatment of Experimental Animals" (2006) Number 398 set by the Ministry of Science and Technology, China.
2.2. Animals and Phenotypes. Six adult individuals (three males and three females, two years old) from KS and TS breeds were randomly selected, respectively, and any two or more individuals with a traceable phylogenetic relationship were avoided in the sampling process. Adipose tissue biopsies were collected from the tails and stored as previously described .
2.3. Small RNA Library Construction and Sequencing. Total RNA from the mixed adipose tissues of six adult sheep was isolated using the RNAiso plus kit (TaKaRa, Dalian, China) according to the manufacturer's protocol. The small RNA libraries were prepared following Illumina[R] TruSeq[TM] Small RNA Sample Preparation protocol. The small RNA libraries were sequenced using an Illumina/Solexa 1 G Genome Analyzer System (BGI, Shenzhen). The sequencing reads were deposited to the Sequence Read Archive, under accession code SRP093866.
2.4. Expression Profiling. The raw miRNA sequencing reads were filtered and only the unique mapping reads were used for subsequent bioinformatics analysis, such as annotation and gene expression. The expression of miRNA in two samples (TS and KS) was normalized by transcripts per million (TPM) as previous description . After normalization, the expression value of a miRNA, more than one in the two samples, was kept for the differential expression analysis, and we revised 0 value to be 0.01 when we calculated fold change of the two breeds . The differentially expressed genes (DEGs) were determined according to its absolute fold change ([log.sub.2](TS/KS)) of Th and P value of [??] 0.05. Scatter plots were used to demonstrate differentially expressed miRNA between the two sheep phenotypes.
2.5. Quantitative RT-PCR. Total RNA from the mixed sheep tails (six from KS and six from TS) was isolated using a miRNA extraction kit (CWBIO, Beijing, China). First strand cDNA was synthesised with reference to RevertAid First Strand cDNA Synthesis kit (Thermo Scientific #K1622, Fermentas) manufacturer's instructions, using special stem-loop RT primers  (Additional file 1 in Supplementary Material available online at https://doi.org/10.1155/2017/1857580). PCR-iQ5 system (Bio-Rad, Hercules, CA, USA) was used to analyse miRNA genes expression level. The reaction mixture volume is 25 [micro]L, containing 2 pL of cDNA (25 ng), 12.5 [micro]L SYBR Premix Ex TaqTM II (TaKaRa, Dalian, China), 1.0 [micro]L specific forward primer, 1 [micro]L universal primer, and 8.5 [micro]L water, the reaction program referenced to previous description . U6 snRNA was used as an endogenous control. All reactions were performed in triplicate, and [2.sup.-[DELTA][DELTA]Ct] method was used to analyse the data. t-tests (parametric test, unpaired test) was used to determine significance of the results using SPSS21. Graphpad Prism 6 was used to plot the graphs. All data were presented as mean [+ or -] SD and values of P < 0.05 were considered statistically significant.
2.6. mRNA Profiling Bioinformatics. The sequences of reads with poly(A) tracts were deposited to the Sequence Read Archive, under accession code SRP093860. The FastQC package was used to assess the quality of raw reads and then mapped to the sheep genome assembly Oar_v3.1 (GCA_000298735.1) using TopHat2 . Quality control was performed on the two samples. The high-quality (HQ) clean reads were presented after a series of procedures (remove adapter, ploy N, and low-quality reads). All the following analyses were based on the HQ clean reads. The unannotated HQ clean reads were assembled by Cufflinks . Orthonormal expression levels were calculated as FPKM (Fragments Per Kilobase of transcript per Million mapped reads) , and FPKM distribution of the samples was presented as gene's coverage. The edgeR package for bioconductor was applied to identify differentially expressed genes between TS and KS . Genes were defined as differentially expressed with the threshold FDR < 0.05 and [absolute value of [log.sub.2](fold change)] > 1.
2.7. Integration Analysis of miRNA and mRNA Sequencing Data. MicroRNA-mRNA integration analysis was performed on the mRNAs that were differentially expressed between the KS and TS breeds in our previous study , and the miRNAs were identified and classified according to the miRNA family again, in order to meet the current study. Individual miRNAs determined to be differentially expressed were submitted to each of the three databases ((RNAhybridv2.1.2)+svm_light(v6.01), Miranda(v3.3a), and TargetScan(Version: 7.0)) to determine the predicted mRNA targets, and the parameter was taken as the default. We selected only mRNA targets that were differentially expressed in the same samples. Next, the correlation between the expression levels of each miRNA and any of the predicted, differentially expressed target mRNA from any of the three databases was computed. According to the analysis results and the prior research, we chose miRNA-125a-5p and one of the target genes (ESRRa) for validation analysis.
2.8. Functional Enrichment of Target Genes. For each miRNA with a significant association with its set of mRNA targets (based on the gene set enrichment test), the GO and KEGG terms were extracted for all of its mRNA targets. We reveal the functions significantly related to predicted target gene candidates of miRNAs used gene ontology (GO) analysis.
The primary pathways of the candidate target genes were determined by KEGG pathway analysis. The formula used to calculate this is the same as that used in the GO analysis. Genes with P < 0.05 are considered as significantly enriched in target gene candidates. The KEGG indicates the main pathways in which the candidate genes are involved.
3.1. Small RNA Diversity in Sheep Adipose Tissues. The sequencing of two small RNA libraries from the KS and TS yielded 12 M counts of sequencing reads, respectively (Table 1). All of the clean reads were aligned with the sheep genome (Oar_v3.1) sequence using SOAP software. The length distribution and the distribution of miRNAs on the genome were analysed. The high-quality reads in both groups exhibited the canonical size range distribution that is common to mammalian miRNAs  (Figure 1(a)). The vast majority of the reads were approximately 21~23 nucleotides (nt) in length, and ~22 nt reads accounted for 50% and 39% for each of high-quality reads from the KS and TS breeds, respectively. This is in agreement with previous founding in other adipose tissues . These results demonstrated the high reliability of using the small RNA-sequencing approach to obtain miRNA reads for further studies.
There were 8,000,311 (67.53%) and 6,720,288 (56.68%) reads matched to the sheep genome sequence in KS and TS libraries, respectively. All clean reads were classified and annotated by tag2annotation software (developed by BGI), aligned against the Rfam database (Ver. 10.1) and the miRBase (Ver. 19.0). There were 11,599,356 and 11,603,988 total conserved miRNAs reads, and 0.77% and 0.81% of clean reads were identified as potential novel miRNAs from KS and TS libraries, respectively, and more information was given in the additional file (Additional file 2).
3.2. Differential miRNA Expression in Sheep Adipose Tissues. To identify conserved miRNAs in sheep, all sRNA sequences were mapped to known miRNAs in the miRBase 19.0 database. We presented a total of 815 miRNAs found in abundance in adipose tissues of sheep (Additional file 3). Of 815 identified miRNAs, 539 miRNAs were expressed in both breeds, whereas 179 and 97 miRNAs were specifically expressed in KS and TS, respectively. The top ten abundant miRNAs from KS, TS, and both groups are listed in Table 2. Of the shared miRNAs, 93 were upregulated and 33 were downregulated in the adult adipose tissues of sheep (Figure 1(b)).
We next analysed the differentially expressed miRNAs in the two breeds. One hundred and seventy-five miRNAs exhibited similar expression levels (nonsignificant), 33 miRNAs showed significantly reduced expression (P < 0.05) and 276 miRNAs showed an increase expression (P < 0.05) in KS, compared to TS (Figure 1(b)). MicroRNAs with similar expression patterns in different sample pairs were clustered together. Clustering analysis was based on the sample difference model by using Cluster software, and the results were viewed with Java Treeview. All differentially expressed miRNAs clustered together after six rounds of clustering (Figure 2).
Novel miRNAs could be predicted through the characteristic hairpin structure of miRNA precursors. We predicted novel miRNAs by exploring secondary structure, the Dicer cleavage site, and the minimum free energy of the unannotated small RNA reads, which could be mapped to sheep genome sequences by using Mireap software. In total, 186,755 unannotated sequences were used to predict novel miRNAs. Thirty-five potential novel miRNAs were identified; 7 miRNAs were expressed in the two samples, whereas 18 and 10 miRNAs were specifically expressed in KS and TS, respectively (Additional file 4).
3.3. Validation of miRNA by qPCR. In order to validate the reliability of miRNA expression data, quantitative RTPCR was conducted to assess the expression of six upregulated miRNAs (miR-27c, miR-183, miR-224, miR-339a, miR-2070-3p, and miR-101) and four downregulated ones (miR-27d, miR-106a, miR-204, and miR-301) among the differentially expressed miRNAs, using adipose tissues from both KS and TS groups. U6 was used as an internal control. Stable miRNA expression levels in both KS and TS tissues correlated well between qPCR and miRNA-seq data (Figure 3).
3.4. Differential Expression of mRNAs in TS and KS Sheep Adipose Tissue. In total, 55,354,542 and 51,410,666 clean reads were obtained from the TS and KS adipose tissues, respectively. The percentages of reads containing HQ clean reads, adaptors, N, low-quality reads, poly A, and clean reads were calculated. Approximately 99% of the clean reads passed the filter (Additional file 5). After alignment by TopHat, 86.78% and 86.65% unique reads were found in Tibetan and Kazakhstan sheep, respectively (Additional file 6). We identified 15330 (76.13%) known genes and 458 novel genes in Tibetan sheep and 15162 (75.29%) known genes and 459 novel genes in Kazakhstan sheep. Those expression levels were quantified by FPKM > 0.05 (Additional file 7).
DEGs were selected out by FDR < 0.05 and FPKM fold change (|[log.sub.2](KS/TS)|). Compared with KS, 1931 DEGs were identified. Seven hundred and sixty-six of these were upregulated expression, and 1165 were downregulated (Additional file 8). The scatter plot showed significant DEGs in the above comparisons between TS and KS (Figure 4).
3.5. GO and Pathway Enrichment Analyses of Differentially Expressed Genes. To determine the function of differentially expressed genes, all of the DEGs in this study were mapped to terms in the GO database (http://www.geneontology.org/) using Blast2GO . In order to assess the general functional characteristics of genes activity in sheep adipose tissues, we performed functional enrichment analysis of GO using DAVID software and pathway enrichment using specific KEGG terms . A total of 1522 genes were categorized into the three main categories of GO classification. Figure 5 shows the GO classification of DEGs.
KEGG pathway annotation showed that 4899 genes and 537 DEGs genes were annotated for 209 biological functions (Additional file 9). Figure 6 showed the top 20 of pathway enrichment. We found that several significantly overrepresented categories of GO biological processes were associated with fat metabolism and deposition. We found that the cAMP signaling pathway had a critical role in both adipogenesis and lipid partitioning in white adipose tissue and well-characterized mechanisms controlling adipocyte differentiation [23, 24].
3.6. Identification of Potential mRNA Targets of miRNAs. MicroRNA recognizes its target mRNA through binding to a seed sequence, which localizes on the 3' untranslated regions (37UTR) of target mRNAs . In an attempt to determine potential genes whose mRNAs might be targeted by particular miRNAs, we conducted an integrated analysis between the expression of specific miRNA and the expression of all the predicted miRNA target genes (mRNA) using TargetScan. We identified potential miRNA target genes that are downregulated at the transcriptional level and are inversely correlated with the miRNA expression in the same KS and TS adipose tissues. Table 3 lists the top 10 differentially expressed miRNAs and their target genes determined by miRNAmRNA correlation analysis, and the details were presented in Additional file 10.
3.7. Validation of the Target Gene of miR-125a-5p. For validation, we selected the a priori gene ESRR[alpha], which inhibits preadipocyte differentiation [26, 27]. Using RNAhybrid (http://bibiserv.techfak.uni-bielefeld.de/rnahybrid) and miRanda (http://www.microrna.org) and combining the results of the analyses, we presumed that miR-125a-5p regulates ESRRa translation. MiR-125a-5p was found to directly target the ESRRa 3' UTR sequence (Figure 7(a)). To verify this finding, the 3' UTR containing miR-125a-5p targeted site was cloned and inserted into the psiCHECK[TM]-2 reporter plasmid (psiCHECK[TM]-2-ESRR[alpha]-3'UTR WT) and also the mutated 3' UTR (psiCHECK[TM]-2-ESRRa-3'UTR-mutant). The two types were cotransfected with the miR-125a-5p mimics or negative control sample into 293T cells, respectively, and luciferase activity was measured after 48 h later. The luciferase activity of the miR-125a-5p group was significantly lower than that of the NC group (P < 0.01), whereas the mutated ESRRa 3'UTR exhibited increased luciferase expression (Figure 7(c)). Thus, the miR-125a-5p was confirmed as one of the ESRR[alpha] translation targeted regulation sequences and the primer sequences used for this are listed in Additional file 1.
3.8. Functional Enrichment Based on miRNA-Associated Target Genes. In order to assess the general functional characteristics of the miRNA-related activity in sheep adipose tissues, we selected the putative miRNA target genes derived from miRNA-mRNA correlation analysis and performed functional enrichment analysis of gene ontology (GO) using DAVID software and pathway enrichment using specific KEGG terms . We found that several significantly overrepresented categories of GO biological processes were associated with biological functions, including environmental information processing, cellular processes, metabolism, and genetic information processing (Additional file 11). KEGG pathway annotation showed that 4,687 target genes were annotated for 226 biological functions (Additional file 12). Most of these genes were involved in metabolic processes, regulation of biological process, and signal transduction. Several of these KEGG terms were related to fat deposition and fatty acid metabolism, including MAPK signaling pathway, focal adhesion, Pyruvate metabolism, FoxO signaling pathway, and TNF signaling pathway. This suggests that miRNAs are functionally relevant to lipid metabolism and fat deposition processes.
Next generation sequencing approaches have become powerful tools to predict and identify the novel genes and small RNAs of livestock species and other organisms. Compared with other meat-producing livestock such as cattle and pig, fewer studies have focused on the role of miRNAs in the process of fat deposition in sheep, though there are a few omics studies screening miRNAs in sheep skin  and muscles . In an effort to gain insights into the molecular mechanisms underlying fat deposition of adipose tissue governed by miRNA regulation in a fat-tailed and a thin-tailed sheep breeds, we evaluated whole transcriptome profiles of miRNAs and mRNA expression on the same adipose tissue from sheep tails .
In the present study, we were able to discover 35 novel miRNAs which have eluded previous efforts and identify 239 miRNAs that have not been previously annotated in sheep. Our analysis detected 815 differentially expressed miRNAs between the fat deposits of the two sheep breeds used here. Not surprisingly, a proportion of the miRNAs was found to be differentially expressed in our study had been previously reported to play a role in adiposity, adipocyte development/differentiation, and other metabolic disturbances in livestock, as well as adiposity and related disease risk in humans. This further highlights the importance of miRNAs in adipose tissue development and metabolism. Previous studies have shown an increase in miR-143 during human and murine preadipocyte differentiation, and the following experiment verified that it inhibited preadipocyte differentiation [30-32]. miR-378 may promote bovine adipogenesis in white adipose tissue through targeting MAPK1 and PPARy . The interactions demonstrated effects on the MAPK signaling pathway, focal adhesion, and Pyruvate metabolism pathways. miR-103 and miR-30 have been reported to enhance adipogenesis [32-35], while miR-27 and miR-138 could inhibit adipogenesis [36, 37]. miR-122, miR-370, and miR-378 have been demonstrated to have a key role in lipid metabolism [7, 38, 39]. miR-148a might modulate fat deposition by targeting MAPKAPK5, MAPK3, and MAP2K2, and miR-148 has been shown to affect obesity and modulates adipocyte differentiation [40, 41]. miR-125b-5p inhibited proliferation and promoted adipogenic differentiation in 3T3-L1 preadipocytes , and other studies have shown that miR-125b impaired brite adipocyte formation and function . These reports are consistent with our findings reported.
In addition to microRNA expression profiling, we performed miRNA-mRNA integration analysis using the miRNA and mRNA expression data from corresponding tissues in distinct sheep breeds and determined a limited number of putative target genes for the differentially expressed miRNAs. The integration analysis removed a large portion of false positive miRNA target genes and increased the precision of our predictions . Therefore, the identified candidates could be taken as reference genes to validate the differentially expressed miRNAs really affecting the mRNA expression levels during adiposeness. Here, we showed that miR-125a inhibits porcine preadipocytes differentiation by targeting ESRRa . This kind of targeted effect also exists in sheep, so we speculate that the function of miR-125a in ovine preadipocytes is similar to that in porcine preadipocytes.
The DAVID gene annotation analysis of the miRNA-associated targets indicated that the GO biological processes were enriched to modification processes, metabolic processes, signal transduction, and regulation of biological processes. This demonstrates the distinct roles of adipocyte development in the tail areas as the reserves of energy, in response to external or internal environmental conditions during the dry and cold seasons. KEGG pathway annotation of their putative miRNA-related targets revealed that the MAPK signaling pathway, focal adhesion, Pyruvate metabolism, FoxO signaling pathway, and TNF signaling pathway contribute to the formation of sheep fat/thin-tails phenotypes. Various tail phenotypes are unique for sheep compared to other animals in the word, interestingly, a mass of adipose depot in the fat-tail sheep breed and less in the thin-tail sheep breed. In a manner, this phenomenon is a form of energy storage in the fight against extreme conditions, and the function of tail adipose tissue is similar to subcutaneous adipose tissue of other animals according to the enrich pathways, which is associated with metabolism risk [20,45].
In conclusion, we identified a number of miRNAs that are differentially expressed between the fat-tailed and short-tailed sheep breeds. We further highlighted gene targets of related miRNAs that maybe involved in regulating fat deposition and adiposeness, in sheep and other livestock. This occurs via the key signaling pathways including focal adhesion, Pyruvate metabolism, and the MAPK, FoxO, and TNF signaling pathway. Further studies are needed to verify the correlation between key miRNAs and their target genes by in vitro approach and elucidate the functional impacts that miRNAs serve during adiposeness. Our results also provide evidence for the interaction of miRNAs and genes in the regulation of obesity and metabolic syndromes, which suggests that this may serve as an animal model for human' obesity and metabolic syndromes researches.
The authors declare that they have no competing interests.
Guangxian Zhou, Xiaolong Wang, and Chao Yuan analysed the data and wrote the paper. Guangxian Zhou performed all experiments, and Danju Kang participated in the animal experiments and miRNA data analysis. Xiaochun Xu and Rongqing Geng participated in data analysis and manuscript revision. Jiping Zhou, Yuxin Yang, and Zhaoxia Yang gave some advice on the experiment design and pathologic analysis. Yulin Chen supervised and provided continuous guidance for the experiment and the data analysis. Guangxian Zhou and Xiaolong Wang contributed equally to this work.
This work was supported by grants from the NSFC (nos. 31372279 and 31402038), Natural Science Foundation (2014JM3068), and International Sci-Tech Cooperation Program (2014KW14-01), of Shaanxi province. Thanks are due to Gene Denovo Co., Ltd., for providing technical support. The authors also thank the 34th Conference of the International Society of Animal Genetics (ISAG) for giving them an opportunity to present this work.
 A. Davidson, Oxford Companion to Food, Oxford University Press, Oxford, UK, 1999.
 X. Wang, G. Zhou, X. Xu et al., "Transcriptome profile analysis of adipose tissues from fat and short-tailed sheep," Gene, vol. 549, no. 2, pp. 252-257, 2014.
 D. P. Bartel, "MicroRNAs: genomics, biogenesis, mechanism, and function," Cell, vol. 116, no. 2, pp. 281-297, 2004.
 H. M. Heneghan, N. Miller, and M. J. Kerin, "Role of microRNAs in obesity and the metabolic syndrome," Obesity Reviews, vol. 11, no. 5, pp. 354-361, 2010.
 R. Alexander, H. Lodish, and L. Sun, "MicroRNAs in adipogenesis and as therapeutic targets for obesity," Expert Opinion on Therapeutic Targets, vol. 15, no. 5, pp. 623-636, 2011.
 R. A. McGregor and M. S. Choi, "microRNAs in the regulation of adipogenesis and obesity," Current Molecular Medicine, vol. 11, no. 4, pp. 304-316, 2011.
 C. Hilton, M. J. Neville, and F. Karpe, "MicroRNAs in adipose tissue: their role in adipogenesis and obesity," International Journal of Obesity, vol. 37, no. 3, pp. 325-332, 2013.
 H.-C. Liu, J. A. Hicks, N. Trakooljul, and S.-H. Zhao, "Current knowledge of microRNA characterization in agricultural animals," Animal Genetics, vol. 41, no. 3, pp. 225-231, 2010.
 A. Fatima and D. G. Morris, "MicroRNAs in domestic livestock," Physiological Genomics, vol. 45, no. 16, pp. 685-696, 2013.
 W. Jin, M. V. Dodson, S. S. Moore, J. A. Basarab, and L. L. Guan, "Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development," BMC Molecular Biology, vol. 11, article 29, 2010.
 A. Eirin, S. M. Riester, X.-Y. Zhu et al., "MicroRNA and mRNA cargo of extracellular vesicles from porcine adipose tissue-derived mesenchymal stem cells," Gene, vol. 551, no. 1, pp. 55-64, 2014.
 H. Wang, Y. Zheng, G. Wang, and H. Li, "Identification of microRNA and bioinformatics target gene analysis in beef cattle intramuscular fat and subcutaneous fat," Molecular BioSystems, vol. 9, no. 8, pp. 2154-2162, 2013.
 B. Li and C. N. Dewey, "RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome," BMC Bioinformatics, vol. 12, article 323, 2011.
 C. Yuan, X. Wang, R. Geng, X. He, L. Qu, and Y. Chen, "Discovery of cashmere goat (Capra hircus) microRNAs in skin and hair follicles by Solexa sequencing," BMC Genomics, vol. 14, article 511, 2013.
 C. Chen, D. A. Ridzon, A. J. Broomer et al., "Real-time quantification of microRNAs by stem-loop RT-PCR," Nucleic Acids Research, vol. 33, no. 20, article no. e179, 2005.
 D. Kim, G. Pertea, C. Trapnell, H. Pimentel, R. Kelley, and S. L. Salzberg, "TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions," Genome Biology, vol. 14, no. 4, article R36, 2013.
 C. Trapnell, B. A. Williams, G. Pertea et al., "Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation," Nature Biotechnology, vol. 28, no. 5, pp. 511-515, 2010.
 A. Mortazavi, B. A. Williams, K. McCue, L. Schaeffer, and B. Wold, "Mapping and quantifying mammalian transcriptomes by RNA-Seq," Nature Methods, vol. 5, no. 7, pp. 621-628, 2008.
 M. D. Robinson, D. J. McCarthy, and G. K. Smyth, "edgeR: a bioconductor package for differential expression analysis of digital gene expression data," Bioinformatics, vol. 26, no. 1, pp. 139-140, 2010.
 J. Ma, S. Yu, F. Wang et al., "MicroRNA transcriptomes relate intermuscular adipose tissue to metabolic risk," International Journal of Molecular Sciences, vol. 14, no. 4, pp. 8611-8624, 2013.
 A. Conesa, S. Gotz, J. M. Garcia-Gomez, J. Terol, M. Talon, and M. Robles, "Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research," Bioinformatics, vol. 21, no. 18, pp. 3674-3676, 2005.
 M. Kanehisa and S. Goto, "KEGG: kyoto encyclopedia of genes and genomes," Nucleic Acids Research, vol. 28, no. 1, pp. 27-30, 2000.
 M. Chen, H. Chen, A. Nguyen et al., "GSa deficiency in adipose tissue leads to a lean phenotype with divergent effects on cold tolerance and diet-induced thermogenesis," Cell Metabolism, vol. 11, no. 4, pp. 320-330, 2010.
 K. Ravnskjaer, A. Madiraju, and M. Montminy, "Role of the cAMP pathway in glucose and lipid metabolism," Handbook of Experimental Pharmacology, vol. 233, pp. 29-49, 2015.
 M. W. Jones-Rhoades and D. P. Bartel, "Computational identification of plant MicroRNAs and their targets, including a stress-induced miRNA," Molecular Cell, vol. 14, no. 6, pp. 787-799, 2004.
 H.-L. Ji, C.-C. Song, Y.-F. Li et al., "miR-125a inhibits porcine preadipocytes differentiation by targeting ERRa," Molecular and Cellular Biochemistry, vol. 395, no. 1-2, pp. 155-165, 2014.
 K. Lapid, A. Lim, D. J. Clegg, D. Zeve, and J. M. Graff, "Oestrogen signalling in white adipose progenitor cells inhibits differentiation into brown adipose and smooth muscle cells," Nature Communications, vol. 5, article 5196, 2014.
 G. Liu, R. Liu, Q. Li et al., "Identification of microRNAs in wool follicles during anagen, catagen, and telogen phases in tibetan sheep," PLoS ONE, vol. 8, no. 10, Article ID e77801, 2013.
 S. Zhang, F. Zhao, C. Wei et al., "Identification and characterization of the miRNA transcriptome of Ovis aries," PLoS ONE, vol. 8, no. 3, Article ID e58905, 2013.
 C. Esau, X. Kang, E. Peralta et al., "MicroRNA-143 regulates adipocyte differentiation," Journal of Biological Chemistry, vol. 279, no. 50, pp. 52361-52365, 2004.
 R. Takanabe, K. Ono, Y. Abe et al., "Up-regulated expression of microRNA-143 in association with obesity in adipose tissue of mice fed high-fat diet," Biochemical and Biophysical Research Communications, vol. 376, no. 4, pp. 728-732, 2008.
 H. Xie, B. Lim, and H. F. Lodish, "MicroRNAs induced during adipogenesis that accelerate fat cell development are downregulated in obesity," Diabetes, vol. 58, no. 5, pp. 1050-1057, 2009.
 H. Enomoto, T. Furuichi, A. Zanma et al., "Runx2 deficiency in chondrocytes causes adipogenic changes in vitro," Journal of Cell Science, vol. 117, no. 3, pp. 417-425, 2004.
 M. Karbiener, C. Neuhold, P. Opriessnig, A. Prokesch, J. G. Bogner-Strauss, and M. Scheideler, "MicroRNA-30c promotes human adipocyte differentiation and co-represses PAI-1 and ALK2," RNA biology, vol. 8, no. 5, pp. 850-860, 2011.
 L.-E. Zaragosi, B. Wdziekonski, K. L. Brigand et al., "Small RNA sequencing reveals miR-642a-3p as a novel adipocytespecific microRNA and miR-30 as a key regulator of human adipogenesis," Genome Biology, vol. 12, no. 7, article R64, 2011.
 Q. Lin, Z. Gao, R. M. Alarcon, J. Ye, and Z. Yun, "A role of miR-27 in the regulation of adipogenesis," FEBS Journal, vol. 276, no. 8, pp. 2348-2358, 2009.
 M. Trajkovski, J. Hausser, J. Soutschek et al., "MicroRNAs 103 and 107 regulate insulin sensitivity," Nature, vol. 474, no. 7353, pp. 649-653, 2011.
 I. Gerin, G. T. Bommer, C. S. McCoin, K. M. Sousa, V. Krishnan, and O. A. MacDougald, "Roles for miRNA-378/378* in adipocyte gene expression andlipogenesis," American Journal of Physiology - Endocrinology and Metabolism, vol. 299, no. 2, pp. E198-E206, 2010.
 C. Fernandez-Hernando, Y. Suarez, K. J. Rayner, and K. J. Moore, "MicroRNAs in lipid metabolism," Current Opinion in Lipidology, vol. 22, no. 2, pp. 86-92, 2011.
 C. Shi, M. Zhang, M. Tong et al., "MiR-148a is associated with obesity and modulates adipocyte differentiation of mesenchymal stem cells through wnt signaling," Scientific Reports, vol. 5, article 9930, 2015.
 D. Pan, C. Mao, B. Quattrochi et al., "MicroRNA-378 controls classical brown fat expansion to counteract obesity," Nature Communications, vol. 5, article 4725, 2014.
 D. Ouyang, Y. Ye, D. Guo et al., "MicroRNA-125b-5p inhibits proliferation and promotes adipogenic differentiation in 3T3-L1 preadipocytes," Acta Biochimica et Biophysica Sinica, vol. 47, no. 5, pp. 355-361, 2015.
 M. Giroud, D. F. Pisani, M. Karbiener et al., "miR-125b affects mitochondrial biogenesis and impairs brite adipocyte formation and function," Molecular Metabolism, vol. 5, no. 8, pp. 615-625, 2016.
 J. C. Huang, T. Babak, T. W. Corson et al., "Using expression profiling data to identify human microRNA targets," Nature Methods, vol. 4, no. 12, pp. 1045-1049, 2007
 M. M. Ibrahim, "Subcutaneous and visceral adipose tissue: structural and functional differences," Obesity Reviews, vol. 11, no. 1, pp. 11-18, 2010.
Guangxian Zhou, (1) Xiaolong Wang, (1) Chao Yuan, (1,2) Danju Kang, (1) Xiaochun Xu, (1,3) Jiping Zhou, (4) Rongqing Geng, (1,5) Yuxin Yang, (1) Zhaoxia Yang, (1) and Yulin Chen (1)
(1) College of Animal Science and Technology, Northwest A&F University, Yangling, China
(2) Lanzhou Institute of Husbandry and Pharmaceutical Sciences of CAAS, Lanzhou, China
(3) College of Biology Sciences and Engineering, Beifang University of Nationalities, Yinchuan, Ningxia 750021, China
(4) Qinghai Academy of Animal Science and Veterinary Medicine, Xining, China
(5) College of Pharmacy, Yancheng Teachers University, Yancheng, China
Correspondence should be addressed to Yulin Chen; email@example.com
Received 21 August 2016; Revised 5 December 2016; Accepted 15 December 2016; Published 15 February 2017
Academic Editor: Gang Liu
Caption: Figure 1: (a) Length distribution of the clean reads based on total abundance and distinct sequences. (b) Differential expression of miRNAs between the adipose tissues from KS and TS breed. Each point represents a single miRNA. The X and Y axes show the expression level of miRNAs in these two samples, respectively. Red points represent miRNAs with a ratio > 2, blue points represent miRNAs with 1/2 < ratio [??] 2, and green points represent miRNAs with ratio [??] 1/2; ratio = normalized expression in treatment/normalized expression in control.
Caption: Figure 2: Heat map of miRNAs differentially expressed in sheep adipose tissues. Red indicates that the miRNA has a higher expression level in the treatment samples; green indicates that the miRNA has a higher expression in the control samples and gray indicates that the miRNA has no expression in at least one sample. Each row in the figure represents one miRNA, and each column shows one sample pair. Each cell shows the differential expression of a miRNA in one sample pair.
Caption: Figure 3: qPCR validation of miRNAs identified in sheep adipose tissues. Red represents TS and blue represents KS (0.01 < * p < 0.05, ** p < 0.01).
Caption: Figure 4: Differentially expressed gene between TS and KS DEGs. Red spots represented upregulated genes, and green spots indicated downregulated genes. Back spots represented genes that did not show obvious changes between the TS and KS.
Caption: Figure 5: GO classification of DEG. The x-axis indicated the subcategories; the right 7-axis indicated the number of DEG. The height of columns represented the percentage (number) of upregulated expression (red) and downregulated expression (green) genes.
Caption: Figure 6: Top 20 of pathway enrichment. The circle size represented gene number. Q-value was shown by color gradient.
Caption: Figure 7: miR-125a-5p downregulates the expression of ESRR[alpha] by targeting its 3' UTR site. (a) The pairing schematic of ESRR[alpha] 3' UTR and miR-125a-5p. The nucleotide in blue represents "seed sequence" of miR-125a-5p, and the mutation nucleotides are presented in red. (b) The insertion site of ESRR[alpha] 3' UTR or its mutation luciferase reporter vector map. (c) The two types were cotransfected with miR-125a-5p mimics (or negative control) into 293T cells. Luciferase assay was performed 48 h after transfection. Results are presented as relative luciferase activity (mean [+ or -] SD; n = 4; ANOVA; * P < 0.01).
Table 1: Summary of total small RNA tags by Solexa sequencing. Type KS Counts Percentage (%) total_reads 12,000,000 high_quality 11,977,968 100 3'adapter_null 54,706 0.46 insert_null 4,294 0.04 5'adapter_contaminants 5,406 0.05 smaller_than_18nt 66,241 0.55 polyA 4 0.00 clean_reads 11,847,317 98.91 Type TS Counts Percentage (%) total_reads 12,000,000 high_quality 11,973,801 100 3'adapter_null 61,670 0.52 insert_null 5,122 0.04 5'adapter_contaminants 5,610 0.05 smaller_than_18nt 45,410 0.38 polyA 3 0.00 clean_reads 11,855,986 99.02 Table 2: Highly expressed miRNAs (top 10) in sheep adipose tissues. In fat-tailed tissue (KS) miRNA P value miR-2070-3p 0 miR-222 2.86E-189 miR-502-3p 1.24R-179 miR-7446-3p 8.35B-122 miR-6724-5p 4.92B-104 miR-126 2.53E-101 miR-378e 0 miR-30c 9.26B-60 miR-193 3.71B-59 miR-2312 2.13B-44 In thin-tailed tissue (TS) miRNA P value miR-4792 3.55B-18 miR-106a 8.87B-19 miR-3649 4.44B-19 miR-4152-5p 2.17.E-22 miR-125a-3p 3.34B-302 miR-4492 5.80B-48 miR-993-5p 5.57B-54 miR-7930-3p 3.35B-61 miR-7475-5p 2.33R-115 miR-6238 5.41R-162 In both tissues miRNA P value miR-25-3p 2.07B-01 miR-25 1.09R-02 miR-10a-5p 2.43B- 18 miR-10 3.94B- 18 miR-10a 1.77E-23 miR-92a-3p 1.01R-26 miR-10-5p 5.49B-72 miR-92b-3p 4.72B-29 let-7 0 miR-486-5p 0 Table 3: Lists of miRNA and their targets with functions related to adipogenesis and/or fat metabolism. miRNA Target Partial target genes numbers miR-2070-3p 1352 SH3D21, BCL7C, ACTR3B, EPC1 miR-222 624 RGS6, HMG20A, RBM15, NFE2 miR-502-3p 462 CRTC1, FGD1, CCL8, STARD8 miR-6238 391 Mcph1, PDZK1 miR-7446-3p 414 KLF13, SIAH2, TUB miR-7475-5p 517 LDB1, DVL3, PEG3, LRP1, LATS2, EFHD2 miR-125a-5p 2246 ESRRa, SENP2, BCL2L12, SREBP-1, ABCA2, NNMT miR-126 2438 TNKS2, PTPRU, RGS14, NAP1L5 miR-378e 1044 IGF1R, CACNB2, RASIP1, API5, SCD5, SLC25A29 miR-7930-3p 1793 CABIN1, PCDHA2, PLXNA4
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Zhou, Guangxian; Wang, Xiaolong; Yuan, Chao; Kang, Danju; Xu, Xiaochun; Zhou, Jiping; Geng, Rongqing|
|Publication:||BioMed Research International|
|Date:||Jan 1, 2017|
|Previous Article:||Protective Effect of a Polyherbal Aqueous Extract Comprised of Nigella sativa (Seeds), Hemidesmus indicus (Roots), and Smilax glabra (Rhizome) on...|
|Next Article:||Corrigendum to "Instructors' Perceptions of Mostly Seated Exercise Classes: Exploring the Concept of Chair Based Exercise".|