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

RNA-seq analysis of small intestine transcriptional changes induced by starvation stress in piglets

Piglets may experience a variety of stress injuries, but the molecular regulatory mechanisms underlying these injuries are not well understood. In this study, we analysed the ileum of Large White (LW) and Mashen (MS) piglets at different times of starvation using chemical staining and transcriptome analysis. The intestinal barrier of piglets was damaged after starvation stress, but the intestinal antistress ability of MS piglets was stronger than LW piglets. A total of 8021 differentially expressed genes (DEGs) were identified in two breeds. Interestingly, the immune capacity (CHUK, TLR3) of MS piglets increased significantly after short-term starvation stress, while energy metabolism (NAGS, PLA2G12B, AGCG8) was predominant in LW piglets. After long-term starvation stress, the level of energy metabolism (PLIN5, PLA2G12B) was significantly increased in MS piglets. The expression of immune (HLA-DQB1, IGHG4, COL3A1, CD28, LAT) and disease (HSPA1B, MINPPI, ADH1C, GAL3ST1) related genes were significantly increased in two breeds of piglets. These results suggest that short-term stress mainly enhances immunity and energy metabolism in piglets, while long-term starvation produces greater stress on piglets, making it difficult for them to compensate for the damage to their bodies through self-regulation. This information can help improve the stress resistance of piglets through molecular breeding.

Whole-genome resequencing of native and imported dairy goat identifies genes associated with productivity and immunity

Understanding the differences in genetic variation between local Chinese dairy goat breeds and imported breeds can help germplasm innovation and molecular breeding. However, the research is limited in this area. In this study, whole-genome resequencing data from 134 individuals of both local and imported dairy goat breeds were analyzed, and their differences in genomic genetic variation, genetic diversity, and population structure were subsequently identified. We also screened candidate genes associated with important traits of dairy goats such as milk production (STK3, GHR, PRELID3B), reproduction (ATP5E), growth and development (CTSZ, GHR), and immune function (CTSZ, NELFCD). Furthermore, we examined allele frequency distributions for the genes of interest and found significant differences between the two populations. This study provides valuable resources for the study of genetic diversity in dairy goats and lays the foundation for the selective breeding of dairy goats in the future.

Single-nucleus RNA sequencing and lipidomics reveal characteristics of transcriptional and lipid composition in porcine longissimus dorsi muscle

BACKGROUND

Global per capita meat consumption continues to rise, especially pork. Meat quality is influenced by the content of intramuscular fat (IMF) as a key factor. The longissimus dorsi muscle of Dahe pigs (DHM, IMF: 7.98% ± 1.96%) and Dahe black pigs (DHBM, IMF: 3.30% ± 0.64%) was studied to explore cellular heterogeneity and differentially expressed genes (DEGs) associated with IMF deposition using single-nucleus RNA sequencing (snRNA-seq). The lipid composition was then analyzed using non-targeted lipidomics.

RESULTS

A total of seven cell subpopulations were identified, including myocytes, fibroblast/fibro/adipogenic progenitors (FAPs), satellite cells, endothelial cells, macrophages, pericytes, and adipocytes. Among them, FAPs and adipocytes were more focused because they could be associated with lipid deposition. 1623 DEGs in the FAPs subpopulation of DHBM were up-regulated compared with DHM, while 1535 were down-regulated. These DEGs enriched in the glycolysis/gluconeogenesis pathway. 109 DEGs were up-regulated and 806 were down-regulated in the adipocyte subpopulation of DHBM compared with DHM, which were mainly enriched in the PPAR signaling pathway and fatty acid (FA) biosynthesis. The expression level of PPARG, ABP4, LEP, and ACSL1 genes in DHM was higher than that in DHBM. Lipidomics reveals porcine lipid composition characteristics of muscle tissue. A total of 41 lipid classes and 2699 lipid species were identified in DHM and DHBM groups. The top ten relative peak areas of lipid classes in DHM and DHBM were triglyceride (TG), phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylserine (PS), diglyceride (DG), cardiolipin (CL), ceramides (Cer), Simple Glc series (Hex1Cer), sphingomyelin (phSM), and phosphatidylinositol (PI). The relative peak areas of 35 lipid species in DHM were lower than DHBM, and 28 lipid species that were higher. There was a significant increase in the TG fatty acyl chains C6:0, C17:0, and C11:4, and a significant decrease in C16:0, C18:1, C18:2, and C22:4 in DHBM (p < 0.05).

CONCLUSIONS

C16:0 FA may downregulate the expression level of PPARG gene, which leads to the downregulation of fat metabolism-related genes such as ACSL, PLIN2, and FABP4 in DHBM compared with DHM. This may be the reason that the lipid deposition ability of Dahe pigs is stronger than that of Dahe black pigs, which need further investigation.

Skeletal muscle transcriptomics of sheep acclimated to cold desert and tropical regions identifies genes and pathways accentuating their diversity

The current study attempts to investigate the differences in gene expression in longissimus thoracis muscles between sheep breeds acclimated to diverse environments. Changthangi sheep inhabits the cold arid plateau of Ladakh, at an altitude above 3000 m with prevalence of rarefied atmosphere. Muzzafarnagri sheep, on the other hand is found in the sub-tropical hot and humid plains at an altitude of about 250 m. Comparative transcriptomics was used to provide a molecular perspective of the differential adaptation of the two breeds. RNA sequencing data was generated from four biological replicates of the longissimus thoracis muscles from both breeds. The common genes expressed in both breeds were involved in muscle contraction and muscle fibre organization. The most significant pathways enriched in Changthangi muscles were glycogen metabolism, reduction of cytosolic Ca++ levels and NFE2L2 regulating anti-oxidant, while those in Muzzafarnagri were extracellular matrix organization and collagen formation. The hub genes identified in Changthangi were involved in hematopoiesis and HIF signaling pathway, suggesting the molecular acclimatization of Changthangi to the high altitude cold desert of Ladakh. The nodal genes discovered in Muzzafarnagri sheep were associated with the extracellular matrix which accentuates its significance in the development, growth and repair of muscles. The observed transcriptomic differences underscore the morphological and adaptive disparity between the two breeds. The candidate genes and pathways identified in this study will form the basis for future research on adaptation to high altitude and body size in small ruminants.

Genomic evidence for human-mediated introgressive hybridization and selection in the developed breed

BACKGROUND

The pig (Sus Scrofa) is one of the oldest domesticated livestock species that has undergone extensive improvement through modern breeding. European breeds have advantages in lean meat development and highly-productive body type, whereas Asian breeds possess extraordinary fat deposition and reproductive performance. Consequently, Eurasian breeds have been extensively used to develop modern commercial breeds for fast-growing and high prolificacy. However, limited by the sequencing technology, the genome architecture of some nascent developed breeds and the human-mediated impact on their genomes are still unknown.

RESULTS

Through whole-genome analysis of 178 individuals from an Asian locally developed pig breed, Beijing Black pig, and its two ancestors from two different continents, we found the pervasive inconsistent gene trees and species trees across the genome of Beijing Black pig, which suggests its introgressive hybrid origin. Interestingly, we discovered that this developed breed has more genetic relationships with European pigs and an unexpected introgression from Asian pigs to this breed, which indicated that human-mediated introgression could form the porcine genome architecture in a completely different type compared to native introgression. We identified 554 genomic regions occupied 63.30 Mb with signals of introgression from the Asian ancestry to Beijing Black pig, and the genes in these regions enriched in pathways associated with meat quality, fertility, and disease-resistant. Additionally, a proportion of 7.77% of genomic regions were recognized as regions that have been under selection. Moreover, combined with the results of a genome-wide association study for meat quality traits in the 1537 Beijing Black pig population, two important candidate genes related to meat quality traits were identified. DNAJC6 is related to intramuscular fat content and fat deposition, and RUFY4 is related to meat pH and tenderness.

CONCLUSIONS

Our research provides insight for analyzing the origins of nascent developed breeds and genome-wide selection remaining in the developed breeds mediated by humans during modern breeding.

Dietary Betaine Improves Glucose Metabolism in Obese Mice

BACKGROUND

Obesity caused by the overconsumption of energy-dense foods high in fat and sugar has contributed to the growing prevalence of type 2 diabetes. Betaine, found in food or supplements, has been found to lower blood glucose concentrations, but its exact mechanism of action is not well understood.

OBJECTIVES

A comprehensive evaluation of the potential mechanisms by which betaine supplementation improves glucose metabolism.

METHODS

Hyperglycemic mice were fed betaine to measure the indexes of glucose metabolism in the liver and muscle. To explore the mechanism behind the regulation of betaine on glucose metabolism, Ribonucleic Acid-Seq was used to analyze the livers of the mice. In vitro, HepG2 and C2C12 cells were treated with betaine to more comprehensively evaluate the effect of betaine on glucose metabolism.

RESULTS

Betaine was added to the drinking water of high-fat diet-induced mice, and it was found to reduce blood glucose concentrations and liver triglyceride concentrations without affecting body weight, confirming its hypoglycemic effect. To investigate the specific mechanism underlying its hypoglycemic effect, protein-protein interaction enrichment analysis of the liver revealed key nodes associated with glucose metabolism, including cytochrome P450 family activity, insulin sensitivity, glucose homeostasis, and triglyceride concentrations. The Kyoto Encyclopedia of Genes and Genomes and gene ontogeny enrichment analyses showed significant enrichment of the Notch signaling pathway. These results provided bioinformatic evidence for specific pathways through which betaine regulates glucose metabolism. Key enzyme activities involved in glucose uptake, glycogen synthesis, and glycogenolysis pathways of the liver and muscle were measured, and improvements were observed in these pathways.

CONCLUSIONS

This study provides new insight into the mechanisms by which betaine improves glucose metabolism in the liver and muscle and supports its potential as a drug for the treatment of metabolic disorders related to glucose.

Transcriptome profiling of longissimus dorsi during different prenatal stages to identify genes involved in intramuscular fat deposition in lean and obese pig breeds

BACKGROUND

There was significant difference in muscle development between fat-type and lean-type pig breeds.

METHODS AND RESULTS

In current study, transcriptome analysis and bioinformatics analysis were used to compare the difference in longissimus dorsi (LD) muscle at three time-points (38 days post coitus (dpc), 58 dpc, and 78 dpc ) between Huainan (HN) and Large white (LW) pig breeds. A total of 24500 transcripts were obtained in 18 samples, and 2319, 2799, and 3713 differently expressed genes (DEGs) were identified between these two breeds at 38 dpc, 58 dpc, and 78 dpc, respectively. And the number and foldchange of DEGs were increased, the alternative splice also increased. The cluster analysis of DEGs indicated the embryonic development progress of LD muscle between these two breeds was different. There were 539 shared DEGs between HN and LW at three stages, and the top-shared DEGs were associated with muscle development and lipid deposition, such as KLF4, NR4A1, HSP70, ZBTB16 and so on.

CONCLUSIONS

The results showed DEGs between Huainan (HN) and Large white (LW) pig breeds, and contributed to the understanding the muscle development difference between HN and LW, and provided basic materials for improvement of meat quality.

Integrated analysis of lncRNA, miRNA and mRNA expression profiles reveals regulatory pathways associated with pig testis function

Long noncoding RNA (lncRNA) and microRNA (miRNA) are known to play pivotal roles in mammalian testicular function and spermatogenesis. However, their impact on porcine male reproduction has yet to be well unraveled. Here, we sequenced and identified lncRNA and miRNA expressed in the testes of Chinese indigenous Banna mini-pig inbred line (BMI) and introduced Western Duroc (DU) and Large White (LW) pigs. By pairwise comparison (BMI vs DU, BMI vs LW, and DU vs LW), we found the gene expression differences in the testes between Chinese local pigs and introduced Western commercial breeds were more striking than those between introduced commercial breeds. Furthermore, we found 1622 co-differentially expressed genes (co-DEGs), 122 co-differentially expressed lncRNAs (co-DELs), 39 co-differentially expressed miRNAs (co-DEMs) in BMI vs introduced commercial breeds (DU and LW). Functional analysis revealed that these co-DEGs and co-DELs/co-DEMs target genes were enriched in male sexual function pathways, including MAPK, AMPK, TGF-β/Smad, Hippo, NF-kappa B, and PI3K/Akt signaling pathways. Additionally, we established 10,536 lncRNA-mRNA, 11,248 miRNA-mRNA pairs, and 62 ceRNA (lncRNA-miRNA-mRNA) networks. The ssc-miR-1343 had the most interactive factors in the ceRNA network, including 20 mRNAs and 3 lncRNAs, consisting of 56 ceRNA pairs. These factors played extremely important roles in the regulation of testis function as key nodes in the interactive regulatory network. Our results provide insight into the functional roles of lncRNAs and miRNAs in porcine testis and offer a valuable resource for understanding the differences between Chinese indigenous and introduced Western pigs.

Identification of genomic characteristics and selective signals in Guizhou black goat

BACKGROUND

Guizhou black goat is one of the indigenous black goat breeds in the southwest region of Guizhou, China, which is an ordinary goat for mutton production. They are characterized by moderate body size, black coat, favorite meat quality with tender meat and lower odor, and tolerance for cold and crude feed. However, little is known about the genetic characteristics or variations underlying their important economic traits.

RESULTS

Here, we resequenced the whole genome of Guizhou black goat from 30 unrelated individuals breeding in the five core farms. A total of 9,835,610 SNPs were detected, and 2,178,818 SNPs were identified specifically in this breed. The population structure analysis revealed that Guizhou black goat shared a common ancestry with Shaanbei white cashmere goat (0.146), Yunshang black goat (0.103), Iran indigenous goat (0.054), and Moroccan goat (0.002). However, Guizhou black goat showed relatively higher genetic diversity and a lower level of linkage disequilibrium than the other seven goat breeds by the analysis of the nucleotide diversity, linkage disequilibrium decay, and runs of homozygosity. Based on FST and θπ values, we identified 645, 813, and 804 selected regions between Guizhou black goat and Yunshang black goat, Iran indigenous goat, and cashmere goats. Combined with the results of XP-EHH, there were 286, 322, and 359 candidate genes, respectively. Functional annotation analysis revealed that these genes are potentially responsible for the immune response (e.g., CD28, CD274, IL1A, TLR2, and SLC25A31), humility-cold resistance (e.g., HBEGF, SOSTDC1, ARNT, COL4A1/2, and EP300), meat quality traits (e.g., CHUK, GAB2, PLAAT3, and EP300), growth (e.g., GAB2, DPYD, and CSF1), fertility (e.g., METTL15 and MEI1), and visual function (e.g., PANK2 and NMNAT2) in Guizhou black goat.

CONCLUSION

Our results indicated that Guizhou black goat had a high level of genomic diversity and a low level of linkage disequilibrium in the whole genome. Selection signatures were detected in the genomic regions that were mainly related to growth and development, meat quality, reproduction, disease resistance, and humidity-cold resistance in Guizhou black goat. These results would provide a basis for further resource protection and breeding improvement of this very local breed.

Whole-genome resequencing reveals the uniqueness of Subei yak

Subei yak is an essential local yak in the Gansu Province, which genetic resource has recently been discovered. It is a meat-milk dual-purpose variety with high fecundity and relatively stable population genetic structure. However, its population genetic structure and genetic diversity are yet to be reported. Therefore, this study aimed to identify molecular markers of Subei yak genome by whole-genome resequencing, and to analyze the population structure and genetic diversity of Subei yak. This study screened 12,079,496 single nucleotide polymorphism (SNP) molecular markers in the 20 Subei yaks genome using whole-genome resequencing technology. Of these SNPs, 32.09% were located in the intronic region of the genome. Principal component analysis, phylogenetic analysis, and population structure analysis revealed that the Subei yak belonged to an independent group in the domestic yak population. A selective clearance analysis was carried out on Subei yak and other domestic yaks, and the genes under positive selection were annotated. The functional enrichment analysis showed that Subei yak possessed prominent selection characteristics in terms of external environment perception, hypoxia adaptation, and muscle development. Furthermore, Subei yak showed excellent muscle fat deposition and meat quality traits. Thus, this study will serve as a reference for discovering population structure, genetic evolution, and other unique traits of Subei yak and for expanding the genetic variation catalog of yaks.

Comparative muscle transcriptome of Mali and Hampshire breeds of pigs: a preliminary study

Muscle development is an important priority of pig breeding programs. There is a considerable variation in muscularity between the breeds, but the regulation mechanisms of genes underlying myogenesis are still unclear. Transcriptome data from two breeds of pigs with divergent muscularity (Mali and Hampshire) were integrated with histology, immunofluorescence and meat yield to identify differences in myogenesis during the early growth phase. The muscle transcriptomics analysis revealed 17,721 common, 1413 and 1115 unique transcripts to Hampshire and Mali, respectively. This study identified 908 differentially expressed genes (p < 0.05; log2FC > ±1) in the muscle samples, of which 550 were upregulated and 358 were downregulated in Hampshire pigs, indicating differences in physiological process related to muscle function and development. Expression of genes related to myoblast fusion (MYMK), skeletal muscle satellite cell proliferation (ANGPT1, CDON) and growth factors (HGF, IGF1, IGF2) were higher in Hampshire than Mali, even though transcript levels of several other myogenesis-related genes (MYF6, MYOG, MSTN) were similar. The number of fibers per fascicle and the expression of myogenic marker proteins (MYOD1, MYOG and PAX7) were more in Hampshire as compared to Mali breed of pig, supporting results of transcriptome studies. The results suggest that differences in muscularity between breeds could be related to the regulation of myoblast fusion and myogenic activities. The present study will help to identify genes that could be explored for their utility in the selection of animals with different muscularities.

Transcriptome Analysis Revealed Potential Genes of Skeletal Muscle Thermogenesis in Mashen Pigs and Large White Pigs under Cold Stress

Pigs are susceptible to cold stress due to the absence of brown fat caused by the partial deletion of uncoupling protein 1 during their evolution. Some local pig breeds in China exhibit potential cold adaptability, but research has primarily focused on fat and intestinal tissues. Skeletal muscle plays a key role in adaptive thermogenesis in mammals, yet the molecular mechanism of cold adaptation in porcine skeletal muscle remains poorly understood. This study investigated the cold adaptability of two pig breeds, Mashen pigs (MS) and Large White pigs (LW), in a four-day cold (4 °C) or normal temperature (25 °C) environment. We recorded phenotypic changes and collected blood and longissimus dorsi muscle for transcriptome sequencing. Finally, the PRSS8 gene was randomly selected for functional exploration in porcine skeletal muscle satellite cells. A decrease in body temperature and body weight in both LW and MS pigs under cold stress, accompanied by increased shivering frequency and respiratory frequency, were observed. However, the MS pigs demonstrated stable physiological homeostasis, indicating a certain level of cold adaptability. The LW pigs primarily responded to cold stress by regulating their heat production and glycolipid energy metabolism. The MS pigs exhibited a distinct response to cold stress, involving the regulation of heat production, energy metabolism pathways, and robust mitochondrial activity, as well as a stronger immune response. Furthermore, the functional exploration of PRSS8 in porcine skeletal muscle satellite cells revealed that it affected cellular energy metabolism and thermogenesis by regulating ERK phosphorylation. These findings shed light on the diverse transcriptional responses of skeletal muscle in LW and MS pigs under cold stress, offering valuable insights into the molecular mechanisms underlying cold adaptation in pigs.

A genome-wide association study for loin depth and muscle pH in pigs from intensely selected purebred lines

BACKGROUND

Genome-wide association studies (GWAS) aim at identifying genomic regions involved in phenotype expression, but identifying causative variants is difficult. Pig Combined Annotation Dependent Depletion (pCADD) scores provide a measure of the predicted consequences of genetic variants. Incorporating pCADD into the GWAS pipeline may help their identification. Our objective was to identify genomic regions associated with loin depth and muscle pH, and identify regions of interest for fine-mapping and further experimental work. Genotypes for ~ 40,000 single nucleotide morphisms (SNPs) were used to perform GWAS for these two traits, using de-regressed breeding values (dEBV) for 329,964 pigs from four commercial lines. Imputed sequence data was used to identify SNPs in strong ([Formula: see text] 0.80) linkage disequilibrium with lead GWAS SNPs with the highest pCADD scores.

RESULTS

Fifteen distinct regions were associated with loin depth and one with loin pH at genome-wide significance. Regions on chromosomes 1, 2, 5, 7, and 16, explained between 0.06 and 3.55% of the additive genetic variance and were strongly associated with loin depth. Only a small part of the additive genetic variance in muscle pH was attributed to SNPs. The results of our pCADD analysis suggests that high-scoring pCADD variants are enriched for missense mutations. Two close but distinct regions on SSC1 were associated with loin depth, and pCADD identified the previously identified missense variant within the MC4R gene for one of the lines. For loin pH, pCADD identified a synonymous variant in the RNF25 gene (SSC15) as the most likely candidate for the muscle pH association. The missense mutation in the PRKAG3 gene known to affect glycogen content was not prioritised by pCADD for loin pH.

CONCLUSIONS

For loin depth, we identified several strong candidate regions for further statistical fine-mapping that are supported in the literature, and two novel regions. For loin muscle pH, we identified one previously identified associated region. We found mixed evidence for the utility of pCADD as an extension of heuristic fine-mapping. The next step is to perform more sophisticated fine-mapping and expression quantitative trait loci (eQTL) analysis, and then interrogate candidate variants in vitro by perturbation-CRISPR assays.

Identification and characterization of circRNAs related to meat quality during embryonic development of the longissimus dorsi muscle in two pig breeds

Meat quality, an important economic trait, is regulated by many factors, especially by genetic factors, including coding genes, miRNAs, and lncRNAs. Recent studies have elucidated that circRNAs also play a key role in muscle development and lipid deposition. However, the functions and regulatory mechanisms of circRNAs in meat quality remain mostly unknown. The circRNA expression profiles between Huainan pigs (Chinese indigenous pigs, fat-type, Huainan HN) and Large White pigs (Western commercial pigs, lean-type, LW) in the longissimus dorsi (LD) muscle at 38, 58, and 78 days post conception (dpc) were compared by sequencing. In total, 39,887 circRNAs were identified in 18 samples, and 60, 78, and 86 differentially expressed circRNAs (DECs) were found at the three stages mentioned above between these two breeds. The parent genes of DECs were enriched in myogenesis, proliferation, adipogenesis and muscle fiber-type transition. The circRNA-miRNA interaction networks included 38 DECs and 47 miRNAs, and these miRNAs were involved in muscle development and lipid metabolism. Two shared DECs (circ_0030593 and circ_0032760) of these three stages were selected, their head-to-tail junction sites were validated by Sanger sequencing, and RT‒qPCR results suggested that these two DECs might be involved in intramuscular fat deposition. These findings provide a basis for understanding the role of circRNAs in meat quality.

Effects of Low-Ambient-Temperature Stimulation on Modifying the Intestinal Structure and Function of Different Pig Breeds

Simple Summary

Low ambient temperature resulted in the body’s cold stress response, while local wild boars in the middle-temperate zone performed better than commercial pigs. Therefore, three breeds—Large White (LW) pigs, a local Mashen (MS) pig breed and Jinfen White (JFW) pigs, a hybrid breed from wild boar—were investigated in an artificial climate chamber. The results implicated that low-ambient-temperature stimulation increased trypsin activity in duodenal chyme and promoted inflammatory response in Mashen pigs. The cold-resistance mechanism of MS pigs should be explored to reduce hogs’ stress caused by low-ambient-temperature stimulation.

Abstract

Ambient temperature (Ta) fluctuation is a key factor affecting the growth performance and economic returns of pigs. However, whether the response of intestinal structure and function are related to pig breeds in low Ta has not been investigated yet. In this study, Large White (LW) pigs, Jinfen White (JFW) pigs and Mashen (MS) pigs were raised in artificial climate chambers under normal Ta (25 °C) and low Ta (4 °C) for 96 h. Afterwards, the decrease in body temperature and complete blood counts (CBC) of all pigs were measured. Hematoxylin–eosin, immunohistochemical staining, qPCR and ELISA were used to investigate their intestinal mucosa integrity and inflammatory response. The results showed that MS pigs could maintain a normal body temperature and villus structure after 4 °C stimulation compared with those of LW and JFW pigs. Villus height and villus height/crypt depth of MS pigs were significantly higher than those of LW and JFW pigs at 4 °C. Low-Ta stimulation increased the digestion of carbohydrates of all pigs. Meanwhile, low Ta enhanced the activity of lipase in LW pigs and increased trypsin activity in MS and JFW pigs. Furthermore, low-Ta stimulation significantly downregulated the protein of tight junction and upregulated the mRNA expression of inflammatory cytokines in MS pigs. MS pigs also showed stronger spleen immune function at 4 °C. These results indicated that the local MS pig breed had stronger intestinal function in low Ta by producing a stronger inflammatory response, which lays the foundation for further study on the mechanism of cold tolerance in pigs.

Genome-Wide Expression Profiling of mRNAs, lncRNAs and circRNAs in Skeletal Muscle of Two Different Pig Breeds

Simple Summary

Variation exists in muscle-related traits, such as muscle growth and meat quality, between obese and lean pigs. In this study, the transcriptome profiles of skeletal muscle between Beijing Blackand Yorkshire pigs were characterized to explore the molecular mechanism underlying skeletal muscle-relatedtraits. Gene Ontology (GO) and KEGG pathway enrichment analyses showed that differentially expressed mRNAs, lncRNAs, and circRNAs involved in skeletal muscle development and fatty acid metabolism played a key role in the determination of muscle-related traits between different pig breeds. These results provide candidate genes responsible for muscle phenotypic variation and are valuable for pig breeding.

Abstract

RNA-Seq technology is widely used to analyze global changes in the transcriptome and investigate the influence on relevant phenotypic traits. Beijing Black pigs show differences in growth rate and meat quality compared to western pig breeds. However, the molecular mechanisms responsible for such phenotypic differences remain unknown. In this study, longissimus dorsi muscles from Beijing Black and Yorkshire pigs were used to construct RNA libraries and perform RNA-seq. Significantly different expressions were observed in 1051 mRNAs, 322 lncRNAs, and 82 circRNAs. GO and KEGG pathway annotation showed that differentially expressed mRNAs participated in skeletal muscle development and fatty acid metabolism, which determined the muscle-related traits. To explore the regulatory role of lncRNAs, the cis and trans-target genes were predicted and these lncRNAswere involved in the biological processes related to skeletal muscle development and fatty acid metabolismvia their target genes. CircRNAs play a ceRNA role by binding to miRNAs. Therefore, the potential miRNAs of differentially expressed circRNAs were predicted and interaction networks among circRNAs, miRNAs, and key regulatory mRNAs were constructed to illustrate the function of circRNAs underlying skeletal muscle development and fatty acid metabolism. This study provides new clues for elucidating muscle phenotypic variation in pigs.

Integrative ATAC-seq and RNA-seq Analysis of the Longissimus Muscle of Luchuan and Duroc Pigs

Luchuan pig is a typical obese pig breed in China, and the diameter and area of its longissimus dorsi muscle fibers are significantly smaller than those of Duroc (lean) pig. Skeletal muscle fiber characteristics are related to meat quality of livestock. There is a significant correlation between the quality of different breeds of pork and the characteristics of muscle fiber, which is an important factor affecting the quality of pork. The diameter and area of muscle fibers are related to muscle growth and development. Therefore, we used the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) analysis to investigate the potential mechanism underlying the difference in skeletal muscle growth and development between the two types of pigs. First, transposase-accessible chromatin was analyzed to map the landscape of open chromatin regions and transcription factor binding sites. We identified several transcription factors that potentially affected muscle growth and development, including TFAP4, MAX, NHLH1, FRX5, and TGIF1. We also found that transcription factors with basic helix-loop-helix structures had a preference for binding to genes involved in muscle development. Then, by integrating ATAC-seq and RNA-seq, we found that the Wnt signaling pathway, the mTOR signaling pathway, and other classical pathways regulate skeletal muscle development. In addition, some pathways that might regulate skeletal muscle growth, such as parathyroid hormone synthesis, secretion, and action, synthesis and degradation of ketone bodies, and the thyroid hormone signaling pathway, which were significantly enriched. After further study, we identified a number of candidate genes (ASNS, CARNS1, G0S2, PPP1R14C, and SH3BP5) that might be associated with muscle development. We also found that the differential regulation of chromatin openness at the level of some genes was contrary to the differential regulation at the level of transcription, suggesting that transcription factors and transcriptional repressors may be involved in the regulation of gene expression. Our study provided an in-depth understanding of the mechanism behind the differences in muscle fibers from two species of pig and provided an important foundation for further research on improving the quality of pork.

In Utero Fetal Weight in Pigs Is Regulated by microRNAs and Their Target Genes

Impaired skeletal muscle growth in utero can result in reduced birth weight and poor carcass quality in pigs. Recently, we showed the role of microRNAs (miRNAs) and their target genes in prenatal skeletal muscle development and pathogenesis of intrauterine growth restriction (IUGR). In this study, we performed an integrative miRNA-mRNA transcriptomic analysis in longissimus dorsi muscle (LDM) of pig fetuses at 63 days post conception (dpc) to identify miRNAs and genes correlated to fetal weight. We found 13 miRNAs in LDM significantly correlated to fetal weight, including miR-140, miR-186, miR-101, miR-15, miR-24, miR-29, miR-449, miR-27, miR-142, miR-99, miR-181, miR-199, and miR-210. The expression of these miRNAs decreased with an increase in fetal weight. We also identified 1315 genes significantly correlated to fetal weight at 63 dpc, of which 135 genes were negatively correlated as well as identified as potential targets of the above-listed 13 miRNAs. These miRNAs and their target genes enriched pathways and biological processes important for fetal growth, development, and metabolism. These results indicate that the transcriptomic profile of skeletal muscle can be used to predict fetal weight, and miRNAs correlated to fetal weight can serve as potential biomarkers of prenatal fetal health and growth.

Evidence Against the Causal Relationship Between a Putative Cis-Regulatory Variant of MYH3 and Intramuscular Fat Content in Pigs

Improving meat quality has become the main goal of modern pig breeding. Intramuscular fat content (IMF) is an important trait influencing meat quality of livestock, but the molecular mechanism behind this trait is still unclear. Recently, Cho et al. reported the discovery of the first causal mutation affecting IMF and red flesh color (a*) in pigs, namely XM_013981330.2:g.−1805_−1810del, a 6-bp deletion variant in the porcine MYH3 promoter region. The objective of this study was to reassess the causality of this mutation for its potential commercial application. By Sanger sequencing, we firstly identified several new variants (including a 4-bp deletion) at or near the 6-bp deletion site, which formed four haplotypes in multiple breeds. Unexpectedly, the 6-bp deletion allele, previously determined as the MYH3 Q allele because of its significantly positive effect on IMF and a*, was found not only in Chinese indigenous breeds, but also in four western commercial breeds with relatively lower IMF levels, including Duroc, Large White, Landrace and Pietrain. More surprisingly, we found that the MYH3 Q allele and the haplotypes harboring it had no significant effects on IMF, marbling and color score in three large-scale divergent pig populations: the heterogeneous F6 and F7 pigs and commercial crossbred Duroc × (Landrace × Yorkshire) pigs. Transient transfection analysis in porcine satellite cells showed that the 6-bp deletion variants had a negligible effect on transcription of reporter gene, but could attenuate the MRF (myogenesis regulatory factors)-induced increase in luciferase activity of the MYH3 promoter vector. The MYH3 protein level in muscle did not differ significantly among the haplotype groups. Therefore, our results cannot support the causal relationship between the 6-bp deletion in MYH3 and IMF trait, suggesting that the causal mutation for the IMF QTL on SSC12 needs to be further identified.