Background: The genetic mechanisms that underlie phenotypic differentiation in breeding animals have important implications in evolutionary biology and agriculture. However, the contribution of cis-regulatory variants to pig phenotypes is poorly understood. Therefore, our aim was to elucidate the molecular mechanisms by which non-coding variants cause phenotypic differences in pigs by combining evolutionary biology analyses and functional genomics. Results: We obtained a high-resolution phased chromosome-scale reference genome with a contig N50 of 18.03 Mb for the Luchuan pig breed (a representative eastern breed) and profiled potential selective sweeps in eastern and western pigs by resequencing the genomes of 234 pigs. Multi-tissue transcriptome and chromatin accessibility analyses of these regions suggest that tissue-specific selection pressure is mediated by promoters and distal cis-regulatory elements. Promoter variants that are associated with increased expression of the lysozyme ( LYZ ) gene in the small intestine might enhance the immunity of the gastrointestinal tract and roughage tolerance in pigs. In skeletal muscle, an enhancer-modulating single-nucleotide polymorphism that is associated with up-regulation of the expression of the troponin C1, slow skeletal and cardiac type ( TNNC1 ) gene might increase the proportion of slow muscle fibers and affect meat quality. Conclusions: Our work sheds light on the molecular mechanisms by which non-coding variants shape phenotypic differences in pigs and provides valuable resources and novel perspectives to dissect the role of gene regulatory evolution in animal domestication and breeding.

The genetic basis of complex traits and phenotypic differentiation remains unclear in pigs. Using nine genomes—seven of which were newly generated, high‐quality de novo assembled genomes—and 1081 resequencing genomes, we built a pan‐genome and identified 134.24 Mb nonredundant nonreference sequences, 1099 novel protein‐coding genes, 187,927 structural variations (SVs) and 30,143,962 single‐nucleotide polymorphisms (SNPs). Analysis of selective domestication revealed BRCA1 associated with enhanced adipocyte growth and fat deposition, and ABCA3 linked to an alleviated immune response and reduced lung injury. Integrating 162 transcriptomes and 162 methylomes of skeletal muscle across 27 developmental stages revealed the regulatory mechanism of phenotypic differentiation between Eastern and Western breeds. Artificial selection reshaped local DNA methylation status and imparted regulatory effects on the progression patterns of heterochronic genes such as GHSR and BDH1 , particularly during embryonic development. Altogether, our work provides valuable resources for understanding molecular mechanisms behind phenotypic variations and enhancing the genetic improvement programs in pigs.

Myogenesis is a central step in prenatal myofiber formation, postnatal myofiber hypertrophy, and muscle damage repair in adulthood. RNA-Seq technology has greatly helped reveal the molecular mechanism of myogenesis, but batch effects in different experiments inevitably lead to misinterpretation of differentially expressed genes (DEGs). We previously applied the robust rank aggregation (RRA) method to effectively circumvent batch effects across multiple RNA-Seq datasets from 3T3-L1 cells. Here, we also used the RRA method to integrate nine RNA-Seq datasets from C2C12 cells and obtained 3140 robust DEGs between myoblasts and myotubes, which were then validated with array expression profiles and H3K27ac signals. The upregulated robust DEGs were highly enriched in gene ontology (GO) terms related to muscle cell differentiation and development. Considering that the cooperative binding of transcription factors (TFs) to enhancers to regulate downstream gene expression is a classical epigenetic mechanism, differentially expressed TFs (DETFs) were screened, and potential novel myogenic factors (MAF, BCL6, and ESR1) with high connection degree in protein–protein interaction (PPI) network were presented. Moreover, KLF5 cooperatively binds with the three key myogenic factors (MYOD, MYOG, and MEF2D) in C2C12 cells. Motif analysis speculates that the binding of MYOD and MYOG is KLF5-independent, while MEF2D is KLF5-dependent. It was revealed that KLF5-binding sites could be exploited to filter redundant MYOD-, MYOG-, and MEF2D-binding sites to focus on key enhancers for myogenesis. Further functional annotation of KLF5-binding sites suggested that KLF5 may regulate myogenesis through the PI3K-AKt signaling pathway, Rap1 signaling pathway, and the Hippo signaling pathway. In general, our study provides a wealth of untapped candidate targets for myogenesis and contributes new insights into the core regulatory mechanisms of myogenesis relying on KLF5-binding signal.

Silencers, the yin to enhancers’ yang, play a pivotal role in fine-tuning gene expression throughout the genome. However, despite their recognized importance, comprehensive identification of these regulatory elements in the genome is still in its early stages. We developed a method called Ss-STARR-seq to directly determine the activity of silencers in the whole genome. In this study, we applied Ss-STARR-seq to human cell lines K562, LNCaP, and 293 T, and identified 134,171, 137,753, and 125,307 silencers on a genome-wide scale, respectively, these silencers function in various cells in a cell-specific manner. Silencers exhibited a substantial enrichment of transcriptional-inhibitory motifs, including REST, and demonstrated overlap with the binding sites of repressor transcription factors within the endogenous environment. Interestingly, H3K27me3 did not reflect silencer activity but facilitated the silencer’s inhibitory role on gene expression. Additionally, the silencer did not have any significant histone markers at the genome-wide level. Our findings unveil that aspect-silencers not only transition into enhancers throughout diverse cell lines but also achieve functional conversion with insulators. Regarding to biological effects, knockout experiments underscored the functional redundancy and specificity of silencers in regulating gene expression and cell proliferation. In summary, this study pioneers the elucidation of the genome-wide silencer landscape in human cells, delineates their global regulatory features, and identifies specific silencers influencing cancer cell proliferation.

The Bama Xiang pig (BM) is a unique pig species in Guangxi Province, China. Compared to other breeds of domestic pig, such as the Debao pig (DB), it is smaller in size, better in meat quality, resistant to rough feeding and strong in stress resistance. These unique advantages of Bama Xiang pigs make them of great edible value and scientific research value. However, the differences in muscle metabolites between Bama Xiang pigs (BM) and Debao pigs (DB) are largely unexplored. Here, we identified 214 differential metabolites between these two pig breeds by LC-MS. Forty-one such metabolites are enriched into metabolic pathways, and these metabolites correspond to 11 metabolic pathways with significant differences. In Bama pigs, the abundance of various metabolites such as creatine, citric acid, L-valine and hypoxanthine is significantly higher than in Debao pigs, while the abundance of other metabolites, such as carnosine, is significantly lower. Among these, we propose six differential metabolites: L-proline, citric acid, ribose 1-phosphate, L-valine, creatine, and L-arginine, as well as four potential differential metabolites (without the KEGG pathway), alanyl-histidine, inosine 2′-phosphate, oleoylcarnitine, and histidinyl hydroxyproline, as features for evaluating the meat quality of Bama pigs and for differentiating pork from Bama pigs and Debao pigs. This study provides a proof-of-concept example of distinguishing pork from different pig breeds at the metabolite level and sheds light on elucidating the biological processes underlying meat quality differences. Our pork metabolites data are also of great value to the genomics breeding community in meat quality improvement.