The Meishan pig genome reveals structural variation-mediated gene expression and phenotypic divergence underlying Asian pig domestication

LncRNA RMG controls liquid-liquid phase separation of MEIS2 to regulate myogenesis

Long non-coding RNAs (lncRNAs) regulate liquid-liquid phase separation (LLPS), driving the formation of biomolecular condensates essential for cellular function. However, this regulatory mechanism is yet to be reported in skeletal muscles. In this study, we comprehensively analyzed lncRNAs in skeletal muscle across multiple pig breeds, developmental stages, and tissues. Our analysis identified over 10,000 novel lncRNAs. We found that the lnc-regulator of muscle growth (lnc-RMG) regulates myogenesis by modulating the LLPS of Meis homeobox 2 (MEIS2). Lnc-RMG was specifically expressed in the skeletal muscle, with significantly higher expression in the fetal stage than in the embryonic stage. Notably, lnc-RMG was highly conserved between pigs and humans and exhibits similar biological functions in myogenesis. Furthermore, lnc-RMG knockdown promoted skeletal muscle regeneration. Mechanistically, lnc-RMG produces mature microRNA (miR)-133a-3p, which targets and inhibits MEIS2 expression, thereby inhibiting MEIS2 LLPS. This inhibition promoted the transcription of transforming growth factor-β receptor II (TGFβR2), ultimately regulating myogenesis. Overall, our findings revealed a novel lnc-RMG/miR-133a-3p/MEIS2/TGFβR2 axis that regulated myogenesis through LLPS and provided new insights into the molecular mechanisms that drive muscle development and regeneration. These findings highlight potential therapeutic targets for muscle-related diseases and novel strategies for livestock improvement.

Dietary Short-Chain Fatty Acids Supplementation Improves Reproductive Performance and Gut Microbiota in Gilts

BACKGROUND

Short-chain fatty acids (SCFAs) have emerged as critical modulators of female reproductive function and host gut microbiota.

OBJECTIVES

This study aimed to investigate the impact of dietary SCFAs supplementation on reproductive performance and gut microbiota in gilts, and to elucidate the underlying mechanisms.

METHODS

Eighty gilts (95 d old) were randomly assigned to either a control group (Ctrl, 40 gilts) receiving a basal diet, or a SCFAs treatment group (SCFAs, 40 gilts) receiving a basal diet supplemented with 0.13% sodium acetate, 0.11% sodium propionate, and 0.09% sodium butyrate. At third estrus, 13 gilts (6 from Ctrl and 7 from SCFAs) were killed for follicular development and gut microbiota analysis, whereas the remaining gilts completed gestation for reproductive performance assessment.

RESULTS

SCFAs group had higher total number born (11.97 compared with 9.44) and total number born alive (11.28 compared with 9.34) compared with Ctrl group (P < 0.05). SCFAs group had increased counts of secondary follicles (36.14 compared with 26.83), antral follicles (10.29 compared with 6.67), and corpus luteum (25.09 compared with 19.33), alongside had reduced atretic follicles (15.32 compared with 20.67) compared with Ctrl group (P < 0.05). Proteomic analysis revealed that SCFAs-induced differentially expressed proteins (DEPs) were significantly enriched in the follicular development-related pathways (P < 0.05). Apoptosis-related DEPs positively correlated with follicular development indices (P < 0.05), consistent with the reduced apoptosis observed in ovarian granulosa cells of the SCFAs group. Additionally, SCFAs supplementation improved both the composition and alpha-diversity (P < 0.05) of gilts' gut microbiota. Furthermore, both the SCFAs-enriched bacteria and plasma SCFAs concentrations showed positive associations with gilts' follicular development indices (P < 0.05).

CONCLUSIONS

Dietary SCFAs supplementation enhances reproductive performance in gilts by promoting ovarian follicular maturation and optimizing gut microbiota composition.

Chromosome-Level Genome Assembly of the Meishan Pig and Insights into Its Domestication Mechanisms

Pigs are essential agricultural animals, and among the various breeds, the Meishan pig, a native breed of China, is renowned for its high reproductive performance. This breed has been introduced to many countries to enhance local pig breeding programs. However, there have been limited genomic and population genetics studies focusing on Meishan pigs. We created a chromosomal-level genomic assembly using high-depth PacBio sequencing and Illumina sequencing data collected from a Meishan pig. Additionally, we analyzed whole-genome sequencing (WGS) data from Chinese boars and Meishan pigs to identify domestication selection signals within the Meishan breed. The assembled genome of the Meishan pig (MSjxau) was found to be 2.45 Gb in size, with a scaffold length of 139.17 Mb. The quality value was 37.06, and the BUSCO score was 96.2%, indicating good completeness, continuity, and accuracy. We annotated transposable elements, segmental duplication, and genes in the MSjxau genome. By combining these data with 28 publicly available genomes, we provide a high-quality structural variants resource for pigs. Furthermore, we identified 716 selective sweep intervals between Chinese wild pigs and Meishan pigs, where the selected gene PGR may be linked to the high fertility observed in Meishan pigs. Our study offers valuable genomic and variation resources for pig breeding and identifies several genes associated with the domestication of the Meishan pig. This lays the groundwork for further investigation into the genetic mechanisms behind complex traits in pigs.

Muscle Characteristics Comparison Analysis Reveal Differences in the Meat Quality and Nutritional Components of Three Shanghai Local Pig Breeds

To study the differences in meat quality and nutritional components between the local Shanghai pig breeds Meishan pig (MS), Shawutou pig (SWT), Fengjing pig (FJ), and the commercial Duroc × Landrace × Yorkshire (DLY) crossbred pigs, and to provide data support for the selection and breeding of superior pig breeds, this study selected 30 piglets each of three local pig breeds and DLY with similar birth ages and weights around 25 kg, fed them the same daily ration with uniform nutritional components, and slaughtered ten of them at around 100 kg weight for evaluation of differences in meat quality indicators (primarily intramuscular fat content, tenderness value, texture, etc.) and amino acid content among the varieties. The results indicated significant differences among the four pig breeds in intramuscular fat content, with MS having the highest content and significant differences in tenderness value compared to the other three breeds (p < 0.05). In terms of texture indicators, MS and SWT differed significantly from FJ and DLY in terms of hardness and fracturability indicators (p < 0.05), with certain differences existing among the four breeds in other indicators. In amino acid content, the differences in total amino acid content among the three Shanghai local pig breeds were not significant (p > 0.05), but all were significantly higher than DLY. Further analysis revealed significant differences in amino acid content between Shanghai local pig breeds and DLY, with Shanghai local pigs showing markedly higher levels of serine, proline, isoleucine, leucine, and histidine compared to DLY (p < 0.05). Regarding nucleotides, the cytidine monophosphate (CMP) indicator of MS differed significantly from the other three breeds (p < 0.05), SWT's uridine monophosphate (UMP) indicator differed significantly from FJ and DLY, and FJ and DLY's inosine monophosphate (IMP) indicator was significantly higher than MS and SWT (p < 0.05), while SWT's adenosine monophosphate (AMP) indicator was significantly higher than the other three breeds (p < 0.05). The results of this study suggest that the meat quality and nutritional composition of Shanghai local pigs are significantly superior to DLY, with MS exhibiting significantly better meat quality and nutrition compared to SWT and FJ among the three local pig breeds.

Genome Wide Identification of Structure Variations in Five Italian Turkey Populations

Structural variants (SVs) are one of the main sources of genetic variants and have a significant impact on phenotype evolution, disease susceptibility, and environmental adaptations. We used 73 whole genome sequencing (12x) to apply a mapping approach to identify SVs in five turkey populations. A notable degree of genetic isolation was observed between the Basilicata and Apulian populations, as indicated by principal component analysis and admixture results. A total of 11,733 SVs were detected, including 6712 deletions, 2671 duplications, 1430 inversions, and 920 translocations. The Variant Effect Predictor (VEP) analysis predicted various consequences of filtered SVs as follows: intron variants (35.8%), intergenic variants (9.6%), coding sequence variants (8.3%), downstream gene variants (7.5%), and transcript ablations (7.3%). Our functional annotation of genes overlapping with SVs was mainly enriched in recognized pathways governing positive regulation of nucleoplasm, protein binding, mitochondrion, negative regulation of cell population proliferation, identical protein binding, and calcium signaling. We produced a comprehensive SV catalog utilizing unique whole-genome turkey data. This SV catalog not only increases our understanding of genetic diversity in turkeys but also enhances our knowledge of the role of SVs in their phenotypic traits.

Pangenome Reveals Gene Content Variations and Structural Variants Contributing to Pig Characteristics

Pigs are one of the most essential sources of high-quality proteins in human diets. Structural variants (SVs) are a major source of genetic variants associated with diverse traits and evolutionary events. However, the current linear reference genome of pigs restricts the accurate presentation of position information for SVs. In this study, we generated a pangenome of pigs and a genome variation map of 599 deeply sequenced genomes across Eurasia. Additionally, we established a section-wide gene repertoire, revealing that core genes are more evolutionarily conserved than variable genes. Furthermore, we identified 546,137 SVs, their enrichment regions, and relationships with genomic features and found significant divergence across Eurasian pigs. More importantly, the pangenome-detected SVs could complement heritability estimates and genome-wide association studies based only on single nucleotide polymorphisms. Among the SVs shaped by selection, we identified an insertion in the promoter region of the TBX19 gene, which may be related to the development, growth, and timidity traits of Asian pigs and may affect the gene expression. The constructed pig pangenome and the identified SVs in this study provide rich resources for future functional genomic research on pigs.

A near telomere-to-telomere genome assembly of the Jinhua pig: enabling more accurate genetic research

BACKGROUND

Pigs are crucial sources of meat and protein, valuable animal models, and potential donors for xenotransplantation. However, the existing reference genome for pigs is incomplete, with thousands of segments and centromeres and telomeres missing, which limits our understanding of the important traits in these genomic regions.

FINDINGS

We present a near-complete genome assembly for the Jinhua pig (JH-T2T) and provide a set of diploid Jinhua reference genomes, constructed using PacBio HiFi, ONT long reads, and Hi-C reads. This assembly includes all 18 autosomes and the X and Y sex chromosomes, with only 6 gaps. It features annotations of 46.90% repetitive sequences, 33 telomeres, 17 centromeres, and 23,924 high-confident genes. Compared to the Sscrofa11.1, JH-T2T closes nearly all gaps, extends sequences by 177 Mb, predicts more intact telomeres and centromeres, and gains 799 more genes and loses 114 genes. Moreover, it enhances the mapping rate for both Western and Chinese local pigs, outperforming Sscrofa11.1 as a reference genome. Additionally, this comprehensive genome assembly will facilitate large-scale variant detection.

CONCLUSIONS

This study produced a near-gapless assembly of the pig genome and provides a set of haploid Jinhua reference genomes. Our findings represent a significant advance in pig genomics, providing a robust resource that enhances genetic research, breeding programs, and biomedical applications.

Suppressing IGF2R mitigates hypoxia-induced apoptosis by reducing the expression of pro-apoptotic factor BAX

Oxidative stress caused by hypoxia can lead to serious bodily damage and functional degradation. Our previous study in pigs showed that the insulin-like growth factor II receptor (IGF2R) gene might participate in the process of hypoxia adaptability. To investigate the function and mechanism of IGF2R in cellular hypoxia tolerance, we analyze the effect of IGF2R on cell survival capacity under hypoxia conditions in intestinal porcine enterocyte cell line (IPEC-J2) cells. The results show that under hypoxia condition (3% O2), cell viability is significantly reduced, the expression of IGF2R and cell apoptosis are significantly increased. Functional analysis suggests that suppressing IGF2R expression under hypoxia does not affect cell cycle and cell proliferation but increases cellular viability. Meanwhile, the expression of the pro-apoptotic gene BAX is reduced, the hypoxia-induced apoptosis is rescued, and cell survival is significantly improved. Transcriptome analysis suggests that global gene expression changes in knockdown IGF2R under hypoxia, IGF2R may regulate apoptosis through oxidative phosphorylation. Our findings demonstrate that suppressing IGF2R expression under hypoxia can rescue hypoxia-induced cell injury by reducing the expression of BAX, highlighting the potential ability of IGF2R regulation for the treatment of hypoxia stress.

Chromosome-level genome assembly of Huai pig (Sus scrofa)

Although advances in long-read sequencing technology and genome assembly techniques have facilitated the study of genomes, little is known about the genomes of unique Chinese indigenous breeds, including the Huai pig. Huai pig is an ancient domestic pig breed and is well-documented for its redder meat color and high forage tolerance compared to European domestic pigs. In the present study, we sequenced and assembled the Huai pig genome using PacBio, Hi-C, and Illumina sequencing technologies. The final highly contiguous chromosome-level Huai pig genome spans 2.53 Gb with a scaffold N50 of 138.92 Mb. The Benchmarking Universal Single-Copy Orthologs (BUSCO) completeness score for the assembled genome was 95.33%. Remarkably, 23,389 protein-coding genes were annotated in the Huai-pig genome, along with 45.87% repetitive sequences. Overall, this study provided new foundational resources for future genetic research on Chinese domestic pigs.

Genome-wide detection of structural variation in some sheep breeds using whole-genome long-read sequencing data

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep.

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.

A chromosome-level genome of Chenghua pig provides new insights into the domestication and local adaptation of pigs

As a country with abundant genetic resources of pigs, the domestication history of pigs in China and the adaptive evolution of Chinese pig breeds at different latitudes have rarely been elucidated at the genome-wide level. To fill this gap, we first assembled a high-quality chromosome-level genome of the Chenghua pig and used it as a benchmark to analyse the genomes of 272 samples from three genera of three continents. The divergence of the three species belonging to three genera, Phacochoerus africanus, Potamochoerus porcus, and Sus scrofa, was assessed. The introgression of pig breeds redefined that the migration routes were basically from southern China to central and southwestern China, then spread to eastern China, arrived in northern China, and finally reached Europe. The domestication of pigs in China occurred ∼12,000 years ago, earlier than the available Chinese archaeological domestication evidence. In addition, FBN1 and NR6A1 were identified in our study as candidate genes related to extreme skin thickness differences in Eurasian pig breeds and adaptive evolution at different latitudes in Chinese pig breeds, respectively. Our study provides a new resource for the pig genomic pool and refines our understanding of pig genetic diversity, domestication, migration, and adaptive evolution at different latitudes.

Mapping and functional characterization of structural variation in 1060 pig genomes

BACKGROUND

Structural variations (SVs) have significant impacts on complex phenotypes by rearranging large amounts of DNA sequence.

RESULTS

We present a comprehensive SV catalog based on the whole-genome sequence of 1060 pigs (Sus scrofa) representing 101 breeds, covering 9.6% of the pig genome. This catalog includes 42,487 deletions, 37,913 mobile element insertions, 3308 duplications, 1664 inversions, and 45,184 break ends. Estimates of breed ancestry and hybridization using genotyped SVs align well with those from single nucleotide polymorphisms. Geographically stratified deletions are observed, along with known duplications of the KIT gene, responsible for white coat color in European pigs. Additionally, we identify a recent SINE element insertion in MYO5A transcripts of European pigs, potentially influencing alternative splicing patterns and coat color alterations. Furthermore, a Yorkshire-specific copy number gain within ABCG2 is found, impacting chromatin interactions and gene expression across multiple tissues over a stretch of genomic region of ~200 kb. Preliminary investigations into SV's impact on gene expression and traits using the Pig Genotype-Tissue Expression (PigGTEx) data reveal SV associations with regulatory variants and gene-trait pairs. For instance, a 51-bp deletion is linked to the lead eQTL of the lipid metabolism regulating gene FADS3, whose expression in embryo may affect loin muscle area, as revealed by our transcriptome-wide association studies.

CONCLUSIONS

This SV catalog serves as a valuable resource for studying diversity, evolutionary history, and functional shaping of the pig genome by processes like domestication, trait-based breeding, and adaptive evolution.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Assembly of novel sequences for Chinese domestic pigs reveals new genes and regulatory variants providing new insights into their diversity

There is an increasing understanding that a reference genome representing an individual cannot capture all the gene repertoire of a species. Here, we conduct a population-scale missing sequences detection of Chinese domestic pigs using whole-genome sequencing data from 534 individuals. We identify 132.41 Mb of sequences absent in the reference assembly, including eight novel genes. In particular, the breeds spread in Chinese high-altitude regions perform significantly different frequencies of new sequences in promoters than other breeds. Furthermore, we dissect the role of non-coding variants and identify a novel sequence inserted in the 3'UTR of the FMO3 gene, which may be associated with the intramuscular fat phenotype. This novel sequence could be a candidate marker for meat quality. Our study provides a comprehensive overview of the missing sequences in Chinese domestic pigs and indicates that this dataset is a valuable resource for understanding the diversity and biology of pigs.

Pig pangenome graph reveals functional features of non-reference sequences

BACKGROUND

The reliance on a solitary linear reference genome has imposed a significant constraint on our comprehensive understanding of genetic variation in animals. This constraint is particularly pronounced for non-reference sequences (NRSs), which have not been extensively studied.

RESULTS

In this study, we constructed a pig pangenome graph using 21 pig assemblies and identified 23,831 NRSs with a total length of 105 Mb. Our findings revealed that NRSs were more prevalent in breeds exhibiting greater genetic divergence from the reference genome. Furthermore, we observed that NRSs were rarely found within coding sequences, while NRS insertions were enriched in immune-related Gene Ontology terms. Notably, our investigation also unveiled a close association between novel genes and the immune capacity of pigs. We observed substantial differences in terms of frequencies of NRSs between Eastern and Western pigs, and the heat-resistant pigs exhibited a substantial number of NRS insertions in an 11.6 Mb interval on chromosome X. Additionally, we discovered a 665 bp insertion in the fourth intron of the TNFRSF19 gene that may be associated with the ability of heat tolerance in Southern Chinese pigs.

CONCLUSIONS

Our findings demonstrate the potential of a graph genome approach to reveal important functional features of NRSs in pig populations.

Identification of a Novel Long Non-Coding RNA G8110 That Modulates Porcine Adipogenic Differentiation and Inflammatory Responses

Long non-coding RNAs (lncRNAs) have been extensively studied, and their crucial roles in adipogenesis, lipid metabolism, and gene expression have been revealed. However, the exact regulatory or other mechanisms by which lncRNAs influence the functioning of mesenteric adipose tissue (MAT) remain largely unknown. In this paper, we report the identification of a new lncRNA, named G8110, from the MAT of Bama pigs. The coordinated expression levels of lncRNA G8110 and NFE2L1 were significantly decreased in the MAT of obese Bama pigs compared with those in the MAT of lean pigs. Using a bone mesenchymal stem cell adipogenic differentiation model, we found that lncRNA G8110 played a role in adipocyte differentiation by positively regulating NFE2L1. We also found that lncRNA G8110 inhibited the formation of intracellular lipid synthesis, promoted lipid metabolism, and inhibited the expression of inflammatory cytokines. Our findings regarding lipid synthesis may further promote the role of lncRNAs in driving adipose tissue remodeling and maintaining metabolic health.

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.

Molecular Mechanism of MYL4 Regulation of Skeletal Muscle Development in Pigs

The processes of muscle growth and development, including myoblast proliferation, migration, differentiation, and fusion, are modified by a variety of regulatory factors. MYL4 plays an important role in atrial development, atrial cardiomyopathy, muscle-fiber size, and muscle development. The structural variation (SV) of MYL4 was found via the de novo sequencing of Ningxiang pigs, and the existence of SV was verified in the experiments. The genotype distribution of Ningxiang pigs and Large White pigs was detected, and it was found that Ningxiang pigs were mainly of the BB genotype and that Large White pigs were mainly of the AB genotype. However, the molecular mechanisms behind the MYL4-mediated regulation of skeletal muscle development need to be deeply explored. Therefore, RT-qPCR, 3'RACE, CCK8, EdU, Western blot, immunofluorescence, flow cytometry, and bioinformation analysis were used to explore the function of MYL4 in myoblast development. The cDNA of MYL4 was successfully cloned from Ningxiang pigs, and its physicochemical properties were predicted. The expression profiles in six tissues and four stages of Ningxiang pigs and Large White pigs were found to be the highest in the lungs and 30 days after birth. The expression of MYL4 increased gradually with the extension of the myogenic differentiation time. The myoblast function test showed that the overexpression of MYL4 inhibited proliferation and promoted apoptosis and differentiation. The knockdown of MYL4 showed the opposite result. These results enhance our understanding of the molecular mechanisms of muscle development and provide a solid theoretical foundation for further exploring the role of the MYL4 gene in muscle development.

Phosphoproteomic analysis identifies differentially expressed phosphorylation sites that affect muscle fiber type in pigs

Skeletal muscle of livestock is composed of both fast- and slow-twitch muscle fibers, which are key factors in their meat quality. However, the role of protein phosphorylation in muscle fiber type is not completely understood. Here, a fast-twitch (biceps femoris, BF) and slow-twitch (soleus, SOL) muscle tissue sample was collected from three male offspring of Duroc and Meishan pigs. We demonstrate that the meat quality of SOL muscle is significantly better than that of BF muscle. We further used phosphoproteomic profiling of BF and SOL muscles to identify differences between these muscle types. A total of 2,327 phosphorylation sites from 770 phosphoproteins were identified. Among these sites, 287 differentially expressed phosphorylation sites (DEPSs) were identified between BF and SOL. GO and KEGG enrichment analysis of proteins containing DEPSs showed that these phosphorylated proteins were enriched in the glycolytic process GO term and the AMPK signaling pathway. A protein-protein interaction (PPI) analysis reveals that these phosphorylated proteins interact with each other to regulate the transformation of muscle fiber type. These analyses reveal that protein phosphorylation modifications are involved in porcine skeletal muscle fiber type transformation. This study provides new insights into the molecular mechanisms by which protein phosphorylation regulates muscle fiber type transformation and meat quality in pigs.

Long-range interaction within the chromatin domain determines regulatory patterns in porcine skeletal muscle

Spatial chromatin structure is crucial for understanding the early growth and development of porcine skeletal muscle. However, its characteristic of 3D architecture and elaborate regulation of gene transcription remains unclear. In this study, ChIA-PET method is used to study the changes of early chromatin three-dimensional structure in skeletal muscle of lean type Yorkshire pig and fat type Meishan pig. Integrating the in situ Hi-C data revealed the 3D architecture and long-range interaction of the porcine muscle. The results showed the CTCF/RNAPII mediated long-range interaction shapes the different chromatin architecture and dominates the unique regulation of enhancers. In addition, the results revealed that key myogenic genes like ssc-mir-1 had a unique enhancer regulation function in myogenesis. Interestingly, the FGF6 gene is of breed-specific regulation, implying the difference between two breeds in skeletal muscle development. Our research thus may provide a clue for the porcine genetic improvement of skeletal muscle.

Large Fragment InDels Reshape Genome Structure of Porcine Alveolar Macrophage 3D4/21 Cells

The porcine monomyeloid cell line, or 3D4/21 cells, is an effective tool to study the immune characteristics and virus infection mechanism of pigs. Due to the introduction of the neomycin resistance gene and the SV40 large T antigen gene, its genome has undergone essential changes, which are still unknown. Studying the variation in genome structure, especially the large fragments of insertions and deletions (InDels), is one of the proper ways to reveal these issues. In this study, an All-seq method was established by combining Mate-pair and Shotgun sequencing methods, and the detection and verification of large fragments of InDels were performed on 3D4/21 cells. The results showed that there were 844 InDels with a length of more than 1 kb, of which 12 regions were deletions of more than 100 kb in the 3D4/21 cell genome. In addition, compared with porcine primary alveolar macrophages, 82 genes including the CD163 had lost transcription in 3D4/21 cells, and 72 genes gained transcription as well. Further referring to the Hi-C structure, it was found that the fusion of the topologically associated domains (TADs) caused by the deletion may lead to abnormal gene function. The results of this study provide a basis for elaborating the genome structure and functional variation in 3D4/21 cells, provide a method for rapid and convenient detection of large-scale InDels, and provide useful clues for the study of the porcine immune function genome and the molecular mechanism of virus infection.

Genomic Analysis of Porcine Reproductive and Respiratory Syndrome Virus 1 Revealed Extensive Recombination and Potential Introduction Events in China

Simple Summary

Porcine reproductive and respiratory syndrome, caused by the porcine reproductive and respiratory syndrome virus, is considered one of the most devastating swine diseases worldwide. Porcine reproductive and respiratory syndrome virus 1 was first isolated in China in 2006, and there have been few reports concerning its genetic characteristics in China. We hope to find out the regularity of genetic diversity, recombination, and evolution of the virus by analyzing all available genomic sequences during 1991–2018. We found that high-frequency recombination regions were concentrated in non-structural protein 2 and structural proteins 2 to 4 and extensive deletions in non-structural protein 2; phylogenetic analysis revealed four independent introductions in China. Our results suggest that attention should be paid to the prevention and control of porcine reproductive and respiratory syndrome virus 1 and the rational use of vaccine strains. These results will help us to understand the recombination of porcine reproductive and respiratory syndrome virus and strengthen viral inspection before mixing herds of swine to reduce the probability of novel recombinant variants. Moreover, our study might form the basis of monitoring and control measures to prevent the spread of this economically important virus.

Abstract

Porcine reproductive and respiratory syndrome (PRRS), caused by the PRRS virus (PRRSV), is considered one of the most devastating swine diseases worldwide. PRRSV-1 was first isolated in China in 2006. However, there were few reports concerning the genetic characteristics of PRRSV-1 in China. In this study, three PRRSV-1 strains (HL85, HeB3, and HeB47) were detected by a general RT-qPCR method from clinical samples in 2018. HeB47 was identified as a recombinant between the BJEU06-1 and CReSA228-like strains. To further analyze the recombination and deletion features of PRRSV-1, all the available 88 complete genome sequences (isolated in 19 countries) from 1991 to 2018 in GenBank were analyzed. The high-frequency recombination regions were concentrated in NSP2 and GP2 to GP4. More importantly, phylogenetic analysis of PRRSV-1 revealed four independent introductions in China. Therefore, it is necessary to strengthen the important monitoring of breeding pigs and pork products and epidemiological surveys on pig farms to prevent the further spread of PRRSV-1.

Origin, evolution, and tissue-specific functions of the porcine repetitive element 1

Background

The porcine repetitive element 1 (PRE1) is the most abundant short interspersed nuclear element (SINE) in the Sus scrofa genome and it has been suggested that some PRE1 can have regulatory functions. The million copies of PRE1 in the porcine genome have accumulated abundant CpG dinucleotides and unique structural variations, such as direct repeats and patterns of sequence degeneration. The aims of this study were to analyse these structural variations to trace the origin and evolutionary pattern of PRE1 and to investigate potential methylation-related functions of PRE1 based on methylation patterns of PRE1 CpG dinucleotides in different tissues.

Results

We investigated the evolutionary trajectory of PRE1 and found that PRE1 originated from the ancestral CHRS-S1 family through three main successive partial duplications. We found that the partial duplications and deletions of PRE1 were likely due to RNA splicing events during retrotransposition. Functionally, correlation analysis showed that the methylation levels of 103 and 261 proximal PRE1 were, respectively, negatively and positively correlated with the expression levels of neighboring genes (Spearman correlation, P < 0.01). Further epigenomic analysis revealed that, in the testis, demethylation of proximal PRE1 in the HORMAD1 and HACD3 genes had tissue-specific enhancer and promoter functions, while in the muscle, methylation of proximal PRE1 repeats in the TCEA3 gene had an enhancer function.

Conclusions

The characteristic sequences of PRE1 reflect unique patterns of origin and evolution and provide a structural basis for diverse regulatory functions.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-022-00745-3.

Initial Analysis of Structural Variation Detections in Cattle Using Long-Read Sequencing Methods

Structural variations (SVs), as a great source of genetic variation, are widely distributed in the genome. SVs involve longer genomic sequences and potentially have stronger effects than SNPs, but they are not well captured by short-read sequencing owing to their size and relevance to repeats. Improved characterization of SVs can provide more advanced insight into complex traits. With the availability of long-read sequencing, it has become feasible to uncover the full range of SVs. Here, we sequenced one cattle individual using 10× Genomics (10 × G) linked read, Pacific Biosciences (PacBio) continuous long reads (CLR) and circular consensus sequencing (CCS), as well as Oxford Nanopore Technologies (ONT) PromethION. We evaluated the ability of various methods for SV detection. We identified 21,164 SVs, which amount to 186 Mb covering 7.07% of the whole genome. The number of SVs inferred from long-read-based inferences was greater than that from short reads. The PacBio CLR identified the most of large SVs and covered the most genomes. SVs called with PacBio CCS and ONT data showed high uniformity. The one with the most overlap with the results obtained by short-read data was PB CCS. Together, we found that long reads outperformed short reads in terms of SV detections.

The de novo assembly of a European wild boar genome revealed unique patterns of chromosomal structural variations and segmental duplications

The rapid progress of sequencing technology has greatly facilitated the de novo genome assembly of pig breeds. However, the assembly of the wild boar genome is still lacking, hampering our understanding of chromosomal and genomic evolution during domestication from wild boars into domestic pigs. Here, we sequenced and de novo assembled a European wild boar genome (ASM2165605v1) using the long‐range information provided by 10× Linked‐Reads sequencing. We achieved a high‐quality assembly with contig N50 of 26.09 Mb. Additionally, 1.64% of the contigs (222) with lengths from 107.65 kb to 75.36 Mb covered 90.3% of the total genome size of ASM2165605v1 (~2.5 Gb). Mapping analysis revealed that the contigs can fill 24.73% (93/376) of the gaps present in the orthologous regions of the updated pig reference genome (Sscrofa11.1). We further improved the contigs into chromosome level with a reference‐assistant scaffolding method. Using the ‘assembly‐to‐assembly’ approach, we identified intra‐chromosomal large structural variations (SVs, length >1 kb) between ASM2165605v1 and Sscrofa11.1 assemblies. Interestingly, we found that the number of SV events on the X chromosome deviated significantly from the linear models fitting autosomes (R² > 0.64, p < 0.001). Specifically, deletions and insertions were deficient on the X chromosome by 66.14 and 58.41% respectively, whereas duplications and inversions were excessive on the X chromosome by 71.96 and 107.61% respectively. We further used the large segmental duplications (SDs, >1 kb) events as a proxy to understand the large‐scale inter‐chromosomal evolution, by resolving parental‐derived relationships for SD pairs. We revealed a significant excess of SD movements from the X chromosome to autosomes (p < 0.001), consistent with the expectation of meiotic sex chromosome inactivation. Enrichment analyses indicated that the genes within derived SD copies on autosomes were significantly related to biological processes involving nervous system, lipid biosynthesis and sperm motility (p < 0.01). Together, our analyses of the de novo assembly of ASM2165605v1 provides insight into the SVs between European wild boar and domestic pig, in addition to the ongoing process of meiotic sex chromosome inactivation in driving inter‐chromosomal interaction between the sex chromosome and autosomes.