Next Article in Journal
A Novel Mechanism of bta-miR-210 in Bovine Early Intramuscular Adipogenesis
Next Article in Special Issue
Application of Genetically Engineered Pigs in Biomedical Research
Previous Article in Journal
Genome-Wide Diversity, Population Structure and Demographic History of Dromedaries in the Central Desert of Iran
Previous Article in Special Issue
Genome-Wide Transcriptome and Metabolome Analyses Provide Novel Insights and Suggest a Sex-Specific Response to Heat Stress in Pigs
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 6
10.3390/genes11060600
genes-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

Identification of Molecular Mechanisms Related to Pig Fatness at the Transcriptome and miRNAome Levels

by
Katarzyna Ropka-Molik
1,*,
Klaudia Pawlina-Tyszko
1,
Kacper Żukowski
2,
Mirosław Tyra
3,
Natalia Derebecka
4,
Joanna Wesoły
4,
Tomasz Szmatoła
1,5 and
Katarzyna Piórkowska
1
1
Department of Animal Molecular Biology, National Research Institute of Animal Production, Krakowska 1, 32-083 Balice, Poland
2
Department of Cattle Breeding, National Research Institute of Animal Production, Krakowska 1, 32-083 Balice, Poland
3
Department of Pig Breeding, National Research Institute of Animal Production, Krakowska 1, 32-083 Balice, Poland
4
Laboratory of High Throughput Technologies, Institute of Molecular Biology and Biotechnology, Faculty of Biology, Uniwersytetu Poznanskiego street 6, 61-614 Poznań, Poland
5
University Centre of Veterinary Medicine, University of Agriculture in Kraków, Al. Mickiewicza 24/28, 30-059 Kraków, Poland
*
Author to whom correspondence should be addressed.
Genes 2020, 11(6), 600; https://doi.org/10.3390/genes11060600
Submission received: 13 April 2020 / Revised: 13 May 2020 / Accepted: 27 May 2020 / Published: 29 May 2020
(This article belongs to the Special Issue Pig Genomics and Genetics)
Browse Figures
Versions Notes

Abstract

:
Fat deposition and growth rate are closely related to pork quality and fattening efficiency. The next-generation sequencing (NGS) approach for transcriptome and miRNAome massive parallel sequencing of adipocyte tissue was applied to search for a molecular network related to fat deposition in pigs. Pigs were represented by three breeds (Large White, Pietrain, and Hampshire) that varied in fat content within each breed. The obtained results allowed for the detection of significant enrichment of Gene Ontology (GO) terms and pathways associated directly and indirectly with fat deposition via regulation of fatty acid metabolism, fat cell differentiation, inflammatory response, and extracellular matrix (ECM) organization and disassembly. Moreover, the results showed that adipocyte tissue content strongly affected the expression of leptin and other genes related to a response to excessive feed intake. The findings indicated that modification of genes and miRNAs involved in ECM rearrangements can be essential during fat tissue growth and development in pigs. The identified molecular network within genes and miRNAs that were deregulated depending on the subcutaneous fat level are proposed as candidate factors determining adipogenesis, fatness, and selected fattening characteristics in pigs.

Graphical Abstract

1. Introduction

The quantity of fat is a critical factor influencing meat quality and fattening efficiency in pig production. The fat level and mass of adipocyte tissue control food intake and body weight via leptin secretion [1,2]. Leptin is considered an essential hormone regulating adipose tissue metabolism and energy expenditure. Leptin acts as an appetite regulation factor (hunger inhibition and satiety stimulation) through interaction with hypothalamic receptors [3] and controls glucose/lipid metabolism and body weight homeostasis through gluconeogenesis regulation [4]. Moreover, pigs can be a suitable animal model for studies of fat accumulation and obesity [5]. The identification of genetic factors determining fatness levels in pigs can be a valuable resource for medical research in humans.
In humans, genetic predisposition is a key contributing factor in obesity [6]. It has been estimated that the genetic basis of phenotypic variations in obesity can range from 40% to 70% [6,7]. In pigs, the high heritability coefficient (h2 = 0.5) indicates the genetic background of fatness traits and the ability to improve them by selection [5]. The identification of genome regions, genes, or single nucleotide polymorphisms (SNPs) related to fatness in pigs has been conducted using different methodological approaches. In 2011, the usage of expression and single nucleotide polymorphism microarrays allowed for the detection of the panel of porcine genes relevant to fatness-associated traits [8]. In 2015, the genome-wide association study (GWAS) method was applied to identify QTLs (quantitative trait loci) located on Sus scrofa chromosome 7 (SSC7) and SSC4 strongly associated with growth and fatness traits based on Chinese pig breeds [9]. Guo et al. [10] using the GWAS method also indicated the involvement of the SSC7 region in the regulation of fatness and growth features. Another approach that can allow for the identification of molecular networks related to fat metabolism and growth traits is transcriptomic research concentrated on examining the hypothalamus–pituitary axis. The hypothalamus is a key brain structure controlling food intake and fat accumulation [11]. The pituitary is a crucial gland involved in fetal adipose metabolism via controlling the fatty acid synthesis and receptor-mediated lipolytic response [12]. Pérez-Montarelo et al. [13] confirmed that the usage of network interaction based on hypothalamic transcriptome analysis could pinpoint candidate genes for fatness in pigs. The transcriptomic profiling of the swine pituitary showed the strong association of this gland with postnatal growth and development [14], and also with fatness-associated traits [15].
A significant number of studies have confirmed an important involvement of adipocyte tissue in food intake [16] and body-weight regulation [17]. Such findings indicated the necessity of the analysis of adipocyte tissue to obtain the full view of metabolic processes related to fatness in pigs. To date, in pigs, transcriptomic research has been applied to identify the regulation processes determining feed efficiency [18] and to examine the body’s response to different nutritional treatments [19]. Processes based on adipocyte tissue gene expression profiling in two Portuguese local pig breeds to determine the fat deposition via controlling of de novo fatty acid synthesis have been proposed, [20]. Furthermore, gene expression profiling in combination with genome resequencing has allowed for the detection of genes potentially related to backfat thickness [21].
The present study aimed to perform comprehensive whole transcriptome and miRNAome analyses [20] of adipocyte tissue to identify the molecular networks strongly related to the most crucial production traits, namely, fatness and feed intake in pigs.

2. Materials and Methods

2.1. Animals, Phenotype Data Collection, and Tissue Sampling

Animals used for gene and miRNA expression analyses were selected from a large pig population, based on fatness phenotypic traits measured after dissection. All pigs, representing three breeds (Pietrain, Pi; Hampshire, Hp; and Large White, LW) were unrelated (at least three generations back), females, and were maintained under the same housing and feeding conditions according to a procedure previously described [22]. The pigs were kept in individual pens and fed to a weight of 105 kg (±2.5 kg) (based on the Pig Test Station (SKURTCh) of the National Research Institute of Animal Production in Chorzelów). Half-carcasses were dissected 24 h after slaughter and several carcass characteristics, including fatness traits, were evaluated. The traits included average backfat thickness (ABT), calculated from five measurement points (cm) (at the thickest point over the shoulder; on the back above the joint between the last thoracic and first lumbar vertebrae; and at three points over the loin: above the rostral edge (loin I), above the middle (loin II), and the caudal edge (loin III) of gluteal muscle cross-section), and the weight of peritoneal fat (kg), using a procedure described by Tyra and Żak [23].
Pigs for next-generation sequencing (NGS) analysis were selected based on the most important fatness parameter, i.e., ABT. In each analyzed breed (Pietrain, 8, Hampshire, 8, and Large White, 8), pigs with different fatness phenotypes were selected (4 pigs with high and 4 with low fatness). The separation of the two groups with extreme fatness within each breed allowed for avoidance of any potential breed effect on obtained genes or miRNA expression profiles (Table 1). Immediately after slaughter, the samples of fat tissue (subcutaneous fat) were collected into the RNAlater solution (Ambion, Thermo Fisher, USA) and stored at −20 °C.
The research did not require the approval of the Animal Experimentation Committee due to the fact that all biological material was collected from animals which were slaughtered and dissected, and, after carcass evaluation, meat was intended for consumption. The pigs were maintained and slaughtered according to Directive 2010/63/EU of the European Parliament and the Council of 22 September 2010.

2.2. Whole Transcriptome Sequencing (RNA-seq)

Extracted total RNA (Direct-zol RNA Mini Prep kit (Zymo Research, Orange, CA, USA)) was subjected to quantity and quality controls using a NanoDrop 2000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) and a TapeStation 2200 instrument (Agilent, Santa Clara, CA, USA), respectively. Samples with RIN value (RNA integrity number) above 7 were chosen for further analyses. A quantity of 300 ng of total RNA was used for cDNA library preparation according to the TruSeq RNA Kit v2 kit protocol (Illumina, San Diego, CA, USA). The quality and quantity of cDNA libraries were checked using Qubit 2.0 and TapeStation 2200. The individual cDNA libraries were ligated with adaptors with different indexes and pooled. After pooling, each cDNA library was loaded into at least four flow cell lanes as four technical replicates. Sequencing of RNA was performed on the HiScanSQ system (Illumina, San Diego, CA, USA) using the TruSeq SR Cluster Kit v3- CBOT-HS kit and TruSeq SBS Kit v 3 - HS chemistry (Illumina, San Diego, CA, USA) according to the standard protocols and with 81 single-end reads.
The quality of raw data was analyzed using FastQC software (overrepresented sequences, sequence duplication levels, adapter contents, sequence length distribution, sequence quality scores, and GC content). The adapter’s sequence reads shorter than 36 bp and/or reads with a quality score lower than 20 were removed using Flexbar software v2.5 [24]. The estimation of transcript abundance was conducted using RSEM software (v1.2.22) [25] supported by STAR aligner software (STAR_2.4.2a) [26]. The reads were aligned to the Sus scrofa reference genome (assembly Sscrofa11.1, Ensemble release 94). The whole procedure was followed by the alignment parameter evaluation using ENCODE3’s STAR-RSEM pipeline. Differentially-expressed genes (DEGs) were identified according to the following criteria: fold change ≥ |1.0| and adjusted p-value < 0.05 separately for each breed using DESeq2 software (version 1.12.4) [26].

2.3. MicroRNA Sequencing

MicroRNA libraries were prepared with the use of the NEBNext Multiplex Small RNA Library Prep Set for Illumina (New England Biolabs, Ipswich, MA, USA) according to the protocol. The libraries were barcoded with different indexed primers to allow multiplexing of the samples during NGS sequencing. The amplified libraries were purified and size-selected on a Novex 6% TBE PAGE gel (Invitrogen, Thermo Fisher Scientific, Waltham, MA, USA). The quality and quantity of obtained libraries were measured with a Qubit 2.0 Fluorometer (Thermo Fisher Scientific, USA) and size-assessed with a 2200 TapeStation instrument (Agilent, Santa Clara, CA, USA). The whole miRNA profile sequencing was performed in 36 cycles on a HiScanSQ (Illumina, San Diego, CA, USA) instrument according to the manufacturer’s protocol.
The raw miRNA reads were quality controlled using FastQC software (parameters as for the transcriptome) [27]. The adaptor trimming and filtering were performed with the UEA sRNA Workbench v3.2 helper tool [28], while the miRCat tool v1.0 was used to detect known and novel miRNA sequences. The identification was carried out on the basis of the Sus scrofa genome (assembly Sscrofa 10.2) and miRBase v21.0 (Griffiths-Jones lab at the Faculty of Biology, Medicine, and Health, University of Manchester, USA) [29,30], due to the lack of 11.1 assembly in miRbase dataset. The obtained miRNAs genome localization based on 10.2 reference were converted to 11.1 reference using NCBI Genome Remapping tool. The default animal parameters were set except for minimum abundance (6 reads), minimum length (17 nt), and maximum length (25 nt). In the next step, potential other non-coding RNAs were identified and excluded employing the RNAcentral database (the RNAcentral Consortium, 2017). Moreover, isomiR-SEA software (version 1.60) [31] with the default settings was used to identify microRNA length and sequence variants (isomiRs). miRNAs differentially expressed between pigs with low and high fatness in each breed were detected using DESeq2 software (v.1.12.4) [26]. The identification of experimentally-validated target genes and enriched biological processes (KEGG, GO) was performed using the mirPath v.3.0 DIANA Tools web application (DIANA-Lab, Athens, Greece) [32], including TarBase v7.0 (DIANA-Lab, Athens, Greece) as a reference [33]. The analysis was performed based on human homologues deposited in miRBase (21.0) due to the lack of data for pig microRNAs.

2.4. Validation of NGS Results

Validation of the RNA-seq results was carried out using a real-time PCR method, for 10 selected DEGs. The primers for selected DEGs were designed based on reference sequences showed in Supplementary Table S1 and using Primer3 Input (version 0.4.0) The detailed information about the chosen genes is presented in Supplementary Table S1. The cDNA was prepared from 300 ng of total RNA using a Maxima First Strand cDNA Synthesis Kit for RT-qPCR (Thermo Fisher Scientific, Waltham, MA, USA) according to the protocol. The exact gene expression levels were measured on a QuantStudio 7 flex instrument (Applied Biosystems, Thermo Fisher Scientific, Waltham, MA, USA) using Sensitive RT HS-PCR EvaGreen Mix (A&A Biotechnology, Gdynia, Poland). The reaction was carried out in three replications for each gene. The expression was calculated using the delta-delta CT method, according to Pfaffl [34] and based on two reference controls OAZ1 [35] and RPS29 [36].
Validation of the detected expression levels of 10 selected microRNAs was carried out using the qPCR method. A quantity of 10 ng of total RNA was reverse transcribed to cDNA using a TaqMan Advanced miRNA cDNA Synthesis Kit (Thermo Fisher Scientific, USA) following the protocol. The quantitative estimation of miRNAs was performed on a QuantStudio 7 flex instrument using TaqMan® Advanced miRNA Assays labelled with VIC or FAM fluorescent dye and with TaqMan® Fast Advanced Master Mix (Thermo Fisher Scientific, USA) according to the standard protocol (Supplementary Table S1). The relative expression levels of selected miRNAs were calculated using the delta-delta CT method according to Pfaffl [34] and based on mir-451a as a reference control. The comparison between NGS data (RNA-seq, miRNA-seq) and relative quantity obtained by the real-time PCR method was performed using the Pearson correlation (SAS software, version 8.02).

3. Results

3.1. RNA Sequencing Results—Differentially-Expressed Genes

The whole transcriptome sequencing of adipocyte tissue allowed for the identification of differentially-expressed genes (DEGs) between pigs varying in backfat thickness. The investigation of three pig breeds enabled the detection of genes related to fatness traits regardless of breed factor. According to the NGS approach, 167 DEGs were identified between Large White pigs with different fatness phenotypes (85 up-regulated and 82 down-regulated genes in pigs with thicker backfat). A total of 32 of these were recognized as novel genes. In turn, in the Pietrain breed, 247 DEGs were detected (53 up-regulated and 194 down-regulated), of which 50 were novel genes or uncharacterized proteins. The highest number of DEGs was identified for the Hampshire breed: 128 up-regulated and 173 down-regulated genes, for a total of 301 genes (46 novel).
The comparison of the detected sets of DEGs across breeds enabled the identification of eight known genes (SCD, SFRP2, MMRN1, PCK1, TNC, C2, CXCl10, and CXCL9) and one novel gene common to all three breeds (Figure 1). The largest group of common genes identified as differentially-expressed was found for Large White and Pietrain breeds (44 DEGs and 12 novel genes). The lowest similarity in the identified gene panel was observed for Hampshire pigs and the other two breeds.
The Gene Ontology (GO) analysis performed for the DEGs identified in at least two breeds showed significantly-overrepresented GO terms (Table 2). The highest number of genes was assigned to innate immune response (GO:0045087), 14 genes; inflammatory response (GO:0006954), 11 genes; and positive regulation of apoptotic process (GO:0043065), 10 genes. The DEGs were also included in GO terms related to lipid metabolism: fatty acid biosynthetic process, long-chain fatty acid biosynthetic process, positive regulation of fat cell differentiation, positive regulation of mast cell degranulation, and response to excessive food intake (Table 2; David software). Furthermore, GO and pathways analyses using over-representation analysis (WebGestalt software) confirmed that detected DEGs were involved in the biosynthesis of unsaturated fatty acids, fatty acid metabolism, and regulation of fat cell differentiation and brown fat differentiation (not significant GO) (Figure 2). The DEGs associated with several identified GO terms were as follows: ACACA (acetyl-CoA Carboxylase α), SCD (stearoyl-CoA desaturase), SCD5 (stearoyl-CoA desaturase 5), FASN (fatty acid synthase), LEP (leptin), and CTGF (connective tissue growth factor) (Table 3; Figure 3).

3.2. MiRNA Sequencing Results—Differentially-Expressed miRNAs

The comparison of whole miRNAome profiles between groups of pigs with significant-different fatness traits allowed for the identification of differentially-expressed (DE) miRNAs: 46 for LW (21 up-regulated and 25 down-regulated genes in pigs with thicker backfat), 61 for Pietrain (36 up-regulated and 25 down-regulated), and 41 for Hampshire (8 up-regulated and 31 down-regulated). As for DEGs, analogous comparison across breeds was also carried out for DE miRNAs and showed the presence of 14 miRNAs common to all three breeds (Figure 4). The lowest number of common miRNAs was detected for Hampshire and Large White pigs. Detailed data on the chromosomal localization of the identified miRNAs according to Sscrofa10.2 and Sscrofa11.1 assemblies are in the Supplementary Table S2.
The GO analysis performed using the mirPath v.3 DIANA tool showed that differentially-expressed miRNAs were involved in extracellular matrix organization (GO:0030198), extracellular matrix disassembly (GO:0022617), and innate immune response (GO:0045087), where the highest number of predicted targeted genes for miRNAs was detected. Eight miRNAs were associated with cellular lipid metabolic process (GO:0044255) (Table 4). Moreover, miR-100-5p, miR-143-3p, miR-10b-5p, and miR-24-3p were involved in fatty acid metabolism and fatty acid biosynthesis (Figure 5). The predicted targeted genes for differentially-expressed miRNAs are summarized in Supplementary Table S3.

3.3. Analysis of Pathways Common for DEGs and DE miRNAs

The pathways analysis was performed for both DEGs and DE miRNAs identified in at least two breeds. Furthermore, the predicted genes regulated by the identified miRNAs, and involved in selected pathways, detected using mirPath v.3 DIANA tool are shown in Table 5.
An enriched pathway significant for both DEGs and DE miRNAs was fatty acid metabolism (KEGG ID: hsa01212 for miRNA/ ssc01212 for DEGs) for which four miRNAs and five genes were detected. Three of all DEGs involved in fatty acid metabolism overlapped with predicted genes regulated by detection miRNAs, namely: SCD, FASN, and ACACA. A similar situation was observed for the fatty acid biosynthesis pathway (p-value for miRNAs < 1.0 × 10−325; not significant for DEGs); FASN and ACACA differentially-expressed genes overlapped with miRNA-predicted targets. The pathways identified as significant and common for DE miRNAs and DEGs were also the Hippo-signaling pathway, ECM−receptor interaction, cell cycle, and P53-signaling pathways.
The pathway analysis of DE miRNAs also showed significant enrichment for the set of pathways related to toll-like receptors (toll-like receptor 2, 3, 4, 5, 9, and 10 signaling pathways; toll-like receptor TLR1:TLR2 signaling pathways; and toll-like receptor TLR6:TLR2 signaling pathways), as well as TRIF-dependent toll-like receptor signaling pathways (Figure S1).

3.4. Validation of Obtained Data

The validation of RNA-seq results was performed for eight DEGs and six miRNAs (Table 6). For the majority of analyzed DEGs, the RNA-seq data and RQ results showed high and significant correlation coefficients (from 0.66 to 0.95), which confirmed the reliability of the NGS results. The highest similarity of results was obtained for PCK1, ACACA, and LEP genes. For four miRNAs, high correlation coefficients were detected, but these were significant only for miR-26a-5p, mir-100-5p, and mir-103a-3p.

4. Discussion

In pigs, fatness traits are one of the important production features due to their strong relationship with meat quality and fattening performance. The inverse correlation between lean meat content and fatness traits results in the decrease of the fatness level in carcasses as meatiness increases [37]. Moreover, the proper fat content, including of intramuscular fat (IMF), is critical to achieving the desired meat quality parameters. The adipocyte tissue, as a secretory organ, also determines the regulation of food intake; thus, the body fat content can influence food intake and body weight [2]. In turn, the food intake and feed conversion ratio (FCR) are primary factors which determine the economic efficiency of pig production.
The high phenotypic variability of fatness characteristics, as well as the heritability coefficient (about 0.5), strongly indicate the genetic background of this group of traits and the possibility of their improvement or modification using genetic markers [5]. From a breeding perspective, the most important step would be to establish such a genetic marker that enables the controlling of pig fatness without losing the already-achieved level of meatiness. The latest successes of high throughput genetic methods, such as NGS sequencing technology and the GWAS approach, have become more applicable and create new possibilities in the search for and identification of the molecular background of phenotypic variations.
In the presented report, comprehensive whole mRNA and miRNA profiling of adipocyte tissue was applied to detect the molecular network related to fat deposition in pigs. The use of three pig breeds representing different usage types (maternal line, Large White; sire lines, Pietrain and Hampshire) allowed indication of a universal genetic basis of fatness characteristics that was not associated with any breed. The transcriptome sequencing of fat tissue delivered from pigs diverse in terms of backfat thickness enabled the pinpointing of significant, enriched Gene Ontology terms and pathways possibly related to fat deposition. For both DEGs and DE miRNAs, the fatty acid metabolism pathway was detected as significant. This pathway involved three DEGs, namely, ACACA, FASN, and SCD, and differentially-expressed miRNAs that recognized these genes as their targets. Moreover, the stearoyl-CoA desaturase gene (SCD) was identified as significant for all three analyzed breeds, and the SCD expression level was higher in pig groups with thicker fat cover. The stearoyl-CoA desaturase plays a vital role during the biosynthesis of monounsaturated fatty acids, within palmitoyl- and oleoyl-CoA, making them readily available to the body [38]. The higher SCD expression observed in this study, in pigs with higher fat content was in accord with previously-found human data. Hulver et al. [39] showed that an increased SCD level is correlated to obesity. Similarly, decrease of SCD transcription leads to a reduction of adiposity [40]. The deficiency of stearoyl-CoA desaturase results in leanness and leads to an increase of metabolic rate as well as insulin sensitivity [38]. Moreover, a study performed on ob/ob mice indicated that SCD is involved in metabolic response to leptin and down-regulation of SCD can be a key element of the metabolic actions of leptin [41].
Reports of many authors have indicated the association of the SCD gene and selected fatness traits in pigs. Using the GWAS approach, Yang et al. [42] showed a possible relationship between the SCD locus and C18:0 content. Similarly, results from the expression Genome-Wide Association Study (eGWAS) indicated that the SCD gene is related to intramuscular fat composition in pigs [43]. Xing et al. [21] compared transcriptomes of subcutaneous adipose tissue of Landrace pigs depending on variable backfat levels, which showed significantly-higher SCD expression in animals with increased fatness. The authors proposed stearoyl-CoA desaturase as a candidate gene for fat deposition. Our results, based on three different pig breeds, confirmed previous observations, thus strongly supporting an essential role of the SCD gene in the determination of fat-associated traits.
Moreover, the present miRNAome sequencing identified that expression of miR-24-3p is also dependent on subcutaneous fat level. In turn, the mirPath v.3 and TarBase v7.0 DIANA bioinformatic tools found that this microRNA modulates expression of the SCD gene, as well as fatty acid synthase (FASN), acetyl-CoA carboxylase α (ACACA), acetyl-coA acyltransferase 1 and 2 (ACAA1 and ACAA2), and malonyl coA-acyl carrier protein transacylase (MCAT). miR-24-3p has been established to play a role in adipogenesis via activating adipocyte differentiation by targeting genes engaged in MAPK7-signaling pathways [44,45]. To date, only a few reports have been published concerning the possible involvement of miR-24-3p in lipid metabolism processes and obesity development. The present study indicates that this miRNA can be related to the regulation of fat deposition processes by targeting the key genes of fatty acid metabolism.
The analogues gene expression profile—increased transcript level in all fatty pigs independent of breed—was also observed for ACACA and FASN genes. A recent study applied the RNA-seq method to analyze the whole transcriptome profile of subcutaneous fat in native Puławska pigs that differed in backfat tissue [46]. The authors identified the higher transcript level of all three genes—ACACA, FASN, and SCD—in pigs with increased backfat thickness. Previous research showed the significant association of the FASN gene with meat quality and the ACACA gene with IMF content. In Iberian pigs, the ACACA gene was proposed as a candidate gene responsible for intramuscular content of saturated fatty acids and monounsaturated fatty acids fatty acids [47]. Furthermore, the up-regulation of FASN and SCD genes were reported in Alentejano pigs and related to the higher fat deposition observed in this breed [20]. The research performed on Iberian × Landrace pig cross lines indicated a significant association of the ACACA gene and IMF palmitic fatty acid percentage [48]. Similarly, the FASN gene was previously correlated to meat quality, fatty acid content, and composition [49,50]. Our research, supported by previous studies, strongly confirms that the identified DEGs determine fat deposition in pigs.

Enriched Metabolic Process and Pathways Associated with Fat Deposition

Genes whose expression varied between thick and thin backfat groups belonged to several processes and pathways directly related to fatty lipid metabolism, within the fatty acid biosynthetic process, the long-chain fatty acid biosynthetic process, positive regulation of fat cell differentiation, and response to dietary excess. The molecular network of food intake controlled by adipocyte tissues has been widely investigated and described. One of the main factors determining the energy balance by hunger inhibition is leptin [51]. The DEGs set comparison showed that LEP gene expression was increased in fat tissue of more obese pigs independent of breed. These findings are in accordance with observations made in humans, where the overexpression of the leptin gene in subcutaneous and omental adipose tissues in obese patients was detected [52]. Adipocyte cells accumulate triglyceride and increase the synthesis of leptin during their growth. Leptin acts via the hypothalamus to control energy balance and feeding behavior [53] Thus, leptin, also called an anti-obesity hormone, reduces food intake and increases energy expenditure [54]. Likewise, the other two DEGs identified in the present study—PPARGC1A (peroxisome proliferator-activated receptor-γ co-activator-1alpha) and PCSK1N (proprotein convertase subtilisin/kexin type 1 inhibitor)—have been associated with obesity in humans. A polymorphism in the PPARGC1A gene is related to obesity and type 2 diabetes [55], while the down-regulation of PCSK1N gene expression is also associated with obesity [56]. In the present research, PPARGC1A expression was up-regulated in thicker backfat tissue in all three investigated pig breeds, whereas expression of the PCSK1N gene was lowered in obese pigs represented by Hampshire and Large White breeds. This suggests that all identified DEGs related to mast cell (MC) degranulation can affect the increase of fat tissue mass and lipid metabolism.
Research performed in humans and animals strongly indicates that MCs are involved in activation of adipocytes and recruitment of inflammatory cells [57]. The positive regulation of the MC degranulation-GO term was identified as significantly deregulated in fat tissue dependent on fat cover thickness. The research allowed for the detection of FGR, FCER1A, FCER1G, and ZAP70 as DEGs associated with MC degranulation. The function of MC degranulation is associated with the release of multiple enzymes which can influence adipocyte tissue modification. In vitro studies on human white adipocyte tissue showed that MCs can promote adipose initiation in response to cold [58]. Furthermore, MC cells induced remodeling of adipose tissue extracellular matrix [57].
Interestingly, the obtained data showed the significant enrichment of ECM organization (GO:0030198) and ECM−receptor interaction pathways. The panel of DEGs was identified including genes coding for collagens, thrombospondins, and laminin, which belong to extracellular matrix remodeling. A recent report highlighted that ECM plays a critical role in adipose tissue expansion through controlling cell migration during body growth and development [59]. It has been proven that adipose tissue expansion is strongly related to ECM remodeling at the level of matrix individual components—collagens, thrombospondins, metalloproteinases, and their inhibitors (Tissue Inhibitor of Metalloproteinase - TIMPs) [60]. The present research allowed the identification of DE miR-145-5p and its target gene, tenascin C (TNC), both of which are involved in the ECM−receptor interaction pathway. Tenascin C, as a glycoprotein member of ECM, is related to tissue developmental processes, while miR-145 is involved in abdominal obesity, insulin resistance, and lipid metabolism [61]. The exact role of miR-145 has not been well established but the increased expression of miR-145 was observed in adipose and liver tissues in obese patients [62]. It has been proposed that this miRNA plays a role in lipolysis and our results supports this thesis, and also confirm that miR-145 overexpression in fat tissue is associated with obese individuals, which was independently observed in the investigated pig breeds.
Another interesting observation is the down-regulation of the TNC gene identified in pigs with thicker fat cover. Tenascin C determines the formation of the collagen network in adipocyte tissue controlling cell migration and proliferation [63]. On the other hand, the excessive cross-linking of adipocytes by fibrotic components can lead to the reduction of their metabolic activity [64]. We hypothesized that the increased expression of the TNC gene in thinner backfat can be related to significant modification of extracellular matrix components, which may result in reduced adipocyte activity and a decreased rate of fat tissue development. Tenascin C is also related to activation of the toll-like receptor 4 (TLR4) signaling pathway, which triggers the obesity-induced inflammatory response [65]. Our results also showed the significant enrichment of both toll-like receptor TLR1:TLR2 signaling and toll-like receptor TLR6:TLR2 signaling pathways for DE miRNAs. Gene expression profiling of human peripheral blood mononuclear cells (PBMCs) showed the increased transcript level of TLR2 and TLR4 genes in obese patients [66]. The increase in expression of the TLR4 gene, reported in differentiating adipocytes in db/db mice suggests that this gene is critical during obesity development processes [67]. In the present study, in addition to the detection of GO terms related to TLR2 and TLR6, the set of DE miRNAs targeting the TLR4 gene was identified, namely, the members of the let-7 family (let-7a-5p, let-7i-5p, and let-7d-5p). To date, many reports have indicated a significant role of let-7 miRNAs in regulation of adipogenesis, lipolysis, and insulin resistance [68]. The differential expression of several let-7 family members may confirm that, beside TLR genes, this miRNA family plays a key role in the fat deposition rate and lipid metabolism.
The other proven role of ECM receptor pathways during fat tissue development is the regulation of adipocyte apoptosis and/or necrosis processes [60]. The results obtained in our study also indicated expression changes of genes involved in the positive regulation of apoptotic process GO term and cell cycle pathways. Other members of the EC matrix—thrombospondins—such as those identified in our study (THBS2, THBS3, and THBS4), are mainly detected in visceral adipose tissues and their increased expression level is associated with obesity. Moreover, thrombospondin 1 has been proposed as a novel marker related to obesity and metabolic syndrome [69]. The identification of the panel of genes involved in ECM receptor pathways and extracellular matrix assembly confirmed that such signaling is essential for fat deposition processes in pigs irrespective of a breed factor.

5. Conclusions

The backfat thickness and growth rate are closely related to pork quality and fattening efficiency, as well as reproductive and immune characteristics. Taking into account whole transcriptome profiling of pigs varying in fat content across three pig breeds, the significant enrichment of Gene Ontology terms and pathways associated directly and indirectly with fat deposition was detected (ACACA, ACOX3, FASN, SCD5, SCD, and PLPP1). The differentially-expressed miRNAs and genes were involved in fatty acid metabolism, positive regulation of fat cell differentiation, and the inflammatory response. Moreover, the results showed that adipocyte tissue content regulated the expression of leptin and other genes related to a response to dietary excess. Our results confirmed previous findings in humans that showed ECM organization and disassembly are fundamental for fat tissue growth and development. The modification of genes and miRNAs involved in ECM rearrangements can also be essential during fat tissue growth and development in pigs. We also suggest that mast cell degranulation can be one of the significant processes associated with adipocyte tissue development. The pinpointed molecular networks within genes and miRNAs deregulated by subcutaneous fat level are proposed as candidate factors determining adipogenesis, fatness, and selected fattening characteristics in pigs. Identified DEGs and DE miRNAs should be investigated in the future in terms of their use as genetic markers associated with pig production traits.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/6/600/s1, Table S1. Detailed information about DEGs and miRNAs whose expression was estimated using real-time PCR methods. Table S2. Detailed data on the chromosomal localization of the identified miRNAs according to Sscrofa10.2 and Sscrofa11.1 (“LW” stands for Large White samples, “Pi” denotes Pietrain samples, and “Hp” stands for Hampshire samples; “NA” denotes potentially novel miRNA identified in this study). Table S3. TarBase experimentally supported interactions—the predicted targeted genes (N) for each of the 14 detected miRNAs. Figure S1. The significant enrichment molecular pathways in which 14 microRNAs commonly detected in all three pig breeds were detected (DIANA-miRPath v3.0, DIANA-Lab, Athens, Greece). Figure S2. The significant enrichment Hippo-signaling pathway (KEGG database, hsa04390) involved genes targeted for identified differentially-expressed miRNAs. Targeted genes are marked orange (genes included in more than one pathway) and yellow (genes included only in one pathway). Figure S3. The significant enrichment cell cycle (KEGG database; hsa04110) involved genes targeted for identified differentially-expressed miRNAs. Targeted genes are marked orange (genes included in more than one pathway) and yellow (genes included only in one pathway).

Author Contributions

For research articles with several authors, a short paragraph specifying their individual contributions must be provided. The following statements should be used “Conceptualization, K.R.-M., K.P.-T., and K.P.; methodology, K.R.-M. and K.P.-T.; software, K.Ż. and T.S.; validation, K.P. and K.R.-M.; formal analysis, K.R.-M., K.P.-T., N.D., and J.W.; investigation, K.R.-M., K.P.-T., N.D., and J.W.; resources, M.T.; data curation, K.Ż. and T.S.; writing—original draft preparation, K.R.-M.; writing—review and editing, K.P.-T., K.P., K.Ż., T.S., M.T., N.D., and J.W.; visualization, K.R.-M.; supervision, K.R.-M.; project administration, K.R.-M. and K.P.; funding acquisition, K.R.-M. and K.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by statutory activity of the National Research Institute of Animal Production, grant number 04.-0.14.1 (cost of biological material collection and storage) and grant number 01-18-05-21 (cost of genetic analyses).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Friedman, J.M.; Halaas, J.L. Leptin and the regulation of body weight in mammals. Nature 1998, 395, 763–770. [Google Scholar] [CrossRef]
  2. Hussain, Z.; Khan, J.A. Food intake regulation by leptin: Mechanisms mediating gluconeogenesis and energy expenditure. Asian Pac. J. Trop. Med. 2017, 10, 940–944. [Google Scholar] [CrossRef] [PubMed]
  3. Elmquist, J.K.; Elias, C.F.; Saper, C.B. From lesions to leptin: Hypothalamic control of food intake and body weight. Neuron 1999, 22, 221–232. [Google Scholar] [CrossRef] [Green Version]
  4. Coelho, M.; Oliveira, T.; Fernandes, R. Biochemistry of adipose tissue: An endocrine organ. Arch. Med. Sci. 2013, 9, 191–200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Switonski, M.; Stachowiak, M.; Cieslak, J.; Bartz, M.; Grzes, M. Genetics of fat tissue accumulation in pigs: A comparative approach. J. Appl. Genet. 2010, 51, 153–168. [Google Scholar] [CrossRef]
  6. Allison, D.B.; Neale, M.C.; Kezis, M.I.; Alfonso, V.C.; Heshka, S.; Heymsfield, S.B. Assortative mating for relative weight: Genetic implications. Behav. Genet. 1996, 26, 103–111. [Google Scholar] [CrossRef]
  7. Ashrafi, K. Obesity and the Regulation of Fat Metabolism. WormBook. 2007. Available online: https://www.ncbi.nlm.nih.gov/books/NBK19757/ (accessed on 28 May 2020).
  8. Ponsuksili, S.; Murani, E.; Brand, B.; Schwerin, M.; Wimmers, K. Integrating expression profiling and whole-genome association for dissection of fat traits in a porcine model. J. Lipid Res. 2011, 52, 668–678. [Google Scholar] [CrossRef] [Green Version]
  9. Qiao, R.; Gao, J.; Zhang, Z.; Li, L.; Xie, X.; Fan, Y.; Cui, L.; Ma, J.; Ai, H.; Ren, J.; et al. Genome-wide association analyses reveal significant loci and strong candidate genes for growth and fatness traits in two pig populations. Genet. Sel. Evol. 2015, 47, 17. [Google Scholar] [CrossRef] [Green Version]
  10. Guo, Y.; Huang, Y.; Hou, L.; Ma, J.; Chen, C.; Ai, H.; Huang, L.; Ren, J. Genome-wide detection of genetic markers associated with growth and fatness in four pig populations using four approaches. Genet. Sel. Evol. 2017, 49, 21. [Google Scholar] [CrossRef] [Green Version]
  11. Timper, K.; Brüning, J.C. Hypothalamic circuits regulating appetite and energy homeostasis: Pathways to obesity. DMM Dis. Model. Mech. 2017, 10, 679–689. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Hausman, D.B.; Hausman, G.J.; Martin, R.J. Influence of the pituitary on lipolysis and lipogenesis in fetal pig adipose tissue. Horm. Metab. Res. 1993, 25, 17–20. [Google Scholar] [CrossRef] [PubMed]
  13. Pérez-Montarelo, D.; Madsen, O.; Alves, E.; Rodríguez, M.C.; Folch, J.M.; Noguera, J.L.; Groenen, M.A.M.; Fernández, A.I. Identification of genes regulating growth and fatness traits in pig through hypothalamic transcriptome analysis. Physiol. Genom. 2014, 46, 195–206. [Google Scholar] [CrossRef] [Green Version]
  14. Shan, L.; Wu, Q.; Li, Y.; Shang, H.; Guo, K.; Wu, J.; Wei, H.; Zhao, J.; Yu, J.; Li, M.H. Transcriptome profiling identifies differentially expressed genes in postnatal developing pituitary gland of miniature pig. DNA Res. 2014, 21, 207–216. [Google Scholar] [CrossRef] [Green Version]
  15. Piórkowska, K.; Żukowski, K.; Tyra, M.; Szyndler-Nędza, M.; Szulc, K.; Skrzypczak, E.; Ropka-Molik, K. The pituitary transcriptional response related to feed conversion in pigs. Genes 2019, 10, 712. [Google Scholar] [CrossRef] [Green Version]
  16. Henry, B.A.; Clarke, I.J. Adipose Tissue Hormones and the Regulation of Food Intake. J. Neuroendocrinol. 2008, 20, 842–849. [Google Scholar] [CrossRef] [PubMed]
  17. Havel, P.J. Role of adipose tissue in body-weight regulation: Mechanisms regulating leptin production and energy balance. Proc. Nutr. Soc. 2000, 59, 359–371. [Google Scholar] [CrossRef] [Green Version]
  18. Xu, Y.; Qi, X.; Hu, M.; Lin, R.; Hou, Y.; Wang, Z.; Zhou, H.; Zhao, Y.; Luan, Y.; Zhao, S.; et al. Transcriptome analysis of adipose tissue indicates that the cAMP signaling pathway affects the feed efficiency of pigs. Genes 2018, 9, 336. [Google Scholar] [CrossRef] [Green Version]
  19. Benítez, R.; Trakooljul, N.; Núñez, Y.; Isabel, B.; Murani, E.; De Mercado, E.; Gómez-Izquierdo, E.; García-Casco, J.; López-Bote, C.; Wimmers, K.; et al. Breed, diet, and interaction effects on adipose tissue transcriptome in iberian and duroc pigs fed different energy sources. Genes 2019, 10, 589. [Google Scholar] [CrossRef] [Green Version]
  20. Albuquerque, A.; Óvilo, C.; Núñez, Y.; Benítez, R.; López-Garcia, A.; García, F.; Félix, M.D.R.; Laranjo, M.; Charneca, R.; Martins, J.M. Comparative transcriptomic analysis of subcutaneous adipose tissue from local pig breeds. Genes 2020, 11, 422. [Google Scholar] [CrossRef] [Green Version]
  21. Xing, K.; Zhu, F.; Zhai, L.; Chen, S.; Tan, Z.; Sun, Y.; Hou, Z.; Wang, C. Identification of genes for controlling swine adipose deposition by integrating transcriptome, whole-genome resequencing, and quantitative trait loci data. Sci. Rep. 2016, 6, 1–10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Ropka-Molik, K.; Dusik, A.; Piórkowska, K.; Tyra, M.; Oczkowicz, M.; Szmatoła, T. Polymorphisms of the membrane-associated ring finger 4, ubiquitin protein ligase gene (MARCH4) and its relationship with porcine production traits. Livest. Sci. 2015, 178, 18–26. [Google Scholar] [CrossRef]
  23. Tyra, M.; Żak, G. Analysis of the possibility of improving the indicators of pork quality through selection with particular consideration of intramuscular fat (imf) conntenttent*. Ann. Anim. Sci. 2013, 13, 33–44. [Google Scholar] [CrossRef] [Green Version]
  24. Dodt, M.; Roehr, J.T.; Ahmed, R.; Dieterich, C. FLEXBAR-flexible barcode and adapter processing for next-generation sequencing platforms. Biology 2012, 1, 895–905. [Google Scholar] [CrossRef] [Green Version]
  25. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics 2011, 12, 323. [Google Scholar] [CrossRef] [Green Version]
  26. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15. [Google Scholar] [CrossRef] [Green Version]
  27. Andrews, S. FastQC A Quality Control tool for High Throughput Sequence Data 2010. Available online: https://www.bioinformatics.babraham.ac.uk/projects/fastqc/ (accessed on 28 May 2020).
  28. Stocks, M.B.; Moxon, S.; Mapleson, D.; Woolfenden, H.C.; Mohorianu, I.; Folkes, L.; Schwach, F.; Dalmay, T.; Moulton, V. The UEA sRNA workbench: A suite of tools for analysing and visualizing next generation sequencing microRNA and small RNA datasets. Bioinformatics 2012, 28, 2059–2061. [Google Scholar] [CrossRef]
  29. Griffiths-Jones, S.; Saini, H.; van Dongen, S.; Enright, A. miRBase: Tools for microRNA genomics. Nucleic Acids Res. 2008, 36, D154–D158. [Google Scholar] [CrossRef] [Green Version]
  30. Griffiths-Jones, S. miRBase: MicroRNA sequences, targets and gene nomenclature. Nucleic Acids Res. 2006, 34, D140–D144. [Google Scholar] [CrossRef]
  31. Urgese, G.; Paciello, G.; Acquaviva, A.; Ficarra, E. IsomiR-SEA: An RNA-Seq analysis tool for miRNAs/isomiRs expression level profiling and miRNA-mRNA interaction sites evaluation. BMC Bioinformatics 2016, 17, 148. [Google Scholar] [CrossRef] [Green Version]
  32. Vlachos, I.S.; Paraskevopoulou, M.D.; Karagkouni, D.; Georgakilas, G.; Vergoulis, T.; Kanellos, I.; Anastasopoulos, I.-L.; Maniou, S.; Karathanou, K.; Kalfakakou, D. DIANA-TarBase v7.0: Indexing more than half a million experimentally supported miRNA:mRNA interactions. Nucleic Acids Res. 2015, 43, D153–D159. [Google Scholar] [CrossRef] [PubMed]
  33. Vlachos, I.S.; Zagganas, K.; Paraskevopoulou, M.D.; Georgakilas, G.; Karagkouni, D.; Vergoulis, T.; Dalamagas, T.; Hatzigeorgiou, A.G. DIANA-miRPath v3.0: Deciphering microRNA function with experimental support. Nucleic Acids Res. 2015, 43, W460–W466. [Google Scholar] [CrossRef] [PubMed]
  34. Pfaffl, M.W.; Tichopad, A.; Prgomet, C.; Neuvians, T.P. Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper—Excel-based tool using pair-wise correlations. Biotechnol. Lett. 2004, 26, 509–515. [Google Scholar] [CrossRef] [PubMed]
  35. Piórkowska, K.; Oczkowicz, M.; Rózycki, M.; Ropka-Molik, K.; Kajtoch, A.P.K. Novel porcine housekeeping genes for real-time RT-PCR experiments normalization in adipose tissue: Assessment of leptin mRNA quantity in different pig breeds. Meat Sci. 2011, 87, 191–195. [Google Scholar] [CrossRef] [PubMed]
  36. Piórkowska, K.; Tyra, M.; Ropka-Molik, K.; Podbielska, A. Evolution of peroxisomal trans-2-enoyl-CoA reductase (PECR) as candidate gene for meat quality. Livest. Sci. 2017, 201, 85–91. [Google Scholar] [CrossRef]
  37. Jacyno, E.; Pietruszka, A.; Kawecka, M.; Biel, W.; Kołodziej-Skalska, A. Phenotypic correlations of backfat thickness with meatiness traits, intramuscular fat, Longissimus muscle cholesterol and fatty acid composition in pigs. South Afr. J. Anim. Sci. 2015, 45, 122–128. [Google Scholar] [CrossRef] [Green Version]
  38. Paton, C.M.; Ntambi, J.M. Biochemical and physiological function of stearoyl-CoA desaturase. Am. J. Physiol. Endocrinol. Metab. 2009, 297, E28. [Google Scholar] [CrossRef] [Green Version]
  39. Hulver, M.W.; Berggren, J.R.; Carper, M.J.; Miyazaki, M.; Ntambi, J.M.; Hoffman, E.P.; Thyfault, J.P.; Stevens, R.; Dohm, G.L.; Houmard, J.A.; et al. Elevated stearoyl-CoA desaturase-1 expression in skeletal muscle contributes to abnormal fatty acid partitioning in obese humans. Cell Metab. 2005, 2, 251–261. [Google Scholar] [CrossRef] [Green Version]
  40. Flowers, M.T.; Ntambi, J.M. Role of stearoyl-coenzyme A desaturase in regulating lipid metabolism. Curr. Opin. Lipidol. 2008, 19, 248–256. [Google Scholar] [CrossRef] [Green Version]
  41. Cohen, P.; Miyazaki, M.; Socci, N.D.; Hagge-Greenberg, A.; Liedtke, W.; Soukas, A.A.; Sharma, R.; Hudgins, L.C.; Ntambi, J.M.; Friedman, J.M. Role for stearoyl-CoA desaturase-1 in leptin-mediated weight loss. Science 2002, 297, 240–243. [Google Scholar] [CrossRef]
  42. Yang, B.; Zhang, W.; Zhang, Z.; Fan, Y.; Xie, X.; Ai, H.; Ma, J.; Xiao, S.; Huang, L.; Ren, J. Genome-Wide Association Analyses for Fatty Acid Composition in Porcine Muscle and Abdominal Fat Tissues. PLoS ONE 2013, 8, e65554. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Puig-Oliveras, A.; Revilla, M.; Castelló, A.; Fernández, A.I.; Folch, J.M.; Ballester, M. Expression-based GWAS identifies variants, gene interactions and key regulators affecting intramuscular fatty acid content and composition in porcine meat. Sci. Rep. 2016, 6, 31803. [Google Scholar] [CrossRef] [PubMed]
  44. Matoušková, P.; Hanousková, B.; Skálová, L. Micrornas as potential regulators of glutathione peroxidases expression and their role in obesity and related pathologies. Int. J. Mol. Sci. 2018, 19, 1199. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Jin, M.; Wu, Y.; Wang, J.; Chen, J.; Huang, Y.; Rao, J.; Feng, C. MicroRNA-24 promotes 3T3-L1 adipocyte differentiation by directly targeting the MAPK7 signaling. Biochem. Biophys. Res. Commun. 2016, 474, 76–82. [Google Scholar] [CrossRef] [PubMed]
  46. Piórkowska, K.; Małopolska, M.; Ropka-Molik, K.; Nędza, M.S.; Wiechniak, A.; Żukowski, K.; Lambert, B.; Tyra, M. Evaluation of scd, acaca and fasn mutations: Effects on pork quality and other production traits in pigs selected based on RNA-seq results. Animals 2020, 10, 123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Pena, R.N.; Noguera, J.L.; García-Santana, M.J.; González, E.; Tejeda, J.F.; Ros-Freixedes, R.; Ibáñez-Escriche, N. Five genomic regions have a major impact on fat composition in Iberian pigs. Sci. Rep. 2019, 9, 1–9. [Google Scholar] [CrossRef] [Green Version]
  48. Muñoz, M.; Alves, E.; Corominas, J.; Folch, J.M.; Casellas, J.; Noguera, J.L.; Silió, L.; Fernández, A.I. Survey of SSC12 regions affecting fatty acid composition of intramuscular fat using high-density SNP data. Front. Genet. 2012, 2, 101. [Google Scholar] [CrossRef] [Green Version]
  49. Kim, S.-W.; Choi, Y.-I.; Choi, J.-S.; Kim, J.-J.; Choi, B.-H.; Kim, T.-H.; Kim, K.-S. Porcine Fatty Acid Synthase Gene Polymorphisms Are Associated with Meat Quality and Fatty Acid Composition. Korean J. Food Sci. Anim. Resour. 2011, 31, 356–365. [Google Scholar] [CrossRef]
  50. Grzes, M.; Sadkowski, S.; Rzewuska, K.; Szydlowski, M.; Switonski, M. Pig fatness in relation to FASN and INSIG2 genes polymorphism and their transcript level. Mol. Biol. Rep. 2016, 43, 381–389. [Google Scholar] [CrossRef] [Green Version]
  51. Brennan, A.M.; Mantzoros, C.S. Drug Insight: The role of leptin in human physiology and pathophysiology—Emerging clinical applications. Nat. Clin. Pract. Endocrinol. Metab. 2006, 2, 318–327. [Google Scholar] [CrossRef]
  52. Lönnqvist, F.; Arner, P.; Nordfors, L.; Schalling, M. Overexpression of the obese (ob) gene in adipose tissue of human obese subjects. Nat. Med. 1995, 1, 950–953. [Google Scholar] [CrossRef] [PubMed]
  53. Aslam, M. Leptin: Fights against obesity! Pak. J. Physiol. 2006, 2, 54–60. [Google Scholar]
  54. Hill, J.O.; Wyatt, H.R.; Peters, J.C. Energy balance and obesity. Circulation 2012, 126, 126–132. [Google Scholar] [CrossRef] [PubMed]
  55. Ridderstråle, M.; Johansson, L.E.; Rastam, L.; Lindblad, U. Increased risk of obesity associated with the variant allele of the PPARGC1A Gly482Ser polymorphism in physically inactive elderly men. Diabetologia 2006, 49, 496–500. [Google Scholar] [CrossRef] [Green Version]
  56. Liu, T.; Zhao, Y.; Tang, N.; Feng, R.; Yang, X.; Lu, N.; Wen, J.; Li, L. Pax6 Directly Down-Regulates Pcsk1n Expression Thereby Regulating PC1/3 Dependent Proinsulin Processing. PLoS ONE 2012, 7, e46934. [Google Scholar] [CrossRef]
  57. Shi, M.A.; Shi, G.P. Different roles of mast cells in obesity and diabetes: Lessons from experimental animals and humans. Front. Immunol. 2012, 3, 7. [Google Scholar] [CrossRef] [Green Version]
  58. Finlin, B.S.; Confides, A.L.; Zhu, B.; Boulanger, M.C.; Memetimin, H.; Taylor, K.W.; Johnson, Z.R.; Westgate, P.M.; Dupont-Versteegden, E.E.; Kern, P.A. Adipose Tissue Mast Cells Promote Human Adipose Beiging in Response to Cold. Sci. Rep. 2019, 9, 8658. [Google Scholar] [CrossRef]
  59. Ruiz-Ojeda, F.J.; Méndez-Gutiérrez, A.; Aguilera, C.M.; Plaza-Díaz, J. Extracellular matrix remodeling of adipose tissue in obesity and metabolic diseases. Int. J. Mol. Sci. 2019, 20, 4888. [Google Scholar] [CrossRef] [Green Version]
  60. Lin, D.; Chun, T.H.; Kang, L. Adipose extracellular matrix remodelling in obesity and insulin resistance. Biochem. Pharmacol. 2016, 119, 8–16. [Google Scholar] [CrossRef] [Green Version]
  61. Huang, Y.; Yan, Y.; Xv, W.; Qian, G.; Li, C.; Zou, H.; Li, Y. A New Insight into the Roles of MiRNAs in Metabolic Syndrome. BioMed Res. Int. 2018. [Google Scholar] [CrossRef] [Green Version]
  62. Heneghan, H.M.; Miller, N.; McAnena, O.J.; O’brien, T.; Kerin, M.J. Differential miRNA expression in omental adipose tissue and in the circulation of obese patients identifies novel metabolic biomarkers. J. Clin. Endocrinol. Metab. 2011, 96, E846–E850. [Google Scholar] [CrossRef] [PubMed]
  63. Divoux, A.; Tordjman, J.; Le Lacasa, D.; Veyrie, N.; Hugol, D.; Aissat, A.; Basdevant, A.; Le Guerre-Millo, M.; Poitou, C.; Zucker, J.-D.; et al. Fibrosis in Human Adipose Tissue: Composition, Distribution, and Link With Lipid Metabolism and Fat Mass Loss. Am. Diabetes Assoc. 2010, 59, 2817–2825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Pellegrinelli, V.; Heuvingh, J.; Du Roure, O.; Rouault, C.; Devulder, A.; Klein, C.; Lacasa, M.; Clément, E.; Lacasa, D.; Clément, K. Human adipocyte function is impacted by mechanical cues. J. Pathol. 2014, 233, 183–195. [Google Scholar] [CrossRef] [PubMed]
  65. Rogero, M.M.; Calder, P.C. Obesity, Inflammation, Toll-Like Receptor 4 and Fatty Acids. Nutrients 2018, 10, 432. [Google Scholar] [CrossRef] [Green Version]
  66. Ahmad, R.; Al-Mass, A.; Atizado, V.; Al-Hubail, A.; Al-Ghimlas, F.; Al-Arouj, M.; Bennakhi, A.; Dermime, S.; Behbehani, K. Elevated expression of the toll like receptors 2 and 4 in obese individuals: Its significance for obesity-induced inflammation. J. Inflamm. 2012, 9, 48. [Google Scholar] [CrossRef] [Green Version]
  67. Song, M.J.; Kim, K.H.; Yoon, J.M.; Kim, J.B. Activation of Toll-like receptor 4 is associated with insulin resistance in adipocytes. Biochem. Biophys. Res. Commun. 2006, 346, 739–745. [Google Scholar] [CrossRef]
  68. Deiuliis, J. MicroRNAs as regulators of metabolic disease: Pathophysiologic significance and emerging role as biomarkers and therapeutics. Int. J. Obes. Vol. 2016, 40, 88–101. [Google Scholar] [CrossRef] [Green Version]
  69. Matsuo, Y.; Tanaka, M.; Yamakage, H.; Metabolism, Y.S.; Satoh-Asahara, N. Thrombospondin 1 as a novel biological marker of obesity and metabolic syndrome. Metabolism 2015, 64, 1490–1499. [Google Scholar] [CrossRef] [Green Version]
Genes 11 00600 g001 550
Figure 1. Comparison of differentially-expressed genes (DEGs) between subcutaneous fat tissue varying in thickness in each of three breeds (LW, Large White; Pi, Pietrain; and Hp, Hampshire).
Figure 1. Comparison of differentially-expressed genes (DEGs) between subcutaneous fat tissue varying in thickness in each of three breeds (LW, Large White; Pi, Pietrain; and Hp, Hampshire).
Genes 11 00600 g001
Genes 11 00600 g002 550
Figure 2. The significant enrichment GO terms and pathways identified based on backfat thickness differences and related to lipid metabolism (WebGestalt software).
Figure 2. The significant enrichment GO terms and pathways identified based on backfat thickness differences and related to lipid metabolism (WebGestalt software).
Genes 11 00600 g002
Genes 11 00600 g003 550
Figure 3. The detected differentially-expressed genes involved in GO terms related to lipid metabolism processes (String software based on Sus scrofa reference; detected genes showed no more than five interactions. Line shape indicates the predicted mode of action: red, interactions that were experimentally determined; blue, interactions from curated databases; black, co-expression; green, text mining associations and interactions based on relevant publications mentioning a transfer from other organisms; yellow, transcriptional regulation).
Figure 3. The detected differentially-expressed genes involved in GO terms related to lipid metabolism processes (String software based on Sus scrofa reference; detected genes showed no more than five interactions. Line shape indicates the predicted mode of action: red, interactions that were experimentally determined; blue, interactions from curated databases; black, co-expression; green, text mining associations and interactions based on relevant publications mentioning a transfer from other organisms; yellow, transcriptional regulation).
Genes 11 00600 g003
Genes 11 00600 g004 550
Figure 4. Venn diagram of differentially-expressed miRNAs in subcutaneous fat tissue dependent on fat content in three analyzed breeds (LW, Large White; Pi, Pietrain; and Hp, Hampshire).
Figure 4. Venn diagram of differentially-expressed miRNAs in subcutaneous fat tissue dependent on fat content in three analyzed breeds (LW, Large White; Pi, Pietrain; and Hp, Hampshire).
Genes 11 00600 g004
Genes 11 00600 g005 550
Figure 5. Venn diagram of differentially-expressed miRNAs in subcutaneous fat tissue dependent on fat content in three analyzed breeds (LW, Large White; Pi, Pietrain; and Hp, Hampshire).
Figure 5. Venn diagram of differentially-expressed miRNAs in subcutaneous fat tissue dependent on fat content in three analyzed breeds (LW, Large White; Pi, Pietrain; and Hp, Hampshire).
Genes 11 00600 g005
Table
Table 1. The differences in fatness traits detected in all three breeds and obtained pig groups used in genetic analyses.
Table 1. The differences in fatness traits detected in all three breeds and obtained pig groups used in genetic analyses.
Backfat Thickness (cm)Weight of Peritoneal Fat (kg)
L H LH
Pietrain0.71±0.0.1 b1.29±0.10 a0.15±0.001 b0.31±0.051 a
Hampshire0.99±0.03 b1.54±0.20 a0.30±0.0160.29±0.067
Large White1.04±0.061.40±0.140.23±0.0700.31±0.122
L, low fatness; H, high fatness. Data are presented as means ± standard error; means with different letters (a,b differ significantly, p-value = 0.05).
Table
Table 2. Significant, enriched Gene Ontology terms for differentially-expressed genes (DEGs) in relation to backfat thickness.
Table 2. Significant, enriched Gene Ontology terms for differentially-expressed genes (DEGs) in relation to backfat thickness.
Gene Ontology/Accession NumberFDR *NGene Name
Extracellular matrix organization (GO:0030198)9.2 × 10−36CYR61; ELN; GFAP; HPSE; POSTN; VTN
Positive regulation of mast cell degranulation (GO:0043306)2.4 × 10−34FGR; FCER1A; FCER1G; ZAP70
Cell adhesion (GO:0007155)3.5 × 10−312TNFAIP6; WISP1; CTGF; CYR61; HAS1; LYVE1;
NOV; POSTN; SELL; TNC; THBS2; THBS3
Innate immune response (GO:0045087)5.9 × 10−214FGR; FCER1G; MX1; MX2; S100A8; B2M; BST2; FGB; IFIH1; JCHAIN; LCN2; PML; TLR2; ZAP70;
Fatty acid biosynthetic process (GO:0006633)8.4 × 10−34ACACA; SCD; FASN; ACLY
Positive regulation of apoptotic process (GO:0043065)5.6 × 10−310BMF; ALDH1A2; CLU; CYP1B1; GADD45G; IGFBP3; SFRP2; TOP2A; TGM2; ZBTB16
Long-chain fatty acid biosynthetic process (GO:0042759)6.0 × 10−33PLPP1; SCD; SCD5
Response to dietary excess (GO:0002021)8.6 × 10−33PPARGC1A; LEP; PCSK1N
Positive regulation of fat cell differentiation (GO:0045600)6.8 × 10−35CCDC3; MEDAG; SFRP2; ZBTB16; ZNF385A
Inflammatory response (GO:0006954)8.6 × 10−311CCL5; CCR5; CD180; FAS; CXCL10; C5AR1; PLP1; TSPAN2; TBXA2R; TLR2; ZAP70
* FDR (false discovery rate); p-values are shown after Benjamini correction and according to David software; N, number of identified DEGs.−
Table
Table 3. The fold-change values of DEGs related to lipid metabolism and detected as significant for at least two breeds.
Table 3. The fold-change values of DEGs related to lipid metabolism and detected as significant for at least two breeds.
GeneAccession NumberPietrainLarge WhiteHampshire
FCFDRFCFDRFCFDR
LEPENSSSCG000000404641.410.032.140.0011.10ns
ACACAENSSSCG000000176941.620.031.670.0011.52ns
SCDENSSSCG000000105542.950.0011.660.0012.230.04
SCD5ENSSSCG00000009245−1.320.04−1.09ns−2.280.001
FASNENSSSCG000000299441.58ns1.470.023.070.005
ACOX3ENSSSCG000000087241.400.051.32ns2.060.001
C2ENSSSCG00000001422−3.180.001−2.000.0012.310.01
ACLYENSSSCG000000174211.650.052.830.0011.13ns
TNCENSSSCG00000005494−3.230.0004−2.550.001−4.17ns
PPARGC1AENSSSCG000000292751.16ns13.450.00013.270.05
PCSK1NENSSSCG000000213281.070.05−1.52ns−2.990.05
TLR2ENSSSCG00000009002−2.630.001−1.590.05−1.34ns
FGRENSSSCG000000035781.130.051.620.001−1.10ns
FCER1AENSSSCG000000064131.400.001−1.190.01−2.69ns
FCER1GENSSSCG00000006357−2.400.001−1.400.05−1.22ns
FDR (false discovery rate); FC, fold change; ns, not significant.
Table
Table 4. Significant, enriched Gene Ontology terms for differentially-expressed miRNAs in relation to backfat thickness.
Table 4. Significant, enriched Gene Ontology terms for differentially-expressed miRNAs in relation to backfat thickness.
Gene Ontology/Accession NumberFDR *N Target GenesN miRNAsmiRNAs
Extracellular matrix organization (GO:0030198)<1.0 × 10−3251135hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; hsa-miR-145-5p; hsa-let-7d-5p
Extracellular matrix disassembly (GO:0022617)<1.0 × 10−325466hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-26a-5p; hsa-miR-145-5p; hsa-let-7d-5p
Cellular lipid metabolic process (GO:0044255)3.1 × 10−13638hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-26a-5p; has-miR-143-3p; has-miR-142-5p; has-miR-145-5p; hsa-let-7d-5p
Cell junction organization (GO:0034330)<1.0 × 10−325666hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-143-3p; hsa-miR-145-5p; hsa-let-7d-5p
MyD88-independent toll-like receptor signaling pathway (GO:0002756)<1.0 × 10−325447hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-26a-5p; hsa-miR-145-5p; hsa-let-7d-5p; has-miR-378a-3p
Cellular component disassembly involved in execution phase of apoptosis (GO:0006921)<1.0 × 10−325287hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-26a-5p; hsa-miR-145-5p; hsa-miR-145-5p; hsa-miR-378a-3p
Toll-like receptor TLR1:TLR2 signaling pathway (GO:0038123)Toll-like receptor TLR6:TLR2 signaling pathway (GO:0038124) <1.0 × 10−325397hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-26a-5p; hsa-miR-139-5p; has-miR-378a-3p; hsa-let-7d-5p
Innate immune response (GO:0045087)<1.0 × 10−3252137hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-24-3p; has-miR-26a-5p; hsa-miR-139-5p; has-miR-142-5p; has-miR-145-5p; hsa-let-7d-5p
* p-values are presented after FDR correction and according to the mirPath v.3 tool; N target genes, the number of predicted target genes of differentially-expressed miRNAs; N miRNAs, number of DE miRNAs.
Table
Table 5. Significant pathways detected for both DEGs (significant for at least two breeds) and DE miRNAs (14 miRNAs common for all breeds).
Table 5. Significant pathways detected for both DEGs (significant for at least two breeds) and DE miRNAs (14 miRNAs common for all breeds).
DE miRNAsDEGs
PathwaysFDR *NmiRNAsNN Target GenesFDR*NGenes
Fatty acid metabolism (hsa01212/ssc01212)6.6 × 10−164hsa-miR-100-5p; hsa-miR-10b-5p; hsa-miR-143-3p; hsa-miR-24-3p9FASN; ACACA;MCAT; ACAA2; ACAA1; CPT1A; CPT2; SCD; ELOVL55.4 × 10−35ACACA;ACOX3; FASN; SCD5; SCD
ECM−receptor interaction (hsa04512/ssc04512)<1.0 × 10−3254hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-143-3p; hsa-miR-145-5p35SPP1; CD44; LAMC1; LAMC3; ITGA3; ITGA5; FN1; TNC; COL1A2; COL4A2; COL5A1; COL5A2; COL6A1; COL27A1; ITGB1; GP5; THBS12.4 × 10−211CHAD; COL1A1; COL5A3; COL6A6; COL11A1; LAMC2; TNC; THBS2; THBS3; THBS4; VTN
Hippo-signaling pathway (hsa04390/ssc04390)<1.0 × 10−3257hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-145-5p; hsa-miR-24-3p; hsa-miR-139-5p; hsa-let-7d-5p; hsa-miR-26a-5p87PPP2CA; BMP5; BMP2; BMP7; FGF1; ACTG1; PPP1CB; CTGF;
Figure S2		0.017YWHAH; CTGF; FZD2; FZD4; NKD1; TGFB3; WNT2B
Fatty acid biosynthesis (hsa00061/ssc00061)<1.0 × 10−3254hsa-miR-100-5p; hsa-miR-10b-5p; hsa-miR-143-3p; hsa-miR-24-3p3FASN; ACACA; MCATns2FASN; ACACA
Cell cycle (hsa04110/ssc04110)<1.0 × 10−3259hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-10b-5p; hsa-miR-143-3p; hsa-miR-24-3p; hsa-miR-26a-5p; hsa-miR-142-5p; hsa-let-7d-5p78ESPL1; CDC6; GSK3B; CCNB2; RBL2; PCNA; E2F1E2F2;YWHAE; MCM6; MCM2
Figure S3		0.0035E2F2; GADD45G; MCM2;
TGFB3; YWHAH;
P53-signaling pathways (hsa04115/ ssc04151)2.5 × 10−9 7 hsa-let-7a-5p; hsa-let-7i-5p; hsa-miR-26a-5p; hsa-miR-143-3p; hsa-miR-10b-5p; hsa-miR-142-5p; hsa-miR-378-3p60ZMAT3; CCNB2; CCNB1;
CDK4; BID; THBS1; CDK2
CCND2; PERP; RRM2B; CDK1; CDKN2A; CDK6
CHEK1; TP53; APAF1; PMAIP1; CD82; CASP3		0.000213CHAD; COL1A1; COL6A6; LAMC2; PCK1; THBS2; THBS3; THBS4; TLR2; TNC; VTN; YWHAH
Accession numbers are presented for DE miRNAs/DEGs from the mirPath v.3 and KEGG database, respectively; * p-values are presented after FRD correction and according to the mirPath v.3 tool for DE miRNAs and after Benjamini correction and according to David software for DEGs; N, number of identified DEGs; ns, not significant; highlighted genes were identified both as differentially-expressed and as predicted genes regulated by miRNAs. The bold and underline were used to highlighted the genes identified in both–DEGs and targeted genes groups.
Table
Table 6. Correlation coefficients for NGS results and qPCR data for both miRNAs and DEGs.
Table 6. Correlation coefficients for NGS results and qPCR data for both miRNAs and DEGs.
DEGsCorrelationmiRNAsCorrelation
ROCK10.56hsa-miR-26a-5p0.81 *
LRP120.66 *hsa-let-7a-5p0.50
ACACA0.91 **hsa-mir-100-5p0.75 *
HK20.87 **hsa-mir-378a-3p0.40
LRP60.75 *hsa-mir-103a-3p0.88 *
LEP0.94 **hsa-miR-125b-5p0.73
TNC0.83 *
PCK10.95 **
* p-value < 0.05; ** p < 0.01.

Share and Cite

MDPI and ACS Style

Ropka-Molik, K.; Pawlina-Tyszko, K.; Żukowski, K.; Tyra, M.; Derebecka, N.; Wesoły, J.; Szmatoła, T.; Piórkowska, K. Identification of Molecular Mechanisms Related to Pig Fatness at the Transcriptome and miRNAome Levels. Genes 2020, 11, 600. https://doi.org/10.3390/genes11060600

AMA Style

Ropka-Molik K, Pawlina-Tyszko K, Żukowski K, Tyra M, Derebecka N, Wesoły J, Szmatoła T, Piórkowska K. Identification of Molecular Mechanisms Related to Pig Fatness at the Transcriptome and miRNAome Levels. Genes. 2020; 11(6):600. https://doi.org/10.3390/genes11060600

Chicago/Turabian Style

Ropka-Molik, Katarzyna, Klaudia Pawlina-Tyszko, Kacper Żukowski, Mirosław Tyra, Natalia Derebecka, Joanna Wesoły, Tomasz Szmatoła, and Katarzyna Piórkowska. 2020. "Identification of Molecular Mechanisms Related to Pig Fatness at the Transcriptome and miRNAome Levels" Genes 11, no. 6: 600. https://doi.org/10.3390/genes11060600

APA Style

Ropka-Molik, K., Pawlina-Tyszko, K., Żukowski, K., Tyra, M., Derebecka, N., Wesoły, J., Szmatoła, T., & Piórkowska, K. (2020). Identification of Molecular Mechanisms Related to Pig Fatness at the Transcriptome and miRNAome Levels. Genes, 11(6), 600. https://doi.org/10.3390/genes11060600

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