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

Transcriptomic regulatory analysis of skeletal muscle development in landrace pigs

The development of pig skeletal muscle is a complex dynamic regulation process, which mainly includes the formation of primary and secondary muscle fibers, the remodeling of muscle fibers, and the maturation of skeletal muscle; However, the regulatory mechanism of the entire developmental process remains unclear. This study analyzed the whole-transcriptome data of skeletal muscles at 27 developmental nodes (E33-D180) in Landrace pigs, and their key regulatory factors in the development process were identified using the bioinformatics method. Firstly, we constructed a transcriptome expression map of skeletal muscle development from embryo to adulthood in Landrace pig. Subsequently, due to drastic change in gene expression, the perinatal periods including E105, D0 and D9, were focused, and the genes related to the process of muscle fiber remodeling and volume expansion were revealed. Then, though conjoint analysis with miRNA and lncRNA transcripts, a ceRNA network were identified, which consist of 11 key regulatory genes (such as CHAC1, RTN4IP1 and SESN1), 7 miRNAs and 43 lncRNAs, and they potentially play an important role in the process of muscle fiber differentiation, muscle fiber remodeling and volume expansion, intramuscular fat deposition, and other skeletal muscle developmental events. In summary, we reveal candidate genes and underlying molecular regulatory networks associated with perinatal skeletal muscle fiber type remodeling and expansion. These data provide new insights into the molecular regulation of mammalian skeletal muscle development and diversity.

Genome-wide characteristics and potential functions of circular RNAs from the embryo muscle development in Chengkou mountain chicken

The Chengkou mountain chicken, a native Chinese poultry breed, holds significant importance in the country's poultry sector due to its delectable meat and robust stress tolerance. Muscle growth and development are pivotal characteristics in poultry breeding, with muscle fiber development during the embryonic period crucial for determining inherent muscle growth potential. Extensive evidence indicates that non-coding RNAs (ncRNAs) play a regulatory role in muscle growth and development. Among ncRNAs, circular RNAs (circRNAs), characterized by a closed-loop structure, have been shown to modulate biological processes through the regulation of microRNAs (miRNAs). This study seeks to identify and characterize the spatiotemporal-specific expression of circRNAs during embryonic muscle development in Chengkou mountain chicken, and to construct the potential regulatory network of circRNAs-miRNA-mRNAs. The muscle fibers of HE-stained sections became more distinct, and their boundaries were more defined over time. Subsequent RNA sequencing of 12 samples from four periods generated 9,904 novel circRNAs, including 917 differentially expressed circRNAs. The weighted gene co-expression network analysis (WGCNA)-identified circRNA source genes significantly enriched pathways related to cell fraction, cell growth, and muscle fiber growth regulation. Furthermore, a competitive endogenous RNA (ceRNA) network constructed using combined data of present and previous differentially expressed circRNAs, miRNA, and mRNA revealed that several circRNA transcripts regulate MYH1D, MYH1B, CAPZA1, and PERM1 proteins. These findings provide insight into the potential pathways and mechanisms through which circRNAs regulate embryonic muscle development in poultry, a theoretical support for trait improvement in domestic chickens.

Advancing skeletal health and disease research with single-cell RNA sequencing

Orthopedic conditions have emerged as global health concerns, impacting approximately 1.7 billion individuals worldwide. However, the limited understanding of the underlying pathological processes at the cellular and molecular level has hindered the development of comprehensive treatment options for these disorders. The advent of single-cell RNA sequencing (scRNA-seq) technology has revolutionized biomedical research by enabling detailed examination of cellular and molecular diversity. Nevertheless, investigating mechanisms at the single-cell level in highly mineralized skeletal tissue poses technical challenges. In this comprehensive review, we present a streamlined approach to obtaining high-quality single cells from skeletal tissue and provide an overview of existing scRNA-seq technologies employed in skeletal studies along with practical bioinformatic analysis pipelines. By utilizing these methodologies, crucial insights into the developmental dynamics, maintenance of homeostasis, and pathological processes involved in spine, joint, bone, muscle, and tendon disorders have been uncovered. Specifically focusing on the joint diseases of degenerative disc disease, osteoarthritis, and rheumatoid arthritis using scRNA-seq has provided novel insights and a more nuanced comprehension. These findings have paved the way for discovering novel therapeutic targets that offer potential benefits to patients suffering from diverse skeletal disorders.

Divergent Regulatory Roles of Transcriptional Variants of the Chicken LDB3 Gene in Muscle Shaping

LIM domain binding 3 (LDB3) serves as a striated muscle-specific Z-band alternatively spliced protein that plays an important role in mammalian skeletal muscle development, but its regulatory role and molecular mechanism in avian muscle development are still unclear. In this study, we reanalyzed RNA sequencing data sets of 1415 samples from 21 chicken tissues published in the NCBI GEO database. First, three variants (LDB3-X, LDB3-XN1, and LDB3-XN2) generated by alternative splicing of the LDB3 gene were identified in chicken skeletal muscle, among which LDB3-XN1 and LDB3-XN2 are novel variants. LDB3-X and LDB3-XN1 are derived from exon skipping in chicken skeletal muscle at the E18-D7 stage and share three LIM domains, but LDB3-XN2 lacks a LIM domain. Our results preliminarily suggest that the formation of three variants of LDB3 is regulated by RBM20. The three splice isomers have divergent functions in skeletal muscle according to in vitro and in vivo assays. Finally, we identified the mechanism by which different variants play different roles through interactions with IGF2BP1 and MYHC, which promote the proliferation and differentiation of chicken myoblasts, in turn regulating chicken myogenesis. In conclusion, this study revealed the divergent roles of three LDB3 variants in chicken myogenesis and muscle remodeling and demonstrated their regulatory mechanism through protein-protein interactions.

Biomaterials-Based Technologies in Skeletal Muscle Tissue Engineering

For many clinically prevalent severe injuries, the inherent regenerative capacity of skeletal muscle remains inadequate. Skeletal muscle tissue engineering (SMTE) seeks to meet this clinical demand. With continuous progress in biomedicine and related technologies including micro/nanotechnology and 3D printing, numerous studies have uncovered various intrinsic mechanisms regulating skeletal muscle regeneration and developed tailored biomaterial systems based on these understandings. Here, the skeletal muscle structure and regeneration process are discussed and the diverse biomaterial systems derived from various technologies are explored in detail. Biomaterials serve not merely as local niches for cell growth, but also as scaffolds endowed with structural or physicochemical properties that provide tissue regenerative cues such as topographical, electrical, and mechanical signals. They can also act as delivery systems for stem cells and bioactive molecules that have been shown as key participants in endogenous repair cascades. To achieve bench-to-bedside translation, the typical effect enabled by biomaterial systems and the potential underlying molecular mechanisms are also summarized. Insights into the roles of biomaterials in SMTE from cellular and molecular perspectives are provided. Finally, perspectives on the advancement of SMTE are provided, for which gene therapy, exosomes, and hybrid biomaterials may hold promise to make important contributions.

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.

Tankyrase-1 regulates RBP-mediated mRNA turnover to promote muscle fiber formation

Poly(ADP-ribosylation) (PARylation) is a post-translational modification mediated by a subset of ADP-ribosyl transferases (ARTs). Although PARylation-inhibition based therapies are considered as an avenue to combat debilitating diseases such as cancer and myopathies, the role of this modification in physiological processes such as cell differentiation remains unclear. Here, we show that Tankyrase1 (TNKS1), a PARylating ART, plays a major role in myogenesis, a vital process known to drive muscle fiber formation and regeneration. Although all bona fide PARPs are expressed in muscle cells, experiments using siRNA-mediated knockdown or pharmacological inhibition show that TNKS1 is the enzyme responsible of catalyzing PARylation during myogenesis. Via this activity, TNKS1 controls the turnover of mRNAs encoding myogenic regulatory factors such as nucleophosmin (NPM) and myogenin. TNKS1 mediates these effects by targeting RNA-binding proteins such as Human Antigen R (HuR). HuR harbors a conserved TNKS-binding motif (TBM), the mutation of which not only prevents the association of HuR with TNKS1 and its PARylation, but also precludes HuR from regulating the turnover of NPM and myogenin mRNAs as well as from promoting myogenesis. Therefore, our data uncover a new role for TNKS1 as a key modulator of RBP-mediated post-transcriptional events required for vital processes such as myogenesis.

Long non-coding RNAs and their role in muscle regeneration

In mammals, most of the genome is transcribed to generate a large and heterogeneous variety of non-protein coding RNAs, that are broadly grouped according to their size. Long noncoding RNAs include a very large and versatile group of molecules. Despite only a minority of them has been functionally characterized, there is emerging evidence indicating long noncoding RNAs as important regulators of expression at multiple levels. Several of them have been shown to be modulated during myogenic differentiation, playing important roles in the regulation of skeletal muscle development, differentiation and homeostasis, and contributing to neuromuscular diseases. In this chapter, we have summarized the current knowledge about long noncoding RNAs in skeletal muscle and discussed specific examples of long noncoding RNAs (lncRNAs and circRNAs) regulating muscle stem cell biology. We have also discussed selected long noncoding RNAs involved in the most common neuromuscular diseases.

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.

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.

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.

Circular RNA circIGF2BP3 Promotes the Proliferation and Differentiation of Chicken Primary Myoblasts

The quality and quantity of animal meat are closely related to the development of skeletal muscle, which, in turn, is determined by myogenic cells, including myoblasts and skeletal muscle satellite cells (SMSCs). Circular RNA, an endogenous RNA derivative formed through specific reverse splicing in mRNA precursors, has the potential to influence muscle development by binding to miRNAs or regulating gene expression involved in muscular growth at the transcriptional level. Previous high-throughput sequencing of circRNA in chicken liver tissue revealed a circular transcript, circIGF2BP3, derived from the gene encoding insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3). In this study, we confirmed the presence of the natural circular molecule of circIGF2BP3 through an RNase R enzyme tolerance assay. RT-qPCR results showed high circIGF2BP3 expression in the pectoral and thigh muscles of Yuexi frizzled feather chickens at embryonic ages 14 and 18, as well as at 7 weeks post-hatch. Notably, its expression increased during embryonic development, followed by a rapid decrease after birth. As well as using RT-qPCR, Edu, CCK-8, immunofluorescence, and Western blot techniques, we demonstrated that overexpressing circIGF2BP3 could promote the proliferation and differentiation of chicken primary myoblasts through upregulating genes such as proliferating cell nuclear antigen (PCNA), cyclin D1 (CCND1), cyclin E1 (CCNE1), cyclin dependent kinase 2 (CDK2), myosin heavy chain (MyHC), myoblast-determining 1 (MyoD1), myogenin (MyoG), and Myomaker. In conclusion, circIGF2BP3 promotes the proliferation and differentiation of myoblasts in chickens. This study establishes a foundation for further investigation into the biological functions and mechanisms of circIGF2BP3 in myoblasts proliferation and differentiation.

Evolutionary conserved circular MEF2A RNAs regulate myogenic differentiation and skeletal muscle development

Circular RNAs (circRNAs) have been recognized as critical regulators of skeletal muscle development. Myocyte enhancer factor 2A (MEF2A) is an evolutionarily conserved transcriptional factor that regulates myogenesis. However, it remains unclear whether MEF2A produces functional circRNAs. In this study, we identified two evolutionarily conserved circular MEF2A RNAs (circMEF2As), namely circMEF2A1 and circMEF2A2, in chicken and mouse muscle stem cells. Our findings revealed that circMEF2A1 promotes myogenesis by regulating the miR-30a-3p/PPP3CA/NFATC1 axis, whereas circMEF2A2 facilitates myogenic differentiation by targeting the miR-148a-5p/SLIT3/ROBO2/β-catenin signaling pathway. Furthermore, in vivo experiments demonstrated that circMEF2As both promote skeletal muscle growth. We also discovered that the linear MEF2A mRNA-derived MEF2A protein binds to its own promoter region, accelerating the transcription of MEF2A and upregulating the expression of both linear MEF2A and circMEF2As, forming a MEF2A autoregulated positive feedback loop. Moreover, circMEF2As positively regulate the expression of linear MEF2A by adsorbing miR-30a-3p and miR-148a-5p, which directly contribute to the MEF2A autoregulated feedback loop. Importantly, we found that mouse circMEF2As are essential for the myogenic differentiation of C2C12 cells. Collectively, our results demonstrated the evolution, function, and underlying mechanisms of circMEF2As in animal myogenesis, which may provide novel insight for both the farm animal meat industry and human medicine.

Genomic Scan for Runs of Homozygosity and Selective Signature Analysis to Identify Candidate Genes in Large White Pigs

Large White pigs are extensively utilized in China for their remarkable characteristics of rapid growth and the high proportion of lean meat. The economic traits of pigs, comprising reproductive and meat quality traits, play a vital role in swine production. In this study, 2295 individuals, representing three different genetic backgrounds Large White pig populations were used: 500 from the Canadian line, 295 from the Danish line, and 1500 from the American line. The GeneSeek 50K GGP porcine HD array was employed to genotype the three pig populations. Firstly, genomic selective signature regions were identified using the pairwise fixation index (FST) and locus-specific branch length (LSBL). By applying a top 1% threshold for both parameters, a total of 888 candidate selective windows were identified, harbouring 1571 genes. Secondly, the investigation of regions of homozygosity (ROH) was performed utilizing the PLINK software. In total, 25 genomic regions exhibiting a high frequency of ROHs were detected, leading to the identification of 1216 genes. Finally, the identified potential functional genes from candidate genomic regions were annotated, and several important candidate genes associated with reproductive traits (ADCYAP1, U2, U6, CETN1, Thoc1, Usp14, GREB1L, FGF12) and meat quality traits (MiR-133, PLEKHO1, LPIN2, SHANK2, FLVCR1, MYL4, SFRP1, miR-486, MYH3, STYX) were identified. The findings of this study provide valuable insights into the genetic basis of economic traits in Large White pigs and may have potential use in future pig breeding programs.

Genome-wide analysis of circular RNA-mediated ceRNA regulation in porcine skeletal muscle development

BACKGROUND

As a diverse and abundant class of endogenous RNAs, circular RNAs (circRNAs) participate in various biological processes including cell proliferation and apoptosis. Nevertheless, few researchers have investigated the role of circRNAs in muscle development in cultivated pigs.

RESULTS

In this study, we used RNA-seq to construct circRNA expression profiles in skeletal muscle of Jinfen White pigs at the age of 1, 90, and 180 days. Among the 16,990 identified circRNAs, 584 circRNAs were differentially expressed. Moreover, the enrichment analysis of DE circRNA host genes showed that they were mainly involved in muscle contraction, muscle organ development and muscle system processes, as well as AMPK and cAMP-related signal pathways. We also constructed a circRNA-miRNA-mRNA co-expression network to find key circRNAs which many involved in the regulation of porcine skeletal muscle development through the competitive endogenous RNA (ceRNA) mechanism. It is noteworthy that circ_0018595/miR-1343/PGM1 axis may play a regulatory role in the development of porcine skeletal muscle.

CONCLUSIONS

This study identified the circRNAs and present the circRNA expression profile in the development of pigs, revealed that DE circRNA host genes participate in different cell fates and enriched the porcine ceRNA network. Thus, this work will become a valuable resource for further in-depth study of the regulatory mechanism of circRNA in the development of porcine skeletal muscle.

CircPLXNA2 Affects the Proliferation and Apoptosis of Myoblast through circPLXNA2/gga-miR-12207-5P/MDM4 Axis

CircRNAs are newly identified special endogenous RNA molecules that covalently close a loop by back-splicing with pre-mRNA. In the cytoplasm, circRNAs would act as molecular sponges to bind with specific miRNA to promote the expression of target genes. However, knowledge of circRNA functional alternation in skeletal myogenesis is still in its infancy. In this study, we identified a circRNA–miRNA–mRNA interaction network in which the axis may be implicated in the progression of chicken primary myoblasts’ (CPMs) myogenesis by multi-omics (i.e., circRNA-seq and ribo-seq). In total, 314 circRNA–miRNA–mRNA regulatory axes containing 66 circRNAs, 70 miRNAs, and 24 mRNAs that may be relevant to myogenesis were collected. With these, the circPLXNA2-gga-miR-12207-5P-MDM4 axis aroused our research interest. The circPLXNA2 is highly differentially expressed during differentiation versus proliferation. It was demonstrated that circPLXNA2 inhibited the process of apoptosis while at the same time stimulating cell proliferation. Furthermore, we demonstrated that circPLXNA2 could inhibit the repression of gga-miR-12207-5p to MDM4 by directing binding to gga-miR-12207-5p, thereby restoring MDM4 expression. In conclusion, circPLXNA2 could function as a competing endogenous RNA (ceRNA) to recover the function of MDM4 by directing binding to gga-miR-12207-5p, thereby regulating the myogenesis.

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.

Integrative transcriptomics analysis and experimental validation reveal immunomodulatory patterns in keratoconus

Keratoconus is a progressive disorder of the cornea and is typically considered a noninflammatory disease. However, increasing evidence indicates that immune disorders play an essential role in keratoconus progression, but the immune-related etiology remains elusive. Here, we comprehensively utilized bioinformatics approaches and experimental methods to explore the potential immunoregulatory mechanism of keratoconus progression. Transcriptomics data containing two keratoconus patient groups was derived from the public dataset GSE151631. The intersection of genes and known immunological genes was used to obtain differentially expressed immune-related genes. We utilized various protein clustering algorithms to screen out and validated the hub immune-related genes, and further explored their potential biological functions via gene annotation and pathway enrichment analyses. Moreover, the underlying immune landscape and drug targets were predicted by immune cell infiltration analysis and drug-gene interaction analysis. Furthermore, keratoconus-related immunoregulatory competitive endogenous RNA networks were constructed and experimentally validated. After filtering and experimental validation, nine keratoconus-associated immune-related genes were credible. Infiltrated monocytes might play an essential role in the progression of keratoconus. Moreover, eleven intersecting drugs targeting four genes, CCR2, CCR5, F2RL1, and ADORA1, were considered as potential druggable molecular targets for keratoconus. Furthermore, in the competitive endogenous RNA network, we identified several lncRNAs and miRNAs as critical noncoding RNAs regulating the hub genes. Overall, our data indicated that the immunomodulatory patterns had undergone changes in the pathogenesis of keratoconus, which might facilitate the understanding of keratoconus-related immune processes and provide novel insights into developing new immunotherapies for keratoconus.

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.

Pyruvate dehydrogenase B regulates myogenic differentiation via the FoxP1-Arih2 axis

BACKGROUND

Sarcopenia, the age-related decline in skeletal muscle mass and function, diminishes life quality in elderly people. Improving the capacity of skeletal muscle differentiation is expected to counteract sarcopenia. However, the mechanisms underlying skeletal muscle differentiation are complex, and effective therapeutic targets are largely unknown.

METHODS

The human Gene Expression Omnibus database, aged mice and primary skeletal muscle cells were used to assess the expression level of pyruvate dehydrogenase B (PDHB) in human and mouse aged state. d-Galactose (d-gal)-induced sarcopenia mouse model and two classic cell models (C2C12 and HSkMC) were used to assess the myogenic effect of PDHB and the underlying mechanisms via immunocytochemistry, western blotting, quantitative real-time polymerase chain reaction, RNA interference or overexpression, dual-luciferase reporter assay, RNA sequencing and untargeted metabolomics.

RESULTS

We identified that a novel target PDHB promoted myogenic differentiation. PDHB expression decreased in aged mouse muscle relative to the young state (-50% of mRNA level, P < 0.01) and increased during mouse and primary human muscle cell differentiation (+3.97-fold, P < 0.001 and +3.79-fold, P < 0.001). Knockdown or overexpression of PDHB modulated the expression of genes related to muscle differentiation, namely, myogenic factor 5 (Myf5) (-46%, P < 0.01 and -27%, P < 0.05; +1.8-fold, P < 0.01), myogenic differentiation (MyoD) (-55%, P < 0.001 and -34%, P < 0.01; +2.27-fold, P < 0.001), myogenin (MyoG) (-60%, P < 0.001 and -70%, P < 0.001; +5.46-fold, P < 0.001) and myosin heavy chain (MyHC) (-70%, P < 0.001 and -69%, P < 0.001; +3.44-fold, P < 0.001) in both C2C12 cells and HSkMC. Metabolomic and transcriptomic analyses revealed that PDHB knockdown suppressed pyruvate metabolism (P < 0.001) and up-regulated ariadne RBR E3 ubiquitin protein ligase 2 (Arih2) (+7.23-fold, P < 0.001) in cellular catabolic pathways. The role of forkhead box P1 (FoxP1) (+4.18-fold, P < 0.001)-mediated Arih2 transcription was the key downstream regulator of PDHB in muscle differentiation. PDHB overexpression improved d-gal-induced muscle atrophy in mice, which was characterized by significant increases in grip strength, muscle mass and mean muscle cross-sectional area (1.19-fold to 1.5-fold, P < 0.01, P < 0.05 and P < 0.001).

CONCLUSIONS

The comprehensive results show that PDHB plays a sarcoprotective role by suppressing the FoxP1-Arih2 axis and may serve as a therapeutic target in sarcopenia.

Noncoding RNAs: Master Regulator of Fibroblast to Myofibroblast Transition in Fibrosis

Myofibroblasts escape apoptosis and proliferate abnormally under pathological conditions, especially fibrosis; they synthesize and secrete a large amount of extracellular matrix (ECM), such as α-SMA and collagen, which leads to the distortion of organ parenchyma structure, an imbalance in collagen deposition and degradation, and the replacement of parenchymal cells by fibrous connective tissues. Fibroblast to myofibroblast transition (FMT) is considered to be the main source of myofibroblasts. Therefore, it is crucial to explore the influencing factors regulating the process of FMT for the prevention, treatment, and diagnosis of FMT-related diseases. In recent years, non-coding RNAs, including microRNA, long non-coding RNAs, and circular RNAs, have attracted extensive attention from scientists due to their powerful regulatory functions, and they have been found to play a vital role in regulating FMT. In this review, we summarized ncRNAs which regulate FMT during fibrosis and found that they mainly regulated signaling pathways, including TGF-β/Smad, MAPK/P38/ERK/JNK, PI3K/AKT, and WNT/β-catenin. Furthermore, the expression of downstream transcription factors can be promoted or inhibited, indicating that ncRNAs have the potential to be a new therapeutic target for FMT-related diseases.

Transcriptomic analysis of human ALS skeletal muscle reveals a disease-specific pattern of dysregulated circRNAs

Circular RNAs are abundant, covalently closed transcripts that arise in cells through back-splicing and display distinct expression patterns across cells and developmental stages. While their functions are largely unknown, their intrinsic stability has made them valuable biomarkers in many diseases. Here, we set out to examine circRNA patterns in amyotrophic lateral sclerosis (ALS). By RNA-sequencing analysis, we first identified circRNAs and linear RNAs that were differentially abundant in skeletal muscle biopsies from ALS compared to normal individuals. By RT-qPCR analysis, we confirmed that 8 circRNAs were significantly elevated and 10 were significantly reduced in ALS, while the linear mRNA counterparts, arising from shared precursor RNAs, generally did not change. Several of these circRNAs were also differentially abundant in motor neurons derived from human induced pluripotent stem cells (iPSCs) bearing ALS mutations, and across different disease stages in skeletal muscle from a mouse model of ALS (SOD1G93A). Interestingly, a subset of the circRNAs significantly elevated in ALS muscle biopsies were significantly reduced in the spinal cord samples from ALS patients and ALS (SOD1G93A) mice. In sum, we have identified differentially abundant circRNAs in ALS-relevant tissues (muscle and spinal cord) that could inform about neuromuscular molecular programs in ALS and guide the development of therapies.

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.

Role of noncoding RNAs in pancreatic ductal adenocarcinoma associated cachexia

Cachexia is an acute syndrome that is very commonly observed in patients with cancer. Cachexia is the number one cause of death in patients with metastatic disease and is also the major factor for physical toxicity and financial burden. More importantly, the majority of patients with advanced-stage pancreatic ductal adenocarcinoma (PDAC) cancer undergo cachexia. Pancreatic cancer causes deaths of ∼50,000 Americans and about 400,000 people worldwide every year. The high mortality rates in metastatic PDAC are due to systemic pathologies and cachexia, which quickens death in these patients. About 90% of all patients with PDAC undergo wasting of muscle causing mobility loss and leading to a number of additional pathological conditions. PDAC-associated cancer cachexia emanates from complex signaling cues involving both mechanical and biological signals. Tumor invasion is associated with the loss of pancreatic function-induced digestive disorders and malabsorption, which causes subsequent weight loss and eventually promotes cachexia. Besides, systemic inflammation of patients with PDAC could release chemical cues (e.g., cytokine-mediated Atrogin-1/MAFbx expression) that participate in muscle wasting. Our understanding of genes, proteins, and cytokines involved in promoting cancer cachexia has evolved considerably. However, the role of epigenetic factors, particularly the role of noncoding RNAs (ncRNAs) in regulating PDAC-associated cachexia is less studied. In this review article, the most updated knowledge on the various ncRNAs including microRNAs (miRs), long noncoding RNA (lncRNAs), piwi interacting RNAs (PiwiRNAs), small nucleolar RNA (snoRNAs), and circular RNAs (circRNA) and their roles in cancer cachexia are described.

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.

CircRNAs Related to Breast Muscle Development and Their Interaction Regulatory Network in Gushi Chicken

Circular RNAs (circRNAs) play a significant regulatory role during skeletal muscle development. To identify circRNAs during postnatal skeletal muscle development in chickens, we constructed 12 cDNA libraries from breast muscle tissues of Chinese Gushi chickens at 6, 14, 22, and 30 weeks and performed RNA sequencing. In total, 2112 circRNAs were identified, and among them 79.92% were derived from exons. CircRNAs are distributed on all chromosomes of chickens, especially chromosomes 1-9 and Z. Bioinformatics analysis showed that each circRNA had an average of 38 miRNA binding sites, 61.32% of which have internal ribosomal entry site (IRES) elements. Furthermore, in total 543 differentially expressed circRNAs (DE-circRNAs) were identified. Functional enrichment analysis revealed that DE-circRNAs source genes are engaged in biological processes and muscle development-related pathways; for example, cell differentiation, sarcomere, and myofibril formation, mTOR signaling pathway, and TGF-β signaling pathway, etc. We also established a competitive endogenous RNA (ceRNA) regulatory network associated with skeletal muscle development. The results in this report indicate that circRNAs can mediate the development of chicken skeletal muscle by means of a complex ceRNA network among circRNAs, miRNAs, genes, and pathways. The findings of this study might help increase the number of known circRNAs in skeletal muscle tissue and offer a worthwhile resource to further investigate the function of circRNAs in chicken skeletal muscle development.

Identification and Functional Prediction of Circular RNAs Related to Growth Traits and Skeletal Muscle Development in Duroc pigs

Porcine skeletal muscle is a highly heterogeneous tissue type, and the Longissimus Dorsi muscle (LDM), as the most economical and physiologically metabolized skeletal muscle in pigs, has always been the focus of research and improvement in pig molecular breeding. Circular RNA, as an important new member of regulatory non-coding RNA after microRNA, has become a frontier hot spot in life science research. This study aims to explore candidate circRNAs related to growth, meat quality, and skeletal muscle development among Duroc pigs with different average daily gain (ADG). Eight pigs were selected and divided into two groups: H group (high-ADG) and L group (low-ADG), followed by RNA-Seq analysis to identify circRNAs. The results showed that backfat at 6-7 rib (BF) and Intramuscular fat (IMF) content in the H group was significantly lower than L group, but ribeye area (REA) in the H group was higher than in the L group. In RNA-seq, 296 Differentially expressed (DE) circRNAs (157 upregulated and 139 downregulated) were identified and exons flanking long introns are easier to circularize to produce circRNAs. Most of the DE circRNAs were enriched in Quantitative trait locus (QTL) regions related to meat quality and growth traits. In addition, a gene can produce one or more circRNA transcripts. It was also found that the source genes of DE circRNAs were enriched in MAPK, FoXO, mTOR, PI3K-Akt, and Wnt signaling pathways. The results showed that different ADG, carcass, and meat quality traits among half-sibling Duroc pigs with the same diet may be due to the DE circRNAs related to skeletal muscle growth and development.

CircARID1A binds to IGF2BP3 in gastric cancer and promotes cancer proliferation by forming a circARID1A-IGF2BP3-SLC7A5 RNA–protein ternary complex

Background

Gastric cancer (GC) is one of the most common malignant tumors in China. Circular RNAs (circRNAs) are novel non-coding RNAs with important regulatory roles in cancer progression. IGF2BP3 has been found to play oncogenic roles in various cancers including GC, while the exact mechanism of IGF2BP3 is largely unknown.

Methods

The expression of IGF2BP3 in GC was evaluated by Western Blot and bioinformatics analysis. CircRNA expression profiles were screened via IGF2BP3 RIP-seq in GC. Sanger sequencing, RNase R digestion, nucleo-plasmic separation and RNA-FISH assays were used to detect the existence and expression of circARID1A. RNA ISH assay was employed to test the expression of circARID1A in paraffin-embedded GC tissues. Moreover, the function of circARID1A on cellular proliferation was assessed by CCK-8, plate colony formation, EdU assays and GC xenograft mouse model in vivo. Furthermore, the location or binding of circARID1A, IGF2BP3 protein and SLC7A5 in GC was evaluated by RNA-FISH/IF or RNA pull-down assays.

Results

We identified a novel circRNA, circARID1A, that can bind to IGF2BP3 protein. CircARID1A was significantly upregulated in GC tissues compared with noncancerous tissues and positively correlated with tumor length, tumor volume, and TNM stage. CircARID1A knockdown inhibited the proliferation of GC cells in vitro and in vivo and circARID1A played an important role in the oncogenic function of IGF2BP3. Mechanistically, circARID1A served as a scaffold to facilitate the interaction between IGF2BP3 and SLC7A5 mRNA, finally increasing SLC7A5 mRNA stability. Additionally, circARID1A was able to directly bind SLC7A5 mRNA through complementary base-pairing and then formed the circARID1A-IGF2BP3-SLC7A5 RNA–protein ternary complex and promoted the proliferation of GC via regulating AKT/mTOR pathway.

Conclusions

Altogether, our data suggest that circARID1A is involved in the function of IGF2BP3 and GC proliferation, and the circARID1A-IGF2BP3-SLC7A5 axis has the potential to serve as a novel therapeutic target for GC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02466-3.

Integration of RNA-seq and ATAC-seq identifies muscle-regulated hub genes in cattle

As the main product of livestock, muscle itself plays an irreplaceable role in maintaining animal body movement and regulating metabolism. Therefore, it is of great significance to explore its growth, development and regeneration to improve the meat yield and quality of livestock. In this study, we attempted to use RNA-seq and ATAC-seq techniques to identify differentially expressed genes (DEGs) specifically expressed in bovine skeletal muscle as potential candidates for studying the regulatory mechanisms of muscle development. Microarray data from 8 tissue samples were selected from the GEO database for analysis. First, we obtained gene modules related to each tissue through WGCNA analysis. Through Gene Ontology (GO) functional annotation, the module of lightyellow (ME_lightyellow) was closely related to muscle development, and 213 hub genes were screened as follow-up research targets. Further, the difference analysis showed that, except for PREB, all other candidate hub genes were up-regulated (muscle group vs. other-group). ATAC-seq analysis showed that muscle-specific accessible chromatin regions were mainly located in promoter of genes related to muscle structure development (GO:0061061), muscle cell development (GO:0055001) and muscle system process (GO:0003012), which were involved in cAMP, CGMP-PKG, MAPK, and other signaling pathways. Next, we integrated the results of RNA-seq and ATAC-seq analysis, and 54 of the 212 candidate hub genes were identified as key regulatory genes in skeletal muscle development. Finally, through motif analysis, 22 of the 54 key genes were found to be potential target genes of transcription factor MEF2C. Including CAPN3, ACTN2, MB, MYOM3, SRL, CKM, ALPK3, MAP3K20, UBE2G1, NEURL2, CAND2, DOT1L, HRC, MAMSTR, FSD2, LRRC2, LSMEM1, SLC29A2, FHL3, KLHL41, ATXN7L2, and PDRG1. This provides a potential reference for studying the molecular mechanism of skeletal muscle development in mammals.