Next Article in Journal
Glyphosate Herbicide Induces Changes in the Growth Pattern and Somatic Indices of Crossbred Red Tilapia (O. niloticus × O. mossambicus)
Next Article in Special Issue
Expression Characteristics of microRNA in Pig Umbilical Venous Blood and Umbilical Arterial Blood
Previous Article in Journal
Modelling the Spatial Distribution of ASF-Positive Wild Boar Carcasses in South Korea Using 2019–2020 National Surveillance Data
Previous Article in Special Issue
Long Non-Coding RNAs in Insects
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Animals
Volume 11
Issue 5
10.3390/ani11051207
animals-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

RNA-Seq Reveals Function of Bta-miR-149-5p in the Regulation of Bovine Adipocyte Differentiation

by
Hongfang Guo
1,†,
Rajwali Khan
2,3,*,†,
Sayed Haidar Abbas Raza
3,
Syed Muhammad Suhail
2,
Hamayun Khan
4,
Sher Bahadar Khan
4,
Ayman Hassan Abd El-Aziz
5 and
Linsen Zan
3,6,*
1
Medical College, Xuchang University, Xuchang, Henan 461000, China
2
Livestock Management, Breeding and Genetics, The University of Agriculture, Peshawar 25120, Khyber Pakhtunkhwa, Pakistan
3
College of Animal Science and Technology, Northwest A&F University, Yangling 712100, China
4
College of Veterinary Sciences, The University of Agriculture, Peshawar 25120, Khyber Pakhtunkhwa, Pakistan
5
Department of Animal Husbandry and Animal Wealth Development, Faculty of Veterinary Medicine, Damanhour University, Damanhour 22511, Egypt
6
National Beef Cattle Improvement Research Center, Yangling 712100, China
*
Authors to whom correspondence should be addressed.
†
These authors contributed equally to this manuscript.
Animals 2021, 11(5), 1207; https://doi.org/10.3390/ani11051207
Submission received: 18 February 2021 / Revised: 16 April 2021 / Accepted: 18 April 2021 / Published: 22 April 2021
(This article belongs to the Special Issue Non-Coding RNAs Roles in Animals)
Browse Figures
Versions Notes

Abstract

:

Simple Summary

This is a novel study that explored the function of bta-miR-149-5p in the lipid metabolism of bovine adipocytes. Previously, we explored the role of bta-miR-149-5p gene in proliferation and differentiation of bovine adipocytes; however, the underlying mechanism of bta-miR-149-5p function in the regulation of lipid metabolism and adipogenesis in bovine adipocytes is unexplored. A transcriptomic study with integration of bioinformatics study was performed to fully explore the function of bta-miR-149-5p microRNA in bovine adipogenesis. This study will help the scientific community deal with lipogenesis for the breed improvement program of cattle for the provision of healthy meat to consumers.

Abstract

Intramuscular fat is a real challenge for the experts of animal science to improve meat quality traits. Research on the mechanism of adipogenesis provides invaluable information for the improvement of meat quality traits. This study investigated the effect of bta-miR-149-5p and its underlying mechanism on lipid metabolism in bovine adipocytes. Bovine adipocytes were differentiated and transfected with bta-miR-149-5p mimics or its negative control (NC). A total of 115 DEGs including 72 upregulated and 43 downregulated genes were identified in bovine adipocytes. The unigenes and GO term biological processes were the most annotated unigene contributor parts at 80.08%, followed by cellular component at 13.4% and molecular function at 6.7%. The KEGG pathways regulated by the DEGs were PI3K-Akt signaling pathway, calcium signaling pathway, pathways in cancer, MAPK signaling pathway, lipid metabolism/metabolic pathway, PPAR signaling pathway, AMPK signaling pathway, TGF-beta signaling pathway, cAMP signaling pathway, cholesterol metabolism, Wnt signaling pathway, and FoxO signaling pathway. In addition to this, the most important reactome enrichment pathways were R−BTA−373813 receptor CXCR2 binding ligands CXCL1 to 7, R−BTA−373791 receptor CXCR1 binding CXCL6 and CXCL8 ligands, R−BTA−210991 basigin interactions, R−BTA−380108 chemokine receptors binding chemokines, R−BTA−445704 calcium binding caldesmon, and R−BTA−5669034 TNFs binding their physiological receptors. Furthermore, the expression trend of the DEGs in these pathways were also exploited. Moreover, the bta-miR-149-5p significantly (p < 0.01) downregulated the mRNA levels of adipogenic marker genes such as CCND2, KLF6, ACSL1, Cdk2, SCD, SIK2, and ZEB1 in bovine adipocytes. In conclusion, our results suggest that bta-miR-149-5p regulates lipid metabolism in bovine adipocytes. The results of this study provide a basis for studying the function and molecular mechanism of the bta-miR-149-5p in regulating bovine adipogenesis.

1. Introduction

Adipose tissue stores energy in the form of neutral triglycerides and performs a vital function in maintaining energy homeostasis [1]. Adipocytes are the primary cell types of adipose tissue. There are two ways of adipose tissue development: one is hypertrophy (increase in adipocyte cell mass) and other is hyperplasia (increase in adipocyte numbers). In the case of surplus energy, the adipocytes deposits triglycerides (lipid droplets) through the lipogenic pathway, and in the case of energy deficiency, it releases lipids through the lipolysis pathway [2,3]. Mammalian body contains four important fat depots, namely, visceral, subcutaneous, intermuscular, and intramuscular (IM) fat. In livestock species, the IM fat is considered one of the most important factors that determines carcass quality traits. However, selective improvement of IM fat without altering the other three adipose depots is very challenging in the meat industry [4]. The IM fat deposition contributes to the development of meat sensory qualities such as marbling, juiciness, flavor, and tenderness [5]. Hence, approaches for the improvement of IM fat deposition are crucial for the development of meat quality. Therefore, research needs to be conducted to explore the underlying mechanisms of adipogenesis for the improvement of meat quality [6].
miRNAs are a class of endogenous small non-coding, approximately 22 nt, single-stranded RNAs, and play a role in post-transcriptional regulation by targeting mRNA degradation or inhibiting their translation [7]. A single miRNA can target hundreds of potential genes that can integrate complex molecular regulatory networks [8,9]. miR-143 was the first identified, with a proven role in adipocyte differentiation by regulating their target gene ERK5 (extracellular signal-regulated kinase 5) [10]. Previously, a large number of miRNAs were identified with their proven roles in regulation of adipose proliferation and differentiation. For instance, miR-2400 promoted bovine preadipocyte proliferation by directly targeting PR/SET domain 11 (PRDM11) [11]. The regulatory roles of miR-145 were determined in bovine adipogenesis through PI3K/Akt and MAPK signaling pathways [12]. Moreover, miR-144-3p suppresses preadipocyte proliferation [13]. The role of miR-424 in bovine adipogenesis through direct regulation of serine/threonine kinase 11 (STK11 also known as LKB1) was also identified through regulation of AMP-activated protein kinase (AMPK) [14]. Similarly, miR-17-92, miR-27, miR-93, miR-130, and miR-378 are key regulators of adipogenesis through adipocyte proliferation or differentiation [15,16,17,18]. Previously, our group exploited the negative regulatory role of miR-23a–miR27a–24-2 cluster in bovine adipocyte differentiation [19], and bta-miR-130a-b was found to affect bovine adipocyte differentiation by targeting PPARγ and CYP2UI [20]. Recently, studies have shown that miR-149 plays a pivotal role in the pathogenesis of various types of carcinomas including cancers of the digestive system [21], colorectal cancer [22], hepatocellular cancer [23,24], cancer of the respiratory system [25], breast cancer [26], and sebaceous gland carcinoma of the eyelid [27]. Therefore, this microRNA may be used as a potential diagnostic marker and therapeutic target against cancers. The miR-149-3p regulates energy homeostasis through the PRDM16 gene in visceral adipocytes [28]. Remarkably, miR-149 suppresses osteosarcoma through inhibiting TWEAK-Fn14 axis and is used as a potential therapeutic agent in osteosarcoma patients [29]. Interestingly, miR-149 enhances mitochondrial biogenesis and its activity through PGC1α activation by suppressing PARP-2 poly (ADP-ribose) polymerase-2 gene and enhancing cellular NAD+ levels through activation of the SIRT-1 gene. In obesity, miR-149 regulates mitochondrial biogenesis and function in SIRT-1/PGC-1α pathway through PARP-2 gene, and also regulates proliferation and differentiation of follicular B cells [30]. These findings proved the regulatory role of miR-149 in type-2 diabetes, high-fed diet-induced obesity, cancers, and cardiovascular diseases [31]. Previously, we proved the role of bta-miR-149-5p in bovine adipocyte proliferation and differentiation through regulation of TORCs [32]. It is intriguing that the regulatory function of bta-miR-149-5p will expand our understanding about the molecular mechanisms of bovine adipogenesis, which can be applied to enhance meat quality grading.

2. Materials and Methods

2.1. Ethical Statement

The experimental animals were dealt with as per standard operating procedures (SOPs) formulated by the Chinese Council of Animals Care, and further approved by “Experimental Animal Management Committee (EAMC)” of the Northwest Agricultural and Forestry University, notified vide notification no. EAMC/20-23 dated 20.04.2013. The sampling was performed after human killing of the experimental animals at the National Beef Cattle Improvement Research Center, Yangling, Shaanxi, China.

2.2. Isolation of Bovine Primary Preadipocytes

Bovine preadipocytes were isolated from healthy newborn calves (5 days old) of Qinchuan cattle by using previously described methods [19,33,34,35,36,37]. The longissimus dorsi muscle area was exposed and extracted under aseptic conditions, and the muscle was first cleaned with 75% ethanol and then incised with the help of a sharp sterile surgical curved scissors. The adipose tissue samples were washed with PBS containing 10% penicillin/streptomycin and immediately transferred to the cell culture room for further processing. The tissues were separated from the blood vessels and connective tissues under a stereo dissecting microscope with the help of sterile forceps. The adipose tissues were then crushed and subjected to enzyme digestion with collagenase I, 0. 25% (Sigma, Shanghai, China) for 1 h at 37 °C in a shaking water bath. The digestive mixture sample was then neutralized with equal volume of 10% FBS (Invitrogen, Carlsbad, CA, USA) and strained by using 100 μm and then 40 μm strainers. The filtrate was centrifuged at 1500 × g for 10 min, and the cell pellet was washed 2 times with DMEM-F/12 medium (Gibco, Grand Island, NY, USA) without serum. The isolated cells pellet was suspended within DMEM-F/12 medium supplemented with 10% FBS, then seeded in 60 mm collagen-coated cell culture plates, and finally incubated at 37 °C and 5% CO2 for 1 h. The medium was changed and then washed 3 times with PBS to remove the debris and free-floating dead cells.

2.3. Transient Transfection, Cell Differentiation, and Staining for Lipid Droplets

The bovine adipocyte cells were seeded in 24-well plates and cultured in DMEM-F/12 (Dulbecco’s modified Eagle’s medium; Gibco, Grand Island, NY, USA) with 10% FBS (Invitrogen, USA) and 1% antibiotics (100 IU/mL penicillin and 100 μg/mL streptomycin), and incubated at 37 °C with 5% CO2. The cells were transiently transfected with bta-miR-149-5p mimic (50 nM) and mimic negative control (50 nM) at 70–90% confluence and density of 1.2 × 105 cells with growth medium without antibiotics and FBS using lipofectamine 3000 transfection reagent (Invitrogen, CA, USA). Briefly, the RNA oligonucleotides and the transfection reagent were diluted separately in Opti-MEM medium (Gibco) and incubated for 10 min at RT. Both the reagents were mixed, combined, and incubated for another 15 min at RT to form the transfection reagent–RNA complexes. The complexes were then added dropwise to the cells and incubated for 24 h before changing to fresh medium.
At 48 h post-transfection, the cells were differentiated with first differentiation media including 0.5 mM hydro cortisol, 1μM dexamethasone, IBMX (3. isobut-1-methylxanthine), and 167 nM insulin, as previously described [38]. After 2 days, the media was replaced with second differentiation medium (DMEM-F/12) supplanted with 10% fetal bovine serum and 167 nM insulin. The cells were fixed and stained with red oil O staining at the ninth day of differentiation. Briefly, the stock solution was prepared by dissolving red oil O staining powder in isopropanol under dark. The solution was filtered and diluted with 60% deionized water and re-filtered to prepare working solution. Before staining, the adipocytes were washed 3 times with 1X PBS and then fixed with 4% paraformaldehyde solution for 30 min at room temperature, then washed again with 1XPBS for 3 times. Then, the cells were stained with working solution of red oil O staining for 30 min at room temperature. The images were captured using an Olympus IX71 microscope (OLYMPUS), and lipid droplets were observed. The lipid droplets were also analyzed using BODIPY staining. The adipocytes were fixed and then incubated for 30 min at room temperature in a 1:1000 dilution of 1 mg/mL BODIPY 493/503 (Invitrogen) diluted with 1XPBS at room temperature. 4,6-Diamidino-2-phenylindole (DAPI) was used for the identification of the nuclei. Images were captures with a Live Cell Imaging System (Nikon Instruments, Europe BV, Kingston, Surrey, UK). ImageJ software (National Institutes of Health, Bethesda, MA, USA) was used for the quantification of lipid accumulation in adipocytes from the fluorescence intensity of the BODIPY stain (green fluorescence).

2.4. RNA Isolation, Construction of cDNA Library, and qRT-PCR

After 48 h of adipocyte treatment with bta-miR-149-5p mimics or mimics NC (negative control), the adipocytes were collected and RNA was extracted through RNAiso Plus kit (Takara Beijing, China). The RNA integrity both quantity and quality were evaluated through optical density (OD) of 260 and ratio of the OD of 260/280 with Nano Quant plate TM (TECAN, Infinite M200 PRO). The RNA quality was further analyzed through 1% agarose gel. Then, the RNA was subjected to reverse transcription for cDNA synthesis through PrimeScriptTM RT reagent kit (perfect real time) with a gDNA eraser (Takara, Beijing, China). The qRT-PCR was performed using SYBR® Premix Ex Taq II kit (Takara, Beijing, China). For the bta-miR-149-5p expression analysis, we used a mircute miRNA isolation kit (TIANGEN, Beijing, China). The first strand cDNA library was constructed using miRcute Plus miRNA First-Strand cDNA kit (TIANGEN, Beijing, China). The bta-miR-149-5p expression was measured through mircute plus miRNA qPCR detection kit (TIANGEN, Beijing, China). For gene expression analysis, GAPDH and β-actin were used, while for miRNA expression, U6 was used as housekeeping genes. The 2−ΔΔCt method was used for the identification of relative mRNA and miRNA expression levels [39]. The PCR primer sequences are shown in Table S1.

2.5. RNA-Seq Library Construction, Quality Control, and Sequencing

The total RNA was extracted with Trizol reagent kit (Invitrogen, Carlsbad, CA, USA), and the RNA quality was analyzed through an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). The quality was further assessed through RNase free agarose gel electrophoresis. After extraction of total RNA, Oligo (dT) beads were used for the enrichment of mRNA, and fragmentation buffer was used for making short fragments, and then reverse-transcribed into cDNA using random primers. The DNA polymerase I, dNTP, RNase H, and buffer were used for the synthesis of second-strand cDNA, and their fragments were purified using QiaQuick PCR extraction kit (Qiagen, Venlo, the Netherlands); poly(A) was added for end repair and then ligated into Illumina sequencing adapters. The ligation products were size-selected by agarose gel electrophoresis, PCR amplified, and sequenced using Illumina HiSeq2500 by Gene Denovo Biotechnology Co. (Guangzhou, China).

2.6. Bioinformatics Analysis

2.6.1. Filtering of Clean Reads

After sequencing, the raw reads including adapters and low quality bases (reads containing more than 10% of unknown nucleotides “N” and the reads containing more than 50% of low quality bases “Q-value ≤ 20”) were filtered through fastp software version 0.18 [40].

2.6.2. Alignment with Ribosome RNA (rRNA) and Reference Genome

The reads were mapped to ribosomal RNA (rRNA) database using the short reads alignment tool Bowtie-2 software version 2.2.8 [41]. Then, the rRNA-mapped reads were removed, and clean reads were used in gene abundance calculation. The reference genome index was created, and clear reads were mapped with reference genome through HISAT2. 2.4 software using “rna_strandness RF” and with all other parameters set as a default [42].

2.6.3. Quantification of Gene Abundance

The fragment per kilobase of transcript per million mapped reads (FPKM) value was calculated for quantification and variations using StringTie software [43,44]. The FPKM formula is shown as follows:
F P k M = 10 6 c N L / 10 3
where FPkM (A) refers to the expression of gene A, “c” denotes the number of fragments mapped to gene A, “N” is the total number of fragments mapped with reference genes, and “L” is number of bases in gene A.

2.6.4. Deferentially Expressed Genes (DEG) Identification

The DEGs between the 2 groups was analyzed with DESeq2 software [45], and edgeR software was used for variation between the 2 samples [46]. The transcripts with FDR (false discovery rate) below 0.05 and absolute fold change with ≥2 were considered to be DEGs.

2.6.5. GO Enrichment Analysis

The DEGs were mapped into GO terms through GO enrichment analysis using the GO database available at http://www.geneontology.org/ (accessed on 20 April 2021). The significantly enriched GO terms in the DEGs correspond to the respective biological functions. Hypergeometric test was used for the calculation of gene numbers for each term in comparison with the reference genome using the following formula:
P = 1 − ∑ i = 0 m − 1 m i N − M n − i N n
where “N” denotes the number of all transcripts within GO annotation, “n” is the number of DEGs in N, “M” is the number of all transcripts annotated without respective GO terms, and “m” is the number of DEGs in M. The p-value is from FDR correction (FDR ≤ 0.05) and was used as a threshold. The GO terms within these conditions were considered as significantly enriched GO terms in DEGs.

2.6.6. KEGG Pathway and Reactome Enrichment Analysis

The pathway enrichment analysis identified significantly enriched metabolic or signal transduction pathways in DEGs compared with the reference whole genome using the following formula:
P = 1 − ∑ i = 0 m − 1 m i N − M n − i N n
where “N” denotes the number of all transcripts within KEGG annotation, “n” is the number of DEGs in N, “M” is the number of all transcripts annotated without respective KEGG pathways, and “m” is the number of DEGs in M. The p-value is from FDR correction (FDR ≤ 0.05) and was used as a threshold. The pathways within these conditions were considered as significantly enriched KEGG pathways in DEGs. The reactome is a free online database of biological pathways [47,48]. The core unit of the reactome data model is the reaction. The nucleic acids, proteins, complexes, and small molecules that participate in reactions form a network of biological interactions and are grouped into pathways. Examples of biological pathways in the reactome include signaling, innate and acquired immune function, transcriptional regulation, translation, apoptosis, and classical intermediary metabolism. Reactome enrichment analysis identified significantly enriched reactions in DEGs compared with the whole genome background.
The same formula was used for reactome enrichment analysis.

2.7. Protein–Protein Interaction

Protein–protein interaction network was identified using String v10 [49], which determined genes as nodes and interaction as lines in a network. The network file was visualized using Cytoscape (v3.7.1) (Maryland, USA) software to present a core and hub gene biological interaction [50].

2.8. Statistical Analysis

The ANOVA was used for the statistical analysis through SAS version 8.1 software (SAS Institute Inc., Cary, NC, USA). The graphs were designed through Graph-Pad Prism6 (GraphPad, San Diego, CA, USA) software, with mean ± SEM. The variations at * p < 0.05 and ** p < 0.01 were considered to be statistically significant.

3. Results

3.1. Transfection Efficiency and Quality Evaluation of the Samples

To explicate regulatory function of bovine bta-miR-5p in adipogenesis, we analyzed the mimic efficiency of bta-miR-149-5p in bovine adipocytes. The expression level of bta-miR-149-5p in cells transfected with bta-miR-149-5p mimic was significantly (p < 0.01) increased as compared to mimic NC (Figure 1A). Before analyzing the mimic efficiency, we detected the transfection efficiency of the miRNA through FAM-labeled negative control miRNA (gene pharma, Shanghai, China). The transfection efficiency was evaluated at 6-, 12-, and 24-h intervals. The transfected adipocytes were observed under fluorescence microscope, and the photographs were captured though Olympus IX71 microscope (OLYMPUS, Tokyo, Japan) (Figure 1B–D). These findings confirmed the accuracy and success of the transfection experiment performed in the current study and authenticated the reliability of data in subsequent experiments.
The RNA quality shows RNA integrity number (RIN) ranged from 7.5 to 10 (Figure 2A). These figures show integrity of the RNA samples used in the current study. Moreover, sequenced clean reads were filtered (Figure 2C) and aligned with ribosome RNA (rRNA) and reference genome (Figure 2D) for purification and normality of the data for subsequent analysis. The expression abundance also confirmed the significance of the transcriptomic analysis of the samples (Figure 2E).

3.2. Bta-miR-149-5p Regulated Adipogenesis

The lipid droplets (LDs) on day 9 of differentiation were assessed using BODIPY staining (Figure 3). The adipocytes transfected with bta-miR-149-5p mimics exhibited decreased level of LDs as compared with adipocytes transfected with mimic NC. The histogram pixel intensity also showed significant (p < 0.01) reduction in LDs in adipocytes transfected with bta-miR-149-5p mimics as compared with mimic NC. These results illustrate the fact that bta-miR-149-5p inhibited differentiation of bovine adipocytes.
Deep RNA-sequencing analysis revealed that bta-miR-149-5p significantly affected expression of adipogenic marker genes in differentiated bovine adipocytes transfected with bta-miR-149-5p mimic and mimic NC (Figure 4A–C). In total, there were 115 DEGs (differentially expressed genes), with 1.5 ± |FC and (log2 (fc) > 0.58496/log 2 (fc) > 0.58496 and p-value (<0.05). These 115 DEGs contain 72 upregulated and 43 downregulated genes in bovine differentiated adipocytes (Table S2).
The bta-miR-149-5p significantly changed the GO terms in bovine adipocyte transfected with bta-miR-149-5p mimics and mimics NC (Figure 5A). The unigenes and GO term biological processes were the most annotated unigene contributor part with 80.08%, followed by cellular component (13.4%) and molecular function (6.7%). The top 20 GO terms due to DEGs are shown in Figure 5B. The DEGs significantly altered the GO term GO:0044421 extracellular region part and GO:0005576 extracellular region in cellular component GO ontologies. The GO terms altered in molecular function ontologies in the top 20 GO terms were GO:0005515 protein binding and GO:0005488 binding. The most important GO terms altered in cellular component ontologies in the top 20 GO terms were GO:0031012 extracellular matrix, GO:0005125 cytokine activity, GO:0005102 receptor binding, GO:0016787 hydrolase activity, GO:0022610 biological adhesion, GO:0065008 regulation of biological quality, GO:0031589 cell−substrate adhesion, GO:0031323 regulation of cellular metabolic process, GO:0042326 negative regulation of phosphorylation, and GO:0042325 regulation of phosphorylation in bovine adipocytes.
Different expression trends of DEGs were found in different GO terms (Figure 6, Table S3). In GO terms of biological processes, 13 genes were upregulated and 6 were downregulated in GO:0031323 regulation of cellular metabolic process; 12 genes were upregulated and 9 were downregulated in GO:0065008 regulation of biological quality; 9 genes were upregulated and 9 were downregulated in GO:0022610 biological adhesion; and 3 genes were upregulated and 1 gene was downregulated each in GO:0031589 cell-substrate adhesion, GO:0008643 carbohydrate transport, GO:0015749 monosaccharide transport, and GO:0048645 organ formation, respectively. A single gene was up- and downregulated each in GO:0072502 cellular trivalent inorganic anion homeostasis, GO:0002639 positive regulation of immunoglobulin production, GO:0002891 positive regulation of immunoglobulin-mediated immune response, and GO:0045830-positive regulation of isotype switching. In biological function ontologies, 22 genes were upregulated and 13 were downregulated in GO:0005515 protein binding, 9 genes were upregulated and 3 were downregulated in GO:0005102 receptor binding, and all 4 genes were upregulated in GO:0005125 cytokine activity. In the category of cellular component ontologies, 14 genes were upregulated and 4 were downregulated, each in GO:0005576 extracellular region and GO:0044421 extracellular region part, respectively. All six genes were upregulated in GO:0031012 extracellular matrix.
The pathway analysis of the DEGs was performed to explore the pathways and molecular interactions in bovine adipocytes transfected with bta-miR-149-5p mimics or mimics NC. The top 20 significantly enriched KEGG pathways, based on the Q < 0.05 value calculated by hypergeometric test, demonstrated the roles of DEGs in adipocyte differentiation (Figure 7). DEGs altered the top KEGG_A_class enrichment pathways including metabolism (global and overview maps, nucleotide, amino acids, carbohydrates, and lipid metabolism), genetic information processes (signal transduction, signaling molecules, and interaction), cellular process (cellular community eukaryotes, cell growth and death, transport and catabolism, cell motility, cell motility), and organismal systems (immune system, endocrine system, nervous system, development, digestive system, circulatory system); under diseases, it also regulated infectious diseases, cancers, immune diseases, endocrine and metabolic diseases, and cardiovascular diseases. These KEGG_A class pathways perform crucial roles in adipogenesis.
The top 20 pathways affected by bta-miR-149-5p through DEGs within differentiated bovine adipocytes were ko04060 cytokine−cytokine receptor interaction, ko04668 TNF signaling pathway, ko04015 Rap1 signaling pathway, ko04392 Hippo signaling pathway—multiple species, ko04530 tight junction, diseases, ko04923 regulation of lipolysis in adipocyte, ko04062 chemokine signaling pathway, and ko04911 insulin secretion.
The different KEGG_A_class pathways were significantly altered in adipocytes transfected with bta-miR-149-5p mimics and its negative control NC (Figure 8, Table 1 and Table S4). The highly effected KEGG_A_class was the human disease that altered the most important KEGG_B_class pathways including immune diseases, infectious diseases, cardiovascular diseases, and cancers. The second most affected KEGG_A_class pathway was organismal systems, which significantly altered the KEGG_B_class pathways including immune systems, endocrine system, and excretory system; the most important corresponding pathways altered by these KEGG_B class were IL-17 signaling pathway, regulation of lipolysis in adipocyte, chemokine signaling pathway, insulin secretion, and aldosterone-regulated sodium reabsorption. The role of these pathways in adipogenesis is well documented. The environmental information processing and cellular processes were also the most significantly changed KEGG_A_class pathways by bta-miR-149-5p in bovine adipocytes, the former changed signal transduction, signaling molecules, and interaction, and the latter changed cell growth and death, cell motility, and cellular community—eukaryotes KEGG_B_class pathways. These KEGG_B classes altered the most important pathways, including TNF signaling pathway, Hippo signaling pathway—multiple species, Rap1 signaling pathway, tight junction, and focal adhesion. Moreover, metabolism was also found to be one of the mostly affected pathways in the KEGG_A_class, and the KEGG_B_class pathways in this category included lipid metabolism, carbohydrate metabolism, amino acid metabolism, metabolism of cofactors and vitamins, nucleotide metabolism, and global and overview maps.
Moreover, we explored the relationship between the DEGs and signaling pathways through Cytoscape software (Figure 9). The KEGG pathways regulated by these DEGs were PI3K-Akt signaling pathway, calcium signaling pathway, pathways in cancer, MAPK signaling pathway, lipid metabolism/metabolic pathway, PPAR signaling pathway, AMPK signaling pathway, TGF-beta signaling pathway, cAMP signaling pathway, cholesterol metabolism, Wnt signaling pathway, and FoxO signaling pathway. Furthermore, the expression trend of the DEGs in these pathways were also exploited. These data indicate the role of bta-miR-149-5p in adipogenesis through interactive regulation of DEGs and their respective KEGG pathways in bovine adipocytes.
In addition to this, reactome biological pathways were also studied in terms of the DEGs. As shown in Figure 10, the most important pathways in top 20 reactome enrichment pathways were R−BTA−373813 receptor CXCR2 binding ligands CXCL1 to 7, R−BTA−373791 receptor CXCR1 binding CXCL6 and CXCL8 ligands, R−BTA−210991 basigin interactions, R−BTA−380108 chemokine receptors binding chemokines, R−BTA−445699, ATP hydrolysis by myosin, R−BTA−445704 calcium binding caldesmon, and R−BTA−5669034 TNFs binding their physiological receptors. These reactome pathways play crucial roles in adipogenesis.

3.3. Validation and Identification of Key DEGs during Adipogenesis through qRT-PCR

The RNA sequence results were validated through qRT-PCR analysis of key marker genes selected from the list of DEGs (Figure 11). The mRNA levels of CCND2, KLF6, ACSL1 Cdk2, SCD, SIK2, and ZEB1 were significantly (p < 0.01) downregulated in adipocytes transfected with bta-miR-149-5p mimics. However, upregulation of bta-miR-149-5p significantly (p < 0.01) enriched relative mRNA expression of KLF12 and LIPG genes in bovine adipocytes. These results are in line with the results of sequencing data.

4. Discussion

Improvement of intramuscular fat is a challenge for the experts of animal science for the improvement of meat quality traits, and obesity is a global issue—a root cause of certain fatal human disorders [51]. Adipogenesis is a complex biochemical process, characterized by adipocyte proliferation, maturation, and differentiation of preadipocytes. The preadipocytes that are derived from the existing pool of adipocytes differentiate in response to proper signals [52]. Therefore, it is imperative to illuminate the underlying molecular mechanism of adipogenesis. In the recent past, a series of studies proved the crucial functions of microRNAs in the proliferation, differentiation, and function of adipocytes. High-throughput sequencing and microarray analysis revealed the roles of novel miRNAs in animals [53]. Previously, the roles of let-7 [54], miR-21a [55], miR-23a [56,57], miR-27a/b [58,59], miR-130 [60], miR-185 [61], miR-448 [62], and miR-709 [63] in lipogenesis negatively regulated adipocyte differentiation, whereas miR-21 [64,65], miR-26b [66], miR-143 [67], miR-148 [68], miR-183 [69], miR-424 [14], and miR-425 [70] positively regulated adipocyte differentiation. Previously, we confirmed the role of bta-miR-149-5p bovine adipocyte proliferation and differentiation [32]. These findings imply that bta-miR-149-5p is a novel hallmark for dealing with the problem of low intramuscular fat in cattle breeding programs.
In the present study, deep RNA-sequencing analysis revealed that bta-miR-149-5p significantly affected expression of adipogenic marker genes in differentiated bovine adipocytes transfected with bta-miR-149-5p mimic and mimic NC. The unigenes and GO term biological processes were the most annotated unigene contributor part with 80.08%, followed by cellular component of 13.4%, and molecular function of 6.7%. Interestingly, in the present study, the DEGs were upregulated in all three GO categories, namely, in biological processes, biological function, and cellular component ontologies. These findings show that bta-miR-149-5p regulates the functional roles of DEGs in bovine adipocytes through regulation of GO terms. Furthermore, the pathway analysis of the DEGs was performed to exhibit the pathways and molecular interactions in bovine adipocytes transfected with bta-miR-149-5p mimics or mimics NC. Strikingly, the pathway enrichment analysis explored the KEGG_A_classes, including disease and metabolism regulated by the DEGs in adipocytes transfected with bta-miR-149-5p, which specifically categorized the most important KEGG_B class pathways with proven roles in metabolism and adipogenesis. The pathway-based studies helped to exploit the gene roles in biological functions [46,71]. The adipocytes perform various functions, particularly energy storage in the form of triglycerides [72]. Adipose is also an endocrine organ, secreting various bioactive molecules such as adipokines, which regulate metabolic functions in body system. Recently, intensive research has been conducted to show miRNA’s potential regulatory roles in metabolism [73]. Moreover, the expression trend of DEGs in these pathways exhibited an upregulated trend, and the number of regulated DEGs were more than the number of downregulated genes in these pathways.
In addition to this, the specific KEGG pathways regulated by these DEGs were PI3K-Akt signaling pathway, calcium signaling pathway, pathways in cancer, MAPK signaling pathway, lipid metabolism/metabolic pathway, cholesterol metabolism, PPAR signaling pathway, AMPK signaling pathway, TGF-beta signaling pathway, cAMP signaling pathway, Wnt signaling pathway, and FoxO signaling pathway. These data indicate the role of bta-miR-149-5p in adipogenesis through interactive regulation of DEGs and their respective KEGG pathways in bovine adipocytes. All these pathways have proven roles in adipogenesis. Previously, it has been reported that PI3K-Akt signaling pathway plays an important role in adipose tissue function and energy metabolism through regulation of the genes expressed in lipid metabolism [74,75,76]. Previous literature has shown the association of FoxO signaling, PI3K-Akt signaling, and cAMP signaling pathways with lipid metabolism [77,78,79]. The cAMP and calcium signaling pathways act as the second messengers in eukaryotic cells [80] and regulate the intracellular lipolysis in adipocytes through calcium ions and cAMP levels [81]. There is substantial evidence that MAPK signaling pathway [82,83], lipid metabolism, insulin signaling [84], AMPK signaling pathway, Wnt signaling pathway [85,86,87,88], and TGF-beta signaling pathway [89,90] perform crucial roles in adipogenesis. The results of the present study support the interactive roles of the signaling pathways in the regulation of adipogenesis.
Interestingly, the pathway enrichment analysis in the present study also exhibited the role of bta-miR-149-5p in cancer. Previously, the roles of miR-149 in osteosarcoma [29], cancers of the digestive system [21], colorectal cancer [22], hepatocellular cancer [23,24], respiratory system cancer [25], breast cancer [26], and sebaceous gland carcinoma of the eyelid [27] have been well documented. Furthermore, in the present study, overexpression of bta-miR-149-5p significantly (p < 0.05) upregulated LIPG (endothelial lipase) and KLF12 expression, which perform vital roles in cancer [91,92]. However, cancer study is beyond the scope of the present study, yet there is a very close relation of adipocytes and cancers. Adipocytes are the primary stromal cells, performing a dynamic role in the cancer microenvironment. Cancer-associated adipocytes (CAAs) are found next to the cancer cells and communicate with the tumor cells via secreting several factors that mediate local and systemic effects. Recently, miRNAs have been considered mediators of intercellular communication in the cancer microenvironment. Tumor cells coordinate with stromal cells and other extracellular matrix (ECM) elements and establish a comfortable niche to grow in in order to escape the immune system. Inside the cancer microenvironment of the cells, the miRNAs influence and take over the physiological processes of surrounding cells and thus promote cancer progression [93]. However, these findings need further research in order to fully explore this molecular mechanism of microRNA in cancer-associated adipocytes (CAAs).
Moreover, the bta-miR-149-5p significantly (p < 0.01) downregulated the mRNA levels of adipogenic marker genes such as CCND2, KLF6, ACSL1, Cdk2, SCD, SIK2, and ZEB1 in bovine adipocytes. These findings validated the transcriptomic data and also proved the functional roles of bta-miR-149-5p in bovine adipogenesis through regulation of these adipogenic marker genes’ expression. The cyclin D2 (CCND2) and Cdk2 are well-known proliferation marker genes, enriched in early differentiation and then again overexpressed in the later stage of adipocyte differentiation [94,95]. Kruppel-like factor 6 (KLF6) is a transcription factor and is a potent regulator of adipocyte differentiation [96]. ACSL1 is a proven mediator of lipid metabolism, cellular lipid content, and insulin sensitivity in adipocyte differentiation [97]. The crucial role of the SCD gene is well established in the differentiation adipocyte [98] and fatty acid composition during adipogenesis [99]. SIK2 plays a critical role in the regulation of lipid homeostasis and adipogenesis via cAMP signaling pathway [100]. Previously, our research group validated the role of ZEB1 transcription factor in bovine adipocytes, with it found to be potential regulatory motif in the promoter of an adipogenic marker gene (ABHD5) that plays vital role in bovine adipogenesis [101]. These results of qRT-PCR are in line with the results of RNA-Seq data, which further validate and authenticate the present study. Taken together, these findings suggest that the DEGs in the present study play pivotal roles in bovine adipogenesis. Generally, the present study provides valuable evidence for expanding our understanding about the molecular mechanism of bta-miR-149-5p in the regulation of lipid metabolism and adipogenesis.

5. Conclusions

These data indicate the role of bta-miR-149-5p in adipogenesis through interactive regulation of DEGs and their respective KEGG pathways in bovine adipocytes. In conclusion, our results suggest that bta-miR-149-5p regulates lipid metabolism in bovine adipocytes. The results of this study provide a basis for studying the function and molecular mechanism of the bta-miR-149-5p in regulating adipogenesis.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ani11051207/s1, Tables S1–S4.

Author Contributions

Conceptualization, H.G., R.K. and L.Z.; data curation, H.G. and R.K.; formal analysis, R.K., S.M.S. and H.K.; funding acquisition, H.G. and L.Z.; investigation, H.G. and R.K.; methodology, R.K. and S.H.A.R.; project administration, and L.Z.; resources, and L.Z.; software, R.K.; supervision, L.Z.; writing—original draft, H.G. and R.K.; writing—review and editing, S.H.A.R., S.M.S., H.K., S.B.K. and A.H.A.E.-A. All authors have read and agreed to the published version of the manuscript.

Funding

The study was supported by the National Key Research and Development Program of China (no. 2018YFD0501700), National Science and Technology Support Projects (no. 2015BAD03B04), and the National Modern Agricultural Industry Special Program (no. CARS-37).

Institutional Review Board Statement

The experimental animals were dealt as per standard operating procedures (SOPs) formulated by the Chinese Council of Animals Care, and further approved by “Experimental Animal Management Committee (EAMC)” of the Northwest Agricultural and Forestry University. The sampling was performed after human killing of the experimental animals at the National Beef Cattle Improvement Research Center, Yangling, Shaanxi, China.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are available with the manuscript as supplementary files.

Acknowledgments

All the authors of the manuscript are immensely grateful to their respective universities and institutes for their technical assistance and valuable support in the completion of this research project.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Sethi, J.K.; Vidal-Puig, A.J. Thematic review series: Adipocyte Biology. Adipose tissue function and plasticity orchestrate nutritional adaptation. J. Lipid Res. 2007, 48, 1253–1262. [Google Scholar] [CrossRef] [Green Version]
  2. Rutkowski, J.M.; Stern, J.H.; Scherer, P.E. The cell biology of fat expansion. J. Cell Biol. 2015, 208, 501–512. [Google Scholar] [CrossRef] [Green Version]
  3. Ghaben, A.L.; Schere, P.E. Adipogenesis and metabolic health. Nature reviews. Mol. Cell Biol. 2019, 20, 242–258. [Google Scholar]
  4. Martins, T.S.; Sanglard, L.M.P.; Silva, W.; Chizzotti, M.L.; Rennó, L.N.; Serão, N.V.L.; Silva, F.F.; Guimarães, S.E.F.; Ladeira, M.M.; Dodson, M.V.; et al. Molecular Factors Underlying the Deposition of Intramuscular Fat and Collagen in Skeletal Muscle of Nellore and Angus Cattle. PloS ONE 2015, 10, e0139943. [Google Scholar] [CrossRef] [PubMed]
  5. Cai, H.; Li, M.; Sun, X.; Plath, M.; Li, C.; Lan, X.; Lei, C.; Huang, Y.; Bai, Y.; Qi, X.; et al. Global Transcriptome Analysis During Adipogenic Differentiation and Involvement of Transthyretin Gene in Adipogenesis in Cattle. Front. Genet. 2018, 9, 463–473. [Google Scholar] [CrossRef] [PubMed]
  6. Fu, Y.-Y.; Chen, K.-L.; Li, H.-X.; Zhou, G.-H. The adipokine Chemerin induces lipolysis and adipogenesis in bovine intramuscular adipocytes. Mol. Cell. Biochem. 2016, 418, 39–48. [Google Scholar] [CrossRef]
  7. He, L.; Hannon, G.J. MicroRNAs: Small RNAs with a big role in gene regulation. Nat. Rev. Genet. 2004, 5, 522–531. [Google Scholar] [CrossRef]
  8. Thomas, M.; Lieberman, J.; Lal, A. Desperately seeking microRNA targets. Nat. Struct. Mol. Biol. 2010, 17, 1169–1174. [Google Scholar] [CrossRef] [PubMed]
  9. Raza, S.H.A.; Kaster, N.; Khan, R.; Abdelnour, S.A.; El-Hack, M.E.A.; Khafaga, A.F.; Taha, A.; Ohran, H.; Swelum, A.A.; Schreurs, N.M.; et al. The Role of MicroRNAs in Muscle Tissue Development in Beef Cattle. Genes 2020, 11, 295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Esau, C.; Kang, X.; Peralta, E.R.; Hanson, E.; Marcusson, E.G.; Ravichandran, L.V.; Sun, Y.; Koo, S.; Perera, R.; Jain, R. MicroRNA-143 Regulates Adipocyte Differentiation. J. Biol. Chem. 2004, 279, 52361–52365. [Google Scholar] [CrossRef] [Green Version]
  11. Wei, Y.; Cui, Y.F.; Tong, H.L.; Zhang, W.W.; Yan, Y.Q. MicroRNA-2400 promotes bovine preadipocyte proliferation. Biochem. Biophys. Res. Commun. 2016, 478, 1054–1059. [Google Scholar] [CrossRef]
  12. Wang, L.; Zhang, S.; Cheng, G.; Mei, C.; Li, S.; Zhang, W.; Junjvlieke, Z.; Zan, L. MiR-145 reduces the activity of PI3K/Akt and MAPK signaling pathways and inhibits adipogenesis in bovine preadipocytes. Genomics 2020, 112, 2688–2694. [Google Scholar] [CrossRef] [PubMed]
  13. Shen, L.; Li, Q.; Wang, J.; Zhao, Y.; Niu, L.; Bai, L.; Shuai, S.; Li, X.; Zhang, S.; Zhu, L. miR-144-3p Promotes Adipogenesis Through Releasing C/EBPalpha From Klf3 and CtBP2. Front. Genet. 2018, 9, 677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Wang, L.; Zhang, S.; Zhang, W.; Cheng, G.; Khan, R.; Junjvlieke, Z.; Li, S.; Zan, L. miR-424 Promotes Bovine Adipogenesis Through an Unconventional Post-Transcriptional Regulation of STK11. Front. Genet. 2020, 11, 145. [Google Scholar] [CrossRef] [Green Version]
  15. Cioffi, M.; Vallespinosserrano, M.; Trabulo, S.; Fernandezmarcos, P.J.; Firment, A.N.; Vazquez, B.N.; Vieira, C.R.; Mulero, F.; Camara, J.A.; Cronin, U.P. MiR-93 Controls Adiposity via Inhibition of Sirt7 and Tbx3. Cell Rep. 2015, 12, 1594–1605. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. 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] [CrossRef] [Green Version]
  17. Kulyte, A.; Kwok, K.H.M.; De Hoon, M.J.L.; Carninci, P.; Hayashizaki, Y.; Arner, P.; Arner, E. MicroRNA-27a/b-3p and PPARG regulate SCAMP3 through a feed-forward loop during adipogenesis. Sci. Rep. 2019, 9, 1–9. [Google Scholar] [CrossRef] [Green Version]
  18. Liu, S.-Y.; Zhang, Y.-Y.; Gao, Y.; Zhang, L.-J.; Chen, H.-Y.; Zhou, Q.; Chai, M.-L.; Li, Q.-Y.; Jiang, H.; Yuan, B.; et al. MiR-378 Plays an Important Role in the Differentiation of Bovine Preadipocytes. Cell. Physiol. Biochem. 2015, 36, 1552–1562. [Google Scholar] [CrossRef]
  19. Wang, Y.; Zhang, Y.; Su, X.; Wang, H.; Yang, W.; Zan, L. Cooperative and Independent Functions of the miR-23a~27a~24-2 Cluster in Bovine Adipocyte Adipogenesis. Int. J. Mol. Sci. 2018, 19, 3957. [Google Scholar] [CrossRef] [Green Version]
  20. Ma, X.; Wei, D.; Cheng, G.; Li, S.; Wang, L.; Wang, Y.; Wang, X.; Zhang, S.; Wang, H.; Zan, L. Bta-miR-130a/b regulates preadipocyte differentiation by targeting PPARG and CYP2U1 in beef cattle. Mol. Cell. Probes 2018, 42, 10–17. [Google Scholar] [CrossRef]
  21. Wang, N.; Zhou, P.; Chen, Y.; Qu, H.; Lu, K.; Xia, J. MicroRNA-149: A review of its role in digestive system cancers. Pathol. Res. Pract. 2020, 216, 153266. [Google Scholar] [CrossRef]
  22. Liu, K.; Mou, Y.; Shi, X.; Liu, T.; Chen, Z.; Zuo, X. Circular RNA 100146 Promotes Colorectal Cancer Progression by the MicroRNA 149/HMGA2 Axis. Mol. Cell. Biol. 2021, 41, e00445. [Google Scholar] [CrossRef]
  23. Zhang, Y.; Guo, X.; Xiong, L.; Yu, L.; Li, Z.; Guo, Q.; Li, Z.; Li, B.; Lin, N. Comprehensive analysis of microRNA-regulated protein interaction network reveals the tumor suppressive role of microRNA-149 in human hepatocellular carcinoma via targeting AKT-mTOR pathway. Mol. Cancer 2014, 13, 253. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Yu, Y.-X.; Ge, T.-W.; Zhang, P. Circular RNA circGFRA1 promotes angiogenesis, cell proliferation and migration of hepatocellular carcinoma by combining with miR-149. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 11058–11064. [Google Scholar] [PubMed]
  25. Li, H.; Cui, Z.; Lv, X.; Li, J.; Gao, M.; Yang, Z.; Bi, Y.; Zhang, Z.; Wang, S.; Li, S.; et al. Long Non-coding RNA HOTAIR Function as a Competing Endogenous RNA for miR-149-5p to Promote the Cell Growth, Migration, and Invasion in Non-small Cell Lung Cancer. Front. Oncol. 2020, 10, 528520. [Google Scholar] [CrossRef]
  26. Qi, C.; Qin, X.; Zhou, Z.; Wang, Y.; Yang, Q.; Liao, T. Circ_0072995 Promotes Cell Carcinogenesis via Up-Regulating miR-149-5p-Mediated SHMT2 in Breast Cancer. Cancer Manag. Res. 2020, 12, 11169–11181. [Google Scholar] [CrossRef]
  27. Hirano, T.; Yunoki, T.; Furusawa, Y.; Tabuchi, Y.; Hayashi, A. Bioinformatics analysis of the microRNA-mRNA network in sebaceous gland carcinoma of the eyelid. Mol. Med. Rep. 2020, 23, 1. [Google Scholar] [CrossRef] [PubMed]
  28. Ding, H.; Zheng, S.; Garcia-Ruiz, D.; Hou, D.; Wei, Z.; Liao, Z.; Li, L.; Zhang, Y.; Han, X.; Zen, K.; et al. Fasting induces a subcutaneous-to-visceral fat switch mediated by microRNA-149-3p and suppression of PRDM. Nat. Commun. 2016, 7, 11533. [Google Scholar] [CrossRef]
  29. Xu, R.-D.; Feng, F.; Yu, X.-S.; Liu, Z.-D.; Lao, L.-F. miR-149-5p inhibits cell growth by regulating TWEAK/Fn14/PI3K/AKT pathway and predicts favorable survival in human osteosarcoma. Int. J. Immunopathol. Pharmacol. 2018, 32, 2058738418786656. [Google Scholar] [CrossRef] [PubMed]
  30. Malpeli, G.; Barbi, S.; Zupo, S.; Tosadori, G.; Scardoni, G.; Bertolaso, A.; Sartoris, S.; Ugel, S.; Vicentini, C.; Fassan, M.; et al. Identification of microRNAs implicated in the late differentiation stages of normal B cells suggests a central role for miRNA targets ZEB1 and TP. Oncotarget 2017, 8, 11809–11826. [Google Scholar] [CrossRef]
  31. Junaith, S.M.; Ameena, H.; Patricia, S.P.; Aladin, M.B. MicroRNA-149 inhibits PARP-2 and promotes mitochondrial bio-genesis via SIRT-1/PGC-1α network in skeletal muscle. Diabetes 2014, 63, 1546–1559. [Google Scholar]
  32. Khan, R.; Raza, S.H.A.; Junjvlieke, Z.; Wang, X.; Wang, H.; Cheng, G.; Mei, C.; Elnour, I.E.; Zan, L. Bta-miR-149-5p inhibits proliferation and differentiation of bovine adipocytes through targeting CRTCs at both transcriptional and posttranscriptional levels. J. Cell. Physiol. 2020, 235, 5796–5810. [Google Scholar] [CrossRef] [PubMed]
  33. Meissburger, B.; Perdikari, A.; Moest, H.; Müller, S.; Geiger, M.; Wolfrum, C. Regulation of adipogenesis by paracrine factors from adipose stromal-vascular fraction-a link to fat depot-specific differences. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2016, 1861, 1121–1131. [Google Scholar] [CrossRef] [Green Version]
  34. Wei, S.; Zhang, M.; Zheng, Y.; Yan, P. ZBTB16 Overexpression Enhances White Adipogenesis and Induces Brown-Like Adipocyte Formation of Bovine White Intramuscular Preadipocytes. Cell. Physiol. Biochem. 2018, 48, 2528–2538. [Google Scholar] [CrossRef] [PubMed]
  35. Khan, R.; Raza, S.H.A.; Schreurs, N.; Xiaoyu, W.; Hongbao, W.; Ullah, I.; Rahman, A.; Suhail, S.M.; Khan, S.; Linsen, Z. Bioinformatics analysis and transcriptional regulation of TORC1 gene through transcription factors NRF1 and Smad3 in bovine preadipocytes. Genomics 2020, 112, 1575–1587. [Google Scholar] [CrossRef]
  36. Wang, X.; Khan, R.; Raza, S.H.A.; Li, A.; Zhang, Y.; Liang, C.; Yang, W.; Wu, S.; Zan, L. Molecular characterization of ABHD5 gene promoter in intramuscular preadipocytes of Qinchuan cattle: Roles of Evi1 and C/EBPα. Gene 2019, 690, 38–47. [Google Scholar] [CrossRef] [PubMed]
  37. Khan, R.; Raza, S.H.A.; Junjvlieke, Z.; Wang, H.; Cheng, G.; Smith, S.B.; Jiang, Z.; Li, A.; Zan, L. RNA-seq reveal role of bovine TORC2 in the regulation of adipogenesis. Arch. Biochem. Biophys. 2020, 680, 108236. [Google Scholar] [CrossRef]
  38. Junjvlieke, Z.; Mei, C.; Khan, R.; Zhang, W.; Hong, J.; Wang, L.; Li, S.; Zan, L. Transcriptional regulation of bovine elongation of very long chain fatty acids protein 6 in lipid metabolism and adipocyte proliferation. J. Cell. Biochem. 2019, 120, 13932–13943. [Google Scholar] [CrossRef] [PubMed]
  39. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  40. Chen, S.; Zhou, Y.; Chen, Y.; Gu, J. fastp: An ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 2018, 34, i884–i890. [Google Scholar] [CrossRef]
  41. Langmead, B.; Salzberg, S.L. Fast gapped-read alignment with Bowtie. Nat. Methods 2012, 9, 357–359. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kim, D.; Langmead, B.; Salzberg, S.L. HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 2015, 12, 357–360. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Pertea, M.; Pertea, G.M.; Antonescu, C.M.; Chang, T.-C.; Mendell, J.T.; Salzberg, S.L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 2015, 33, 290–295. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Pertea, M.; Kim, D.; Pertea, G.M.; Leek, J.T.; Salzberg, S.L. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 2016, 11, 1650–1667. [Google Scholar] [CrossRef] [PubMed]
  45. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [Green Version]
  46. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  47. Croft, D.; O’Kelly, G.; Wu, G.; Haw, R.; Gillespie, M.; Matthews, L.; Caudy, M.; Garapati, P.; Gopinath, G.; Jassal, B.; et al. Reactome: A database of reactions, pathways and biological processes. Nucleic Acids Res. 2010, 39, D691–D697. [Google Scholar] [CrossRef]
  48. Fabregat, A.; Jupe, S.; Matthews, L.; Sidiropoulos, K.; Gillespie, M.; Garapati, P.; Haw, R.; Jassal, B.; Korninger, F.; May, B.; et al. The Reactome Pathway Knowledgebase. Nucleic Acids Res. 2018, 46, D649–D655. [Google Scholar] [CrossRef]
  49. Van Der Auwera, G.A.; Carneiro, M.O.; Hartl, C.; Poplin, R.; Del Angel, G.; Levy-Moonshine, A.; Jordan, T.; Shakir, K.; Roazen, D.; Thibault, J.; et al. From FastQ Data to High-Confidence Variant Calls: The Genome Analysis Toolkit Best Practices Pipeline. Curr. Protoc. Bioinform. 2013, 43, 11.10.1–11.10.33. [Google Scholar] [CrossRef]
  50. Szklarczyk, D.; Franceschini, A.; Wyder, S.; Forslund, K.; Heller, D.; Huerta-Cepas, J.; Simonovic, M.; Roth, A.; Santos, A.; Tsafou, K.P.; et al. STRING v10: Protein–protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015, 43, D447–D452. [Google Scholar] [CrossRef]
  51. Guilherme, A.; Virbasius, J.V.; Puri, V.; Czech, M.P. Adipocyte dysfunctions linking obesity to insulin resistance and type 2 diabetes. Nat. Rev. Mol. Cell Biol. 2008, 9, 367–377. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Ahn, J.; Lee, H.; Jung, C.H.; Jeon, T.I.; Ha, T.Y. Micro RNA-146b promotes adipogenesis by suppressing the SIRT 1-FOXO 1 cascade. EMBO Mol. Med. 2013, 5, 1602–1612. [Google Scholar] [CrossRef] [PubMed]
  53. Ambros, V. The functions of animal microRNAs. Nature 2004, 431, 350. [Google Scholar] [CrossRef] [PubMed]
  54. Sun, T.; Fu, M.; Bookout, A.L.; Kliewer, S.A.; Mangelsdorf, D.J. MicroRNA let-7 Regulates 3T3-L1 Adipogenesis. Mol. Endocrinol. 2009, 23, 925–931. [Google Scholar] [CrossRef] [Green Version]
  55. Xie, X.; Song, J.; Li, G. MiR-21a-5p suppresses bisphenol A-induced pre-adipocyte differentiation by targeting map2k3 through MKK3/p38/MAPK. Biochem. Biophys. Res. Commun. 2016, 473, 140–146. [Google Scholar] [CrossRef]
  56. Huang, Y.; Huang, J.; Qi, R.; Wang, Q.; Wu, Y.; Wang, J. Effects of MicroRNA-23a on Differentiation and Gene Expression Profiles in 3T3-L1 Adipocytes. Genes 2016, 7, 92. [Google Scholar] [CrossRef] [Green Version]
  57. Shen, L.; Zhang, Y.; Du, J.; Chen, L.; Luo, J.; Li, X.; Li, M.; Tang, G.; Zhang, S.; Zhu, L. MicroRNA-23a regulates 3T3-L1 adipocyte differentiation. Gene 2016, 575, 761–764. [Google Scholar] [CrossRef]
  58. Kang, T.; Lu, W.; Xu, W.; Anderson, L.; Bacanamwo, M.; Thompson, W.; Chen, Y.E.; Liu, D. MicroRNA-27 (miR-27) Targets Prohibitin and Impairs Adipocyte Differentiation and Mitochondrial Function in Human Adipose-derived Stem Cells. J. Biol. Chem. 2013, 288, 34394–34402. [Google Scholar] [CrossRef] [Green Version]
  59. 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] [CrossRef]
  60. Lee, E.K.; Lee, M.J.; Abdelmohsen, K.; Kim, W.; Kim, M.M.; Srikantan, S.; Martindale, J.L.; Hutchison, E.R.; Kim, H.H.; Marasa, B.S.; et al. miR-130 Suppresses Adipogenesis by Inhibiting Peroxisome Proliferator-Activated Receptor Expression. Mol. Cell. Biol. 2010, 31, 626–638. [Google Scholar] [CrossRef] [Green Version]
  61. Ning, C.; Li, G.; You, L.; Ma, Y.; Jin, L.; Ma, J.; Li, X.; Li, M.; Liu, H. MiR-185 inhibits 3T3-L1 cell differentiation by targeting SREBP-1. Biosci. Biotechnol. Biochem. 2017, 81, 1747–1754. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. 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] [CrossRef] [Green Version]
  63. Chen, H.; Mo, D.; Li, M.; Zhang, Y.; Chen, L.; Zhang, X.; Li, M.; Zhou, X.; Chen, Y. miR-709 inhibits 3T3-L1 cell differentiation by targeting GSK3β of Wnt/β-catenin signaling. Cell. Signal. 2014, 26, 2583–2589. [Google Scholar] [CrossRef]
  64. An, X.; Ma, K.; Zhang, Z.; Zhao, T.; Zhang, X.; Tang, B.; Li, Z. miR-17, miR-21, and miR-143 Enhance Adipogenic Differentiation from Porcine Bone Marrow-Derived Mesenchymal Stem Cells. DNA Cell Biol. 2016, 35, 410–416. [Google Scholar] [CrossRef] [PubMed]
  65. Kang, M.; Yan, L.-M.; Zhang, W.-Y.; Li, Y.-M.; Tang, A.-Z.; Ou, H.-S. Role of microRNA-21 in regulating 3T3-L1 adipocyte differentiation and adiponectin expression. Mol. Biol. Rep. 2013, 40, 5027–5034. [Google Scholar] [CrossRef]
  66. Li, G.; Ning, C.; Ma, Y.; Jin, L.; Tang, Q.; Li, X.; Li, M.; Liu, H. miR-26b Promotes 3T3-L1 Adipocyte Differentiation Through Targeting PTEN. DNA Cell Biol. 2017, 36, 672–681. [Google Scholar] [CrossRef] [PubMed]
  67. Yi, C.; Xie, W.-D.; Li, F.; Lv, Q.; He, J.; Wu, J.; Gu, D.; Xu, N.; Zhang, Y. MiR-143 enhances adipogenic differentiation of 3T3-L1 cells through targeting the coding region of mouse pleiotrophin. FEBS Lett. 2011, 585, 3303–3309. [Google Scholar] [CrossRef] [Green Version]
  68. Shi, C.; Zhang, M.; Tong, M.; Yang, L.; Pang, L.; Chen, L.; Xu, G.; Chi, X.; Hong, Q.; Ni, Y.; et al. miR-148a is Associated with Obesity and Modulates Adipocyte Differentiation of Mesenchymal Stem Cells through Wnt Signaling. Sci. Rep. 2015, 5, srep09930. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Chen, C.; Xiang, H.; Peng, Y.-L.; Peng, J.; Jiang, S.-W. Mature miR-183, negatively regulated by transcription factor GATA3, promotes 3T3-L1 adipogenesis through inhibition of the canonical Wnt/β-catenin signaling pathway by targeting LRP6. Cell. Signal. 2014, 26, 1155–1165. [Google Scholar] [CrossRef] [PubMed]
  70. Qi, R.; Wang, J.; Wang, Q.; Qiu, X.; Yang, F.; Liu, Z.; Huang, J. MicroRNA-425 controls lipogenesis and lipolysis in adipocytes. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 2019, 1864, 744–755. [Google Scholar] [CrossRef]
  71. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef] [PubMed]
  72. Bartelt, A.; Heeren, J. Adipose tissue browning and metabolic health. Nat. Rev. Endocrinol. 2014, 10, 24–36. [Google Scholar] [CrossRef]
  73. Arner, P.; Kulyté, A. MicroRNA regulatory networks in human adipose tissue and obesity. Nat. Rev. Endocrinol. 2015, 11, 276–288. [Google Scholar] [CrossRef] [PubMed]
  74. Xiao, J.; Bai, X.; Liao, L.; Zhou, M.; Peng, J.; Xiang, Q.; Ren, Z.; Wen, H.; Jiang, Z.; Tang, Z.; et al. Hydrogen sulfide inhibits PCSK9 expression through the PI3K/Akt-SREBP-2 signaling pathway to influence lipid metabolism in HepG2 cells. Int. J. Mol. Med. 2019, 43, 2055–2063. [Google Scholar] [CrossRef]
  75. Sun, H.; Liu, X.; Long, S.R.; Wang, T.; Ge, H.; Wang, Y.; Yu, S.; Xue, Y.; Zhang, Y.; Li, X.; et al. Antidiabetic effects of pterostilbene through PI3K/Akt signal pathway in high fat diet and STZ-induced diabetic rats. Eur. J. Pharmacol. 2019, 859, 172526. [Google Scholar] [CrossRef]
  76. Junjvlieke, Z.; Khan, R.; Mei, C.; Cheng, G.; Wang, S.; Raza, S.H.A.; Hong, J.; Wang, X.; Yang, W.; Zan, L. Effect of ELOVL6 on the lipid metabolism of bovine adipocytes. Genomics 2020, 112, 2282–2290. [Google Scholar] [CrossRef]
  77. Choi, J.-W.; Lee, H.-J.; Oh, J.W. Signature Genes in Macrodactyly through Transcriptome Network Analysis Reveal their Association of Lipid Metabolism. Korean J. Phys. Anthr. 2017, 30, 77–85. [Google Scholar] [CrossRef] [Green Version]
  78. Song, Y.-F.; Gao, Y.; Hogstrand, C.; Li, D.-D.; Pan, Y.-X.; Luo, Z. Upstream regulators of apoptosis mediates methionine-induced changes of lipid metabolism. Cell. Signal. 2018, 51, 176–190. [Google Scholar] [CrossRef] [Green Version]
  79. Yang, X.; Fu, Y.; Hu, F.; Luo, X.; Hu, J.; Wang, G. PIK3R3 regulates PPARalpha expression to stimulate fatty acid be-ta-oxidation and decrease hepatosteatosis. Exp. Mol. Med. 2018, 50, e431. [Google Scholar] [CrossRef] [Green Version]
  80. Katona, D.; Rajki, A.; Di Benedetto, G.; Pozzan, T.; Spät, A. Calcium-dependent mitochondrial cAMP production enhances aldosterone secretion. Mol. Cell. Endocrinol. 2015, 412, 196–204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Kissebah, A.; Vydelingum, N.; Tulloch, B.; Hope-Gill, H.; Fraser, T. The Role of Calcium in Insulin Action-I. Purification and properties of enzymes regulating lipolysis in human adipose tissue: Effects of cyclic-AMP and calcium ions. Horm. Metab. Res. 1974, 6, 247–255. [Google Scholar] [CrossRef]
  82. Makdissy, N.; Haddad, K.; Mouawad, C.; Popa, I.; Younsi, M.; Valet, P.; Brunaud, L.; Ziegler, O.; Quilliot, D. Regulation of SREBPs by Sphingomyelin in Adipocytes via a Caveolin and Ras-ERK-MAPK-CREB Signaling Pathway. PLoS ONE 2015, 10, e0133181. [Google Scholar] [CrossRef] [Green Version]
  83. Al Dow, M.; Silveira, M.A.D.; Poliquin, A.; Tribouillard, L.; Fournier, É.; Trébaol, E.; Secco, B.; Villot, R.; Tremblay, F.; Bilodeau, S.; et al. Control of adipogenic commitment by a STAT3-VSTM2A axis. Am. J. Physiol. Metab. 2021, 320, E259–E269. [Google Scholar] [CrossRef]
  84. Gao, Z.; Ding, R.; Zhai, X.; Wang, Y.; Chen, Y.; Yang, C.-X.; Du, Z.-Q. Common Gene Modules Identified for Chicken Adiposity by Network Construction and Comparison. Front. Genet. 2020, 11, 537. [Google Scholar] [CrossRef]
  85. Ross, S.E.; Hemati, N.; Longo, K.A.; Bennett, C.N.; Lucas, P.C.; Erickson, R.L.; MacDougald, O.A. Inhibition of adipogen-esis by Wnt signaling. Science 2000, 289, 950–953. [Google Scholar] [CrossRef] [PubMed]
  86. Liu, R.; Liu, X.; Bai, X.; Xiao, C.; Dong, Y. Different expression of lipid metabolism-related genes in Shandong black cattle and Luxi cattle based on transcriptome analysis. Sci. Rep. 2020, 10, 1–14. [Google Scholar] [CrossRef] [PubMed]
  87. Zhao, W.; Pu, M.; Shen, S.; Yin, F. Geniposide improves insulin resistance through AMPK-mediated Txnip protein deg-radation in 3T3-L1 adipocytes. Acta Biochim. Biophys. Sin. 2021, 53, 160–169. [Google Scholar] [CrossRef] [PubMed]
  88. Wang, G.; Wu, B.; Xu, W.; Jin, X.; Wang, K.; Wang, H. The Inhibitory Effects of Juglanin on Adipogenesis in 3T3-L1 Adipocytes. Drug Des. Dev. Ther. 2020, 14, 5349–5357. [Google Scholar] [CrossRef]
  89. Zhang, Z.; Meng, Y.; Gao, F.; Xiao, Y.; Zheng, Y.; Wang, H.-Q.; Gao, Y.; Jiang, H.; Yuan, B.; Zhang, J.-B. TGF-β1-Mediated FDNCR1 Regulates Porcine Preadipocyte Differentiation via the TGF-β Signaling Pathway. Animals 2020, 10, 1399. [Google Scholar] [CrossRef]
  90. Taylor, B.; Shah, A.; Bielczyk-Maczyńska, E. TGF-β is insufficient to induce adipocyte state loss without concurrent PPARγ downregulation. Sci. Rep. 2020, 10, 1–13. [Google Scholar] [CrossRef]
  91. Justine, E.; Yu, J.; Han, S.-Y.; Wolfson, B.; Zhou, Q. The role of endothelial lipase in lipid metabolism, inflammation, and cancer. Histol. Histopathol. 2017, 33, 1–10. [Google Scholar] [CrossRef]
  92. Hong, C.; Deng, R.; Wang, P.; Lu, X.; Zhao, X.; Wang, X.; Cai, R.; Lin, J. LIPG: An inflammation and cancer modulator. Cancer Gene Ther. 2021, 28, 27–32. [Google Scholar] [CrossRef]
  93. Conti, I.; Varano, G.; Simioni, C.; LaFace, I.; Milani, D.; Rimondi, E.; Neri, L.M. miRNAs as Influencers of Cell–Cell Communication in Tumor Microenvironment. Cells 2020, 9, 220. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Imagawa, M. Molecular Mechanisms of Early-stage Adipocyte Differentiation and Multi-functional Roles of Newly Isolated Adipogenic Factors. Yakugaku Zasshi J. Pharm. Soc. Jpn. 2016, 136, 649–658. [Google Scholar] [CrossRef] [Green Version]
  95. Reichert, M.; Eick, D. Analysis of cell cycle arrest in adipocyte differentiation. Oncogene 1999, 18, 459–466. [Google Scholar] [CrossRef] [Green Version]
  96. Lin, Y.; Wen-Jie, Z.; Chang-Qing, L.; Sheng-Xiang, A.; Yue, Z. mir-22-3p/KLF6/MMP14 axis in fibro-adipogenic progenitors regulates fatty infiltration in muscle degeneration. FASEB J. Off. Publ. Fed. Am. Soc. Exp. Biol. 2020, 34, 12691–12701. [Google Scholar] [CrossRef]
  97. Joseph, R.; Poschmann, J.; Sukarieh, R.; Too, P.G.; Julien, S.G.; Xu, F.; Teh, A.L.; Holbrook, J.D.; Ng, K.L.; Chong, Y.S.; et al. ACSL1 Is Associated With Fetal Programming of Insulin Sensitivity and Cellular Lipid Content. Mol. Endocrinol. 2015, 29, 909–920. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  98. Liu, L.; Wang, Y.; Liang, X.; Wu, X.; Liu, J.; Yang, S.; Tao, C.; Zhang, J.; Tian, J.; Zhao, J.; et al. Stearoyl-CoA Desaturase Is Essential for Porcine Adipocyte Differentiation. Int. J. Mol. Sci. 2020, 21, 2446. [Google Scholar] [CrossRef] [Green Version]
  99. Tan, Z.; Du, J.; Shen, L.; Liu, C.; Ma, J.; Bai, L.; Jiang, Y.; Tang, G.; Li, M.; Li, X.; et al. miR-199a-3p affects adipocytes differentiation and fatty acid composition through targeting SCD. Biochem. Biophys. Res. Commun. 2017, 492, 82–88. [Google Scholar] [CrossRef]
  100. Park, J.; Yoon, Y.-S.; Han, H.-S.; Kim, Y.-H.; Ogawa, Y.; Park, K.-G.; Lee, C.-H.; Kim, S.-T.; Koo, S.-H. SIK2 Is Critical in the Regulation of Lipid Homeostasis and Adipogenesis In Vivo. Diabetes 2014, 63, 3659–3673. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Wang, X.; Li, A.; Raza, S.H.A.; Liang, C.; Zhang, S.; Mei, C.; Yang, W.; Zan, L. Transcription Factors ZEB1 and CREB Promote the Transcription of Bovine ABHD5 Gene. DNA Cell Biol. 2021, 40, 219–230. [Google Scholar] [CrossRef] [PubMed]
Animals 11 01207 g001 550
Figure 1. Transfection and interference efficiencies of miR-149-5p: (A) qPCR results showing expression levels of bta-miR-149-5p in adipocytes transfected with bta-miR-149-5p mimic and its negative control. (B–D) Transfection efficiency was checked through FAM/labeled miR-NC in bovine adipocytes; after 6-, 12-, and 24-h intervals, the transfection efficiency was measured through an Olympus I × 71 microscope (OLYMPUS) in bovine adipocytes. U6 was used as a house-keeping gene, and the relative expression level in adipocytes transfected with bta-miR-149-5p mimics was normalized with adipocytes transfected with mimic NC. The values presented with mean ± SEM, n = 3, and ** p < 0.01 were considered statistically significant.
Figure 1. Transfection and interference efficiencies of miR-149-5p: (A) qPCR results showing expression levels of bta-miR-149-5p in adipocytes transfected with bta-miR-149-5p mimic and its negative control. (B–D) Transfection efficiency was checked through FAM/labeled miR-NC in bovine adipocytes; after 6-, 12-, and 24-h intervals, the transfection efficiency was measured through an Olympus I × 71 microscope (OLYMPUS) in bovine adipocytes. U6 was used as a house-keeping gene, and the relative expression level in adipocytes transfected with bta-miR-149-5p mimics was normalized with adipocytes transfected with mimic NC. The values presented with mean ± SEM, n = 3, and ** p < 0.01 were considered statistically significant.
Animals 11 01207 g001
Animals 11 01207 g002 550
Figure 2. Quality and integrity of the RNA samples and data analysis. (A) RNA quality and integrity of adipocytes transfected with miR-149-5p mimic and mimic NC. (B) The reliability of the sequenced data collected from adipocytes transfected with miR-149-5p mimic and mimic NC. (C) The filtering of clean reads: the blue color represents the high-quality reads after removal of the adapters and low-quality bases. (D) The clean reads were aligned and represented with different colors, which exhibited different types of genomic regions: intergenic region (blue), intronic region (green), and exonic region (red). (E) The distribution of FPKM (fragment per kilobase of transcript per million mapped reads) in each sample of the 2 groups (adipocytes transfected with miR-149-5p mimic and mimic NC). The different colors represent different samples processed as mentioned below each sample.
Figure 2. Quality and integrity of the RNA samples and data analysis. (A) RNA quality and integrity of adipocytes transfected with miR-149-5p mimic and mimic NC. (B) The reliability of the sequenced data collected from adipocytes transfected with miR-149-5p mimic and mimic NC. (C) The filtering of clean reads: the blue color represents the high-quality reads after removal of the adapters and low-quality bases. (D) The clean reads were aligned and represented with different colors, which exhibited different types of genomic regions: intergenic region (blue), intronic region (green), and exonic region (red). (E) The distribution of FPKM (fragment per kilobase of transcript per million mapped reads) in each sample of the 2 groups (adipocytes transfected with miR-149-5p mimic and mimic NC). The different colors represent different samples processed as mentioned below each sample.
Animals 11 01207 g002
Animals 11 01207 g003 550
Figure 3. Bta-miR-149-5p represented bovine adipocyte differentiation. The photos were captured through an Olympus 1 × 71 microscope (Olympus, Tokyo, Japan). (A) The pictures of BODIPY-stained LDs in bovine adipocytes are represented in green, and the nuclei stained with 4, 6-diamidino-2-phenylindole (DAPI) are represented in blue. (B) The histogram pixel intensity values of adipocytes transfected with bta-miR-149-5p mimics or mimic NC. The values are presented with mean ± SEM, n = 3; ** p < 0.01 indicates statistical significance.
Figure 3. Bta-miR-149-5p represented bovine adipocyte differentiation. The photos were captured through an Olympus 1 × 71 microscope (Olympus, Tokyo, Japan). (A) The pictures of BODIPY-stained LDs in bovine adipocytes are represented in green, and the nuclei stained with 4, 6-diamidino-2-phenylindole (DAPI) are represented in blue. (B) The histogram pixel intensity values of adipocytes transfected with bta-miR-149-5p mimics or mimic NC. The values are presented with mean ± SEM, n = 3; ** p < 0.01 indicates statistical significance.
Animals 11 01207 g003
Animals 11 01207 g004 550
Figure 4. Quantification of DEGs in differentiated bovine adipocytes transfected with bta-miR-149-5p mimic and mimic NC. (A) The heatmap gene expression chart for the hierarchical cluster analysis of DEGs between 2 different groups of adipocytes transfected with bta-miR-149-5p mimic (n = 3) and mimic NC (n = 3). (B) The bar graph represents the number of DEGs between the 2 groups of differentiated bovine adipocytes transfected with bta-miR-149-5p mimic (n = 3) and mimic NC (n = 3). The red shows the number of upregulated genes and blue represents the number of downregulated genes.
Figure 4. Quantification of DEGs in differentiated bovine adipocytes transfected with bta-miR-149-5p mimic and mimic NC. (A) The heatmap gene expression chart for the hierarchical cluster analysis of DEGs between 2 different groups of adipocytes transfected with bta-miR-149-5p mimic (n = 3) and mimic NC (n = 3). (B) The bar graph represents the number of DEGs between the 2 groups of differentiated bovine adipocytes transfected with bta-miR-149-5p mimic (n = 3) and mimic NC (n = 3). The red shows the number of upregulated genes and blue represents the number of downregulated genes.
Animals 11 01207 g004
Animals 11 01207 g005 550
Figure 5. The Gene Ontology (GO) classification of DEGs in 2 groups of bovine adipocytes, transfected with bta-miR-149-5p (n = 3) mimics and mimics NC (n = 3). (A) GO terms of all unigenes identified in the transcriptome of differentiated bovine adipocytes within bta-miR-149-5p mimics (n = 3) and mimics NC (n = 3). The bars in each GO term were represented with green (downregulated genes) and red bars (upregulated unigenes). (B) The z-score within each category of GO term is represented with blue (biological process), green (cellular component), and yellow (molecular function).
Figure 5. The Gene Ontology (GO) classification of DEGs in 2 groups of bovine adipocytes, transfected with bta-miR-149-5p (n = 3) mimics and mimics NC (n = 3). (A) GO terms of all unigenes identified in the transcriptome of differentiated bovine adipocytes within bta-miR-149-5p mimics (n = 3) and mimics NC (n = 3). The bars in each GO term were represented with green (downregulated genes) and red bars (upregulated unigenes). (B) The z-score within each category of GO term is represented with blue (biological process), green (cellular component), and yellow (molecular function).
Animals 11 01207 g005
Animals 11 01207 g006 550
Figure 6. The Gene Ontology (GO) classification of DEGs in 2 groups of bovine adipocytes, transfected with bta-miR-149-5p (n = 3) mimics and mimics NC (n = 3). Each category of GO term was represented with different colors: yellow reflects biological process, green shows cellular component, and blue depicts molecular function.
Figure 6. The Gene Ontology (GO) classification of DEGs in 2 groups of bovine adipocytes, transfected with bta-miR-149-5p (n = 3) mimics and mimics NC (n = 3). Each category of GO term was represented with different colors: yellow reflects biological process, green shows cellular component, and blue depicts molecular function.
Animals 11 01207 g006
Animals 11 01207 g007 550
Figure 7. The KEGG pathway analysis of DEGs in 2 groups of bovine adipocytes transfected with bta-miR-149-5p (n = 3) mimics and mimics NC (n = 3). The KO enrichment difference bubble chart: the ordinate is -log10 (Qvalue), the abscissa is the z-score value (the difference between the number of upregulated genes and the number of downregulated genes accounts for the proportion of the total differential genes), and the yellow line represents Qvalue = 0.05 threshold. On the right is the list of pathways with the top 20 Q values. Different colors represent different A classes.
Figure 7. The KEGG pathway analysis of DEGs in 2 groups of bovine adipocytes transfected with bta-miR-149-5p (n = 3) mimics and mimics NC (n = 3). The KO enrichment difference bubble chart: the ordinate is -log10 (Qvalue), the abscissa is the z-score value (the difference between the number of upregulated genes and the number of downregulated genes accounts for the proportion of the total differential genes), and the yellow line represents Qvalue = 0.05 threshold. On the right is the list of pathways with the top 20 Q values. Different colors represent different A classes.
Animals 11 01207 g007
Animals 11 01207 g008 550
Figure 8. KO enrichment circle diagram of the KEGG_A_class: The first circle: the first 20 pathways of enrichment; outside the circle is the coordinate ruler of the number of genes. Different colors represent different A classes. The second circle: the number of pathways in the background genes and the Q value. The more genes, the longer the bar, the smaller the Q value, the redder the color. The third circle: the bar graph of upregulated gene ratio, wherein dark purple represents the upregulated gene ratio and light purple represents the downregulated gene ratio; the specific value is shown below. The four circles: Rich factor value of each pathway (the number of differences in the pathway divided by all the numbers); background grid lines each grid represent 0.1.
Figure 8. KO enrichment circle diagram of the KEGG_A_class: The first circle: the first 20 pathways of enrichment; outside the circle is the coordinate ruler of the number of genes. Different colors represent different A classes. The second circle: the number of pathways in the background genes and the Q value. The more genes, the longer the bar, the smaller the Q value, the redder the color. The third circle: the bar graph of upregulated gene ratio, wherein dark purple represents the upregulated gene ratio and light purple represents the downregulated gene ratio; the specific value is shown below. The four circles: Rich factor value of each pathway (the number of differences in the pathway divided by all the numbers); background grid lines each grid represent 0.1.
Animals 11 01207 g008
Animals 11 01207 g009 550
Figure 9. The distribution of up- and downregulated DEG enrichment in different KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. The pathway nodes are represented by pathways (yellow), downregulated genes (green), and upregulated genes (red), which are shown as circles.
Figure 9. The distribution of up- and downregulated DEG enrichment in different KEGG (Kyoto Encyclopedia of Genes and Genomes) pathways. The pathway nodes are represented by pathways (yellow), downregulated genes (green), and upregulated genes (red), which are shown as circles.
Animals 11 01207 g009
Animals 11 01207 g010 550
Figure 10. Top 20 reacrome enrichment pathways. (A) Reactome enrichment bar graph: the top 20 reactome channels with the smallest Q value to plot; the ordinate is the reactome channel, and the abscissa is the percentage of the number of reactome channels in all the DEGs. The darker the color, the smaller the Q value. The value on the column is the number of reactome channels and the Q value. (B) Reactome enrichment circle diagram: the first circle: the top 20 reactome pathways of the enrichment; outside the circle is the coordinate ruler of the number of genes. The second circle: the number of the reactome pathway in the background gene and the Q value. The more the genes, the longer the bars, the smaller the Q value, and the redder the color. The third circle: the bar graph of upregulated gene ratio, wherein dark purple represents the upregulated gene ratio and light purple represents the downregulated gene ratio; the specific value is shown below. The fourth circle: each reactome pathway Rich factor value (the number of differences in the reactome channel divided by all the numbers); in terms of background grid lines, each grid represents 0.1.
Figure 10. Top 20 reacrome enrichment pathways. (A) Reactome enrichment bar graph: the top 20 reactome channels with the smallest Q value to plot; the ordinate is the reactome channel, and the abscissa is the percentage of the number of reactome channels in all the DEGs. The darker the color, the smaller the Q value. The value on the column is the number of reactome channels and the Q value. (B) Reactome enrichment circle diagram: the first circle: the top 20 reactome pathways of the enrichment; outside the circle is the coordinate ruler of the number of genes. The second circle: the number of the reactome pathway in the background gene and the Q value. The more the genes, the longer the bars, the smaller the Q value, and the redder the color. The third circle: the bar graph of upregulated gene ratio, wherein dark purple represents the upregulated gene ratio and light purple represents the downregulated gene ratio; the specific value is shown below. The fourth circle: each reactome pathway Rich factor value (the number of differences in the reactome channel divided by all the numbers); in terms of background grid lines, each grid represents 0.1.
Animals 11 01207 g010
Animals 11 01207 g011 550
Figure 11. Validation of RNA-Seq data through qRT-PCR. Expression (A–I) relative mRNA levels of CCND2, KLF12, LIPG, KLF6, ACSL1, Cdk2, SCD, SIK2, and ZEB1 in bovine adipocytes transfected with bta-miR-149-5p mimics or mimics negative control (NC). The values exhibit the mean ± SEM (n = 3), and statistical significance is shown by * p < 0.05 and ** p < 0.01.
Figure 11. Validation of RNA-Seq data through qRT-PCR. Expression (A–I) relative mRNA levels of CCND2, KLF12, LIPG, KLF6, ACSL1, Cdk2, SCD, SIK2, and ZEB1 in bovine adipocytes transfected with bta-miR-149-5p mimics or mimics negative control (NC). The values exhibit the mean ± SEM (n = 3), and statistical significance is shown by * p < 0.05 and ** p < 0.01.
Animals 11 01207 g011
Table
Table 1. Description of KEGG pathways classes altered by DEGs in bovine adipocytes transfected with bta-miR-149-5p mimics or mimics NC.
Table 1. Description of KEGG pathways classes altered by DEGs in bovine adipocytes transfected with bta-miR-149-5p mimics or mimics NC.
KEGG_A_ClassKEGG_B_ClassPathwayDEGsPvalueQvaluePathway ID
Environmental information processingSignaling molecules and interactionCytokine–cytokine receptor interaction90.0003930.031419ko04060
Signal transductionTNF signaling pathway40.0100950.164541ko04668
Signal transductionHippo signaling pathway—multiple species20.0154180.179513ko04392
Signal transductionRap1 signaling pathway50.0137270.172113ko04015
Organismal systemsImmune systemComplement and coagulation cascades50.0004610.031419ko04610
Immune systemIL-17 signaling pathway40.0049230.143156ko04657
Endocrine systemOvarian steroidogenesis30.0058530.143156ko04913
Endocrine systemRegulation of lipolysis in adipocyte30.006770.143156ko04923
Immune systemChemokine signaling pathway50.0070260.143156ko04062
Endocrine systemInsulin secretion30.017680.192128ko04911
Excretory systemAldosterone-regulated sodium reabsorption20.0257110.24652ko04960
Human diseasesImmune diseasesRheumatoid arthritis50.0005780.031419ko05323
Infectious diseasesPertussis40.0018390.074941ko05133
Infectious diseasesLegionellosis30.0088470.160228ko05134
Infectious diseasesEpithelial cell signaling in Helicobacter pylori infection30.0125950.171085ko05120
CancersMelanoma30.0125950.171085ko05218
Cardiovascular diseasesDilated cardiomyopathy (DCM)30.0224610.228825ko05414
Infectious diseasesSalmonella infection30.0278810.252474ko05132
CancersBladder cancer20.0323440.26341ko05219
Cellular processesCellular community—eukaryotesTight junction40.032526 0.26341 ko04530
Cell growth and deathOocyte meiosis30.0451150.301575ko04114
Cellular community—eukaryotesFocal adhesion40.0475410.301575ko04510
KEGG: Kyoto Encyclopedia of Genes and Genomes. DEG: Deferentially Expressed Genes.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Guo, H.; Khan, R.; Abbas Raza, S.H.; Suhail, S.M.; Khan, H.; Khan, S.B.; Abd El-Aziz, A.H.; Zan, L. RNA-Seq Reveals Function of Bta-miR-149-5p in the Regulation of Bovine Adipocyte Differentiation. Animals 2021, 11, 1207. https://doi.org/10.3390/ani11051207

AMA Style

Guo H, Khan R, Abbas Raza SH, Suhail SM, Khan H, Khan SB, Abd El-Aziz AH, Zan L. RNA-Seq Reveals Function of Bta-miR-149-5p in the Regulation of Bovine Adipocyte Differentiation. Animals. 2021; 11(5):1207. https://doi.org/10.3390/ani11051207

Chicago/Turabian Style

Guo, Hongfang, Rajwali Khan, Sayed Haidar Abbas Raza, Syed Muhammad Suhail, Hamayun Khan, Sher Bahadar Khan, Ayman Hassan Abd El-Aziz, and Linsen Zan. 2021. "RNA-Seq Reveals Function of Bta-miR-149-5p in the Regulation of Bovine Adipocyte Differentiation" Animals 11, no. 5: 1207. https://doi.org/10.3390/ani11051207

APA Style

Guo, H., Khan, R., Abbas Raza, S. H., Suhail, S. M., Khan, H., Khan, S. B., Abd El-Aziz, A. H., & Zan, L. (2021). RNA-Seq Reveals Function of Bta-miR-149-5p in the Regulation of Bovine Adipocyte Differentiation. Animals, 11(5), 1207. https://doi.org/10.3390/ani11051207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop