Transcriptome RNA Sequencing Reveals That Circular RNAs Are Abundantly Expressed in Embryonic Breast Muscle of Duck

Signaling pathways and regulatory networks in quail skeletal muscle development: insights from whole transcriptome sequencing

Quail, as an advantageous avian model organism due to its compact size and short reproductive cycle, holds substantial potential for enhancing our understanding of skeletal muscle development. The quantity of skeletal muscle represents a vital economic trait in poultry production. Unraveling the molecular mechanisms governing quail skeletal muscle development is of paramount importance for optimizing meat and egg yield through selective breeding programs. However, a comprehensive characterization of the regulatory dynamics and molecular control underpinning quail skeletal muscle development remains elusive. In this study, through the application of HE staining on quail leg muscle sections, coupled with preceding fluorescence quantification PCR of markers indicative of skeletal muscle differentiation, we have delineated embryonic day 9 (E9) and embryonic day 14 (E14) as the start and ending points, respectively, of quail skeletal muscle differentiation. Then, we employed whole transcriptome sequencing to investigate the temporal expression profiles of leg muscles in quail embryos at the initiation of differentiation (E9) and upon completion of differentiation (E14). Our analysis revealed the expression patterns of 12,012 genes, 625 lncRNAs, 14,457 circRNAs, and 969 miRNAs in quail skeletal muscle samples. Differential expression analysis between the E14 and E9 groups uncovered 3,479 differentially expressed mRNAs, 124 lncRNAs, 292 circRNAs, and 154 miRNAs. Furthermore, enrichment analysis highlighted the heightened activity of signaling pathways related to skeletal muscle metabolism and intermuscular fat formation, such as the ECM-receptor interaction, focal adhesion, and PPAR signaling pathway during E14 skeletal muscle development. Conversely, the E9 stage exhibited a prevalence of pathways associated with myoblast proliferation, exemplified by cell cycle processes. Additionally, we constructed regulatory networks encompassing lncRNA‒mRNA, miRNA‒mRNA, lncRNA‒miRNA-mRNA, and circRNA-miRNA‒mRNA interactions, thus shedding light on their putative roles within quail skeletal muscle. Collectively, our findings illuminate the gene and non-coding RNA expression characteristics during quail skeletal muscle development, serving as a foundation for future investigations into the regulatory mechanisms governing non-coding RNA and quail skeletal muscle development in poultry production.

A novel protein encoded by circKANSL1L regulates skeletal myogenesis via the Akt-FoxO3 signaling axis

Skeletal muscle is one the largest organs of the body and is involved in animal production and human health. Circular RNAs (circRNAs) have been implicated in skeletal myogenesis through largely unknown mechanisms. Herein, we report the phenotypic and metabolomic analysis of porcine longissimus dorsi muscles in Lantang and Landrace piglets, revealing a high-content of slow-oxidative fibers responsible for high-quality meat product in Lantang piglets. Using single-cell transcriptomics, we identified four myogenesis-related cell types, and the Akt-FoxO3 signaling axis was the most significantly enriched pathway in each subpopulation in the different pig breeds, as well as in fast-twitch glycolytic fibers. Using the multi-dimensional bioinformatic tools of circRNAome-seq and Ribo-seq, we identified a novel circRNA, circKANSL1L, with a protein-coding ability in porcine muscles, whose expression level correlated with myoblast proliferation and differentiation in vitro, as well as the transformation between distinct mature myofibers in vivo. The protein product of circKANSL1L could interact with Akt to decrease the phosphorylation level of FoxO3, which subsequently promoted FoxO3 transcriptional activity to regulate skeletal myogenesis. Our results established the existence of a protein encoded by circKANSL1L and demonstrated its potential functions in myogenesis.

Whole-Transcriptome RNA Sequencing Uncovers the Global Expression Changes and RNA Regulatory Networks in Duck Embryonic Myogenesis

Duck meat is pivotal in providing high-quality protein for human nutrition, underscoring the importance of studying duck myogenesis. The regulatory mechanisms governing duck myogenesis involve both coding and non-coding RNAs, yet their specific expression patterns and molecular mechanisms remain elusive. To address this knowledge gap, we performed expression profiling analyses of mRNAs, lncRNAs, circRNAs, and miRNAs involved in duck myogenesis using whole-transcriptome RNA-seq. Our analysis identified 1733 differentially expressed (DE)-mRNAs, 1116 DE-lncRNAs, 54 DE-circRNAs, and 174 DE-miRNAs when comparing myoblasts and myotubes. A GO analysis highlighted the enrichment of DE molecules in the extracellular region, protein binding, and exocyst. A KEGG analysis pinpointed pathways related to ferroptosis, PPAR signaling, nitrogen metabolism, cell cycle, cardiac muscle contraction, glycerolipid metabolism, and actin cytoskeleton. A total of 51 trans-acting lncRNAs, including ENSAPLT00020002101 and ENSAPLT00020012069, were predicted to participate in regulating myoblast proliferation and differentiation. Based on the ceRNAs, we constructed lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA ceRNA networks involving five miRNAs (miR-129-5p, miR-133a-5p, miR-22-3p, miR-27b-3p, and let-7b-5p) that are relevant to myogenesis. Furthermore, the GO and KEGG analyses of the DE-mRNAs within the ceRNA network underscored the significant enrichment of the glycerolipid metabolism pathway. We identified five different DE-mRNAs, specifically ENSAPLG00020001677, ENSAPLG00020002183, ENSAPLG00020005019, ENSAPLG00020010497, and ENSAPLG00020017682, as potential target genes that are crucial for myogenesis in the context of glycerolipid metabolism. These five mRNAs are integral to ceRNA networks, with miR-107_R-2 and miR-1260 emerging as key regulators. In summary, this study provides a valuable resource elucidating the intricate interplay of mRNA-lncRNA-circRNA-miRNA in duck myogenesis, shedding light on the molecular mechanisms that govern this critical biological process.

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.

Cloning of the RNA m6A Methyltransferase 3 and Its Impact on the Proliferation and Differentiation of Quail Myoblasts

Methyltransferase 3 (METTL3), which has been demonstrated to play a crucial role in a variety of biological processes, is the key enzyme for catalyzing m⁶A modification in RNA. However, the complete protein sequence of METTL3 in quail has not been annotated, and its function in skeletal muscle of quails remains unknown. In the current study, the full-length coding sequence of the quail METTL3 was obtained through the 3' rapid amplification of cDNA ends (3' RACE) and its homology with that of other species was predicted based on a generated phylogenetic tree. A Cell Counting Kit-8 assay and flow cytometry in a quail myoblast cell line (QM7) demonstrated that METTL3 promotes myoblast proliferation. The overexpression of METTL3 in QM7 cells significantly increased the expression levels of the myoblast differentiation markers myogenin (MYOG), myogenic differentiation 1 (MYOD1), and myocyte enhancer factor 2C (MEF2C), further demonstrating that METTL3 promotes myoblast differentiation. Additionally, transcriptome sequencing following METTL3 overexpression revealed that METTL3 controls the expression of various genes involved in RNA splicing and the regulation of gene expression, as well as pathways such as the MAPK signaling pathway. Taken together, our findings demonstrated that METTL3 plays a vital function in quail myoblast proliferation and differentiation and that the METTL3-mediated RNA m⁶A modification represents an important epigenetic regulatory mechanism in poultry skeletal muscle development.