Longissimus dorsi transcriptome analysis of purebred and crossbred Iberian pigs differing in muscle characteristics

Multi-omics reveals the mechanism of quality discrepancy between Gayal (Bos frontalis) and yellow cattle beef

BACKGROUND

Producing high-quality beef with enhanced muscle composition and reduced fat content is critical for meeting consumer preferences and supporting a balanced diet. Given the substantial variability in beef quality across cattle breeds, this study aimed to identify key determinants of meat quality by examining Gayal (Bos frontalis) and yellow cattle (Bos taurus) through a multi-disciplinary approach.

RESULTS

The results demonstrated that Gayal cattle exhibited superior meat quality, characterized by higher levels of protein, flavor-enhancing and essential amino acids, total amino acids, and polyunsaturated fatty acids (PUFAs), alongside reduced fat content, with similar trends observed in serum hormone and amino acid profiles. Distinct differences in gut microbial composition, enzymatic activities, and metabolites were observed between the breeds. Gayal displayed increased abundances of key bacterial taxa such as Akkermansia, Paeniclostridium, Escherichia-Shigella, and Clostridium sensu stricto 1, which were associated with enhanced volatile fatty acids (VFAs), ammoniacal nitrogen, and enzymatic activity in the colon. Transcriptomic analysis of the psoas major (PM) muscle revealed significant changes in genes linked to muscle development, amino acid metabolism, and lipid metabolism. Genes related to intestinal amino acid absorption were upregulated in Gayal, while those connected to short-chain fatty acid absorption were downregulated. Correlation analyses underscored the role of gut microbiota and metabolic profiles in modulating gene expression associated with lipid and amino acid metabolism, ultimately influencing meat flavor and quality.

CONCLUSIONS

These findings provide actionable insights into the genetic and microbial factors underlying beef quality, offering a foundation for enhancing local cattle resources, optimizing breeding programs, and advancing the production of premium beef to meet both market and dietary needs.

IGFBP7 promotes the proliferation and differentiation of primary myoblasts and intramuscular preadipocytes in chicken

Though it is well known that insulin-like growth factor (IGF) binding protein 7 (IGFBP7) plays an important role in myogenesis and adipogenesis in mammals, its impact on the proliferation, differentiation, and lipid deposition in chicken primary myoblasts (CPM) and intramuscular preadipocytes remains unexplored. In the present study, we firstly examined the correlation between SNPs within the genomic sequence of the IGFBP7 gene and carcass and blood chemical traits in a F2 resource population by genetic association analysis, and found that a significant correlation between the SNP (4_49499525) located in the intron region of IGFBP7 and serum high-density lipoproteins (HDL). We then examined the expression patterns of IGFBP7 across different stages of proliferation and differentiation in CPMs and intramuscular preadipocytes via qPCR, and explored the biological functions of IGFBP7 through gain- and loss-of-function experiments and a range of techniques including qPCR, CCK-8, EdU, flow cytometry, Western blot, immunofluorescence, and Oil Red O staining to detect the proliferation, differentiation, and lipid deposition in CPMs and intramuscular preadipocytes. We ascertained that the expression levels of the IGFBP7 gene increased as cell differentiation progresses in CPMs and intramuscular preadipocytes, and that IGFBP7 promotes the proliferation and differentiation of these cells, as well as facilitates intracellular lipid deposition. Furthermore, we investigated the regulatory mechanism of IGFBP7 expression by using co-transfection strategy and dual-luciferase reporter assay, and discovered that the myogenic transcription factors (MRF), myoblast determination factor (MyoD) and myogenin (MyoG), along with the adipocyte-specific transcription factor (TF) CCAAT/enhancer-binding protein α (C/EBPα), can bind to the core transcription activation region of the IGFBP7 promoter located 500 bp upstream from the transcription start site, thereby promoting IGFBP7 transcription and expression. Taken together, our study underscores the role of IGFBP7 as a positive regulator for myogenesis and adipogenesis, while also elucidating the functional and transcriptional regulatory mechanisms of IGFBP7 in chicken skeletal muscle development and intramuscular adipogenesis.

Profiling and functional analysis of circular RNAs in yaks intramuscular fat

Circular RNAs (circRNAs) are a new class of endogenous RNA regulating gene expression. However, the regulatory mechanisms of lipid metabolism in yaks involved in circRNAs remain poorly understood. The IMF plays a crucial role in the quality of yak meat, to greatly improve the meat quality. In this study, the fatty acid profiles of yak IMF were determined and circRNAs were sequenced. The results showed that the total of polyunsaturated fatty acid (PUFA) content of adult yak muscle was significantly higher than that in yak calves (p < 0.05). A total of 29,021 circRNAs were identified in IMF tissue, notably, 99 differentially expressed (DE) circRNAs were identified, to be associated with fat deposition, the most significant of which were circ_12686, circ_6918, circ_3582, ci_106 and ci_123 (A circRNA composed of exons is labelled 'circRNA' and a circRNA composed of introns is labelled 'ciRNA'). KEGG pathway enrichment analysis showed that the differential circRNAs were enriched in four pathways associated with fat deposition (e.g., the peroxisome proliferator-activated receptor signalling, fatty acid degradation, sphingolipid metabolism and sphingolipid signalling pathways). We also constructed co-expression networks of DE circRNA-miRNA using high-throughput sequencing in IMF deposition, from which revealed that ci_106 target binding of bta-miR-130b, bta-miR-148a, bta-miR-15a, bta-miR-34a, bta-miR-130a, bta-miR-17-5p and ci_123 target binding of bta-miR-150 were involved in adipogenesis. The study revealed the role of the circRNAs in the IMF deposition in yak and its influence on meat quality the findings demonstrated the circRNA differences in the development of IMF with the increase of age, thus providing a theoretical basis for further research on the molecular mechanism of IMF deposition in yaks.

Identification and characterization of structural variants related to meat quality in pigs using chromosome-level genome assemblies

BACKGROUND

Many studies have been performed to identify various genomic loci and genes associated with the meat quality in pigs. However, the full genetic architecture of the trait still remains unclear in part because of the lack of accurate identification of related structural variations (SVs) which resulted from the shortage of target breeds, the limitations of sequencing data, and the incompleteness of genome assemblies. The recent generation of a new pig breed with superior meat quality, called Nanchukmacdon, and its chromosome-level genome assembly (the NCMD assembly) has provided new opportunities.

RESULTS

By applying assembly-based SV calling approaches to various genome assemblies of pigs including Nanchukmacdon, the impact of SVs on meat quality was investigated. Especially, by checking the commonality of SVs with other pig breeds, a total of 13,819 Nanchukmacdon-specific SVs (NSVs) were identified, which have a potential effect on the unique meat quality of Nanchukmacdon. The regulatory potentials of NSVs for the expression of nearby genes were further examined using transcriptome- and epigenome-based analyses in different tissues.

CONCLUSIONS

Whole-genome comparisons based on chromosome-level genome assemblies have led to the discovery of SVs affecting meat quality in pigs, and their regulatory potentials were analyzed. The identified NSVs will provide new insights regarding genetic architectures underlying the meat quality in pigs. Finally, this study confirms the utility of chromosome-level genome assemblies and multi-omics analysis to enhance the understanding of unique phenotypes.

Characteristics of Transcriptome and Metabolome Concerning Intramuscular Fat Content in Beijing Black Pigs

To study the characteristics of genes and metabolites related to intramuscular fat (IMF) content with less influence by breed background and individual differences, the skeletal muscle samples from 40 Beijing black pigs with either high or low IMF content were used to perform transcriptome and metabolome analyses. About 99 genes (twofold-change) were differentially expressed. Up-regulated genes in the high IMF pigs were mainly related to fat metabolism. The key genes in charge of IMF deposition are ADIPOQ, CIDEC, CYP4B1, DGAT2, LEP, OPRL1, PLIN1, SCD, and THRSP. KLHL40, TRAFD1, and HSPA6 were novel candidate genes for the IMF trait due to their high abundances. In the low IMF pigs, the differentially expressed genes involved in virus resistance were up-regulated. About 16 and 18 differential metabolites (1.5 fold-change) were obtained in the positive and negative modes, respectively. Pigs with low IMF had weaker fatty acid oxidation due to the down-regulation of various carnitines. Differentially expressed genes were more important in determining IMF deposition than differential metabolites because relatively few differential metabolites were obtained, and they were merely the products under the physiological status of diverged IMF content. This study provided valuable information for further studies on IMF deposition.

A Comprehensive Genomic Analysis of Chinese Indigenous Ningxiang Pigs: Genomic Breed Compositions, Runs of Homozygosity, and Beyond

Ningxiang pigs are a renowned indigenous pig breed in China, known for their meat quality, disease resistance, and environmental adaptability. In recent decades, consumer demand for meats from indigenous breeds has grown significantly, fueling the selection and crossbreeding of Ningxiang pigs (NXP). The latter has raised concerns about the conservation and sustainable use of Ningxiang pigs as an important genetic resource. To address these concerns, we conducted a comprehensive genomic study using 2242 geographically identified Ningxiang pigs. The estimated genomic breed composition (GBC) suggested 2077 pigs as purebred Ningxiang pigs based on a ≥94% NXP-GBC cut-off. The remaining 165 pigs were claimed to be crosses, including those between Duroc and Ningxiang pigs and between Ningxiang and Shaziling pigs, and non-Ningxiang pigs. Runs of homozygosity (ROH) were identified in the 2077 purebred Ningxiang pigs. The number and length of ROH varied between individuals, with an average of 32.14 ROH per animal and an average total length of 202.4 Mb per animal. Short ROH (1-5 Mb) was the most abundant, representing 66.5% of all ROH and 32.6% of total ROH coverage. The genomic inbreeding estimate was low (0.089) in purebred Ningxiang pigs compared to imported western pig breeds. Nine ROH islands were identified, pinpointing candidate genes and QTLs associated with economic traits of interest, such as reproduction, carcass and growth traits, lipid metabolism, and fat deposition. Further investigation of these ROH islands and candidate genes is anticipated to better understand the genomics of Ningxiang pigs.

Differences in Muscle Lipogenic Gene Expression, Carcass Traits and Fat Deposition among Three Iberian Pig Strains Finished in Two Different Feeding Systems

The Iberian pig breed includes several well-differentiated strains. The present study evaluated carcass traits, fat deposition and muscle expression of important lipogenic genes (SCD, ME1, ACACA, FASN, EGR1, ACOX and ACLY) using 65 male pigs of 3 Iberian strains (20 Lampiño, 23 Torbiscal, and 22 Retinto) finished either in a conventional, concentrate-based system (CF) or in montanera (MF), a traditional free-range system with acorn feeding. Torbiscal had the highest ham, Longissimus thoracis and prime cuts yields, and the thinnest subcutaneous adipose tissue (SAT). Retinto had the highest monounsaturated fatty acids (MUFA) and percentage of intramuscular fat (IMF), while Lampiño had the greatest content of saturated fatty acids (SFA), polyunsaturated fatty acids (PUFA), atherogenic (AI) and thrombogenic (TI) indexes in SAT. Conventionally finished pigs had the highest ham, L. thoracis and prime cuts yields, and SFA. Montanera-finished animals had the highest PUFA and MUFA contents, and the lowest AI, TI and n6/n3 ratio in SAT. In relation to gene expression, Retinto had the greatest SCD, FASN and ACLY levels. Most studied genes were overexpressed in CF pigs. In conclusion, MF pigs had healthier fat than CF pigs, and Retinto had the healthiest fat and the greatest lipogenic trend in muscle, supported by IMF and lipogenic gene expression.

Hepatic Transcriptome Analysis Reveals Genes, Polymorphisms, and Molecules Related to Lamb Tenderness

Tenderness is a key meat quality trait that determines the public acceptance of lamb consumption, so genetic improvement toward lamb with higher tenderness is pivotal for a sustainable sheep industry. However, unravelling the genomics controlling the tenderness is the first step. Therefore, this study aimed to identify the transcriptome signatures and polymorphisms related to divergent lamb tenderness using RNA deep sequencing. Since the molecules and enzymes that control muscle growth and tenderness are metabolized and synthesized in the liver, hepatic tissues of ten sheep with divergent phenotypes: five high- and five low-lamb tenderness samples were applied for deep sequencing. Sequence analysis identified the number of reads ranged from 21.37 to 25.37 million bases with a mean value of 22.90 million bases. In total, 328 genes are detected as differentially expressed (DEGs) including 110 and 218 genes that were up- and down-regulated, respectively. Pathway analysis showed steroid hormone biosynthesis as the dominant pathway behind the lamb tenderness. Gene expression analysis identified the top high (such as TP53INP1, CYP2E1, HSD17B13, ADH1C, and LPIN1) and low (such as ANGPTL2, IGFBP7, FABP5, OLFML3, and THOC5) expressed candidate genes. Polymorphism and association analysis revealed that mutation in OLFML3, ANGPTL2, and THOC5 genes could be potential candidate markers for tenderness in sheep. The genes and pathways identified in this study cause variation in tenderness, thus could be potential genetic markers to improve meat quality in sheep. However, further validation is needed to confirm the effect of these markers in different sheep populations so that these could be used in a selection program for lamb with high tenderness.

Comparative and Functional Analysis of miRNAs and mRNAs Involved in Muscle Fiber Hypertrophy of Juvenile and Adult Goats

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate several pathway intermediates and affect the skeletal muscle development in mice, pigs, sheep, and cattle. However, to date, only a small number of miRNAs have been reported in the muscle development of goats. In this report, the longissimus dorsi transcripts of one- and ten-month-old goats were analyzed by sequencing RNAs and miRNAs. The results showed that the ten-month-old Longlin goats had 327 up- and 419 down-regulated differentially expressed genes (DEGs) compared with the one-month-old. In addition, 20 co-up-regulated and 55 co-down-regulated miRNAs involved in the muscle fiber hypertrophy of goats were identified in ten-month-old Longlin and Nubian goats compared with one-month-old. Five miRNA-mRNA pairs (chi-let-7b-3p-MIRLET7A, chi-miR193b-3p-MMP14, chi-miR-355-5p-DGAT2, novel_128-LOC102178119, novel_140-SOD3) involved in the goat skeletal muscle development were identified by miRNA-mRNA negative correlation network analysis. Our results provided new insight into the functional roles of goat muscle-associated miRNAs, allowing a deeper understanding of the transformation of miRNA roles during mammalian muscle development.

Global analysis of the association between pig muscle fatty acid composition and gene expression using RNA-Seq

Fatty acids (FAs) play an essential role as mediators of cell signaling and signal transduction, affecting metabolic homeostasis and determining meat quality in pigs. However, FAs are transformed by the action of several genes, such as those encoding desaturases and elongases of FAs in lipogenic tissues. The aim of the current work was to identify candidate genes, biological processes, and pathways involved in the modulation of intramuscular FA profile from longissimus dorsi muscle. FA profile by gas chromatography of methyl esters and gene expression by RNA-Seq were determined in 129 Iberian × Duroc backcrossed pigs. An association analysis between the muscle transcriptome and its FA profile was performed, followed by a concordance and functional analysis. Overall, a list of well-known (e.g., PLIN1, LEP, ELOVL6, SC5D, NCOA2, ACSL1, MDH1, LPL, LGALS12, TFRC, GOT1, and FBP1) and novel (e.g., TRARG1, TANK, ENSSSCG00000011196, and ENSSSCG00000038429) candidate genes was identified, either in association with specific or several FA traits. Likewise, several of these genes belong to biological processes and pathways linked to energy, lipid, and carbohydrate metabolism, which seem determinants in the modulation of FA compositions. This study can contribute to elucidate the complex relationship between gene expression and FA profile in pig muscle.

Characterization of fatty acid compositions in longissimus thoracis muscle and identification of candidate gene and SNPs related to polyunsaturated fatty acid in Hu sheep

Fatty acid (FA) composition contributes greatly to the quality and nutritional value of lamb meat. In the present study, FA was measured in longissimus thoracis (LT) muscles of 1,085 Hu sheep using gas chromatography. Comparative transcriptomic analysis was conducted in LT muscles to identify differentially expressed genes (DEGs) between six individuals with high polyunsaturated fatty acids (H-PUFA, 15.27% ± 0.42%) and six with low PUFA (L-PUFA, 5.22% ± 0.25%). Subsequently, the single nucleotide polymorphisms (SNPs) in a candidate gene PLIN2 were correlated with FA traits. The results showed a total of 29 FA compositions and 8 FA groups were identified, with the highest content of monounsaturated fatty acids (MUFA, 46.54%, mainly C18:1n9c), followed by saturated fatty acids (SFA, 44.32%, mainly C16:0), and PUFA (8.72%, mainly C18:2n6c), and significant correlations were observed among the most of FA traits. Transcriptomic analyses identified 110 upregulated and 302 downregulated DEGs between H-PUFA and L-PUFA groups. The functional enrichment analysis revealed three significant Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways and 17 gene ontology (GO) terms, in which regulation of lipolysis in adipocytes, the AMPK signaling pathway, and the PPAR signaling pathway may play important roles in FA metabolism and biosynthesis. In addition, weighted gene co-expression network analysis (WGCNA) identified 37 module genes associated with PUFA-related traits. In general, PLIN1, LIPE, FABP4, LEP, ACACA, ADIPOQ, SCD, PCK2, FASN, PLIN2, LPL, FABP3, THRSP, and ACADVL may have a great impact on PUFA metabolism and lipid deposition. Four SNPs within PLIN2 were significantly associated with FA. Of those, SNP1 (g.287 G>A) was significantly associated with C18:1n9c and MUFA, and SNP4 (g.7807 T>C) was significantly correlated with PUFA (C18:3n3). In addition, the combined genotype of SNP1 (g.287 G>A), SNP3 (g.7664 T>C), and SNP4 (g.7807 T>C) were significantly correlated with C16:1, C17:0, C18:1C6, PUFA (C18:3n3, C22:6n3), and n-6/n-3 PUFA. These results contribute to the knowledge of the biological mechanisms and genetic markers involved in the composition of FA in Hu sheep.

VMP1 Regulated by chi-miR-124a Effects Goat Myoblast Proliferation, Autophagy, and Apoptosis through the PI3K/ULK1/mTOR Signaling Pathway

The production of goat meat is determined by the growth speed of muscle fibers, and the autophagy and apoptosis of myoblast cells is a crucial process in the growth of muscle fibers. The rapid growth of muscle fibers occurs from one month old to nine months old in goats; however, the mechanisms of myoblast cells’ autophagy and apoptosis in this process are still unknown. To identify candidate genes and signaling pathway mechanisms involved in myoblast apoptosis and autophagy, we compared the expression characteristics of longissimus dorsi tissues from Wu’an goats—a native goat breed of China—at 1 month old (mon1 group) and 9 months old (mon9 group). Herein, a total of 182 differentially expressed mRNAs (DEGs) in the mon1 vs. mon9 comparison, along with the KEGG enrichments, showed that the PI3K-Akt pathway associated with autophagy and apoptosis was significantly enriched. Among these DEGs, expression of vacuole membrane protein 1 (VMP1)—a key gene for the PI3K-Akt pathway—was significantly upregulated in the older goats relative to the 1-month-old goats. We demonstrated that VMP1 promotes the proliferation and autophagy of myoblasts, and inhibits their apoptosis. The integration analysis of miRNA–mRNA showed that miR-124a was a regulator of VMP1 in muscle tissue, and overexpression and inhibition of miR-124a suppressed the proliferation and autophagy of myoblasts. The PI3K/Akt/mTOR pathway was an important pathway for cell autophagy. Additionally, the activator of the PI3K/Akt/mTOR pathway, the expression of VMP1, and ULK1 were higher than the negative control, and the expression of mTOR was depressed. The expression of VMP1, ULK1, and mTOR was the opposite when the inhibitor was added to the myoblasts. These results show that the PI3K/Akt/mTOR pathway promoted the expression of VMP1 and ULK1. By using adenovirus-mediated apoptosis and proliferation assays, we found that that miR-124a inhibits myoblast proliferation and autophagy, and promotes their apoptosis by targeting VMP1. In conclusion, our results indicated that VMP1 was highly expressed in the LD muscle tissues of nine-month-old goats, and that it was regulated by miR-124a to inhibit myoblast cells’ apoptosis through the PI3K/Akt/mTOR pathway, and to promote proliferation and autophagy. These findings contribute to the understanding of the molecular mechanisms involved in myoblast proliferation, autophagy, and apoptosis.

Modeling microRNA-driven post-transcriptional regulation using exon-intron split analysis in pigs

The contribution of microRNAs (miRNAs) to mRNA post-transcriptional regulation has often been explored by the post hoc selection of downregulated genes and determining whether they harbor binding sites for miRNAs of interest. This approach, however, does not discriminate whether these mRNAs are also downregulated at the transcriptional level. Here, we have characterized the transcriptional and post-transcriptional changes in mRNA expression in two porcine tissues: gluteus medius muscle of fasted and fed Duroc gilts and adipose tissue of lean and obese Duroc-Göttingen minipigs. Exon-intron split analysis of RNA-seq data allowed us to identify downregulated mRNAs with high post-transcriptional signals in fed or obese states, and we assessed whether they harbor binding sites for upregulated miRNAs in any of these two physiological states. We found 26 downregulated mRNAs with high post-transcriptional signals in the muscle of fed gilts and 21 of these were predicted targets of miRNAs upregulated in fed pigs. For adipose tissue, 44 downregulated mRNAs in obese minipigs displayed high post-transcriptional signals, and 25 of these were predicted targets of miRNAs upregulated in the obese state. These results suggest that the contribution of miRNAs to mRNA repression is more prominent in the skeletal muscle system. Finally, we identified several genes that may play relevant roles in the energy homeostasis of the pig skeletal muscle (DKK2 and PDK4) and adipose (SESN3 and ESRRG) tissues. By differentiating transcriptional from post-transcriptional changes in mRNA expression, exon-intron split analysis provides a valuable view of the regulation of gene expression, complementary to canonical differential expression analyses.

Impact of dietary Lactobacillus supplementation on intramuscular fat deposition and meat quality of Sunit sheep

To investigate the impacts of dietary Lactobacillus supplementation on meat quality such as edible quality and nutritional value of Sunit sheep, a 90-day feeding experiment (Lactobacillus dietary group, R group; non-Lactobacillus dietary group, C group) using twelve 3-month-old Sunit sheep was conducted. The deposition of intramuscular fat (IMF) was increased (p < .05) while the share force and cooking loss were decreased (p < .05) in the R group compared with the C group. The proportions of seven kinds of fatty acids (FAs) have changed significantly (p < .05), especially with higher functional FAs and lower trans-FA in the R group. Metabonomics analysis showed that the metabolites and pathway-related lipid syntheses, such as carnitine cycle, tricarboxylic acid cycle, and glycerophosphocholine metabolic pathway, have significantly changed in the R group. The Lactobacillus dietary supplements impacted the variation of IMF deposition and FAs composition by altering the lipid metabolism pathways of Sunit sheep and then changed the edible quality and nutritional value. PRACTICAL APPLICATIONS: It is well known that the intramuscular fat (IMF) and fatty acids composition in livestock is positively correlated with various aspects of meat quality such as edible quality and nutritional value, which are related to consumer preference. The present study analyzed the effects of Lactobacillus supplement on the intramuscular fat deposition and meat quality of Sunit sheep, which resulted in the increase of IMF, and the differences of fatty acids composition, especially the functional fatty acids. It was explored the mechanism of Lactobacillus affect the variation of lipid metabolism pathways and key metabolites in sheep, which suggested that altering the feeding regimen could improve the meat quality of agri-animals.

Understanding the role of myoglobin content in Iberian pigs fattened in an extensive system through analysis of the transcriptome profile

Meat color is the first perceived sensory feature and one of the most important quality traits. Myoglobin is the main pigment in meat, giving meat its characteristic cherry‐red color, highly appreciated by the consumers. In the current study, we used the RNA‐seq technique to characterize the longissimus dorsi muscle transcriptome in two groups of Iberian pigs with divergent breeding values for myoglobin content. As a result, we identified 57 differentially expressed genes and transcripts (DEGs). Moreover, we have validated the RNA‐seq expression of a set of genes by quantitative PCR (qPCR). Functional analyses revealed an enrichment of DEGs in biological processes related to oxidation (HBA1), lipid metabolism (ECH1, PLA2G10, PLD2), inflammation (CHST1, CD209, PLA2G10), and immune system (CD209, MX2, LGALS3, LGALS9). The upstream analysis showed a total of five transcriptional regulatory factors and eight master regulators that could moderate the expression of some DEGs, highlighting SPI1 and MAPK1, since they regulate the expression of DEGs involved in immune defense and inflammatory processes. Iberian pigs with high myoglobin content also showed higher expression of the HBA1 gene and both molecules, myoglobin and hemoglobin, have been described as having a protective effect against oxidative and inflammatory processes. Therefore, the HBA1 gene is a very promising candidate gene to harbor polymorphisms underlying myoglobin content, whereby further studies should be carried out for its potential use in an Iberian pig selection program.

Changes in Biceps femoris Transcriptome along Growth in Iberian Pigs Fed Different Energy Sources and Comparative Analysis with Duroc Breed

Simple Summary

The genetic mechanisms that regulate biological processes, such as skeletal muscle development and growth, or intramuscular fat deposition, have attracted great interest, given their impact on production traits and meat quality. In this sense, a comparison of the transcriptome of skeletal muscle between phenotypically different pig breeds, or along growth, could be useful to improve the understanding of the molecular processes underlying the differences in muscle metabolism and phenotypic traits, potentially driving the identification of causal genes, regulators and metabolic pathways involved in their variability.

Abstract

This experiment was conducted to investigate the effects of developmental stage, breed, and diet energy source on the genome-wide expression, meat quality traits, and tissue composition of biceps femoris muscle in growing pure Iberian and Duroc pigs. The study comprised 59 Iberian (IB) and 19 Duroc (DU) animals, who started the treatment at an average live weight (LW) of 19.9 kg. The animals were kept under identical management conditions and fed two diets with different energy sources (6% high oleic sunflower oil or carbohydrates). Twenty-nine IB animals were slaughtered after seven days of treatment at an average LW of 24.1 kg, and 30 IB animals plus all the DU animals were slaughtered after 47 days at an average LW of 50.7 kg. The main factors affecting the muscle transcriptome were age, with 1832 differentially expressed genes (DEGs), and breed (1055 DEGs), while the effect of diet on the transcriptome was very small. The results indicated transcriptome changes along time in Iberian animals, being especially related to growth and tissue development, extracellular matrix (ECM) composition, and cytoskeleton organization, with DEGs affecting relevant functions and biological pathways, such as myogenesis. The breed also affected functions related to muscle development and cytoskeleton organization, as well as functions related to solute transport and lipid and carbohydrate metabolism. Taking into account the results of the two main comparisons (age and breed effects), we can postulate that the Iberian breed is more precocious than the Duroc breed, regarding myogenesis and muscle development, in the studied growing stage.

Genome-Wide Phylogenetic Analysis, Expression Pattern, and Transcriptional Regulatory Network of the Pig C/EBP Gene Family

The CCAAT/enhancer binding protein (C/EBP) transcription factors (TFs) regulate many important biological processes, such as energy metabolism, inflammation, cell proliferation etc. A genome-wide gene identification revealed the presence of a total of 99 C/EBP genes in pig and 19 eukaryote genomes. Phylogenetic analysis showed that all C/EBP TFs were classified into 6 subgroups named C/EBPα, C/EBPβ, C/EBPδ, C/EBPε, C/EBPγ, and C/EBPζ. Gene expression analysis showed that the C/EBPα, C/EBPβ, C/EBPδ, C/EBPγ, and C/EBPζ genes were expressed ubiquitously with inconsistent expression patterns in various pig tissues. Moreover, a pig C/EBP regulatory network was constructed, including C/EBP genes, TFs and miRNAs. A total of 27 feed-forward loop (FFL) motifs were detected in the pig C/EBP regulatory network. Based on the RNA-seq data, gene expression patterns related to FFL sub-network were analyzed in 27 adult pig tissues. Certain FFL motifs may be tissue specific. Functional enrichment analysis indicated that C/EBP and its target genes are involved in many important biological pathways. These results provide valuable information that clarifies the evolutionary relationships of the C/EBP family and contributes to the understanding of the biological function of C/EBP genes.

Influence of genetic background and dietary oleic acid on gut microbiota composition in Duroc and Iberian pigs

BACKGROUND

Phenotypic variability for productive and meat quality traits has been largely studied in Iberian pigs, especially in genetic selection and nutritional experiments. Complex interactions among genetic background, diet composition and gut microbiota hinder the correct assessment of each factor's contribution on phenotypes. In order to disentangle these interactions, we evaluated changes in gut microbiota composition comparing 48 Iberian and Duroc pigs fed diets with different energy source (standard diet with carbohydrates vs sunflower oil-enriched diet with high oleic acid content).

RESULTS

A higher richness was observed for Iberian pigs (p < 0.05) and compositional analysis was applied for beta-diversity, differential abundance and pairwise log-ratio analyses. We found significant differences in overall microbiota composition between breeds, and also between diets inside breeds, to a lesser extent. Differential abundance analysis revealed that Duroc animals have more proportion of Actinobacteria and Prevotella, while Iberian replace those microorganisms with other more variable taxa. According to dietary differences, high-oleic fed animals were richer in Prevotella. We also found microbial ratios capable of separating animals by breeds and diets, mostly related to Actinobacteria.

CONCLUSION

This study reveals that both genetic background and diet composition might have a relevant impact in gut microbiota composition. The application of compositional data analysis has facilitated the identification of microorganisms and ratios as possibly related to metabolic changes due to genetic background and, to a lower extent, to dietary changes. This may lead to a relevant progress in the knowledge of interactions between pig genetics, environment and gut microbiota.

Muscle Transcriptome Analysis Reveals Molecular Pathways Related to Oxidative Phosphorylation, Antioxidant Defense, Fatness and Growth in Mangalitsa and Moravka Pigs

This work was aimed at evaluating loin transcriptome and metabolic pathway differences between the two main Serbian local pig breeds with divergent characteristics regarding muscle growth and fatness, as well as exploring nutrigenomic effects of tannin supplementation in Mangalitsa (MA) pigs. The study comprised 24 Mangalitsa and 10 Moravka (MO) males, which were kept under identical management conditions. Mangalitsa animals were divided in two nutritional groups (n = 12) receiving a standard (control) or tannin-supplemented diet (1.5%; MAT). Moravka pigs were fed the standard mixture. All animals were slaughtered at a similar age; 120 kg of average live weight (LW) and loin tissue was used for RNA-seq analysis. Results showed 306 differentially expressed genes (DEGs) according to breed, enriched in genes involved in growth, lipid metabolism, protein metabolism and muscle development, such as PDK4, FABP4, MYOD1 and STAT3, as well as a relevant number of genes involved in mitochondrial respiratory activity (MT-NDs, NDUFAs among others). Oxidative phosphorylation was the most significantly affected pathway, activated in Mangalitsa muscle, revealing the basis of a different muscle metabolism. Also, many other relevant pathways were affected by breed and involved in oxidative stress response, fat accumulation and development of skeletal muscle. Results also allowed the identification of potential regulators and causal networks such as those controlled by FLCN, PPARGC1A or PRKAB1 with relevant regulatory roles on DEGs involved in mitochondrial and lipid metabolism, or IL3 and TRAF2 potentially controlling DEGs involved in muscle development. The Tannin effect on transcriptome was small, with only 23 DEGs, but included interesting ones involved in lipid deposition such as PPARGC1B. The results indicate a significant effect of the breed on muscle tissue gene expression, affecting relevant biological pathways and allowing the identification of strong regulatory candidate genes to underlie the gene expression and phenotypic differences between the compared groups.

Key Genes Regulating Skeletal Muscle Development and Growth in Farm Animals

Simple Summary

Skeletal muscle mass is an important economic trait, and muscle development and growth is a crucial factor to supply enough meat for human consumption. Thus, understanding (candidate) genes regulating skeletal muscle development is crucial for understanding molecular genetic regulation of muscle growth and can be benefit the meat industry toward the goal of increasing meat yields. During the past years, significant progress has been made for understanding these mechanisms, and thus, we decided to write a comprehensive review covering regulators and (candidate) genes crucial for muscle development and growth in farm animals. Detection of these genes and factors increases our understanding of muscle growth and development and is a great help for breeders to satisfy demands for meat production on a global scale.

Abstract

Farm-animal species play crucial roles in satisfying demands for meat on a global scale, and they are genetically being developed to enhance the efficiency of meat production. In particular, one of the important breeders’ aims is to increase skeletal muscle growth in farm animals. The enhancement of muscle development and growth is crucial to meet consumers’ demands regarding meat quality. Fetal skeletal muscle development involves myogenesis (with myoblast proliferation, differentiation, and fusion), fibrogenesis, and adipogenesis. Typically, myogenesis is regulated by a convoluted network of intrinsic and extrinsic factors monitored by myogenic regulatory factor genes in two or three phases, as well as genes that code for kinases. Marker-assisted selection relies on candidate genes related positively or negatively to muscle development and can be a strong supplement to classical selection strategies in farm animals. This comprehensive review covers important (candidate) genes that regulate muscle development and growth in farm animals (cattle, sheep, chicken, and pig). The identification of these genes is an important step toward the goal of increasing meat yields and improves meat quality.

Comparative Transcriptome Profile Analysis of Longissimus dorsi Muscle Tissues From Two Goat Breeds With Different Meat Production Performance Using RNA-Seq

Carcass weight, meat quality and muscle components are important traits economically and they underpin most of the commercial return to goat producers. In this study, the Longissimus dorsi muscle tissues were collected from five Liaoning cashmere (LC) goats and five Ziwuling black (ZB) goats with phenotypic difference in carcass weight, some meat quality traits and muscle components. The histological quantitative of collagen fibers and the transcriptome profiles in the Longissimus dorsi muscle tissues were investigated using Masson-trichrome staining and RNA-Seq, respectively. The percentage of total collagen fibers in the Longissimus dorsi muscle tissues from ZB goats was less than those from LC goats, suggesting that these ZB goats had more tender meat. An average of 15,919 and 15,582 genes were found to be expressed in Longissimus dorsi muscle tissues from LC and ZB goats, respectively. Compared to LC goats, the expression levels of 78 genes were up-regulated in ZB goats, while 133 genes were down-regulated. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses revealed that the differentially expressed genes (DEGs) were significantly enriched in GO terms related to the muscle growth and development and the deposition of intramuscular fat and lipid metabolism, hippo signaling pathway and Jak-STAT signaling pathway. The results provide an improved understanding of the genetic mechanisms regulating meat production performance in goats, and will help us improve the accuracy of selection for meat traits in goats using marker-assisted selection based on these differentially expressed genes obtained.

Comparative Transcriptome Analyses of Longissimus thoracis Between Pig Breeds Differing in Muscle Characteristics

The two breeds, Mashen (MS; a northern China breed) and Large White (LW; a western lean breed) pigs, show important phenotypic differences in growth, adaptability, intramuscular fat (IMF) content, and energy metabolism since early developmental stage. The main aim of this study was the evaluation of longissimus thoracis muscle transcriptome profile of both genetic types to identify genes, pathways responsible for their differentiated phenotype. Longissimus thoracis of MS and LW pigs were sampled at 0 day (early stage), 90 days (middle stage), and 180 days (later stage) after birth. A total of 3,487 differentially expressed genes (DEGs) were identified at the three time points. Sample clustering analysis revealed the slower growth rate of MS than LW pigs. Gene expression pattern analysis revealed that multicellular organism growth genes (GHSR, EZR, FOXS1, DRD2, SH3PXD2B, CSF1, and TSHR) were involved in the fast growth rate of LW pigs at early stage. Furthermore, DEGs (ACACA, ACSF3, OXSM, CBR4, and HSD17B8) functionally related to fatty acid synthesis revealed that IMF accumulation occurred around 90 and up to 180 days. These DEGs provided valuable resource to study the phenotypic difference in longissimus thoracis muscle between MS and LW pigs.

Major quantitative trait loci influencing milk production and conformation traits in Guernsey dairy cattle detected on Bos taurus autosome 19

The goal of this study was to identify potential quantitative trait loci (QTL) for 27 production, fitness, and conformation traits of Guernsey cattle through genome-wide association (GWA) analyses, with extra emphasis on BTA19, where major QTL were observed for several traits. Animals' de-regressed predicted transmitting abilities (PTA) from the December 2018 traditional US evaluation were used as phenotypes. All of the Guernsey cattle included in the QTL analyses were predictor animals in the reference population, ranging from 1,077 to 1,685 animals for different traits. Single-trait GWA analyses were carried out by a mixed-model approach for all 27 traits using imputed high-density genotypes. A major QTL was detected on BTA19, influencing several milk production traits, conformation traits, and livability of Guernsey cattle, and the most significant SNP lie in the region of 26.2 to 28.3 Mb. The myosin heavy chain 10 (MYH10) gene residing within this region was found to be highly associated with milk production and body conformation traits of dairy cattle. After the initial GWA analyses, which suggested that many significant SNP are in linkage with one another, conditional analyses were used for fine mapping. The top significant SNP on BTA19 were fixed as covariables in the model, one at a time, until no more significant SNP were detected on BTA19. After this fine-mapping approach was applied, only 1 significant SNP was detected on BTA19 for most traits, but multiple, independent significant SNP were found for protein yield, dairy form, and stature. In addition, the haplotype that hosts the major QTL on BTA19 was traced to a US Guernsey born in 1954. The haplotype is common in the breed, indicating a long-term influence of this QTL on the US Guernsey population.

Differences in the Loin Tenderness of Iberian Pigs Explained through Dissimilarities in Their Transcriptome Expression Profile

Simple Summary

The Iberian pig is the most representative autochthonous breed of the Mediterranean region with unique genetic and phenotypic characteristics. The breed has been successfully preserved by its high-quality meat and high-priced products. Tenderness is one of the most relevant meat quality traits, and meat tenderization is influenced by genetic and environmental effects such as pre-slaughter handling and post-mortem conditions. Tenderness could be included in Iberian pig breeding programs, mainly focused on the improvement of premium-cuts percentage, in order to avoid the meat quality decline. A better biological understanding of this trait is needed. In the current study, we analyze the transcriptome of pigs divergent for Warner–Bratzler shear force through RNA-seq technique for the identification, characterization and quantification of candidate genes involved in biological pathways, networks and functions affecting meat tenderness.

Abstract

Tenderness is one of the most important meat quality traits and it can be measured through shear force with the Warner–Bratzler test. In the current study, we use the RNA-seq technique to analyze the transcriptome of Longissimus dorsi (LD) muscle in two groups of Iberian pigs (Tough and Tender) divergent for shear force breeding values. We identified 200 annotated differentially expressed genes (DEGs) and 245 newly predicted isoforms. The RNAseq expression results of 10 genes were validated with quantitative PCR (qPCR). Functional analyses showed an enrichment of DE genes in biological processes related to proteolysis (CTSC, RHOD, MYH8, ACTC1, GADD45B, CASQ2, CHRNA9 and ANKRD1), skeletal muscle tissue development (ANKRD1, DMD, FOS and MSTN), lipid metabolism (FABP3 and PPARGC1A) and collagen metabolism (COL14A1). The upstream analysis revealed a total of 11 transcription regulatory factors that could regulate the expression of some DEGs. Among them, IGF1, VGLL3 and PPARG can be highlighted since they regulate the expression of genes involved in biological pathways that could affect tenderness. The experiment revealed a set of candidate genes and regulatory factors suggestive to search polymorphisms that could be incorporated in a breeding program for improving meat tenderness.

Transcriptome profiling of longissimus lumborum in Holstein bulls and steers with different beef qualities

Previous research regarding Holstein cows has mainly focused on increasing milk yield. However, in order to maximize the economical profits of Holstein cattle farming, it is necessary to fully take advantage of Holstein bulls to produce high-grade beef. The present study aims to investigate different transcriptomic profiling of Holstein bulls and steers, via high-throughput RNA-sequencing (RNA-seq). The growth and beef quality traits of Holstein steers and bulls were characterized via assessment of weight, rib eye area, marbling score, shear force and intramuscular fat percentage of the longissimus lumborum (LL) muscle. The results indicated that castration improved the meat quality, yet reduced the meat yield. Subsequently, RNA-seq of the LL muscle from Holstein steers and bulls revealed a total of 56 differentially expressed genes (DEGs). We performed the functional enrichment analysis in Gene Ontology (GO) annotations of the DEGs using GOseq R package software and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis using KOBAS tool. Through the integrated analysis of DEGs with reported QTLs and SNPs, seven promising candidate genes potentially affecting the beef quality of LL muscle following castration were discovered, including muscle structural protein coding genes (MYH1, MYH4, MYH10) and functional protein coding genes (GADL1, CYP2R1, EEPD1, SHISA3). Among them, MYH10, GADL1, CYP2R1, EEPD1 and SHISA3 were novel candidate genes associated with beef quality traits. Notably, EEPD1 was associated with both meat quality and reproduction traits, thus indicating its overlapping role in responding to hormone change, and subsequently inducing beef quality improvement. Our findings provide a complete dataset of gene expression profile of LL in Holstein bulls and steers, and will aid in understanding how castration influence meat yield and quality.

Differential Effects of Litter Size and Within-Litter Birthweight on Postnatal Traits of Fatty Pigs

Simple Summary

The current work evaluated the relative effects of the number of piglets in the litter of origin (high vs. low litter size) and the birthweight (low (LBW) vs. normal (NBW) in large litters) on postnatal performance and quality of carcasses and meat in purebred Iberian pigs. The results indicated that NBW piglets born in large litters had disparities in developmental patterns in weight and size, back-fat deposition, and fatty acid composition of viscerae compared to NBW piglets from small litters, which again were different from those found in LBW piglets when compared to their NBW counterparts. However, both growth patterns were altered and might indicate previous phases of metabolic disorders.

Abstract

Fatty pigs are characterized by a thrifty genotype, adapted to harsh environments based on changes in metabolism and energy saving. Thus, we hypothesized that feto-maternal energy partitioning in large litters might have postnatal effects that might be independent of intrauterine growth restriction (IUGR) processes. Hence, the current work reported the influence of two effects on postnatal performance and carcass and meat quality of purebred Iberian pigs: (a) the effects of the number of piglets in the litter (high vs. low litter size), and (b) the effects of birthweight (low (LBW) vs. normal (NBW)) in large litters. The results confirmed that NBW piglets born in large litters had differences in developmental patterns of weight, back-fat deposition, and fatty acid (FA) composition when compared to NBW piglets from small litters. These results were different from those found in LBW piglets when compared to their NBW counterparts, which showed an initial asymmetrical growth and altered muscle FA composition at slaughtering. The assessment of FA composition indicated better metabolic status in NBW piglets from large litters than in LBW piglets. These data support the concept that the prenatal environment, even when the individual may cope with it, inescapably affects postnatal life.

The whole-transcriptome landscape of muscle and adipose tissues reveals the ceRNA regulation network related to intramuscular fat deposition in yak

Background

The Intramuscular fat (IMF) content in meat products, which is positively correlated with meat quality, is an important trait considered by consumers. The regulation of IMF deposition is species specific. However, the IMF-deposition-related mRNA and non-coding RNA and their regulatory network in yak (Bos grunniens) remain unknown. High-throughput sequencing technology provides a powerful approach for analyzing the association between transcriptome-related differences and specific traits in animals. Thus, the whole transcriptomes of yak muscle and adipose tissues were screened and analyzed to elucidate the IMF deposition-related genes. The muscle tissues were used for IMF content measurements.

Results

Significant differences were observed between the 0.5- and 2.5-year-old yaks. Several mRNAs, miRNAs, lncRNAs and circRNAs were generally expressed in both muscle and adipose tissues. Between the 0.5- and 2.5-year-old yaks, 149 mRNAs, 62 miRNAs, 4 lncRNAs, and 223 circRNAs were differentially expressed in muscle tissue, and 72 mRNAs, 15 miRNAs, 9 lncRNAs, and 211 circRNAs were differentially expressed in adipose tissue. KEGG annotation revelved that these differentially expressed genes were related to pathways that maintain normal biological functions of muscle and adipose tissues. Moreover, 16 mRNAs, 5 miRNAs, 3 lncRNAs, and 5 circRNAs were co-differentially expressed in both types of tissue. We suspected that these co-differentially expressed genes were involved in IMF-deposition in the yak. Additionally, LPL, ACADL, SCD, and FASN, which were previously shown to be associated with the IMF content, were identified in the competing endogenous RNA (ceRNA) regulatory network that was constructed on the basis of the IMF deposition-related genes. Three ceRNA subnetworks also revealed that TCONS-00016416 and its target SIRT1 “talk” to each other through the same miR-381-y and miR-208 response elements, whereas TCONS-00061798 and its target PRKCA, and TCONS-00084092 and its target LPL “talk” to each other through miR-122-x and miR-499-y response elements, respectively.

Conclusion

Taken together, our results reveal the potential mRNA and noncoding RNAs involved in IMF deposition in the yak, providing a useful resource for further research on IMF deposition in this animal species.

Comparative Transcriptomic Analysis of Subcutaneous Adipose Tissue from Local Pig Breeds

When compared to modern lean-type breeds, Portuguese local Alentejano (AL) and Bísaro (BI) pig breeds present a high potential for subcutaneous and intramuscular fat (IMF) deposition which contributes for better meat quality. The aim of this work was to explore the genome function to better understand the underlying physiological mechanisms associated with body fat accretion. Dorsal subcutaneous fat samples were collected at slaughter from adult animals (n = 4 for each breed) with ~150 kg body weight. Total RNA was obtained and sequenced for transcriptome analysis using DESeq2. A total of 458 differentially expressed (DE) genes (q-value < 0.05) were identified, with 263 overexpressed in AL and 195 in BI. Key genes involved in de novo fatty acid biosynthesis, elongation and desaturation were upregulated in AL such as ACLY, FASN, ME1, ELOVL6 and SCD. A functional enrichment analysis of the DE genes was performed using Ingenuity Pathway Analysis. Cholesterol synthesis is suggested to be higher in AL via SREBF2, SCAP and PPARG, while lipolytic activity may be more active in BI through GH and AMPK signalling. Increased signalling of CD40 together with the predicted activation of INSIG1 and INSIG2 in BI suggests that this breed is more sensitive to insulin whereas the AL is less sensitive like the Iberian breed.

Enhanced identification of significant regulators of gene expression

Background

Diseases like cancer will lead to changes in gene expression, and it is relevant to identify key regulatory genes that can be linked directly to these changes. This can be done by computing a Regulatory Impact Factor (RIF) score for relevant regulators. However, this computation is based on estimating correlated patterns of gene expression, often Pearson correlation, and an assumption about a set of specific regulators, normally transcription factors. This study explores alternative measures of correlation, using the Fisher and Sobolev metrics, and an extended set of regulators, including epigenetic regulators and long non-coding RNAs (lncRNAs). Data on prostate cancer have been used to explore the effect of these modifications.

Results

A tool for computation of RIF scores with alternative correlation measures and extended sets of regulators was developed and tested on gene expression data for prostate cancer. The study showed that the Fisher and Sobolev metrics lead to improved identification of well-documented regulators of gene expression in prostate cancer, and the sets of identified key regulators showed improved overlap with previously defined gene sets of relevance to cancer. The extended set of regulators lead to identification of several interesting candidates for further studies, including lncRNAs. Several key processes were identified as important, including spindle assembly and the epithelial-mesenchymal transition (EMT).

Conclusions

The study has shown that using alternative metrics of correlation can improve the performance of tools based on correlation of gene expression in genomic data. The Fisher and Sobolev metrics should be considered also in other correlation-based applications.

Trend analysis of the role of circular RNA in goat skeletal muscle development

Background

Circular RNA (circRNA) is produced during the splicing of mRNA (in addition to linear splicing) and is part of the gene regulatory network. The temporal expression patterns the different developmental stages were inseparable from these molecules’ function.

Results

Skeletal muscles of Anhui white goat (AWG) across seven fetal to postnatal development stages were sequenced and 21 RNA sequencing libraries were constructed. We thereby identified 9090 circRNAs and analyzed their molecular properties, temporal expression patterns, and potential functions at the different stages. CircRNAs showed complexities and diversity of formation as the same host gene produces multiple isoforms of these nucleic acids with different expression profiles. The differential expression of 2881 circRNAs (DECs, P < 0.05) was identified and four were randomly selected and validated by qPCR. Moreover, 1118 DECs under strict selected (SDECs, |log2^FC| > 2 and P-adj value < 0.01) showed 4 expression trends (Clusters 0, 19, 16 and 18). Cluster 0 molecules had increasing expression at all stages with effects on muscle through metabolism, regulation of enzyme activity, and biosynthesis. Cluster 16 circRNAs had high expression in the early and late stages and are involved in “Wnt signaling pathway”, “AMPK signaling pathway” and others. Cluster 18 molecules were mainly expressed at F120 and participate in “cytoskeletal protein binding”, “Notch signaling pathway” and so on. Cluster 19 circRNAs were down-regulated at all stages and related to muscle structure and development. Lastly, the SDECs divided the period of skeletal muscle development into three transitional stages: stage 1 (F45 to F90), which related to muscle satellite cell proliferation and muscle fiber structure; stage 2 (F90 to B1), in which the attachment of the cytoplasmic surface to the actin cytoskeleton initiates; and stage 3, which involved the “cGMP-PKG signaling pathway”. Moreover, the paraffin sections messages also validated that there are three transitional stages of skeletal muscle development.

Conclusion

Our current study provides a catalog of goat muscle-related circRNAs that can stratify skeletal muscle development fetus 45 days to newborn 90 days into three developmental stages. These findings better our understanding of functional transitions during mammalian muscle development.

Lipid Deposition and Metabolism in Local and Modern Pig Breeds: A Review

Modern pig breeds, which have been genetically improved to achieve fast growth and a lean meat deposition, differ from local pig breeds with respect to fat deposition, fat specific metabolic characteristics and various other properties. The present review aimed to elucidate the mechanisms underlying the differences between fatty local and modern lean pig breeds in adipose tissue deposition and lipid metabolism, taking into consideration morphological, cellular, biochemical, transcriptomic and proteomic perspectives. Compared to modern breeds, local pig breeds accumulate larger amounts of fat, which generally contains more monounsaturated and saturated fatty acids; they exhibit a higher adipocyte size and higher activity of lipogenic enzymes. Studies using transcriptomic and proteomic approaches highlighted several processes like immune response, fatty-acid turn-over, oxidoreductase activity, mitochondrial function, etc. which differ between local and modern pig breeds.

Impact of genotype, body weight and sex on the prenatal muscle transcriptome of Iberian pigs

Growth is dependent on genotype and diet, even at early developmental stages. In this study, we investigated the effects of genotype, sex, and body weight on the fetal muscle transcriptome of purebred Iberian and crossbred Iberian x Large White pigs sharing the same uterine environment. RNA sequencing was performed on 16 purebred and crossbred fetuses with high body weight (340±14g and 415±14g, respectively) and 16 with low body weight (246±14g and 311±14g, respectively), on gestational day 77. Genotype had the greatest effect on gene expression, with 645 genes identified as differentially expressed (DE) between purebred and crossbred animals. Functional analysis showed differential regulation of pathways involved in energy and lipid metabolism, muscle development, and tissue disorders. In purebred animals, fetal body weight was associated with 35 DE genes involved in development, lipid metabolism and adipogenesis. In crossbred animals, fetal body weight was associated with 60 DE genes involved in muscle development, viability, and immunity. Interestingly, the results suggested an interaction genotype*weight for some DE genes. Fetal sex had only a modest effect on gene expression. This study allowed the identification of genes, metabolic pathways, biological functions and regulators related to fetal genotype, weight and sex, in animals sharing the same uterine environment. Our findings contribute to a better understanding of the molecular events that influence prenatal muscle development and highlight the complex interactions affecting transcriptional regulation during development.

Co-expression network analysis predicts a key role of microRNAs in the adaptation of the porcine skeletal muscle to nutrient supply

BACKGROUND

The role of non-coding RNAs in the porcine muscle metabolism is poorly understood, with few studies investigating their expression patterns in response to nutrient supply. Therefore, we aimed to investigate the changes in microRNAs (miRNAs), long intergenic non-coding RNAs (lincRNAs) and mRNAs muscle expression before and after food intake.

RESULTS

We measured the miRNA, lincRNA and mRNA expression levels in the gluteus medius muscle of 12 gilts in a fasting condition (AL-T0) and 24 gilts fed ad libitum during either 5 h. (AL-T1, N = 12) or 7 h. (AL-T2, N = 12) prior to slaughter. The small RNA fraction was extracted from muscle samples retrieved from the 36 gilts and sequenced, whereas lincRNA and mRNA expression data were already available. In terms of mean and variance, the expression profiles of miRNAs and lincRNAs in the porcine muscle were quite different than those of mRNAs. Food intake induced the differential expression of 149 (AL-T0/AL-T1) and 435 (AL-T0/AL-T2) mRNAs, 6 (AL-T0/AL-T1) and 28 (AL-T0/AL-T2) miRNAs and none lincRNAs, while the number of differentially dispersed genes was much lower. Among the set of differentially expressed miRNAs, we identified ssc-miR-148a-3p, ssc-miR-22-3p and ssc-miR-1, which play key roles in the regulation of glucose and lipid metabolism. Besides, co-expression network analyses revealed several miRNAs that putatively interact with mRNAs playing key metabolic roles and that also showed differential expression before and after feeding. One case example was represented by seven miRNAs (ssc-miR-148a-3p, ssc-miR-151-3p, ssc-miR-30a-3p, ssc-miR-30e-3p, ssc-miR-421-5p, ssc-miR-493-5p and ssc-miR-503) which putatively interact with the PDK4 mRNA, one of the master regulators of glucose utilization and fatty acid oxidation.

CONCLUSIONS

As a whole, our results evidence that microRNAs are likely to play an important role in the porcine skeletal muscle metabolic adaptation to nutrient availability.

TMT-labeled quantitative proteomic analyses on the longissimus dorsi to identify the proteins underlying intramuscular fat content in pigs

The Laiwu pig is famous for its excessively extremely high level of intramuscular fat content (IMF), however, the exact regulatory mechanism underlying intramuscular fat deposition in skeletal muscle is still unknown. As an economically important trait in pigs, IMF is controlled by multiple genes and biological pathways. In this study, we performed an integrated transcriptome-assisted TMT-labeled quantitative proteomic analysis of the longissimus dorsi (LD) muscle in Laiwu pigs at the fastest IMF deposition stage and identified 5074 unique proteins and 52 differentially abundant proteins (DAPs) (>1.5-fold cutoff, p < .05). These DAPs were hierarchically clustered in the LD muscle over two developmental stages from 120 d to 240 d. A comparison between transcriptomic (mRNA) and proteomic data revealed two differentially expressed genes corresponding to the DAPs. Changes in the levels of the nine proteins were further analyzed using RT-qPCR and parallel reaction monitoring (PRM). The proteins identified in this study could serve as candidates for elucidating the molecular mechanism of IMF deposition in pigs. SIGNIFICANCE: The intramuscular fat content (IMF) refers to the amount of fat within muscles and plays an important role in meat quality by affecting meat quality-related traits, such as tenderness, juiciness and flavor. Using the integrated transcriptome-assisted TMT-labeled quantitative proteomic approach to characterize changes in the proteomic profile of the longissimus dorsi muscle, we identified differentially abundant proteins, such as ALDH1B1, OTX2, AnxA6 and Zfp512, that are associated with intramuscular fat deposition and fat biosynthesis in pigs. These proteins could serve as candidates for elucidating the molecular mechanism of IMF deposition in pigs.

Genomic diversity, linkage disequilibrium and selection signatures in European local pig breeds assessed with a high density SNP chip

Genetic characterization of local breeds is essential to preserve their genomic variability, to advance conservation policies and to contribute to their promotion and sustainability. Genomic diversity of twenty European local pig breeds and a small sample of Spanish wild pigs was assessed using high density SNP chips. A total of 992 DNA samples were analyzed with the GeneSeek Genomic Profiler (GGP) 70 K HD porcine genotyping chip. Genotype data was employed to compute genetic diversity, population differentiation and structure, genetic distances, linkage disequilibrium and effective population size. Our results point out several breeds, such as Turopolje, Apulo Calabrese, Casertana, Mora Romagnola and Lithuanian indigenous wattle, having the lowest genetic diversity, supported by low heterozygosity and very small effective population size, demonstrating the need of enhanced conservation strategies. Principal components analysis showed the clustering of the individuals of the same breed, with few breeds being clearly isolated from the rest. Several breeds were partially overlapped, suggesting genetic closeness, which was particularly marked in the case of Iberian and Alentejana breeds. Spanish wild boar was also narrowly related to other western populations, in agreement with recurrent admixture between wild and domestic animals. We also searched across the genome for loci under diversifying selection based on F_ST outlier tests. Candidate genes that may underlie differences in adaptation to specific environments and productive systems and phenotypic traits were detected in potentially selected genomic regions.

Breed, Diet, and Interaction Effects on Adipose Tissue Transcriptome in Iberian and Duroc Pigs Fed Different Energy Sources

In this study, we analyzed the effects of breed, diet energy source, and their interaction on adipose tissue transcriptome in growing Iberian and Duroc pigs. The study comprised 29 Iberian and 19 Duroc males, which were kept under identical management conditions except the nutritional treatment. Two isoenergetic diets were used with 6% high oleic sunflower oil (HO) or carbohydrates (CH) as energy sources. All animals were slaughtered after 47 days of treatment at an average live weight of 51.2 kg. Twelve animals from each breed (six fed each diet) were employed for ham subcutaneous adipose tissue RNA-Seq analysis. The data analysis was performed using two different bioinformatic pipelines. We detected 837 and 1456 differentially expressed genes (DEGs) according to breed, depending on the pipeline. Due to the strong effect of breed on transcriptome, the effect of the diet was separately evaluated in the two breeds. We identified 207 and 57 DEGs depending on diet in Iberian and Duroc pigs, respectively. A joint analysis of both effects allowed the detection of some breed–diet interactions on transcriptome, which were inferred from RNA-Seq and quantitative PCR data. The functional analysis showed the enrichment of functions related to growth and tissue development, inflammatory response, immune cell trafficking, and carbohydrate and lipid metabolism, and allowed the identification of potential regulators. The results indicate different effects of diet on adipose tissue gene expression between breeds, affecting relevant biological pathways.

Effects of lysine deficiency on carcass composition and activity and gene expression of lipogenic enzymes in muscles and backfat adipose tissue of fatty and lean piglets

The purpose of this study was to investigate potential mechanisms involved in fat deposition promoted by dietary lysine deficiency, particularly intramuscular fat (IMF), and differential responses between fatty and lean pigs. Carcass traits and lipogenic enzyme activities and gene expression levels in muscles and adipose tissue were investigated in Iberian (fatty) and Landrace × Large White (LDW) pigs under identical feeding level (based on body weight (BW)) and management conditions. Twenty-eight barrows of 10 kg initial BW, 14 per breed, were fed two isoproteic (200 g CP /kg DM) and isocaloric (14.7 MJ metabolizable energy/kg DM) diets with identical composition except for the lysine content (1.09% for diet adequate in lysine and 0.52% for diet deficient in lysine). At a BW of 25 kg, pigs were slaughtered. Compared with pigs fed the lysine-adequate diet, in both genotypes lysine-deficient diet led to lower carcass protein concentration, lower relative proportions of leaner components (loin, ham and shoulder; P < 0.01), and higher carcass fatty components and carcass lipid concentration (P < 0.001). Irrespective of diet, the activity and gene expression of lipogenic enzymes (fatty acid synthase (FAS), malic enzyme (ME) and glucose-6-phosphate dehydrogenase (G6PDH)) were greater in Iberian than in LDW pigs, particularly in adipose tissue where transcriptional regulators involved in the control of adipogenesis and lipogenesis were also upregulated in Iberian animals. In backfat tissue, there was a small decrease induced by or no effects of lysine-deficient diet on the activity and gene expression of lipogenic enzymes, nor in gene expression levels of upstream regulators of lipogenesis and adipogenesis. In longissimus muscle, the activity of FAS, G6PDH and ME increased with lysine deficiency in both genotypes (P < 0.01) and an upregulation of gene expression of lipogenic enzymes was specifically observed in Iberian pigs (P < 0.05 to P < 0.001). In biceps femoris muscle of lysine-deficient pigs, the activity of FAS and ME enzymes increased, ME1 gene was upregulated (added to FASN gene in the case of Iberian pigs; P < 0.01 to P < 0.001) and PPARA gene was downregulated (P < 0.05). The results show that in both fatty and lean pigs, the effect of lysine deficiency on lipid metabolism was tissue-specific, with an activation of lipogenesis in longissimus and biceps femoris muscle but no apparent stimulation in backfat adipose tissue. Suitable feeding protocols including lysine-deficient diets should be designed for each pig type in order to increase intramuscular lipids without penalizing the growth of lean carcass components.

Identification of Candidate Genes and Regulatory Factors Underlying Intramuscular Fat Content Through Longissimus Dorsi Transcriptome Analyses in Heavy Iberian Pigs

One of the most important determinants of meat quality is the intramuscular fat (IMF) content. The development of high-throughput techniques as RNA-seq allows identifying gene pathways and networks with a differential expression (DE) between groups of animals divergent for a particular trait. The Iberian pig is characterized by having an excellent meat quality and a high content of intramuscular fat. The objectives of the present study were to analyze the longissimus dorsi transcriptome of purebred Iberian pigs divergent for their IMF breeding value to identify differential expressed genes and regulatory factors affecting gene expression. RNA-seq allowed identifying ∼10,000 of the 25,878 annotated genes in the analyzed samples. In addition to this, 42.46% of the identified transcripts corresponded to newly predicted isoforms. Differential expression analyses revealed a total of 221 DE annotated genes and 116 DE new isoforms. Functional analyses identified an enrichment of overexpressed genes involved in lipid metabolism (FASN, SCD, ELOVL6, DGAT2, PLIN1, CIDEC, and ADIPOQ) in animals with a higher content of IMF and an enrichment of overexpressed genes related with myogenesis and adipogenesis (EGR1, EGR2, EGR3, JUNB, FOSB, and SEMA4D) in the animals with a lower content of IMF. In addition to this, potential regulatory elements of these DE genes were identified. Co-expression networks analyses revealed six long non-coding RNAs (lncRNAs) (ALDBSSCG0000002079, ALDBSSCG0000002093, ALDBSSCG0000003455, ALDBSSCG0000004244, ALDBSSCG0000005525, and ALDBSSCG0000006849) co-expressed with SEMA4D and FOSB genes and one (ALDBSSCG0000004790) with SCD, ELOVL6, DGAT2, PLIN1, and CIDEC. Analyses of the regulatory impact factors (RIFs) revealed 301 transcriptionally regulatory factors involved in expression differences, with five of them involved in adipogenesis (ARID5B, CREB1, VDR, ATF6, and SP1) and other three taking part of myogenesis and development of skeletal muscle (ATF3, KLF11, and MYF6). The results obtained provide relevant insights about the genetic mechanisms underlying IMF content in purebred Iberian pigs and a set of candidate genes and regulatory factors for further identification of polymorphisms susceptible of being incorporated in a selection program.

Effects of fetal genotype and sex on developmental response to maternal malnutrition

The present study aimed to determine whether developmental patterns, adiposity level and fatty-acid composition of fetuses exposed to maternal malnutrition are driven by their sex or their genotype, or both, as these may modulate the adaptive response to the intrauterine environment independently of the maternal genotype. We used a single maternal genotype (purebred Iberian (IB) sows), which was inseminated with heterospermic semen (obtained by mixing semen from Iberian and Large White (LW) boars), to obtain four different subsets of fetuses (male and female, purebred (IB×IB) and crossbred (IB×LW)) in Iberian purebred sows. Analysis of fetal phenotypes indicated a better adaptive response of the female offspring, which was modulated by their genotype. When faced with prenatal undernutrition, females prioritised the growth of vital organs (brain, liver, lungs, kidneys and intestine) at the expense of bone and muscle. Moreover, the analysis of fat composition showed a higher availability of essential fatty acids in the female sex than in their male counterparts and also in the Iberian genotype than in crossbred fetuses. These results are of high translational value for understanding ethnic differences in prenatal programming of postnatal health and disease status, and show evidence that prenatal development and metabolic traits are primarily determined by fetal sex and strongly modulated by fetal genotype.

Longissimus dorsi muscle transcriptomic analysis of Yunling and Chinese simmental cattle differing in intramuscular fat content and fatty acid composition

Intramuscular fat (IMF) content and fatty acid (FA) composition vary significantly across beef cattle breeds, which play an important role in taste and nutritional value. However, the molecular mechanisms underlying these phenotypic differences remain unknown. The present study compared meat quality traits between Yunling cattle and Chinese Simmental cattle. Yunling cattle showed a lower IMF content and proportion of monounsaturated fatty acids (MUFA), as well as higher proportions of saturated fatty acids (SFA), polyunsaturated fatty acids (PUFA), and short-chain fatty acids (sc-FA) in the longissimus dorsi (LD) muscle than Chinese Simmental cattle. To further identify the candidate genes and pathways responsible for these phenotypic differences, the transcriptome of LD muscle from the two breeds were measured using RNA-seq. A total of 1347 differentially expressed genes were identified. The major metabolic pathways that were differentially modulated were lipolysis and glycometabolism. Yunling cattle showed a higher expression of lipolysis genes (ALDH9A1, ACSL5, ACADM, ACAT2, ACOT2) and a lower expression of genes related to glycometabolism (PGM1, GALM, PGM1, GPI, LDHA). This research identified candidate genes and pathways for IMF content and FA composition in the LD muscle of beef cattle, which may facilitate the design of new selection strategies to improve meat quality.

Identification and Functional Analysis of Long Intergenic Non-coding RNAs Underlying Intramuscular Fat Content in Pigs

Intramuscular fat (IMF) content is an important trait that can affect pork quality. Previous studies have identified many genes that can regulate IMF. Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in various biological processes. However, lincRNAs related to IMF in pig are largely unknown, and the mechanisms by which they regulate IMF are yet to be elucidated. Here we reconstructed 105,687 transcripts and identified 1,032 lincRNAs in pig longissimus dorsi muscle (LDM) of four stages with different IMF contents based on published RNA-seq. These lincRNAs show typical characteristics such as shorter length and lower expression compared with protein-coding genes. Combined with methylation data, we found that both the promoter and genebody methylation of lincRNAs can negatively regulate lincRNA expression. We found that lincRNAs exhibit high correlation with their protein-coding neighbors in expression. Co-expression network analysis resulted in eight stage-specific modules, gene ontology and pathway analysis of them suggested that some lincRNAs were involved in IMF-related processes, such as fatty acid metabolism and peroxisome proliferator-activated receptor signaling pathway. Furthermore, we identified hub lincRNAs and found six of them may play important roles in IMF development. This work detailed some lincRNAs which may affect of IMF development in pig, and facilitated future research on these lincRNAs and molecular assisted breeding for pig.

Modulatory Effects of Breed, Feeding Status, and Diet on Adipogenic, Lipogenic, and Lipolytic Gene Expression in Growing Iberian and Duroc Pigs

Meat quality depends on tissue composition which is in turn influenced by different factors, such as diet, genotype, age, or sex. We evaluated the effects of breed, 24 h fasting, and dietary energy source (HO: oleic acid versus CH: carbohydrates) on the expression of candidate genes involved in adipogenesis, lipogenesis, and lipolysis in the adipose tissue from Iberian and Duroc growing pigs. The Iberian pigs showed greater feed intake, backfat thickness, and saturated fatty acids (SFA) content in the subcutaneous fat, whereas the Duroc pigs had greater ham weight and polyunsaturated fatty acids (PUFA) content. In both breeds, the diet induced changes in the fatty acid (FA) composition of subcutaneous fat samples. The HO group had higher monounsaturated fatty acids (MUFA) and oleic acid, and lower SFA than the CH group. Regarding gene expression, breed and feeding status (fasting versus postprandial) had significant effects on gene expression, with quantitative interactions between them, while diet showed negligible effects. In general, adipogenic and lipogenic genes were upregulated in the Iberian pigs and in postprandial samples. In contrast, the expression of lipolytic genes showed complex interaction effects. Our results agree with the phenotypic differences between the Iberian and Duroc breeds and with the inhibition of lipogenesis by fasting. Quantitative interactions between breed and feeding status effects were observed, which indicates a different response to fasting of the two breeds, with the obese Iberian breed showing a more stable expression of lipogenic genes. These results highlight the complexity of lipid metabolism regulation, especially in relation to lipolysis processes.

Identification of differentially expressed genes in longissimus dorsi muscle between Wei and Yorkshire pigs using RNA sequencing

Intramuscular fat (IMF) content is an important trait closely related to meat quality, which is highly variable among pig breeds from diverse genetic backgrounds. High-throughput sequencing has become a powerful technique for analyzing the whole transcription profiles of organisms. In order to elucidate the molecular mechanism underlying porcine meat quality, we adopted RNA sequencing to detect transcriptome in the longissimus dorsi muscle of Wei pigs (a Chinese indigenous breed) and Yorkshire pigs (a Western lean-type breed) with different IMF content. For the Wei and Yorkshire pig libraries, over 57 and 64 million clean reads were generated by transcriptome sequencing, respectively. A total of 717 differentially expressed genes (DEGs) were identified in our study (false discovery rate < 0.05 and fold change > 2), with 323 up-regulated and 394 down-regulated genes in Wei pigs compared with Yorkshire pigs. Gene Ontology analysis showed that DEGs significantly related to skeletal muscle cell differentiation, phospholipid catabolic process, and extracellular matrix structural constituent. Pathway analysis revealed that DEGs were involved in fatty acid metabolism, steroid biosynthesis, glycerophospholipid metabolism, and protein digestion and absorption. Quantitative real time PCR confirmed the differential expression of 11 selected DEGs in both pig breeds. The results provide useful information to investigate the transcriptional profiling in skeletal muscle of different pig breeds with divergent phenotypes, and several DEGs can be taken as functional candidate genes related to lipid metabolism (ACSL1, FABP3, UCP3 and PDK4) and skeletal muscle development (ASB2, MSTN, ANKRD1 and ANKRD2).

Genome-wide analysis of circular RNAs in prenatal and postnatal muscle of sheep

Circular RNAs (circRNAs), a type of non-coding RNA with circular structure, were generated by back splicing and widely expressed in animals and plants. Recent studies have shown that circRNAs extensively participate in cell proliferation, cell differentiation, cell autophagy and other biological processes. However, the role and expression of circRNAs in the development and growth of muscle have not been studied in sheep. In our study, we first used RNA-seq to study the circRNAs in prenatal and postnatal longissimus dorsi muscle of sheep. A total of 6113 circRNAs were detected from the RNA-seq data. Several circRNAs were identified using reverse transcription PCR, DNA sequencing and RNase R digestion experiments. The expression levels of several circRNAs in prenatal and postnatal muscle were confirmed by Real-Time RT-PCR. The gene ontology (GO) and KEGG enrichment analysis of the host gene of the circRNAs showed that these circRNAs were mainly involved in the growth and development of muscle related signaling pathways. These circRNAs might sponge microRNAs (miRNAs) in predicted circRNA-miRNA-mRNA networks. The circRNAs expression profiles in muscle provided an important reference for the study of circRNAs in sheep.

Dynamic transcriptome and DNA methylome analyses on longissimus dorsi to identify genes underlying intramuscular fat content in pigs

Background

The intramuscular fat content (IMF) refers to the amount of fat within muscles, including the sum of phospholipids mainly found in cell membranes, triglycerides and cholesterol, and is determined both by hyperplasia and hypertrophy of adipocyte during the development of pigs. The IMF content is an important economic trait that is genetically controlled by multiple genes. The Laiwu pig is an indigenous fatty pig breed distributed in North China, characterized by excessively higher level of IMF content (9%~12%), therefore, is suitable for the identification of genes controlling IMF variations. To identify genes underlying IMF deposition, we performed genome-wide transcriptome and methylome analyses on longissimus dorsi (LD) muscle in Laiwu pigs across four developmental stages.

Results

A total of 22,524 expressed genes were detected and 1158 differentially expressed genes (DEGs) were hierarchically clustered in the LD muscle over four developmental stages from 60 d to 400 d. These genes were significantly clustered into four temporal expression profiles, and genes participating in fat cell differentiation and lipid biosynthesis processes were identified. From 120 d to 240 d, the period with the maximum IMF deposition rate, the lipid biosynthesis related genes (FOSL1, FAM213B and G0S2), transcription factors (TFs) (EGR1, KLF5, SREBF2, TP53 and TWIST1) and enriched pathways (steroid biosynthesis and fatty acid biosynthesis) were revealed; and fat biosynthesis relevant genes showing differences in DNA methylation in gene body or intergenic region were detected, such as FASN, PVALB, ID2, SH3PXD2B and EGR1.

Conclusions

This study provides a comprehensive landscape of transcriptome of the LD muscle in Laiwu pigs ranging from 60 to 400 days old, and methylome of the LD muscle in 120 d and 240 d Laiwu pigs. A set of candidate genes and TFs involved in fat biosynthesis process were identified, which were probably responsible for IMF deposition. The results from this study would provide a reference for the identification of genes controlling IMF variation, and for exploring molecular mechanisms underlying IMF deposition in pigs.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4201-9) contains supplementary material, which is available to authorized users.

Transcriptomic profiling in muscle and adipose tissue identifies genes related to growth and lipid deposition

Growth performance and meat quality are important traits for the pig industry and consumers. Adipose tissue is the main site at which fat storage and fatty acid synthesis occur. Therefore, we combined high-throughput transcriptomic sequencing in adipose and muscle tissues with the quantification of corresponding phenotypic features using seven Chinese indigenous pig breeds and one Western commercial breed (Yorkshire). We obtained data on 101 phenotypic traits, from which principal component analysis distinguished two groups: one associated with the Chinese breeds and one with Yorkshire. The numbers of differentially expressed genes between all Chinese breeds and Yorkshire were shown to be 673 and 1056 in adipose and muscle tissues, respectively. Functional enrichment analysis revealed that these genes are associated with biological functions and canonical pathways related to oxidoreductase activity, immune response, and metabolic process. Weighted gene coexpression network analysis found more coexpression modules significantly correlated with the measured phenotypic traits in adipose than in muscle, indicating that adipose regulates meat and carcass quality. Using the combination of differential expression, QTL information, gene significance, and module hub genes, we identified a large number of candidate genes potentially related to economically important traits in pig, which should help us improve meat production and quality.

Transcriptome analysis of the pectoral muscles of local chickens and commercial broilers using Ribo-Zero ribonucleic acid sequencing

BACKGROUND

The molecular mechanisms underlying meat quality and muscle growth are not clear. The meat quality and growth rates of local chickens and commercial broilers are very different. The Ribo-Zero RNA-Seq technology is an effective means of analyzing transcript groups to clarify molecular mechanisms. The aim of this study was to provide a reference for studies of the differences in the meat quality and growth of different breeds of chickens.

RESULTS

Ribo-Zero RNA-Seq technology was used to analyze the pectoral muscle transcriptomes of Gushi chickens and AA broilers. Compared with AA broilers, 1649 genes with annotated information were significantly differentially expressed (736 upregulated and 913 downregulated) in Gushi chickens with Q≤0.05 (Q is the P-value corrected by multiple assumptions test) at a fold change ≥2 or ≤0.5. In addition, 2540 novel significantly differentially expressed (SDE) genes (1405 upregulated and 1135 downregulated) were discovered. The results showed that the main signal transduction pathways that differed between Gushi chickens and AA broilers were related to amino acid metabolism. Amino acids are important for protein synthesis, and they regulate key metabolic pathways to improve the growth, development and reproduction of organisms.

CONCLUSION

This study showed that differentially expressed genes in the pectoral tissues of Gushi chickens and AA broilers were related to fat metabolism, which affects meat. Additionally, a large number of novel genes were found that may be involved in fat metabolism and thus may affect the formation of meat, which requires further study. The results of this study provide a reference for further studies of the molecular mechanisms of meat formation.