Global DNA methylation pattern involved in the modulation of differentiation potential of adipogenic and myogenic precursors in skeletal muscle of pigs

Integrated analysis of the DNA methylome and RNA transcriptome during the development of skeletal muscle in Duroc pigs

BACKGROUND

Skeletal muscle development plays a crucial role in yield and quality of pork; however, this process is influenced by various factors. In this study, we employed whole-genome bisulfite sequencing (WGBS) and transcriptome sequencing to comprehensively investigate the longissimus dorsi muscle (LDM), aiming to identify key genes that impact the growth and development of Duroc pigs with different average daily gains (ADGs).

RESULTS

Eight pigs were selected and divided into two groups based on ADGs: H (774.89 g) group and L (658.77 g) group. Each pair of the H and L groups were half-siblings. The results of methylation sequencing revealed 2631 differentially methylated genes (DMGs) involved in metabolic processes, signalling, insulin secretion, and other biological activities. Furthermore, a joint analysis was conducted on these DMGs and the differentially expressed genes (DEGs) obtained from transcriptome sequencing of the same individual. This analysis identified 316 differentially methylated and differentially expressed genes (DMEGs), including 18 DMEGs in promoter regions and 294 DMEGs in gene body regions. Finally, LPAR1 and MEF2C were selected as candidate genes associated with muscle development. Bisulfite sequencing PCR (BSP) and quantitative real-time PCR (qRT-PCR) revealed that the promoter region of LPAR1 exhibited significantly lower methylation levels (P < 0.05) and greater expression levels (P < 0.05) in the H group than in the L group. Additionally, hypermethylation was observed in the gene body region of MEF2C, as was a low expression level, in the H group (P < 0.05).

CONCLUSIONS

These results suggest that the differences in the ADGs of Duroc pigs fed the same diet may be influenced by the methylation levels and expression levels of genes related to skeletal muscle development.

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.

Genomic structure analysis and construction of DNA fingerprint for four sheep populations

The Yongdeng Qishan sheep (QS) is a sheep population found locally in China. To gain in-depth knowledge of its population characteristics, three control groups were chosen, comprising the Lanzhou fat-tailed sheep (LFT), TAN sheep (TAN), and Minxian black fur sheep (MBF), inhabiting the nearby environments. This study genotyped a total of 120 individuals from four sheep populations: QS, LFT, TAN, and MBF. Using Specific-Locus Amplified Fragment Sequencing, we conducted genetic diversity, population structure, and selective sweep analysis, and constructed the fingerprint of each population. In total, there were 782 535 single nucleotide polymorphism (SNP) variations identified, with most being situated within regions that are intergenic or intronic. The genetic diversity analysis revealed that the QS population exhibited lower genetic diversity compared to the other three populations. Consistent results were obtained from the principal component, phylogenetic tree, and population structure analysis, indicating significant genetic differences between QS and the other three populations. However, a certain degree of differentiation was observed within the QS population. The linkage disequilibrium (LD) patterns among the four populations showed clear distinctions, with the QS group demonstrating the most rapid LD decline. Kinship analysis supported the findings of population structure, dividing the 90 QS individuals into two subgroups consisting of 23 and 67 individuals. Selective sweep analysis identified a range of genes associated with reproduction, immunity, and adaptation to high-altitude hypoxia. These genes hold potential as candidate genes for marker-assisted selection breeding. Additionally, a total of 86 523 runs of homozygosity (ROHs) were detected, showing non-uniform distribution across chromosomes, with chromosome 1 having the highest coverage percentage and chromosome 26 the lowest. In the high-frequency ROH islands, 79 candidate genes were associated with biological processes such as reproduction and fat digestion and absorption. Furthermore, a DNA fingerprint was constructed for the four populations using 349 highly polymorphic SNPs. In summary, our research delves into the genetic diversity and population structure of QS population. The construction of DNA fingerprint profiles for each population can provide valuable references for the identification of sheep breeds both domestically and internationally.

The Integration of Genome-Wide DNA Methylation and Transcriptomics Identifies the Potential Genes That Regulate the Development of Skeletal Muscles in Ducks

DNA methylation is a pivotal epigenetic regulatory mechanism in the development of skeletal muscles. Nonetheless, the regulators responsible for DNA methylation in the development of embryonic duck skeletal muscles remain unknown. In the present study, whole genome bisulfite sequencing (WGBS) and transcriptome sequencing were conducted on the skeletal muscles of embryonic day 21 (E21) and day 28 (E28) ducks. The DNA methylation pattern was found to fall mainly within the cytosine-guanine (CG) context, with high methylation levels in the intron, exon, and promoter regions. Overall, 7902 differentially methylated regions (DMRs) were identified, which corresponded to 3174 differentially methylated genes (DMGs). By using integrative analysis of both WGBS with transcriptomics, we identified 1072 genes that are DMGs that are negatively associated with differentially expressed genes (DEGs). The gene ontology (GO) analysis revealed significant enrichment in phosphorylation, kinase activity, phosphotransferase activity, alcohol-based receptors, and binding to cytoskeletal proteins. The Kyoto Encyclopedia of Genes and Genomes (KEGGs) analysis showed significant enrichment in MAPK signaling, Wnt signaling, apelin signaling, insulin signaling, and FoxO signaling. The screening of enriched genes showed that hyper-methylation inhibited the expression of Idh3a, Got1, Bcl2, Mylk2, Klf2, Erbin, and Klhl38, and hypo-methylation stimulated the expression of Col22a1, Dnmt3b, Fn1, E2f1, Rprm, and Wfikkn1. Further predictions showed that the CpG islands in the promoters of Klhl38, Klf2, Erbin, Mylk2, and Got1 may play a crucial role in regulating the development of skeletal muscles. This study provides new insights into the epigenetic regulation of the development of duck skeletal muscles.

Identification of a 5-Methylcytosine Site (mC-7) That May Inhibit CXCL11 Expression and Regulate E. coli F18 Susceptibility in IPEC-J2 Cells

The primary pathogen causing post-weaning diarrhea in piglets is Escherichia coli F18 (E. coli F18), hence it is essential to investigate the mechanism governing E. coli F18 resistance in native pig breeds. Based on the previous RNA-seq results of the duodenum from E. coli F18-resistant and -susceptible Meishan piglets, CXCL11, an important functional gene, was preliminarily screened. In this investigation, in order to further examine the expression regulation mechanism of E. coli F18 in intestinal porcine epithelial cells (IPEC-J2) against E. coli F18 infection, CXCL11 gene expression on IPEC-J2 cells infected by E. coli F18 was detected, which was significantly downregulated (p < 0.01). Secondly, the overexpression on the IPEC-J2 cell line was successfully structured, and a relative quantification method of the PILIN, bacteria enumeration, and immunofluorescence assay indicated that the CXCL11 overexpression significantly reduced the ability of E. coli F18 to interact with IPEC-J2 in vitro. The promoter region of the CXCL11 gene was predicted to contain a CpG island (−619 ~ −380 bp) of which 13 CpG sites in the sequencing region were methylated to varying degrees, and the methylation level of one CPG site (mC-7) positively linked negatively with the expression of the CXCL11 gene (p < 0.05). Meanwhile, a dual luciferase assay detected the mutation of the mC-7 site that significantly inhibited the luciferase activity of the CXCL11 gene promoter (p < 0.01). Transcription factor prediction and expression verification indicated that mC-7 is located in the OSR1-binding domain, and that its expression level is related to E. coli F18 susceptibility. We speculated that methylation modification of the mC-7 site of the CpG island in the promoter region of the CXCL11 gene might inhibit the binding of transcription factor OSR1 with the mC-7 site, and then affect its expression level to regulate the susceptibility to E. coli F18.

Integration of ATAC-seq and RNA-seq analysis identifies key genes affecting intramuscular fat content in pigs

Meat quality is one of the most important economic traits in pig breeding and production, and intramuscular fat (IMF) content is the major factor in improving meat quality. The IMF deposition in pigs is influenced by transcriptional regulation, which is dependent on chromatin accessibility. However, how chromatin accessibility plays a regulatory role in IMF deposition in pigs has not been reported. Xidu black is a composite pig breed with excellent meat quality, which is an ideal research object of this study. In this study, we used the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) analysis to identify the accessible chromatin regions and key genes affecting IMF content in Xidu black pig breed with extremely high and low IMF content. First, we identified 21,960 differential accessible chromatin peaks and 297 differentially expressed genes. The motif analysis of differential peaks revealed several potential cis-regulatory elements containing binding sites for transcription factors with potential roles in fat deposition, including Mef2c, CEBP, Fra1, and AP-1. Then, by integrating the ATAC-seq and RNA-seq analysis results, we found 47 genes in the extremely high IMF (IMF_H) group compared with the extremely low IMF (IMF_L) group. For these genes, we observed a significant positive correlation between the differential gene expression and differential ATAC-seq signal (r² = 0.42). This suggests a causative relationship between chromatin remodeling and the resulting gene expression. We identified several candidate genes (PVALB, THRSP, HOXA9, EEPD1, HOXA10, and PDE4B) that might be associated with fat deposition. Through the PPI analysis, we found that PVALB gene was the top hub gene. In addition, some pathways that might regulate fat cell differentiation and lipid metabolism, such as the PI3K-Akt signaling pathway, MAPK signaling pathway, and calcium signaling pathway, were significantly enriched in the ATAC-seq and RNA-seq analysis. To the best of our knowledge, our study is the first to use ATAC-seq and RNA-seq to examine the mechanism of IMF deposition from a new perspective. Our results provide valuable information for understanding the regulation mechanism of IMF deposition and an important foundation for improving the quality of pork.

Exosome loaded hydroxyapatite (HA) scaffold promotes bone regeneration in calvarial defect: an in vivo study

In this study, hydroxyapatite (HA) scaffolds were synthesized and characterized, following the osteogenic and angiogenic effects of HA scaffolds with or without endometrial mesenchymal stem stromal cells (hEnSCs) derived Exosomes were investigated in rat animal model with calvaria defect. The X-ray diffraction (XRD) analysis of HA powder formation was confirmed with Joint Corporation of Powder Diffraction Standards (JCPDS) files numbers of 34-0010 and 24-0033A and Ball mill, and sintering manufactured Nano-size particles. Obtained results containing FE-SEM images presented that the surface of scaffolds has a rough and porous structure, which makes them ideal and appropriate for tissue engineering. Additionally, the XRD showed that these scaffolds exhibited a crystallized structure without undergoing phase transformation; meanwhile, manufactured scaffolds consistently release exosomes; moreover, in vivo findings containing hematoxylin-eosin staining, immunohistochemistry, Masson's trichrome staining, and histomorphometric analysis confirmed that our implant has an osteogenic and angiogenic characteristic. So prepared scaffolds containing exosomes can be proposed as a promising substitute in tissue engineering.

circRNAome profiling reveals circFgfr2 regulates myogenesis and muscle regeneration via a feedback loop

BACKGROUND

Circular RNAs (circRNAs) represent a novel class of non-coding RNAs formed by a covalently closed loop and play crucial roles in many biological processes. Several circRNAs associated with myogenesis have been reported. However, the dynamic expression, function, and mechanism of circRNAs during myogenesis and skeletal muscle development are largely unknown.

METHODS

Strand-specific RNA-sequencing (RNA-seq) and microarray datasets were used to profile the dynamic circRNAome landscape during skeletal muscle development and myogenic differentiation. Bioinformatics analyses were used to characterize the circRNAome and identify candidate circRNAs associated with myogenesis. Bulk and single-cell RNA-seq were performed to identify the downstream genes and pathways of circFgfr2. The primary myoblast cells, C2C12 cells, and animal model were used to assess the function and mechanism of circFgfr2 in myogenesis and muscle regeneration in vitro or in vivo by RT-qPCR, western blotting, dual-luciferase activity assay, RNA immunoprecipitation, RNA fluorescence in situ hybridization, and chromatin immunoprecipitation.

RESULTS

We profiled the dynamic circRNAome in pig skeletal muscle across 27 developmental stages and detected 52 918 high-confidence circRNAs. A total of 2916 of these circRNAs are conserved across human, mouse, and pig, including four circRNAs (circFgfr2, circQrich1, circMettl9, and circCamta1) that were differentially expressed (|log₂ fold change| > 1 and adjusted P value < 0.05) in various myogenesis systems. We further focused on a conserved circRNA produced from the fibroblast growth factor receptor 2 (Fgfr2) gene, termed circFgfr2, which was found to inhibit myoblast proliferation and promote differentiation and skeletal muscle regeneration. Mechanistically, circFgfr2 acted as a sponge for miR-133 to regulate the mitogen-activated protein kinase kinase kinase 20 (Map3k20) gene and JNK/MAPK pathway. Importantly, transcription factor Kruppel like factor 4 (Klf4), the downstream target of the JNK/MAPK pathway, directly bound to the promoter of circFgfr2 and affected its expression via an miR-133/Map3k20/JNK/Klf4 auto-regulatory feedback loop. RNA binding protein G3BP stress granule assembly factor 1 (G3bp1) inhibited the biogenesis of circFgfr2.

CONCLUSIONS

The present study provides a comprehensive circRNA resource for skeletal muscle study. The functional and mechanistic analysis of circFgfr2 uncovered a circRNA-mediated auto-regulatory feedback loop regulating myogenesis and muscle regeneration, which provides new insight to further understand the regulatory mechanism of circRNAs.

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.