Goat Skeletal Muscle Research Articles

Single-nucleus transcriptional profiling reveals TCF7L2 as a key regulator in adipogenesis in goat skeletal muscle development

Intramuscular adipogenesis plays an important role in muscle development, which determines the quality of goat meat. However, its underlying cellular and molecular mechanisms remain poorly understood. In this study, we provided detailed cellular atlases of goat longissimus dorsi during muscle development at single-nucleus resolution, and identified the subpopulations of fibroblasts/fibro-adipogenic progenitors (FAPs) and muscle satellite cell (MuSC), as well as the differentiation trajectory of FAPs subpopulations. Cellular ligand-receptor interaction analysis revealed enriched BMP and IGF pathways implicated in within-tissue crosstalk centered around FAPs. Through single-nucleus gene regulatory network analysis and in vitro interference verification, we found that TCF7L2 was a critical transcriptional factor (TF) in early adipogenesis in skeletal muscle. Overall, our work reveals the cellular intricacies and diversity of goat longissimus dorsi during muscle development, implementing insights into the critical roles of BMP, IGF pathways and TCF7L2 TF in intramuscular adipogenesis.

ADAR1 Promotes Myogenic Proliferation and Differentiation of Goat Skeletal Muscle Satellite Cells.

As one of the most important economic traits for domestic animal husbandry, skeletal muscle is regulated by an intricate molecular network. Adenosine deaminase acting on RNA (ADAR1) involves various physiological processes and diseases, such as innate immunity and the development of lung adenocarcinoma, breast cancer, gastric cancer, etc. However, its role in skeletal muscle growth requires further clarification. Here, we explored the functions of ADAR1 in the myogenic process of goat skeletal muscle satellite cells (MuSCs). The ADAR1 transcripts were noticeably enriched in goat visceral tissues compared to skeletal muscle. Additionally, its levels in slow oxidative muscles like the psoas major and minor muscles were higher than in the fast oxidative glycolytic and fast glycolytic muscles. Among the two common isoforms from ADAR1, p110 is more abundant than p150. Moreover, overexpressing ADAR1 enhanced the proliferation and myogenic differentiation of MuSCs. The mRNA-seq performed on MuSCs' knockdown of ADAR1 obtained 146 differentially expressed genes (DEGs), 87 upregulated and 59 downregulated. These DEGs were concentrated in muscle development and process pathways, such as the MAPK and cAMP signaling pathways. Furthermore, many DEGs as the key nodes defined by protein-protein interaction networks (PPI), including STAT3, MYH3/8, TGFβ2, and ACTN4, were closely related to the myogenic process. Finally, RNA immunoprecipitation combined with qPCR (RIP-qPCR) showed that ADAR1 binds to PAX7 and MyoD mRNA. This study indicates that ADAR1 promotes the myogenic development of goat MuSCs, which provides a useful scientific reference for further exploring the ADAR1-related regulatory networks underlying mammal skeletal muscle growth.

Regulatory role of lncMD1 in goat skeletal muscle satellite cell differentiation via miR-133a-3p and miR-361-3p targeting

Skeletal muscle satellite cells (SMSCs) are pivotal in skeletal muscle development and are influenced by numerous regulatory factors. This study focuses on the regulatory and functional mechanism roles of lncMD1, a muscle-specific long non-coding RNA, in the proliferation and differentiation of goat SMSCs. Employing in vitro cultured goat SMSCs, this study demonstrated that lncMD1, functions as a decoy for miR-133a-3p and miR-361-3p. This interaction competitively binds these microRNAs to modulate the expression of dynactin subunit 2 (DCTN2) and dynactin subunit 1 (DCTN1), thereby affects SMSCs proliferation and differentiation. These findings enhance the understanding of non-coding RNAs in goat SMSCs growth cycles and offer a theoretical foundation for exploring the molecular processes of goat skeletal muscle myogenic development.

Fat Mass- and Obesity-Associated Protein (FTO) Promotes the Proliferation of Goat Skeletal Muscle Satellite Cells by Stabilizing DAG1 mRNA in an IGF2BP1-Related m6A Manner.

Skeletal muscle development is spotlighted in mammals since it closely relates to animal health and economic benefits to the breeding industry. Researchers have successfully unveiled many regulatory factors and mechanisms involving myogenesis. However, the effect of N6-methyladenosine (m6A) modification, especially demethylase and its regulated genes, on muscle development remains to be further explored. Here, we found that the typical demethylase FTO (fat mass- and obesity-associated protein) was highly enriched in goats' longissimus dorsi (LD) muscles. In addition, the level of m6A modification on transcripts was negatively regulated by FTO during the proliferation of goat skeletal muscle satellite cells (MuSCs). Moreover, a deficiency of FTO in MuSCs significantly retarded their proliferation and promoted the expression of dystrophin-associated protein 1 (DAG1). m6A modifications of DAG1 mRNA were efficiently altered by FTO. Intriguingly, the results of DAG1 levels and its m6A enrichment from FB23-2 (FTO demethylase inhibitor)-treated cells were consistent with those of the FTO knockdown, indicating that the regulation of FTO on DAG1 depended on m6A modification. Further experiments showed that interfering FTO improved m6A modification at site DAG1-122, recognized by Insulin-like growth factor 2 mRNA-binding protein 1 (IGF2BP1) and consequently stabilized DAG1 transcripts. Our study suggests that FTO promotes the proliferation of MuSCs by regulating the expression of DAG1 through m6A modification. This will extend our knowledge of the m6A-related mechanism of skeletal muscle development in animals.

Role of circPAPD7 in regulating proliferation and differentiation of goat skeletal muscle satellite cells

The circular RNA (circRNA) plays a crucial role in various biological processes, particularly posttranscriptional regulation. However, the role of circRNA in the development of goat skeletal muscle has not been thoroughly explored. Here, we identified circPAPD7, which is a novel circular RNA that is preferentially expressed in the skeletal muscle. Functional assays demonstrated that circPAPD7 promoted proliferation and inhibited differentiation in goat skeletal muscle satellite cells (MuSCs). Mechanistically, it was discovered that circPAPD7 interacts with miR-26a-5p. Moreover, the rescue experiments indicated that the overexpression of circPAPD7 may reverse the inhibitory impact of miR-26a-5p on myoblast proliferation and the accelerated effects on differentiation. Furthermore, we provided evidence that circPAPD7 functions as a sponge for miR-26a-5p, thereby facilitating the upregulation of EZH2 expression in goat MuSCs. Together, the results revealed that circPAPD7 promote proliferation and inhibit differentiation of goat MuSCs via the miR-26a-5p/EZH2 pathway.

Evaluation of Reference Gene Stability in Goat Skeletal Muscle Satellite Cells during Proliferation and Differentiation Phases.

The process of skeletal muscle development is intricate and involves the regulation of a diverse array of genes. Accurate gene expression profiles are crucial for studying muscle development, making it essential to choose the right reference genes for real-time quantitative PCR (RT-qPCR). In the present study, eight candidate reference genes were identified from our previous transcriptome sequencing analysis of caprine skeletal muscle satellite cells (MuSCs), and two traditional reference genes (ACTB and GAPDH) were assessed. The quantitative levels of the candidate reference genes were determined through the RT-qPCR technique, while the stability of their expression was evaluated utilizing the GeNorm, NormFinder, BestKeeper, and RefFinder programs. Furthermore, the chosen reference genes were utilized for the normalization of the gene expression levels of PCNA and Myf5. It was determined that conventional reference genes, including ACTB and GAPDH, were not appropriate for normalizing target gene expression. Conversely, RPL14 and RPS15A, identified through RNA sequencing analysis, exhibited minimal variability and were identified as the optimal reference genes for normalizing gene expression during the proliferation and differentiation of goat MuSCs. Our research offers a validated panel of optimal reference genes for the detection of differentially expressed genes in goat muscle satellite cells using RT-qPCR.

Role and Regulatory Mechanism of circRNA_14820 in the Proliferation and Differentiation of Goat Skeletal Muscle Satellite Cells.

Skeletal muscle satellite cells (SMSCs), a type of myogenic stem cell, play a pivotal role in postnatal muscle regeneration and repair in animals. Circular RNAs (circRNAs) are a distinct class of non-coding RNA molecules capable of regulating muscle development by modulating gene expression, acting as microRNAs, or serving as protein decoys. In this study, we identified circ_14820, an exonic transcript derived from adenosine triphosphatase family protein 2 (ATAD2), through initial RNA-Seq analysis. Importantly, overexpression of circ_14820 markedly enhanced the proliferation of goat SMSCs while concomitantly suppressing their differentiation. Moreover, circ_14820 exhibited predominant localization in the cytoplasm of SMSCs. Subsequent small RNA and mRNA sequencing of circ_14820-overexpressing SMSCs systematically elucidated the molecular regulatory mechanisms associated with circ_14820. Our preliminary findings suggest that the circ_14820-miR-206-CCND2 regulatory axis may govern the development of goat SMSCs. These discoveries contribute to a deeper understanding of circRNA-mediated mechanisms in regulating skeletal muscle development, thereby advancing our knowledge of muscle biology.

Mechanistic of AMPK/ACC2 regulating myoblast differentiation by fatty acid oxidation of goat

Myoblast differentiation depends on fatty acid oxidation (FAO)，and its rate-limiting enzyme acetyl-CoA carboxylase 2 (ACC2) participate in the regulation skeletal muscle development. However, the precise regulatory mechanism is still unknown. Using previous RNA-sequencing data from our laboratory, we explored the effect of ACC2 on myoblast differentiation, as a candidate gene, since its expression is higher in myoblasts of lamb (first day of age) than that of the fetus (75th day of pregnancy). Our findings show that siACC2 inhibited myoblast proliferation, promoted differentiation, and boosted mitochondrial and fatty acid oxidation activities. The effect of ACC2 on goat muscle cell differentiation was modulated by Etomoxir, a CPT1A inhibitor. Notably, the AMPK/ACC2 pathway was found to regulate fatty acid oxidation and goat muscle cell differentiation. Inhibiting the AMPK/ACC2 pathway significantly reduced CPT1A expression. These findings indicate that AMPK/ACC2 regulate goat myoblast differentiation via fatty acid oxidation, contributing to understanding the mechanism of goat skeletal muscle development.

CircUSP13 facilitates the fast-to-slow myofiber shift via the MAPK/ERK signaling pathway in goat skeletal muscles.

Understanding how skeletal muscle fiber proportions are regulated is essential for understanding muscle function and improving the quality of mutton. While circular RNA (circRNA) has a critical function in myofiber type transformation, the specific mechanisms are not yet fully understood. Prior evidence indicates that circular ubiquitin-specific peptidase 13 (circUSP13) can promote myoblast differentiation by acting as a ceRNA, but its potential role in myofiber switching is still unknown. Herein, we found that circUSP13 enhanced slow myosin heavy chain (MyHC-slow) and suppressed MyHC-fast expression in goat primary myoblasts (GPMs). Meanwhile, circUSP13 evidently enhanced the remodeling of the mitochondrial network while inhibiting the autophagy of GPMs. We obtained fast-dominated myofibers, via treatment with rotenone, and further demonstrated the positive role of circUSP13 in the fast-to-slow transition. Mechanistically, activation of the mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) pathway significantly impaired the slow-to-fast shift in fully differentiated myotubes, which was restored by circUSP13 or IGF1 overexpression. In conclusion, circUSP13 promoted the fast-to-slow myofiber type transition through MAPK/ERK signaling in goat skeletal muscle. These findings provide novel insights into the role of circUSP13 in myofiber type transition and contribute to a better understanding of the genetic mechanisms underlying meat quality.

An Identification of Functional Genetic Variants in B4GALNT2 Gene and Their Association with Growth Traits in Goats.

β-1,4-N-acetylgalactosamine transferase 2 (B4GALNT2) is a vital candidate gene that affects the growth traits in sheep. However, whether it has the same function in goats remains to be investigated further. This study selected 348 Nanjiang Yellow goats, screened all exons, and conserved non-coding regions of the B4GALNT2 gene for single-nucleotide polymorphisms (SNPs). Our results revealed the presence of a synonymous mutation, rs672215506, within the exon of the B4GALNT2 gene in the Nanjiang Yellow goat population. The mutation resulted in a decrease in the mRNA stability of the B4GALNT2 gene. The results of SNP detection of the conserved non-coding region of the B4GALNT2 gene showed five potential regulatory SNPs in the Nanjiang Yellow goat population. Except for rs66095343, the ~500 bp fragments of the other four SNPs (rs649127714, rs649573228, rs652899012, and rs639183528) significantly increased the luciferase activity both in goat skeletal muscle satellite cells (MuSCs) and 293T cells. The genetic diversity indexes indicated low or intermediate levels for all six SNPs analyzed, and the genotype frequencies were in Hardy-Weinberg equilibrium. Association analysis showed that rs660965343, rs649127714, and rs649573228 significantly correlate with growth traits in the later stage of growth and development of Nanjiang Yellow goats. The haplotype combinations of H2H3 and H2H2 had higher body weight and greater body size. Moreover, H2H2 haplotype combinations significantly correlated with the litter size of the Nanjiang Yellow goats. The results of our study demonstrate the potential role of the B4GALNT2 gene as a functional genetic marker in the breeding programs of Nanjiang Yellow goats.

METTL3 Promotes the Differentiation of Goat Skeletal Muscle Satellite Cells by Regulating MEF2C mRNA Stability in a m6A-Dependent Manner.

The development of mammalian skeletal muscle is a highly complex process involving multiple molecular interactions. As a prevalent RNA modification, N6-methyladenosine (m6A) regulates the expression of target genes to affect mammalian development. Nevertheless, it remains unclear how m6A participates in the development of goat muscle. In this study, methyltransferase 3 (METTL3) was significantly enriched in goat longissimus dorsi (LD) tissue. In addition, the global m6A modification level and differentiation of skeletal muscle satellite cells (MuSCs) were regulated by METTL3. By performing mRNA-seq analysis, 8050 candidate genes exhibited significant changes in expression level after the knockdown of METTL3 in MuSCs. Additionally, methylated RNA immunoprecipitation sequencing (MeRIP-seq) illustrated that myocyte enhancer factor 2c (MEF2C) mRNA contained m6A modification. Further experiments demonstrated that METTL3 enhanced the differentiation of MuSCs by upregulating m6A levels and expression of MEF2C. Moreover, the m6A reader YTH N6-methyladenosine RNA binding protein C1 (YTHDC1) was bound and stabilized to MEF2C mRNA. The present study reveals that METTL3 enhances myogenic differentiation in MuSCs by regulating MEF2C and provides evidence of a post-transcriptional mechanism in the development of goat skeletal muscle.

MiR-188-5p regulates the proliferation and differentiation of goat skeletal muscle satellite cells by targeting calcium/calmodulin dependent protein kinase II beta.

The aim of this study was to reveal the role and regulatory mechanism of miR-188-5p in the proliferation and differentiation of goat muscle satellite cells. Goat skeletal muscle satellite cells isolated in the pre-laboratory were used as the test material. First, the expression of miR-188-5p in goat muscle tissues at different developmental stages was detected by quantitative reverse transcription polymerase chain reaction (qRT-PCR). In addition, miR-188-5p was transfected into goat skeletal muscle satellite cells by constructing mimics and inhibitors of miR-188-5p, respectively. The changes of differentiation marker gene expression were detected by qPCR method. It was highly expressed in adult goat latissimus dorsi and leg muscles, goat fetal skeletal muscle, and at the differentiation stage of muscle satellite cells. Overexpression and interference of miR-188-5p showed that miR-188-5p inhibited the proliferation and promoted the differentiation of goat muscle satellite cells. Target gene prediction and dual luciferase assays showed that miR-188-5p could target the 3'untranslated region of the calcium/calmodulin dependent protein kinase II beta (CAMK2B) gene and inhibit luciferase activity. Further functional studies revealed that CAMK2B promoted the proliferation and inhibited the differentiation of goat muscle satellite cells, whereas si-CAMK2B restored the function of miR-188-5p inhibitor. These results suggest that miR-188-5p inhibits the proliferation and promotes the differentiation of goat muscle satellite cells by targeting CAMK2B. This study will provide a theoretical reference for future studies on the molecular mechanisms of skeletal muscle development in goats.

HuR Promotes the Differentiation of Goat Skeletal Muscle Satellite Cells by Regulating Myomaker mRNA Stability.

Human antigen R (HuR) is an RNA-binding protein that contributes to a wide variety of biological processes and diseases. HuR has been demonstrated to regulate muscle growth and development, but its regulatory mechanisms are not well understood, especially in goats. In this study, we found that HuR was highly expressed in the skeletal muscle of goats, and its expression levels changed during longissimus dorsi muscle development in goats. The effects of HuR on goat skeletal muscle development were explored using skeletal muscle satellite cells (MuSCs) as a model. The overexpression of HuR accelerated the expression of myogenic differentiation 1 (MyoD), Myogenin (MyoG), myosin heavy chain (MyHC), and the formation of myotubes, while the knockdown of HuR showed opposite effects in MuSCs. In addition, the inhibition of HuR expression significantly reduced the mRNA stability of MyoD and MyoG. To determine the downstream genes affected by HuR at the differentiation stage, we conducted RNA-Seq using MuSCs treated with small interfering RNA, targeting HuR. The RNA-Seq screened 31 upregulated and 113 downregulated differentially expressed genes (DEGs) in which 11 DEGs related to muscle differentiation were screened for quantitative real-time PCR (qRT-PCR) detection. Compared to the control group, the expression of three DEGs (Myomaker, CHRNA1, and CAPN6) was significantly reduced in the siRNA-HuR group (p < 0.01). In this mechanism, HuR bound to Myomaker and increased the mRNA stability of Myomaker. It then positively regulated the expression of Myomaker. Moreover, the rescue experiments indicated that the overexpression of HuR may reverse the inhibitory impact of Myomaker on myoblast differentiation. Together, our findings reveal a novel role for HuR in promoting muscle differentiation in goats by increasing the stability of Myomaker mRNA.

Transcriptome-Wide Study Revealed That N6-Methyladenosine Participates in Regulation Meat Production in Goats.

In mammals, skeletal muscle development is a complex biological process regulated by many factors. N6-methyladenosine (m6A) RNA modification plays an important role in many biological processes. However, the regulation of m6A on skeletal muscle growth and development in adult goats remains unclear. In this study, Duan goats (DA) and Nubia goats (NBY), both female and 12 months old, were selected as the research objects, and m6A-Seq and RNA-Seq were mainly used to detect the difference of m6A modification and gene expression during the development of the longissimus dorsi (LD) muscle in the two breeds. The results showed that compared with DA, the meat production performance of NBY was better than that of DA, and the modification level of m6A was higher than that of DA in LD. The m6A-Seq of LD indicated m6A peaks were mainly enriched in the coding sequence (CDS) and stop codon. A total of 161 differentially methylated genes (DMGs) and 1294 differentially expressed genes (DEGs) were identified in two breeds. GO and KEGG analysis showed that DMGs were closely related to cellular metabolism, and most of DMGs were enriched in pathways related to energy metabolism, muscle growth and development, mainly MAPK signaling pathway, Wnt signaling pathway and CGMP-PKG signaling pathway. The DEGs were significantly enriched in actin binding, calcium ion binding, angiogenesis, and other biological processes, and most of them were enriched in PI3K-Akt and CGMP-PKG signaling pathways. Combined analysis of m6A-Seq and RNA-Seq data revealed a negative correlation between differentially methylated m6A levels and mRNA abundance, and mRNA expression of the gene with m6A peak near 3'UTR will decrease. In addition, 11 DMGs regulating cell differentiation, muscle growth and development were identified. This study displayed the m6A profiles and distribution patterns in the goat transcriptome, determined the potential role of m6A modification in muscle growth and provided a new reference for the further study of goat skeletal muscle development.

Comparative and Functional Analysis of miRNAs and mRNAs Involved in Muscle Fiber Hypertrophy of Juvenile and Adult Goats.

MicroRNAs (miRNAs) are small non-coding RNAs that post-transcriptionally regulate several pathway intermediates and affect the skeletal muscle development in mice, pigs, sheep, and cattle. However, to date, only a small number of miRNAs have been reported in the muscle development of goats. In this report, the longissimus dorsi transcripts of one- and ten-month-old goats were analyzed by sequencing RNAs and miRNAs. The results showed that the ten-month-old Longlin goats had 327 up- and 419 down-regulated differentially expressed genes (DEGs) compared with the one-month-old. In addition, 20 co-up-regulated and 55 co-down-regulated miRNAs involved in the muscle fiber hypertrophy of goats were identified in ten-month-old Longlin and Nubian goats compared with one-month-old. Five miRNA-mRNA pairs (chi-let-7b-3p-MIRLET7A, chi-miR193b-3p-MMP14, chi-miR-355-5p-DGAT2, novel_128-LOC102178119, novel_140-SOD3) involved in the goat skeletal muscle development were identified by miRNA-mRNA negative correlation network analysis. Our results provided new insight into the functional roles of goat muscle-associated miRNAs, allowing a deeper understanding of the transformation of miRNA roles during mammalian muscle development.

MiR-193b-3p Promotes Proliferation of Goat Skeletal Muscle Satellite Cells through Activating IGF2BP1

As a well-known cancer-related miRNA, miR-193b-3p is enriched in skeletal muscle and dysregulated in muscle disease. However, the mechanism underpinning this has not been addressed so far. Here, we probed the impact of miR-193b-3p on myogenesis by mainly using goat tissues and skeletal muscle satellite cells (MuSCs), compared with mouse C2C12 myoblasts. miR-193b-3p is highly expressed in goat skeletal muscles, and ectopic miR-193b-3p promotes MuSCs proliferation and differentiation. Moreover, insulin-like growth factor-2 mRNA-binding protein 1 (IGF2BP1) is the most activated insulin signaling gene when there is overexpression of miR-193b-3p; the miRNA recognition element (MRE) within the IGF1BP1 3' untranslated region (UTR) is indispensable for its activation. Consistently, expression patterns and functions of IGF2BP1 were similar to those of miR-193b-3p in tissues and MuSCs. In comparison, ectopic miR-193b-3p failed to induce PAX7 expression and myoblast proliferation when there was IGF2BP1 knockdown. Furthermore, miR-193b-3p destabilized IGF2BP1 mRNA, but unexpectedly promoted levels of IGF2BP1 heteronuclear RNA (hnRNA), dramatically. Moreover, miR-193b-3p could induce its neighboring genes. However, miR-193b-3p inversely regulated IGF2BP1 and myoblast proliferation in the mouse C2C12 myoblast. These data unveil that goat miR-193b-3p promotes myoblast proliferation via activating IGF2BP1 by binding to its 3' UTR. Our novel findings highlight the positive regulation between miRNA and its target genes in muscle development, which further extends the repertoire of miRNA functions.

CircTCF4 Suppresses Proliferation and Differentiation of Goat Skeletal Muscle Satellite Cells Independent from AGO2 Binding.

The proliferation and differentiation of mammalian skeletal muscle satellite cells (MuSCs) are highly complicated. Apart from the regulatory signaling cascade driven by the protein-coding genes, non-coding RNAs such as microRNAs (miRNA) and circular RNAs (circRNAs) play essential roles in this biological process. However, circRNA functions in MuSCs proliferation and differentiation remain largely to be elucidated. Here, we screened for an exonic circTCF4 based on our previous RNA-Seq data, specifically expressed during the development of the longest dorsal muscle in goats. Subsequently, the circular structure and whole sequence of circTCF4 were verified using Sanger sequencing. Besides, circTCF4 was spatiotemporally expressed in multiple tissues from goats but strikingly enriched in muscles. Furthermore, circTCF4 suppressed MuSCs proliferation and differentiation, independent of AGO2 binding. Finally, we conducted Poly(A) RNA-Seq using cells treated with small interfering RNA targeting circTCF4 and found that circTCF4 would affect multiple signaling pathways, including the insulin signaling pathway and AMPK signaling pathway related to muscle differentiation. Our results provide additional solid evidence for circRNA regulating skeletal muscle formation.

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

Simple SummaryThe skeletal muscle growth and development affects the production of goat meat. Here, we performed transcriptome sequencing of longissimus dorsi muscle of goats at two development stages (namely, 1-month-old and 9-month-old) by RNA sequencing, assessing the lncRNA expression profile during skeletal muscle development. This study identified regulatory lncRNAs and pathways related to the development of skeletal muscle, providing a reference for future study on the molecular mechanism that regulates the skeletal growth and development.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.

Characteristics and Expression of circ_003628 and Its Promoted Effect on Proliferation and Differentiation of Skeletal Muscle Satellite Cells in Goats.

Simple SummaryCircular RNAs (circRNAs) are new regulators of the development of skeletal muscle in mammals. Herein, circ_003628 in Longissimus dorsi muscle tissue of goats, previously found by RNA-seq, was selected to construct an expression profile in different caprine tissues, and investigate the effect on proliferation and differentiation of caprine skeletal muscle satellite cells (SMSCs), using RT-qPCR, EdU, CCK-8 and immunofluorescence assays. The results showed that circ_003628 had the highest expression level both in the longissimus dorsi muscle among nine caprine tissues collected, and on day 6 after differentiation during SMSCs differentiation periods. The interfering of circ_003628 inhibited the viability, proliferation, and differentiation of goat SMSCs.In our previous a study, circ_003628 was one of the most highly expressed circular RNAs (circRNAs) in the Longissimus dorsi muscle of goats found by RNA-seq, suggesting that the circRNA may be important for caprine muscle growth and development. However, there have been no reports describing the molecular mechanisms by which circ_003628 regulates the activities of goat skeletal muscle satellite cells (SMSCs). In this study, reverse transcriptase-PCR (RT-PCR) and DNA sequencing were used to validate the authenticity of circ_003628, and its characteristics, expression profile and effect on goat SMSCs were also studied using real-time quantitative-PCR (RT-qPCR), EdU, CCK-8 and immunofluorescence assays. Circ_003628 is partially originated from 13 exons, 12 introns and 3′-untranslated regions (UTR) of caprine Myosin Heavy Chain 1 (MYH1), and 25 exons and 5′ UTR of Myosin Heavy Chain 4 (MYH4), as well as intergenic sequences between the two genes. A total of 77.07% of circ_003628 were located in the nuclei of goat SMSCs, while 22.93% were expressed in the cytoplasm. The circRNAs were only expressed in triceps brachii, quadriceps femoris and longissimus dorsi muscle tissues in nine caprine tissues investigated, with the highest expression level in longissimus dorsi muscle. The expression level of circ_003628 gradually increased during differentiation periods of goat SMSCs and reached the maximum on day 6 after differentiation. The small interfering RNA of circ_003628 (named si-circ_003628) inhibited the viability and proliferation of goat SMSCs, and also decreased the expression of four cell proliferation marker genes: paired box 7 (Pax7), cyclin-dependent kinase 2 (CDK2), CDK4 and CyclinD1 in goat SMSCs. Transfection of si-circ_003628 significantly decreased the area of MyHC-labeled myotubes of goat SMSCs, as well as the expression levels of three differentiation marker genes: myosin heavy chain (MyHC), myogenin (MyoG), and myocyte enhancer factor 2C (MEF2C). These results suggest that circ_003628 promotes the viability, proliferation, and differentiation of goat SMSCs, and they also provide an improved understanding of the roles of circ_003628 in skeletal muscle growth and development in goats.

Chi-miR-99b-3p Regulates the Proliferation of Goat Skeletal Muscle Satellite Cells In Vitro by Targeting Caspase-3 and NCOR1.

Simple SummaryIn this study, we investigated the role of chi-miR-99b-3p in goat skeletal muscle satellite cells (SMSCs). We confirmed that chi-miR-99b-3p promoted proliferation through targeting Caspase-3 and nuclear receptor corepressor 1, and inhibiting the intrinsic apoptosis-related genes in SMSCs, whereas inhibition of chi-miR-99b-3p had the opposite effect. Furthermore, integrative transcriptomic analysis revealed that overexpression of chi-miR-99b-3p induced various differentially expressed (DE)-gene-associated processes. In addition, 47 dysregulated miRNAs, including 16 upregulated and 31 downregulated miRNAs, which were involved in the biological pathways associated with the DE genes, were identified. Our study demonstrated that chi-miR-99b-3p was an effective facilitator of goat SMSCs. This study is helpful for future studies in skeletal muscle development.We previously found that chi-miR-99b-3p was highly expressed in the skeletal muscle of 7-month-old (rapid growth period) goats and speculated that it may be associated with muscle development. To further investigate the role of chi-miR-99b-3p in goats, we found that chi-miR-99b-3p acted as a myogenic miRNA in the regulation of skeletal muscle development. Dual-luciferase reporter assays, qRT-PCR, and Western blot results confirmed that Caspase-3 and nuclear receptor corepressor 1 were direct targets for chi-miR-99b-3p as their expression was inhibited by this miR. Cell proliferation and qRT-PCR assays showed that chi-miR-99b-3p promoted proliferation through relevant targets and intrinsic apoptosis-related genes in goat skeletal muscle satellite cells (SMSCs), whereas inhibition of chi-miR-99b-3p had the opposite effect. Furthermore, integrative transcriptomic analysis revealed that overexpression of chi-miR-99b-3p induced various differentially expressed (DE) genes mainly associated with the cell cycle, relaxin signaling pathway, DNA replication, and protein digestion and absorption. Notably, most of the cell-cycle-related genes were downregulated in SMSCs after miR-99b-3p upregulation, including the pro-apoptosis-related gene BCL2. In addition, 47 DE miRNAs (16 upregulated and 31 downregulated) were determined by Small RNA-sequencing in SMSCs after chi-miR-99b-3p overexpression. Based on the KEGG enrichment analysis, we found that these DE miRNAs were involved in the biological pathways associated with the DE genes. Our study demonstrated that chi-miR-99b-3p was an effective facilitator of goat SMSCs and provided new insights into the mechanisms by which miRNAs regulate skeletal muscle growth in goats.