Isolation, proliferation and differentiation for the PSCs and PADSCs.
Most of the muscles of large yellow croaker belong to the white muscle of vertebrates. It is mainly composed of skeletal muscle and fat. Muscle satellite cells are considered as the main stem cells and are crucial for muscle growth and regeneration. They are generally quiescent until to be activated and become myogenic cells, i.e., myoblasts, which will undergo proliferation and myogenic differentiation to form new muscle fibers.
The culture of large yellow croaker muscle cells has been reported recently, where tissue block culture was used12. To isolate more pure muscle stem cells, we followed an enzymolysis approach from higher vertebrates and other fishes13. The epaxial muscle of large yellow croaker was physically decomposed. The satellite cells were released by collagenase and trypsin treatments and collected by filtration (Fig. 1a). By optimizing the combination of different media, serum ratio and growth factors (data not shown), the proliferation medium of PSCs was finally determined as high-glucose DMEM with 15% fetal bovine serum (FBS) and 10 ng/ml bFGF. At 9 days post-incubation, the primary cultured cells reached about 80% confluence, and the cell morphology began to appear as fibroblasts. After 3 passages, these cells gradually displayed typical fibroblasts morphology (Fig. 1b). Cells from Passage 3 were used to evaluate their proliferation ability with CCK8 assay kit. The data showed that the number of cells doubled less than 2 days after seeding, and the proliferation plateau began to appear on the third day (Fig. 1c). This proliferation capacity was maintained for at least 13 passages (data not shown here). The cells from Passage 3 were also characterized by immunofluorescent staining with the specific myoblast markers such as Pax7 and MyoD1. The results showed that about 43% of the cells were Pax7 positive, and 99% of the cells were MyoD1 positive (Fig. 1d, e). These data suggested that the isolated cells did possess the properties of PSCs.
To evaluate the differentiation potentials of these PSCs, even after multiple passages, the serum starvation myogenic differentiation protocol for mammalian satellite cells was initially adopted. Although myotubes were induced within 2–3 days with medium containing 2% horse serum (HS), they detached shortly and died soon after. To improve differentiation efficiency and survival rate, different combinations of basal medium, serum ratios and myogenic factors were evaluated. Finally, the F12 medium containing 8% HS, 10 ng/ml IGF-1, 50 nM necrosulfonamide and 200 µM ascorbic acid was identified as an effective recipe for myogenic differentiation. On day 3 of culture in this differentiation medium, elongated myotubes began to appear, and on day 6, striations similar to skeletal muscle were observed in some of myotubes (Fig. 1f). Immunofluorescent staining of Desmin (a muscle-specific protein) showed that Desmin-positive myotubes with multiple-nuclei were detected on both days 3 and 6 (Fig. 1f, g). For unexpected reasons, the fusion index (the proportion of nuclei in multinucleated myotubes) was higher on day 3 (6.39%) than on day 6 (1.30%) (Fig. 1g). Although the average number of nuclei in individual myotubes was higher on day 6, the total number of myotubes on day 6 was much less than on day 3 (Fig. 1f).
The isolation, proliferation and differentiation of large yellow croaker preadipocytes have already been reported14. We followed the protocol to isolate PADSCs from ventral coelomic fat tissues. Similar to PSCs, the culture conditions for PADSC proliferation and differentiation were optimized. The proliferation medium was high-glucose DMEM containing 8% FBS, which was reduced by half compared with that previously reported14. The cell growth rate maintained high in the first two days and reached a plateau after three days (Fig. 1c). Immunofluorescence staining of adipose tissue marker HoxC9 showed that almost all cultured cells were positive (Fig. 1i), suggesting that these cells originated from abdomen fat tissue were PADSCs. The feature of PADSCs was confirmed by the adipogenic differentiation. DMEM/F12 containing 10% HS, 10 µg/mL insulin, 0.5 µm IBMX, 0.25 µm dexamethasone and 1% Lipid Mixture Solution, was used as the differentiation medium. Staining with Oil Red O showed that lipid droplets began to appear in almost all cultured cells on day 1, and the number and size of lipid droplets increased obviously on day 6 (Fig. 1j). These results demonstrated that the isolated cells were PADSCs.
Improving The Myogenic Differentiation Efficiency By Transcriptomic Analysis
In addition to the low efficiency of myogenic differentiation, another major concern in our PSC differentiation culture was the decrease of the fusion index from day 3 to day 6 (Fig. 1g). To address these questions, we performed transcriptomic profiling by bulk RNA-seq at different time points (day 0, 3, 6) of differentiation. Clustering analysis showed that the expression profile changes on day 3 were more dramatical than those on day 6, compared with day 0 (Fig. 2a), which was consistent with the fusion index. The differentially expressed genes (DEGs) analysis showed that 1,650 of 22,093 genes were significantly upregulated (log2 FC > 1, FDR < 0.05), and 1,063 were significantly downregulated (log2 FC < − 1, FDR < 0.05) from day 0 to day 3 (Fig. 2b). Gene ontology (GO) enrichment revealed that the upregulated genes were predominantly related to terms associated with muscle development and function, whereas most of the downregulated genes were enriched under terms such as negative regulation of cell differentiation, developmental process and some metabolic process, suggesting that serum starvation induced myogenic differentiation (Fig. 2c). However, we noted that genes related to cell cycle regulation were not in the downregulated gene category, which was different from that previously described in mice and bovine15. There were at least two reasons for phenomenon. The first one was the low frequency of myogenic differentiation of PSCs. The second reason was the low viability of fish myocytes in the medium lacking serum (data not shown).
The results of DEGs from comparison showed that 664 genes were significantly upregulated, and 1,165 were significantly downregulated between day 6 and day 3 (Supplementary Fig. 1a), wherein the term of mitotic cell cycle process was also not enriched into the downregulated gene category (Supplementary Fig. 2a). While, the term of mitotic cell cycle process was enriched into the downregulated gene category form DEGs between day 6 and day 0 (Supplementary Fig. 2b), where 1,636 genes were significantly upregulated, and 1,428 were significantly downregulated (Supplementary Fig. 1b). These data indicated that cell proliferation gradually decreased during the differentiation from day 0 to day 6 (Fig. 2e).
Gene expression profiles of myogenic factors and differentiation marker genes showed that the myogenic stem/progenitor marker pax7 and the early myogenic determination marker myod1 were continuously downregulated from day 0 to day 6 of differentiation; the late myogenic determination marker myog, the myogenic differentiation markers desmin and the myoblast fusion factor mymk were upregulated on day 3 and then downregulated on day 6; the skeletal-muscle-specific markers myh11 and mylk steadily increased from day 0 to day 6 (Fig. 2d, Supplementary Fig. 1c). Interestingly, the myogenic regulator myf5 was not detected in the whole process. This might be because the PSCs isolated from epaxial muscle of large yellow croaker were Myf5−. The vertebrate epaxial muscle originates from embryonic dermomyotome, in which there is developmental heterogeneity of the cell lineage. At least two populations of myoblast progenitors coexist in developing dermomyotome. The major population is the fast-cycling population, Pax7+/Myf5+, which possesses limited differentiation capacity and will participate in terminal differentiation more quickly16. The minor population is the slow-cycling population, Pax7+/Myf5−, with self-renew capability, which is important to maintain muscle homeostasis and regeneration17, 18, 19. Our results demonstrated that the PSCs isolated from fish epaxial muscle belong to Pax7+/Myf5− slow-cycling population.
In order to find out what signals might be responsible for low myogenic efficiency, we performed KEGG analysis. The results showed that a number of important signaling pathways, including MAPK, TGF-β, Wnt, Notch, hippo, HIF-1, etc., were enriched in the DEGs (Supplementary Fig. 1c). Among them, two signaling pathways, Notch and TGF-β, attracted our attention in particular. In muscle tissue, Notch signaling is required for the maintenance of quiescence of satellite cells20. TGF-β signaling plays a negative role in myogenesis of adult muscle tissue (reviewed by Burks & Cohn21). Balance between Notch and TGF-β signaling regulates proliferation of satellite-cells, and attenuation of TGF-β restores regeneration in old, injured muscle22. In the differentiation process of PSCs, the expression of main members of Notch pathway on day 3 such as lfng, hes4, hey1 and hes1, did not significantly decrease, and the expression of notch1 even increased, though all of these genes decreased on day 6 (Fig. 2f). Since the downregulation of Notch signaling at the early stage is crucial for PSCs to transition from stem cell fate into differentiation state, relatively high level of Notch signal on day 3 might be one of obstacles for myogenic differentiation.
The changes of TGF-β signaling were more complicated in myogenic differentiation of PSCs (Fig. 2g). Compared to day 0, on day 3, the expression of some members of TGF-β signaling (dcn, tgffbr3, tmoda and smad3) increased, while the expression of some members (smad6, smad2, tgfb1a, tgfb5 and tgfbr1b) decreased; on day 6, the expression of dcn, tgfbr2 and smad7 increased, while the expression of smad6 and smad2 decreased (Fig. 2g). The data indicated that TGF-β signaling was dynamically changed during the induction of myogenic differentiation. A previous study reported that TGF-β induced myogenic cells to transform into myofibroblastic cells after injury23. Therefore, we wondered whether the low fusion index was due to TGF-β signal induced transformation of PSCs into myofibroblasts in the serum starvation culture. The expression of α-smooth muscle actin (α-sma) (a specific marker for myofibroblast) and most extracellular-matrix (ECM) components was upregulated on day 3 and upregulated further on day 6 (Fig. 2h, i), indicating that myofibroblastic differentiation occurred in the serum starvation culture.
To improve the myogenic efficiency and prevent myofibrosis, we inhibited Notch and TGF-β signals with two small molecule drugs: LY411575 (a Notch inhibitor) and RepSox (a TGFβR-1/ALK5 inhibitor) respectively. Two chemicals were added to the basal medium (F12 medium containing 8%HS, 1×PS) for the culture of first two days and DMEM medium (containing 8% FBS, 200nM dexamethasone and 1×PS) for the culture of later stage. On day 6, the numbers of myotubes were increased slightly in the treatments with either LY411575 or RepSox alone, and increased much obviously in the treatment with both chemicals together (Fig. 2k). The fusion index in the control was about 1%, which was increased to 8% and 3.5% by the treatment of RepSox or LY411575 respectively, and synergistically increased to 32% by the treatment of two chemicals together (Fig. 2k). The results demonstrated that the myogenic efficiency of PSCs was greatly improved by inhibition of Notch and TGF-β signals.
Effects Of Supporting Materials Of 3d-scaffold On Fish Cell Adhesion And Growth
To select the supporting materials for 3D-culture of PSCs, we compared various materials from different resources including porcine and fish gelatins (PG/FG)24, hyaluronic acid (HA)25, silk fibroin (SF)26, and chitosan27. The evaluation of material-to-cell compatibility showed that gelatin-based gel had the highest cell viabilities (Fig. 3a and Supplementary Fig. 3). As a partial hydrolysate of collagen, gelatin has natural biocompatibility and can provide an extracellular microenvironment suitable for the survival of fish cells28. The best temperature for fish cell growth is about 27°C29. To maintain the proliferation capacity of fish cells, the gel point of 3D culture materials for large yellow croaker cells should be lower than ~ 27°C, which is different from conventional bio-inks for human or other mammalian cells (about 37°C). Therefore, we made cell culture gels by using porcine- and fish-derived gelatin combined with sodium alginate (SA) to measure their rheological properties. The gelation temperature of FG material was lower than that of PG material. The gel points (i.e. G’=G’’) of 5% and 10% FG were 22.15 and 25.51°C, respectively (Fig. 3b). PG inherently contains higher amounts of hydrophobic amino acids and forms stronger hydrophobic interactions than FG24, resulting in higher gel points (31.25 and 34.39°C for 5% and 10% PG, respectively) which could not be used for fish cell culture. The gel with 10% FG/1% SA had relatively stable value of viscosity, though similar trends of viscosity decrease were observed for both PG and FG materials with SA at 0.1–100 rad/s of shear rate (Fig. 3c).
Subsequently, the gel (10% FG/1% SA) were mixed with PSCs and treated with CaCl2 at different concentrations for different times to pre-fabricate scaffold materials with altering stiffness of elastic modulus. We found that 1% CaCl2 solidification for 10 min or 5% CaCl2 solidification for 5 min induced most PSCs to adhere (Fig. 3d). It is known that large changes in matrix stiffness influence the focal-adhesion and cytoskeleton structures for differentiated cells30, 31. Mechanosensing receptors on the PSC membrane activated actomyosin-mediated contraction and cytoskeletal rearrangement after sensing matrix stiffness in a narrow range of ~ 9–14 kPa (Fig. 3e, f). The cell adhesion rates could reach up to 80.69% for the PSCs treated with 1% CaCl2 for 10 min (Fig. 3f). The activated PSCs became spindle-shaped and gradually proliferated, migrated and differentiated in the 3D culture materials, while non adherent cells died after a few days.
We observed that PSCs had low cell proliferation rate in the 3D culture environment. To deal with this problem, a p53 (a well-known tumor repressor) inhibitor (Pifithrin-α) and a Yap (a positive regulator in cell survival and proliferation) activator (XMU-MP-1) were applied in the culture. The results showed that the cell density in the control medium was about 8.02×104 cells/µL, which was increased to 8.31×, 10.73× and 11.57×104 cells/µL by the treatment of Pifithrin or XMU-MP-1 or two chemicals together (Fig. 3g). The data demonstrated that the cell proliferation rate of PSCs was significantly improved by the combination of two chemicals.
Myogenic differentiation experiments showed that the fusion index in the 3D culture with the basal medium (F12 medium containing 8%HS, 1×PS) was 9.94%, which was much higher than that (about 0.91%) in 2D culture, suggesting that cells in 3D culture inherently have a higher fibrogenic tendency than cells in 2D culture due to the stimulation of supporting material stiffness. Consistent with 2D culture, combination of LY411575 (Notch inhibitor) and RepSox (TGFβR-1/ALK5 inhibitor) also significantly increased the fusion index in 3D differentiation culture (from 9.94–41.79% on day 7) (Fig. 2k and Fig. 3h).
Muscle Tissue Formation Under Guidance Of Biomimetic Model Of Scaffold
The reconstruction of the microstructure of cultured meat, especially based on its real biological structure, has been a challenge due to the reproduction of tissue-like texture. In addition to the supporting role, layered porous scaffolds can also enhance cell adhesion and proliferation, which indicates the importance of microstructure to the biological function of scaffolds32, 33. On the purpose of reconstruction of a printable optimized model of scaffolds, we used a digital analysis of micro-CT scanning and computerized tools to analyze the textures of large yellow croaker muscle tissue, and conducted a 3D bioprinting method to manufacture the biomimic scaffolds (Fig. 4a).
The 3D structure of the micro-CT slice images of large yellow croaker epaxial muscle tissue showed that fish muscle fibers were an unregular polygonal shape and were surrounded by a layer of myocommata composed of collagen and syndesm (Fig. 4b and Supplementary Fig. 4a). Fat tissues as intermuscular filler tissues distributed between muscle fibers (Fig. 4b). 3D volume rendering gave a brand-new view of large yellow croaker muscle compared to conventional imaging (data not shown). Muscle fibers and myocommata seemed to be anisotropic diffusion in space, forming a relatively gradual change in 3D structure. However, for high-precision scaffold construction, it was very difficult to adjust the irregular trajectory and needle diameter for 3D bioprinting of myofibrils34. Therefore, in order to simulate muscle fibers, we had to simplify the digital information of fish tissue structure to have a periodic superposition of two-dimensional plane with proper size (see below). The simplified flows were showed in a case of cross section of epaxial muscle tissue (Supplementary Fig. 4b), which guides the definition of longitudinal biomimic scaffolds shown in Fig. 4b. Due to the limitations in methods and equipments of 3D bioprinting for biomimetic and edible cell scaffolds, especially the bionic structure with spatial anisotropy, we designed a regular lattice model based on the calculation of fish micro-CT slice images (Fig. 4b and Supplementary Fig. 4c, d). According to the image processing and numerical transformation35, the digital muscle fibers were extracted in parallel or interlaced (with angles of 30°), and labeled as muscle scaffold model (Fig. 4b). The anisotropic lattice structure was rendered blue representing fat and syndesm. The top view of intermuscular filler model was shown in Fig. 4b, whose inverse shape was similar to the microstructure of real large yellow croaker muscle (front view in orange). By this, we obtained the ratio of muscle and fat (approximately 6:4) for modeling.
Then a 3D extrusion-printer was used to construct large yellow croaker muscle scaffold with the gel compound mixed with PSCs at 24℃. Each layer of the muscle scaffold consisted of 17 gel fibers, six of which were inclined at 30° to the other gel fibers (Fig. 4b). The printing was repeated ten times to form a muscle model with a length of 20 mm, a width of 12 mm and a thickness of 4 mm. The muscle scaffold was manufactured in batch process and placed in the proliferation medium. On day 3, spindle shaped PSCs were observed without ordered growth direction on biomimic muscle scaffold; On day 10, proliferated PSCs almost filled the muscle scaffold, and then the culture was replaced with differentiation medium (Fig. 4c). On day 7 of differentiation culture, myotubes were formed and integrated into scaffolds with aligned rearrangement (Fig. 4d). These aligned muscle fibers consisted of up to 5.35×107 of PSCs, with the cell growth rate maintained high from day 4 to day 10, and reached a plateau after seven days of differentiation (Fig. 4e).
Characterization Of Tissue-like Cultured Fish Fillets
The differentiated PADSCs of adipocytes were filled into the muscle scaffolds to form cell cultured fish fillets (Fig. 5a). To evaluate cultured meats, first, we compared the numbers and proportions of muscle cells and adipocytes in the cultured fish meats and muscle tissues of raw large yellow croaker (Fig. 5b, c). The results showed that the average numbers of muscle cells and adipocytes were respectively 5.67×107 and 4.02×107 in about 0.96 cm3 cultured fish fillet, which were slightly higher than those in the same volume of native muscle tissue (4.61×107 and 3.22×107, respectively).
Next, we used histological staining and SEM to analyze the microstructure of cultured fish meats. Hematoxylin-eosin and immunofluorescence section staining showed that after 17 days culture, the muscle fibers exhibited a similar arrangement to the native muscle fibers (Fig. 5d and Supplementary Fig. 5a). The muscle fibers in the scaffold were composed of myotubes with a diameter of about 400 µm. The ratio of muscle to fat in cell cultured fish fillets calculated from the cross-sections was 1.73:1, which was similar to that of native muscle tissue (1.88:1).
Furthermore, we performed textural analysis on empty scaffolds and scaffolds seeded with muscles/adipocytes. Textural measurements showed that the hardness of cultured fish fillet was 5.37 ± 0.42 N, which was in the same range as that of native yellow croaker fillet (5.72 ± 0.80 N) (Fig. 5e). There were also no obvious differences in gumminess, resilience and springiness between cultured meats and native muscle tissues. The hardness, gumminess and chewiness were increased, and springiness and resilience were decreased in scaffolds seeded with muscles/adipocytes, compared to empty scaffolds. However, chewiness of cultured fish fillets was lower than that of native muscle tissues. It might be due to the relatively higher free water content in cultured fish fillets (see below).
To explore the distribution and variation of water in cultured fish fillet, we analyzed the mobility and proportion of water molecules by low-field NMR relaxation behavior. By multi-exponential fitting analysis36, the cultured meats showed three peaks at 0.25 ± 0.00 ms, 2.83 ± 0.62 ms, and 335.87 ± 26.40 ms, representing T2b, T22, and T23 respectively, which accounted for tightly bound water, immobilized water and free water (Fig. 5f). Although empty scaffolds and native muscle tissues also had these three peaks as cultured meats, the contents of free water and bound water were different in the three samples. The free water content in cultured meats decreased by 8.23%, while the tightly bound water content increased by 9.52%, compared to those in empty scaffold, indicating that the muscle fibers formed in the cultured fish fillets transformed some of the free water into bound water. However, the water distribution in cultured fish fillets (11.10% tightly bound water, 0.95% immobilized water and 87.95% free water) still differed from that in native muscle tissue (3.33% tightly bound water, 96.48% immobilized water and 0.19% free water) (Supplementary Fig. 5c). Thus, the simulation of ECM in cultured fish meats could be focused in future work.
Taken together, the appearance and many characteristics (including total cell numbers, ratio of muscle cells and adipocytes, hardness, gumminess, resilience, springiness, etc.) of our cultured fish fillets were similar to those of native muscle tissues of large yellow croaker.