Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition

Integration of Multiomics Data Reveals Selection Characteristics of ITGB1 That Are Associated with Size Differentiation in Pigs

Min pigs, a prominent local breed from Northeast China, have diverged into two distinct breeds, Ermin (EM) pigs and Hebao (HB) pigs, through prolonged natural and artificial selection. Although these two breeds exhibit distinct differences in body size, the genetic mechanisms underlying this variation remain poorly understood. In this study, we performed whole-genome resequencing and transcriptome analysis on EM and HB pigs to elucidate the genetic basis of body size variation in Min pigs through genome-wide selection signal analysis and the identification of differentially expressed genes (DEGs). The analysis of genetic diversity and population genetic structure across 14 pig breeds revealed that, compared with other breeds, Min pigs present relatively high genetic diversity and a unique genetic structure. Notably, EM pigs exhibited significant genetic differentiation from HB pigs. Integrated analysis of whole-genome resequencing and transcriptome data revealed candidate genes associated with body size variation in Min pigs, including ENPP1, ENPP3, SPP1, CLU, ITGA11, ITGB1, IQGAP2, BMP7, and F2RL2. These genes are enriched primarily in pathways related to ECM-receptor interactions; pantothenate and CoA biosynthesis; starch and sucrose metabolism; nicotinate and nicotinamide metabolism; pyrimidine metabolism; nucleotide metabolism; cellular responses to lipids; biomineral tissue development; biomineralization; and other pathways related to cell signaling, metabolic responses, lipid deposition, and skeletal development. Notably, ITGB1 on chromosome 10 showed strong positive selection in EM pigs, with an SNP locus exhibiting a significant G/A allele frequency difference between EM pigs (G = 52.94%, A = 47.06%) and HB pigs (G = 0%, A = 100%). Our findings suggest that Min pigs potentially modulate lipid metabolism efficiency in adipose tissue through variations in the expression of the ITGB1 gene, potentially contributing to body size differences. These results provide new insights into the genetic mechanisms underlying body size variation in domestic pigs and serve as a valuable reference for identifying and breeding pig breeds with distinct body sizes.

Integrated transcriptomic and metabolomic analyses reveal heterosis for meat quality of Neijiang pigs

Obese pig breeds have excellent meat quality, while lean pig breeds have high lean meat percentage and feed conversion rate. However, due to their respective shortcomings, obese pig and lean pig breeds are unable to balance production and consumption needs. Therefore, this study crossbred the obese Chinese pig breed Neijiang (NJ) with lean type Large White pigs (LW) to produce Neijiang × Large White(NL) pigs. This study compared the differences in carcass and meat quality traits between NJ pigs and NL pigs, and for the first time comprehensively analyzed the longissimus dorsi muscle of NJ pigs and NL pigs using transcriptomics and metabolomics. The results of slaughter and meat quality testing indicate that the carcass performance of NL pigs was significantly higher than that of NJ pigs, and the excellent meat quality characteristics of NJ pigs were also retained on NL pigs. The results of transcriptomics and metabolomics showed that there were 635 differentially expressed genes (DEGs) and 11 significantly different metabolites (SDM) in the longissimus dorsi muscle of NJ and NL pigs. The results of multi omics joint analysis showed that betaine, uridine triphosphate, glycerol 3-phosphate, and glutathione in SDMs were enriched in the shared KEGG pathway and significantly correlated with C1QTNF12, GGA3, SLC16A6, and RXRG in DEGs. In general, it is feasible to enhance the production performance of NJ pigs through crossbreeding with LW pigs. The hybrid offspring inherit the advantages of these two varieties, maintaining excellent meat quality while also having better carcass performance.

Integrated analyses of the intestinal microbiome and transcriptome in Ningxiang piglets

Ningxiang (NX) pig has been recognized as one of the most famous Chinese indigenous breeds due to its characteristics in stress resistance. However, intestinal microbial feature and gene profiling in NX piglets have not been studied. Here, we compared the intestinal microbiome and transcriptome between NX and Duroc × Landrace × Large white (DLY) piglets and found the high enrichment of several colonic Bacteroides, Prevotella and Clostridium species in NX piglets. Further functional analyses revealed their predominant function in methane, glycolysis and gluconeogenesis metabolism. Our mRNA-sequencing data unraveled the distinct colonic gene expression between these two breeds. In particular, we showed that the improved intestinal function in NX piglets may be determined by enhanced intestinal barrier gene expression and varied immune gene expression through modulating the composition of the gut microbes. Together, our study revealed the intestinal characteristics of NX piglets, providing their potential application in improving breeding strategies and developing dietary interventions.

H3K27ac modification and transcription characteristics of adipose and muscle tissues in Chuxiang Black pig

The establishment of high-quality pork breeds for improving meat quality in the pig industry is needed. The Chuxiang Black (CX) pig is a new breed developed from Chinese local pigs and Western lean pigs that has a high proportion of lean meat and excellent meat quality. However, the characteristics of cis-regulatory elements in CX pigs are still unknown. In this study, cis-regulatory elements of muscle and adipose tissues in CX pigs were investigated using ChIP-seq and RNA sequencing. Compared with the reported cis-regulatory elements of muscle and adipose tissues, 1768 and 1012 highly activated enhancers and 433 and 275 highly activated promoters in CX muscle and adipose tissues were identified, respectively. Motif analysis showed that transcription factors, such as MEF2A and MEF2C, were core regulators of highly activated enhancers and promoters in muscle. Similarly, the transcription factors JUNB and CUX1 were identified as essential for highly activated enhancers and promoters in CX adipose tissue. These results enrich the resources for the analysis of cis-regulatory elements in the pig genome and provide new basic data for further meat quality improvement through breeding in pigs.

Genome-Wide Transcriptome Profiling Reveals the Mechanisms Underlying Hepatic Metabolism under Different Raising Systems in Yak

Yak meat is nutritionally superior to beef cattle but has a low fat content and is slow-growing. The liver plays a crucial role in lipid metabolism, and in order to determine whether different feeding modes affect lipid metabolism in yaks and how it is regulated, we employed RNA sequencing (RNA-seq) technology to analyze the genome-wide differential gene expression in the liver of yaks maintained under different raising systems. A total of 1663 differentially expressed genes (DEGs) were identified (|log2FC| ≥ 0 and p-value ≤ 0.05), including 698 down-regulated and 965 up-regulated genes. According to gene ontology (GO) and KEGG enrichment analyses, these DEGs were significantly enriched in 13 GO terms and 26 pathways (p < 0.05). Some DEGs were enriched in fatty acid degradation, PPAR, PI3K-Akt, and ECM receptor pathways, which are associated with lipid metabolism. A total of 16 genes are well known to be related to lipid metabolism (e.g., APOA1, FABP1, EHHADH, FADS2, SLC27A5, ACADM, CPT1B, ACOX2, HMGCS2, PLIN5, ACAA1, IGF1, FGFR4, ALDH9A1, ECHS1, LAMA2). A total of 11 of the above genes were significantly enriched in the PPAR signaling pathway. The reliability of the transcriptomic data was verified using qRT-PCR. Our findings provide new insights into the mechanisms regulating yak meat quality. It shows that fattening improves the expression of genes that regulate lipid deposition in yaks and enhances meat quality. This finding will contribute to a better understanding of the various factors that determine yak meat quality and help develop strategies to improve yield and quality.

Visfatin Affects the Transcriptome of Porcine Luteal Cells during Early Pregnancy

Visfatin/NAMPT (VIS), the hormone exerting a pleiotropic effect, is also perceived as an important factor in the regulation of reproductive processes and pregnancy maintenance. Previous studies confirmed its involvement in the control of porcine pituitary and ovary function. In this study, we hypothesized that VIS may affect the global transcriptome of luteal cells and thus regulate the functioning of the ovaries. Illumina's NovaSeq 6000 RNA sequencing was performed to investigate the differentially expressed genes (DEGs) and long non-coding RNAs (DELs) as well as the occurrence of differential alternative splicing events (DASs) in the porcine luteal cells exposed to VIS (100 ng/mL) during the implantation period. The obtained results revealed 170 DEGs (99 up- and 71 downregulated) assigned to 45 functional annotations. Moreover, we revealed 40 DELs, of which 3 were known and 37 were described for the first time. We identified 169 DASs events. The obtained results confirmed a significant effect of VIS on the transcriptome and spliceosome of luteal cells, including the genes involved in the processes crucial for successful implantation and pregnancy maintenance as angiogenesis, steroidogenesis, inflammation, cell development, migration, and proliferation.

Identification of candidate regulatory genes for intramuscular fatty acid composition in pigs by transcriptome analysis

BACKGROUND

Intramuscular fat (IMF) content and its fatty acid (FA) composition are typically controlled by several genes, each with a small effect. In the current study, to pinpoint candidate genes and putative regulators involved in FA composition, we performed a multivariate integrative analysis between intramuscular FA and transcriptome profiles of porcine longissimus dorsi (LD) muscle. We also carried out a combination of network, regulatory impact factor (RIF), in silico prediction of putative target genes, and functional analyses to better support the biological relevance of our findings.

RESULTS

For this purpose, we used LD RNA-Seq and intramuscular FA composition profiles of 129 Iberian × Duroc backcrossed pigs. We identified 378 correlated variables (13 FA and 365 genes), including six FA (C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, C18:0, and C16:1n-7) that were among the most interconnected variables in the predicted network. The detected FA-correlated genes include genes involved in lipid and/or carbohydrate metabolism or in regulation of IMF deposition (e.g., ADIPOQ, CHUK, CYCS, CYP4B1, DLD, ELOVL6, FBP1, G0S2, GCLC, HMGCR, IDH3A, LEP, LGALS12, LPIN1, PLIN1, PNPLA8, PPP1R1B, SDR16C5, SFRP5, SOD3, SNW1, and TFRC), meat quality (GALNT15, GOT1, MDH1, NEU3, PDHA1, SDHD, and UNC93A), and transport (e.g., EXOC7 and SLC44A2). Functional analysis highlighted 54 over-represented gene ontology terms, including well-known biological processes and pathways that regulate lipid and carbohydrate metabolism. RIF analysis suggested a pivotal role for six transcription factors (CARHSP1, LBX1, MAFA, PAX7, SIX5, and TADA2A) as putative regulators of gene expression and intramuscular FA composition. Based on in silico prediction, we identified putative target genes for these six regulators. Among these, TADA2A and CARHSP1 had extreme RIF scores and present novel regulators in pigs. In addition, the expression of TADA2A correlated (either positively or negatively) with C20:4n-6, C18:2n-6, C20:3n-6, C18:1n-9, and that of CARHSP1 correlated (positively) with the C16:1n-7 lipokine. We also found that these two transcription factors share target genes that are involved in lipid metabolism (e.g., GOT1, PLIN1, and TFRC).

CONCLUSIONS

This integrative analysis of muscle transcriptome and intramuscular FA profile revealed valuable information about key candidate genes and potential regulators for FA and lipid metabolism in pigs, among which some transcription factors are proposed to control gene expression and modulate FA composition differences.

Identification of Dezhou donkey muscle development-related genes and long non-coding RNA based on differential expression analysis

Long non-coding RNAs (lncRNAs) play a critical role in the development of muscles. However, the role of lncRNAs in regulating skeletal muscle development has not been studied systematically in the donkey. In this study, we performed the RNA sequencing for different stages of muscles in donkeys, and investigate their expression profile, which showed that 3215 mRNAs (p-adjust <0.05) and 471 lncRNAs (p-value <0.05) were significantly differently expressed (DE) verified by RT-qPCR. GO and KEGG enrichment analysis indicated that DE genes and target genes of DE lncRNAs were associated with muscle development in the donkey. We also found these four target genes (DCN, ITM2A, MUSTN1, ARRDC2) involved in skeletal muscle growth and development. Combined with transcriptome data, network, and RT-qPCR results showed that four co-expression networks of DCN and lnc-008278, ITM2A and lnc_017247, MUSTN1 and lnc_030153, and ARRDC2 and lnc_033914, which may play an important role in the formation and development of muscle in the donkey.

Transcriptome analysis of divergent residual feed intake phenotypes in the M. longissimus thoracis et lumborum of Wannan Yellow rabbits

Introduction: Feed efficiency is an important economic trait in rabbit meat production. The identification of molecular mechanisms and candidate genes for feed efficiency may improve the economic and environmental benefits of the rabbit meat industry. As an alternative to the conventional feed conversion ratio, residual feed intake (RFI) can be used as an accurate indicator of feed efficiency. Methods: RNA sequencing was used to identify the differentially expressed genes (DEGs) in the M. longissimus thoracis et lumborum of eight Wannan Yellow rabbits with excessively high or low RFIs (HRFI or LRFI, respectively). Thereafter, Gene Ontology (GO) analysis, enrichment using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, gene set enrichment analysis (GSEA), and protein-protein interaction (PPI) network analysis was conducted. Results: In total, 445 DEGs were identified in the M. longissimus thoracis et lumborum of rabbits with high and low RFIs. The significantly enriched GO terms identified in these two groups were primarily involved in energy and mitochondrial metabolism and oxidation-reduction processes. KEGG analysis identified 11 significantly enriched pathways, including oxidative phosphorylation, PI3K-Akt signaling, and extracellular matrix-receptor interaction pathways. According to GSEA, the expressions of genes and pathways related to mitochondrial function were upregulated in HRFI rabbits, whereas genes with upregulated expressions in LRFI rabbits were related to immune response and energy metabolism. Additionally, PPI network analysis revealed five potential candidate genetic markers. Conclusion: Comparative analysis of the M. longissimus thoracis et lumborum transcriptomes in HRFI and LRFI rabbits revealed FOS, MYC, PRKACB, ITGA2, and FN1 as potential candidate genes that affect feed efficiency in rabbits. In addition, key signaling pathways involved in oxidative phosphorylation and PI3K-Akt and ECM-receptor interaction signaling impact rabbit feed efficiency. These findings will aid in breeding programs to improve feed efficiency and optimize RFI selection of rabbits for meat production.

Transcriptome profile of skeletal muscle using different sources of dietary fatty acids in male pigs

Pork is of great importance in world trade and represents the largest source of fatty acids in the human diet. Lipid sources such as soybean oil (SOY), canola (CO), and fish oil (FO) are used in pig diets and influence blood parameters and the ratio of deposited fatty acids. In this study, the main objective was to evaluate changes in gene expression in porcine skeletal muscle tissue resulting from the dietary oil sources and to identify metabolic pathways and biological process networks through RNA-Seq. The addition of FO in the diet of pigs led to intramuscular lipid with a higher FA profile composition of C20:5 n-3, C22:6 n-3, and SFA (C16:0 and C18:0). Blood parameters for the FO group showed lower cholesterol and HDL content compared with CO and SOY groups. Skeletal muscle transcriptome analyses revealed 65 differentially expressed genes (DEG, FDR 10%) between CO vs SOY, and 32 DEG for CO vs FO, and 531 DEG for SOY vs FO comparison. Several genes, including AZGP1, PDE3B, APOE, PLIN1, and LIPS, were found to be down-regulated in the diet of the SOY group compared to the FO group. The enrichment analysis revealed DEG involved in lipid metabolism, metabolic diseases, and inflammation between the oil groups, with specific gene functions in each group and altered blood parameters. The results provide mechanisms to help us understand the behavior of genes according to fatty acids.

RNA-seq transcriptome profiling of pigs' liver in response to diet with different sources of fatty acids

Pigs (Sus scrofa) are an animal model for metabolic diseases in humans. Pork is an important source of fatty acids (FAs) in the human diet, as it is one of the most consumed meats worldwide. The effects of dietary inclusion of oils such as canola, fish, and soybean oils on pig gene expression are mostly unknown. Our objective was to evaluate FA composition, identify changes in gene expression in the liver of male pigs fed diets enriched with different FA profiles, and identify impacted metabolic pathways and gene networks to enlighten the biological mechanisms' variation. Large White male pigs were randomly allocated to one of three diets with 18 pigs in each; all diets comprised a base of corn and soybean meal to which either 3% of soybean oil (SOY), 3% canola oil (CO), or 3% fish oil (FO) was added for a 98-day trial during the growing and finishing phases. RNA sequencing was performed on the liver samples of each animal by Illumina technology for differential gene expression analyses, using the R package DESeq2. The diets modified the FA profile, mainly in relation to polyunsaturated and saturated FAs. Comparing SOY vs. FO, 143 differentially expressed genes (DEGs) were identified as being associated with metabolism, metabolic and neurodegenerative disease pathways, inflammatory processes, and immune response networks. Comparing CO vs. SOY, 148 DEGs were identified, with pathways related to FA oxidation, regulation of lipid metabolism, and metabolic and neurodegenerative diseases. Our results help explain the behavior of genes with differential expression in metabolic pathways resulting from feeding different types of oils in pig diets.

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.

Exon-intron split analysis reveals posttranscriptional regulatory signals induced by high and low n-6/n-3 polyunsaturated fatty acid ratio diets in piglets

Polyunsaturated fatty acids (PUFA), such as omega-6 (n-6) and omega-3 (n-3), play a vital role in nutrient metabolism, inflammatory response, and gene regulation. microRNAs (miRNA), which can potentially degrade targeted messenger RNAs (mRNA) and/or inhibit their translation, might play a relevant role in PUFA-related changes in gene expression. Although differential expression analyses can provide a comprehensive picture of gene expression variation, they are unable to disentangle when in the mRNA life cycle the regulation of expression is taking place, including any putative functional miRNA-driven repression. To capture this, we used an exon-intron split analysis (EISA) approach to account for posttranscriptional changes in response to extreme values of n-6/n-3 PUFA ratio. Longissimus dorsi muscle samples of male and female piglets from sows fed with n-6/n-3 PUFA ratio of 13:1 (SOY) or 4:1 (LIN), were analyzed in a bidirectional contrast (LIN vs. SOY, SOY vs. LIN). Our results allowed the identification of genes showing strong posttranscriptional downregulation signals putatively targeted by significantly upregulated miRNA. Moreover, we identified genes primarily involved in the regulation of lipid-related metabolism and immune response, which may be associated with the pro- and anti-inflammatory functions of the n-6 and n-3 PUFA, respectively. EISA allowed us to uncover regulatory networks complementing canonical differential expression analyses, thus providing a more comprehensive view of muscle metabolic changes in response to PUFA concentration.

Effective Quality Breeding Directions-Comparison and Conservative Analysis of Hepatic Super-Enhancers between Chinese and Western Pig Breeds

The transcriptional initiation of genes is closely bound to the functions of cis-regulatory elements, including promoters, typical enhancers (TEs), and recently-identified super-enhancers (SEs). In this study, we identified these cis-regulatory elements in the livers of two Chinese (Meishan and Enshi Black) and two Western (Duroc and Large White) pig breeds using ChIP-seq data, then explored their similarities and differences. In addition, we analyzed the conservation of SEs among different tissues and species (pig, human, and mouse). We observed that SEs were more significantly enriched by transcriptional initiation regions, TF binding sites, and SNPs than other cis-elements. Western breeds included fewer SEs in number, while more growth-related QTLs were associated with these SEs. Additionally, the SEs were highly tissue-specific, and were conserved in the liver among humans, pigs, and mice. We concluded that intense selection could concentrate functional SEs; thus, SEs could be applied as effective detection regions in genomic selection breeding.

Analysis of lncRNA in the skeletal muscle of rabbits at different developmental stages

It is universally acknowledged that lncRNA plays an important role in the regulation of animal skeletal muscle development regulation. However, there is a lack of relevant research on lncRNA in rabbit skeletal muscle development. Thus, we explored the expression profiles of lncRNA in rabbits at three growth stages (2-week-old fetus, 6-week-old post-weaning, and 6-month-old adult) using RNA-seq. A total of 554 differentially expressed lncRNAs (235 up- and 319 down-regulated) were found between the post-weaning and fetus groups and 19 (7 up- and 12 down-regulated) between the post-weaning and adult groups and 429 (115 up- and 314 down-regulated) between the fetus and adult. The enrichment pathways in the post-weaning and fetus groups were mainly concentrated at AMPK and PI3K-Akt signaling pathways, and the co-expression results revealed that LINC-2903, LINC-2374, LINC-8591 plays a role in early maintenance of skeletal muscle development. The enriched pathways in the fetus and adult groups were mainly involved in PI3K-Akt signaling pathways with a strong association found in mTOR signaling pathways. Analysis of the co-expression results suggests that LINC-5617 may be involved in the proliferation of embryonic skeletal muscle cells, and that LINC-8613 and LINC-8705 may provide energy for postnatal skeletal muscle development. The specific roles of different lncRNAs in different developmental stages of New Zealand White rabbits obtained. This will contribute to the subsequent study on the regulatory mechanism of muscle development in New Zealand White rabbits.

Chimeric RNA TNNI2-ACTA1-V1 Regulates Cell Proliferation by Regulating the Expression of NCOA3

Chimeric RNA is a crucial target for tumor diagnosis and drug therapy, also having its unique biological role in normal tissues. TNNI2-ACTA1-V1 (TA-V1), a chimeric RNA discovered by our laboratory in porcine muscle tissue, can inhibit the proliferation of Porcine Skeletal Muscle Satellite Cells (PSCs). The regulatory mechanism of TA-V1 in PSCs remains unclear, but we speculate that NCOA3, DDR2 and RDX may be the target genes of TA-V1. In this study, we explored the effects of NCOA3, DDR2 and RDX on cell viability and cell proliferation by CCK-8 assay, EdU staining and flow cytometry. Furthermore, the regulatory pathway of proliferation in PSCs mediated by TA-V1 through NCOA3 or CyclinD1 was elucidated by co-transfection and co-immunoprecipitation (Co-IP). The results revealed that overexpression of NCOA3 significantly increased cell viability and the expression level of CyclinD1, and also promotes cell proliferation by changing cells from the G1 phase to the S phase. In addition, inhibiting the expression of NCOA3 substantially reduced cell viability and inhibited cell proliferation. Overexpression of DDR2 and RDX had no significant effect on cell viability and proliferation. Co-transfection experiments showed that NCOA3 could rescue the proliferation inhibition of PSCs caused by TA-V1. Co-IP assay indicated that TA-V1 directly interacts with NCOA3. Our study explores the hypothesis that TA-V1 directly regulates NCOA3, indirectly regulating CyclinD1, thereby regulating PSCs proliferation. We provide new putative mechanisms of porcine skeletal muscle growth and lay the foundation for the study of chimeric RNA in normal tissues.

Differential Gene Expression Associated with Soybean Oil Level in the Diet of Pigs

The aim of this study was to identify the differentially expressed genes (DEG) from the skeletal muscle and liver samples of animal models for metabolic diseases in humans. To perform the study, the fatty acid (FA) profile and RNA sequencing (RNA-Seq) data of 35 samples of liver tissue (SOY1.5, n = 17 and SOY3.0, n = 18) and 36 samples of skeletal muscle (SOY1.5, n = 18 and SOY3.0, n = 18) of Large White pigs were analyzed. The FA profile of the tissues was modified by the diet, mainly those related to monounsaturated (MUFA) and polyunsaturated (PUFA) FA. The skeletal muscle transcriptome analysis revealed 45 DEG (FDR 10%), and the functional enrichment analysis identified network maps related to inflammation, immune processes, and pathways associated with oxidative stress, type 2 diabetes, and metabolic dysfunction. For the liver tissue, the transcriptome profile analysis revealed 281 DEG, which participate in network maps related to neurodegenerative diseases. With this nutrigenomics study, we verified that different levels of soybean oil in the pig diet, an animal model for metabolic diseases in humans, affected the transcriptome profile of skeletal muscle and liver tissue. These findings may help to better understand the biological mechanisms that can be modulated by the diet.

Transcriptome Sequencing-Based Mining of Genes Associated With Pubertal Initiation in Dolang Sheep

Improving the fertility of sheep is an important goal in sheep breeding as it greatly increases the productivity. Dolang sheep is a typical representative breed of lamb in Xinjiang and is the main local sheep breed and meat source in the region. To explore the genes associated with the initiation of puberty in Dolang sheep, the hypothalamic tissues of Dolang sheep prepubertal, pubertal, and postpubertal periods were collected for RNA-seq analysis on the Illumina platform, generating 64.08 Gb clean reads. A total of 575, 166, and 648 differentially expressed genes (DEGs) were detected in prepuberty_vs._puberty, postpuberty_vs._prepuberty, and postpuberty_vs._puberty analyses, respectively. Based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses, the related genes involved in the initiation of puberty in Dolang sheep were mined. Ten genes that have direct or indirect functions in the initiation of puberty in Dolang sheep were screened using the GO and KEGG results. Additionally, quantitative real-time PCR was used to verify the reliability of the RNA-Seq data. This study provided a new approach for revealing the mechanism of puberty initiation in sheep and provided a theoretical basis and candidate genes for the breeding of early-pubertal sheep by molecular techniques, and at the same time, it is also beneficial for the protection, development, and utilization of the fine genetic resources of Xinjiang local sheep.

Association of IFNA16 and TNFRSF19 Polymorphisms with Intramuscular Fat Content and Fatty Acid Composition in Pigs

Interferon-alpha-16 (IFNA16) and tumor necrosis factor receptor superfamily member 19 (TNFRSF19) are cytokines that may play a role in adipogenesis and fatness. Single nucleotide polymorphisms (SNPs) of the porcine IFNA16 and TNFRSF19 genes were verified and their association with intramuscular fat (IMF) content and fatty acid (FA) composition were evaluated in commercial crossbred pigs. Two non-synonymous SNPs of the porcine IFNA16 c.413G > A and TNFRSF19 c.860G > C loci were detected in commercial crossbred pigs. The porcine IFNA16 c.413G >A polymorphism was significantly associated with stearic acid, total saturated FAs (SFAs), and the ratio of monounsaturated FAs (MUFAs) to SFAs (p < 0.05). Furthermore, the porcine TNFRSF19 c.860G > C polymorphism was found to be significantly associated with IMF content and arachidic acid levels (p < 0.05). The results revealed that porcine IFNA16 and TNFRSF19 polymorphisms are related to IMF content and/or FA composition and affirmed the importance of these cytokine genes as potential candidate genes for lipid deposition and FA composition in the muscle tissue of pigs.

Hepatic transcriptome analysis identifies genes, polymorphisms and pathways involved in the fatty acids metabolism in sheep

Fatty acids (FA) in ruminants, especially unsaturated FA (USFA) have important impact in meat quality, nutritional value, and flavour quality of meat, and on consumer’s health. Identification of the genetic factors controlling the FA composition and metabolism is pivotal to select sheep that produce higher USFA and lower saturated (SFA) for the benefit of sheep industry and consumers. Therefore, this study was aimed to investigate the transcriptome profiling in the liver tissues collected from sheep with divergent USFA content in longissimus muscle using RNA deep-sequencing. From sheep (n = 100) population, liver tissues with higher (n = 3) and lower (n = 3) USFA content were analysed using Illumina HiSeq 2500. The total number of reads produced for each liver sample were ranged from 21.28 to 28.51 million with a median of 23.90 million. Approximately, 198 genes were differentially regulated with significance level of p-adjusted value <0.05. Among them, 100 genes were up-regulated, and 98 were down-regulated (p<0.01, FC>1.5) in the higher USFA group. A large proportion of key genes involved in FA biosynthesis, adipogenesis, fat deposition, and lipid metabolism were identified, such as APOA5, SLC25A30, GFPT1, LEPR, TGFBR2, FABP7, GSTCD, and CYP17A. Pathway analysis revealed that glycosaminoglycan biosynthesis- keratan sulfate, adipokine signaling, galactose metabolism, endocrine and other factors-regulating calcium metabolism, mineral metabolism, and PPAR signaling pathway were playing important regulatory roles in FA metabolism. Importantly, polymorphism and association analyses showed that mutation in APOA5, CFHR5, TGFBR2 and LEPR genes could be potential markers for the FA composition in sheep. These polymorphisms and transcriptome networks controlling the FA variation could be used as genetic markers for FA composition-related traits improvement. However, functional validation is required to confirm the effect of these SNPs in other sheep population in order to incorporate them in the sheep breeding program.

Identification of Chimeric RNAs in Pig Skeletal Muscle and Transcriptomic Analysis of Chimeric RNA TNNI2-ACTA1 V1

Chimeric RNA was considered a special marker of cancer. However, recent studies have demonstrated that chimeric RNAs also exist in non-cancerous cells and tissues. Here, we analyzed and predicted jointly 49 chimeric RNAs by Star-Fusion and FusionMap. One chimeric RNA, we named TNNI2-ACTA1, and its eight transcript variants were identified by reverse transcriptase–polymerase chain reaction. The overexpression of TNNI2-ACTA1 V1 inhibited the proliferation of porcine skeletal muscle satellite cells through down-regulating the mRNA expression levels of cell cycle–related genes cyclinD1. However, as parental genes, there is no such effect in the TNNI2 and ACTA1. To explore the underlying mechanism for this phenomenon, we used RNA-seq to profile the transcriptomes of PSCs with overexpression. Compared with the negative control group, 1,592 differentially expressed genes (DEGs) were upregulated and 1,077 DEGs downregulated in TNNI2 group; 1,226 DEGs were upregulated and 902 DEGs downregulated in ACTA1 group; and 13 DEGs were upregulated and 16 DEGs downregulated in TNNI2-ACTA1 V1 group, respectively. Compared with the parental gene groups, three specific genes were enriched in the TNNI2-ACTA1 V1 group (NCOA3, Radixin, and DDR2). These three genes may be the key to TNNI2-ACTA1 V1 regulating cell proliferation. Taken together, our study explores the role of chimeric RNAs in normal tissues. In addition, our study as the first research provides the foundation for the mechanism of chimeric RNAs regulating porcine skeletal muscle growth.

Transcriptomics analysis of Daheng broilers reveals that PLIN2 regulates chicken preadipocyte proliferation, differentiation and apoptosis

BACKGROUND

Intramuscular fat content, an important meat quality trait, strongly affects flavor, juiciness, and tenderness. Sex hormones regulate lipid metabolism, and female hormones stimulate fat deposition, thereby making the female chickens always fatter than males. In this study, the effect of sex on IMF deposition was screened following transcriptomics in chickens.

METHODS AND RESULTS

Results confirmed significantly higher IMF content of 150-day female chickens as compared to the male chickens. The female chickens manifested higher serum TG, LDL-C, and VLDL, and significantly lower HDL-C contents than male chickens. Moreover, differential expression of genes involved in lipid metabolism were obtained in the muscle and liver between female and male chicken, which could partly interpret the possible reasons for the sex-mediated differences of IMF content. Cellular results revealed that inhibition of PLIN2 significantly inhibited chicken preadipocyte proliferation and induces apoptosis of preadipocytes, as well as promoted adipocyte differentiation.

CONCLUSIONS

According to our results, PLIN2 may be considered as a molecular marker for poultry meat quality and applying this gene in early breed selection.

The Impact of Maternal and Piglet Low Protein Diet and Their Interaction on the Porcine Liver Transcriptome around the Time of Weaning

Maternal diet during early gestation affects offspring phenotype, but it is unclear whether maternal diet during late gestation influences piglet metabolism. We evaluated the impact of two dietary protein levels in sow late gestation diet and piglet nursery diet on piglet metabolism. Diets met or exceeded the crude protein and amino acid requirements. Sows received either 12% (Lower, L) or 17% (Higher, H) crude protein (CP) during the last five weeks of gestation, and piglets received 16.5% (L) or 21% (H) CP from weaning at age 3.5 weeks. This resulted in a 2 × 2 factorial design with four sow/piglet diet treatment groups: HH and LL (match), HL and LH (mismatch). Piglet hepatic tissues were sampled and differentially expressed genes (DEGs) were determined by RNA sequencing. At age 4.5 weeks, 25 genes were downregulated and 22 genes were upregulated in the mismatch compared to match groups. Several genes involved in catabolic pathways were upregulated in the mismatch compared to match groups, as were genes involved in lipid metabolism and inflammation. The results show a distinct interaction effect between maternal and nursery diets, implying that sow late gestation diet could be used to optimize piglet metabolism.

Effects of Chinese wolfberry and astragalus extracts on growth performance, pork quality, and unsaturated fatty acid metabolism regulation in Tibetan fragrant pigs

We studied the effects of wolfberry and astragalus extract on the growth performance, carcass traits, and meat quality of Tibetan fragrant pigs, and we want to explain the mechanism of the difference from the level of RNA Seq. Twelve healthy 120-day-old Tibetan fragrant pigs weighing 35 ± 2 kg were divided randomly into two groups, each with six pigs. The control group was fed a basal diet, and the wolfberry and astragalus extract (WAE) group was fed a basal diet +1‰ of WAE. The experimental period was 90 days. Compared with the control group, the growth performance of the WAE group was significantly improved (p < .05), pork marble score significantly improved (p < .05), vitamin E content significantly increased (p < .05), unsaturated fatty acid content significantly increased (p < .05). A total of 256 differentially expressed genes were obtained by transcriptome sequencing, among which 114 were up-regulated and 142 were down-regulated. GO analysis showed that the differentially expressed genes were related to biological functions, such as monounsaturated fatty acid biosynthesis, fatty acid metabolism, lipoprotein decomposition, and lipase activity. Pathway analysis showed that these differentially expressed genes were mainly involved in unsaturated fatty acid biosynthesis regulation, glycerin metabolism, and lipopolysaccharide regulation in fat. WAE improved Tibetan fragrant pigs growth performance. By intervening in key genes related to fatty acid metabolism, the unsaturated fatty acid contents in pork were regulated, which improved the nutritional value of the pork.

Transcriptomic Profiling of Skeletal Muscle Reveals Candidate Genes Influencing Muscle Growth and Associated Lipid Composition in Portuguese Local Pig Breeds

Simple Summary

Screening and interpretation of differentially expressed genes and associated biological pathways was conducted among experimental groups with divergent phenotypes providing valuable information about the metabolic events occurring and identification of candidate genes with major regulation roles. This comparative transcriptomic analysis includes the first RNA-seq analysis of the Longissimus lumborum muscle tissue from two Portuguese autochthonous pig breeds with different genetic backgrounds, Alentejano and Bísaro. Moreover, a complementary candidate gene approach was employed to analyse, by real time qPCR, the expression profile of relevant genes involved in lipid metabolism, and therefore with potential impacts on meat composition. This study contributes to explaining the biological basis of phenotypical differences occurring between breeds, particularly the ones related to meat quality traits that affect consumer interest.

Abstract

Gene expression is one of the main factors to influence meat quality by modulating fatty acid metabolism, composition, and deposition rates in muscle tissue. This study aimed to explore the transcriptomics of the Longissimus lumborum muscle in two local pig breeds with distinct genetic background using next-generation sequencing technology and Real-Time qPCR. RNA-seq yielded 49 differentially expressed genes between breeds, 34 overexpressed in the Alentejano (AL) and 15 in the Bísaro (BI) breed. Specific slow type myosin heavy chain components were associated with AL (MYH7) and BI (MYH3) pigs, while an overexpression of MAP3K14 in AL may be associated with their lower loin proportion, induced insulin resistance, and increased inflammatory response via NFkB activation. Overexpression of RUFY1 in AL pigs may explain the higher intramuscular (IMF) content via higher GLUT4 recruitment and consequently higher glucose uptake that can be stored as fat. Several candidate genes for lipid metabolism, excluded in the RNA-seq analysis due to low counts, such as ACLY, ADIPOQ, ELOVL6, LEP and ME1 were identified by qPCR as main gene factors defining the processes that influence meat composition and quality. These results agree with the fatter profile of the AL pig breed and adiponectin resistance can be postulated as responsible for the overexpression of MAP3K14′s coding product NIK, failing to restore insulin sensitivity.

Effects of Chinese wolfberry and Astragalus extract on the antioxidant capacity of Tibetan pig liver

The objective of this study is to determine the effect of Chinese wolfberry (Lycium barbarum) and Astragalus (Astragalus membranaceus) extract (WAE) on the antioxidant capacity of Tibetan pig liver, and discussed the regulatory effect of WAE on the liver antioxidant mechanism. Twelve healthy 120-day-old Tibetan black pigs (35±2 kg) were divided randomly into two groups. The WAE group was fed a basal diet supplemented with 1% WAE for 90 days. The control group was fed the same diet, but without the WAE. We found that liver superoxide dismutase 1 (SOD1) activity (P<0.05), total antioxidative capacity (T-AOC) (P<0.05), and catalase (CAT) activity (P<0.01) significantly increased in the WAE group compared with the control group; malondialdehyde (MDA) content decreased, but this was not significant (P >0.05). Transcriptome sequencing analysis detected 106 differentially expressed genes (DEGs) related to oxidative stress. GO enrichment analysis showed these DEGs were involved in the positive regulation of reactive oxygen metabolism and biosynthesis, process regulation, and regulation of the oxidative stress response. KEGG Pathway enrichment analysis showed they were enriched in the PI3K-Akt, AMPK, Rap1, and peroxisome signaling pathways. The expression levels of key peroxisome biosynthesis genes (e.g., PEX3 and PEX11B) and key antioxidant genes (e.g., CAT and SOD1) were significantly higher in the WAE group than in the control group. The PRDX1 and PRDX5 content also was significantly higher in the WAE group. This study showed that the WAE regulated the antioxidant and anti-stress ability of Tibetan pig liver through a "peroxisome antioxidant-oxidant stress" signaling pathway.

Expression analysis of porcine miR-33a/b in liver, adipose tissue and muscle and its potential role in fatty acid metabolism

mir-33a and mir-33b are co-transcribed with the SREBF2 and SREBF1 transcription factors, respectively. The main role of SREBF1 is the regulation of genes involved in fatty acid metabolism, while SREBF2 regulates genes participating in cholesterol biosynthesis and uptake. Our objective was to study the expression of both miR-33a and miR-33b, together with their host SREBF genes, in liver, adipose tissue and muscle to better understand the role of miR-33a/b in the lipid metabolism of pigs. In our study, the expression of miR-33a, miR-33b and SREBF2 in liver, adipose tissue, and muscle was studied in 42 BC1_LD (25% Iberian x 75% Landrace backcross) pigs by RT-qPCR. In addition, the expression of in-silico predicted target genes and fatty acid composition traits were correlated with the miR-33a/b expression. We observed different tissue expression patterns for both miRNAs. In adipose tissue and muscle a high correlation between miR-33a and miR-33b expression was found, whereas a lower correlation was observed in liver. The expression analysis of in-silico predicted target-lipid related genes showed negative correlations between miR-33b and CPT1A expression in liver. Conversely, positive correlations between miR-33a and PPARGC1A and USF1 gene expression in liver were observed. Lastly, positive and negative correlations between miR-33a/b expression and saturated fatty acid (SFA) and polyunsaturated fatty acid (PUFA) content, respectively, were identified. Overall, our results suggested that both miRNAs are differentially regulated and have distinct functions in liver, in contrast to muscle and adipose tissue. Furthermore, the correlations between miR-33a/b expression both with the expression of in-silico predicted target-lipid related genes and with fatty acid composition, opens new avenues to explore the role of miR33a/b in the regulation of lipid metabolism.

Profiling and Functional Analysis of Long Noncoding RNAs and mRNAs during Porcine Skeletal Muscle Development

Gene transcripts or mRNAs and long noncoding RNAs (lncRNAs) are differentially expressed during porcine skeletal muscle development. However, only a few studies have been conducted on skeletal muscle transcriptome in pigs based on timepoints according to the growth curve for porcine. Here, we investigated gene expression in Qingyu pigs at three different growth stages: the inflection point with the maximum growth rate (MGI), the inflection point of the gradually increasing stage to the rapidly increasing stage (GRI), and the inflection point of the rapidly increasing stage to the slowly increasing stage (RSI). Subsequently, we explored gene expression profiles during muscle development at the MGI, GRI and RSI stages by Ribo-Zero RNA sequencing. Qingyu pigs reached the MGI, GRI and RSI stages at 156.40, 23.82 and 288.97 days of age with 51.73, 3.14 and 107.03 kg body weight, respectively. A total of 14,530 mRNAs and 11,970 lncRNAs were identified at the three stages, and 645, 323 differentially expressed genes (DEGs) and 696, 760 differentially expressed lncRNAs (DELs) were identified in the GRI vs. MGI, and RSI vs. MGI, comparisons. Functional enrichment analysis revealed that genes involved in immune system development and energy metabolism (mainly relate to amino acid, carbohydrate and lipid) were enriched at the GRI and MGI stages, respectively, whereas genes involved in lipid metabolism were enriched at the RSI stage. We further characterized G1430, an abundant lncRNA. The full-length sequence (316 nt) of lncRNA G1430 was determined by rapid amplification of cDNA ends (RACE). Subcellular distribution analysis by quantitative real-time PCR (qRT-PCR) revealed that G1430 is a cytoplasmic lncRNA. Binding site prediction and dual luciferase assay showed that lncRNA G1430 directly binds to microRNA 133a (miR-133a). Our findings provide the basis for further investigation of the regulatory mechanisms and molecular genetics of muscle development in pigs.

The effect of Acot2 overexpression or downregulation on the preadipocyte differentiation in Chinese Red Steppe cattle

The quality and nutritional value of beef is closely linked to its content of intramuscular fat (IMF). The differentiation of preadipocytes and the deposition of lipid droplets in the adipocytes are the key to regulate the IMF content. The differentiation of adipocytes is regulated by a series of transcription factors and genes. Acyl-CoA thioesterase 2 (Acot2) hydrolyzes the acyl-coenzyme A (CoA) into free fatty acids and CoA and has the potential to maintain the free fatty acids and acyl CoA at the cellular level. In this work, we detected the expression of the Acot2 gene during the adipocyte differentiation in Chinese Red Steppe cattle, and then interfered and overexpressed the Acot2 gene in the preadipocytes to explore its regulatory role in the adipocyte differentiation. The results showed that the expression and regulation of Acot2 mainly occurred at the later stage of the adipocyte differentiation. The interference with the Acot2 gene significantly inhibited the lipid droplets accumulation and triglyceride content, while its overexpression significantly promoted both of them. The results of this study show that the Acot2 gene is a positive regulator of the adipocyte differentiation and may become a new target to improve the beef quality.

Profiling of skeletal muscle tissue for long non-coding RNAs related to muscle metabolism in the QingYu pig at the growth inflection point

OBJECTIVE

Investigation of muscle growth at different developmental stages is an appropriate strategy for studying the mechanisms underlying muscle development and differences in phenotypes. In particular, the muscle development mechanisms and the difference between the fastest and slowest growth rates.

METHODS

In this study, we used a growth curve model to fit the growth inflection point (IP) of QingYu pigs and compared differences in the long non-coding RNA (lncRNA) transcriptome of muscle both at the growth IP and plateau phase (PP).

RESULTS

The growth curve of the QingYu pig had a good fit (R2 = 0.974) relative to a typical S-curve and reached the IP at day 177.96. At the PP, marbling, intramuscular fat, and monounsaturated fatty acids had increased significantly and the percentage of lean muscle and polyunsaturated fatty acids had decreased. A total of 1,199 mRNAs and 62 lncRNAs were differentially expressed at the IP compared with the PP. Additional to gene ontology and Kyoto encyclopedia of genes and genomes pathway analyses, these differentially expressed protein coding genes were principally related to muscle growth and lipid metabolism.

CONCLUSION

Our results suggest that the identified differentially expressed lncRNAs, could play roles in muscle growth, fat deposition and regulation of fatty acid composition at the IP and PP.

Revisiting the Pig IGHC Gene Locus in Different Breeds Uncovers Nine Distinct IGHG Genes

IgG subclass diversification is common in placental mammals. It has been well documented in humans and mice that different IgG subclasses, with diversified functions, synergistically regulate humoral immunity. However, our knowledge on the genomic and functional diversification of IgG subclasses in the pig, a mammalian species with high agricultural and biomedical importance, is incomplete. Using bacterial artificial chromosome sequencing and newly assembled genomes generated by the PacBio sequencing approach, we characterized and mapped the IgH C region gene locus in three indigenous Chinese breeds (Erhualian, Xiang, and Luchuan) and compared them to that of Duroc. Our data revealed that IGHG genes in Chinese pigs differ from the Duroc, whereas the IGHM, IGHD, IGHA, and IGHE genes were all single copy and highly conserved in the pig breeds examined. Most striking were differences in numbers of IGHG genes: there are seven genes in Erhualian pigs, six in the Duroc, but only five in Xiang pigs. Phylogenetic analysis suggested that all reported porcine IGHG genes could be classified into nine subclasses: IGHG1, IGHG2a, IGHG2b, IGHG2c, IGHG3, IGHG4, IGHG5a, IGHG5b, and IGHG5c. Using sequence information, we developed a mouse mAb specific for IgG3. This study offers a starting point to investigate the structure-function relationship of IgG subclasses in pigs.

Transcriptome analysis of differential gene expression in the longissimus dorsi muscle from Debao and landrace pigs based on RNA-sequencing

RNA-seq analysis was used to identify differentially expressed genes (DEGs) at the genetic level in the longissimus dorsi muscle from two pigs to investigate the genetic mechanisms underlying the difference in meat quality between Debao pigs and Landrace pigs. Then, these DEGs underwent functional annotation, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, and protein–protein interaction (PPI) analyses. Finally, the expression levels of specific DEGs were assessed using qRT-PCR. The reference genome showed gene dosage detection of all samples which showed that the total reference genome comprised 22342 coding genes, including 14743 known and 190 unknown genes. For detection of the Debao pig genome, we obtained 14168 genes, including 13994 known and 174 unknown genes. For detection of the Landrace pig genome, we obtained 14404 genes, including 14223 known and 181 unknown genes. GO analysis and KEGG signaling pathway analysis show that DEGs are significantly related to metabolic regulation, amino acid metabolism, muscular tissue, muscle structure development etc. We identified key genes in these processes, such as FOS, EGR2, and IL6, by PPI network analysis. qRT-PCR confirmed the differential expression of six selected DEGs in both pig breeds. In conclusion, the present study revealed key genes and related signaling pathways that influence the difference in pork quality between these breeds and could provide a theoretical basis for improving pork quality in future genetic thremmatology.

Transcriptome analysis reveals the genetic basis of skeletal muscle glycolytic potential based on a pig model

Glycolytic potential (GP) calculated based on glucose, glycogen, glucose-6-phosphate, and lactate contents is a critical factor for multiple meat quality characteristics. However, the genetic basis of glycolytic metabolism is still unclear. In this study, we constructed six RNA-Seq libraries using longissimus dorsi (LD) muscles from pigs divergent for GP phenotypic values and generated the whole genome-wide gene expression profiles. Furthermore, we identified 25,880 known and 220 novel genes from these skeletal muscle libraries, and 222 differentially expressed genes (DEGs) between the higher and lower GP groups. Notably, we found that the Lactate dehydrogenase B (LDHB) and Fructose-2, 6-biphosphatase 3 (PFKFB3) expression levels were higher in the higher GP group than the lower GP group, and positively correlated with GP and lactic acid (LA), and reversely correlated with pH value at 45 min postmortem (pH_45min). Besides, LDHB and PFKFB3 expression were positively correlated with drip loss measured at 48 h postmortem (DL_48h) and drip loss measured at 24 h postmortem (DL_24h). Collectively, we identified a serial of DEGs as the potential key candidate genes affecting GP and found that LDHB and PFKFB3 are closely related to GP and GP-related traits. Our results lay a solid basis for in-depth studies of the regulatory mechanisms on GP and GP-related traits in pigs.

Comprehensive analysis of the differential expression profile of microRNAs in rats with spinal cord injury treated by electroacupuncture

Abnormal microRNA (miRNA) expression has been implicated in spinal cord injury (SCI), but the underlying mechanisms are poorly understood. To observe the effect of electroacupuncture (EA) on miRNA expression profiles in SCI rats and investigate the potential mechanisms involved in this process, Sprague-Dawley rats were divided into sham, SCI and SCI+EA groups (n=6 each). Basso, Beattie and Bresnahan (BBB) scoring and hematoxylin-eosin staining of cortical tissues were used to evaluate spinal cord recovery with EA treatment 21 days post-surgery across the three groups. To investigate miRNA expression profiles, 6 Sprague-Dawley rats were randomly divided into SCI and SCI+EA groups (n=3 in each group) and examined using next-generation sequencing. Integrated miRNA-mRNA-pathway network analysis was performed to elucidate the interaction network of the candidate miRNAs, their target genes and the involved pathways. Behavioral scores suggested that hindlimb motor functions improved with EA treatments. Apoptotic indices were lower in the SCI+EA group compared with the SCI group. It was also observed that 168 miRNAs were differentially expressed between the SCI and SCI+EA groups, with 29 upregulated and 139 downregulated miRNAs in the SCI+EA group. Changes in miRNA expression are involved in SCI physiopathology, including inflammation and apoptosis. Reverse transcription-quantitative PCR measurement of the five candidate miRNAs, namely rno-miR-219a-5p, rno-miR-486, rno-miR-136-5p, rno-miR-128-3p, and rno-miR-7b, was consistent with RNA sequencing data. Integrated miRNA-mRNA-pathway analysis suggested that the MAPK, Wnt and NF-κB signaling pathways were involved in EA-mediated recovery from SCI. The present study evaluated the miRNA expression profiles involved in EA-treated SCI rats and demonstrated the potential mechanism and functional role of miRNAs in SCI in rats.

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.

METTL14 promotes the migration and invasion of breast cancer cells by modulating N6‑methyladenosine and hsa‑miR‑146a‑5p expression

Breast cancer (BC) is the most frequently diagnosed cancer and the leading cause of cancer-related death among women worldwide. Evidence indicates that posttranscriptional N6-methyladenosine (m6A) modification modulates BC development. In the present study, we assessed BC and normal tissues to investigate this connection. RNA m6A levels were determined by methylation quantification assay. The effects of methyltransferase-like 14 (METTL14) gain-of-expression or co-transfection with an m6A inhibitor on cell migration and invasion abilities were determined by Transwell assays. The levels of differentially expressed (DE) miRNAs were verified by real-time quantitative PCR (RT-qPCR). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses (KEGG) were performed to analyze potential function of target genes of the DE miRNAs. The effects of candidate miRNAs modulated by METTL14 on cell migration and invasion abilities were confirmed by Transwell assays. We demonstrated that m6A methyltransferase METTL14 was significantly upregulated in BC tissues compared with normal tissues. METTL14 gain- and loss-of-expression regulated m6A levels in MCF-7 and MDA-MB-231 cells. The cell function assays revealed that METTL14 overexpression enhanced the migration and invasion capacities of BC cells. Moreover, treatment with the m6A inhibitor suppressed this enhanced cell migration and invasion. Additionally, aberrant expression of METTL14 reshaped the miRNA profile in BC cell lines. The remodeled DE miRNA/mRNA network was found to be most enriched in cancer pathways, and DE miRNAs were enriched in cell adhesion terms. hsa-miR-146a-5p modulated by METTL14 promoted cell migration and invasion. METTL14 modulates m6A modification and hsa-miR-146a-5p expression, thereby affecting the migration and invasion of breast cancer cells.

Identification of Metabonomics Changes in Longissimus Dorsi Muscle of Finishing Pigs Following Heat Stress through LC-MS/MS-Based Metabonomics Method

Simple Summary

Limited research exists on muscle metabolomics of finishing pigs under heat stress. In this study, nine different metabolites in the longissimus dorsi (LD) muscle of finishing pigs under heat stress were screened and identified. Through quantitative verification, it was concluded that the content of L-carnitine in the LD muscles of the finishing pigs could be significantly decreased due to heat stress, which might be a biomarker for monitoring the animal health status and muscle quality under heat stress.

Abstract

Heat stress (HS) negatively affects meat quality by affecting material and energy metabolism, and exploring the mechanism underlying the muscle response to chronic HS in finishing pigs is important for the global pork industry. This study investigated changes in the metabolic profiles of the longissimus dorsi (LD) muscle of finishing pigs under high temperature using ultra-performance liquid chromatography coupled with electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC-ESI-QTOF-MS) and multivariate data analysis (MDA). Castrated male DLY pigs (Duroc × Landrance × Yorkshire pigs, n = 24) from 8 litters were divided into three treatment groups: constant optimal ambient temperature at 22 °C and ad libitum feeding (CR, n = 8); constant high ambient temperature at 30 °C and ad libitum feeding (HS, n = 8); and constant optimal ambient temperature 22 °C and pair-feeding to the control pigs (PF, n = 8). The metabolic profile data from LD muscle samples were analyzed by MDA and external search engines. Nine differential metabolites (L-carnosine, acetylcholine, inosinic acid, L-carnitine, L-anserine, L-α-glycerylphosphorylcholine, acetylcarnitine, thiamine triphosphate, and adenosine thiamine diphosphate) were involved in antioxidant function, lipid metabolism, and cell signal transduction, which may decrease post mortem meat quality and play important roles in anti-HS. Four metabolites (L-carnosine, acetylcholine, inosinic acid, and L-carnitine) were verified, and it was indicated that the muscle L-carnitine content was significantly lower in HS than in CR (p < 0.01). The results show that constant HS affects the metabolites in the LD muscle and leads to coordinated changes in the endogenous antioxidant defense and meat quality of finishing pigs. These metabonomics results provide a basis for researching nutritional strategies to reduce the negative effects of heat stress on livestock and present new insights for further research.

Transcriptome profiling of the liver among the prenatal and postnatal stages in chickens

The liver is an important organ that has pivotal functions in the synthesis of several vital proteins, the metabolism of various biologically useful materials, the detoxification of toxic substances, and immune defense. Most liver functions are not mature at a young age and many changes happen during postnatal liver development, which lead to differential functions of the liver at different developmental stages. However, the transcriptome details of what changes occur in the liver after birth and the molecular mechanisms for the regulation of the developmental process are not clearly known in chickens. Here, we used RNA-sequencing to analyze the transcriptome of chicken liver from the prenatal (at an embryonic day of 13) to the postnatal stages (at 5 wk and 42 wk of age). A total of approximately 161.17 Gb of raw data were obtained, with 4,127 putative and 539 differentially expressed lncRNAs, and with 13,949 putative and 6,370 differentially expressed mRNAs. Coexpression of lncRNAs-mRNAs in hepatic transcriptome analysis showed that the liver plays important roles in providing energy for organisms through the mitochondrial respiratory chain in chickens, meanwhile, acting as a crucial part of antioxidant stress. The developmental transcriptome date revealed that antioxidant defenses are likely to act on chicken embryo development and that significant functional changes during postnatal liver development are associated with the liver maturation of chickens. These results provide a timeline for the functional transcriptome transition from the prenatal to adult stages in chickens and will be helpful to reveal the underlying molecular mechanisms of liver development.

Transcriptome analysis of miRNA and mRNA in the livers of pigs with highly diverged backfat thickness

Fat deposition is very important in pig production, and its mechanism is not clearly understood. MicroRNAs (miRNAs) play critical roles in fat deposition and energy metabolism. In the current study, we investigated the mRNA and miRNA transcriptome in the livers of Landrace pigs with extreme backfat thickness to explore miRNA-mRNA regulatory networks related to lipid deposition and metabolism. A comparative analysis of liver mRNA and miRNA transcriptomes from pigs (four pigs per group) with extreme backfat thickness was performed. We identified differentially expressed genes from RNA-seq data using a Cufflinks pipeline. Seventy-one differentially expressed genes (DEGs), including twenty-eight well annotated on the porcine reference genome genes, were found. The upregulation genes in pigs with higher backfat thickness were mainly involved in fatty acid synthesis, and included fatty acid synthase (FASN), glucokinase (GCK), phosphoglycerate dehydrogenase (PHGDH), and apolipoprotein A4 (APOA4). Cytochrome P450, family 2, subfamily J, polypeptide 34 (CYP2J34) was lower expressed in pigs with high backfat thickness, and is involved in the oxidation of arachidonic acid. Moreover, 13 differentially expressed miRNAs were identified. Seven miRNAs were associated with fatty acid synthesis, lipid metabolism, and adipogenic differentiation. Based on comprehensive analysis of the transcriptome of both mRNAs and miRNAs, an important regulatory network, in which six DEGs could be regulated by differentially expressed miRNAs, was established for fat deposition. The negative correlate in the regulatory network including, miR-545-5p and GRAMD3, miR-338 and FASN, and miR-127, miR-146b, miR-34c, miR-144 and THBS1 indicate that direct suppressive regulation may be involved in lipid deposition and energy metabolism. Based on liver mRNA and miRNA transcriptomes from pigs with extreme backfat thickness, we identified 28 differentially expressed genes and 13 differentially expressed miRNAs, and established an important miRNA-mRNA regulatory network. This study provides new insights into the molecular mechanisms that determine fat deposition in pigs.

Comparative transcriptomic analysis of hypothalamus-pituitary-liver axis in bighead carp (Hypophthalmichthys nobilis) with differential growth rate

Background

Growth rate is one of the most important features for aquaculture species and deciphering its regulation mechanism has great significance both in genetics and in economics. Hypothalamus-pituitary growth axis (HP growth axis) or neuro-endocrine axis plays a vital role in growth regulation in different aquaculture animals.

Results

In this study, the HP and liver transcriptomes of two female groups (H and L) with phenotypically extreme growth rate were sequenced using RNA-Seq. A total of 30,524 and 22,341 genes were found expressed in the two tissues, respectively. The average expression levels for the two tissues were almost the same, but the median differed significantly. A differential expression analysis between H and L groups identified 173 and 204 differentially expressed genes (DEGs) in HP and liver tissue, respectively. Pathway analysis revealed that DEGs in HP tissue were enriched in regulation of cell proliferation and angiogenesis while in liver tissue these genes were overrepresented in sterol biosynthesis and transportation. Genomic overlapping analyses found that 4 and 5 DEGs were within growth-related QTL in HP and liver tissue respectively. A deeper analysis of these 9 genes revealed 3 genes were functionally linked to the trait of interest. The expression of 2075 lncRNAs in HP tissue and 1490 in liver tissue were also detected, and some of lncRNAs were highly expressed in the two tissues.

Conclusions

Above all, the results of the present study greatly contributed to the knowledge of the regulation of growth and then assisted the design of new selection strategies for bighead carp with improved growth-related traits.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5691-4) contains supplementary material, which is available to authorized users.

Integrated Analysis of miRNA and mRNA Expression Profiles Reveals Functional miRNA-Targets in Development Testes of Small Tail Han Sheep

Small Tail Han Sheep is a highly valued local breed in China because of their precocity, perennial estrus, and high fecundity. The average annual lambing rate of ewes is as high as 180–270%, the semen of ram has characteristics of high yield, high density, and good motility. To reveal the key miRNAs and miRNA-targets underlying testis development and spermatogenesis in male Small Tail Han Sheep, integrated analysis of miRNA and mRNA expression profiles in 2-, 6-, and 12-month-old testes was performed by RNA-seq technology and bioinformatics methods. The results showed that total of 153 known sheep miRNAs and 2712 novel miRNAs were obtained in 2-,6 - and 12-month-old Small Tail Han Sheep testes; 5, 1, and 4 differentially expressed (DE) known sheep miRNAs, and 132, 105, and 24 DE novel miRNAs were identified in 2- vs. 6-, 6- vs. 12-, and 2- vs. 12-month-old testes, respectively. We combined miRNA results of this study and the mRNA results obtained in our previous study to predict the target mRNAs of DE known sheep miRNAs; 131, 10, and 15 target mRNAs of DE known sheep miRNAs and 76, 1, and 11 DE miRNA–targets were identified in the three groups, respectively. GO and KEGG analyses showed that: in 2- vs. 6-month-olds, the target genes of DE known sheep miRNAs were involved in 100 biological processes and 11 signaling pathways; in 6- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 4 biological processes; and in 2- vs. 12-month-olds, the target genes of DE known sheep miRNAs were involved in 17 biological processes and 4 signaling pathways. Three miR–target regulatory networks were constructed based on these DE miRNA–targets. The key miRNA-Targets involved in testis development and spermatogenesis were screened. 6 known sheep miRNAs and 6 novel miRNAs were selected to validate the accuracy of miRNA sequencing data by qRT-PCR. The binding sites of oar-miR-379-5p with WNT8A was validated by a dual luciferase reporter gene detection system.

Systems genetics identifies a macrophage cholesterol network associated with physiological wound healing

Among other cells, macrophages regulate the inflammatory and reparative phases during wound healing but genetic determinants and detailed molecular pathways that modulate these processes are not fully elucidated. Here, we took advantage of normal variation in wound healing in 1,378 genetically outbred mice, and carried out macrophage RNA-sequencing profiling of mice with extreme wound healing phenotypes (i.e., slow and fast healers, n = 146 in total). The resulting macrophage coexpression networks were genetically mapped and led to the identification of a unique module under strong trans-acting genetic control by the Runx2 locus. This macrophage-mediated healing network was specifically enriched for cholesterol and fatty acid biosynthetic processes. Pharmacological blockage of fatty acid synthesis with cerulenin resulted in delayed wound healing in vivo, and increased macrophage infiltration in the wounded skin, suggesting the persistence of an unresolved inflammation. We show how naturally occurring sequence variation controls transcriptional networks in macrophages, which in turn regulate specific metabolic pathways that could be targeted in wound healing.

Transcriptome Analyses Reveal Adult Metabolic Syndrome With Intrauterine Growth Restriction in Pig Models

Epidemiological data have indicated that intrauterine growth retardation (IUGR) is a risk factor for the adult metabolic syndrome in pigs. However, the causative genetic mechanism leading to the phenotype in adulthood has not been well characterized. In the present study, both normal and IUGR adult pigs were used as models to survey the differences in global gene expression in livers through transcriptome sequencing. The transcriptome libraries generated 104.54 gb of data. In normal and IUGR pigs, 16,948 and 17,078 genes were expressed, respectively. A total of 1,322 differentially expressed genes (DEGs) were identified. Enrichment analysis of the DEGs revealed that the top overrepresented gene ontology (GO) terms and pathways were related to oxidoreductase activity, ATPase activity, amino catabolic process, glucose metabolism, and insulin signaling pathway. The increased gluconeogenesis (GNG) and decreased glycogen synthesis in the liver contributed to the glucose intolerance observed in IUGR. The reduced expression of insulin signaling genes (such as PI3K and AKT) indicated an elevated risk of diabetes in adulthood. Together, these findings provide a comprehensive understanding of the molecular mechanisms of adult IUGR pigs and valuable information for future studies of therapeutic intervention in IUGR metabolic syndrome.

Identification of Long Non-Coding RNAs Related to Skeletal Muscle Development in Two Rabbit Breeds with Different Growth Rate

Skeletal muscle development plays an important role in muscle quality and yield, which decides the economic value of livestock. Long non-coding RNAs (lncRNAs) have been reported to be associated with skeletal muscle development. However, little is revealed about the function of lncRNAs in rabbits' muscle development. LncRNAs and mRNAs in two rabbit breeds (ZIKA rabbits (ZKR) and Qixin rabbits (QXR)) with different growth rates at three developmental stages (0 day, 35 days, and 84 days after birth) were researched by transcriptome sequencing. Differentially expressed lncRNAs and mRNAs were identified for two rabbit breeds at the same stages by DESeq package. Co-expression correlation analysis of differentially expressed lncRNAs and mRNAs were performed to construct lncRNA⁻mRNA pairs. To explore the function of lncRNAs, Gene Ontology (GO) analysis of co-expression mRNAs in lncRNA⁻mRNA pairs were performed. In three comparisons, there were 128, 109, and 115 differentially expressed lncRNAs, respectively. LncRNAs TCONS_00013557 and XR_518424.2 differentially expressed in the two rabbit breeds might play important roles in skeletal muscle development, for their co-expressed mRNAs were significantly enriched in skeletal muscle development related GO terms. This study provides potentially functional lncRNAs in skeletal muscle development of two rabbit breeds and might be beneficial to the production of rabbits.