Single-cell RNA-seq reveals novel interaction between muscle satellite cells and fibro-adipogenic progenitors mediated with FGF7 signalling.

BACKGROUND: Muscle satellite cells (MuSCs) exert essential roles in skeletal muscle adaptation to growth, injury and ageing, and their functions are extensively modulated by microenvironmental factors. However, the current knowledge about the interaction of MuSCs with niche cells is quite limited. METHODS: A 10× single-cell RNA sequencing (scRNA-seq) was performed on porcine longissimus dorsi and soleus (SOL) muscles to generate a single-cell transcriptomic dataset of myogenic cells and other cell types. Sophisticated bioinformatic analyses, including unsupervised clustering analysis, marker gene, gene set variation analysis (GSVA), AUCell, pseudotime analysis and RNA velocity analysis, were performed to explore the heterogeneity of myogenic cells. CellChat analysis was used to demonstrate cell-cell communications across myogenic cell subpopulations and niche cells, especially fibro-adipogenic progenitors (FAPs). Integrated analysis with human and mice datasets was performed to verify the expression of FGF7 across diverse species. The role of FGF7 on MuSC proliferation was evaluated through administering recombinant FGF7 to porcine MuSCs, C2C12, cardiotoxin (CTX)-injured muscle and d-galactose ( d-gal)-induced ageing model. RESULTS: ScRNA-seq totally figured out five cell types including myo-lineage cells and FAPs, and myo-lineage cells were further classified into six subpopulations, termed as RCN3+, S100A4+, ID3+, cycling (MKI67+), MYF6+ and MYMK+ satellite cells, respectively. There was a higher proportion of cycling and MYF6+ cells in the SOL population. CellChat analysis uncovered a particular impact of FAPs on myogenic cells mediated by FGF7, which was relatively highly expressed in SOL samples. Administration of FGF7 (10 ng/mL) significantly increased the proportion of EdU+ porcine MuSCs and C2C12 by 4.03 ± 0.81% (P < 0.01) and 6.87 ± 2.17% (P < 0.05), respectively, and knockdown of FGFR2 dramatically abolished the pro-proliferating effects (P < 0.05). In CTX-injured muscle, FGF7 significantly increased the ratio of EdU+/Pax7+ cells by 15.68 ± 5.45% (P < 0.05) and elevated the number of eMyHC+ regenerating myofibres by 19.7 ± 4.25% (P < 0.01). Under d-gal stimuli, FGF7 significantly reduced γH2AX+ cells by 17.19 ± 3.05% (P < 0.01) in porcine MuSCs, induced EdU+ cells by 4.34 ± 1.54% (P < 0.05) in C2C12, and restored myofibre size loss and running exhaustion in vivo (all P < 0.05). CONCLUSIONS: Our scRNA-seq reveals a novel interaction between muscle FAPs and satellite cells mediated by FGF7-FGFR2. Exogenous FGF7 augments the proliferation of satellite cells and thus benefits muscle regeneration and counteracts age-related myopathy.

MicroRNA-542-3p targets Pten to inhibit the myoblasts proliferation but suppresses myogenic differentiation independent of targeted Pten.

BACKGROUND: Non-coding RNA is a key epigenetic regulation factor during skeletal muscle development and postnatal growth, and miR-542-3p was reported to be conserved and highly expressed in the skeletal muscle among different species. However, its exact functions in the proliferation of muscle stem cells and myogenesis remain to be determined. METHODS: Transfection of proliferative and differentiated C2C12 cells used miR-542-3p mimic and inhibitor. RT-qPCR, EdU staining, immunofluorescence staining, cell counting kit 8 (CCK-8), and Western blot were used to evaluate the proliferation and myogenic differentiation caused by miR-542-3p. The dual luciferase reporter analysis and rescued experiment of the target gene were used to reveal the molecular mechanism. RESULTS: The data shows overexpression of miR-542-3p downregulation of mRNA and protein levels of proliferation marker genes, reduction of EdU+ cells, and cellular vitality. Additionally, knocking it down promoted the aforementioned phenotypes. For differentiation, the miR-542-3p gain-of-function reduced both mRNA and protein levels of myogenic genes, including MYOG, MYOD1, et al. Furthermore, immunofluorescence staining immunized by MYHC antibody showed that the myotube number, fluorescence intensity, differentiation index, and myotube fusion index all decreased in the miR-542-3p mimic group, compared with the control group. Conversely, these phenotypes exhibited an increased trend in the miR-542-3p inhibitor group. Mechanistically, phosphatase and tensin homolog (Pten) was identified as the bona fide target gene of miR-542-3p by dual luciferase reporter gene assay, si-Pten combined with miR-542-3p inhibitor treatments totally rescued the promotion of proliferation by loss-function of miR-542-3p. CONCLUSIONS: This study indicates that miR-542-3p inhibits the proliferation and differentiation of myoblast and Pten is a dependent target gene of miR-542-3p in myoblast proliferation, but not in differentiation.

The long noncoding RNA lnc-H19 is important for endurance exercise by maintaining slow muscle fiber types.

Skeletal muscle consists of different muscle fiber types whose heterogeneity is characterized by different metabolic patterns and expression of MyHC isomers. The transformation of muscle fiber types is regulated by a complex molecular network in which long noncoding (lnc) RNAs play an important role. In this study, we found that lnc-H19 is more enriched in slow muscle fibers. In vitro, interference of lnc-H19 by siRNA significantly promoted the expression of fast muscle fiber gene MyHC IIB and inhibited the expression of the slow muscle fiber gene MyHC I, thereby leading to a fast muscle fiber phenotype. In addition, interference of lnc-H19 significantly inhibited mRNA expression of the mitochondrial genes, such as COX5A, COX-2, UQCRFSL, FABP3, and CD36. Overexpression of lnc-H19 resulted in an opposite result. In vivo, knockdown of lnc-H19 by AAV-shRNA-H19 suppressed the mRNA expression of the slow muscle fiber gene MyHC I and the protein expression of slow-MyHC. Simultaneously, mitochondria were reduced in number, swollen, and vacuolated. The activities of succinate dehydrogenase, lactic dehydrogenase, and superoxide dismutase were significantly inhibited, and malondialdehyde content was significantly increased, indicating that deficiency of lnc-H19 leads to decreased oxidative metabolism and antioxidant capacity in muscle. Furthermore, inhibition of lnc-H19 decreased the weight-bearing swimming time and limb suspension time of mice. In conclusion, our results revealed the role of lnc-H19 in maintaining slow muscle fiber types and maintaining exercise endurance, which may help to further improve the regulatory network of lnc-H19 in muscle function.

Lnc-Malat1 promotes slow myofiber-type transformation through sponging miR-129-5p in C2C12 myotubes.

Long non-coding metastasis-associated lung adenocarcinoma transcript (lnc-Malat1) emerges as a novel regulator in skeletal muscle development, while its function and the related mechanism is not fully revealed yet. In this study, knockdown of lnc-Malat1 by siRNA significantly inhibited the expression of myoblast marker genes (MyHC, MyoD, and MyoG) and slow muscle fiber marker genes (MyHC I), together with repressed expression of mitochondria-related genes COX5A, ACADM, CPTA1, FABP3, and NDUFA1. Overexpression of lnc-Malat1 exerted an opposite effect, promoting myoblast differentiation and slow muscle fiber formation. Dual luciferase reporter assay revealed a direct interaction between lnc-Malat1 and miR-129-5p, and overexpression of lnc-Malat1 significantly inhibited miR-129-5p expression, thereby elevating the expression of Mef2a, miR-129-5p target protein. In addition, enforced expression of lnc-Malat1 restored the inhibitory effect of miR-129-5p on myoblast differentiation and MyHC I expression. Taken together, our results suggest that lnc-Malat1 promotes myoblast differentiation, and maintains the slow muscle fiber phenotype via adsorbing miR-129-5p.

miR-424(322)-5p targets Ezh1 to inhibit the proliferation and differentiation of myoblasts.

The proliferation and differentiation of myoblasts are considered the key biological processes in muscle development and muscle-related diseases, in which the miRNAs involved remain incompletely understood. Previous research reported that miR-424(322)-5p is highly expressed in mouse skeletal muscle. Therefore, C2C12 cells are used as a model to clarify the effect of miR-424(322)-5p on the proliferation and differentiation of myoblasts. The data show that miR-424(322)-5p exhibits a decreasing trend upon myogenic differentiation. Overexpression of miR-424(322)-5p inhibits the proliferation of myoblasts, manifested by downregulation of proliferation marker genes ( CCNB1, CCND2, and CDK4), decreased percentage of EdU + cells, and reduced cell viability. In contrast, these phenotypes are promoted in myoblasts treated with an inhibitor of miR-424(322)-5p. Interestingly, its gain of function inhibits the expression of myogenic regulators, including MyoD, MyoG, MyHC, and Myf5. Additionally, immunofluorescence staining of MyHC and MyoD shows that overexpression of miR-424(322)-5p reduces the number of myotubes and decreases the myotube fusion index. Consistently, inhibition of its function mediated by an inhibitor promotes the expressions of myogenic markers and myotube fusion. Mechanistically, gene prediction and dual-luciferase reporter experiments confirm that enhancer of zeste homolog 1 ( Ezh1) is one of the targets of miR-424(322)-5p. Furthermore, knockdown of Ezh1 inhibits the proliferation and differentiation of myoblasts. Compared with NC and inhibitor treatment, inhibitor+si- EZH1 treatment rescues the phenotypes of proliferation and differentiation mediated by the miR-424(322)-5p inhibitor. Taken together, these data indicate that miR-424(322)-5p targets Ezh1 to negatively regulate the proliferation and differentiation of myoblasts.

Molecular characterization of TRIB1 gene and its role in regulation of steroidogenesis in bos grunniens granulosa cells.

To explore the expression pattern of the TRIB1 gene in yak follicles and its effect on the steroidogenesis of granulosa cells (GCs). Here, 4-5 years old female yaks were treated as the subjects. Immunohistochemically assay found that TRIB1 protein was expressed in different developmental follicles. Among different cell types of follicles, including cumulus cells (CCs), granulosa cells (GCs) and theca cells (TCs), the TRIB1 protein was most abundant in GCs (P < 0.0001). In addition, we cloned the coding sequence (CDS) of the yak TRIB1 gene, which is 1119 bp, encoding 372 amino acids (AA). The amino acid sequence homology of TRIB1 is >80% to those of other species, except for zebrafish. To further explore the function of TRIB1 in steroidogenesis, the pcDNA3.1(+)-TRIB1 eukaryotic expression vector was constructed and then transfected into GCs. The data showed that overexpression of TRIB1 significantly reduced the progesterone (P4) secretion of granulosa cells measured by ELISA assay (P < 0.05), but not Estradiol (E2) secretion. Consistently, TRIB1 gain-of-function downregulated the mRNA levels of steroidogenesis related genes steroidogenic acute regulatory protein (StAR), cytochrome P450 family 11 subfamily A member 1 (CYP11A1) and 3ß-hydroxysteroid dehydrogenase (3ß-HSD) (P < 0.01), while cytochrome P450 family 17 subfamily A member 1 (CYP17A1) and cytochrome P450 family 19 subfamily A member 1 (CYP19A1) had no significant difference (P > 0.05). Interestingly, mito-tracker staining showed that mitochondrial number significantly decreased in TRIB1 overexpressed GCs (P < 0.01). Further, overexpression of TRIB1 inhibited mRNA levels of mitochondrial biogenesis related genes, including Mitochondrial transcription factor (TFAM) and Peroxisome proliferator-activated receptor alpha co-activator (PPARGC1A) (P < 0.05). Conclusively, this work indicates that TRIB1 inhibited progesterone synthesis of GCs might be involved in the reduction of the mitochondria number.

LKB1 Regulates Goat Intramuscular Adipogenesis Through Focal Adhesion Pathway.

Intramuscular fat (IMF) deposition is one of the most important factors to affect meat quality in livestock and induce insulin resistance and adverse metabolic phenotypes for humans. However, the key regulators involved in this process remain largely unknown. Although liver kinase B1 (LKB1) was reported to participate in the development of skeletal muscles and classical adipose tissues. Due to the specific autonomic location of intramuscular adipocytes, deposited between or within muscle bundles, the exact roles of LKB1 in IMF deposition need further verified. Here, we cloned the goat LKB1 coding sequence with 1,317 bp, encoding a 438 amino acid peptide. LKB1 was extensively expressed in detected tissues and displayed a trend from decline to rise during intramuscular adipogenesis. Functionally, knockdown of LKB1 by two individual siRNAs enhanced the intramuscular preadipocytes differentiation, accompanied by promoting lipid accumulation and inducing adipogenic transcriptional factors and triglyceride synthesis-related genes expression. Conversely, overexpression of LKB1 restrained these biological signatures. To further explore the mechanisms, the RNA-seq technique was performed to compare the difference between siLKB1 and the control group. There were 1,043 differential expression genes (DEGs) were screened, i.e., 425 upregulated genes and 618 downregulated genes in the siLKB1 group. The Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis predicted that the DEGs were mainly enriched in the focal adhesion pathway and its classical downstream signal, the PI3K-Akt signaling pathway. Specifically, knockdown of LKB1 increased the mRNA level of focal adhesion kinase (FAK) and vice versa in LKB1-overexpressed cells, a key component of the activated focal adhesion pathway. Convincingly, blocking this pathway by a specific FAK inhibitor (PF573228) rescued the observed phenotypes in LKB1 knockdown adipocytes. In conclusion, LKB1 inhibited goat intramuscular adipogenesis through the focal adhesion pathway. This work expanded the genetic regulator networks of IMF deposition and provided theoretical support for improving human health and meat quality from the aspect of IMF deposition.

Ablation of KDM2A Inhibits Preadipocyte Proliferation and Promotes Adipogenic Differentiation.

Epigenetic signals and chromatin-modifying proteins play critical roles in adipogenesis, which determines the risk of obesity and which has recently attracted increasing interest. Histone demethylase 2A (KDM2A) is an important component of histone demethylase; however, its direct effect on fat deposition remains unclear. Here, a KDM2A loss of function was performed using two unbiased methods, small interfering RNA (siRNA) and Cre-Loxp recombinase systems, to reveal its function in adipogenesis. The results show that the knockdown of KDM2A by siRNAs inhibited the proliferation capacity of 3T3-L1 preadipocytes. Furthermore, the promotion of preadipocyte differentiation was observed in siRNA-treated cells, manifested by the increasing content of lipid droplets and the expression level of adipogenic-related genes. Consistently, the genetic deletion of KDM2A by Adipoq-Cre in primary adipocytes exhibited similar phenotypes to those of 3T3-L1 preadipocytes. Interestingly, the knockdown of KDM2A upregulates the expression level of Transportin 1(TNPO1), which in turn may induce the nuclear translocation of PPARγ and the accumulation of lipid droplets. In conclusion, the ablation of KDM2A inhibits preadipocyte proliferation and promotes its adipogenic differentiation. This work provides direct evidence of the exact role of KDM2A in fat deposition and provides theoretical support for obesity therapy that targets KDM2A.