Analysis of dynamic and widespread lncRNA and miRNA expression in fetal sheep skeletal muscle

Understanding lncRNAs: key regulators of myogenesis and lipogenesis in farm animals

Long non-coding RNAs (lncRNAs) are RNA molecules exceeding 200 nucleotides in length. Recent studies have demonstrated their involvement in regulating gene expression and various biological processes. Among these, myogenesis and lipogenesis are particularly important because of their direct effects on muscle development and fat deposition in farm animals. These processes are crucial for determining meat quality, growth rates, and overall economic value in animal husbandry. Although the specific mechanisms through which lncRNAs influence these pathways are still under investigation, further research into their roles in muscle and fat development is crucial for optimizing farm animal breeding strategies. Here, we review the characteristics of lncRNAs, including their biogenesis, localization, and structures, with a particular focus on their association with myogenesis and adipogenesis. This review seeks to establish a theoretical foundation for enhancing farm animal production. In particular, focusing on lncRNAs may reveal how these molecules can enhance the economic traits of farm animals, thereby contributing to the optimization of farm animal breeding processes.

Lnc-MEG8 regulates yak myoblast differentiation via the miR-22-3p/RTL1 axis

BACKGROUND

The yak (Bos grunniens) is essential to the livelihoods of Tibetan people on the Qinghai-Tibet Plateau; however, its growth and productivity are constrained by the region's harsh climate and high altitude. Yak skeletal muscle myoblasts, which have evolved to thrive under these challenging conditions, offer a valuable model for investigating muscle development. In this study, we performed transcriptome profiling of yak longissimus dorsi muscle at different growth stages, identifying a key long non-coding RNA, LncRNA-XR_314844 (Lnc-MEG8), with a potential role in muscle development.

RESULTS

We developed a novel technique to isolate high-quality yak myoblasts, enabling detailed analysis of Lnc-MEG8. Our results indicated that Lnc-MEG8's subcellular localization varies during muscle cell growth: it is found in both the nucleus and cytoplasm during proliferation but shifts mainly to the cytoplasm during differentiation. Functional experiments showed that Lnc-MEG8 promotes cell proliferation and inhibits differentiation, while its silencing had the opposite effect. Further analysis revealed that both Lnc-MEG8 and the gene RTL1 share miR-22-3p as a common target. Dual-luciferase assays confirmed miR-22-3p directly targets both Lnc-MEG8 and RTL1 mRNA. Co-transfection of Lnc-MEG8 and a miR-22-3p mimic restored RTL1 expression, highlighting Lnc-MEG8's regulatory role. Lnc-MEG8 also counteracts miR-22-3p's suppression of key muscle genes such as MyF5 and MyoG, facilitating myotube formation.

CONCLUSION

These findings demonstrate that the Lnc-MEG8-miR-22-3p-RTL1 axis plays a crucial role in yak muscle development, providing insights that could advance muscle tissue engineering and enhance yak meat quality.

Integrated analysis of mRNAs and lncRNAs reveals candidate marker genes and potential hub lncRNAs associated with growth regulation of the Pacific Oyster, Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is an economically important shellfish around the world. Great efforts have been made to improve its growth rate through genetic breeding. However, the candidate marker genes, pathways, and potential lncRNAs involved in oyster growth regulation remain largely unknown. To identify genes, lncRNAs, and pathways involved in growth regulation, C. gigas spat was cultured at a low temperature (15 ℃) to yield a growth-inhibited model, which was used to conduct comparative transcriptome analysis with spat cultured at normal temperature (25 ℃).

RESULTS

In total, 8627 differentially expressed genes (DEGs) and 1072 differentially expressed lncRNAs (DELs) were identified between the normal-growth oysters (cultured at 25 ℃, hereinafter referred to as NG) and slow-growth oysters (cultured at 15 ℃, hereinafter referred to as SG). Functional enrichment analysis showed that these DEGs were mostly enriched in the AMPK signaling pathway, MAPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process. LncRNAs analysis identified 265 cis-acting pairs and 618 trans-acting pairs that might participate in oyster growth regulation. The expression levels of LNC_001270, LNC_003322, LNC_011563, LNC_006260, and LNC_012905 were inducible to the culture temperature and food abundance. These lncRNAs were located at the antisense, upstream, or downstream of the SREBP1/p62, CDC42, CaM, FAS, and PIK3CA genes, respectively. Furthermore, the expression of the trans-acting lncRNAs, including XR_9000022.2, LNC_008019, LNC_015817, LNC_000838, LNC_00839, LNC_011859, LNC_007294, LNC_006429, XR_002198885.1, and XR_902224.2 was also significantly associated with the expression of genes enriched in AMPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process.

CONCLUSIONS

In this study, we identified the critical growth-related genes and lncRNAs that could be utilized as candidate markers to illustrate the molecular mechanisms underlying the growth regulation of Pacific oysters.

Global Long Noncoding RNA Expression Profiling of MSTN and FGF5 Double-Knockout Sheep Reveals the Key Gatekeepers of Skeletal Muscle Development

Improving livestock and poultry growth rates and increasing meat production are urgently needed worldwide. Previously, we produced a myostatin (MSTN) and fibroblast growth factor 5 (FGF5) double-knockout (MF^-/-) sheep by CRISPR Cas9 system to improve meat production, and also wool production. Both MF^-/- sheep and the F1 generation (MF^+/-) sheep showed an obvious "double-muscle" phenotype. In this study, we identified the expression profiles of long noncoding RNAs (lncRNAs) in wild-type and MF^+/- sheep, then screened out the key candidate lncRNAs that can regulate myogenic differentiation and skeletal muscle development. These key candidate lncRNAs can serve as critical gatekeepers for muscle contraction, calcium ion transport and skeletal muscle cell differentiation, apoptosis, autophagy, and skeletal muscle inflammation, further revealing that lncRNAs play crucial roles in regulating muscle phenotype in MF^+/- sheep. In conclusion, our newly identified lncRNAs may emerge as novel molecules for muscle development or muscle disease and provide a new reference for MSTN-mediated regulation of skeletal muscle development.

Identification and Characterization of lncRNAs Expression Profile Related to Goat Skeletal Muscle at Different Development Stages

LncRNAs are essential for regulating skeletal muscle. However, the expression profile and function of lncRNAs in goat muscle remains unclear. Here, an average of ~14.58 Gb high-quality reads were obtained from longissimus dorsi tissues of 1-month-old (n = 3) and 9-month-old (n = 3) Wu'an black goats using RNA sequencing. Of a total of 3441 lncRNAs, 1281 were lincRNAs, 805 were antisense lncRNAs, and 1355 were sense_overlapping lncRNAs. These lncRNAs shared some properties with goats, such as fewer exons, shorter transcript, and open reading frames (ORFs) length. Among them, 36 differentially expressed lncRNAs (DE lncRNA) were identified, and then 10 random lncRNAs were validated by RT-qPCR. Furthermore, 30 DE lncRNAs were neighboring 71 mRNAs and several genes were functionally enriched in muscle development-related pathways, such as APC, IFRD1, NKX2-5, and others. Additionally, 36 DE lncRNAs and 2684 mRNAs were included in co-expression interactions. A lncRNA-miRNA-mRNA network containing 4 lncRNAs, 3 miRNAs, and 8 mRNAs was finally constructed, of which XR_001296113.2 might regulate PDLIM7 expression by interaction with chi-miR-1296 to affect skeletal muscle development. This study revealed the expression profile of goat lncRNAs for further investigative studies and provides a fuller understanding of skeletal muscle development.

Integrative analysis of transcriptomics and proteomics of longissimus thoracis of the Hu sheep compared with the Dorper sheep

Meat quality is becoming more important for sheep breeding programs. Meat quality is a complex trait affected by genetic and environmental factors. In the present study, an integrative analysis of the longissimus thoracis tissue transcriptome and proteome was conducted to identify genes, proteins, and pathways related to meat quality in sheep. The sheep breeds Hu and Dorper were considered. These breeds were compared for the differences in muscle fiber structure, chemical composition, and amino acid composition. In the Hu sheep vs. Dorper sheep comparison, 22 DEGs/DEPs showed the same mRNA and protein expression trends. These genes are associated with lipid transport, lipid metabolism, and muscular system development. Moreover, some pathways such as "lipid transport", "lipoprotein metabolic process", "Alanine, aspartate and glutamate metabolism", and "Arginine biosynthesis" were significantly enriched in this study. The reliability of the RNA-Seq results was verified by qRT-PCR. These findings provide new insights into the molecular mechanisms of meat quality in sheep.

Deep Small RNA Sequencing Reveals Important miRNAs Related to Muscle Development and Intramuscular Fat Deposition in Longissimus dorsi Muscle From Different Goat Breeds

MicroRNAs (miRNAs) are a class of small non-coding RNAs that have been shown to play important post-transcriptional regulatory roles in the growth and development of skeletal muscle tissues. However, limited research into the effect of miRNAs on muscle development in goats has been reported. In this study, Liaoning cashmere (LC) goats and Ziwuling black (ZB) goats with significant phenotype difference in meat production performance were selected and the difference in Longissimus dorsi muscle tissue expression profile of miRNAs between the two goat breeds was then compared using small RNA sequencing. A total of 1,623 miRNAs were identified in Longissimus dorsi muscle tissues of the two goat breeds, including 410 known caprine miRNAs, 928 known species-conserved miRNAs and 285 novel miRNAs. Of these, 1,142 were co-expressed in both breeds, while 230 and 251 miRNAs were only expressed in LC and ZB goats, respectively. Compared with ZB goats, 24 up-regulated miRNAs and 135 miRNAs down-regulated were screened in LC goats. A miRNA-mRNA interaction network showed that the differentially expressed miRNAs would target important functional genes associated with muscle development and intramuscular fat deposition. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the target genes of differentially expressed miRNAs were significantly enriched in Ras, Rap 1, FoxO, and Hippo signaling pathways. This study suggested that these differentially expressed miRNAs may be responsible for the phenotype differences in meat production performance between the two goat breeds, thereby providing an improved understanding of the roles of miRNAs in muscle tissue of goats.

Expression Profile Analysis to Identify Circular RNA Expression Signatures in Muscle Development of Wu'an Goat Longissimus Dorsi Tissues

The growth and development of skeletal muscle is a physiological process regulated by a variety of genes and signaling pathways. As a posttranscriptional regulatory factor, circRNA plays a certain regulatory role in the development of animal skeletal muscle in the form of a miRNA sponge. However, the role of circRNAs in muscle development and growth in goats is still unclear. In our study, apparent differences in muscle fibers in Wu'an goats of different ages was firstly detected by hematoxylin-eosin (HE) staining, the circRNA expression profiles of longissimus dorsi muscles from 1-month-old (mon1) and 9-month-old (mon9) goats were screened by RNA-seq and verified by RT-qPCR. The host genes of differentially expressed (DE) circRNAs were predicted, and Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes analyses (KEGG) of host genes with DE circRNAs were performed to explore the functions of circRNAs. The circRNA-miRNA-mRNA networks were then constructed using Cytoscape software. Ten significantly differentially expressed circRNAs were also verified in the mon1 and mon9 groups by RT-qPCR. Luciferase Reporter Assay was used to verify the binding site between circRNA and its targeted miRNA. The results showed that a total of 686 DE circRNAs were identified between the mon9 and mon1 groups, of which 357 were upregulated and 329 were downregulated. Subsequently, the 467 host genes of DE circRNAs were predicted using Find_circ and CIRI software. The circRNA-miRNA-mRNA network contained 201 circRNAs, 85 miRNAs, and 581 mRNAs; the host mRNAs were associated with "muscle fiber development" and "AMPK signaling pathway" and were enriched in the FoxO signaling pathway. Competing endogenous RNA (ceRNA) network analysis showed that novel_circ_0005314, novel_circ_0005319, novel_circ_0009256, novel_circ_0009845, novel_circ_0005934 and novel_circ_0000134 may play important roles in skeletal muscle growth and development between the mon9 and mon1 groups. Luciferase Reporter Assay confirmed the combination between novel_circ_0005319 and chi-miR-199a-5p, novel_circ_0005934 and chi-miR-450-3p and novel_circ_0000134 and chi-miR-655. Our results provide specific information related to goat muscle development and a reference for the goat circRNA profile.

RNA-Seq Reveals miRNA and mRNA Co-regulate Muscle Differentiation in Fetal Leizhou Goats

Muscle differentiation is an essential link in animal growth and development, and microRNA and mRNA are indispensable in skeletal muscle differentiation. To improve the meat quality and production of the Leizhou goat, it is vital to understand the molecular mechanism by which its skeletal muscle differentiates. By RNA sequencing (RNA-SEQ), we established miRNA-mRNA profiles of Leizhou goats at three stages: fetal day 70, 90, and 120. There were 991 differently expressed mRNAs and 39 differentially expressed miRNAs found, with the differentially expressed mRNAs mainly enriched in calcium ion binding, ECM-receptor interaction, and Focal adhesion. CKM and MYH3, two muscle differentiation markers, were significantly differentially expressed during this period. In addition, we found that chi-miR-129-5p, chi-miR-433, and chi-miR-24-3p co-regulate muscle differentiation with their target genes. Finally, we can confirm that muscle differentiation occurred in Leizhou goat between 90 and 120 days of the fetus. This study is helpful to better explore the molecular mechanism of goat muscle differentiation.

MiR-22-3p Inhibits Proliferation and Promotes Differentiation of Skeletal Muscle Cells by Targeting IGFBP3 in Hu Sheep

The growth and development of skeletal muscle require a series of regulatory factors. MiRNA is a non-coding RNA with a length of about 22 nt, which can inhibit the expression of mRNA and plays an important role in the growth and development of muscle cells. The role of miR-22-3p in C2C12 cells and porcine skeletal muscle has been reported, but it has not been verified in Hu sheep skeletal muscle. Through qPCR, CCK-8, EdU and cell cycle studies, we found that overexpression of miR-22-3p inhibited proliferation of skeletal muscle cells (p < 0.01). The results of qPCR and immunofluorescence showed that overexpression of miR-22-3p promoted differentiation of skeletal muscle cells (p < 0.01), while the results of inhibiting the expression of miR-22-3p were the opposite. These results suggested that miR-22-3p functions in growth and development of sheep skeletal muscle cells. Bioinformatic analysis with mirDIP, miRTargets, and RNAhybrid software suggested IGFBP3 was the target of miR-22-3p, which was confirmed by dual-luciferase reporter system assay. IGFBP3 is highly expressed in sheep skeletal muscle cells. Overexpression of IGFBP3 was found to promote proliferation of skeletal muscle cells indicated by qPCR, CCK-8, EdU, and cell cycle studies (p < 0.01). The results of qPCR and immunofluorescence experiments proved that overexpression of IGFBP3 inhibited differentiation of skeletal muscle cells (p < 0.01), while the results of interfering IGFBP3 with siRNA were the opposite. These results indicate that miR-22-3p is involved in proliferation and differentiation of skeletal muscle cells by targeting IGFBP3.

Comparison of MicroRNA Transcriptomes Reveals the Association between MiR-148a-3p Expression and Rumen Development in Goats

Simple Summary

In ruminants, the rumen epithelium plays an important role in nutrient absorption, metabolism and transport. MicroRNAs (miRNAs) have been reported to regulate the proliferation of diverse epithelial cells. In this study, we profiled the miRNA transcriptomes of goat rumens at four development stages and screened for candidate miRNAs related to rumen development. MiR-148a-3p was found to be highly expressed in the rumen tissues and induced the proliferation of GES-1 cells by targeting QKI. Our findings provide some insights into the functional roles of miRNAs in rumen growth and functional development in ruminants.

Abstract

The rumen is an important digestive organ of ruminants. From the fetal to adult stage, the morphology, structure and function of the rumen change significantly. However, the knowledge of the intrinsic genetic regulation of these changes is still limited. We previously reported a genome-wide expression profile of miRNAs in pre-natal goat rumens. In this study, we combined and analyzed the transcriptomes of rumen miRNAs during pre-natal (E60 and E135) and post-natal (D30 and D150) stages. A total of 66 differentially expressed miRNAs (DEMs) were identified in the rumen tissues from D30 and D150 goats. Of these, 17 DEMs were consistently highly expressed in the rumens at the pre-weaning stages (E60, E135 and D30), while down-regulated at D150. Noteworthy, annotation analysis revealed that the target genes regulated by the DEMs were mainly enriched in MAPK signaling pathway, Jak-STAT signaling pathway and Ras signaling pathway. Interestingly, the expression of miR-148a-3p was significantly high in the embryonic stage and down-regulated at D150. The potential binding sites of miR-148a-3p in the 3′-UTR of QKI were predicted by the TargetScan and verified by the dual luciferase report assay. The co-localization of miR-148a-3p and QKI through in situ hybridization was observed in the rumen tissues but not in the intestinal tracts. Moreover, the expression of miR-148a-3p in the epithelium was significantly higher than that in the other layers of the rumen, suggesting that miR-148a-3p is involved in the development of the rumen epithelial cells by targeting QKI. Subsequently, miR-148a-3p inhibitor was found to induce the proliferation of GES-1 cells. Taken together, our study identified DEMs involved in the development of the rumen and provides insights into the regulation mechanism of rumen development in goats.