http://orcid.org/0000-0001-5842-1703
Figures
Abstract
In this study, we examined the role of mRNAs and miRNAs in variations in intramuscular fat content in the longissimus dorsi muscle in Xinjiang brown cattle. Two groups of Xinjiang brown cattle with extremely different intramuscular fat content in the longissimus dorsi were selected for combined of miRNA and mRNA analysis using an RNA-Seq. In total, 296 mRNAs and 362 miRNAs were significantly differentially expressed, including 155 newly predicted miRNAs, 275 significantly upregulated genes, 252 significantly upregulated miRNAs, 21 significantly downregulated genes and 110 significantly downregulated miRNAs. The combined miRNA and mRNA analysis identified 96 differentially expressed miRNAs and 27 differentially expressed mRNAs. In all, 47 upregulated miRNAs had a regulatory effect on 14 differentially downregulated target genes, and 49 downregulated miRNAs had a regulatory effect on 13 upregulated target genes. To verify the sequencing results, 10 differentially expressed genes (DEGs) and 10 differentially expressed miRNAs were selected for qRT-PCR. The qRT-PCR results confirmed the sequencing results. The results of this study shed light on the molecular regulation of bovine adipose tissue, which might help with the development of new strategies for improving meat quality and animal productivity in beef cattle to provide healthier meat products for consumers.
Citation: Li N, Zhang Y, Li H-P, Han L, Yan X-M, Li H-B, et al. (2018) Differential expression of mRNA-miRNAs related to intramuscular fat content in the longissimus dorsi in Xinjiang brown cattle. PLoS ONE 13(11): e0206757. https://doi.org/10.1371/journal.pone.0206757
Editor: Yun Zheng, Kunming University of Science and Technology, CHINA
Received: June 15, 2018; Accepted: October 18, 2018; Published: November 9, 2018
Copyright: © 2018 Li et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All relevant data are within the paper and its Supporting Information files.
Funding: This work was supported by the National Beef Cattle Industrial Technology System (CARS-37), Key Research of Autonomous Region and Plan Projects (2017B01001) and the National Nature Science Foundation (31460647).
Competing interests: The authors have declared that no competing interests exist.
Introduction
Xinjiang brown cattle are an improved breed with an independent intellectual property right. Numbering over 1.7 million, Xinjiang brown cattle are the main cattle breed farmed in Xinjiang. However, no systemic breed selection has been performed since the breed was established [1], leading to large fluctuations in the quality of snow beef and greatly impairing the international competitiveness of Xinjiang brown cattle products. The molecular mechanisms underlying fat deposition in Xinjiang brown cattle have become a popular topic of research.
MicroRNAs (miRNAs) constitute a class of endogenous, small non-coding RNAs of approximately 21–25 nucleotides (nt) in length [2]. In mammals, miRNAs regulate signaling pathways related to fat production by inhibiting transcriptional factors involved in adipogenesis. miRNAs have been reported to negatively regulate fat production by acting on insulin signaling and C/EBPs [3]. In 3T3-L1 cells, miR-320 and miR-29 act directly on phosphatidyl inositol kinase (PIK), which is an important factor in insulin signaling, to reduce insulin activity and inhibit adipocyte formation [4]. miR-27a and miR-27b act on the prohibitin (PHB) gene to negatively regulate adipocyte formation by inhibiting PHB [5].
However, the identification of miRNA target genes remains a considerable challenge because the effective target sequences of miRNAs may not be conserved [6]. In addition, to regulate mRNAs, the bases in the seed sequence of miRNAs bind 6 complementary bases at the target site of the mRNA [7], which increases the complexity of the regulatory effect of miRNAs on mRNAs. The following strategy has been shown to be practically effective: first, small RNA (sRNA) sequencing is performed to identify and predict differentially expressed miRNA target genes; then, the differentially expressed mRNAs are screened; and finally, a correlation analysis between the target genes and differentially expressed genes (DEGs) is conducted to determine the regulatory relationship between the miRNAs and mRNAs. This method has been widely applied to identify miRNA target genes in different tissues [8, 9]. In the present study, sRNA sequencing, prediction of differentially expressed miRNA target genes and screening for differentially expressed mRNAs were performed via RNA-Seq. Then, a correlation analysis of the DEGs was performed to determine the key factors that regulate fat deposition in Xinjiang brown cattle. The study findings provide a theoretical basis for screening effective biomarkers of lipid metabolism.
Materials and methods
Animals
The cattle were provided by the Yili Tianxi Breeding Company. Thirty-three healthy castrated Xinjiang brown cattle aged 10–12 months were selected for a 22-month fattening experiment. All animal procedures were approved by the IACUC of Institute of Research Livertock, Xinjiang Academy of Animal Science (approval number: 20150501).
Slaughtering, segmentation and sampling
Cattle were slaughtered after being fasted for 12 h. The longissimus dorsi at the 12th and 13th ribs was harvested, labeled and preserved in liquid nitrogen. The fat content in the longissimus dorsi was measured. The two cattle with the lowest fat content in the longissimus dorsi were chosen as control group A, and the two cattle with the highest fat content were chosen as experimental group B. See Fig 1.
A: Control group A, with the lowest fat content in the longissimus dorsi. B: Group B, with the highest fat content in the longissimus dorsi.
mRNA library construction and sequencing
Total RNA was isolated using a mirVana miRNA Isolation Kit-AM1561 (Ambion, Austin, TX). An enriched cDNA library was constructed using the isolated RNA. Polyacrylamide gel electrophoresis (PAGE) was performed, and the RNA bands were then cut from the gel. Then, 18–30 nt RNA was obtained, and miRNAs were constructed. Finally, miRNA and mRNA sequencing was performed using a BGISEQ-50 sequencing system (BGI).
Data processing
Impurities were filtered to obtain clean reads. Gene quantification was performed using the RSEM program [10]. DEGs with an FDR ≤0.001 and a fold-change ≥2 were chosen using the Poisson process.
The clean reads obtained by the sRNA sequencing were aligned to reference sequences using HISAT [11] and Bowtie2 [12]. sRNA expression was normalized using TPM [13].
Prediction of miRNA target genes and correlation analysis between miRNAs and mRNAs
The miRNA target genes were predicted using RNAhybrid, miRanda and TargetScan [14, 15]. The differentially expressed mRNAs and miRNAs were used to construct DEG libraries. Then, a correlation analysis was performed between the DEGs in the two libraries, and differentially expressed target genes that were negatively correlated with miRNAs were selected. An integrated network analysis was conducted. Differentially expressed miRNAs and their target genes were chosen from the different fat-content groups.
GO and pathway enrichment analyses
GO and pathway enrichment analyses were conducted using ‘GO: TermFinder’ (http://www.yeastgenome.org/help/analyze/go-term-finder). The significantly enriched pathways (Q value ≤0.05) among the DEGs were identified by performing an alignment against the whole genome.
Construction of miRNA-target interaction network
GO and KEGG pathway analyses were conducted for the DEGs that were negatively correlated with the differentially expressed miRNAs. The genes enriched in both key pathways and GO terms were selected to construct the DE-miRNA network.
Quantitative validation of the differentially expressed mRNAs and miRNAs
Using GAPDH and U6 as internal reference genes, the relevant differentially expressed mRNAs and miRNAs and their downstream functions were chosen for qRT-PCR validation. The primers were synthesized by Sangon Biotech (Shanghai) Co., Ltd. (S1 Table). The qRT-PCR validation of the mRNAs was performed using ABI Power SYBR Green PCR Master Mix according to the manufacturer’s instructions. The qRT-PCR validation of the miRNAs was performed using a miRcute miRNA qPCR Detection Kit (SYBR Green) according to the manufacturer’s instructions. The qRT-PCR parameters were as follows: 40 cycles of incubation at 37°C for 60 min, 95°C for 10 min, 95°C for 15 s, and 60°C for 1 min. The melting curve was plotted.
Results
Analysis of the reads obtained from mRNA and miRNA sequencing
The samples were sequenced using RNA-Seq. This generated 24,136,913–24,137,316 raw reads. After removing low-quality reads, 24,090,787–24,103,791 clean reads were retained. After alignment, the clean reads accounted for 99.80% of the reads in Xinjiang brown cattle.
In total, 24,151,865–71,722,737 raw reads were obtained from 4 sRNA libraries. After filtering out impurities, 22,650,477–64,436,993 clean reads were retained, accounting for 93.84% of the total. The reads were aligned to the known sRNA database to identify sRNAs. Sequences with a length of 18–30 nt accounted for the largest proportion, and this satisfied the requirement for experimental sRNAs, as shown in Fig 2A.
A: The x-axis represents the length of the sRNA, and the y-axis represents the proportion of sRNAs with different lengths. The sequence length distribution reflects the types of sRNAs. B: Smaller correlations are closer to white; while larger correlations (closer to 1) are closer to blue.
A correlation analysis of the mRNAs among the four samples was then performed, as shown in Fig 2B. The correlation coefficients were all above 0.98, suggesting that the expression patterns were highly similar among the samples.
Differential mRNA and miRNA expression in longissimus dorsi with different fat contents
Hierarchical clustering of the DEGs was performed using Cluster software, and the results are shown in Fig 3.
A: Hierarchical clustering of the DEGs between the groups. B: Hierarchical clustering of the miRNAs. Each row represents a DEG. Different types of DEGs are represented by different colors. Red, upregulated DEGs; blue, downregulated DEGs.
In total, 296 DEGs, including 275 upregulated and 21 downregulated genes, were identified. In all, 362 differentially expressed miRNAs, including 155 newly predicted miRNAs, 252 upregulated genes and 110 downregulated genes, were revealed. The target genes of the differentially expressed miRNAs were predicted, and 209,960 target genes, including 108,256 novel miRNAs, were involved.
In the GO enrichment analysis of the DEGs, 214 of the 296 DEGs were annotated. In all, 712 enriched biological process GO terms, 110 enriched cellular component GO terms, and 111 enriched molecular function GO terms were revealed. Pathway enrichment analysis identified 39 significant pathways, 8 of which were related to lipid metabolism. The enriched lipid metabolism pathways included the biosynthesis of unsaturated fatty acids, fatty acid elongation and fatty acid biosynthesis. After removal of repetitive genes, 11 gene clusters remained. (S2 and S3 Tables).
GO and pathway enrichment analyses of the target genes of the differentially expressed sRNAs
GO functional analysis was performed on the 22,590 target genes of miRNAs that were differentially expressed between the two groups. In all, 401 cellular component GO terms, 783 molecular function GO terms, and 3,231 biological process GO terms were identified. Pathway enrichment analysis identified 44 significant pathways. The predicted target genes were involved in metabolic pathways, regulation of the actin cytoskeleton, and the MAPK signaling pathway. There were 14 significant pathways. The enriched lipid metabolism pathways included fatty acid degradation, ether lipid metabolism, and glycerophospholipid metabolism. After removal of repetitive target genes, 460 clusters remained, as shown in S3 Table.
Correlation analysis of the target genes of the differentially expressed miRNAs and differentially expressed mRNAs
DEG and target gene libraries were constructed considering that miRNAs post-transcriptionally regulate, degrade or inhibit translation [16]. A correlation analysis of the genes from the two libraries was conducted. DEGs that were negatively correlated with the miRNAs and were related to lipid metabolism were chosen, and a regulatory network was constructed. In all, 96 differentially expressed miRNAs in groups A and B were correlated to 27 target genes. There were 47 upregulated miRNAs and 49 downregulated miRNAs. The 47 upregulated miRNAs had 14 downregulated target genes, and the 49 downregulated miRNAs had 13 upregulated target genes. According to the miRNA-target network (Fig 4), the miRNAs novel-mir2, novel-mir240, miR-130a, miR-125a, miR-181a and miR-30b-5p were located at the center.
A: Upregulated miRNAs and their downregulated target genes. B: Downregulated miRNAs and their upregulated target genes. Red lines represent the intensity of the interaction between the upregulated miRNAs and their downregulated target genes; green lines represent the intensity of the interaction between the downregulated miRNAs and their upregulated target genes.
GO and pathway enrichment analyses of the differentially expressed target genes of the miRNAs
GO enrichment analysis revealed that the differentially targeted genes belonged to 71 biological processes, including fatty acid metabolic process and neutral lipid metabolic process; 41 cellular components, including extracellular region part, extracellular region, and cytoplasmic region; and 33 molecular functions, including fatty acid synthase activity, lipase inhibitory activity, and lipid binding. Pathway enrichment analysis showed that the upregulated DEGs belonged to 16 pathways, including biosynthesis of unsaturated fatty acids, fatty acid metabolism, and fatty acid elongation, and the downregulated DEGs belonged to 8 pathways, including the MAPK signaling pathway, glutamatergic synapses, and the glucagon signaling pathway (S4 Table).
qRT-PCR validation of the differentially expressed miRNAs and mRNAs
Based on the combined analysis of the miRNAs and mRNAs, we selected 10 DEGs for fluorescent quantitative PCR. FABP4, LEP, KLF6, AEBP1, CXCL9, CAV1, ANGPTL8, DGAT2, SCD and MAP2K6 were associated with adipogenesis and adipocyte differentiation. Let-7i-5p, miR-125a, miR-199a-5p, miR-200b, miR-143, miR-122, miR-375, miR-204-5p, miR-181a and miR-429 affect adipocyte differentiation. The qRT-PCR validation confirmed that the differentially expressed miRNAs and mRNAs obtained by sequencing were accurate, as shown in Fig 5.
A: Comparison of the mRNA sequencing and qRT-PCR results. B: Comparison of the miRNA sequencing and qRT-PCR results.
Discussion
Meat quality is determined on the basis of marbling, texture, fat color, meat color and overall maturity [17]. Fat deposition is positively related to beef flavor and palatability. The deposition of intramuscular adipose tissue influences both meat quality and animal productivity [18]. Understanding the molecular regulation of adipose tissue is a key step toward improving beef quality. Bovine adipose tissue may be regulated by upstream mechanisms such as miRNAs [19] and transcriptional factors [20]. However, few studies have investigated the differential expression of miRNAs related to intramuscular fat content in the longissimus dorsi in Xinjiang brown cattle. Garin et al. (2014) found that FABP4, also known as aP2, suppresses PPARγ and adipogenesis [21]. Bojlul et al. (2013) showed that leptin plays an important role in energy homeostasis [22]. Lee et al. (2012) found that KLF8 is an upstream regulator of PPARγ and C/EBPα and is a key component during adipogenesis [23]. Boqué et al. (2013) revealed that apple polyphenol (AP) intake led to decreased adipocyte enhancer-binding protein (AEBP1) mRNA levels in epididymal adipocytes [24]. Gerst et al. (2017) showed that palmitate and fetuin-A induced CXCL8 mRNA expression in differentiated adipocytes [25]. Sara et al. (2015) found that Caveolin-1 (Cav-1) is a key component of adipocyte caveolae and plays an important role in regulating insulin signaling [26]. Bruno et al. (2016) revealed that Angiopoietin-like protein 8 (ANGPTL8) is implicated in lipid and glucose homeostasis [27]. Ryo et al. (2005) showed that diacylglycerol acyltransferase (DGAT) enzymes facilitated the final step in mammalian triglyceride synthesis and that DGAT was associated with obesity, insulin resistance, and leptin resistance [28]. Niclas et al. (2011) found that Stearoyl-CoA desaturase-1 (SCD-1) is a key enzyme in lipogenesis and that inhibiting its activity impairs lipid synthesis [29]. Jihye et al. (2012) revealed that MAP2K6 affects carcass weight (CW), back fat thickness (BFT) and marbling score (MS) [30].
Giroud et al. (2016) showed that Let-7i-5p was downregulated and repressed brite adipocyte function in human brite adipocytes [31]. Ji et al. (2014) revealed that miR-125a suppresses porcine preadipocytes differentiation by downregulating ERRa [32]. Shi et al. (2014) found that miR-199a-5p facilitated cell proliferation and inhibited lipid deposition in porcine adipocytes by directly targeting Cav-1 [33]. Shen et al. (2018) reported that miR-200b inhibited preadipocyte proliferation by attenuating KLF4 [34]. Kim et al. (2015) showed that miR-143 facilitated adipogenesis by directly targeting pref-1 in adipocytes [35]. Liao et al. (2018) revealed that miR-122 regulated PPAR-γ signaling and adipocyte differentiation in vitro and in human adipose tissue [36]. Liu et al. (2016) found that miR-375 inhibited porcine preadipocyte differentiation by directly targeting BMPR2 [37]. Du et al. (2018) reported that miR-204-5p promoted adipocyte differentiation by directly targeting KLF3 [38]. Li et al. (2013) found that miR-181a accelerated lipid droplet accumulation and increased triglyceride levels by negatively regulating TNF-α [39]. Jiang et al. (2018) showed that miR-499 represses adipogenic differentiation by directly targeting PRDM16 in skeletal muscle satellite cells (SMSCs) [40].
Along with the increased depth provided by next-generation sequencing, RNA-Seq has become a powerful tool for performing comprehensive transcriptome analyses [41]. Compared with microarray and SAGE technologies, RNA-Seq is more sensitive and more widely used [42]. The massive data obtained from RNA-Seq can be analyzed using Bowtie2 and TargetScan software [15]. Here, we applied RNA-Seq to sequence the longissimus dorsi in Xinjiang brown cattle, and the data were analyzed using Bowtie2 software. The lean reads obtained from the transcriptome library accounted for over 95% of the original reads. Moreover, the clean reads obtained from the sRNA library accounted for over 93%, indicating that most known miRNAs are expressed in the muscles of Xinjiang brown cattle. However, not all miRNAs were detected, indicating the tissue specificity of the miRNAs [43]. From another perspective, more tissue-specific miRNAs will likely be discovered in cattle in the future.
Establishing biological replicates is required in any biological experiment, and high-throughput sequencing is no exception [44]. Correlating gene expression between different samples is an important indicator of experimental reliability and the reasonability of sample selection. The closer the correlation coefficient is to 1, the higher the similarity between the two samples. We used only two biological replicates in this study, and this may have led to false-positive results. According to FPKM, the correlation coefficients between the two groups were all above 0.985, which satisfied the experimental requirement. Furthermore, 10 miRNAs and 10 miRNAs were selected from 6 samples for qRT-PCR validation, and the qRT-PCR results were consistent with the sequencing results, supporting the reliability of the DEGs obtained by RNA-Seq. Thus, using only one or two biological replicates was sufficient to identify key DEGs.
miRNAs are widely present in the genomes of animals and plants and play regulatory roles in lipid metabolism in cattle. Jin detected miRNAs in adipose tissues on the backs of three types of crossbred cattle with different thicknesses of back fat. Eighty-nine miRNAs were detected, including 42 differentially expressed miRNAs, with the expression of miR-378 being the most significantly different [45]. This finding was confirmed in our study. Li observed that miR-143 was upregulated in intramuscular adipocytes that differentiated into mature adipocytes in cattle [46]. Once the expression of miR-143 was inhibited, the differentiation of the precursors of the intramuscular adipocytes was also inhibited. In the present study, miR-143 expression was significantly higher in the experimental group than in the control group, suggesting that miR-143 plays a potential role in the formation of intramuscular fat in cattle.
To further explore the potential regulatory role of the aforementioned miRNAs, the target genes of the differentially expressed miRNAs were predicted. In all, 209,960 target genes were involved, and the miRNA-target interaction network indicated that one miRNA may regulate several dozens or hundreds of target genes. In addition, each target gene was regulated by several miRNAs. In this study, 1,074 target genes were predicted for miR-143, which exhibited the largest fold-change difference between the two groups.
The post-transcriptional regulation of target genes by miRNAs has been proved in numerous studies [47–49]. Here, RNA-Seq was innovatively applied to perform a combined analysis of the target genes of differentially expressed miRNAs and differentially expressed mRNAs. Ninety-six negatively correlated miRNAs and 27 target genes were identified. These 27 target genes are key genes responsible for the variability in the fat content and meat quality traits of the longissimus dorsi in Xinjiang brown cattle. The related pathways were the focus of this study. These differentially expressed target genes belonged to 71 biological processes, 41 cellular components and 33 molecular functions. Pathway enrichment analysis indicated that the upregulated DEGs belonged to 16 pathways, including the biosynthesis of unsaturated fatty acids, fatty acid metabolism, the pentose phosphate pathway and fatty acid elongation, while the downregulated DEGs belonged to 8 pathways, including the MAPK signaling pathway, glutamatergic synapses and hypertrophic cardiomyopathy (HCM). The MAPK signaling pathway, which was the most significantly enriched pathway, is involved in the conversion of extracellular signals into intracellular signals in eukaryotes. Studies have shown that miRNAs can target certain genes in the MAPK signaling pathway, thus influencing signal transduction and participating in adipocyte metabolism [50–52]. Lin et al. found that PHB was a target gene of miR-27a and miR-27b [53]. PHB is highly expressed in adipocytes and is closely associated with adipocyte differentiation. PHB overexpression in 3T3-L1 cells inhibits insulin-induced adipogenic differentiation. However, many pathways may influence adipocyte differentiation, and MAPK signaling is but one such key pathway. Moreover, the pRB-E2F, Wnt, biosynthesis of unsaturated fatty acids and fatty acid metabolism pathways identified in this study also play a role in adipocyte metabolism. These pathways, together with miRNAs, form a large network that regulates cell differentiation.
Conclusion
We applied RNA-Seq to the longissimus dorsi of Xinjiang brown cattle and identified 24,098,031 clean reads. The average alignment rate of the clean reads against the reference genome was 95.35%. sRNA sequencing produced 46,007,468 clean reads, which accounted for 93.33% of the total reads. Thus, the sequencing had a high yield and effectively detected mRNAs and miRNAs in the tissues.
Four mRNA and sRNA libraries were constructed for the longissimus dorsi of Xinjiang brown cattle, and 296 differentially expressed mRNAs and 362 differentially expressed miRNAs were identified as follows: 155 newly predicted miRNAs, 275 upregulated genes, 252 upregulated miRNAs, 21 downregulated genes and 110 downregulated miRNAs. The combined analysis of the miRNAs and mRNAs revealed 96 differentially expressed miRNAs and 27 differentially expressed mRNAs. We found that 47 upregulated miRNAs had a regulatory effect on 14 downregulated target genes and that 49 downregulated miRNAs had a regulatory effect on 13 upregulated target genes. An miRNA-target interaction network was constructed. Of the DEGs, those in metabolic pathways, especially pathways related to lipid metabolism, were enriched by the largest proportions.
Supporting information
S1 Table. Primers for reverse transcription and quantitative real-time PCR.
https://doi.org/10.1371/journal.pone.0206757.s001
(DOCX)
S2 Table. Differential expression of miRNAs in longissimus dorsi with different fat content.
https://doi.org/10.1371/journal.pone.0206757.s002
(XLS)
S3 Table. Differential expression of mRNAs in longissimus dorsi with different fat content.
https://doi.org/10.1371/journal.pone.0206757.s003
(XLS)
S4 Table. Differentially expressed target genes of miRNAs.
https://doi.org/10.1371/journal.pone.0206757.s004
(XLS)
Acknowledgments
This work was supported by the National Beef Cattle Industrial Technology System (CARS-37), Key Research of Autonomous Region and Plan Projects (2017B01001) and the National Nature Science Foundation (31460647).
References
- 1. Li N, Li HB, Yan XM, Du W, Zhang JS, Yuan LX, et al. Correlation between four metabolism-related genes in different adipose tissues and adipocyte morphology in Xinjiang brown cattle. Int J Clin Exp Med. 2016;9(3):5912–21.
- 2. Bartel DP. MicroRNAs: Genomics, biogenesis, mechanism, and function (Reprinted from Cell, vol 116, pg 281–297, 2004). Cell. 2007;131(4):11–29.
- 3. Ling HY, Wen GB, Feng SD, Tuo QH, Ou HS, Yao CH, et al. MicroRNA-375 promotes 3T3-L1 adipocyte differentiation through modulation of extracellular signal-regulated kinase signalling. Clin Exp Pharmacol P. 2011;38(4):239–46. pmid:21291493
- 4. He AB, Zhu LB, Gupta N, Chang YS, Fang F. Overexpression of micro ribonucleic acid 29, highly up-regulated in diabetic rats, leads to insulin resistance in 3T3-L1 adipocytes. Mol Endocrinol. 2007;21(11):2785–94. pmid:17652184
- 5. Karbiener M, Fischer C, Nowitsch S, Opriessnig P, Papak C, Ailhaud G, et al. microRNA miR-27b impairs human adipocyte differentiation and targets PPAR gamma. Biochem Bioph Res Co. 2009;390(2):247–51. pmid:19800867
- 6. Tay Y, Zhang JQ, Thomson AM, Lim B, Rigoutsos I. MicroRNAs to Nanog, Oct4 and Sox2 coding regions modulate embryonic stem cell differentiation. Nature. 2008;455(7216):1124–U12. pmid:18806776
- 7. Lim LP, Lau NC, Garrett-Engele P, Grimson A, Schelter JM, Castle J, et al. Microarray analysis shows that some microRNAs downregulate large numbers of target mRNAs. Nature. 2005;433(7027):769–73. pmid:15685193
- 8. Dmitriev P, Barat A, Polesskaya A, O’Connell MJ, Robert T, Dessen P, et al. Simultaneous miRNA and mRNA transcriptome profiling of human myoblasts reveals a novel set of myogenic differentiation-associated miRNAs and their target genes. Bmc Genomics. 2013;14.
- 9. Miao XY, Qin QLX. Genome-wide transcriptome analysis of mRNAs and microRNAs in Dorset and Small Tail Han sheep to explore the regulation of fecundity. Mol Cell Endocrinol. 2015;402(C):32–42. pmid:25573241
- 10. Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. Bmc Bioinformatics. 2011;12. pmid:21816040
- 11. Kim D, Landmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–U121. pmid:25751142
- 12. Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3). pmid:19261174
- 13. t Hoen PAC, Ariyurek Y, Thygesen HH, Vreugdenhil E, Vossen RHAM, de Menezes RX, et al. Deep sequencing-based expression analysis shows major advances in robustness, resolution and inter-lab portability over five microarray platforms. Nucleic Acids Res. 2008;36(21). pmid:18927111
- 14. Kruger J, Rehmsmeier M. RNAhybrid: microRNA target prediction easy, fast and flexible. Nucleic Acids Res. 2006;34:W451–W4. pmid:16845047
- 15. Agarwal V, Bell GW, Nam JW, Bartel DP. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4. pmid:26267216
- 16. Healy NA, Creighton CJ, Fu X, Tsimelzon A, Hilsenbeck SG, Miller N, et al. Identification of miRNAs and Their Associated Target Genes Involved in Endocrine Resistance in Breast Cancer. Cancer Research. 2011;71.
- 17. Fan YY, Fu GW, Fu CZ, Zan LS, Tian WQ. A missense mutant of the PPAR-gamma gene associated with carcass and meat quality traits in Chinese cattle breeds. Genet Mol Res. 2012;11(4):3781–8. pmid:22930427.
- 18. Romao J. M. H ML, McAllister T. A., and Guan L. L. Effect of age on bovine subcutaneous fat proteome: Molecular mechanisms of physiological variations during beef cattle growth. American Society of Animal Science. 2014.
- 19. Shen L, Chen L, Zhang S, Zhang Y, Wang J, Zhu L. MicroRNA-23a reduces slow myosin heavy chain isoforms composition through myocyte enhancer factor 2C (MEF2C) and potentially influences meat quality. Meat Sci. 2016;116:201–6. pmid:26897085.
- 20. Grindflek E, Hoen N, Sundvold H, Rothschild MF, Plastow G, Lien S. Investigation of a peroxisome proliferator-activated receptor gamma haplotype effect on meat quality and carcass traits in pigs. Anim Genet. 2004;35(3):238–41. pmid:15147398.
- 21. Garin-Shkolnik T, Rudich A, Hotamisligil GS, Rubinstein M. FABP4 attenuates PPARgamma and adipogenesis and is inversely correlated with PPARgamma in adipose tissues. Diabetes. 2014;63(3):900–11. pmid:24319114.
- 22. Bahar B, O’Doherty JV, O’Doherty AM, Sweeney T. Chito-oligosaccharide inhibits the de-methylation of a ’CpG’ island within the leptin (LEP) promoter during adipogenesis of 3T3-L1 cells. PLoS One. 2013;8(3):e60011. pmid:23544120.
- 23. Lee H, Kim HJ, Lee YJ, Lee MY, Choi H, Lee H, et al. Kruppel-like factor KLF8 plays a critical role in adipocyte differentiation. PLoS One. 2012;7(12):e52474. pmid:23285057.
- 24. Boque N, de la Iglesia R, de la Garza AL, Milagro FI, Olivares M, Banuelos O, et al. Prevention of diet-induced obesity by apple polyphenols in Wistar rats through regulation of adipocyte gene expression and DNA methylation patterns. Mol Nutr Food Res. 2013;57(8):1473–8. pmid:23529981.
- 25. Gerst F, Wagner R, Kaiser G, Panse M, Heni M, Machann J, et al. Metabolic crosstalk between fatty pancreas and fatty liver: effects on local inflammation and insulin secretion. Diabetologia. 2017;60(11):2240–51. pmid:28791439.
- 26. Palacios-Ortega S, Varela-Guruceaga M, Algarabel M, Ignacio Milagro F, Alfredo Martinez J, de Miguel C. Effect of TNF-Alpha on Caveolin-1 Expression and Insulin Signaling During Adipocyte Differentiation and in Mature Adipocytes. Cell Physiol Biochem. 2015;36(4):1499–516. pmid:26159107.
- 27. Martinez-Perez B, Ejarque M, Gutierrez C, Nunez-Roa C, Roche K, Vila-Bedmar R, et al. Angiopoietin-like protein 8 (ANGPTL8) in pregnancy: a brown adipose tissue-derived endocrine factor with a potential role in fetal growth. Transl Res. 2016;178:1–12. pmid:27469268.
- 28. Suzuki R, Tobe K, Aoyama M, Sakamoto K, Ohsugi M, Kamei N, et al. Expression of DGAT2 in white adipose tissue is regulated by central leptin action. J Biol Chem. 2005;280(5):3331–7. pmid:15550388.
- 29. Franck N, Gummesson A, Jernas M, Glad C, Svensson PA, Guillot G, et al. Identification of adipocyte genes regulated by caloric intake. J Clin Endocrinol Metab. 2011;96(2):E413–8. pmid:21047925.
- 30. Lee C R J K Y K C K J. Association of bovine carcass phenotypes with genes in an adaptive thermogenesis pathway. Mol Biol Rep. 2012;39(1441–1445). pmid:21617946
- 31. Giroud M, Karbiener M, Pisani DF, Ghandour RA, Beranger GE, Niemi T, et al. Let-7i-5p represses brite adipocyte function in mice and humans. Sci Rep. 2016;6:28613. pmid:27345691.
- 32. Yang J H-l S C-C, L Y-F, H J-j, L Y-l, Z X-l G-s . miR-125a inhibits porcine preadipocytes differentiation by targeting ERRa. Mol Cell Biochem. 2014;395:155–65. pmid:24952481
- 33. Shi XE, Li YF, Jia L, Ji HL, Song ZY, Cheng J, et al. MicroRNA-199a-5p affects porcine preadipocyte proliferation and differentiation. Int J Mol Sci. 2014;15(5):8526–38. pmid:24830555.
- 34. Shen L, Gan M, Li Q, Wang J, Li X, Zhang S, et al. MicroRNA-200b regulates preadipocyte proliferation and differentiation by targeting KLF4. Biomed Pharmacother. 2018;103:1538–44. pmid:29864940.
- 35. Kim Y-J M K TS, S K-S, SH . Expression of pref-1/dlk-1 is regulated by microRNA-143 in 3T3-L1 cells. Mol Biol Rep. 2015;46:617–24. pmid:25366176
- 36. Liao CH, Wang CY, Liu KH, Liu YY, Wen MS, Yeh TS. MiR-122 marks the differences between subcutaneous and visceral adipose tissues and associates with the outcome of bariatric surgery. Obes Res Clin Pract. 2018. pmid:29960868.
- 37. Liu S, Sun G, Yuan B, Zhang L, Gao Y, Jiang H, et al. miR-375 negatively regulates porcine preadipocyte differentiation by targeting BMPR2. FEBS Lett. 2016;590(10):1417–27. Epub 2016/04/10. pmid:27059117.
- 38. Du J, Zhang P, Gan M, Zhao X, Xu Y, Li Q, et al. MicroRNA-204-5p regulates 3T3-L1 preadipocyte proliferation, apoptosis and differentiation. Gene. 2018;668:1–7. pmid:29775748.
- 39. Li H, Chen X, Guan L, Qi Q, Shu G, Jiang Q, et al. MiRNA-181a regulates adipogenesis by targeting tumor necrosis factor-alpha (TNF-alpha) in the porcine model. PLoS One. 2013;8(10):e71568. pmid:24098322.
- 40. Jiang J, Li P, Ling H, Xu Z, Yi B, Zhu S. MiR-499/PRDM16 axis modulates the adipogenic differentiation of mouse skeletal muscle satellite cells. Hum Cell. 2018. pmid:30097922.
- 41. Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8. pmid:18516045
- 42. Liu CG, Calin GA, Meloon B, Gamliel N, Sevignani C, Ferracin M, et al. An oligonucleotide microchip for genome-wide microRNA profiling in human and mouse tissues. P Natl Acad Sci USA. 2004;101(26):9740–4. pmid:15210942
- 43. Aboobaker AA, Tomancak P, Patel N, Rubin GM, Lai EC. Drosophila microRNAs exhibit diverse spatial expression patterns during embryonic development. P Natl Acad Sci USA. 2005;102(50):18017–22. pmid:16330759
- 44. Ye F, Samuels DC, Clark T, Guo Y. High-throughput sequencing in mitochondrial DNA research. Mitochondrion. 2014;17:157–63. pmid:24859348
- 45. Jin WW, Dodson MV, Moore SS, Basarab JA, Guan LL. Characterization of microRNA expression in bovine adipose tissues: a potential regulatory mechanism of subcutaneous adipose tissue development. Bmc Mol Biol. 2010;11. pmid:20423511
- 46. Li HX, Zhang Z, Zhou X, Wang ZY, Wang GL, Han ZY. Effects of MicroRNA-143 in the differentiation and proliferation of bovine intramuscular preadipocytes. Mol Biol Rep. 2011;38(7):4273–80. pmid:21113671
- 47. Cheng DF, Lu YY, Zeng L, Liang GW, He XF. Si-CSP9 regulates the integument and moulting process of larvae in the red imported fire ant, Solenopsis invicta. Sci Rep-Uk. 2015;5. pmid:25784646
- 48. Crea F, Watahiki A, Quagliata L, Xue H, Pikor L, Parolia A, et al. Identification of a long non-coding RNA as a novel biomarker and potential therapeutic target for metastatic prostate cancer. Oncotarget. 2014;5(3):764–74. pmid:24519926
- 49. Han DX, Xiao Y, Wang CJ, Jiang H, Gao Y, Yuan B, et al. Regulation of FSH expression by differentially expressed miR-186-5p in rat anterior adenohypophyseal cells. Plos One. 2018;13(3). pmid:29534107
- 50. Han DX, Sun XL, Xu MQ, Chen CZ, Jiang H, Gao Y, et al. Roles of differential expression of microRNA-21-3p and microRNA-433 in FSH regulation in rat anterior pituitary cells. Oncotarget. 2017;8(22):36553–65. Epub 2017/04/13. pmid:28402262.
- 51. Yang S, Guo LJ, Su YY, Wen J, Du J, Li XY, et al. Nitric oxide balances osteoblast and adipocyte lineage differentiation via the JNK/MAPK signaling pathway in periodontal ligament stem cells. Stem Cell Res Ther. 2018;9.
- 52. Kim KA, Kim JH, Wang Y, Sul HS. Pref-1 (Preadipocyte factor 1) activates the MEK/extracellular signal-regulated kinase pathway to inhibit adipocyte differentiation. Mol Cell Biol. 2007;27(6):2294–308. pmid:17210639
- 53. Lin XZ, Luo J, Zhang LP, Wang W, Shi HB, Zhu JJ. miR-27a suppresses triglyceride accumulation and affects gene mRNA expression associated with fat metabolism in dairy goat mammary gland epithelial cells. Gene. 2013;521(1):15–23. pmid:23537996