Ryanodine receptor RyR1-mediated elevation of Ca2+ concentration is required for the late stage of myogenic differentiation and fusion
Journal of Animal Science and Biotechnology volume 13, Article number: 9 (2022)
Cytosolic Ca2+ plays vital roles in myogenesis and muscle development. As a major Ca2+ release channel of endoplasmic reticulum (ER), ryanodine receptor 1 (RyR1) key mutations are main causes of severe congenital myopathies. The role of RyR1 in myogenic differentiation has attracted intense research interest but remains unclear.
In the present study, both RyR1-knockdown myoblasts and CRISPR/Cas9-based RyR1-knockout myoblasts were employed to explore the role of RyR1 in myogenic differentiation, myotube formation as well as the potential mechanism of RyR1-related myopathies. We observed that RyR1 expression was dramatically increased during the late stage of myogenic differentiation, accompanied by significantly elevated cytoplasmic Ca2+ concentration. Inhibition of RyR1 by siRNA-mediated knockdown or chemical inhibitor, dantrolene, significantly reduced cytosolic Ca2+ and blocked multinucleated myotube formation. The elevation of cytoplasmic Ca2+ concentration can effectively relieve myogenic differentiation stagnation by RyR1 inhibition, demonstrating that RyR1 modulates myogenic differentiation via regulation of Ca2+ release channel. However, RyR1-knockout-induced Ca2+ leakage led to the severe ER stress and excessive unfolded protein response, and drove myoblasts into apoptosis.
Therefore, we concluded that Ca2+ release mediated by dramatic increase in RyR1 expression is required for the late stage of myogenic differentiation and fusion. This study contributes to a novel understanding of the role of RyR1 in myogenic differentiation and related congenital myopathies, and provides a potential target for regulation of muscle characteristics and meat quality.
Ryanodine receptor 1 (RyR1), located on the endoplasmic/sarcoplasmic reticulum (ER/SR) membrane, is highly expressed in the skeletal muscle. It serves as a critical Ca2+ channel in mediating intracellular flux of Ca2+, triggering the contraction of the skeletal muscle. Various mutations or epigenetic changes in the RyR1 gene have been demonstrated to associate with muscle myopathies including malignant hyperthermia and several congenital myopathies [1,2,3,4]. RyR1 crystal structure has been resolved using electron cryomicroscopy as a 6-transmembrane ion channel with an EF-hand domain for Ca2+-mediated allosteric gating and a huge cytoplasmic domain on top of each transmembrane domain [5,6,7]. These studies contribute to the deep understanding of the structure, function and channel activity of RyR1. In particular, the impacts of two malignant hyperthermia-associated mutations on the RyR1 3D structure were recently solved to gain insights into the pathogenesis [8, 9]. Diagnostic gene-sequencing has been employed to avoid RyR1-related congenital myopathies, and several disease-modulating therapeutic strategies and salvage therapies have been developed against RyR1-related myopathies [10,11,12,13,14,15].
It has been well known that cytoplasmic Ca2+ mediates myogenic differentiation and skeletal muscle development . As a major Ca2+ release channel of endoplasmic reticulum (ER), RyRs-mediated Ca2+ release plays a role in the morphogenesis of mammalian skeletal muscle, while most of RyR1 mutations present gain-of-function phenotype and result in leaking of the internal Ca2+ store . Notably, the proportion of pigs carrying the single-base mutation of RyR1 unexpectedly increased during intensive breeding for lean meat in pigs, which leads to the increased yield of abnormal meat characteristics of pale, soft and exudative . However, the block of RyRs activity by ryanodine selectively retards fetal myoblast differentiation, and lead to both physiological and pathological consequences of muscle morphology and functions . It has been demonstrated that homozygous RyR1-null mice died after birth and displayed small limbs and abnormal skeletal muscle organization [20, 21]. Therefore, we hypothesize that RyR1 is not only involved in causing congenital myopathies, but implicated in myogenesis and subsequent muscle development. Therefore, the mechanisms of RyR1 action in myogenesis need to be elucidated.
Endoplasmic/sarcoplasmic reticulum is responsible for proper folding, processing, and trafficking of proteins and plays an important role in cellular Ca2+ homeostasis. Alterations in cellular Ca2+ dynamics directly trigger ER stress and activate the unfolded protein response (UPR) [22, 23]. ER stress and resultant UPR modulation are implicated in various human diseases including sarcopenia [24,25,26], and roles of ER stress and UPR pathways in skeletal muscle health and diseases receive increased research attention . Therefore, we surmise that ER stress signaling is involved in RyR1-related muscle myopathies.
Myogenic differentiation and fusion play a pivotal role in the development and skeletal muscle maturation via forming mature multinuclear myofibers, impacting on meat quality in livestock . In this study, we employed CRISPR/Cas9-based RyR1-KO myoblasts and siRNA-mediated RyR1-knockdown myoblasts to explore the role of RyR1 in myogenesis and formation mechanism of RyR1-related myopathies.
Materials and methods
Cell culture and myogenic differentiation
Mouse skeletal myoblast C2C12 cells were purchased from the National Infrastructure of Cell Line Resource in China. Proliferating myoblasts were maintained in DMEM/high glucose medium (Hyclone, Logan, UT, USA) supplemented with 10% FBS (Gibco, Carlsbad, CA, USA) in a humidified CO2 incubator (5% CO2, 37 °C; HF90, Heal Force, Hongkong, China). For myogenic differentiation, myoblasts with 80% ~ 90% confluence were induced by DMEM/high glucose medium containing 2% horse serum (Hyclone, Logan, UT, USA).
Myogenic cells of pigs (n = 3) were isolated using preplate techniques from the skeletal muscle according to the protocol our lab established previously [29, 30]. The cells were cultured in growth medium in a 100 mm dish coated with collagen I (Sigma-Aldrich, Louis, MO, USA) at 37 °C and 5% CO2. The growth medium was composed of DMEM/F12 (Hyclone), 10% FBS (Gibco-BRL, Carlsbad, CA, USA), 2 mmol/L glutamine (Gibco-BRL), and 5 ng/mL bFGF (Peptech, Burlington, MA, USA). As for myogenic induction, cells were cultured for 5 d in DMEM/F12 medium containing 2% horse serum (Hyclone).
Cellular Ca2+ concentration measurement
Ca2+ concentration in the cytoplasm or ER was measured using flow cytometry. Briefly, cells were collected, washed with PBS (phosphate buffered saline) and HBSS (Hanks balanced salt solutions) subsequently, and then incubated in 5 μg/mL Fluo-3 acetoxymethyl ester (Cayman, Ann Arbor, MI, USA) or Mag-fluo-AM (GENMED, Shanghai, China) for 30 min at 37 °C in dark. After three washes with PBS supplemented with 1% FBS, cells were resuspended in 200 μL PBS containing 1% FBS. Flow cytometry was carried out immediately using a FACS Calibur Cytometer and Image Cytometry software (BD, Franklin, NJ, USA). Calcium-bound Fluo-3 or Mag-fluo-AM has an emission maximum of 526 nm which was quantified by excitation with a 488-nm laser and signals were collected using a 530/30 nm band-pass filter. Each sample generated 20,000 live gated events. Debris, multicellularity, and dead cells were excluded by forward scatter (FSC) and side scatter. For detecting dynamic change of Ca2+ concentration of unfused cells during myogenic differentiation, a blank control combined with a house-keeping control (the proliferative C2C12 cells) was used to correct the deviation caused by the loaded indicator amount and the voltage used in each measurement. Mean fluorescence intensity was determined from the entire cell population and then adjusted by relative cell size calculated according to FSC to represent Ca2+ concentration.
Cell viability and apoptosis assays
Cell viability was measured using Cell Counting Kit-8 (CA1210, Solarbio, Beijing, China) according to the manufacturer’s guideline. Briefly, cells were cultured in 96-well plates for 24 h, and then CCK-8 reagent was added at 100 μL per well. One hour later, the absorbance of culture medium was analyzed by microplate spectrophotometer. In addition, cell proliferation activity was also measured by Cell-Light™ EdU Apollo®488 Cell Tracking Kit (RIBOBIO, Guangzhou, China). After pre-cultured for 24 h in 96-well plates, cells were cultured continuously for another 2 h in new media supplemented with 50 μmol/L EdU reagent. Then cells were fixed by 4% paraformaldehyde, permeabilized by 0.2% Trutib X-100, and fluorescently-tagged with Hoechst3342 using nucleus staining methods. The newly proliferated cells were visualized by an Apollo reaction system. Cell proliferation rate was analyzed by ImageJ (v1.51h, National Institutes of Health, Bethesda, MD, USA).
Apoptosis was tested using an Annexin V-FITC/PI Apoptosis Detection Kit (Gene Protein Link, Beijing, China) according to the manufacturer’s protocol. Briefly, cells were stained with a combination of Annexin V-FITC and propidium iodide in darkness for 15 min at room temperature, and then analyzed by the flow cytometry system.
RNA isolation and qRT-PCR
Cells used for total RNA extraction obtained from 3 separate experiments (different batches of cells and on different days). Total RNA was extracted from cells using HiPure Total RNA Mini Kit (Magen, Beijing, China), and then reverse-transcribed into cDNA using a PrimeScript™ RT reagent Kit with gDNA Eraser (Takara, Osaka, Japan). Synthesized cDNA was used for RT-qPCR analysis by employing a quantitative real-time PCR kit (Takara, Osaka, Japan) with an AJ qTOWER 2.2 Real-Time PCR system (Analytik Jena AG, Jena, Germany) according to standard procedures. All samples were measured in triplicate. The primers used in the experiment were listed in Table S1. To compare the mRNA expression of RyR1 and RyR3 in cells, the amplification efficiency of their primers was used to rectify the qRT-PCR cycle number. GAPDH (glyceraldehyde-3-phosphate dehydrogenase) was used as an internal control. Relative gene expression level was calculated by 2−ΔΔCt method .
Protein extraction and western blot analysis
The relative abundances of proteins concerning ER stress, MAPK signaling pathway, and apoptosis were determined by Western Blot. Cell samples were collected and lysed in RIPA buffer (Huaxingbio, Beijing, China) composed of 50 mmol/L Tris-HCl (pH 7.4), 150 mmol/L NaCl, 1% NP-40, and 0.1% SDS, plus a Halt protease/phosphatase inhibitor cocktail (Thermo Fisher Scientific, Waltham, MA, USA). The homogenate was centrifuged at 14,000 × g for 15 min at 4 °C and the supernatant was isolated for Western Blot analysis. Protein concentrations were determined using a BCA Protein Assay Kit (Huaxingbio, Beijing, China). Equal amounts of protein (30 μg), together with a pre-stained protein ladder (Thermo Fisher Scientific, Waltham, MA, USA), were electrophoresed on SDS polyacrylamide gel, electro- transferred to a polyvinylidene difluoride membrane (Millipore, Bedford, OH, USA), and blocked for 1 h in 5% non-fat dry milk at room temperature in Tris-Buffered saline and Tween-20 (TBST; 20 mmol/L Tris-Cl, 150 mmol/L NaCl, 0.05% Tween 20, pH 7.4). Samples were incubated with corresponding primary antibodies overnight at 4 °C. After washing with TBST (pH 7.4), membranes were incubated with the secondary antibody (DyLight 800, Goat Anti-Rabbit IgG). Protein bands were detected with the Odyssey Clx kit (LI-COR, Lincoln, NE, USA) and quantified using an Alpha Imager 2200 (Alpha InnoTec, CA, USA). GAPDH was taken as an internal standard to calculate relative protein expression. Antibodies used for Western Blot in this study were listed in Table S2.
Cells were fixed with 4% paraformaldehyde (PFA)/PBS for 30 min. After the neutralization of excess formyl group by 2 mg/mL glycine, cells were permeabilized by 0.2% Trutib X-100 in PBS for 10 min. After blocked with 3% BSA/PBS, cells were incubated with primary antibody (anti-RyR1, 1:300, MA3–925, Thermo Fisher Scientific, Waltham, MA, IL, USA; anti-myosin, 1:300, M4276, Sigma-Aldrich, Louis, MO, USA) overnight and then incubated with Fluorescein-Conjugated secondary antibody (ZF-0311, ZSGB-BIO, Beijing, China) at 1:100. Nuclei were stained with DAPI (Thermo Fisher Scientific, Waltham, MA, USA). Finally, myotubes were visualized by an inverted fluorescence microscope.
Chemical blockers of Ca2+ channels
Two kinds of chemical blocker were employed in the experimental treatments. Dantrolene (DAN) (Stock solution: 200 mmol/L in DMSO, HY-12542A, MedChemExpress, South Brunswick, NJ, USA), an inhibitor of RyR1, was added as a final concentration of 10 μmol/L in culture media. Thapsigargin (THA, ab120286, Abcam, Cambridge, UK), an ER stress-inducing agent, was used at a final concentration of 100 nmol/L in culture media. Equal amounts of vehicle (DMSO) were used as the control. During the proliferating period, myoblasts were treated with chemical blockers, respectively, for 48 h and then collected for further analysis. Upon myogenic induction, myoblasts were treated with chemical blockers, respectively, for 5 d and then used for mRNA extraction and immunocytochemistry. Each treatment was conducted in three independently repeated experiments.
Small interfering RNA transfection
RNA interference of RyR1 (mouse, gene ID: 20190) expression was performed using a 21-base pair small interfering RNA (siRNA) duplex (designed and synthesized by IBSBIO, Shanghai, China). The sense strand nucleotide sequence for RyR1 siRNA was 5′-CCUGCUCUAUGAACUUCUAGC-3′ (sense strand) and 5′-UAGAAGUUCAUAGAGCAGGUU-3′ (anti-sense strand). A scrambled siRNA (siControl, sense strand: 5′-UUCUCCGAACGUGUCACGUTT-3′, anti-sense strand: 5′-ACGUGACACGUUCGGAGAATT-3′) with the same nucleotide composition as RyR1 siRNA but lacks significant sequence homology to the RyR1 was also designed as a negative control. Briefly, myoblasts were plated in a 100 mm cell culture dish for 24 h, and then transfected with 100 nmol/L siRNA using 16 μL Lipofectamine 3000 (Invitrogen, Carlsbad, CA, USA) in each dish. After transfection for 24 h, myogenic differentiation was induced in cells.
Gene-edited myoblasts with RyR1-knockout were generated via the Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-CRISPR associated protein 9 (Cas9) system. The plasmid vectors expressing Cas9 protein and guide-RNA (gRNA) were designed and synthesized by Syngentech (Beijing, China). One gRNA among the three gRNA targeting RyR1 was selected for using in further study according to their shearing efficiency. The sequences of gRNA-1, gRNA-2, and gRNA-3 are as follows: 5′-GGCGATGATCTCTATTCTTA-3′, 5′-TACAGCCCCTACCCCGGAGG-3′, and 5′-AGCTCAGGCCACCCACCTGA-3′, respectively. When the myoblasts were grown to 80% ~ 90% confluency, 1 × 106 cells were collected and transfected with 2 μg CRISPR-Cas9 vectors by electroporation using Nucleofector Program B-032 (VCA-1003, Lonza, Basel, Switzerland). Infected myoblasts were selected by incubation with 200 μg/mL hygromycin for 1 week. The stable RyR1-knockout cell line was obtained through single cell clone techniques. Genotype identification and verification of putative off-target sites (Table S3) were conducted via DNA-sequencing technology. The related primers were list in Table S4.
Student’s t-test was used between two groups; one-way or two-way ANOVA with Tukey’s test was used among multiple groups. Data are presented as mean ± SEM. The criterion for statistical significance was set at P < 0.05.
Cytoplasmic Ca2+ dynamics and expression patterns of Ca2+ channels during myogenic differentiation
Upon myogenic induction, myoblasts gradually expressed plenty of myosin from the d 0 to 6 (Fig. S1A). During the entire period (d 1–5) of myogenic differentiation, cytoplasmic Ca2+ concentration (labeled by Fluo-3) of myoblasts was significantly increased (Fig. S1B and C). Myf5 (myogenic factor 5) and MyoD1 (myogenic differentiation 1) expression were significantly increased on d 2 relative to d 0, then sharply decreased on d 6 and 8 and even lower than the initial level before myogenic induction (Fig. 1A and B). Meanwhile, MyoG (myogenin) expression showed continuous increase and reached a plateau on d 4 (Fig. 1C).
As for Ca2+ transporters, CAV1.1 (also known as CACNA1S, calcium voltage-gated channel subunit alpha1 S), CRACR2B (calcium release activated channel regulator 2B), ITPR1 (inositol 1,4,5-trisphosphate receptor type 1), and ORAI2 (ORAI calcium release-activated calcium modulator 2) whose expression patterns showed similar with Myf5 and MyoD1 (Fig. 1D-G). The expression patterns of RyR1 and STIM1 (stromal interaction molecule 1) resembled MyoG well (Fig. 1H and I). Notably, RyR1 mRNA expression increased more than 100-fold in myoblastic C2C12 cells upon myogenic induction. Consistently, RyR1 protein expression was significantly increased upon myogenic induction (Fig. S1D and E). In addition, the mRNA expression of ATP2A2 (ATPase ER/SR Ca2+ transporting 2; Fig. 1J), ATP2B (ATPase plasma membrane Ca2+ transporting 1; Fig. 1K) and CRACR2A (calcium release activated channel regulator 2A; Fig. 1L) were significantly increased in myoblastic C2C12 cells on d 2 during myogenic induction and maintained a plateau for subsequent days.
Since RyR3 mRNA expression level was less than 1/10 that of RyR1 in C2C12 cells, the increase in RyR3 mRNA expression level during myogenic induction was far below that of RyR1 (Fig. S2A and B). From this view, therefore, RyR1 rather than RyR3, another member of RYRs expressed in skeletal muscle, may exert a major role in mediating myogenesis of myoblastic C2C12 cells can be attributed to. Consistently, we also observed that RyR1 mRNA expression showed almost a 25-fold increase during myogenic differentiation, while the mRNA expression of RyR3 was not changed in myogenic cells of pigs (Fig. S2C-E).
RyR1-mediated elevation of cytoplasmic Ca2+ concentration is indispensable for myogenesis
As shown in Fig. 2A and B, cytoplasmic Ca2+ concentration was significantly decreased upon treatment with dantrolene (DAN), an inhibitor of RyR1, which blocks the release of Ca2+ from SR/ER . Moreover, RyR1 knockdown by siRNA significantly also decreased cytoplasmic Ca2+ concentration of C2C12 cells (Fig. 2C and D). Consistently, the protein expression of RyR1 was effectively suppressed by siRNA interference (Fig. 2E and F). Functional constraints of RyR1 by either DAN or siRyR1 dramatically inhibited the formation of multinucleated myotubes (Fig. 2G and H). During myogenic induction, DAN effectively blocked MyoD1 expression on d 2 and both MyoG and Mymk (myomaker, myoblast fusion factor) expression on d 4 (Fig. 2I). At the d 4 during myogenic induction, siRyR1-knockdown significantly reduced RyR1, MyoG, and Mymk expression without impact on Myf5 and MyoD1 expression (Fig. 2J).
Thapsigargin (THA) treatment significantly increased cytoplasmic Ca2+ concentration (Fig. 3A and B), but did not influence the mRNA expressions of myogenic-specific genes including Myf5, MyoD1, MyoG and Mymk (Fig. 3C). RyR1-knockdown by siRNA interference was independent of THA treatment (P < 0.01 for siRNA treatment, Fig. 3D and E). On d 4 during myogenic induction, siRyR1-knockdown significantly increased Myf5 and MyoD1 mRNA expression, while THA effectively eliminated alterations in Myf5 and MyoD1 mRNA expressions induced by siRyR1-knockdown (P < 0.05 for siRNA treatment × THA treatment, Fig. 4A). Accordingly, myotube formation was significantly inhibited by siRyR1-knockdown, which was effectively recovered by THA (P < 0.01 for siRNA treatment × THA treatment, Fig. 4B-D).
Effects of RyR1-knockout on cell proliferation and differentiation
RyR1 was successfully knockout by CRISPR/Cas-9 gene editing system targeting Exon 18 of RyR1 via gRNA-3 which demonstrated the highest shearing efficiency among three designed gRNAs and resulted in the frame-shift mutation of RyR1 (Fig. S3). Homozygote and heterozygote of RyR1-knockout cells, named as RyR1−/− and RyR1+/−, were obtained by monoclonal cultivation and identified via gene sequencing on the target site and putative off-target sites of gRNA-3. Relative to the wild type cells (WT), RyR1−/− showed the higher concentration of Ca2+ in the cytoplasm, but the lower level of Ca2+ in the ER (Fig. 5A-C). In RyR1−/− or RyR1+/−, the mRNA expression of ATP2A2, ATP2B, CRACR2B, and ORAI1 was significantly increased, while the mRNA expression of CAV1.1 and ORAI2 was decreased relative to WT (Fig. 5D). Cell proliferation (the proportion of EdU+ cells) of RyR1−/− or RyR1+/− significantly declined relative to WT (Fig. 5E and F), which was also demonstrated by the CCK-8 test (Fig. 5G). The mRNA expression of Myf5, MyoD1, and MyoG was significantly increased in RyR1−/− or RyR1+/− relative to WT (Fig. 5H).
RyR1-knockout triggered apoptosis
On d 2 during myogenic induction, apoptosis instead of myotube formation was significantly accelerated in both RyR1−/− and RyR1+/− relative to WT (Fig. 6A and B). The protein expression level of cyclin D1 and caspase-3 was significantly lowered in RyR1−/− or RyR1+/−, but the expression of cleaved caspase-3 was significantly increased as compared with WT (Fig. 6C and D).
ER stress level of cells was evaluated, and both phosphorylated and total protein expression of IRE1α and PERK were dramatically elevated in RyR1−/− and RyR1+/−, while phosphorylation of EIF2α was significantly decreased (Fig. 6E and F). As to several apoptosis-related proteins, the protein abundance of CHOP (also known as DDIT3), caspase-9, and caspase-12 in RyR1−/− or RyR1+/− was significantly increased relative to WT, while the protein expression of ERP44 and HSPA5 was not influenced (Fig. 6G and H).
Skeletal muscle mass is maintained by myogenic differentiation of myogenic progenitors and subsequent myoblast fusion , in which cytosolic Ca2+ dynamics plays a vital role . Particularly, myotube formation requires net Ca2+ influx into myoblasts [34,35,36]. RyR1, serving as a major Ca2+ release channel of ER, has attracted the most intense research interests. However, the role of RyR1 in myogenesis remains unclear.
RyR1 protein expression pattern during myogenic differentiation of myoblastic C2C12 cells provides a well-controlled model for investigations of RyR1 function during myogenic differentiation . In the current study, we observed that cytoplasmic Ca2+ concentration of C2C12 myoblasts was significantly elevated during myogenic differentiation, which was consistent with our previously study in primary myogenic cells of pigs . Cytoplasmic Ca2+ dynamics is tightly regulated by various channels and transports in cells . Relative to the extracellular matrix and ER, cytoplasmic Ca2+ concentration is maintained at very low levels (10–100 nmol/L) under resting conditions. Ca2+ is released from the ER, the main storage site of intracellular Ca2+, through the transmembrane channels RyR1 and ITPR1 [39, 40]. Cytoplasmic Ca2+ influx from extracellular matrix occurs through plasma membrane channels, such as CAV1, CAV2, and CAV3 . In maintaining the resting state, excessive amounts of cytoplasmic Ca2+ re-accumulates in the ER by SR/ER Ca2+-ATPase (SERCA, also called ATP2A)  or is expulsed in the external milieu by plasma membrane Ca2+-ATPase (PMCA, also called ATP2B) [43, 44]. Furthermore, store-operated calcium entry (SOCE), mediated by STIM (ER Ca2+ sensors), activates CRAC and ORAI channels located at plasma membrane to maintain cellular Ca2+ homeostasis [45, 46].
It has been well known that lineage commitment and differentiation of myoblasts are governed by the programmed expression and functional activation of myogenic regulatory transcription factors (MRFs) [47, 48]. In this study, we discovered RyR1 mRNA expression was increased more than 100-fold in C2C12 cells and almost 25-fold in myogenic cells of pigs during myogenic differentiation along with the normal expression pattern of Ca2+ transporters as well as MRFs during myogenic differentiation. In addition, the expression patterns of RyR1 showed similar with MyoG, which was increased in the late stage of myogenesis. It strongly indicates that RyR1 is not necessary for initial myogenic commitment, but more for later stage of myogenic differentiation.
Notably, RyR1 restriction via either DAN treatment or siRNA interference significantly decreased cytoplasmic Ca2+ concentration, and then blocked formation of multi-nuclei myotubes and expression of MRFs. As a major Ca2+ release channel of ER, the dramatic increase of RyR1 expression should be responsible for the significant elevation of cytosolic Ca2+ concentration. Therefore, we deduced that dramatic increase of RyR1 expression is required for the myoblast fusion at later stage of myogenesis. In addition, CAV1.1 is a physiological activator of RyR1 in the excitation–contraction coupling of skeletal muscle . However, in this study, the expression of CAV1.1 was not significantly increased as RyR1 during myogenic differentiation, which indicated that RyR1 expression in myogenic cells is independent of CAV1.1.
Thapsigargin, an inhibitor of SERCA, which mediated the reuptake of cytoplasmic Ca2+ into the sarcoplasmic reticulum, was used to increase cytoplasmic Ca2+ concentration. In the current study, the dose of THA was far below that used in previous studies [50, 51] in order to keep cell viability from the adverse effect of THA, and guaranteed that THA did not affect the efficiency of siRyR1-knockdown. As a result, THA treatment significantly increased cytoplasmic Ca2+ concentration but did not affect expression of myogenic specific genes. However, myotube formation was significantly inhibited by siRyR1-knockdown, but recovered by THA treatment, accompanied by the expression of MRFs. Therefore, RyR1-mediated elevation of cytoplasmic Ca2+ concentration was indispensable for myotubes formation.
To further explore the roles of RyR1 in myogenic differentiation and subsequent myoblast fusion, a model of RyR1-KO myoblasts by CRISPR/Cas9 gene-editing was employed in this study. The RyR1 Ca2+ release channel is composed of macromolecular complexes consisting of a homotetramer of 560-kDa RyR1 subunits that form scaffolds for proteins that regulate channel function including protein kinase A (PKA) and the phosphodiesterase PDE4D3 (both of which are targeted to the channel via the anchoring protein mAKAP), PP1 (targeted via spinophilin), and calstabin1 (FKBP12) [52, 53]. In contrast to the result of siRNA interference, we observed higher concentration of Ca2+ in the cytoplasm of RyR1-KO myoblasts, which suggested that considering the role of RyR1 as a calcium channel, RyR1-KO could differ from siRyR1 knockdown in terms of controlling Ca2+ dynamics between ER and cytoplasm. The deletion of RyR1 resulted in the dysfunction of the protein complex, losing the control of Ca2+ flowing out of the ER and accelerated Ca2+ flooding into the cytoplasm from the ER, which also distorted the expression patterns of other Ca2+ transporters. This observation is consistent with previous studies that demonstrated both over-activation and mutations of RyR1 resulted in leakage of Ca2+ from the SR [13, 54]. Cytoplasmic Ca2+ elevation directly activated ER stress [22, 23], and has a strong influence on differentiation through oxidative signaling and G0/G1 cell cycle arrest . In the current study, expression of Cyclin D1 was abolished and cell viability was decreased in RyR1-KO cells, demonstrating that cell cycle and cell viability were suppressed during the enhanced ER stress caused by RyR1-KO mediated cytoplasmic Ca2+ elevation.
ER stress signaling gives rise to apoptosis . In the present study, myogenic differentiation potential of RyR1-KO myoblasts reflected by MRFs expression was dramatically enhanced. However, the differentiation process of RyR1-KO myoblasts was interrupted by apoptosis, which indicated that RyR1 knockout make myoblasts too fragile to undertake the stress of myogenic induction. UPR was a protective response for cells under stress, however, excessive or prolonged UPR can cause apoptosis . Caspase-3, belonging to a highly conserved family of cysteinyl aspartate-specific proteases, is an essential regulator of apoptosis. In the current study, cleaved activation of caspase-3 was significantly stimulated in RyR1-KO cells, which was supported by a previous study of RyRs-mediated ER stress . In the present study, the total and phosphorylated protein expressions of IRE1α and PERK were sharply increased in RyR1-KO cells, indicating that RyR1-knockout activated serious ER stress. Furthermore, aggravated ER stress significantly increased expression of CHOP, caspase-9, and caspase-12. Meanwhile, the positive effects of ERP44 and HSPA5 against ER stress were not enhanced in RyR1-KO myoblasts. In addition, phosphorylation levels of JNK and Erk1/2, whose activations are beneficial for the resistance to ER stress-induced apoptosis [59, 60], were also significantly reduced in RyR1-KO cells. Therefore, we deduced that RyR1 deletion led to serious ER stress and excessive UPR than the tolerance thresholds of cells, and accounted for the apoptosis of RyR1-KO myoblasts upon myogenic induction.
In summary, dramatic increase in RyR1 expression is indispensable for myogenesis, and RyR1-mediated Ca2+ release plays a critical role in myoblast fusion at the later stage of myogenesis. RyR1-KO led to cytoplasmic Ca2+ elevation and enhanced myogenic differentiation potential, while serious ER stress and excessive UPR over the tolerance thresholds of cells resulted by RyR1-KO triggered the process of apoptosis of myoblasts upon myogenic induction. This study contributes to a novel understanding of the role of RyR1 in muscle development and related congenital myopathies, and provides a potential target for regulation of muscle characteristics and meat quality.
Availability of data and materials
All data generated or analyzed during this study are included in this published article [and its supplementary information files].
ATPase ER/SR Ca2+ transporting 2
ATPase plasma membrane Ca2+ transporting 1
CRISPR associated protein 9
Calcium voltage-gated channel subunit alpha1 S
C/EBP homologous protein
Calcium release activated channel regulator 2A
Calcium release activated channel regulator 2B
Clustered Regularly Interspaced Short Palindromic Repeats
DNA-damage inducible transcript 3
Eukaryotic translation initiation factor 2α
Extracellular regulated protein kinases
ER protein 44
Heat shock protein 5
Inositol 1,4,5-trisphosphate receptor type 1
Stress-activated protein kinase/Jun-amino-terminal kinase
Myogenic regulatory transcription factors
Myogenic factor 5
Myogenic differentiation 1
ORAI calcium release-activated calcium modulator 2
Protein kinase A
Ryanodine receptor 1
Store-operated calcium entry
Stromal interaction molecule 1
Tris-Buffered saline and Tween-20
Unfolded protein response
Maclennan DH, Duff C, Zorzato F, Fujii J, Phillips M, Korneluk RG, et al. Ryanodine receptor gene is a candidate for predisposition to malignant hyperthermia. Nature. 1990;343(6258):559–61. https://doi.org/10.1038/343559a0.
Rachel R, Danielle C, Marie-Anne S, Jane H, Philip H. Mutations in RYR1 in malignant hyperthermia and central core disease. Hum Mutat. 2010;27(10):977–89. https://doi.org/10.1002/humu.20356.
Ravenscroft G, Davis MR, Lamont P, Forrest A, Laing NG. New era in genetics of early-onset muscle disease: breakthroughs and challenges. Semin Cell Dev Biol. 2017;64:160–70. https://doi.org/10.1016/j.semcdb.2016.08.002.
Rokach O, Sekulic-Jablanovic M, Voermans N, Wilmshurst J, Pillay K, Heytens L, et al. Epigenetic changes as a common trigger of muscle weakness in congenital myopathies. Hum Mol Genet. 2015;24(16):4636–47. https://doi.org/10.1093/hmg/ddv195.
Efremov RG, Leitner A, Aebersold R, Raunser S. Architecture and conformational switch mechanism of the ryanodine receptor. Nature. 2015;517(7532):39–43. https://doi.org/10.1038/nature13916.
Ran Z, Clarke OB, Amédée DG, Grassucci RA, Steven R, Filippo M, et al. Structure of a mammalian ryanodine receptor. Nature. 2015;517(7532):44–9. https://doi.org/10.1038/nature13950.
Yan Z, Bai XC, Yan CY, Wu JP, Li ZQ, Xie T, et al. Structure of the rabbit ryanodine receptor RyR1 at near-atomic resolution. Nature. 2015;517(7532):50–5. https://doi.org/10.1038/nature14063.
Woll KA, Haji-Ghassemi O, Van Petegem F. Pathological conformations of disease mutant ryanodine receptors revealed by cryo-EM. Nat Commun. 2021;12(1):807. https://doi.org/10.1038/s41467-021-21141-3.
Iyer KA, Hu Y, Nayak AR, Kurebayashi N, Murayama T, Samsó M. Structural mechanism of two gain-of-function cardiac and skeletal RyR mutations at an equivalent site by cryo-EM. Sci Adv. 2020;6(31):eabb2964. https://doi.org/10.1126/sciadv.abb2964.
Brennan S, Garcia-Castañeda M, Michelucci A, Sabha N, Malik S, Groom L, et al. Mouse model of severe recessive RYR1-related myopathy. Hum Mol Genet. 2019;28(18):3024–36. https://doi.org/10.1093/hmg/ddz105.
Capogrosso RF, Mantuano P, Uaesoontrachoon K, Cozzoli A, Giustino A, Dow T, et al. Ryanodine channel complex stabilizer compound S48168/ARM210 as a disease modifier in dystrophin-deficient mdx mice: proof-of-concept study and independent validation of efficacy. FASEB J. 2018;32(2):1025–43. https://doi.org/10.1096/fj.201700182RRR.
Lawal TA, Todd JJ, Meilleur KG. Ryanodine receptor 1-related myopathies: diagnostic and therapeutic approaches. Neurotherapeutics. 2018;15(4):885–99. https://doi.org/10.1007/s13311-018-00677-1.
Mori S, Iinuma H, Manaka N, Ishigami-Yuasa M, Murayama T, Nishijima Y, et al. Structural development of a type-1 ryanodine receptor (RyR1) ca-release channel inhibitor guided by endoplasmic reticulum ca assay. Eur J Med Chem. 2019;179:837–48. https://doi.org/10.1016/j.ejmech.2019.06.076.
Suman M, Sharpe JA, Bentham RB, Kotiadis VN, Menegollo M, Pignataro V, et al. Inositol trisphosphate receptor-mediated Ca2+ signalling stimulates mitochondrial function and gene expression in core myopathy patients. Hum Mol Genet. 2018;27(13):2367–82. https://doi.org/10.1093/hmg/ddy149.
Kushnir A, Todd JJ, Witherspoon JW, Yuan Q, Reiken S, Lin H, et al. Intracellular calcium leak as a therapeutic target for RYR1-related myopathies. Acta Neuropathol. 2020;139(6):1089–104. https://doi.org/10.1007/s00401-020-02150-w.
Nasipak BT, Padilla-Benavides T, Green KM, Leszyk JD, Mao W, Konda S, et al. Opposing calcium-dependent signalling pathways control skeletal muscle differentiation by regulating a chromatin remodelling enzyme. Nat Commun. 2015;6(1):7441. https://doi.org/10.1038/ncomms8441.
Bellinger AM, Reiken S, Dura M, Murphy PW, Deng SX, Landry DW, et al. Remodeling of ryanodine receptor complex causes "leaky" channels: a molecular mechanism for decreased exercise capacity. Proc Natl Acad Sci U S A. 2008;105(6):2198–202. https://doi.org/10.1073/pnas.0711074105.
Fujii J, Otsu K, Zorzato F, de Leon S, Khanna VK, Weiler JE, et al. Identification of a mutation in porcine ryanodine receptor associated with malignant hyperthermia. Science. 1991;253(5018):448–51. https://doi.org/10.1126/science.1862346.
Pisaniello A, Serra C, Rossi D, Vivarelli E, Sorrentino V, Molinaro M, et al. The block of ryanodine receptors selectively inhibits fetal myoblast differentiation. J Cell Sci. 2003;116(Pt 8):1589–97. https://doi.org/10.1242/jcs.00358.
Filipova D, Walter AM, Gaspar JA, Brunn A, Linde NF, Ardestani MA, et al. Gene profiling of embryonic skeletal muscle lacking type I ryanodine receptor ca (2+) release channel. Sci Rep. 2016;6(1):20050. https://doi.org/10.1038/srep20050.
Filipova D, Henry M, Rotshteyn T, Brunn A, Carstov M, Deckert M, et al. Distinct transcriptomic changes in E14.5 mouse skeletal muscle lacking RYR1 or Cav1.1 converge at E18.5. PLoS ONE. 2018;13(3):e0194428. https://doi.org/10.1371/journal.pone.0194428.
Ozcan L, Tabas I. Role of endoplasmic reticulum stress in metabolic disease and other disorders. Annu Rev Med. 2012;63(63):317–28. https://doi.org/10.1146/annurev-med-043010-144749.
Ron D, Walter P. Signal integration in the endoplasmic reticulum unfolded protein response. Nat Rev Mol Cell Biol. 2007;8(7):519–29. https://doi.org/10.1038/nrm2199.
Jheng JR, Chen YS, Ao UI, Chan DC, Huang JW, Hung KY, et al. The double-edged sword of endoplasmic reticulum stress in uremic sarcopenia through myogenesis perturbation. J Cachexia Sarcopenia Muscle. 2018;9(3):570–84. https://doi.org/10.1002/jcsm.12288.
Ozcan U, Cao Q, Yilmaz E, Lee AH, Iwakoshi NN, Ozdelen E, et al. Endoplasmic reticulum stress links obesity, insulin action, and type 2 diabetes. Science. 2004;306(5695):457–61. https://doi.org/10.1126/science.1103160.
Villalobos-Labra R, Subiabre M, Toledo F, Pardo F, Sobrevia L. Endoplasmic reticulum stress and development of insulin resistance in adipose, skeletal, liver, and foetoplacental tissue in diabesity. Mol Asp Med. 2019;66:49–61. https://doi.org/10.1016/j.mam.2018.11.001.
Bohnert KR, McMillan JD, Kumar A. Emerging roles of ER stress and unfolded protein response pathways in skeletal muscle health and disease. J Cell Physiol. 2018;233(1):67–78. https://doi.org/10.1002/jcp.25852.
Günther S, Kim J, Kostin S, Lepper C, Fan CM, Braun T. Myf5-positive satellite cells contribute to Pax7-dependent long-term maintenance of adult muscle stem cells. Cell Stem Cell. 2013;13(5):590–601. https://doi.org/10.1016/j.stem.2013.07.016.
Qiu K, Xu D, Wang L, Zhang X, Jiao N, Gong L, et al. Association analysis of single-cell RNA sequencing and proteomics reveals a vital role of ca signaling in the determination of skeletal muscle development potential. Cells. 2020;9(4):E1045. https://doi.org/10.3390/cells9041045.
Sun W, He T, Qin C, Qiu K, Zhang X, Luo Y, et al. A potential regulatory network underlying distinct fate commitment of myogenic and adipogenic cells in skeletal muscle. Sci Rep. 2017;7(1):44133. https://doi.org/10.1038/srep44133.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) method. Methods. 2001;25(4):402–8. https://doi.org/10.1006/meth.2001.1262.
Ikemoto T, Hosoya T, Aoyama H, Kihara Y, Suzuki M, Endo M. Effects of dantrolene and its derivatives on Ca2+ release from the sarcoplasmic reticulum of mouse skeletal muscle fibres. Br J Pharmacol. 2010;134(4):729–36. https://doi.org/10.1038/sj.bjp.0704307.
Comai G, Tajbakhsh S. Molecular and cellular regulation of skeletal myogenesis. Curr Top Dev Biol. 2014;110:1–73. https://doi.org/10.1016/B978-0-12-405943-6.00001-4.
Bijlenga P, Liu JH, Espinos E, Haenggeli CA, Fischer-Lougheed J, Bader CR, et al. T-type α 1H Ca2+ channels are involved in Ca2+ signaling during terminal differentiation (fusion) of human myoblasts. Proc Natl Acad Sci U S A. 2000;97(13):7627–32. https://doi.org/10.1073/pnas.97.13.7627.
Porter GA Jr, Makuck RF, Rivkees SA. Reduction in intracellular calcium levels inhibits myoblast differentiation. J Biol Chem. 2002;277(32):28942–7. https://doi.org/10.1074/jbc.M203961200.
Qiu K, Zhang X, Wang L, Jiao N, Xu D, Yin J. Protein expression landscape defines the differentiation potential specificity of adipogenic and myogenic precursors in the skeletal muscle. J Proteome Res. 2018;17(11):3853–65. https://doi.org/10.1021/acs.jproteome.8b00530.
Airey JA, Baring MD, Sutko JL. Ryanodine receptor protein is expressed during differentiation in the muscle cell lines BC3H1 and C2C12. Dev Biol. 1991;148(1):365–74. https://doi.org/10.1016/0012-1606(91)90344-3.
Raffaello A, Mammucari C, Gherardi G, Rizzuto R. Calcium at the center of cell signaling: interplay between endoplasmic reticulum, mitochondria, and lysosomes. Trends Biochem Sci. 2016;41(12):1035–49. https://doi.org/10.1016/j.tibs.2016.09.001.
Rink TJ. Receptor-mediated calcium entry. FEBS Lett. 1990;268(2):381–5. https://doi.org/10.1016/0014-5793(90)81290-5.
Mackrill JJ. Protein–protein interactions in intracellular Ca2+−release channel function. Biochem J. 1999;337(Pt 3, 3):345–61.
Villalobos C, García-sancho J. Capacitative Ca2+ entry contributes to the Ca2+ influx induced by thyrotropin-releasing hormone (TRH) in GH3 pituitary cells. Pflügers Arch. 1995;430(6):923–35. https://doi.org/10.1007/BF01837406.
Misquitta CM, Mack DP, Grover AK. Sarco/endoplasmic reticulum ca 2+ (SERCA)-pumps: link to heart beats and calcium waves. Cell Calcium. 1999;25(4):277–90. https://doi.org/10.1054/ceca.1999.0032.
Bruce JIE. Metabolic regulation of the PMCA: role in cell death and survival. Cell Calcium. 2018;69:28–36. https://doi.org/10.1016/j.ceca.2017.06.001.
Brini M, Carafoli E. The plasma membrane Ca2+ ATPase and the plasma membrane sodium calcium exchanger cooperate in the regulation of cell calcium. Cold Spring Harb Perspect Biol. 2011;3(2):487–96. https://doi.org/10.1101/cshperspect.a004168.
Shaw PJ, Qu B, Hoth M, Feske S. Molecular regulation of CRAC channels and their role in lymphocyte function. Cell Mol Life Sci. 2013;70(15):2637–56. https://doi.org/10.1007/s00018-012-1175-2.
Prakriya M, Lewis RS. Store-operated calcium channels. Physiol Rev. 2015;95(4):1383–436. https://doi.org/10.1152/physrev.00020.2014.
Tapscott SJ. The circuitry of a master switch: Myod and the regulation of skeletal muscle gene transcription. Development. 2005;132(12):2685–95. https://doi.org/10.1242/dev.01874.
Saab R, Spunt SL, Skapek SX. Myogenesis and rhabdomyosarcoma : the Jekyll and Hyde of skeletal muscle. Curr Top Dev Biol. 2011;94(94):197–234. https://doi.org/10.1016/B978-0-12-380916-2.00007-3.
Schartner V, Romero NB, Donkervoort S, Treves S, Munot P, Pierson TM, et al. Dihydropyridine receptor (DHPR, CACNA1S) congenital myopathy. Acta Neuropathol. 2017;133(4):517–33. https://doi.org/10.1007/s00401-016-1656-8.
Földi I, Tóth AM, Szabó Z, Mózes E, Berkecz R, Datki ZL, et al. Proteome-wide study of endoplasmic reticulum stress induced by thapsigargin in N2a neuroblastoma cells. Neurochem Int. 2013;62(1):58–69. https://doi.org/10.1016/j.neuint.2012.11.003.
Zhang X, Yuan Y, Jiang L, Zhang J, Gao J, Shen Z, et al. Endoplasmic reticulum stress induced by tunicamycin and thapsigargin protects against transient ischemic brain injury: involvement of PARK2-dependent mitophagy. Autophagy. 2014;10(10):1801–13. https://doi.org/10.4161/auto.32136.
Brillantes AB, Ondrias K, Scott A, Kobrinsky E, Ondriasová E, Moschella MC, et al. Stabilization of calcium release channel (ryanodine receptor) function by FK506-binding protein. Cell. 1994;77(4):513–23. https://doi.org/10.1016/0092-8674(94)90214-3.
Marx SO, Reiken S, Hisamatsu Y, Gaburjakova M, Gaburjakova J, Yang YM, et al. Phosphorylation-dependent regulation of ryanodine receptors: a novel role for leucine/isoleucine zippers. J Cell Biol. 2001;153(4):699–708. https://doi.org/10.1083/jcb.153.4.699.
Guerrero-Hernández A, Ávila G, Rueda A. Ryanodine receptors as leak channels. Eur J Pharmacol. 2014;739:26–38. https://doi.org/10.1016/j.ejphar.2013.11.016.
Thrivikraman G, Madras G, Basu B. Electrically driven intracellular and extracellular nanomanipulators evoke neurogenic/cardiomyogenic differentiation in human mesenchymal stem cells. Biomaterials. 2016;77:26–43. https://doi.org/10.1016/j.biomaterials.2015.10.078.
Nakanishi K, Dohmae N, Morishima N. Endoplasmic reticulum stress increases myofiber formation in vitro. FASEB J. 2007;21(11):2994–3003. https://doi.org/10.1096/fj.06-6408com.
Feng R, Zhai WL, Yang HY, Jin H, Zhang QX. Induction of ER stress protects gastric cancer cells against apoptosis induced by cisplatin and doxorubicin through activation of p38 MAPK. Biochem Biophys Res Commun. 2011;406(2):299–304. https://doi.org/10.1016/j.bbrc.2011.02.036.
Wang H, Dong Y, Zhang J, Xu Z, Wang G, Swain CA, et al. Isoflurane induces endoplasmic reticulum stress and caspase activation through ryanodine receptors. Br J Anaesth. 2014;113(4):695–707. https://doi.org/10.1093/bja/aeu053.
Urano F, Wang X, Bertolotti A, Zhang Y, Chung P, Harding HP, et al. Coupling of stress in the ER to activation of JNK protein kinases by transmembrane protein kinase IRE1. Science. 2000;287(5453):664–6. https://doi.org/10.1126/science.287.5453.664.
Chen Q, Hang Y, Zhang T, Tan L, Li S, Jin Y. USP10 promotes proliferation and migration and inhibits apoptosis of endometrial stromal cells in endometriosis through activating the Raf-1/MEK/ERK pathway. Am J Physiol Cell Physiol. 2018;315(6):C863–C72. https://doi.org/10.1152/ajpcell.00272.2018.
This study was financially supported by the National Natural Science Foundation of China (Grant No. 31790412), National key research and development program of China (Grant No. 2018YFD0500402), and the National Natural Science Foundation of China (Grant No. 31672431).
Ethics approval and consent to participate
All studies using pigs were conducted in accordance with the guidelines of, and approved by, the Institutional Animal Care and Use Committee of China Agricultural University.
Consent for publication
The authors declare no conflict of interest.
Cytoplasmic Ca2+ concentration of C2C12 cells during myogenic differentiation. (A) Immunostaining with myosin antibody on d 0, 2, 4, and 6 during myogenic differentiation. (B) Cytoplasmic Ca2+ signals of unfused C2C12 cells labeled by Fluo-3 on d 0–5 during myogenic differentiation (n = 3). (C) Quantitative results of cytoplasmic Ca2+ signals. (D-E) The proteins expression of RyR1 in C2C12 cells during myogenic differentiation (n = 3). GM: growth medium, representing cells cultured in growth medium before myogenic induction; DM: differentiation medium, representing cells on d 4 during myogenic differentiation cultured in differentiation medium. The data are presented as the mean ± SEM. *Represents significant difference with the value on d 0 (P < 0.05).
The mRNA expression of RyR3 during myogenic differentiation of C2C12 cells. *Represents significant difference with the value at D0 (P < 0.05). (B) The ratio of RyR3 and RyR1 mRNA expression in proliferating C2C12 cells. (C) The mRNA expression of RyR1 and RyR3 during myogenic differentiation of myogenic cells of pigs. (D) The ratio of RyR3 and RyR1 mRNA expression in proliferating myogenic cells of pigs. (E) The ratio of RyR3 and RyR1 mRNA expression of myogenic cells of pigs on d 4 during myogenic differentiation. *Represents significant difference between the two groups (P < 0.05). GM: growth medium, representing cells cultured in growth medium before myogenic induction; DM: differentiation medium, representing cells on d 4 during myogenic differentiation cultured in differentiation medium.
The establishment of RyR1-knockout cells. (A) The plasmid vector used for CRISPR/Cas9 gene-editing. (B) The shearing efficiency of three gRNA targeted to RyR1. (C) The target site of the selected gRNA on RyR1 gene. (D) Homozygote and heterozygote of RyR1-knockout, named as RyR1−/− and RyR1+/−, was verified by gene sequencing.
The primers used for qRT-PCR assays. Table S2 Information of antibodies used for Western Blot. Table S3 The potential off-target sites (OTS) of CRISPR/Cas9 system. Table S4 Primers used for PCR and DNA sequencing.
About this article
Cite this article
Qiu, K., Wang, Y., Xu, D. et al. Ryanodine receptor RyR1-mediated elevation of Ca2+ concentration is required for the late stage of myogenic differentiation and fusion. J Animal Sci Biotechnol 13, 9 (2022). https://doi.org/10.1186/s40104-021-00668-x