Effect of miR-149-5p on intramuscular fat deposition in pigs based on metabolomics and transcriptomics
BMC Genomics volume 24, Article number: 293 (2023)
As one of the important traits in pig production, meat quality has important research significance and value. Intramuscular fat (IMF) content is one of the most important factors affecting pork quality. Many experimental studies have shown that IMF content is closely related to the flavor, tenderness, and juiciness of pork. Therefore, it is of great significance to study the mechanism of porcine IMF deposition. Previous research indicated that miR-149-5p promoted the proliferation of porcine intramuscular (IM) preadipocytes and decreased their ability to differentiate, albeit the exact mechanism of action is unknown. In vitro, foreign pigs showed increased miR-149-5p expression and reduced fat deposition when compared to Queshan Black pigs. This study conducted metabolomics and transcriptomics analyses of porcine IM preadipocytes overexpressing miR-149-5p to verify their effects on lipid formation. According to metabolomics analysis, the overexpression of miR-149-5p has significantly altered the lipid, organic acid, and organic oxygen metabolites of porcine IM preadipocytes. Specially speaking, it has changed 115 metabolites, including 105 up-regulated and 10 down-regulated ones, as well as the composition of lipid, organic acid, and organic oxygen metabolism-related metabolites. RNA-seq analysis showed that overexpression of miR-149-5p significantly altered 857 genes, of which 442 were up-regulated, and 415 were down-regulated, with enrichment to MAPK, IL-17, PI3K-Akt, and ErbB signaling pathways. We found that overexpression of miR-149-5p inhibited adipogenic differentiation by changing cAMP signaling pathway in porcine IM preadipocytes. In addition, the overexpression of miR-149-5p may affect the transport of Cu2+ by targeting ATP7A and inhibiting adipogenic differentiation. These findings elucidate the regulatory function of miR-149-5p in porcine IM preadipocytes, which may be a key target for controlling pork quality.
Pork is an important source of animal protein, energy, and iron for humans as well as the economic backbone of animal husbandry . An increasing number of studies show that pork currently accounts for more than 40% of all meat products consumed worldwide [2, 3]. Because of improving living standards, people now have more stringent demands for pork quality. Thus, improving the quality of pork is the top priority for the pig industry. In essence, the growth of fat and muscle is related to the breed of pig. Pork quality is closely related to the amount of intramuscular fat (IMF) [4,5,6]. Many experimental studies have demonstrated that IMF content is closely related to the flavor, tenderness, and juiciness of pork [7,8,9]. The fat breakdown has a considerable impact on the flavor of meat products. During processing, it can decompose into volatile compounds such as aliphatic hydrocarbons, aldehydes, ketones, alcohols, carboxylic acids, and esters. These chemicals can directly affect the flavor of meat products through volatilization as well as indirectly through mixing these chemicals with other compounds to produce other flavor substances . IMF is mainly distributed in perimysium and endomysium. With the rise in IMF content, the number of muscle fibers per unit area will decrease, thus affecting the tenderness of meat products . Pork quality is a complex trait, which is influenced by genetic factors, environment, management, feeding and slaughter conditions . Therefore, understanding the mechanisms underlying animal fat deposition is essential for improving meat quality.
Preadipocytes develop into white adipocytes, which serve as both humans’ and animals’ main fuel for storing energy. Adipose tissue, one of the largest and busiest organs in both humans and animals, is essential for maintaining energy balance . The two primary mechanisms that contribute to the development of adipose tissue are hypertrophy (differentiation) and adipocyte proliferation [14, 15]. The development of preadipocytes is controlled by several important genes, including lipin 1 (LPIN1) , fatty acid binding protein 4 (FABP4) , and peroxisome proliferator-activated receptor gamma (PPARG) . The primary source of IMF is fibrous adipogenic progenitor cells (FAPs), which can differentiate into adipogenesis and fibrosis. More than 90% of FAPs in mice can differentiate into different forms of fat, and when muscles are damaged, this number decreases . However, only 30% of human FAPs can form fat . Therefore, learning more about the differentiation of porcine IM preadipocytes can help us better understand animal development and eventually, improve meat quality.
Approximately 22-nucleotide non-coding single stranded RNA molecules known as microRNAs (miRNAs) regulate the growth of adipose tissue, lipid metabolism, and adipogenesis in many animals [21,22,23,24,25,26]. Mounting data suggest that miRNA is crucial to the development of porcine adipocytes [27,28,29,30]. The miRNA typically binds to the 3’ untranslated region (3’UTR) sequence of its mRNA target to degrade or inhibit the process of translation , but it can also bind to the unconventional target of the 5’ untranslated region (5’UTR) and CDS to prevent the process of translation . Numerous studies have shown that several miRNAs, such as miR-122 , miR-335 , miR-27 , miR-370 , and miR-30 , act as important regulators of lipid metabolism and fat deposition. A novel miRNA, which is called miR-149, has been found in the phosphatidylinositol glypican 1 (GPC1) gene’s first intron. Its precursor can produce the miR-149-5p and miR-149-3p active chains, which interact with target genes to perform their biological effects. Mohamed’s study found that miR-149 inhibited PARP-2 expression through the SIRT-1/PGC-1 network and promoted mitochondrial biogenesis in skeletal muscle . By targeting FTO, miR-149-3p controls the shift of bone marrow mesenchymal stem cells from adipogenesis to osteogenesis . MiR-149 controls non-alcoholic fatty liver disease by targeting FGF-21 ; it also inhibits the ATF6 pathway to reduce inflammation and apoptosis in nonalcoholic fatty liver disease brought on by endoplasmic reticulum stress . According to Liu et al., miR-149-5p reduces the differentiation of porcine IM preadipocytes and 3T3-L1 cells while promoting their growth . However, it is still unclear how miR-149-5p influences the development of porcine IM preadipocyte. In this experiment, we conducted metabolomics and transcriptomics sequencing analysis on the overexpression of miR-149-5p.
In previous studies, we identified the role of miR-149-5p in porcine IM preadipocytes and 3T3-L1 cells. The results showed that overexpression of miR-149-5p significantly inhibited the adipogenic differentiation of porcine IM preadipocytes and 3T3-L1 cells while interference with miR-149-5p significantly promoted its adipogenic differentiation. However, the regulatory and metabolic roles of miR-149-5p in porcine IM preadipocytes remain unclear. To this end, this study aims to investigate the metabolism and transcription mechanism of miR-149-5p in porcine IM preadipocytes, the regulatory role of miR-149-5p in porcine IM preadipocytes differentiation, as well as the effect and mechanism of miR-149-5p on pork quality.
Breed-specific differences in the expression of miR-149-5p
The study compared the slaughter characteristics of (Duroc×Landrace×Yorkshire) DLY and Queshan Black pigs. Vitro studies showed that compared with DLY, Queshan Black pigs had lower body weight at the same age (P < 0.05, Fig. 1A), higher marbling score (P < 0.05 Fig. 1B), higher IMF content (P < 0.05, Fig. 1C), higher backfat thickness (P < 0.01, Fig. 1D), higher fat ratio (P < 0.01, Fig. 1E), and lower lean ratio (P < 0.01, Fig. 1F). The fat deposition capacity of Queshan Black pigs was significantly higher than that of DLY, and the lean ratio was significantly lower than that of DLY (P < 0.01). The tissue expression profile of Queshan Black pigs showed that miR-149-5p had the highest expression in subcutaneous fat, followed by longissimus dorsi (Fig. 2A). Compared with Yorkshire, the expression of miR-149-5p in the liver, longissimus dorsi, leg muscle, and subcutaneous fat of Queshan Black pigs decreased significantly (P < 0.05 Fig. 2B). In addition, expressions of adipogenic differentiation marker genes (CEBPA, PPARG, FABP4, and LIPE) of the longissimus dorsi in Queshan Black pigs were significantly higher than those in Yorkshire (P < 0.05 Fig. 2C).
MiR-149-5p overexpression altered total metabolites in porcine IM preadipocyte
To investigate the effect of miR-149-5p on metabolism, we transfected porcine IM preadipocytes with miR-149-5p mimics. A total of 115 differential metabolites were found, of which 10 were down-regulated and 105 were up-regulated (Fig. 3A). KEGG results showed that the differential metabolites were highly enriched to glucosinolate biosynthesis, phenylpropanoid biosynthesis, methane metabolism, and other signal pathways (Fig. 3B). We categorized 115 metabolites according to data from the Human Metabolome Database and LIPID MAPS to comprehend the molecular type of metabolites. The amount and composition of metabolites associated with organic acids, their derivatives, and lipid metabolism showed significant alterations in cells overexpressing miR-149-5p (Fig. 3C). Among the metabolites, organic acid and derivatives metabolites (56.5%), lipid and lipid-like molecules metabolites (17.4%), nucleosides, nucleotides, and analogues metabolites (0.9%), organic nitrogen compounds metabolites (3.5%), organic oxygen compounds metabolites (3.5%), organoheterocyclic compounds metabolites (4.3%), and benzenoids metabolites (1.7%) in IM preadipocytes underwent significant change. 17 up-regulated and 3 down-regulated lipids were among the significantly changed lipid and lipid-like molecules metabolites (Fig. 3D), which were then divided into lipid subclasses (Fig. 3E). The majority of the significantly changed lipid metabolites in miR-149-5p overexpressed IM preadipocytes were fatty acids (FAS, 55.0%), glycerophospholipid (GP, 25.0%), prenol lipid (PL, 10.0%), and sphingolipid (ST, 10.0%). Significantly changed lipid metabolites were shown on the heatmap (Fig. 3F).
Overexpression of miR-149-5p changed the expression of genes involved in metabolic pathways
We used RNA-seq to map the transcriptional changes in adipocytes as a result of miR-149-5p overexpression to investigate how the adipocyte metabolites were changed. We found 857 differentially expressed genes, among which 442 were up-regulated, and 415 were down-regulated upon miR-149-5p overexpression (Fig. 4A). Significantly changed genes were displayed in the heatmap (Fig. 4B). Up-regulated genes in the overexpression of miR-149-5p gathered in PPAR and IL-17 signaling pathways, according to KEGG pathway analysis (Fig. 4C). In addition, down-regulated in the overexpression of miR-149-5p gathered in IL-17, MAPK, and ErbB signaling pathways (Fig. 4D). The analysis of differentially up-regulated genes revealed a significant increase in the process of peptide biosynthetic and amide biosynthetic (Fig. 4E). The analysis of differentially down-regulated genes revealed a significant decrease in regulation of protein kinase activity, protein phosphorylation, and interleukin-6 production (Fig. 4F).
Overexpression of miR-149-5p affected lipid metabolism of porcine IM preadipocytes
To determine the signal pathway involved in lipid metabolism, we used bioinformatics analysis methods to select metabolite-related genes from the KEGG pathway (Fig. 5A–D). Heatmap demonstrates that the miR-149-5p overexpression dramatically altered the expression of the followings: MAPK signal pathway-related genes, including interleukin 1 alpha (IL1A), ATPase copper transporting alpha (ATP7A), and epiregulin (EREG) (Fig. 5A); IL-17 signal pathway-related genes, including prostaglandin-endoperoxide synthase 2 (PTGS2), mitogen-activated protein kinase 9 (MAPK9), interleukin 17 receptor C (IL17RC), and CCAAT enhancer binding protein beta (CEBPB) (Fig. 5B); PI3K-Akt signal pathway-related genes, including fms related receptor tyrosine kinase 1 (FLT1), hepatocyte growth factor (HGF), and colony stimulating factor 1 (CSF1) (Fig. 5C); and ErbB signal pathway-related genes, including MYC proto-oncogene, bHLH transcription factor (MYC), cbl proto-oncogene B (CBLB), and amphiregulin (AREG) (Fig. 5D). We performed protein-protein interaction analysis on the genes in these four pathways and found that 40 genes had significant protein-protein interaction, especially (chemokine (C-C motif) ligand 2) CCL2, (mitogen-activated protein kinase 9) MAPK9, (matrix metallopeptidase 3) MMP3, and HGF (Fig. 5E). We constructed a cAMP signaling pathway metabolic network to describe how differentially expressed genes (DEGs) mediate the regulation of miR-149-5p (Fig. 6A). In addition, we performed a heatmap clustering analysis of cAMP signaling pathway-related genes (Fig. 6B).
MiR-149-5p influenced lipid formation by targeting ATP7A
According to the results of RNA-seq, once miR-149-5p was overexpressed, the expression of ATP7A drastically dropped, and since RNA22 v2 suggested that they had binding sites (Fig. 7A), it was assumed that ATP7A was miR-149-5p’s target. To prove that, dual-luciferase assay detection was done on wild-type and mutant vectors. The findings showed that miR-149-5p significantly decrease the activity of the wild-type vector without having impact on the activity of the mutant vector when compared to the control group (P < 0.01, Fig. 7B). ATP7A is a Cu2+ transporter, which may affect the formation of IMF by affecting the transport of Cu2+. MiR-149-5p expression was considerably enhanced upon transfection (P < 0.01, Fig. 7C). MiR-149-5p mimics dramatically reduced the expression of ATP7A in porcine IM preadipocytes (P < 0.01), and the expressions of Cu2+ transporters CTR1 and CTR2 but insignificantly (Fig. 7D). To detect the targeting of miR-149-5p and ATP7A, this study interfered with miR-149-5p (Fig. 7E). MiR-149-5p inhibitors significantly increased the expression of ATP7A and CTR2 in porcine IM preadipocytes (P < 0.05), and the expression of Cu2+ transporter CTR1 but insignificantly (Fig. 7F). Meanwhile, compared with foreign pigs, Queshan Black pigs showed significantly increased expressions of ATP7A, CTR1, and CTR2 in the longissimus dorsi, leg muscles, and subcutaneous fat (Figs. 7G–I).
MiR-149-5p regulates energy metabolism in addition to genes associated with fat storage and obesity. We summarized the metabolic and transcriptional characteristics of porcine IM preadipocyte overexpressed with miR-149-5p. Furthermore, this experiment found that miR-149-5p may affect the transport of Cu2+ by targeting ATP7A, further affecting the formation of IMF. Phospholipids are significant compounds that impact the taste of meat products, making up IMF together with triglycerides. They are mainly composed of C:16 and C:18 fatty acids . Pork quality is determined by the kind and percentage of fatty acids, and this can serve as a benchmark when judging the quality of other meats . Monounsaturated and saturated fatty acids can improve the quality of meat, whereas polyunsaturated fatty acids have the opposite effect . IMF is also closely related to the tenderness of meat products. It not only destroys the cross-linked structure between muscle fiber bundles but also weakens the strength of connective tissue, thus reducing the shearing force of muscle and improving the tenderness of meat products .
Having a long history and a sizable population in Queshan County, Henan Province, the Queshan Black pig is a superior local breed of pigs. They are characterized by strong fecundity, strong adaptability, good meat quality, and stable genetic performance [47,48,49]. In vitro, studies showed that the fat deposition ability of Queshan Black pigs was significantly higher than that of foreign pigs, and the expression of miR-149-5p in subcutaneous fat was significantly decreased, indicating that the expression of miR-149-5p was possibly related to fat accumulation. In addition, compared with foreign pigs, miR-149-5p in Queshan Black pigs was significantly decreased in liver, longissimus dorsi, and leg muscles, indicating that miR-149-5p expression was associated with possible IMF deposition. The liver is an important metabolic organ in animals, which is closely related to lipid deposition. In previous studies, miR-149-5p overexpression inhibited the differentiation of porcine IM preadipocytes and the expression of fat formation-related genes. At the same time, many studies found that miR-149 is strongly correlated to fat production. For instance, miR-149 induces subcutaneous fat change into visceral fat and inhibits PRDM16 . In addition, the decrease of miR-149-3p significantly increased the expression of PRDM16 in SAT, improved the systemic insulin sensitivity of high fat-fed mice, and alleviated SAT inflammation and liver steatosis . In studies on cattle, the expression of miR-149-5p was substantially higher in bovine preadipocytes than in adipocytes . By targeting CRTC genes, miR-149-5p can prevent the proliferation and differentiation of bovine adipocytes . Mouse mesenchymal stem cells promote differentiation by miR-149, which targets Dab2 . Moreover, miR-149-5p inhibits melanoma cell proliferation by targeting LRIG2 .
MAPK is a kind of serine/threonine protein kinase that widely exists in mammals and can be triggered by external signals, such as bacterial complexes, cytokines, growth factors, and physical stress. The regulation of the MAPK signaling pathway, which is connected to the regulation of a variety of cellular activities, has become one of the most active research areas in cell signal transduction in recent years. The MAPK signaling pathway can phosphorylate a wide range of substrates, including a wide range of nuclear transcription factors and protein kinases, thereby controlling the transcription of related genes; it takes part in several physiological processes, such as cell growth, development, and functional synchronization among cells; it is crucial for adipocyte differentiation, which in turn affects obesity . PPI is composed of proteins that interact with each other, to participate in cell cycle regulation, energy and material metabolism, gene expression regulation, biological signal transmission, and other vital processes. CCL2, MAPK9, MMP3, and HGF have significant protein-protein interactions, which are related to adipogenesis. In vitro knockout of CCL2 reduces macrophage migration to EMSC, thereby reducing fat content . In luteolin-treated epididymal white adipose tissue, MAPK9 gene expression and adipogenesis were found to be positively correlated . MMP3 overexpression prevents the differentiation of 3T3-L1 cells . A high fat-fed fatty liver can be greatly improved by overexpressing HGF . Transcriptomics research is an important means of functional genomics research, which studies gene expression from the level of all mRNA in a cell, and metabolomics studies the expression of cell metabolites. The combination of different omics technologies is an important strategy in the era of functional genomics, which has the characteristics of integrity and high throughput. Therefore, transcriptomics and metabolomics technologies are applied to complex metabolic studies such as obesity, and the results at different levels are integrated to obtain the exact mechanism of the occurrence of research objects. Studies on transcriptomics and metabolomics revealed that CRTC3 overexpression promoted adipogenic differentiation by up-regulating Ca2+, and cAMP signaling pathway in IM and SC adipocytes .
ATP7A, also known as Menkes disease protein, is a P-type ATPase that transports Cu2+ across cell membranes. It is well known that lipid metabolism and copper homeostasis are closely related . This experiment discovered that the miR-149-5p target gene is ATP7A. In the meantime, the ATP7A in the longissimus dorsi, leg muscles, and subcutaneous fat of Queshan Black pigs was significantly higher than that of Yorkshire, and the Cu2+ transporters CTR1 and CTR2 were also increased, which further indicated that miR-149-5p may affect the transport of Cu2+ by targeting ATP7A, thereby inhibiting the formation of IMF. The study found that cold-stimulated brown adipose tissue showed an increase in the level of Cu2+ and a significant increase in ATP7A level . Tao et al. found that the level of Cu2+ concentration of ATP7A knockout mice significantly increased, and the fat content significantly decreased . The differential expression of ATP7A in porcine IM preadipocytes possibly influences the effect of the level of Cu2+ on lipid metabolism in adipocytes.
Based on transcriptomics and metabolomics, we studied the differences in miR-149-5p expression in adipose tissue of local and foreign pigs and revealed the effect of miR-149-5p overexpression on adipocyte metabolism. MiR-149-5p overexpression may inhibit fat accumulation in adipocytes by regulating MAPK, IL-17, PI3K-Akt, and ErbB-related signaling pathways. We discovered that Cu2+ transporter ATP7A is the target gene of miR-149-5p. MiR-149-5p may inhibit the transport of Cu2+ and further inhibit IMF formation by targeting ATP7A. This study provides the first detailed analysis and fresh perspectives of how miR-149-5p overexpression affects the metabolism and transcription of porcine IM preadipocytes and enhances our comprehension of the molecular characteristics of fat and energy metabolism regulated by miR-149-5p. And the results of this study provide a theoretical formation for improving pork quality traits.
Materials and methods
Experimental animals and sample collection
Three Queshan Black pigs and three DLY at a weight of about 110 kg were randomly selected for vitro experiments to study the expression pattern of miR-149-5p and its relationship with fat deposition in pigs. These pigs were fed at the same nutritional level and under the same management conditions. Before being slaughtered, the pigs were required to fast for 24 hours while being given free access to water. After slaughter, the lean and fat weights of the pigs were measured, and the lean and fat meat rates were calculated. We used a Vernier caliper to measure the thickness of the three-point backfat. Marbling was scored by using the US NPPC flesh color comparison board. And the IMF content of the longissimus dorsi was measured by using an extraction method.
Cell culture and transfection
Porcine IM preadipocytes were isolated from 3-day-old Queshan Black piglets by using the published method . The improved Eagle medium contains 89% Durbeko (Gibco, Carlsbad, USA), 10% FBS (Gibco, Carlsbad, USA), and 1% penicillin-streptomycin (Solarbio, Beijing, China) growth medium to cultivate porcine IM preadipocytes. Cells were inoculated into 6-well plates and transfected with Lipofectamine 3000 (Thermo Fisher Scientific, USA) when the cell density reached 70–80%.
RNA extraction and reverse transcription
RNA extraction operation steps refer to the instructions of Trizol reagent (Thermo Fisher Scientific, USA), and RNA concentration was measured by the nucleic acid detector. RNA was reverse transcribed into cDNA by a reverse transcription kit (Accurate Biology, Hunan, China).
Adipogenic marker genes such as CEBPA, PPARG, FABP4, and LIPE were selected to identify the adipogenic differentiation ability of the longissimus dorsi in Queshan Black and Yorkshire. ATP7A was selected to verify the reliability of RNA-seq. The target genes were quantitatively detected by CFX96 Real-time PCR instrument. The primers were synthesized by Shanghai Sangon Biotech Company. Based on the internal reference gene GAPDH, the expression of the target genes in each template was calculated by 2-ΔΔCt method. At the same time, the expression of miR-149-5p takes U6 (RiboBio, Guangzhou, China) as an internal reference.
Dual-luciferase reporter assay
To figure out if miR-149-5p targets ATP7A, two versions of the luciferase reporter plasmid vector (psi-CHECK2) were made: one with a wild-type 3′UTR (WT-ATP7A) and the other with a mutant 3′UTR (Mut-ATP7A). Following that, miR-149-5p mimics and psiCHECK2-ATP7A-3′UTR were cotransfected into 293T cells, the control group was used instead. After lysing the cells, the firefly luciferase activities were found after 48 hours’ transfection by using a luciferase reporter assay kit (Promega, Madison, WI, U.S.A.).
After 8 days of differentiation, miR-149-5p and NC cells were collected, treated with Trizol reagent, stored on dry ice, and sent to Shanghai Personal Bio Company (Shanghai, China) for RNA-seq. The first cDNA was synthesized by using random oligonucleotides and Super Script II. Subsequently, the second strand cDNA was synthesized by using DNA polymerase I and RNase h, and the remaining dangling was converted to a blunt end by exonuclease/polymerase activity and the enzyme was removed. After adenylation at the 3’ end of the DNA fragment, the Illumina PE aptamer oligonucleotide was ligated to prepare for hybridization. In order to select the optimal cDNA fragment with a length of 400–500 bp, the library fragment was purified by using the AMPure XP system (Beckman Coulter, Beverly, CA, USA). Agilent 2100 Bioanalyzer (Agilent, 2100) and Agilent High Sensitivity DNA Kit (Agilent, 5067 − 4626) were used to detect the quality of the library. After homogenization, multiplexed DNA libraries were mixed in equal volumes. The mixed library was gradually diluted and quantified, and PE150 mode sequencing was performed on the Illumina sequencer. After the sample was sequenced on the computer, the image file was obtained, and the original data (Raw Data) of FASTQ was generated by the software of the sequencing platform. Download the reference genome and gene model annotation files directly from the genome website. An index of the reference genome is constructed by using HISAT2v2.0.5 and the paired-end clean reads are compared to the reference genome by using HISAT2v2.0.5. HTSeq (0.9.1) was used to compare the Read Count values on each gene as the original expression of the gene. Differential expression analysis of the two comparison combinations was performed by using DESeq2 software (Version 1.12.4). A significant change was referred to as the threshold (|log2 fold change| > log2 1.2; P < 0.05). GO enrichment and KEGG analysis of DEGs were used by the Gene Denovo Cloud Platform (https://www.omicshare.com/) [65,66,67].
After 8 days of differentiation, miR-149-5p and NC cells were collected, stored on dry ice, and sent to Shanghai Personal Bio Company (Shanghai, China) for LC-MS/MS analysis. Take each sample, make it slowly thawed at 4 °C, and add 1mL methanol: acetonitrile: water (2:2:1, v/v) for full vortex mixing, ultrasonic crushing at low temperature, and incubation for 1 h at -20 °C to precipitate protein, centrifuged at 13,000 rpm and 4 °C for 15 min, with the supernatant freeze-dried and stored at -80 °C for later use. During mass spectrometry analysis, 100 µL acetonitrile aqueous solution (acetonitrile: water = 1:1, v/v) was added for reconstitution, vortexed, and centrifuged at 14000 g and 4 °C for 15 min. Use its company website for drawing. The samples were separated by Agilent 1290 Infinity LC ultra-high performance liquid chromatography (UHPLC) HILIC column. To avoid the influence of instrument detection signal fluctuation, the continuous analysis of samples is carried out in random order. QC samples are inserted into the sample queue to monitor and evaluate the stability of the system and the reliability of the experimental data. The samples were separated by UHPLC and analyzed by Triple TOF 6600 mass spectrometer (AB SCIEX).
MiR-149-5p was compared with NC to reveal the changes of metabolites in porcine IM preadipocytes after overexpression of miR-149-5p. Orthogonal partial least squares discriminant analysis (OPLS-DA) was used for statistical analysis to figure out the changes in metabolites between groups. A significant change was referred as the threshold (|log2 fold change| > log2 1; P < 0.05; VIP > 1). GO enrichment and KEGG analysis of differential metabolites was used by the Gene Denovo Cloud Platform (https://www.omicshare.com/).
The binding site of mRNA and miRNA was predicted by using the online tool RNA22 v2 (https://cm.jefferson.edu/rna22/Interactive/).
The Student’s t-test of SPSS 26.0 software was used, and *P < 0.05 and **P < 0.01 was used as the criterion for a significant difference. The statistical data were expressed as ‘mean ± SD (n = 3)’.
Availability of data and materials
All raw data of transcriptomics sequencing have been deposited to the National Genomics Data Center (NGDC, https://bigd.big.ac.cn) with the dataset accession number CRA009692 (https://ngdc.cncb.ac.cn/gsa/s/eg5Yhm4u). The metabolite data set supporting the results of this article is included within the article, and can be found in the Supplemental Information (Additional file 2: Table S2).
Chang KC, da Costa N, Blackley R, Southwood O, Evans G, Plastow G, Wood JD, Richardson RI. Relationships of myosin heavy chain fibre types to meat quality traits in traditional and modern pigs. Meat Sci. 2003;64(1):93–103.
Xiong Z, Sun DW, Pu H, Gao W, Dai Q. Applications of emerging imaging techniques for meat quality and safety detection and evaluation: a review. Crit Rev Food Sci Nutr. 2017;57(4):755–68.
Kim JM, Lim KS, Ko KB, Ryu YC. Estimation of pork quality in live pigs using biopsied muscle fibre number composition. Meat Sci. 2018;137:130–3.
Bertol TM, de Campos RM, Ludke JV, Terra NN, de Figueiredo EA, Coldebella A, dos Santos Filho JI, Kawski VL, Lehr NM. Effects of genotype and dietary oil supplementation on performance, carcass traits, pork quality and fatty acid composition of backfat and intramuscular fat. Meat Sci. 2013;93(3):507–16.
Jung JH, Shim KS, Na CS, Choe HS. Studies on intramuscular Fat percentage in live swine using real-time ultrasound to determine pork quality. Asian-Australas J Anim Sci. 2015;28(3):318–22.
Poleti MD, Regitano LCA, Souza G, Cesar ASM, Simas RC, Silva-Vignato B, Oliveira GB, Andrade SCS, Cameron LC, Coutinho LL. Longissimus dorsi muscle label-free quantitative proteomic reveals biological mechanisms associated with intramuscular fat deposition. J Proteom. 2018;179:30–41.
Hu H, Wang J, Zhu R, Guo J, Wu Y. Effect of myosin heavy chain composition of muscles on meat quality in Laiwu pigs and Duroc. Sci China C Life Sci. 2008;51(2):127–32.
Lu P, Li D, Yin J, Zhang L, Wang Z. Flavour differences of cooked longissimus muscle from chinese indigenous pig breeds and hybrid pig breed (Duroc×Landrace×Large White). Food Chem. 2008;107(4):1529–37.
Gao SZ, Zhao SM. Physiology, affecting factors and strategies for control of pig meat intramuscular fat. Recent Pat Food Nutr Agric. 2009;1(1):59–74.
Mottram DS. Flavour formation in meat and meat products: a review. Food Chem. 1998;62(4):415–24.
van Laack RL, Stevens SG, Stalder KJ. The influence of ultimate pH and intramuscular fat content on pork tenderness and tenderization. J Anim Sci. 2001;79(2):392–7.
Rosenvold K, Andersen HJ. Factors of significance for pork quality-a review. Meat Sci. 2003;64(3):219–37.
Scherer PE. Adipose tissue: from lipid storage compartment to endocrine organ. Diabetes. 2006;55(6):1537–45.
Longo M, Zatterale F, Naderi J, Parrillo L, Formisano P, Raciti GA, Beguinot F, Miele C. Adipose tissue dysfunction as determinant of Obesity-Associated Metabolic Complications. Int J Mol Sci. 2019;20(9):2358.
Fuster JJ, Ouchi N, Gokce N, Walsh K. Obesity-Induced Changes in Adipose tissue microenvironment and their impact on Cardiovascular Disease. Circ Res. 2016;118(11):1786–807.
Carthew RW, Sontheimer EJ. Origins and Mechanisms of miRNAs and siRNAs. Cell. 2009;136(4):642–55.
Gui L, Raza SHA, Ma B, Easa AA, Althobaiti F, Shukry M, Alotaibi MA, Al Hazani TMI, Dawood MAO, Khan R, et al. CEBPβ binding directly to the promoter region drives CEBPɑ transcription and improves FABP4 transcriptional activity in adipose tissue of yak (Bos grunniens). Res Vet Sci. 2021;141:174–9.
Lefterova MI, Haakonsson AK, Lazar MA, Mandrup S. PPARγ and the global map of adipogenesis and beyond. Trends Endocrinol Metab. 2014;25(6):293–302.
Joe AW, Yi L, Natarajan A, Le Grand F, So L, Wang J, Rudnicki MA, Rossi FM. Muscle injury activates resident fibro/adipogenic progenitors that facilitate myogenesis. Nat Cell Biol. 2010;12(2):153–63.
Collao N, Farup J, De Lisio M. Role of metabolic stress and Exercise in regulating Fibro/Adipogenic progenitors. Front Cell Dev Biol. 2020;8:9.
Klöting N, Berthold S, Kovacs P, Schön MR, Fasshauer M, Ruschke K, Stumvoll M, Blüher M. MicroRNA expression in human omental and subcutaneous adipose tissue. PLoS ONE. 2009;4(3):e4699.
McDaneld TG, Smith TP, Doumit ME, Miles JR, Coutinho LL, Sonstegard TS, Matukumalli LK, Nonneman DJ, Wiedmann RT. MicroRNA transcriptome profiles during swine skeletal muscle development. BMC Genomics. 2009;10:77.
Liu HC, Hicks JA, Trakooljul N, Zhao SH. Current knowledge of microRNA characterization in agricultural animals. Anim Genet. 2010;41(3):225–31.
Ortega FJ, Moreno-Navarrete JM, Pardo G, Sabater M, Hummel M, Ferrer A, Rodriguez-Hermosa JI, Ruiz B, Ricart W, Peral B, et al. MiRNA expression profile of human subcutaneous adipose and during adipocyte differentiation. PLoS ONE. 2010;5(2):e9022.
Rantalainen M, Herrera BM, Nicholson G, Bowden R, Wills QF, Min JL, Neville MJ, Barrett A, Allen M, Rayner NW, et al. MicroRNA expression in abdominal and gluteal adipose tissue is associated with mRNA expression levels and partly genetically driven. PLoS ONE. 2011;6(11):e27338.
Guo Y, Mo D, Zhang Y, Zhang Y, Cong P, Xiao S, He Z, Liu X, Chen Y. MicroRNAome comparison between intramuscular and subcutaneous vascular stem cell adipogenesis. PLoS ONE. 2012;7(9):e45410.
Guo Y, Chen Y, Zhang Y, Zhang Y, Chen L, Mo D. Up-regulated miR-145 expression inhibits porcine preadipocytes differentiation by targeting IRS1. Int J Biol Sci. 2012;8(10):1408–17.
Li M, Qiu Y, Liu S, Dong P, Ning X, Li Y, Yang G, Sun S. [Over-expressed MiR-103 promotes porcine adipocyte differentiation]. Sheng Wu Gong Cheng Xue Bao. 2012;28(8):927–36.
Pan S, Yang X, Jia Y, Li R, Zhao R. Microvesicle-shuttled miR-130b reduces fat deposition in recipient primary cultured porcine adipocytes by inhibiting PPAR-g expression. J Cell Physiol. 2014;229(5):631–9.
Taniguchi M, Nakajima I, Chikuni K, Kojima M, Awata T, Mikawa S. MicroRNA-33b downregulates the differentiation and development of porcine preadipocytes. Mol Biol Rep. 2014;41(2):1081–90.
Lagos-Quintana M, Rauhut R, Lendeckel W, Tuschl T. Identification of novel genes coding for small expressed RNAs. Science. 2001;294(5543):853–8.
Lytle JR, Yario TA, Steitz JA. Target mRNAs are repressed as efficiently by microRNA-binding sites in the 5’ UTR as in the 3’ UTR. Proc Natl Acad Sci U S A. 2007;104(23):9667–72.
Esau C, Davis S, Murray SF, Yu XX, Pandey SK, Pear M, Watts L, Booten SL, Graham M, McKay R, et al. miR-122 regulation of lipid metabolism revealed by in vivo antisense targeting. Cell Metab. 2006;3(2):87–98.
Nakanishi N, Nakagawa Y, Tokushige N, Aoki N, Matsuzaka T, Ishii K, Yahagi N, Kobayashi K, Yatoh S, Takahashi A, et al. The up-regulation of microRNA-335 is associated with lipid metabolism in liver and white adipose tissue of genetically obese mice. Biochem Biophys Res Commun. 2009;385(4):492–6.
Lin Q, Gao Z, Alarcon RM, Ye J, Yun Z. A role of miR-27 in the regulation of adipogenesis. Febs j. 2009;276(8):2348–58.
Iliopoulos D, Drosatos K, Hiyama Y, Goldberg IJ, Zannis VI. MicroRNA-370 controls the expression of microRNA-122 and Cpt1alpha and affects lipid metabolism. J Lipid Res. 2010;51(6):1513–23.
Soh J, Iqbal J, Queiroz J, Fernandez-Hernando C, Hussain MM. MicroRNA-30c reduces hyperlipidemia and atherosclerosis in mice by decreasing lipid synthesis and lipoprotein secretion. Nat Med. 2013;19(7):892–900.
Mohamed JS, Hajira A, Pardo PS, Boriek AM. MicroRNA-149 inhibits PARP-2 and promotes mitochondrial biogenesis via SIRT-1/PGC-1α network in skeletal muscle. Diabetes. 2014;63(5):1546–59.
Li Y, Yang F, Gao M, Gong R, Jin M, Liu T, Sun Y, Fu Y, Huang Q, Zhang W, et al. Mir-149-3p regulates the switch between adipogenic and osteogenic differentiation of BMSCs by targeting FTO. Mol Ther Nucleic Acids. 2019;17:590–600.
Xiao J, Lv D, Zhao Y, Chen X, Song M, Liu J, Bei Y, Wang F, Yang W, Yang C. miR-149 controls non-alcoholic fatty liver by targeting FGF-21. J Cell Mol Med. 2016;20(8):1603–8.
Chen Z, Liu Y, Yang L, Liu P, Zhang Y, Wang X. MiR-149 attenuates endoplasmic reticulum stress-induced inflammation and apoptosis in nonalcoholic fatty liver disease by negatively targeting ATF6 pathway. Immunol Lett. 2020;222:40–8.
Liu Y, Dou Y, Qi K, Li C, Song C, Li X, Li X, Qiao R, Wang K, Han X. CircSETBP1 Acts as a MiR-149-5p sponge to promote intramuscular Fat deposition by regulating CRTCs. J Agric Food Chem. 2022;70(40):12841–51.
Van Deckel MJ, Casteels M, Warnants N, Van Damme L, Boucqué CV. Omega-3 fatty acids in pig nutrition: implications for the intrinsic and sensory quality of the meat. Meat Sci. 1996;44(1–2):55–63.
Hugo A, Roodt E. Significance of porcine fat quality in meat technology: a review. Food reviews international. 2007;23(2):175–98.
Cameron ND, Enser M, Nute GR, Whittington FM, Penman JC, Fisken AC, Perry AM, Wood JD. Genotype with nutrition interaction on fatty acid composition of intramuscular fat and the relationship with flavour of pig meat. Meat Sci. 2000;55(2):187–95.
Aaslyng MD, Oksama M, Olsen EV, Bejerholm C, Baltzer M, Andersen G, Bredie WL, Byrne DV, Gabrielsen G. The impact of sensory quality of pork on consumer preference. Meat Sci. 2007;76(1):61–73.
Li X, Qiao R, Ye J, Wang M, Zhang C, Lv G, Wang K, Li X, Han X. Integrated miRNA and mRNA transcriptomes of spleen profiles between Yorkshire and Queshan black pigs. Gene. 2019;688:204–14.
Qi K, Liu Y, Li C, Li X, Li X, Wang K, Qiao R, Han X. Construction of circRNA-related ceRNA networks in longissimus dorsi muscle of Queshan Black and large White pigs. Mol Genet Genomics. 2022;297(1):101–12.
Qiao R, Li X, Han X, Wang K, Lv G, Ren G, Li X. Population structure and genetic diversity of four Henan pig populations. Anim Genet. 2019;50(3):262–5.
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 PRDM16. Nat Commun. 2016;7:11533.
Zheng S, Guo S, Sun G, Shi Y, Jiang X. Gain of metabolic benefit with ablation of mir-149-3p from Subcutaneous Adipose tissue in Diet-Induced obese mice. Mol Therapy - Nucleic Acids. 2019;18:194–203.
Yu X, Fang X, Gao M, Mi J, Zhang X, Xia L, Zhao Z, Albrecht E, Maak S, Yang R. Isolation and identification of bovine Preadipocytes and Screening of MicroRNAs Associated with adipogenesis. Anim (Basel). 2020;10(5):818.
Khan R, Raza SHA, Junjvlieke Z, Wang X, Wang H, Cheng G, Mei C, Elsaeid Elnour I, 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(7–8):5796–810.
Lu M, Xu L, Min W, Tao G, Tao C. miR-149 promotes the myocardial differentiation of mouse bone marrow stem cells by targeting Dab2. Mol Med Rep. 2018;17(6):8502–9.
Chen W, Zhang J, Xu H, Jie D, Zhang X. The negative regulation of miR-149-5p in melanoma cell survival and apoptosis by targeting LRIG2. Am J Translational Res. 2017;9(9):4331.
Bost F, Aouadi M, Caron L, Binétruy B. The role of MAPKs in adipocyte differentiation and obesity. Biochimie. 2005;87(1):51–6.
Borkar R, Wang X, Zheng D, Miao Z, Zhang Z, Li E, Wu Y, Xu RH. Human ESC-derived MSCs enhance fat engraftment by promoting adipocyte reaggregation, secreting CCL2 and mobilizing macrophages. Biomaterials. 2021;272:120756.
Kwon EY, Choi MS. Luteolin Targets the Toll-Like Receptor Signaling Pathway in Prevention of Hepatic and Adipocyte Fibrosis and Insulin Resistance in Diet-Induced Obese Mice. Nutrients. 2018;10(10):1415.
Wu Y, Lee MJ, Ido Y, Fried SK. High-fat diet-induced obesity regulates MMP3 to modulate depot- and sex-dependent adipose expansion in C57BL/6J mice. Am J Physiol Endocrinol Metab. 2017;312(1):E58–e71.
Kosone T, Takagi H, Horiguchi N, Ariyama Y, Otsuka T, Sohara N, Kakizaki S, Sato K, Mori M. HGF ameliorates a high-fat diet-induced fatty liver. Am J Physiol Gastrointest Liver Physiol. 2007;293(1):G204–210.
Liu J, Wang L, Chen W, Li J, Shan T. CRTC3 regulates the lipid metabolism and adipogenic differentiation of Porcine Intramuscular and Subcutaneous Adipocytes by activating the Calcium Pathway. J Agric Food Chem. 2021;69(25):7243–55.
Morrell A, Tallino S, Yu L, Burkhead JL. The role of insufficient copper in lipid synthesis and fatty-liver disease. IUBMB Life. 2017;69(4):263–70.
Wang C, Liang X, Tao C, Yao X, Wang Y, Wang Y, Li K. Induction of copper and iron in acute cold-stimulated brown adipose tissues. Biochem Biophys Res Commun. 2017;488(3):496–500.
Tao C, Wang Y, Zhao Y, Pan J, Fan Y, Liang X, Cao C, Zhao J, Petris MJ, Li K, et al. Adipocyte-specific disruption of ATPase copper transporting α in mice accelerates lipoatrophy. Diabetologia. 2019;62(12):2340–53.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Kanehisa M. Toward understanding the origin and evolution of cellular organisms. Protein science: a publication of the Protein Society. 2019;28(11):1947–51.
Kanehisa M, Furumichi M, Sato Y, Kawashima M, Ishiguro-Watanabe M. KEGG for taxonomy-based analysis of pathways and genomes. Nucleic Acids Res. 2023;51(D1):D587–d592.
This research was funded by the 14th Five-Year National Key R&D Program, grant number 2021YFD1301202; Agricultural Breeds Research Project of Henan Province, grant number 2022020101.
Ethics approval and consent to participate
All of the experiments involving animals were carried out in accordance with the guidelines for the care and use of experimental animals established by the Ministry of Science and Technology of the People’s Republic of China (Approval Number DWLL20211193). The animal study was reviewed and approved by the Henan Agricultural University Animal Care and Use Committee. All experiments were performed in accordance with relevant guidelines and followed the ARRIVE guidelines for the reporting of animals’ experiments.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Liu, Y., Wei, Y., Dou, Y. et al. Effect of miR-149-5p on intramuscular fat deposition in pigs based on metabolomics and transcriptomics. BMC Genomics 24, 293 (2023). https://doi.org/10.1186/s12864-023-09382-6
- Porcine intramuscular preadipocytes