Reorganization of chromatin architecture during prenatal development of porcine skeletal muscle

Fibin is a crucial mitochondrial regulatory gene in skeletal muscle development

Fin bud initiation factor homolog (Fibin) is a secreted protein that is relatively conserved among species. It is closely related to fin bud development and can regulate a variety of cellular processes. In our previous high-throughput chromosome conformation capture (Hi-C) study of pig embryonic muscle development, it was found that Fibin has high expression and activity during the development of pig primary muscle fibers. Therefore, we speculated Fibin participated in myogenesis severely. Specific deletion of Fibin in mouse skeletal muscle resulted in abnormal primary muscle fiber development during the embryonic period and a substantial decrease in skeletal muscle mass in adulthood. In vitro, knocking out Fibin in C2C12 cells promoted cell proliferation; however, after myogenic induction, cells lacking Fibin had almost no ability to differentiate into myotubes. During myogenic differentiation, loss of Fibin disrupts the normal function of mitochondria and impairs oxidative phosphorylation, resulting in decrease of NADH and FADH in the electron transport chain. Transmission electron microscopy also showed that mitochondrial morphology of Fibin-deficient C2C12 was impaired. In conclusion, our research has unveiled a novel mechanism of myogenesis regulation in mitochondrial function and potential target Fibin, and improved understanding of a broad range of mitochondrial muscle diseases.

Comprehensive genome-wide analysis of genetic loci and candidate genes associated with litter traits in purebred Berkshire pigs of Korea

OBJECTIVE

The objective of this study was to identify genomic regions and candidate genes associated with the total number of piglets born (TNB), number of piglets born alive (NBA), and total number of stillbirths (TNS) in Berkshire pigs.

METHODS

This study used a total of 11,228 records and 2,843 single-nucleotide polymorphism (SNP) data obtained from Illumina porcine 60 K and 80 K chips. The estimated genomic breeding values (GEBVs) and SNP effects were estimated using weighted single-step genomic BLUP (WssGBLUP).

RESULTS

The heritabilities of the TNB, NBA, and TNS were determined using single-step genomic best linear unbiased prediction (ssGBLUP). The heritability estimates were 0.13, 0.12, and 0.015 for TNB, NBA, and TNS, respectively. When comparing the accuracy of breeding value estimates, the results using pedigree-based BLUP (PBLUP) were 0.58, 0.60, and 0.31 for TNB, NBA, and TNS, respectively. In contrast, the accuracy increased to 0.67, 0.66, and 0.42 for TNB, NBA, and TNS, respectively, when using WssGBLUP, specifically in the last three iterations. The results of weighted single-step genome-wide association studies (WssGWAS) showed that the highest variance explained for each trait was predominantly located in the Sus scrofa chromosome 5 (SSC5) region. Specifically, the variance exceeded 4% for TNB, 3% for NBA, and 6% for TNS. Within the SSC5 region (12.26 to 12.76 Mb), which exhibited the highest variance for TNB, 20 SNPs were identified, and five candidate genes were identified: TIMP3, SYN3, FBXO7, BPIFC, and RTCB.

CONCLUSION

The identified SNP markers for TNB, NBA, and TNS were expected to provide valuable information for genetic improvement as an understanding of their expression and genetic architecture in Berkshire pigs. With the accumulation of more phenotype and SNP data in the future, it is anticipated that more effective SNP markers will be identified.

Dynamic alternations of three-dimensional chromatin architecture contribute to phenotypic characteristics of breast muscle in chicken

Breast muscle growth rate and intramuscular fat (IMF) content show apparent differences between fast-growing broilers and slow-growing indigenous chickens. However, the underlying genetic basis of these phenotypic characteristics remains elusive. In this study, we investigate the dynamic alterations of three-dimensional genome architecture and chromatin accessibility in breast muscle across four key developmental stages from embryo to starter chick in Arbor Acres (AA) broilers and Yufen (YF) indigenous chickens. The limited breed-specifically up-regulated genes (Bup-DEGs) are embedded in breed-specific A compartment, while a majority of the Bup-DEGs involving myogenesis and adipogenesis are regulated by the breed-specific TAD reprogramming. Chromatin loops allow distal accessible regions to interact with myogenic genes, and those loops share an extremely low similarity between chicken with different growth rate. Moreover, AA-specific loop interactions promote the expression of 40 Bup-DEGs, such as IGF1, which contributes to myofiber hypertrophy. YF-specific loop interactions or distal accessible regions lead to increased expression of 5 Bup-DEGs, including PIGO, PEMT, DHCR7, TMEM38B, and DHDH, which contribute to IMF deposition. These results help elucidate the regulation of breast muscle growth and IMF deposition in chickens.

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.

Novel Genetic and Biochemical Insights into the Spectrum of NEFL-Associated Phenotypes

Background

NEFL encodes for the neurofilament light chain protein. Pathogenic variants in NEFL cause demyelinating, axonal and intermediate forms of Charcot-Marie-Tooth disease (CMT) which present with a varying degree of severity and somatic mutations have not been described yet. Currently, 34 different CMT-causing pathogenic variants in NEFL in 174 patients have been reported. Muscular involvement was also described in CMT2E patients mostly as a secondary effect. Also, there are a few descriptions of a primary muscle vulnerability upon pathogenic NEFL variants.

Objectives

To expand the current knowledge on the genetic landscape, clinical presentation and muscle involvement in NEFL-related neurological diseases by retrospective case study and literature review.

Methods

We applied in-depth phenotyping of new and already reported cases, molecular genetic testing, light-, electron- and Coherent Anti-Stokes Raman Scattering-microscopic studies and proteomic profiling in addition to in silico modelling of NEFL-variants.

Results

We report on a boy with a muscular phenotype (weakness, myalgia and cramps, Z-band alterations and mini-cores in some myofibers) associated with the heterozygous p.(Phe104Val) NEFL-variant, which was previously described in a neuropathy case. Skeletal muscle proteomics findings indicated affection of cytoskeletal proteins. Moreover, we report on two further neuropathic patients (16 years old girl and her father) both carrying the heterozygous p.(Pro8Ser) variant, which has been identified as 15% somatic mosaic in the father. While the daughter presented with altered neurophysiology,neurogenic clump feet and gait disturbances, the father showed clinically only feet deformities. As missense variants affecting proline at amino acid position 8 are leading to neuropathic manifestations of different severities, in silico modelling of these different amino acid substitutions indicated variable pathogenic impact correlating with disease onset.

Conclusions

Our findings provide new morphological and biochemical insights into the vulnerability of denervated muscle (upon NEFL-associated neuropathy) as well as novel genetic findings expanding the current knowledge on NEFL-related neuromuscular phenotypes and their clinical manifestations. Along this line, our data show that even subtle expression of somatic NEFL variants can lead to neuromuscular symptoms.

Advances in the discovery of genetic elements underlying longissimus dorsi muscle growth and development in the pig

As a major source of protein in human diets, pig meat plays a crucial role in ensuring global food security. Key determinants of meat production refer to the chemical and physical compositions or characteristics of muscle fibers, such as the number, hypertrophy potential, fiber-type conversion and intramuscular fat deposition. However, the growth and formation of muscle fibers comprises a complex process under spatio-temporal regulation, that is, the intermingled and concomitant proliferation, differentiation, migration and fusion of myoblasts. Recently, with the fast and continuous development of next-generation sequencing technology, the integration of quantitative trait loci mapping with genome-wide association studies (GWAS) has greatly helped animal geneticists to discover and explore thousands of functional or causal genetic elements underlying muscle growth and development. However, owing to the underlying complex molecular mechanisms, challenges to in-depth understanding and utilization remain, and the cost of large-scale sequencing, which requires integrated analyses of high-throughput omics data, is high. In this review, we mainly elaborate on research advances in integrative analyses (e.g. GWAS, omics) for identifying functional genes or genomic elements for longissimus dorsi muscle growth and development for different pig breeds, describing several successful transcriptome analyses and functional genomics cases, in an attempt to provide some perspective on the future functional annotation of genetic elements for muscle growth and development in pigs.

An Analysis of a Transposable Element Expression Atlas during 27 Developmental Stages in Porcine Skeletal Muscle: Unveiling Molecular Insights into Pork Production Traits

The development and growth of porcine skeletal muscle determine pork quality and yield. While genetic regulation of porcine skeletal muscle development has been extensively studied using various omics data, the role of transposable elements (TEs) in this context has been less explored. To bridge this gap, we constructed a comprehensive atlas of TE expression throughout the developmental stages of porcine skeletal muscle. This was achieved by integrating porcine TE genomic coordinates with whole-transcriptome RNA-Seq data from 27 developmental stages. We discovered that in pig skeletal muscle, active Tes are closely associated with active epigenomic marks, including low levels of DNA methylation, high levels of chromatin accessibility, and active histone modifications. Moreover, these TEs include 6074 self-expressed TEs that are significantly enriched in terms of muscle cell development and myofibril assembly. Using the TE expression data, we conducted a weighted gene co-expression network analysis (WGCNA) and identified a module that is significantly associated with muscle tissue development as well as genome-wide association studies (GWAS) of the signals of pig meat and carcass traits. Within this module, we constructed a TE-mediated gene regulatory network by adopting a unique multi-omics integration approach. This network highlighted several established candidate genes associated with muscle-relevant traits, including HES6, CHRNG, ACTC1, CHRND, MAMSTR, and PER2, as well as novel genes like ENSSSCG00000005518, ENSSSCG00000033601, and PIEZO2. These novel genes hold promise for regulating muscle-related traits in pigs. In summary, our research not only enhances the TE-centered dissection of the genetic basis underlying pork production traits, but also offers a general approach for constructing TE-mediated regulatory networks to study complex traits or diseases.

White-throated sparrow (Zonotrichia albicollis) liver and pectoralis flight muscle transcriptomic changes in preparation for migration

Migratory songbirds undertake challenging journeys to reach their breeding grounds each spring. They accomplish these nonstop flapping feats of endurance through a suite of physiological changes, including the development of substantial fat stores and flight muscle hypertrophy and an increased capacity for fat catabolism. In addition, migratory birds may show large reductions in organ masses during flight, including the flight muscle and liver, which they must rapidly rebuild during their migratory stopover before replenishing their fat stores. However, the molecular basis of this capacity for rapid tissue remodeling and energetic output has not been thoroughly investigated. We performed RNA-sequencing analysis of the liver and pectoralis flight muscle of captive white-throated sparrows in experimentally photostimulated migratory and nonmigratory condition to explore the mechanisms of seasonal change to metabolism and tissue mass regulation that may facilitate these migratory journeys. Based on transcriptional changes, we propose that tissue-specific adjustments in preparation for migration may alleviate the damaging effects of long-duration activity, including a potential increase to the inflammatory response in the muscle. Furthermore, we hypothesize that seasonal hypertrophy balances satellite cell recruitment and apoptosis, while little evidence appeared in the transcriptome to support myostatin-, insulin-like growth factor 1-, and mammalian target of rapamycin-mediated pathways for muscle growth. These findings can encourage more targeted molecular studies on the unique integration of pathways that we find in the development of the migratory endurance phenotype in songbirds.NEW & NOTEWORTHY Migratory songbirds undergo significant physiological changes to accomplish their impressive migratory journeys. However, we have a limited understanding of the regulatory mechanisms underlying these changes. Here, we explore the transcriptomic changes to the flight muscle and liver of white-throated sparrows as they develop the migratory condition. We use these patterns to develop hypotheses about metabolic flexibility and tissue restructuring in preparation for migration.

Transcriptome analysis of mRNAs, lncRNAs, and miRNAs in the skeletal muscle of Tibetan chickens at different developmental stages

Introduction: As a valuable genetic resource, native birds can contribute to the sustainable development of animal production. Tibetan chickens, known for their special flavor, are one of the important local poultry breeds in the Qinghai-Tibet Plateau. However, Tibetan chickens have a slow growth rate and poor carcass traits compared with broilers. Although most of the research on Tibetan chickens focused on their hypoxic adaptation, there were fewer studies related to skeletal muscle development. Methods: Here, we performed the transcriptional sequencing of leg muscles from Tibetan chicken embryos at E (embryonic)10, E14, and E18. Results: In total, 1,600, 4,610, and 2,166 DE (differentially expressed) mRNAs, 210, 573, and 234 DE lncRNAs (long non-coding RNAs), and 52, 137, and 33 DE miRNAs (microRNAs) were detected between E10 and E14, E10 and E18, and E14 and E18, respectively. Functional prediction showed several DE mRNAs and the target mRNAs of DE lncRNAs and DE miRNAs were significantly enriched in sarcomere organization, actin cytoskeleton organization, myofibril, muscle fiber development, and other terms and pathways related to muscle growth and development. Finally, a lncRNA-miRNA-mRNA ceRNA (competing endogenous RNA) network associated with muscle growth and development, which contained 6 DE lncRNAs, 13 DE miRNAs, and 50 DE mRNAs, was constructed based on the screened DE RNAs by Gene Ontology (GO) enrichment. These DE RNAs may play a critical regulatory role in the skeletal muscle development of chickens. Discussion: The results provide a genomic resource for mRNAs, lncRNAs, and miRNAs potentially involved in the skeletal muscle development of chickens, which lay the foundation for further studies of the molecular mechanisms underlying skeletal muscle growth and development in Tibetan chickens.

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

BACKGROUND

N6-methyladenosine (m⁶A) 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 m⁶A 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 m⁶A 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 m⁶A writers and co-occurrence of 5mC and m⁶A modification were revealed from public datasets of C2C12 myoblasts. The protein-protein interactions between the DNA methylase and the m⁶A 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/m⁶A 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 m⁶A and gene expression profiles revealed a genome/transcriptome-wide correlated dynamics and co-occurrence of 5mC and m⁶A modifications as a consequence of 5mC/m⁶A crosstalk in the postnatal myogenesis progress of pigs. Last, we identified a group of myogenesis-related genes collaboratively regulated by both 5mC and m⁶A 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.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Loss of Monoallelic Expression of IGF2 in the Adult Liver Via Alternative Promoter Usage and Chromatin Reorganization

In mammals, genomic imprinting operates via gene silencing mechanisms. Although conservation of the imprinting mechanism at the H19/IGF2 locus has been generally described in pigs, tissue-specific imprinting at the transcript level, monoallelic-to-biallelic conversion, and spatio-temporal chromatin reorganization remain largely uninvestigated. Here, we delineate spatially regulated imprinting of IGF2 transcripts, age-dependent hepatic mono- to biallelic conversion, and reorganization of topologically associating domains at the porcine H19/IGF2 locus for better translation to human and animal research. Whole-genome bisulfite sequencing (WGBS) and RNA sequencing (RNA-seq) of normal and parthenogenetic porcine embryos revealed the paternally hypermethylated H19 differentially methylated region and paternal expression of IGF2. Using a polymorphism-based approach and omics datasets from chromatin immunoprecipitation sequencing (ChIP–seq), whole-genome sequencing (WGS), RNA-seq, and Hi-C, regulation of IGF2 during development was analyzed. Regulatory elements in the liver were distinguished from those in the muscle where the porcine IGF2 transcript was monoallelically expressed. The IGF2 transcript from the liver was biallelically expressed at later developmental stages in both pigs and humans. Chromatin interaction was less frequent in the adult liver compared to the fetal liver and skeletal muscle. The duration of genomic imprinting effects within the H19/IGF2 locus might be reduced in the liver with biallelic conversion through alternative promoter usage and chromatin remodeling. Our integrative omics analyses of genome, epigenome, and transcriptome provided a comprehensive view of imprinting status at the H19/IGF2 cluster.

Reorganization of 3D genome architecture across wild boar and Bama pig adipose tissues

BACKGROUND

A growing body of evidence has revealed that the mammalian genome is organized into hierarchical layers that are closely correlated with and may even be causally linked with variations in gene expression. Recent studies have characterized chromatin organization in various porcine tissues and cell types and compared them among species and during the early development of pigs. However, how chromatin organization differs among pig breeds is poorly understood.

RESULTS

In this study, we investigated the 3D genome organization and performed transcriptome characterization of two adipose depots (upper layer of backfat [ULB] and greater omentum [GOM]) in wild boars and Bama pigs; the latter is a typical indigenous pig in China. We found that over 95% of the A/B compartments and topologically associating domains (TADs) are stable between wild boars and Bama pigs. In contrast, more than 70% of promoter-enhancer interactions (PEIs) are dynamic and widespread, involving over a thousand genes. Alterations in chromatin structure are associated with changes in the expression of genes that are involved in widespread biological functions such as basic cellular functions, endocrine function, energy metabolism and the immune response. Approximately 95% and 97% of the genes associated with reorganized A/B compartments and PEIs in the two pig breeds differed between GOM and ULB, respectively.

CONCLUSIONS

We reported 3D genome organization in adipose depots from different pig breeds. In a comparison of Bama pigs and wild boar, large-scale compartments and TADs were mostly conserved, while fine-scale PEIs were extensively reorganized. The chromatin architecture in these two pig breeds was reorganized in an adipose depot-specific manner. These results contribute to determining the regulatory mechanism of phenotypic differences between Bama pigs and wild boar.