Differences in gene expression profiles for subcutaneous adipose, liver, and skeletal muscle tissues between Meishan and Landrace pigs with different backfat thicknesses

Comparison of Growth Performance, Carcass Properties, Fatty Acid Profile, and Genes Involved in Fat Metabolism in Nanyang and Landrace Pigs

This study compared the growth, carcass properties, fatty acid profile, lipid-producing enzyme activity, and expression pattern of genes involved in fat metabolism in Nanyang and Landrace pigs. In the study, 32 Nanyang (22.16 ± 0.59 kg) and 32 Landrace barrows (21.37 ± 0.57 kg) were selected and divided into two groups, each with eight pens and four pigs per pen. The trial period lasted 90 days. The findings showed that the Nanyang pigs had lower average daily weight gain and lean percentage and higher average backfat thickness and lipogenic enzyme activities, including for acetyl-CoA carboxylase, glucose-6-phosphate dehydrogenase, malic enzyme, and fatty acid synthase, than the Landrace pigs. A total of 14 long-chain fatty acids were detected using HPLC-MS, in which it was found that the levels of C14:0, C18:1n-9, C20:1n-9, C20:4n-6, and MUFA were up-regulated and C18:2n-6, C18:3n-3, PUFA n6, n3/n6, and total PUFA were down-regulated in the Nanyang pigs. Moreover, the mRNA levels for genes involved in fat metabolism, ME1, FAS, and LPL, were higher and the expression of SREBP1 mRNA was lower in the Nanyang pigs. Our results suggest genetic differences between the pig breeds in terms of growth, carcass traits, lipogenic enzyme activities, fatty acid profile, and the mRNA expression of genes involved in fat metabolism in subcutaneous fat tissue, which may provide a basis for high-quality pork production. Further studies are needed to investigate the regulation of lipid metabolism.

An Analysis of a Transposable Element Expression Atlas during 27 Developmental Stages in Porcine Skeletal Muscle: Unveiling Molecular Insights into Pork Production Traits

The development and growth of porcine skeletal muscle determine pork quality and yield. While genetic regulation of porcine skeletal muscle development has been extensively studied using various omics data, the role of transposable elements (TEs) in this context has been less explored. To bridge this gap, we constructed a comprehensive atlas of TE expression throughout the developmental stages of porcine skeletal muscle. This was achieved by integrating porcine TE genomic coordinates with whole-transcriptome RNA-Seq data from 27 developmental stages. We discovered that in pig skeletal muscle, active Tes are closely associated with active epigenomic marks, including low levels of DNA methylation, high levels of chromatin accessibility, and active histone modifications. Moreover, these TEs include 6074 self-expressed TEs that are significantly enriched in terms of muscle cell development and myofibril assembly. Using the TE expression data, we conducted a weighted gene co-expression network analysis (WGCNA) and identified a module that is significantly associated with muscle tissue development as well as genome-wide association studies (GWAS) of the signals of pig meat and carcass traits. Within this module, we constructed a TE-mediated gene regulatory network by adopting a unique multi-omics integration approach. This network highlighted several established candidate genes associated with muscle-relevant traits, including HES6, CHRNG, ACTC1, CHRND, MAMSTR, and PER2, as well as novel genes like ENSSSCG00000005518, ENSSSCG00000033601, and PIEZO2. These novel genes hold promise for regulating muscle-related traits in pigs. In summary, our research not only enhances the TE-centered dissection of the genetic basis underlying pork production traits, but also offers a general approach for constructing TE-mediated regulatory networks to study complex traits or diseases.

Genome wide copy number variations using Porcine 60K SNP Beadchip in Landlly pigs

In the present study, Porcine 60K SNP genotype data from 69 Landlly pigs were used to explore Copy Number Variations (CNVs) across the autosomes. A total of 386 CNVs were identified using Hidden Markov Model (HMM) in PennCNV software, which were subsequently aggregated to 115 CNV regions (CNVRs). Among the total detected CNVRs, 58 gain, 49 were loss type while remaining 8 events were both gain and loss types. Identified CNVRs covered 12.5 Mb (0.55%) of Sus scrofa reference 11.1 genome. Comparison of our results with previous investigations on pigs revealed that approximately 75% CNVRs were novel, which may be due to differences in genetic background, environment and implementation of artificial selection in Landlly pigs. Functional annotation and pathway analysis showed the significant enrichment of 267 well-annotated Sus scrofa genes in CNVRs. These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity. Additionally, KIT and FUT1 were two major genes detected on CNVR in our population. This investigation provided a comprehensive overview of CNV distribution in the Indian porcine genome for the first time, which may be useful for further investigating the association of important quantitative traits in Landlly pigs.Highlights115 CNVRs were identified in 69 Landlly pig population.Approximately 75% detected CNVRs were novel for Landlly population.Significant enrichment of 267 well-annotated Sus scrofa genes observed in these CNVRs.These genes were involved in different biological functions like sensory perception, meat quality traits, back fat thickness and immunity.Comprehensive CNV map in the Indian porcine genome developed for the first time.

Expression and Genetic Effects of GLI Pathogenesis-Related 1 Gene on Backfat Thickness in Pigs

Backfat thickness (BFT) is an important carcass composition trait and regarded as a breeding focus. Our initial transcriptome analysis of pig BFT identified GLI pathogenesis-related 1 (GLIPR1) as one of the promising candidate genes. This study was conducted to identify the expression profiles, polymorphisms, and genetic effects of the GLIPR1 gene on BFT in pigs. The expression of the GLIPR1 gene existed in every detected tissue, and there was a significantly higher expression in spleen and adipose tissue than others (p < 0.05). At the different ages of pig, the expression of the GLIPR1 gene was low at an early age, increased with growth, and reached the highest level at 180 days. Genetic polymorphism analysis was detected in 553 individuals of the Large White × Minzhu F2 population. Four SNPs in the promoter significantly associated with 6−7 rib BFT (p < 0.05) were predicted to alter the transcription factor binding sites (TFBS), and the mutations of g.38758089 T>G and g.38758114 G>C were predicted to change the TFs associated with the regulation of adipogenesis. Haplotypes were formed by the detected SNPs, and one block showed a strong association with BFT (p < 0.05). In summary, our results indicate that the expression profiles and genetic variants of GLIPR1 affected the BFT of pigs. To our knowledge, this study is the first to demonstrate the biological function and genetic effects of the GLIPR1 gene on the BFT of pig and provide genetic markers to optimize breeding for BFT in pigs.

Single-Step Genome Wide Association Study Identifies QTL Signals for Untrimmed and Trimmed Thigh Weight in Italian Crossbred Pigs for Dry-Cured Ham Production

Simple Summary

Along with the traditional traits, swine breeding programs for Italian dry-cured ham production have recently aimed to include novel phenotypes. The identification of the genomic regions underlying such new traits helps to untangle their genetic architecture and may provide useful information to be integrated in genetic selection. With this aim, we estimated genetic parameters and conducted a single step genome wide association studies (GWAS) on untrimmed and trimmed thigh weight considering two pig crossbred lines approved for Italian Protected Designation of Origin ham production. Quantitative trait loci (QTLs) were characterized based on the variance of 10-SNP sliding windows genomic estimated breeding values. In particular, we identified interesting QTL signals on several chromosomes, notably on chromosome 4, 6, 7 and 15. A high heritability and genetic correlation were observed for the two traits under investigation and although independent studies including other pig populations are required to disentangle the possible effects of specific linkage disequilibrium in our population, our findings suggest that such QTL could be investigated in future pig breeding programs to improve the reliability of genomic estimated breeding values for the dry-cured ham production.

Abstract

Protected Designation of Origin (PDO) dry-cured ham is the most important product in the Italian pig breeding industry, mainly oriented to produce heavy pig carcasses to obtain hams of the right weight and maturity. Recently, along with the traditional traits swine breeding programs have aimed to include novel carcass traits. The identification at the genome level of quantitative trait loci (QTLs) affecting such new traits helps to reveal their genetic determinism and may provide information to be integrated in prediction models in order to improve prediction accuracy as well as to identify candidate genes underlying such traits. This study aimed to estimate genetic parameters and perform a single step genome wide association studies (ssGWAS) on novel carcass traits such as untrimmed (UTW) and trimmed thigh weight (TTW) in two pig crossbred lines approved for the ham production of the Italian PDO. With this purpose, phenotypes were collected from ~1800 animals and 240 pigs were genotyped with Illumina PorcineSNP60 Beadchip. The single-step genomic BLUP procedure was used for the heritability estimation and to implement the ssGWAS. QTL were characterized based on the variance of 10-SNP sliding window genomic estimated breeding values. Moderate heritabilities were detected and QTL signals were identified on chromosome 1, 4, 6, 7, 11 and 15 for both traits. As expected, the genetic correlation among the two traits was very high (~0.99). The QTL regions encompassed a total of 249 unique candidate genes, some of which were already reported in association with growth, carcass or ham weight traits in pigs. Although independent studies are required to further verify our findings and disentangle the possible effects of specific linkage disequilibrium in our population, our results support the potential use of such new QTL information in future breeding programs to improve the reliability of genomic prediction.

A compendium and comparative epigenomics analysis of cis-regulatory elements in the pig genome

Although major advances in genomics have initiated an exciting new era of research, a lack of information regarding cis-regulatory elements has limited the genetic improvement or manipulation of pigs as a meat source and biomedical model. Here, we systematically characterize cis-regulatory elements and their functions in 12 diverse tissues from four pig breeds by adopting similar strategies as the ENCODE and Roadmap Epigenomics projects, which include RNA-seq, ATAC-seq, and ChIP-seq. In total, we generate 199 datasets and identify more than 220,000 cis-regulatory elements in the pig genome. Surprisingly, we find higher conservation of cis-regulatory elements between human and pig genomes than those between human and mouse genomes. Furthermore, the differences of topologically associating domains between the pig and human genomes are associated with morphological evolution of the head and face. Beyond generating a major new benchmark resource for pig epigenetics, our study provides basic comparative epigenetic data relevant to using pigs as models in human biomedical research.

miR-32-5p Regulates Lipid Accumulation in Intramuscular Fat of Erhualian Pigs by Suppressing KLF3

Intramuscular fat (IMF) and subcutaneous fat (SCF) are important traits affecting the economics of the pork industry, in which less SCF and more IMF content is desirable. However, the mechanisms that regulate IMF and SCF content are not clear yet. In this study, we demonstrate that KLF3 (Krüppel-like factor 3) was negatively correlated with IMF content in the longissimus dorsi muscle of Erhualian pigs. In addition, the expression level of KLF3 was significantly higher in IMF than SCF. Overexpression and knockdown experiments revealed that KLF3 could suppress adipocyte differentiation in vitro by downregulating adipogenic markers, including PPARG, C/EBPA, and FABP4. Luciferase activity analysis proved that miR-32-5p was able to suppress KLF3. Notably, miR-32-5p level was negatively correlated to KLF3 mRNA level in both IMF and SCF tissues. The same relationship was proved in samples with different IMF content. Further studies showed that miR-32-5p could promote adipocyte differentiation via inhibiting KLF3. Our results suggest that the miR-32-5p-KLF3 pathway is involved in the regulation of differential fat deposition of IMF and SCF tissues.

Transcriptional Profiles of Skeletal Muscle Associated With Increasing Severity of White Striping in Commercial Broilers

Development of the white striping (WS) abnormality adversely impacts overall quality of broiler breast meat. Its etiology remains unclear. This study aimed at exploring transcriptional profiles of broiler skeletal muscles exhibiting different WS severity to elucidate molecular mechanisms underlying the development and progression of WS. Total RNA was isolated from pectoralis major of male 7-week-old Ross 308 broilers. The samples were classified as mild (n = 6), moderate (n = 6), or severe (n = 4), based on number and thickness of the white striations on the meat surface. The transcriptome was profiled using a chicken gene expression microarray with one-color hybridization technique. Gene expression patterns of each WS severity level were compared against each other; hence, there were three comparisons: moderate vs. mild (C1), severe vs. moderate (C2), and severe vs. mild (C3). Differentially expressed genes (DEGs) were identified using the combined criteria of false discovery rate ≤ 0.05 and absolute fold change ≥1.2. Differential expression of 91, 136, and 294 transcripts were identified in C1, C2, and C3, respectively. There were no DEGs in common among the three comparisons. Based on pathway analysis, the enriched pathways of C1 were related with impaired homeostasis of macronutrients and small biochemical molecules with disrupted Ca²⁺-related pathways. Decreased abundance of the period circadian regulator suggested the shifted circadian phase when moderate WS developed. The enriched pathways uniquely obtained in C2 were RNA degradation, Ras signaling, cellular senescence, axon guidance, and salivary secretion. The DEGs identified in those pathways might play crucial roles in regulating cellular ion balances and cell-cycle arrest. In C3, the pathways responsible for phosphatidylinositol 3-kinase-Akt signaling, p53 activation, apoptosis, and hypoxia-induced processes were modified. Additionally, pathways associated with a variety of diseases with the DEGs involved in regulation of [Ca²⁺], collagen formation, microtubule-based motor, and immune response were identified. Eight pathways were common to all three comparisons (i.e., calcium signaling, Ras-associated protein 1 signaling, ubiquitin-mediated proteolysis, vascular smooth muscle contraction, oxytocin signaling, and pathway in cancer). The current findings support the role of intracellular ion imbalance, particularly Ca²⁺, oxidative stress, and impaired programmed cell death on WS progression.

Effects of long-term feeding of rapeseed meal on skeletal muscle transcriptome, production efficiency and meat quality traits in Norwegian Landrace growing-finishing pigs

This study was performed to investigate the effects of dietary inclusion of 20% rapeseed meal (RSM) as an alternative to soybean meal (SBM) in a three-month feeding experiment with growing finishing pigs. Dietary alteration affected growth performance, several carcass traits and transcriptional responses in the skeletal muscle, but did not affect measured meat quality traits. In general, pigs fed the RSM test diet exhibited reduced growth performance compared to pigs on SBM control diet. Significant transcriptional changes in the skeletal muscle of growing pigs fed RSM diet were likely the consequence of an increased amount of fiber and higher polyunsaturated fatty acids, and presence of bioactive phytochemicals, such as glucosinolates. RNAseq pipeline using Tophat2-Cuffdiff identified 57 upregulated and 63 downregulated genes in RSM compared to SBM pigs. Significantly enriched among downregulated pathways was p53-mediated signalling involved in cellular proliferation, while activation of negative growth regulators (IER5, KLF10, BTG2, KLF11, RETREG1, PRUNE2) in RSM fed pigs provided further evidence for reduced proliferation and increased cellular death, in accordance with the observed reduction in performance traits. Upregulation of well-known metabolic controllers (PDK4, UCP3, ESRRG and ESRRB), involved in energy homeostasis (glucose and lipid metabolism, and mitochondrial function), suggested less available energy and nutrients in RSM pigs. Furthermore, several genes supported more pronounced proteolysis (ABTB1, OTUD1, PADI2, SPP1) and reduced protein synthesis (THBS1, HSF4, AP1S2) in RSM muscle tissue. In parallel, higher levels of NR4A3, PDK4 and FGF21, and a drop in adropin, ELOVL6 and CIDEC/FSP27 indicated increased lipolysis and fatty acid oxidation, reflective of lower dressing percentage. Finally, pigs exposed to RSM showed greater expression level of genes responsive to oxidative stress, indicated by upregulation of GPX1, GPX2, and TXNIP.

Comparing pig breeds with genetically low and high backfat thickness: differences in expression of adiponectin, its receptor, and blood metabolites

Here we characterized gene expressions in subcutaneous adipose tissue and blood metabolites of pigs with genetically low backfat (Landrace) and high backfat (Meishan). As pigs aged from 1 wk-to 3-mo old, mRNA levels of adipose-specific genes increased, although their gene expressions coding for major enzymes involved in lipid metabolism (lipoprotein lipase, fatty acid synthase, and hormone-sensitive lipase) did not differ between lean and fat pigs. Instead, there were significant effects for adiponectin and its receptor AdipoR1 mRNA levels between the two breeds of which respective expressions were lower and higher in Meishan by 3 mo of age. Contrary to changes in gene expressions, the concentrations of blood glucose, triglyceride (TG), and NEFA in both breeds decreased during growth, and 3-mo-old Meishan evidenced lower glucose with higher TG than the Landrace. The homeostasis model assessment insulin resistance (HOMA-IR) index was also calculated from the measurements of fasting glucose and insulin concentration, and Meishan showed a higher value than the Landrace. We next examined these differences in Landrace and Meishan crossbreds, which were phenotypically distinguishable by the backfat thickness as the former lean type and the latter fat type. As with the purebreds, high backfat Meishan crosses showed the characteristics of lower glucose and higher TG in circulating levels and also lower adiponectin transcripts in subcutaneous adipose tissue. Collectively, our results demonstrate that levels of adiponectin and its receptor gene expressions, blood glucose, blood lipids, and HOMA-IR in pigs vary between lean and fat. These observations strongly suggest the possibility that overall metabolic differences rather than adipocyte ability itself contribute to the fatness of genetically high backfat pigs.