Myogenic Transdifferentiation of Pig Fibroblasts Using the Porcine MYOD1 Inducible Vector

doi:10.21203/rs.3.rs-102478/v1

Download PDF

Research

Myogenic Transdifferentiation of Pig Fibroblasts Using the Porcine MYOD1 Inducible Vector

https://doi.org/10.21203/rs.3.rs-102478/v1

This work is licensed under a CC BY 4.0 License

Journal Publication

published 13 May, 2021

Read the published version in Journal of Animal Science and Biotechnology →

You are reading this latest preprint version

Background Myogenic transdifferentiation can be accomplished through ectopic MYOD1 expression, which is facilitated by various signaling pathways associated with myogenesis. In this study, we attempted to transdifferentiate pig embryonic fibroblasts (PEFs) myogenically into skeletal muscle through overexpression of the pig MYOD1 gene and modulation of the FGF, TGF-β, WNT, and cAMP signaling pathways.

Results The MYOD1 overexpression vector was constructed based on comparative sequence analysis, demonstrating that pig MYOD1 has evolutionarily conserved domains across various species. Although forced MYOD1 expression through these vectors triggered the expression of endogenous muscle markers, transdifferentiated muscle cells from fibroblasts were not observed. Therefore, various signaling molecules, including FGF2, SB431542, CHIR99021, and forskolin, along with MYOD1 overexpression were applied to enhance the myogenic reprogramming. The modified conditions led to the derivation of myotubes and activation of muscle markers in PEFs, as determined by qPCR and immunostaining. Notably, a sarcomere-like structure was observed, indicating that terminally differentiated skeletal muscle could be obtained from transdifferentiated cells.

Conclusions In summary, we established a protocol for reprogramming MYOD1-overexpressing PEFs into the mature skeletal muscle using signaling molecules. Our myogenic reprogramming can be used as a cell source for muscle disease models in regenerative medicine and the production of cultured meat in cellular agriculture.

Animal Science

Biotechnology and Bioengineering

MYOD1

pig

transdifferentiation

skeletal muscle

sequence analysis

To date, numerous studies have defined various differentiation strategies for the differentiation of PSCs (pluripotent stem cells) into specific cell types of ecto-, meso-, and endodermal lineages, which have advantages as infinite cell sources (1–3). Whereas these directed differentiation approaches recapitulate in vivo developmental processes through modulation of signaling pathways, they reportedly have defects in terms of time consumption and low efficiency because the small molecules used are also involved in differentiation into other cell types (4, 5). For example, the derivation of striated myofibers, a form of mature skeletal muscle, from mouse ESCs (embryonic stem cells) and human iPSCs (induced pluripotent stem cells) reportedly takes 3 ~ 4 weeks, albeit with the use of an advanced protocol with a short experimental period (6, 7). Therefore, direct reprogramming using the activation of the transcriptional program has been applied to facilitate the differentiation of PSCs into the desired cell type (5). However, problems derived from the use of PSCs still remain, including the potential risk of teratoma formation (8) and the requirement of a time-consuming and fine-tuning process for the derivation and maintenance of the pluripotent state. To overcome the above problems, a transdifferentiation approach is required, in which fully differentiated somatic cells, instead of PSCs, are induced to differentiate directly into target cell types via ectopic expression of transcription factors. Transdifferentiation toward the muscle lineage has been widely studied for a long time since the identification of the master transcription factor in myogenesis, Myod1 (9).

Forced Myod1 expression converts cell fates into muscle, which is responsible for Myod1 functioning as a transcriptional and epigenetic regulator, leading to the activation of muscle-specific genes in a feed-forward manner (10, 11). Small molecules can be supplemented under culture conditions to Myod1/MYOD1-overexpressing somatic cells to improve the efficacy of transdifferentiation. The surface of myocytes and the basement membrane of myotubes secrete FGF2 (fibroblast growth factor 2), which play roles in the activation of muscle stem cells and the expansion of myoblasts, inhibiting myogenic terminal differentiation (12, 13). Because TGF-β signaling negatively functions in myogenesis by decreasing Myog activity, inhibition of TGF type I-β receptors, such as SB431542, reportedly enhances the elongation of myotubes through myoblast fusion (13). In fact, the expression of follistatin, which inhibits TGF-β signaling, is detected in the paraxial mesoderm in vivo (14). During in vivo myogenesis, the WNT activator secreted from the notochord, neural tube, and surrounding tissue is involved in a series of specifications along with the presomitic mesoderm, somite, dermomyotome, and myotome, as determined by in vitro-directed differentiation in human and mouse PSCs (3, 6, 7). It demonstrates that WNT signaling induces the commitment of myogenic precursors and has been substantiated by previous research showing that the WNT activator CHIR99021 with SB431542 and FGF2 induces myogenic specification from human iPSCs (14). In a previous report, the adenylyl cyclase activator forskolin was identified as a myoblast proliferation-promoting factor through CREB-mediated WNT, leading to the upregulation of the Pax3, Myf5, and Myod1 genes (15). The combination of FGF2, a WNT activator, and forskolin stimulated skeletal muscle differentiation in human iPSCs and especially forskolin improved satellite cell expansion in mice (16).

In pigs, myogenic reprogramming could be used as a cell source for muscle disease models and the production of cultured meat. Although, in pig iPSCs, skeletal myotubes were reportedly observed through the activation of WNT signaling and ectopically expressed MYOD1 by overexpression vector and 5-azacytidine (5-aza) (17), the aforementioned limitations derived from iPSCs remain. Moreover, it has been demonstrated that muscle stem cells lose the potential for proliferation and differentiation in long-term culture in vitro (18). As an alternative, fibroblasts can be obtained in high yield because of their large proportion in the body and are known to be the most effective cell type to accept expression of Myod1 because they are derived from the mesoderm, which is identical in origin to the muscle (19). Therefore, in this study, we attempted to transdifferentiate pig fibroblasts into skeletal muscle myogenically through overexpression of the pig MYOD1 gene and modulation of the FGF, TGF-β, WNT, and cAMP signaling pathways.

Amino acid sequence analysis

The MYOD1 amino acid sequences of pig (NP_001002824.1), mouse (NP_034996.2), human (NP_002469.2), horse (NP_001304182.1), cow (NP_001035568.2), and sheep (NP_001009390.1) were obtained from the NCBI database (https://www.ncbi.nlm.nih.gov/). Comparative analysis of the amino acid sequences was performed using the multiple sequence alignment program MUSCLE (MUltipleby https://www.ebi.ac.uk/Tools/msa/muscle/), and the similarity of proteins from pigs and the other species was analyzed using Protein BLAST (https://blast.ncbi.nlm.nih.gov/).

Construction of inducible p MYOD1 -overexpression vector

Total RNA was extracted from satellite cells from the biceps femoris of 3-day-old LYD pigs based on our previous study (20), and cDNA was subsequently synthesized. Nested PCR was performed using an outer primer containing 20-bp sequences upstream and downstream of the MYOD1 gene and an inner primer containing the EcoR1 sequence. The following primer set was used: Outer-MYOD1-F, 5'-ATAGAGCAGGGTGGTGGACA-3‘, Outer-MYOD1-R, 5'-CTCAAACTTCTGGGCGCGAG-3‘, Inner-MYOD1-F, 5'-GAATTCTGGGATATGGAGCTGCTGTCGC-3‘, Inner-MYOD1-R, 5'-GAATTCTCAGAGCACCTGGTAGATAGGGGTTGG-3‘. The PCR product was purified using electrophoresis and inserted into a T easy vector (Promega, Madison, WI, USA). DNA was extracted from the selected vectors through TA cloning and was sequenced using the following M-13 primer: F, 5'-GTAAAACGACGGCCAG-3‘, R, 5'-CAGGAAACAGCTATGAC-3‘. Subsequently, a doxycycline (DOX)-inducible FUW-tetO-pMYOD1 vector was constructed using treatment with the restriction enzyme EcoR1 and the FUW-tetO-MCS vector (plasmid #84008; Addgene, Watertown, MA, USA).

Lentiviral vector production

Lentiviral vectors were produced as previously described (21). Briefly, HEK 293 LTV cells (Cell Biolabs, San Diego, CA, USA) were used as the packaging cell line, and five plasmids were used for the production of lentiviral vectors: FUW-tetO-pMYOD1 and FUW-M2rtTA (the transfer plasmid); pLP1 and pLP2 (the packaging plasmids; Invitrogen, Waltham, MA, USA); and pLP/VSVG (the envelope plasmid; Invitrogen). These plasmids were transfected into HEK 293 LTV cells using the calcium phosphate precipitation method. Subsequently, the LTV culture supernatants were filtered and concentrated. The derived virus pellets were stored at -76 °C until use.

Myogenic transdifferentiation of pig fibroblasts

The pig embryonic fibroblasts (PEFs) isolated in previous study were used (22). The cells were cultured in Dulbecco’s modified Eagle’s medium (DMEM; Welgene, Gyeongsan, Korea) supplemented with 10% fetal bovine serum (FBS; collected and processed in the United States; Genedepot, Katy, TX, USA), 1 × Glutamax (Gibco, Gaithersburg, MD, USA), 0.1 mM ß-mercaptoethanol (Gibco), and 1 × antibiotic/antimycotic (Gibco). The exogenous gene was transduced with a lentivirus carrying the FUW-tetO-pMYOD1 or FUW-M2rtTA plasmid. The cells were incubated with lentivirus using 8 mg/ml polybrene (Sigma-Aldrich, St. Louis, MO, USA) for 48 hours.

The virus-infected cells were cultured under various conditions, which consisted of basal medium and signaling molecules. The basal medium was composed of DMEM, 10% FBS, 10% knockout serum replacement (KSR; Gibco), 10% 1 × Glutamax, 0.1 mM ß-mercaptoethanol, and 1 × antibiotic/antimycotic. The tested signaling molecules were 20 ng/ml fibroblast growth factor 2 (FGF2; R&D Systems, Minneapolis, MN, USA), 4 µM TGF β receptor inhibitor SB431542 (Cayman chemical, Ann Arbor, MI, USA), 3 µM GSK3B inhibitor CHIR99021 (Cayman chemical), and 10 µM cAMP activator forskolin (Cayman chemical). DOX (4 ng/ml) was added to activate the inducible vector. In the ‘differentiation step’, each medium was replaced with a serum-free basal medium containing 2% horse serum. The media were changed every day, and all cells were cultured under humidified conditions with 5% CO₂ at 37 °C.

Genomic DNA (gDNA) extraction and polymerase chain reaction (PCR)

Genomic DNA was extracted using the G-spin™ Total DNA Extraction Kit (iNtRON, Seongnam, Korea). Amplifications were performed using the transgene-specific primers listed in Tables 1 and 2 × PCR Master mix solution (iNtRON) containing 5 pmol of each primer set and 50 ng gDNA in a 10 µl reaction volume. PCRs were performed in a thermocycler under the following conditions: 94 °C for 5 min, followed by 35 cycles of denaturation at 95 °C for 30 sec, annealing for 30 sec (annealing temperatures depended on each primer set), and extension at 72 °C for 30 sec, with a final extension at 72 °C for 7 min. The amplified PCR products were visualized using electrophoresis on a 1% agarose gel stained with ethidium bromide.

Table 1

Primer sets used for PCR and qPCR.
Gene	Primer sequence		Annealing Temperature (°C)	Product size (bp)
FUW-TetO-pMYOD1 (ExoMYOD1)	F	CCAGGTGCTCTGAGAATTCGATA	60	114
	R	CCACATAGCGTAAAAGGAGCA
EndoMYOD1	F	AGGGACAGGATAGAGCAGGG	60	199
	R	TCAAATCTACGTCGCGGAGC
PAX7	F	GTGCCCTCAGTGAGTTCGAT	58	152
	R	TCCAGACGGTTCCCTTTGTC
MYF5	F	AGTTCGGGGACGAGTTTGAG	60	232
	R	TCAAACGCCTGGTTGACCTT
MYOG	F	GAGCTGTATGAGACATCCCCC	60	75
	R	GTGGACGGGCAGGTAGTTTT
MHC	F	ACAGTGAAGACGGAAGCAGG	60	153
	R	TGCGTAACGCTCTTTGAGGT
pACTB	F	CCGGGACCTGACCGACTACC	60	126
	R	TCGAAGTCCAGGGCGACGTA
GAPDH	F	TGCTCCTCCCCGTTCGAC	60	100
	R	ATGCGGCCAAATCCGTTC

Table 2

Similarity analysis of the pig MYOD1 amino acid sequence compared to other species.
	Similarity between pig and the other species (%)
	Mouse	Human	Horse	Cow	Sheep
MYOD1 protein	89.03	93.12	94.36	95.92	95.61
Acidic	92.59	92.59	98.15	100.00	100.00
H/C	89.19	97.30	89.19	100.00	100.00
bHLH	100.00	100.00	100.00	100.00	100.00
Helix3	92.86	92.86	92.86	92.86	92.86
Similarity comparison between pig and the other species was confirmed on the score via the Protein BLAST program.

Immunocytochemistry (ICC) analysis

Before staining, all cell samples were preincubated for 5 min at 4 °C and fixed with 4% paraformaldehyde for 15 min. After washing twice with Dulbecco’s phosphate-buffered saline (DPBS; Welgene), the samples were treated for 1 h with 10% goat serum in DPBS to prevent nonspecific binding. Serum-treated cells were incubated overnight at 4 °C with primary antibodies. The primary antibodies used were as follows: rabbit anti-MYOD1 (1:100, Thermo Fisher, Waltham, MA, USA; PA5-23078) and mouse anti-MHC (1:50, Sigma-Aldrich; 05-716). When antibodies against MYOD1 were applied, fixed cells were treated for 15 min with 0.2% Triton-X100 (Sigma-Aldrich) before serum blocking. After incubation with the primary antibody, the cells were treated for 2 h at room temperature with an Alexa Fluor-conjugated secondary antibody. Nuclei were stained with Hoechst 33342 (Molecular Probes, Eugene, OR, USA). Images of stained cells were captured using a TE2000-U inverted microscope (Nikon, Tokyo, Japan).

Quantitative real-time polymerase chain reaction (qPCR)

Total RNA was extracted from the cells using TRIzol^® Reagent (Invitrogen) according to the manufacturer’s instructions. cDNA was synthesized using the High Capacity RNA-to-cDNA Kit (Applied Biosystems, Waltham, MA, USA), producing a final volume of 20 µl. The derived cDNA samples were amplified with PowerSYBR^Ⓡ Green PCR Master Mix (Applied Biosystems) containing 0.5 pmol of each primer set listed in Table 1 in a 10 µl reaction volume. Amplification and detection were conducted using the ABI 7300 Real-Time PCR system (Applied Biosystems) under the following conditions: one cycle of 50 °C for 2 min and 95 °C for 10 min, followed by 40 or 45 cycles of denaturation at 95 °C for 15 sec and annealing/extension for 1 min (annealing/extension temperatures depended on each primer set). The relative expression level was calculated by normalizing the threshold cycle (Ct) values of each gene to that of ACTB (beta-actin) via the Δ–Ct method (23).

Statistical analysis

The data from the qPCR analyses are presented as the mean ± standard error of the mean (SEM) and were analyzed using Prism 6 software (GraphPad Software; San Diego, CA, USA). The significance of differences was determined by one-way analyses of variance followed by Fisher’s least significant difference test. Differences were considered significant at p < 0.05.

Pig MYOD1 overexpression vector construction

Myod1/MYOD1 is reportedly identified as a master transcription factor in myogenesis, thereby inducing myogenic transdifferentiation in nonmuscle cells (9, 24–27). Thus, comparative MYOD1 sequence analysis was performed among various species to assess whether pig MYOD1 contains functionally conserved sequences for myogenesis, as in the other species, before vector construction for ectopic pig MYOD1 (pMYOD1) expression. The pig MYOD1 whole protein is composed of 319 amino acids and has a similar size to MYOD1 protein in other species (Fig. 1A). Additionally, in all the analyzed species, including pig, the MYOD1 protein contains acidic domains, histidine and cysteine-rich (H/C) domains, basic-helix1-loop-helix2 (bHLH) domains, and helix3 domains, which are located at 3–56 a.a, 63–99 a.a, 102–167 a.a, and 246–259 a.a (in the case of mouse: 245–258 a.a), respectively, indicating there are identical locations and sizes among the species (28, 29). Subsequently, the MYOD1 sequences from those species were aligned, and the similarity score was assessed at the amino acid level (Fig. 1B and Table 2). In the bHLH domain, Ala¹¹⁴-Thr¹¹⁵ and Arg¹¹¹ reportedly endow the MYOD1 protein with myogenic activity (30, 31). In fact, these residues (marked in red; Fig. 1B) were confirmed to be conserved in the basic region of all the analyzed species. The MYOD1 whole proteins of pig and ungulates (pig-horse: 94.36%, pig-cow: 95.92%, and pig-sheep: 95.61%) were more similar than those of pig and human (93.12%) or rodent (89.03%). It is noteworthy that the bHLH domain showed a 100.00% similarity between pig and all other species. Altogether, porcine MYOD1, especially the bHLH domain, was identified as an evolutionarily conserved protein, which seems that its role is also conserved in myogenesis across species.

Based on the above analyses, doxycycline (DOX)-inducible pMYOD1 overexpression vectors were generated including the MYOD1 gene isolated from satellite cells in 3-day-old LYD biceps femoris. To verify the function of these constructed vectors, they were introduced into pig embryonic fibroblasts (PEFs) through lentiviral infection for stable transgene expression (4, 27). First, the integration of the exogenous pMYOD1 gene was confirmed by PCR targeting the FUW-tetO-pMYOD1 sequence in gDNA of PEFs infected with a lentivirus carrying the vectors (pMYOD1-PEFs) (Fig. 2A). At day 9, the gene of interest had been inserted stably into the genome of PEFs. Then, MYOD1 expression was confirmed at the protein level using immunostaining in pMYOD1-PEFs (Fig. 2B). The vectors were activated by the addition of DOX, leading to ectopic MYOD1 expression. Finally, the expression pattern of muscle-associated genes (Exo-MYOD1, Endo-MYOD1, PAX7, MYF5, and MYOG) was analyzed in pMYOD1-PEFs using qPCR (Fig. 2C). These genes have been characterized as myogenic lineage-specific markers, such as skeletal muscle progenitor/myoblasts (PAX7, MYF5, and MYOD1) and myocytes (MYOD1 and MYOG) (3). Exo-MYOD1 overexpression by vector activation increased the expression of endogenous muscle-associated genes. During extended cell culture, the expression of these genes was stably maintained. In conclusion, we constructed a DOX-inducible pMYOD1 overexpression vector that triggered the expression of muscle markers in long-term culture, indicating its stable function.

Myogenic transdifferentiation of PEFs by overexpression of MYOD1

The expression of endogenous muscle markers, including Endo-MYOD1, was enhanced via ectopic MYOD1 expression, as shown in Fig. 2C. However, because of the mild changes in each gene, transdifferentiated muscle cells from fibroblasts were not observed. It has been shown that full transdifferentiation is achieved by genetic modulation along with suitable culture conditions for specific cell types (32), suggesting that optimization of culture conditions is required. According to previous studies, a 2-step transdifferentiation protocol was employed with some modification to derive myogenic cells through MYOD1 overexpression (8, 17, 33). In the ‘induction’ step, the myogenic program was activated with the stimulation of transcription factors associated with myogenesis, thereby leading to the commitment into a myogenic lineage. Briefly, pMYOD1-PEFs were treated with various signaling molecules, such as FGF2, SB431542, CHIR99021, and forskolin, which have been shown to be involved in the regulation of myogenesis (8, 14). Then, transdifferentiation was promoted through serum starvation in the ‘differentiation’ step.

To ensure efficient myogenic conversion, we investigated the transition point where exogenous MYOD1 leads to the peak expression of endogenous skeletal muscle-specific genes. In addition, pMYOD1-PEFs were cultured in mitogen-rich media in which myogenic induction was sustained without entering the differentiation process (Fig. 3A). While the FSCHF group had a long cylindrical shape due to elongation of the cytoplasm, a typical fibroblastic and round shape was observed in the control and FGF groups, respectively (Fig. 3B). The alteration of the FSCHF group was maintained by the end of the culture period, resulting in a similar morphology to myoblasts, as previously reported (17, 34). The myogenic genes were upregulated by forced MYOD1 expression (Fig. 3C). Compared to other groups, the relative gene expression of the FSCHF group was higher across all the genes and culture periods. These results suggested that the Exo-MYOD1 effect was enhanced by the combination of four signaling molecules rather than DOX-induced exogenous MYOD1 per se or additional FGF2. While the gene expression of the FGF group was gradually changed, that of the FSCHF group was significantly increased on day 6 and then decreased. The aforementioned expression patterns were observed in all the markers except PAX7, which was upregulated up to day 9 in both groups. Taken together, pMYOD1-PEFs treated with a cocktail of FGF2, SB431542, CHIR99021, and forskolin for 6 days were used for further experiments.

Based on the above observations, we established a myogenic transdifferentiation protocol in which the FSCHF medium for induction into a myogenic lineage was replaced with a low-serum medium for the initiation of differentiation on day 6 (Fig. 4A). The replaced culture condition was classified into two groups distinguished by the addition of DOX (+ DOX and -DOX) to assess whether consistent activation of the myogenic program could enhance transdifferentiation. Notably, multinucleated myotubes via fusion of myoblasts were observed on day 8 in both groups (Fig. 4B). According to the qPCR analyses performed with a sample from day 9, transcripts of the Exo-MYOD1, MYF5, MYOG, and Myosin heavy chain (MHC) genes were upregulated in the +/- DOX groups (Fig. 4C). Unlike that of the other genes, gene expression of Endo-MYOD1 was significantly decreased in the +/- DOX group. Across all of the genes, especially MYF5 and MHC, the + DOX group showed significantly higher expression levels than the -DOX group, demonstrating that continuous MYOD1 overexpression during the differentiation step facilitates myogenic transdifferentiation. The expression of MHC, a marker of late differentiation in myogenesis, was detected in the day 9 sample by immunofluorescence analysis (Fig. 4D). Interestingly, a sarcomere-like structure with a striated pattern was observed, as previously reported (35), indicating that the mature myotube could be assembled. Therefore, the established protocol using ectopic MYOD1 expression and signaling molecules associated with myogenesis, such as FGF2, a TGF-β inhibitor, a WNT activator, and a cAMP activator, enabled fibroblasts to be reprogrammed into skeletal muscle.

Pig MYOD1 has an evolutionarily conserved bHLH domain that regulates myogenesis

Previously, numerous attempts have been made to reprogram nonmuscle cells into skeletal muscle by modulating the expression of the transcription factor Myod1, which plays an important role in myogenesis (10, 11). 5-azacytidine (5-aza), a DNA methyltransferase inhibitor, induces demethylation of the Myod1 locus in differentiated somatic cells and increases Myod1 transcripts, leading to transdifferentiation into the myogenic lineage (9, 36). However, 5-aza does not target site-specific demethylation; it induces genome-wide demethylation because 5-aza is incorporated into genomic DNA as a competitive analog of cytosine and then disintegrates DNA methyltransferase by trapping (37). In fact, when fibroblasts were treated with 5-aza, adipocytes and chondrocytes were also observed (38). Additionally, it has been substantiated that 5-aza exhibits cytotoxicity by causing DNA double-strand breaks and apoptosis (39). Thus, myogenic transdifferentiation using 5-aza is inappropriate in terms of safety and efficiency. As an alternative strategy, forced Myod1/MYOD1 expression without induction by 5-aza has been used to activate the myogenic program in human and mouse fibroblasts (27, 34, 40). In pig, myotubes were differentiated from induced pluripotent stem cells (iPSCs) via a combination of a 5-aza treatment and ectopically expressed MYOD1 (17). Although contractile porcine myotubes with sarcomeres were obtained from iPSCs within eleven days, the entire process would be more time-consuming because of the iPSC derivation period. Here, we established a direct transdifferentiation protocol with pig fibroblasts for the generation of skeletal muscle using porcine MYOD1 overexpression and myogenesis-associated signaling molecules, bypassing the induction of the pluripotent state.

Comparative analysis of the MYOD1 a.a sequence showed that the bHLH domain in MYOD1 was thoroughly conserved across various species. A previous study suggested that conserved sites are less permissive to evolutionary mutation due to their functional or structural importance (41). This finding indicates that the bHLH domain is an essential part of the MYOD1 protein, the function of which is associated with myogenesis, as in other species (34, 42). The basic region in the bHLH domain of Myod1 recognizes and binds the E-box, which is a conserved DNA sequence that is frequently distributed throughout the genome rather than specifically located in the regulatory region of myogenic genes (42). The myogenic specificity of Myod1 is derived from Ala¹¹⁴- Thr¹¹⁵ in the basic region (29, 43). In fact, it has been verified that Myod1 binds to both canonical and noncanonical E-boxes (10). These myogenic codes direct Myod1 to bind to noncanonical E-boxes of Myog, one of the Myod1 target loci, and to interact with Pbx/Meis cofactors associated with myogenic genes (30). In particular, Ala¹¹⁴ leads to an appropriate conformational change allowing myogenic activity by mediating the contact of Arg¹¹¹ and guanine in DNA (31). In addition to the basic region, the HLH motif of Myod1 is also required for myogenesis because the HLH motif dimerizes with other bHLH proteins, whose basic domains are involved in E-box binding (42). The other domains are also highly conserved and have been known to be functional Myod1 domains (11, 44). For example, the acidic domain acts as a transcriptional activation domain (TAD) through additional DNA binding near E boxes, and both the H/C and helix3 domains are involved in chromatin remodeling to allow active transcription of sequential myogenic genes. Accordingly, pig MYOD1 is capable of muscle-specific gene expression through conserved domains, especially the bHLH motif, preventing differentiation into nonmuscle cell fate.

Myogenic transdifferentiation is enhanced by activation of the FGF, WNT, and cAMP signaling pathways and inhibition of the TGF-β signaling pathway

Doxycyline (DOX)-inducible pMYOD1 overexpression vectors were produced with the conserved pig MYOD1 sequence, as described above. Because differentiation is blocked when proliferation is enhanced by mitogen in serum, a low concentration of horse serum has been widely employed for myogenic differentiation in humans and mice (45–47). However, transdifferentiated muscle cells were not observed when the aforementioned conventional culture conditions were applied. This result was consistent with the qPCR data that showed mild upregulation of endogenous muscle genes (Fig. 2C). It was surmised that myogenic reprogramming in pigs requires more supportive culture conditions than in mice and humans. For these reasons, a mixture of small molecules, such as FGF2, a TGF-β inhibitor (SB431542), a WNT activator (CHIR99021), and a cAMP activator (forskolin), was selected to facilitate the conversion of cell fate into myogenic lineage, as applied in previous research (8). While activation of the FGF, WNT, and cAMP signaling pathways is required for specification into the myogenic lineage and proliferation of committed myoblasts, inhibition of the TGF-β signaling pathway is involved in myotube formation (3, 13, 15). These signaling molecules, except forskolin, are secreted during myogenesis in vivo, thus recapitulating the endogenous signaling pathway for muscle formation and regeneration. Forskolin has been reported to be involved in skeletal muscle differentiation from human iPSCs and to improve satellite cell expansion in mice (16). Supporting the function of forskolin, our preliminary study showed that the removal of forskolin resulted in the failure of myogenic transdifferentiation (data not shown).

With treatment of these four signaling molecules, endogenous myogenic genes were highly triggered by exogenous MYOD1 (Exo-MYOD1) (Fig. 3C). Consistent with these observations, it has been demonstrated that ectopically expressed MYOD1 upregulates endogenous MYOD1 (Endo-MYOD1), PAX7, MYF5, and MYOG, which reinforces the notion that MYOD1 is a key regulator in myogenesis. Exo-MYOD1 increased Endo-MYOD1 via an autoregulatory loop, as previously reported (26). Chromatin immunoprecipitation reportedly showed that mesoderm or myoblast markers, including Pax7/PAX7 and Myog/MYOG, were directly activated by Exo-Myod1, leading to an increase in Myf5 (27, 48). Because MYF5 is upstream of MYOD1, it was reported that MYF5 was not expressed in MYOD1-overexpressing human iPSCs (33). However, when treated with signaling molecules, forced Myod1 expression converted mouse fibroblasts into Pax7- and Myf5-positive myogenic progenitor cells (8). The temporal expression of myogenic genes has a hierarchy with stage-specific markers: PAX7 in myogenic progenitor cells and MYF5 and MYOD1 in committed myogenic cells, followed by MYOG in the differentiation phase (3, 49). Therefore, various cell types that belong to the myogenic lineage can be derived in mitogen-rich media. In the C2C12 myoblast cell line, mitogen-rich culture conditions stimulated myoblast growth without differentiation, whereas skeletal muscle was differentiated with the expression of differentiation-specific genes under low-mitogen culture conditions (50). For the enrichment of terminal differentiated myotubes, the culture conditions were switched, and a medium containing 2% horse serum without signaling molecules on day 6 was used, in which endogenous muscle genes peaked by Exo-MYOD1.

As shown in the applied protocol presented in Fig. 4A, multinucleated and elongated myotubes were observed. These myotubes were formed through the fusion of mononucleated myoblasts exiting from the cell cycle, followed by reorganization of the cytoskeleton (3). Based on the higher expression of muscle genes in the + DOX group, continuous MYOD1 overexpression seems to enhance myogenic conversion in a feed-forward mechanism (42). Myotubes and myofiber express myosin heavy chain (MHC) (6), which is a downstream gene of Myod1 (26). MHC provides contractility to eukaryotic cells through filament assembly in the form of striated sarcomeres, such as skeletal and cardiac muscle (51). In fact, the expression of MHC was increased at the RNA level and a sarcomere-like structure was detected by immunostaining for MHC proteins in accordance with previous research (35). Altogether, our protocol enables the transdifferentiated muscle to undergo terminal differentiation and maturation into skeletal muscle.

In summary, fully differentiated somatic cells of pigs can be reprogrammed into mature skeletal muscle by the pig MYOD1 gene. Modulation of the FGF, TGF-β, WNT, and cAMP signaling pathways is required for the cell fate conversion into the myogenic lineage. The transdifferentiated muscle expressed skeletal muscle markers and had the structure of a striated sarcomere, implying that these matured myotubes possess contractile capacity. Given the role of pig as a significant livestock for supplying meat, the myogenic reprogramming of pig cells can be applied to increase agricultural yield and produce cultured meat. Additionally, pigs provide biomedical applications in preclinical studies for human disease because of their anatomical and physiological similarities with humans (52, 53). For example, transdifferentiated pig muscle cells can offer a cell source for skeletal muscle disease modeling and drug screening in regenerative medicine. Finally, this study provides fundamental knowledge for developmental biology in revealing the genetic network and signaling pathways underlying myogenesis.

PEFs: Pig embryonic fibroblasts

PSCs: Pluripotent stem cells

ESCs: Embryonic stem cells

iPSCs: Induced pluripotent stem cells

pMYOD1-PEFs: PEFs infected with a lentivirus carrying the pMYOD1 gene

ICC: Immunocytochemistry

qPCR: Quantitative real-time polymerase chain reaction

gDNA: genomic DNA

5-aza: 5-azacytidine

DOX: Doxycycline

FBS: Fetal bovine serum

KSR: Knockout serum replacement

FGF2: Fibroblast growth factor 2

Cont.: Control

FSCHF: FGF2, SB431542, CHIR99021, and forskolin

pACTB: Porcine beta actin

pMYOD1: Pig MYOD1

Exo-MYOD1: Exogenous MYOD1

Endo-MYOD1: Endogenous MYOD1

MHC: Myosin heavy chain

H/C: histidine and cysteine-rich

bHLH: basic-helix1-loop-helix2

TAD: transcriptional activation domain

Ethics approval and consent to participate

The care and experimental use of pigs were approved by the Institutional Animal Care and Use Committee (IACUC) at Seoul National University (approval nos.: SNU-140328-2 for pig fetal fibroblast isolation).

Consent for publication

Not applicable.

Availability of data and material

Not applicable.

Competing interests

No potential conflicts of interest relevant to this article are reported.

Funding

This work was supported by the BK21 Plus Program and the National Research Foundation of Korea (NRF) Grant funded by the Korea government (NRF-2019R1C1C1004514).

Authors’ contributions

J. Jeong, K.-H. Choi, S.-H. Kim, D.-K. Lee, J.-N. Oh, M. Lee, GC. Choe and C.-K. Lee designed research. J. Jeong and K.-H. Choi performed, analyzed and interpreted all experiments. J. Jeong and C.-K. Lee wrote the paper. C.-K. Lee finally approved the manuscript.

Acknowledgements

Not applicable.

Li K, Kong Y, Zhang M, Xie F, Liu P, Xu S. Differentiation of pluripotent stem cells for regenerative medicine. Biochemical biophysical research communications. 2016;471(1):1-4.
Choi Y-J, Park S-K, Kang H-I, Roh S-H. Effects of Exogenous Insulin-like Growth Factor 2 on Neural Differentiation of Parthenogenetic Murine Embryonic Stem Cells. Reproductive Developmental Biology. 2012;36(1):33-7.
Chal J, Pourquié O. Making muscle: skeletal myogenesis in vivo and in vitro. Development. 2017;144(12):2104-22.
Oh Y, Jang J. Directed differentiation of pluripotent stem cells by transcription factors. Molecules and cells. 2019;42(3):200.
Kodaka Y, Rabu G, Asakura A. Skeletal muscle cell induction from pluripotent stem cells. Stem cells international. 2017;2017.
Chal J, Al Tanoury Z, Hestin M, Gobert B, Aivio S, Hick A, et al. Generation of human muscle fibers and satellite-like cells from human pluripotent stem cells in vitro. Nature protocols. 2016;11(10):1833-50.
Chal J, Oginuma M, Al Tanoury Z, Gobert B, Sumara O, Hick A, et al. Differentiation of pluripotent stem cells to muscle fiber to model Duchenne muscular dystrophy. Nature biotechnology. 2015;33(9):962-9.
Bar-Nur O, Gerli MF, Di Stefano B, Almada AE, Galvin A, Coffey A, et al. Direct reprogramming of mouse fibroblasts into functional skeletal muscle progenitors. Stem cell reports. 2018;10(5):1505-21.
Davis RL, Weintraub H, Lassar AB. Expression of a single transfected cDNA converts fibroblasts to myoblasts. Cell. 1987;51(6):987-1000.
Berkes CA, Tapscott SJ, editors. MyoD and the transcriptional control of myogenesis. Seminars in cell & developmental biology; 2005: Elsevier.
Hernández-Hernández O, Ávila-Avilés RD, Hernández-Hernández JM. Chromatin Landscape During Skeletal Muscle Differentiation. Frontiers in Genetics. 2020;11:1136.
Charge SB, Rudnicki MA. Cellular and molecular regulation of muscle regeneration. Physiological reviews. 2004;84(1):209-38.
Jiwlawat N, Lynch E, Jeffrey J, Van Dyke JM, Suzuki M. Current progress and challenges for skeletal muscle differentiation from human pluripotent stem cells using transgene-free approaches. Stem Cells International. 2018;2018.
Nakajima T, Shibata M, Nishio M, Nagata S, Alev C, Sakurai H, et al. Modeling human somite development and fibrodysplasia ossificans progressiva with induced pluripotent stem cells. Development. 2018;145(16).
Chen AE, Ginty DD, Fan C-M. Protein kinase A signalling via CREB controls myogenesis induced by Wnt proteins. Nature. 2005;433(7023):317-22.
Xu C, Tabebordbar M, Iovino S, Ciarlo C, Liu J, Castiglioni A, et al. A zebrafish embryo culture system defines factors that promote vertebrate myogenesis across species. Cell. 2013;155(4):909-21.
Genovese NJ, Domeier TL, Telugu BPV, Roberts RM. Enhanced development of skeletal myotubes from porcine induced pluripotent stem cells. Scientific reports. 2017;7:41833.
Ding S, Wang F, Liu Y, Li S, Zhou G, Hu P. Characterization and isolation of highly purified porcine satellite cells. Cell death discovery. 2017;3(1):1-11.
Olson EN. MyoD family: a paradigm for development? Genes development. 1990;4(9):1454-61.
Choi KH, Yoon JW, Kim M, Jeong J, Ryu M, Park S, et al. Optimization of Culture Conditions for Maintaining Pig Muscle Stem Cells In Vitro. Food Sci Anim Resour. 2020;40(4):659-67.
Choi K-H, Park J-K, Kim H-S, Uh K-J, Son D-C, Lee C-K. Epigenetic changes of lentiviral transgenes in porcine stem cells derived from embryonic origin. PloS one. 2013;8(8):e72184.
Choi K-H, Park J-K, Son D, Hwang JY, Lee D-K, Ka H, et al. Reactivation of endogenous genes and epigenetic remodeling are barriers for generating transgene-free induced pluripotent stem cells in pig. PloS one. 2016;11(6):e0158046.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. methods. 2001;25(4):402-8.
Taylor SM, Jones PA. Changes in phenotypic expression in embryonic and adult cells treated with 5‐azacytidine. Journal of cellular physiology. 1982;111(2):187-94.
Weintraub H, Tapscott SJ, Davis RL, Thayer MJ, Adam MA, Lassar AB, et al. Activation of muscle-specific genes in pigment, nerve, fat, liver, and fibroblast cell lines by forced expression of MyoD. Proceedings of the National Academy of Sciences. 1989;86(14):5434-8.
Bichsel C, Neeld D, Hamazaki T, Chang L-J, Yang L-J, Terada N, et al. Direct reprogramming of fibroblasts to myocytes via bacterial injection of MyoD protein. Cellular reprogramming. 2013;15(2):117-25.
Kabadi AM, Thakore PI, Vockley CM, Ousterout DG, Gibson TM, Guilak F, et al. Enhanced MyoD-induced transdifferentiation to a myogenic lineage by fusion to a potent transactivation domain. ACS synthetic biology. 2015;4(6):689-99.
Fong AP, Yao Z, Zhong JW, Johnson NM, Farr Iii GH, Maves L, et al. Conversion of MyoD to a neurogenic factor: binding site specificity determines lineage. Cell reports. 2015;10(12):1937-46.
Davis RL, Cheng P-F, Lassar AB, Weintraub H. The MyoD DNA binding domain contains a recognition code for muscle-specific gene activation. Cell. 1990;60(5):733-46.
Heidt AB, Rojas A, Harris IS, Black BL. Determinants of myogenic specificity within MyoD are required for noncanonical E box binding. Molecular cellular biology. 2007;27(16):5910-20.
Huang J, Weintraub H, Kedes L. Intramolecular regulation of MyoD activation domain conformation and function. Molecular cellular biology. 1998;18(9):5478-84.
Xie X, Fu Y, Liu J. Chemical reprogramming and transdifferentiation. Current Opinion in Genetics Development. 2017;46:104-13.
Tanaka A, Woltjen K, Miyake K, Hotta A, Ikeya M, Yamamoto T, et al. Efficient and reproducible myogenic differentiation from human iPS cells: prospects for modeling Miyoshi Myopathy in vitro. PloS one. 2013;8(4):e61540.
Conerly ML, Yao Z, Zhong JW, Groudine M, Tapscott SJ. Distinct activities of Myf5 and MyoD indicate separate roles in skeletal muscle lineage specification and differentiation. Developmental cell. 2016;36(4):375-85.
Jiwlawat S, Lynch E, Glaser J, Smit-Oistad I, Jeffrey J, Van Dyke JM, et al. Differentiation and sarcomere formation in skeletal myocytes directly prepared from human induced pluripotent stem cells using a sphere-based culture. Differentiation. 2017;96:70-81.
Lassar AB, Paterson BM, Weintraub H. Transfection of a DNA locus that mediates the conversion of 10T12 fibroblasts to myoblasts. Cell. 1986;47(5):649-56.
Stresemann C, Lyko F. Modes of action of the DNA methyltransferase inhibitors azacytidine and decitabine. International journal of cancer. 2008;123(1):8-13.
Taylor SM, Jones PA. Multiple new phenotypes induced in 10T12 and 3T3 cells treated with 5-azacytidine. Cell. 1979;17(4):771-9.
Kiziltepe T, Hideshima T, Catley L, Raje N, Yasui H, Shiraishi N, et al. 5-Azacytidine, a DNA methyltransferase inhibitor, induces ATR-mediated DNA double-strand break responses, apoptosis, and synergistic cytotoxicity with doxorubicin and bortezomib against multiple myeloma cells. Molecular cancer therapeutics. 2007;6(6):1718-27.
Ito N, Kii I, Shimizu N, Tanaka H, Takeda Si. Direct reprogramming of fibroblasts into skeletal muscle progenitor cells by transcription factors enriched in undifferentiated subpopulation of satellite cells. Scientific reports. 2017;7(1):1-12.
Capra JA, Singh M. Predicting functionally important residues from sequence conservation. Bioinformatics. 2007;23(15):1875-82.
Tapscott SJ. The circuitry of a master switch: Myod and the regulation of skeletal muscle gene transcription. Development. 2005;132(12):2685-95.
Brennan TJ, Chakraborty T, Olson EN. Mutagenesis of the myogenin basic region identifies an ancient protein motif critical for activation of myogenesis. Proceedings of the National Academy of Sciences. 1991;88(13):5675-9.
Hernández-Hernández JM, García-González EG, Brun CE, Rudnicki MA, editors. The myogenic regulatory factors, determinants of muscle development, cell identity and regeneration. Seminars in cell & developmental biology; 2017: Elsevier.
Kim B, Ko N-Y, Hwang S-s, Kim D-H, Park J-K, Ryoo Z-Y, et al. Growth Factors Supplementation in Culture Medium Leads to Active Proliferation of Porcine Fibroblasts. Reproductive Developmental Biology. 2011;35(3):301-6.
Kimura E, Han JJ, Li S, Fall B, Ra J, Haraguchi M, et al. Cell-lineage regulated myogenesis for dystrophin replacement: a novel therapeutic approach for treatment of muscular dystrophy. Human molecular genetics. 2008;17(16):2507-17.
Lattanzi L, Salvatori G, Coletta M, Sonnino C, De Angelis MC, Gioglio L, et al. High efficiency myogenic conversion of human fibroblasts by adenoviral vector-mediated MyoD gene transfer. An alternative strategy for ex vivo gene therapy of primary myopathies. The Journal of clinical investigation. 1998;101(10):2119-28.
Gianakopoulos PJ, Mehta V, Voronova A, Cao Y, Yao Z, Coutu J, et al. MyoD directly up-regulates premyogenic mesoderm factors during induction of skeletal myogenesis in stem cells. Journal of Biological Chemistry. 2011;286(4):2517-25.
Tomczak KK, Marinescu VD, Ramoni MF, Sanoudou D, Montanaro F, Han M, et al. Expression profiling and identification of novel genes involved in myogenic differentiation. The FASEB journal. 2004;18(2):1-23.
Cao Y, Yao Z, Sarkar D, Lawrence M, Sanchez GJ, Parker MH, et al. Genome-wide MyoD binding in skeletal muscle cells: a potential for broad cellular reprogramming. Developmental cell. 2010;18(4):662-74.
Weiss A, McDonough D, Wertman B, Acakpo-Satchivi L, Montgomery K, Kucherlapati R, et al. Organization of human and mouse skeletal myosin heavy chain gene clusters is highly conserved. Proceedings of the National Academy of Sciences. 1999;96(6):2958-63.
Choi K-H, Lee C-K. Pig pluripotent stem cells as a candidate for biomedical application. J Anim Reprod Biotechnol. 2019;34(3):139-47.
Ezashi T, Yuan Y, Roberts RM. Pluripotent stem cells from domesticated mammals. Annual review of animal biosciences. 2016;4:223-53.

Download PDF

Journal Publication

published 13 May, 2021

Read the published version in Journal of Animal Science and Biotechnology →

Editorial decision: Major revision
06 Dec, 2020
Review #1 received at journal
22 Nov, 2020
Review #2 received at journal
22 Nov, 2020
Reviewer #2 agreed at journal
08 Nov, 2020
Reviewers invited by journal
08 Nov, 2020
Reviewer #1 agreed at journal
08 Nov, 2020
First submitted to journal
02 Nov, 2020
Editor assigned by journal
02 Nov, 2020
Submission checks completed at journal
02 Nov, 2020
Editor invited by journal
02 Nov, 2020

You are reading this latest preprint version

Myogenic Transdifferentiation of Pig Fibroblasts Using the Porcine MYOD1 Inducible Vector

Status:

Journal Publication

Version 1

Abstract

Figures

Background

Methods

Amino acid sequence analysis

Lentiviral vector production

Myogenic transdifferentiation of pig fibroblasts

Genomic DNA (gDNA) extraction and polymerase chain reaction (PCR)

Immunocytochemistry (ICC) analysis

Quantitative real-time polymerase chain reaction (qPCR)

Statistical analysis

Results

Discussion

Pig MYOD1 has an evolutionarily conserved bHLH domain that regulates myogenesis

Conclusions

List of Abbreviations

Declarations

References

Status:

Journal Publication

Version 1