Whole-genome analysis reveals distinct adaptation signatures to diverse environments in Chinese domestic pigs.

BACKGROUND: Long-term natural and artificial selection has resulted in many genetic footprints within the genomes of pig breeds across distinct agroecological zones. Nevertheless, the mechanisms by which these signatures contribute to phenotypic diversity and facilitate environmental adaptation remain unclear. RESULTS: Here, we leveraged whole-genome sequencing data from 82 individuals from 6 domestic pig breeds originating in tropical, high-altitude, and frigid regions. Population genetic analysis suggested that habitat isolation significantly shaped the genetic diversity and contributed to population stratification in local Chinese pig breeds. Analysis of selection signals revealed regions under selection for adaptation in tropical (55.5 Mb), high-altitude (43.6 Mb), and frigid (17.72 Mb) regions. The potential functions of the selective sweep regions were linked to certain complex traits that might play critical roles in different geographic environments, including fat coverage in frigid environments and blood indicators in tropical and high-altitude environments. Candidate genes under selection were significantly enriched in biological pathways involved in environmental adaptation. These pathways included blood circulation, protein degradation, and inflammation for adaptation to tropical environments; heart and lung development, hypoxia response, and DNA damage repair for high-altitude adaptation; and thermogenesis, cold-induced vasodilation (CIVD), and the cell cycle for adaptation to frigid environments. By examining the chromatin state of the selection signatures, we identified the lung and ileum as two candidate functional tissues for environmental adaptation. Finally, we identified a mutation (chr1: G246,175,129A) in the cis-regulatory region of ABCA1 as a plausible promising variant for adaptation to tropical environments. CONCLUSIONS: In this study, we conducted a genome-wide exploration of the genetic mechanisms underlying the adaptability of local Chinese pig breeds to tropical, high-altitude, and frigid environments. Our findings shed light on the prominent role of cis-regulatory elements in environmental adaptation in pigs and may serve as a valuable biological model of human plateau-related disorders and cardiovascular diseases.

Haplotype-resolved assembly of a pig genome using single-sperm sequencing.

Single gamete cell sequencing together with long-read sequencing can reliably produce chromosome-level phased genomes. In this study, we employed PacBio HiFi and Hi-C sequencing on a male Landrace pig, coupled with single-sperm sequencing of its 102 sperm cells. A haplotype assembly method was developed based on long-read sequencing and sperm-phased markers. The chromosome-level phased assembly showed higher phasing accuracy than methods that rely only on HiFi reads. The use of single-sperm sequencing data enabled the construction of a genetic map, successfully mapping the sperm motility trait to a specific region on chromosome 1 (105.40-110.70 Mb). Furthermore, with the assistance of Y chromosome-bearing sperm data, 26.16 Mb Y chromosome sequences were assembled. We report a reliable approach for assembling chromosome-level phased genomes and reveal the potential of sperm population in basic biology research and sperm phenotype research.

CRISPR Screen Identifies the RNA-Binding Protein Eef1a1 as a Key Regulator of Myogenesis.

Skeletal muscle myogenesis hinges on gene regulation, meticulously orchestrated by molecular mechanisms. While the roles of transcription factors and non-coding RNAs in myogenesis are widely known, the contribution of RNA-binding proteins (RBPs) has remained unclear until now. Therefore, to investigate the functions of post-transcriptional regulators in myogenesis and uncover new functional RBPs regulating myogenesis, we employed CRISPR high-throughput RBP-KO (RBP-wide knockout) library screening. Through this approach, we successfully identified Eef1a1 as a novel regulatory factor in myogenesis. Using CRISPR knockout (CRISPRko) and CRISPR interference (CRISPRi) technologies, we successfully established cellular models for both CRISPRko and CRISPRi. Our findings demonstrated that Eef1a1 plays a crucial role in promoting proliferation in C2C12 myoblasts. Through siRNA inhibition and overexpression methods, we further elucidated the involvement of Eef1a1 in promoting proliferation and suppressing differentiation processes. RIP (RNA immunoprecipitation), miRNA pull-down, and Dual-luciferase reporter assays confirmed that miR-133a-3p targets Eef1a1. Co-transfection experiments indicated that miR-133a-3p can rescue the effect of Eef1a1 on C2C12 myoblasts. In summary, our study utilized CRISPR library high-throughput screening to unveil a novel RBP, Eef1a1, involved in regulating myogenesis. Eef1a1 promotes the proliferation of myoblasts while inhibiting the differentiation process. Additionally, it acts as an antagonist to miR-133a-3p, thus modulating the process of myogenesis.

Regulation of myo-miR-24-3p on the Myogenesis and Fiber Type Transformation of Skeletal Muscle.

Skeletal muscle plays critical roles in providing a protein source and contributing to meat production. It is well known that microRNAs (miRNAs) exert important effects on various biological processes in muscle, including cell fate determination, muscle fiber morphology, and structure development. However, the role of miRNA in skeletal muscle development remains incompletely understood. In this study, we observed a critical miRNA, miR-24-3p, which exhibited higher expression levels in Tongcheng (obese-type) pigs compared to Landrace (lean-type) pigs. Furthermore, we found that miR-24-3p was highly expressed in the dorsal muscle of pigs and the quadriceps muscle of mice. Functionally, miR-24-3p was found to inhibit proliferation and promote differentiation in muscle cells. Additionally, miR-24-3p was shown to facilitate the conversion of slow muscle fibers to fast muscle fibers and influence the expression of GLUT4, a glucose transporter. Moreover, in a mouse model of skeletal muscle injury, we demonstrated that overexpression of miR-24-3p promoted rapid myogenesis and contributed to skeletal muscle regeneration. Furthermore, miR-24-3p was found to regulate the expression of target genes, including Nek4, Pim1, Nlk, Pskh1, and Mapk14. Collectively, our findings provide evidence that miR-24-3p plays a regulatory role in myogenesis and fiber type conversion. These findings contribute to our understanding of human muscle health and have implications for improving meat production traits in livestock.

Comprehensive Atlas of Alternative Splicing Reveals NSRP1 Promoting Adipogenesis through CCDC18.

Alternative splicing (AS) plays a crucial role in regulating gene expression, function, and diversity. However, limited reports exist on the identification and comparison of AS in Eastern and Western pigs. Here, we analyzed 243 transcriptome data from eight tissues, integrating information on transcription factors (TFs), selection signals, splicing factors (SFs), and quantitative trait loci (QTL) to comprehensively study alternative splicing events (ASEs) in pigs. Five ASE types were identified, with Mutually Exclusive Exon (MXE) and Skipped Exon (SE) ASEs being the most prevalent. A significant portion of genes with ASEs (ASGs) showed conservation across all eight tissues (63.21-76.13% per tissue). Differentially alternative splicing genes (DASGs) and differentially expressed genes (DEGs) exhibited tissue specificity, with blood and adipose tissues having more DASGs. Functional enrichment analysis revealed coDASG_DEGs in adipose were enriched in pathways associated with adipose deposition and immune inflammation, while coDASG_DEGs in blood were enriched in pathways related to immune inflammation and metabolism. Adipose deposition in Eastern pigs might be linked to the down-regulation of immune-inflammation-related pathways and reduced insulin resistance. The TFs, selection signals, and SFs appeared to regulate ASEs. Notably, ARID4A (TF), NSRP1 (SF), ANKRD12, IFT74, KIAA2026, CCDC18, NEXN, PPIG, and ROCK1 genes in adipose tissue showed potential regulatory effects on adipose-deposition traits. NSRP1 could promote adipogenesis by regulating alternative splicing and expression of CCDC18. Conducting an in-depth investigation into AS, this study has successfully identified key marker genes essential for pig genetic breeding and the enhancement of meat quality, which will play important roles in promoting the diversity of pork quality and meeting market demand.

Characterizing core microbiota and regulatory functions of the pig gut microbiome.

Domestic pigs (Sus scrofa) are the leading terrestrial animals used for meat production. The gut microbiota significantly affect host nutrition, metabolism, and immunity. Hence, characterization of the gut microbial structure and function will improve our understanding of gut microbial resources and the mechanisms underlying host-microbe interactions. Here, we investigated the gut microbiomes of seven pig breeds using metagenomics and 16S rRNA gene amplicon sequencing. We established an expanded gut microbial reference catalog comprising 17 020 160 genes and identified 4910 metagenome-assembled genomes. We also analyzed the gut resistome to provide an overview of the profiles of the antimicrobial resistance genes in pigs. By analyzing the relative abundances of microbes, we identified three core-predominant gut microbes (Phascolarctobacterium succinatutens, Prevotella copri, and Oscillibacter valericigenes) in pigs used in this study. Oral administration of the three core-predominant gut microbes significantly increased the organ indexes (including the heart, spleen, and thymus), but decreased the gastrointestinal lengths in germ-free mice. The three core microbes significantly enhanced intestinal epithelial barrier function and altered the intestinal mucosal morphology, as was evident from the increase in crypt depths in the duodenum and ileum. Furthermore, the three core microbes significantly affected several metabolic pathways (such as "steroid hormone biosynthesis," "primary bile acid biosynthesis," "phenylalanine, tyrosine and tryptophan biosynthesis," and "phenylalanine metabolism") in germ-free mice. These findings provide a panoramic view of the pig gut microbiome and insights into the functional contributions of the core-predominant gut microbes to the host.

Single-cell profiling of African swine fever virus disease in the pig spleen reveals viral and host dynamics.

African swine fever, one of the major viral diseases of swine, poses an imminent threat to the global pig industry. The high-efficient replication of the causative agent African swine fever virus (ASFV) in various organs in pigs greatly contributes to the disease. However, how ASFV manipulates the cell population to drive high-efficient replication of the virus in vivo remains unclear. Here, we found that the spleen reveals the most severe pathological manifestation with the highest viral loads among various organs in pigs during ASFV infection. By using single-cell-RNA-sequencing technology and multiple methods, we determined that macrophages and monocytes are the major cell types infected by ASFV in the spleen, showing high viral-load heterogeneity. A rare subpopulation of immature monocytes represents the major population infected at late infection stage. ASFV causes massive death of macrophages, but shifts its infection into these monocytes which significantly arise after the infection. The apoptosis, interferon response, and antigen-presentation capacity are inhibited in these monocytes which benefits prolonged infection of ASFV in vivo. Until now, the role of immature monocytes as an important target by ASFV has been overlooked due to that they do not express classical monocyte marker CD14. The present study indicates that the shift of viral infection from macrophages to the immature monocytes is critical for maintaining prolonged ASFV infection in vivo. This study sheds light on ASFV tropism, replication, and infection dynamics, and elicited immune response, which may instruct future research on antiviral strategies.

Exploring Multi-Tissue Alternative Splicing and Skeletal Muscle Metabolism Regulation in Obese- and Lean-Type Pigs.

Alternative splicing (AS) is a crucial mechanism in post-transcriptional regulation, contributing significantly to the diversity of the transcriptome and proteome. In this study, we performed a comprehensive AS profile in nine tissues obtained from Duroc (lean-type) and Luchuan (obese-type) pigs. Notably, 94,990 AS events from 14,393 genes were identified. Among these AS events, it was observed that 80% belonged to the skipped exon (SE) type. Functional enrichment analysis showed that genes with more than ten AS events were closely associated with tissue-specific functions. Additionally, the analysis of overlap between differentially alternative splicing genes (DSGs) and differentially expressed genes (DEGs) revealed the highest number of overlapped genes in the heart and skeletal muscle. The novelty of our study is that it identified and validated three genes (PYGM, MAPK11 and CAMK2B) in the glucagon signaling pathway, and their alternative splicing differences were highly significant across two pig breeds. In conclusion, our study offers novel insights into the molecular regulation of diverse tissue physiologies and the phenotypic differences between obese- and lean-type pigs, which are helpful for pig breeding.

Transcription Factor SATB2 Regulates Skeletal Muscle Cell Proliferation and Migration via HDAC4 in Pigs.

Skeletal muscle development remarkably affects meat production and growth rate, regulated by complex regulatory mechanisms in pigs. Specific AT sequence-binding protein 2 (SATB2) is a classic transcription factor and chromatin organizer, which holds a profound effect in the regulation of chromatin remodeling. However, the regulation role of SATB2 concerning skeletal muscle cell fate through chromatin remodeling in pigs remains largely unknown. Here, we observed that SATB2 was expressed higher in the lean-type compared to the obese-type pigs, which also enriched the pathways of skeletal muscle development, chromatin organization, and histone modification. Functionally, knockdown SATB2 led to decreases in the proliferation and migration markers at the mRNA and protein expression levels, respectively, while overexpression SATB2 had the opposite effects. Further, we found histone deacetylase 4 (HDAC4) was a key downstream target gene of SATB2 related to chromatin remodeling. The binding relationship between SATB2 and HDAC4 was confirmed by a dual-luciferase reporter system and ChIP-qPCR analysis. Besides, we revealed that HDAC4 promoted the skeletal muscle cell proliferation and migration at the mRNA and protein expression levels, respectively. In conclusion, our study indicates that transcription factor SATB2 binding to HDAC4 positively contributes to skeletal muscle cell proliferation and migration, which might mediate the chromatin remodeling to influence myogenesis in pigs. This study develops a novel insight into understanding the molecular regulatory mechanism of myogenesis, and provides a promising gene for genetic breeding in pigs.

How changes in landscape patterns affect the carbon emission: a case study in the Chengdu-Chongqing Economic Circle, China.

The construction of low-carbon cities is an optimal means to balance the competing interests of economic growth and carbon emission reduction. This study focuses on the optimization of land use patterns with a low carbon orientation, taking the Chengdu-Chongqing Economic Circle (CCEC), the fourth-largest economic growth pole in China, as an example. The panel data regression analysis is carried out to identify the dynamic correlations between the landscape changes and the carbon emission induced by land use and land cover change (LICE) of each city, each year, for the last 20 years. The results show that the CCEC has witnessed a 142.85% increase in carbon emissions during the period studied, with the growth of built-up land contributing 94% of total carbon emissions from 2000 to 2020. By constructing the panel regression model, this study finds that the intensity of carbon emissions increases significantly as the urban built-up land area and the agglomeration of artificial structures increase. The conversion of cropland, which dominates the landscape pattern, to built-up land has led to further fragmentation of the landscape pattern and a reduction in LPI, thus increasing carbon emissions. And a more complex regional landscape pattern will have a positive impact on carbon emission reduction. Based on the above findings, suggestions are articulated for carbon emission reduction.

Genome-Wide Detection of Copy Number Variations Associated with Miniature Features in Horses.

Copy number variations (CNVs) are crucial structural genomic variants affecting complex traits in humans and livestock animals. The current study was designed to conduct a comprehensive comparative copy number variation analysis among three breeds, Debao (DB), Baise (BS), and Warmblood (WB), with a specific focus on identifying genomic regions associated with miniature features in horses. Using whole-genome next-generation resequencing data, we identified 18,974 CNVs across 31 autosomes. Among the breeds, we found 4279 breed-specific CNV regions (CNVRs). Baise, Debao, and Warmblood displayed 2978, 986, and 895 distinct CNVRs, respectively, with 202 CNVRs shared across all three breeds. After removing duplicates, we obtained 1545 CNVRs from 26 horse genomes. Functional annotation reveals enrichment in biological functions, including antigen processing, cell metabolism, olfactory conduction, and nervous system development. Debao horses have 970 genes overlapping with CNVRs, possibly causing their small size and mountainous adaptations. We also found that the genes GHR, SOX9, and SOX11 may be responsible for the miniature features of the Debao horse by analyzing their overlapping CNVRs. Overall, this study offers valuable insights into the widespread presence of CNVs in the horse genome. The findings contribute to mapping horse CNVs and advance research on unique miniature traits observed in the Debao horse.

Circular RNA transcriptome across multiple tissues reveal skeletal muscle-specific circPSME4 regulating myogenesis.

There are significant differences in meat production, growth rate and other traits between Western commercial pigs and Chinese local pigs. Comparative transcriptome approaches have identified many coding and non-coding candidate genes associated various traits. However, the expression and function of circular RNAs (circRNAs) in different pig tissues are largely unknown. In this study, we conducted a comprehensive analysis of the genome-wide circRNA expression profile across ten tissues in Luchuan (a Chinese local breed) and Duroc (a Western commercial breed) pigs. We identified a total of 56,254 circRNAs, of which 42.9 % were not previously annotated. We found that 33.7 % of these circRNAs were differentially expressed. Enrichment analysis revealed that differentially expressed circRNAs might contribute to the phenotypic differentiation between Luchuan and Duroc pigs. We identified 538 tissue-specific circRNAs, most of which were specifically expressed in the brain and skeletal muscle. Competitive endogenous RNA network analysis suggested that skeletal muscle-specific circPSME4 was co-expressed with MYOD1 and targeted by ssc-miR-181d-3p. Functional analysis revealed that circPSME4 knockdown could promote the proliferation and differentiation of myoblasts. Together, our findings provide valuable resources of circRNAs for animal breeding and biomedical research. We demonstrated that circPSME4 is a novel regulator of skeletal muscle development.

Integrative metabolomic and transcriptomic analysis reveals difference in glucose and lipid metabolism in the longissimus muscle of Luchuan and Duroc pigs.

Luchuan pig, an obese indigenous Chinese porcine breed, has a desirable meat quality and reproductive capacity. Duroc, a traditional western breed, shows a faster growth rate, high feed efficiency and high lean meat rate. Given the unique features these two porcine breeds have, it is of interest to investigate the underlying molecular mechanisms behind their distinctive nature. In this study, the metabolic and transcriptomic profiles of longissimus dorsi muscle from Duroc and Luchuan pigs were compared. A total of 609 metabolites were identified, 77 of which were significantly decreased in Luchuan compared to Duroc, and 71 of which were significantly elevated. Most differentially accumulated metabolites (DAMs) upregulated in Luchuan were glycerophospholipids, fatty acids, oxidized lipids, alcohols, and amines, while metabolites downregulated in Luchuan were mostly amino acids, organic acids and nucleic acids, bile acids and hormones. From our RNA-sequencing (RNA-seq) data we identified a total of 3638 differentially expressed genes (DEGs), 1802 upregulated and 1836 downregulated in Luchuan skeletal muscle compared to Duroc. Combined multivariate and pathway enrichment analyses of metabolome and transcriptome results revealed that many of the DEGs and DAMs are associated with critical energy metabolic pathways, especially those related to glucose and lipid metabolism. We examined the expression of important DEGs in two pathways, AMP-activated protein kinase (AMPK) signaling pathway and fructose and mannose metabolism, using Real-Time Quantitative Reverse Transcription PCR (qRT-PCR). Genes related to glucose uptake, glycolysis, glycogen synthesis, fatty acid synthesis (PFKFB1, PFKFB4, MPI, TPI1, GYS1, SLC2A4, FASN, IRS1, ULK1) are more activated in Luchuan, while genes related to fatty acid oxidation, cholesterol synthesis (CPT1A, HMGCR, FOXO3) are more suppressed. Energy utilization can be a decisive factor to the distinctive metabolic, physiological and nutritional characteristics in skeletal muscle of the two breeds we studied. Our research may facilitate future porcine breeding projects and can be used to reveal the potential molecular basis of differences in complex traits between various breeds.

MicroRNA-200c-5p Regulates Migration and Differentiation of Myoblasts via Targeting Adamts5 in Skeletal Muscle Regeneration and Myogenesis.

Skeletal muscle, as a regenerative organization, plays a vital role in physiological characteristics and homeostasis. However, the regulation mechanism of skeletal muscle regeneration is not entirely clear. miRNAs, as one of the regulatory factors, exert profound effects on regulating skeletal muscle regeneration and myogenesis. This study aimed to discover the regulatory function of important miRNA miR-200c-5p in skeletal muscle regeneration. In our study, miR-200c-5p increased at the early stage and peaked at first day during mouse skeletal muscle regeneration, which was also highly expressed in skeletal muscle of mouse tissue profile. Further, overexpression of miR-200c-5p promoted migration and inhibited differentiation of C2C12 myoblast, whereas inhibition of miR-200c-5p had the opposite effect. Bioinformatic analysis predicted that Adamts5 has potential binding sites for miR-200c-5p at 3'UTR region. Dual-luciferase and RIP assays further proved that Adamts5 is a target gene of miR-200c-5p. The expression patterns of miR-200c-5p and Adamts5 were opposite during the skeletal muscle regeneration. Moreover, miR-200c-5p can rescue the effects of Adamts5 in the C2C12 myoblast. In conclusion, miR-200c-5p might play a considerable function during skeletal muscle regeneration and myogenesis. These findings will provide a promising gene for promoting muscle health and candidate therapeutic target for skeletal muscle repair.

Recent Progress on Circular RNAs in the Development of Skeletal Muscle and Adipose Tissues of Farm Animals.

Circular RNAs (circRNAs) are a highly conserved and specifically expressed novel class of covalently closed non-coding RNAs. CircRNAs can function as miRNA sponges, protein scaffolds, and regulatory factors, and play various roles in development and other biological processes in mammals. With the rapid development of high-throughput sequencing technology, thousands of circRNAs have been discovered in farm animals; some reportedly play vital roles in skeletal muscle and adipose development. These are critical factors affecting meat yield and quality. In this review, we have highlighted the recent advances in circRNA-related studies of skeletal muscle and adipose in farm animals. We have also described the biogenesis, properties, and biological functions of circRNAs. Furthermore, we have comprehensively summarized the functions and regulatory mechanisms of circRNAs in skeletal muscle and adipose development in farm animals and their effects on economic traits such as meat yield and quality. Finally, we propose that circRNAs are putative novel targets to improve meat yield and quality traits during animal breeding.

Coordinated transcriptional and post-transcriptional epigenetic regulation during skeletal muscle development and growth in pigs.

BACKGROUND: N6-methyladenosine (m6A) and DNA 5-methylcytosine (5mC) methylation plays crucial roles in diverse biological processes, including skeletal muscle development and growth. Recent studies unveiled a potential link between these two systems, implicating the potential mechanism of coordinated transcriptional and post-transcriptional regulation in porcine prenatal myogenesis and postnatal skeletal muscle growth. METHODS: Immunofluorescence and co-IP assays were carried out between the 5mC writers and m6A writers to investigate the molecular basis underneath. Large-scale in-house transcriptomic data were compiled for applying weighted correlation network analysis (WGCNA) to identify the co-expression patterns of m6A and 5mC regulators and their potential role in pig myogenesis. Whole-genome bisulfite sequencing (WGBS) and methylated RNA immunoprecipitation sequencing (MeRIP-seq) were performed on the skeletal muscle samples from Landrace pigs at four postnatal growth stages (days 30, 60, 120 and 180). RESULTS: Significantly correlated expression between 5mC writers and m6A writers and co-occurrence of 5mC and m6A modification were revealed from public datasets of C2C12 myoblasts. The protein-protein interactions between the DNA methylase and the m6A methylase were observed in mouse myoblast cells. Further, by analyzing transcriptome data comprising 81 pig skeletal muscle samples across 27 developmental stages, we identified a 5mC/m6A epigenetic module eigengene and decoded its potential functions in pre- or post-transcriptional regulation in postnatal skeletal muscle development and growth of pigs. Following integrative multi-omics analyses on the WGBS methylome data and MeRIP-seq data for both m6A and gene expression profiles revealed a genome/transcriptome-wide correlated dynamics and co-occurrence of 5mC and m6A modifications as a consequence of 5mC/m6A crosstalk in the postnatal myogenesis progress of pigs. Last, we identified a group of myogenesis-related genes collaboratively regulated by both 5mC and m6A modifications in postnatal skeletal muscle growth in pigs. CONCLUSIONS: Our study discloses a potential epigenetic mechanism in skeletal muscle development and provides a novel direction for animal breeding and drug development of related human muscle-related diseases.

Identification of circRNAs Associated with Adipogenesis Based on RNA-seq Data in Pigs.

Adipocytes or fat cells play a vital role in the storage and release of energy in pigs, and many circular RNAs (circRNAs) have emerged as important regulators in various tissues and cell types in pigs. However, the spatio-temporal expression pattern of circRNAs between different adipose deposition breeds remains elusive. In this study, RNA sequencing (RNA-seq) produced transcriptome profiles of Western Landrace (lean-type) and Chinese Songliao black pigs (obese-type) with different thicknesses of subcutaneous fat tissues and were used to identify circRNAs involved in the regulation of adipogenesis. Gene expression analysis revealed 883 circRNAs, among which 26 and 11 circRNAs were differentially expressed between Landrace vs. Songliao pigs and high- vs. low-thickness groups, respectively. We also analyzed the interaction between circRNAs and microRNAs (miRNAs) and constructed their interaction network in adipogenesis; gene ontology classification and pathway analysis revealed two vital circRNAs, with the majority of their target genes enriched in biological functions such as fatty acids biosynthesis, fatty acid metabolism, and Wnt/TGF-ß signaling pathways. These candidate circRNAs can be taken as potential targets for further experimental studies. Our results show that circRNAs are dynamically expressed and provide a valuable basis for understanding the molecular mechanism of circRNAs in pig adipose biology.

Integration of multi-omics data reveals cis-regulatory variants that are associated with phenotypic differentiation of eastern from western pigs.

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.