- Department of Cell Biology, Physiology and Immunology, Faculty of Biology, University of Barcelona, Barcelona, Spain
Fish muscle regeneration is still a poorly known process. In the present study, an injury was done into the left anterior epaxial skeletal muscle of seventy 15 g gilthead sea bream (Sparus aurata) juveniles to evaluate at days 0, 1, 2, 4, 8, 16 and 30 post-wound, the expression of several muscle genes. Moreover, transcripts’ expression in the bone (uninjured tissue) was also analyzed. Histology of the muscle showed the presence of dead tissue the first day after injury and how the damaged fibers were removed and replaced by new muscle fibers by day 16 that kept growing up to day 30. Gene expression results showed in muscle an early upregulation of igf-2 and a downregulation of ghr-1 and igf-1. Proteolytic systems expression increased with capn2 and ctsl peaking at 1 and 2 days post-injury, respectively and mafbx at day 8. A pattern of expression that fitted well with active myogenesis progression 16 days after the injury was then observed, with the recovery of igf-1, pax7, cmet, and cav1 expression; and later on, that of cav3 as well. Furthermore, the first days post-injury, the cytokines il-6 and il-15 were also upregulated confirming the tissue inflammation, while tnfα was only upregulated at days 16 and 30 to induce satellite cells recruitment; overall suggesting a possible role for these molecules as myokines. The results of the bone transcripts showed an upregulation first, of bmp2 and ctsk at days 1 and 2, respectively; then, ogn1 and ocn peaked at day 4 in parallel to mstn2 downregulation, and runx2 and ogn2 increased after 8 days of muscle injury, suggesting a possible tissue crosstalk during the regenerative process. Overall, the present model allows studying the sequential involvement of different regulatory molecules during muscle regeneration, as well as the potential relationship between muscle and other tissues such as bone to control musculoskeletal development and growth, pointing out an interesting new line of research in this group of vertebrates.
1 Introduction
The gilthead sea bream (Sparus aurata, Linnaeus 1758) is one of the most important species in the Mediterranean aquaculture of the last decades as it represents a source of energy and proteins of high quality and thus, it has a high commercial value (1). However, there are still areas of their physiology that are not well described; therefore, more research is needed to optimize their growth and development in terms of farming.
In contrast to other vertebrate groups, some fish species exhibit indeterminate growth, which means that they can grow in length and weight beyond the maturation stage for as long as they live (2). Most of this growth is due to an increase in somatic tissue, specifically skeletal muscle, and implies new muscle cells recruitment to increase the number of muscle fibers (hyperplasia), as well as the increase in size of pre-existing fibers (hypertrophy). These phenomena are known to be an important determinant of the characteristics of flesh texture and, from an aquaculture perspective, muscle quality is of main interest (3, 4).
The growth hormone (Gh) and insulin-like growth factors (Igfs) axis is the most important endocrine system controlling many physiological processes of the body (5). The Gh is the central regulatory molecule of somatic growth in vertebrates, including teleost fish, and although its main target is the liver it can exert significant direct effects on other tissues such as skeletal muscle. However, it is well known that most of the effects of Gh on growth are indirect, through the Igfs (6). The large tissue distribution of Igfs-producing cells including their receptors, together with the extensive distribution of Gh receptors (Ghrs), make this mediation possible. Many studies, including several in fish, have demonstrated systemic effects as well as the existence of local paracrine/autocrine actions of Gh and Igfs (7, 8). In muscle, Igfs directly stimulate a variety of signals including proliferation and differentiation of satellite cells (SCs), hypertrophy of myoblasts, and inhibition of cell atrophy. In gilthead sea bream, it has been shown that Igf-2 alone appears to be an important stimulator of proliferation, while Igf-1 seems to have a synergetic activity with Gh, playing different roles in the control of myogenesis (9). Moreover, in the same species, three splice variants have been described (igf-1a, igf-1b and igf-1c) with differential responses depending on the physiological or stimulatory condition (10). Igfs exert these actions through the binding and activation of their cellular receptors and in turn multiple intracellular signal transduction cascades, including the phosphatidylinositol 3-kinase (Pi3k)-Akt-Tor cascade, which controls an important gene expression program (11).
The process of muscle generation, known as myogenesis, consists in that SCs, the muscle stem cells that can be recognized for the expression of Pax3/7 (12), proliferate and differentiate into myoblasts that fuse together (i.e., hyperplasia) or with pre-existing fibers (i.e., hypertrophy) to form multinucleated myofibers. Myogenesis is regulated by the myogenic regulatory factors (Mrfs), a set of four transcription factors that are expressed sequentially being Myf5 and MyoD the first to be activated controlling SCs activation and proliferation, and later Myogenin and Mrf4 that coordinate myoblasts differentiation and fusion (13). In adulthood, those processes can also occur in response to challenging conditions that threaten muscle homeostasis. Thanks to the ability of muscle cells to modify both, the amount and types of proteins being synthesized within them, defined as muscle plasticity, it is possible to cope with any stimulus that disrupts the normal condition of the organism (14). In this framework, muscle regulation implies a protein turnover to ensure the correct function of the myofibers, i.e., a balance between protein synthesis and degradation. The main catabolic systems that control muscle proteolysis have been also described in gilthead sea bream and are the ubiquitin-proteasome complex (UbP); calpains, a family of Ca2+-dependent proteases formed by catalytic and regulatory small subunits; cathepsins, some of them related to the autophagy-lysosome system; and caspases, the apoptosis protease system (15, 16).
One of the phenomena that can alter normal muscle functions and body homeostasis are plasma membrane disruptions induced by muscle contraction as well as severe damage caused by traumatic injuries or muscle diseases, leading all these processes to fiber rupture. Nevertheless, skeletal muscle presents an intrinsic mechanism capable of regenerating the entire damaged contractile system. The fundamental role in muscle regeneration is played by the SCs resident in the muscle, which are activated in response to various previous stimuli that may come from the tissue itself or from other systems of the organism (17). When necrosis of the affected myofibers occurs because of the trauma, first an inflammatory response is initiated with the release of different cytokines, as the interleukins, from neutrophils, among other immune cells (18, 19). These induce macrophages infiltration and activate the proteolytic systems to remove all the dead material, followed by an activation of SCs to generate new fibers or to fuse with undamaged fibers. Finally, newly formed fibers maturate and remodel the regenerated muscle to recover its normal functions (20).
In addition, there are other factors that can also have critical roles during myogenesis. Proliferating SCs are characterized by the presence of proliferating cell nuclear antigen (Pcna), which is classically known to stimulate the entry of SCs to the cell cycle (21), and of hepatocyte growth factor receptor (cMet) that turns on also cell proliferation in skeletal muscle (22, 23). On the other hand, two novel muscle-specific fusogens, Myomaker and Myomixer, have been described and could play a fundamental role in the stage in which myoblast fusion takes place. In fact, the lack of some of these two muscle-specific membrane proteins results in an incomplete myogenesis since cells cannot fuse and therefore multinucleated myofibers are not formed (24). Other factors with an important role are caveolae, which are invaginations of the plasma membrane that function as messenger centers in tissues, including skeletal muscle, controlling signal transduction. The transmembrane proteins that enable the formation of these structures are members of the caveolin family (25). Caveolin-1 (Cav-1) was the first member described and was identified as being expressed in quiescent SCs but not in mature myofibers and is downregulated when cells begin to proliferate (26). The other two members of the family are Caveolin-2 (Cav-2) and Caveolin-3 (Cav-3). Characterization of the former revealed that colocalizes with Cav-1 and is co-expressed in the same circumstances. Otherwise, Cav-3, unlike Cav-1, is involved in myoblast differentiation and myotube fusion, so that a deficiency of Cav-3 results in cellular immaturity (27). Finally, another family of proteins essential for the maintenance of skeletal muscle homeostasis in adults are the Wnts, although information about the role of these molecules in fish is very limited. Canonical Wnt signaling regulates the differentiation of myoblasts, while the non-canonical pathway controls the self-renewal of SCs and the growth of muscle fibers (28, 29).
In some teleost fishes, like the gilthead sea bream, and as opposed to mammals, the skeletal system is defined as acellular since the skeleton is generated and maintained by chondrocytes, osteoblasts, and osteoclasts but it lacks osteocytes within the calcified extracellular matrix (ECM) (30, 31). Bone lineage determination and differentiation of osteoblasts from mesenchymal stem cells (MSCs) implies various key regulatory actors including specific transcription factors, mineral availability, and environmental conditions (32). The transcription factor expression pattern that controls osteogenesis was elucidated in gilthead sea bream few years ago and it was demonstrated that there is conservation of transcripts compared to mammals (33). Indeed, Runx2 is crucial for the commitment of the osteoprogenitor cells from MSCs and for the differentiation of osteoblasts; besides it coordinates the expression of other genes, such as osteopontin (op), osteonectin (on) and osteocalcin (ocn), among others, which are key factors regulating the ECM mineralization (34). In addition, it is essential to point out that bone is among the only tissues that has a cell type whose main action is to destroy the own tissue, the osteoclast (35). This function is known as bone remodeling and it is crucial to guarantee mineral homeostasis, having regulatory implications also in other tissues such as the skeletal muscle (36).
Bones represent the attachment location for skeletal muscles and therefore this tissue is also important to cope with threatening conditions. Hence, the musculoskeletal system is a complex interconnected structure (37). Together, the musculoskeletal system enables locomotion and has an essential metabolic role, and the interaction of muscle and bone tissues goes beyond the mechanical and involves a great deal of biochemical communication (38, 39). Both tissues are a relevant source of signaling molecules that can act in an autocrine, paracrine, or endocrine way. It is for these reasons that the functions and processes of each individual tissue are regulated by the other one (40, 41).
In this scenario, the main objective of this work is to perform an in vivo study to characterize skeletal muscle regeneration in gilthead sea bream. For this purpose, a mechanical injury will be produced in the muscle and changes in the transcription of key genes involved in the growth regulation of skeletal muscle, including the different systems and molecules mentioned above will be evaluated. In addition, several bone gene markers will be studied to determine the implication of this close-related tissue regulating myogenesis.
2 Materials and methods
2.1 Experimental design
For the muscle regeneration study, 140 gilthead sea bream (Sparus aurata) juveniles (initial body weight = 15.4 ± 3.5 g; initial length = 8.7 ± 0.6 cm (mean ± SEM)) were obtained from a commercial hatchery (Piscimar, Borriana, Spain) and were placed and adapted to the fish facilities of the Faculty of Biology (University of Barcelona). Fish were randomly distributed in three 200 L seawater tanks (45-47 fish/tank). Each tank had a constant flux of 700 L/h in a seawater semi-closed recirculation system with a weekly renewal of 20-30% and a salinity of 35-37 ‰ at a constant temperature of 23 ± 1°C. In the water tank room, there was a photoperiod of 12 h light/dark. Fish were fed ad libitum 3 times per day with a commercial diet (Perla, Skretting, Burgos, Spain) and were kept in the described conditions for their acclimatization during 2 weeks before the experiment. Once the study started, all the animals continued being fed with the same protocol as during the adaptation period. The study was carried out following the EU recommendations and the procedures established by the Spanish and Catalan governments. The protocol was approved by the Ethics and Animal Care Committee of the University of Barcelona (CEEA 37/20).
This study aimed to better understand adult muscle regeneration after an injury. To do that, 140 gilthead sea bream were divided into two groups: injured fish (I) and control fish (C) and randomly distributed in three tanks. Specifically, in each seawater tank, 33 fish from the injured group were housed with 14 fish from the control group. First, all the gilthead sea bream juveniles were anesthetized with MS222 (100 mg/L) and then measured and weighed. To identify the fish, a passive integrated transponder (PIT) tag (ID-100A (1.25) Nano transponder; Trovan Electronic Identification Systems, Madrid, Spain) was inserted subcutaneously into the left anterior epaxial muscle just below the first radius. Subsequently, the injury was performed with a 2.108 mm (14G) diameter needle inserted vertically into the left epaxial muscle below the sixth radius to a depth of 1 cm. To know exactly in the future days where the needle was introduced, the tip of the sixth radius was cut, also to the control fish for comparative purposes. Then, the wound was healed with iodine alcoholic solution and the fish were allowed to recover in a separated small tank before being returned to the 200 L tank.
Tissue samples were obtained at days 0, 1, 2, 4, 8, 16 and 30 after the injury. At each time, 16 fish were randomly selected for the study groups. To carry out the samplings, fish were first anesthetized, identified reading the pit tag, weighed to note the changes on body weight and then, blood was drawn to ensure a cruelty-free death. In injured fish (I), a section of the epaxial muscle was removed from the left loin (injured muscle) and the right loin as a self-control for each fish. The tissue size of the muscle extracted was 0.5 cm wide and 1 cm long just below the sixth cut radius. The spinal column was also excised counting as the bone sample of a fish with injury. Moreover, a sample of the vertebral column was also extracted from the control fish (C) to obtain bone control samples. All tissue samples were placed in RNase-free Eppendorf tubes that were stored in liquid nitrogen during sampling and then at -80 °C until the performance of the gene expression analyses.
2.2 Histology analysis
Epaxial muscle slices were fitted in histological cassettes and were dehydrated in graded ethanol series and embedded in paraffin. Sections of 7-10 µm were cut with a microtome (pfm, ROTARY 3003, Köln, Germany) and stained with Hematoxylin and Eosin and Sirius Red for collagen staining (42). All preparations were observed under a light microscope and photographed (Olympus PM10SP Automatic Photomicrography System). All reagents for histology staining were purchased from Merck (Merck, Mollet del Vallès, Spain).
2.3 Gene expression
2.3.1 RNA extraction and cDNA synthesis
Total RNA was extracted from 100 mg of epaxial muscle or vertebral bone with 1 mL of TRI Reagent Solution® (Applied Biosystems, Alcobendas, Spain) and homogenized on a Precellys® Evolution coupled to a Cryolys cooling unit (Bertin Instruments, Montigny-le-Bretonneux, France). The RNA extraction was performed following the manufacturer’s instructions of the TRI Reagent Solution®. The final concentration of each sample was obtained using the Nanodrop 2000 (Thermo Scientific, Alcobendas, Spain) and RNA integrity was confirmed in a 1% agarose gel (w/v) stained with SYBR-Safe DNA Gel Stain® (Life Technologies, Alcobendas, Spain).
For RNA-complementary DNA (cDNA) synthesis, 2 μg of total RNA were treated with DNase I Amplification Grade® (Life Technologies, Alcobendas, Spain) to remove all genomic DNA. Reverse transcription was carried out with the First Strand cDNA synthesis Transcriptor Kit ® (Roche, Sant Cugat del Valles, Spain) following the manufacturer’s protocol.
2.3.2 Quantitative real time PCR
The mRNA transcript levels of the genes for both tissues were analyzed by real-time quantitative PCR (RT-qPCR) using the CFX384™ Real-Time System (Bio-Rad, El Prat de Llobregat, Spain). The primers used for each tissue are listed in Table S1 and the specificity of the amplification, and the absence of non-specific primer bindings were tested in silico. First, a dilution curve with a pool of samples was run to confirm primer efficiency and to determine the appropriate cDNA dilution for each assay. All the analyses were performed in triplicate wells using 384-well plates with 2.5 µL of iTaq Universal SYBR Green Supermix (Bio-Rad, El Prat de Llobregat, Spain), 0.25 μL of forward (250 nM) plus reverse (250 nM) primers, 1.25 μL of DEPC water and 1 µL of diluted cDNA for each sample, in a final volume of 5 µL following the conditions described by Salmerón et al. (16). Moreover, three negative controls were also included and ran in duplicate: no template control (NTC), no reverse transcriptase control (RTC), and PCR control (PCR, MilliQ water). The expression level of each gene was calculated with the ΔΔCq method considering the efficiency of each primer pair (43) and was analyzed relative to the geometric mean of the reference genes ribosomal protein s18 (rps18), ribosomal protein l27a (rpl27) and elongation factor 1 alpha (ef1a). The reference genes, the most stable under different conditions, were confirmed with the geNorm algorithm. Both, reference genes stability and relative expression of the target genes were determined using the Bio-Rad CFX Manager Software v. 3.1 (Hercules, CA, USA). Finally, to specifically determine that the changes in relative gene expression were due to the injury, the expression results of the injured samples (i.e., left loin and column of injured fish) were normalized by the gene expression values of the corresponding control samples (i.e., right loin of injured fish and column of control fish).
2.4 Statistical analyses
Data were analyzed using IBM® SPSS® Statistics v.25 (IBM, Armonk, NY, USA) and were presented as mean ± standard error of the media (SEM). Normal distribution was analyzed using the Shapiro-Wilk test and homogeneity of the variances (homoscedasticity) was assessed with Levene’s test. Differences were tested by one-way analysis of variance (ANOVA) and the post-hoc Tukey HSD. If necessary, the nonparametric Kruskal Wallis test and the post-hoc Games-Howell were used. Additionally, one-way ANOVA was performed to verify that the tank did not influence the measured parameters. Statistical differences were considered significant when p < 0.05. Data were plotted using GraphPad Prism® v. 7 (GraphPad Software, La Jolla, CA, USA, www.graphpad.com).
3 Results
3.1 Histological characterization during muscle regeneration
The qualitative analysis of the injured muscle tissue samples showed that on day 1 it was possible to observe the presence of dead fibers beginning the process of necrosis in the area where the injury occurred (Figures 1A, B). Additionally, an increase in collagen was observed in the myoseptum of the injured area (Figure 1E). Later, at day 16, it was no longer possible to find any trace of injured fibers and instead new healthy fibers were observed, in addition to a reduction in the amount of collagen in the myoseptum (Figures 1C, F). Finally, in the muscle samples at day 30 neither there were signs of injury nor increased collagen in the myoseptum. Furthermore, it could be seen how the new fibers formed in the injured area were more similar in size to the fibers of unaffected regions showing the entry of the cells into the maturation stage (Figures 1D, G).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g001.jpg)
Figure 1 Histological evaluation of the muscle regeneration process at 1, 16 and 30 days post-injury. Muscle sections (7-10 µm) were stained with hematoxylin/eosin (A-D) and Sirius red (E-G). The arrowheads mark the presence of muscle fibers in the process of necrosis and the asterisks mark the myoseptum.
3.2 Expression of genes from the Gh-Igfs axis during muscle regeneration
Total Igf-1 mRNA levels (igf-1abc) significantly decreased from day 0 to day 1 and, subsequently, increased from day 4 post-injury to recover normal levels that were maintained at day 30 (Figure 2A). Concerning the Igf-1 splice variants, igf-1b mRNA levels showed a significant increase between days 1 and 2 that was recovered at day 8, while the other two variants (igf-1a and igf-1c) did not show any clear pattern of expression or significant differences (Figures 2B–D). Contrarily, the igf-2 gene was significantly upregulated from day 0 to day 1 and after that, the expression decreased by day 2 and specially at day 16 remaining low until day 30 (Figure 2E). Growth hormone receptors (Ghrs) reported a significant fall in the case of ghr-1 expression at day 1 and a notable subsequent recovery 16 days after injury, although it was not possible to determine any significant expression change in the paralogue ghr-2 (Figures 2F, G). With regards to downstream signaling molecules, akt2 gene expression remained constant through the experiment and the mRNA levels of the serine/threonine protein kinase tor showed a non-significant tendency to decrease between day 0 and day 1 (Figures 2H, I).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g002.jpg)
Figure 2 Relative gene expression of the insulin-like growth factors, igf-1abc (A), igf-1a (B), igf-1b (C), igf-1c (D), igf-2 (E), growth hormone receptors, ghr-1 (F), ghr-2 (G), and signaling molecules, akt2 (H), and tor (I) in injured epaxial muscle of fish. Data are shown as means ± SEM, n=10. Statistical analysis was assessed by one-way ANOVA. Different letters indicate significant differences of factor time (days post-injury) (Tukey’s HSD or Games-Howell post-hoc test, p < 0.05).
3.3 Proteolytic systems genes expression during muscle regeneration
Calpain-1 (capn1) and its regulatory subunits (capns1a and capns1b) showed basal mRNA levels without changes throughout the different days of study (Figures 3A–C). On the contrary, capn2 showed a significant increase in mRNA levels at day 1 then decreasing until day 8 (Figure 3D). Then, capn3 was downregulated from day 0 to day 4 after injury then recovering the initial expression levels by day 8 (Figure 3E). Cathepsin L (ctsl) presented a significant transient increase of mRNA levels at day 2, but no changes in cstda gene expression were observed (Figures 3F, G). Regarding the ubiquitin-proteasome system members, mafbx significantly increased from day 4 to day 8 to later decrease by day 16, whereas murf1 showed significantly higher mRNA levels at day 30 compared to days 1 and 4, and the mRNA levels of n3 were not changed in this study (Figures 3H–J).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g003.jpg)
Figure 3 Relative gene expression of the proteolytic systems' markers, capn1 (A), capns1a (B), capns1b (C), capn2 (D), capn3 (E), ctsl (F), ctsda (G), mafbx (H), murf1 (I), and n3 (J) in injured epaxial muscle of fish. Data are shown as means ± SEM, n=10. Statistical analysis was assessed by one-way ANOVA. Different letters indicate significant differences of factor time (days post-injury) (Tukey’s HSD or Games-Howell post-hoc test, p < 0.05).
3.4 Gene expression of other local markers during muscle regeneration
Pax7 and the hepatocyte growth factor receptor (cmet) mRNA levels were low until day 8 to then increase significantly at day 16 maintaining the expression elevated until day 30 (Figures 4A, B). The pcna expression did not present significant changes (Figure 4C). The caveolin gene, cav1 was maintained in a basal expression showing only a significant upregulation at day 16 compared to days 1 and 2 (Figure 4D). Moreover, cav3 showed a similar profile but in this case, the gene expression strongly increased from day 8 to 16 and continued high until 30 days after injury (Figure 4E). Concerning the negative regulators of muscle growth, mstn1 mRNA levels increased progressively from day 4 until day 30, with significant differences at day 16 post-injury (Figure 4F). Similarly, mstn2 mRNA levels showed an upregulation from day 8 to 16 and continued high until the end of the study (Figure 4G). The secreted factor wnt5b, from the Wnt family, increased gradually its mRNA levels from day 4 until day 30, being significant that rise at 16 days onwards (Figure 4H). In the case of the small leucine-rich proteoglycan osteoglycin, muscle mRNA levels of ogn1 increased significantly from day 8 to 30 and similarly ogn2 increased from day 8 to day 16 and then remained high until day 30 (Figures 4I, J). Finally, for the vasoendothelial growth factor (vegfa), minor fluctuations in its mRNA levels were observed, with only a significant increase in expression when comparing between days 4 and 16 (Figure 4K).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g004.jpg)
Figure 4 Relative gene expression of other muscle-related markers involved in muscle regeneration, pax7 (A), cmet (B), pcna (C), cav1 (D), cav3 (E), mstn1 (F), mstn2 (G), wnt5b (H), ogn1 (I), ogn2 (J), and vegfa (K) in injured epaxial muscle of fish. Data are shown as means ± SEM, n=10. Statistical analysis was assessed by one-way ANOVA. Different letters indicate significant differences of factor time (days post-injury) (Tukey’s HSD or Games-Howell post-hoc test, p < 0.05).
3.5 Inflammatory markers gene expression during muscle regeneration
The gene expression of the pro-inflammatory cytokine interleukin 6 (il-6) was incremented significantly at day 1 post-injury and then decreased progressively up to day 8 when the expression returned to significantly lower basal levels (Figure 5A). The il-15 expression was similarly slightly increased at day 1 being those levels significantly higher comparing with days 16 and 30 (Figure 5B). The colony-stimulating factor 1 receptor (csf1r) gene expression raised significantly from day 2 to 4 and remained high until day 30 (Figure 5C). Moreover, tnfα showed a significant upregulation at day 16 and continued high until the end of the study (Figure 5D) and, the il-1β expression did not show any significant change (Figure 5E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g005.jpg)
Figure 5 Relative gene expression of the inflammatory markers, il-6 (A), il-15 (B), csf1r (C), tnfα (D), il-1β (E) in injured epaxial muscle of fish. Data are shown as means ± SEM, n=10. Statistical analysis was assessed by one-way ANOVA. Different letters indicate significant differences of factor time (days post-injury) (Tukey’s HSD or Games-Howell post-hoc test, p < 0.05).
3.6 Bone markers gene expression during muscle regeneration
The bone morphogenetic protein 2 (bmp2) presented a significant mRNA levels reduction from day 1 to day 8 post-injury, followed by an important recovery by day 16 (Figure 6A). Concerning the bone mRNA levels of osteoglycin, ogn1 increased significantly from day 0 to day 4, then its levels decreased in a major way until day 16 and finally, it raised importantly again at day 30 (Figure 6B). The ogn2, showed a parallel but delayed pattern of expression showing its significantly highest mRNA levels at day 8 (Figure 6C). The mRNA levels of the key mineralization marker osteocalcin (ocn) slightly increased from day 2 to 4 and then fell significantly at day 8 to finally being recovered at day 30 (Figure 6D). The analysis of the transcription factor runx2 expression showed a significant increment in the mRNA levels between day 4 and day 8, to then slightly decrease again although not significantly (Figure 6E). Osteonectin (on) mRNA levels increased progressively in a long-term way from day 0 to day 30 presenting significant differences only between those two, initial and final days (Figure 6F). Regarding myostatins in bone, mstn1 showed no expression and mstn2 mRNA levels decreased significantly from day 0 to 4 and then slightly recovered its expression at day 16 (Figure 6G). Finally, the mRNA levels of the proteolytic marker cathepsin k (ctsk) decreased significantly from day 2 to day 8 after the muscle injury being maintained afterwards (Figure 6H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g006.jpg)
Figure 6 Relative gene expression of the bone markers, bmp2 (A), ogn1 (B), ogn2 (C), ocn (D), runx2 (E), on (F), mstn2 (G), and ctsk (H) in the vertebra of injured of fish. Data are shown as means ± SEM, n=10. Statistical analysis was assessed by one-way ANOVA. Different letters indicate significant differences of factor time (days post-injury).
4 Discussion
Juvenile gilthead sea breams were injured in the epaxial white muscle with a 14G needle. Then, muscle histological analysis plus assessment of changes in the transcription of key genes from muscle and bone tissues during a 30-day regeneration period was performed. The histological evaluation revealed at day 1 fibers in the process of necrosis, while 16 days after the injury no traces of damaged or necrotic fibers were observed in the gilthead sea bream. Our results agree with those of Rowlerson et al. (44), who found in the same species that muscle regeneration was weakly detected seven days post-injury, still with some necrotizing fibers, but at 21 days post-lesion many areas were regenerated, and the presence of new young fibers could be observed. In fact, in our study, at day 30 the myotomes showed a complete normal morphology with many new and small muscle fibers present.
The Gh-Igfs axis is the main regulator of muscle growth and contributes to the process of regeneration after tissue injury in vertebrates (20, 45) including fish (46). In this sense, the contribution of Igf-1 during tissue repair has been recognized (47) but the role of Igf-2 is not so well known although its contribution to myogenesis was suggested by demonstrating binding to the Igf-1 receptor, activating the Akt pathway and then MyoD (48, 49). Furthermore, Ge et al. (50) found in mice that igf-2 mRNA levels also increased drastically during early muscle regeneration. In this study, the peak of igf-2 at the beginning of muscle repair (i.e., day 1) is consistent with this early and pre-Igf-1 function at the onset of myogenesis. In fact, it has been previously observed in gilthead sea bream myocytes that Igf-2 contributes to muscle regulation mainly in the early stages of development and before Igf-1 (9, 51). Indeed, Igf-1 action takes place slightly later towards the initial stages of myoblast differentiation (52). In agreement with that, in the present study, total igf-1 expression decreased at day 1, recovering the expression levels at days 4, 8, 16 and 30 after injury, probably mainly due to the upregulation of the splice variant igf-1b at day 2 post-injury. Concerning receptors, ghr-1 expression followed a similar profile to that of igf-1 with a rapid decrease the first days after injury and a progressive and slow increase afterwards, returning to initial levels at day 16 when the animals recovered from the injury, a response that agrees with its anabolic role in this species (8, 53).
The rebuilding of skeletal muscle after an injury requires a fine regulation of proteinases expression. In that sense, calpains, cathepsins and the UbP are the major systems for muscle protein degradation, and they may contribute to the remodeling of skeletal muscle during regeneration (54). Previously, we have found important transcriptional responses of these proteolytic genes in gilthead sea bream subjected to sustained exercise, different diets, or the combination of both factors (15, 55, 56) or in response to 21 days of fasting (57). In this study, we found at least one gene expressed from each one of the three proteolytic systems at different times during regeneration. In fact, the expression of capn2 and ctsl was increased at day 1 and 2 respectively, while mafbx increased at day 8 post-injury and murf1 at day 30 compared to the early days. In a study done in rats that were injected with 0.75% bupivacaine, a drug that at this concentration produces a chemical injury, the enzymatic activity of Capn1 and Capn2 was maximum at day 5 post-injury, that of Cathepsin L and D at day 3, and the proteasome peaked at day 5 to progressively return to control levels by day 21 (54). This profile of proteolytic gene expression can be comparable to the expression profiles obtained in this study in gilthead sea bream and suggests that the sequence of activation of proteolytic genes during muscle regeneration appears to be quite well conserved through vertebrates. Specifically, calpains seemed to have a more relevant role in proliferation of gilthead sea bream cultured myocytes and the UbP system in cells differentiation (58). In addition, it has been demonstrated that mafbx interacts with myod regulating its expression; therefore, in this study mafbx might be inhibiting myod expression prematurely at day 8 in order to ensure the correct activity of Mrfs from day 16 onwards (59). In summary, some of these proteolytic genes could be necessary in our study species, to modulate myogenesis, as observed previously in mammals (60). Thus, these genes could be involved in myoblast proliferation and differentiation (61, 62) but are also required for myoblast fusion (63, 64).
In line with this, as a part of the characterization of the recently identified fusogens Myomaker and Myomixer in gilthead sea bream, aiming to estimate their functions, gene expression of these molecules in addition to the Mrfs were studied under different myogenesis conditions including this model of muscle regeneration (65). All these molecules showed low expression values until day 8, and then started to increase to reach a peak at day 16 post-wound to subsequently slowly decrease to recover initial levels by day 30. Recently, Manneken et al. (66) reviewed the process of muscle regeneration in zebrafish (Danio rerio), focusing on the differences to amniotes myogenesis during both, development, and repair. It is noticeable the similarity of the transcriptional profile of the different Mrfs analyzed in that study with those in gilthead sea bream. Moreover, Rowlerson et al. (44) showed that at 7 days post-lesion a weak reaction to desmin antibody appeared only in bigger fibers but not in presumed newly formed fibers; and at 21 days the regeneration was well advanced in many muscle areas. Forcina et al. (67) and Musarò (14) also reviewed the mechanisms regulating muscle regeneration in mammals, indicating that the MRFs participate in the process between 4 and 7 days post-injury, reaching the maturation and functional recovery of the muscle between days 10 and 15. The profile is again very similar to that observed in our previous study in gilthead sea bream (65) although the faster myogenic regeneration process in mammals could be explained by their higher body temperature compared to fish. In fact, in rainbow trout (Oncorhynchus mykiss), a species that grows at a lower water temperature than gilthead sea bream, the complete muscle regeneration in a similar experiment took longer (i.e., day 30) and myomaker and myomixer indeed peaked later, between days 16 and 30 after injury (68).
Other regulatory factors of muscle growth like pax7, cmet, caveolins, wnt5b, ogn, vegfa or myostatins followed the same tendency, with non-significant small fluctuations up to day 8 but significant increases in expression at 16 and/or 30 days post-lesion. Indeed, the parallelism between the expression profiles of these genes and the Mrfs makes sense due to its recognized role during mammalian myogenesis. For instance, very little is known about the role of caveolins in fish, but the increase observed in cav1 and cav3 suggests their participation in fish myogenesis activating SCs as described in mice and humans (26, 69); following pax7 and cmet, factors that are known to stimulate cell cycle activation and proliferation of SCs (21, 22). Moreover, the progressive increase in the expression of wnt5b was indicative of its potential role at the end of muscle regeneration during myoblast differentiation and fusion although a larger contribution from this or other members of the Wnt family could be expected in the early stages (28). However, again, the information in fish is very scarce and probably their role is not exactly the same as in mammals (29). The expression of ogn1 and ogn2 in the muscle increased from day 8. In a previous study, Costa et al. (70) determined the expression of ogn in primary cultures of gilthead sea bream myocytes and observed a significant increase in their expression at day 8 of culture development, at the stage of differentiation and fusion of myocytes to form myotubes. In the present study, it is possible to consider that ogn transcripts in muscle are also participating in the fusion of the newly formed fibers. Regarding the expression of vegfa, the peak observed at 16 days post-injury is consistent with the main period of new muscle fibers formation plus maturation, and therefore, increased angiogenesis and myogenic factors highest expression; supported by the role previously demonstrated for this factor in muscle remodeling in gilthead sea bream (56). Finally, in the case of mstn1 and 2, their gene expression was increased by day 30, which fits with their known role in gilthead sea bream as negative regulators of myogenesis, therefore helping to the completion of the process in this regeneration scenario (71). Thus, the period between days 8 and 16 after the injury seems to be in this species the main period of activation of myogenic factors to induce muscle regeneration. This is a little later that the timing described in mammals but faster than that observed in rainbow trout (67, 68, 72, 73).
In the healing process of a tissue, acute inflammation and immune cells play a critical role in almost all phases of regeneration. Damaged muscle fibers undergoing necrosis induce the infiltration of many types of immune cells (74, 75). Therefore, the infiltrating cells help clearing dead cells from the injured area as well as secrete different types of cytokines to recruit more immune cells in addition to producing cellular responses to regulate muscle cell activation, proliferation, and differentiation (19). In our study, it has been possible to observe a first wave of il-6 and il-15 increased expression in the early days after injury. In a mouse model of cardiotoxin-induced muscle injury, Zhang et al. (76) showed how high levels of il-6 were detected from the first day after injury (24 h) and demonstrated that this cytokine is crucial for macrophage infiltration, and that the lack of it in a knockout mouse model inhibited myoblast proliferation. Moreover, some studies support the fact that Il-6 also induces proteolysis of muscle cells by activating and modulating proteolytic systems (77, 78). On the other hand, il-15 is highly expressed in mammalian muscle, and it is believed that, in addition to its inflammatory function in a pro-inflammatory context, it may be a myokine as it promotes myogenesis in response to exercise as well as facilitates muscle regeneration (79, 80). The results of our study could be showing such mentioned functions, and in the regeneration in gilthead sea bream both Il-6 and Il-15 could be acting in a pro-inflammatory phase in addition to favoring myogenesis, although we cannot distinguish if these cytokines are specifically expressed in muscle cells, immune cells, or both. Furthermore, Tnfα has been reported to play an important role in muscle regeneration in mammals (81, 82). This cytokine can attract muscle SCs to the damaged site to promote their proliferation and differentiation (83). The increased expression of tnfα in gilthead sea bream muscle 16 days after injury could be indicative of this function since its upregulation occurred at the same time as the entire myogenic program was initiated. In fact, the releasing of a whole continuum of pro-inflammatory cytokines into the injured muscle stimulates the infiltration of monocytes that go on to exhibit an inflammatory profile to carry out phagocytosis and then, switch to anti-inflammatory macrophages to stimulate myogenesis (84–86). The colony stimulating factor 1 (Csf1) controls the production, differentiation, and function of macrophages and its receptor (Csf1R) mediates the biological effects (87). In our study, the strong increase detected in the expression of csf1r from day 4 on could indicate the arrival to the muscle of macrophages after the wave of pro-inflammatory interleukins. Moreover, the expression of this receptor remains elevated up to day 30, which could mean the maintenance of a population of anti-inflammatory and myogenesis-supporting macrophages until the end of the regenerative process.
Referring to genes expressed in bone, the small leucine-rich proteoglycans ogn1 and ogn2 have been characterized in gilthead sea bream during the in vitro development of different cell types including bone and muscle (70). Hence, their significantly increased expression in the bone tissue close to the muscle wound at days 4 and 8 post-injury, respectively, in the current study, agrees with their reported upregulation at the onset of osteoblasts as well as myocytes differentiation. Moreover, Tanaka et al. (88) also demonstrated the regulatory crosstalk role between bone and muscle of ogn, but references in fish are very poor. In this study, the earlier expression of ogn in bone and later in muscle could indicate a stimulus from bone to enhance the recovery of the injured muscle. Similarly, ocn expression modulation during regeneration with increased levels at day 4 versus day 8, support its role as a promising osteokine participating in the regulation of the muscle response after an injury in fish, as already demonstrated in mammals (89). Indeed, Ocn has been proven to be able to increase importantly muscle activity in elderly mice individuals (90) and to directly promote protein synthesis in myotubes, explaining why this hormone is responsible for muscle maintenance during aging (91). Thus, the increase of ocn at day 4 post-injury before the upregulation of Mrfs would support its stimulatory role in fish muscle too. The bmp2 and ctsk genes showed an early activation of expression at days 1 and 2 post-injury, respectively, although bmp2 expression decreased at day 8 recovering basal levels at the end of the experiment. This early upregulation of bmp2 could be related to the highest expression found for this molecule at the beginning of embryonic development in gilthead sea bream (92), but also to situations in which bone (and muscle) growth is stimulated, such as for example in Atlantic salmon (Salmo salar) reared at elevated water temperatures (93). Moreover, the peak of expression of ctsk at day 2 resembled the increase in expression observed after 20 days of fasting (57) suggesting that ctsk and therefore, bone resorption, may have a role in muscle remodeling after stressful conditions. The key osteogenic factor runx2 after a fall at day 4, peaked at day 8, keeping similar levels until the end of the experiment, perhaps suggesting an induction of osteogenesis at this time to support the subsequent muscle mass development. Furthermore, the inter-tissue inhibitory role of Myostatin has been already well characterized in mammals (89, 94). In our model species, a previous study of fasting and refeeding proved how myostatin could act coordinating musculoskeletal growth and suggested its potential role as osteokine in fish (57). In this regeneration study, a regulatory role for myostatin could also be hypothesized, since bone decreased mstn2 expression by day 4 in response to the muscle injury possibly to allow proper muscle remodeling, thus giving more insights of the importance of myostatin between the bone and muscle crosstalk in fish. Finally, the ECM marker on showed a progressive increase during regeneration with the maximum expression reached at day 30 coinciding with the highest expression observed for on at the end of in vitro and in vivo osteogenesis in gilthead sea bream according to its key role regulating the later stages of tissue development (33, 92). In summary, these results open a promising new line of research to study the possible crosstalk between bone and muscle in fish to regulate skeletal muscle development and growth.
5 Conclusions
In gilthead sea bream, the generation of a wound and the presence of damaged and dead muscle fibers firstly activated pro-inflammatory genes as il-6 and il-15 that contributed cleaning the injured tissue and stimulating the muscle regeneration process (Figure 7). Then, proteolytic factors, first capn2 and ctsl, and later on the UbP system member mafbx allowed an effective degradation of the damaged tissue.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/1101356/fendo-14-1101356-HTML/image_m/fendo-14-1101356-g007.jpg)
Figure 7 Graphical summary of the proposed model for fish muscle regeneration after an injury. The figure shows an overview of the different processes occurring after an injury in the gilthead sea bream muscle in a circular way, presenting the progressive appearance of the main molecules transcriptionally upregulated/modulated in both, muscle and bone. The circle is divided into the four main steps that take place during the regeneration process, from the initial inflammation and proteolysis to the satellite cells activation and myogenesis. Note that Mrfs, myomaker and myomixer data has been extracted from (65). The hypothetical muscle and bone interactions suggested are indicated by crossing arrows and question mark symbols. Created with BioRender.com.
At the very early stage of 1-day post-injury, it was noticeable the increase of igf-2 and the decreases of ghr1 and total igf-1 expression. The return to basal levels at day 8 of igf-2 was parallel to the increase of igf-1 and simultaneous to the induction of many other genes’ expression involved in the regulation of myogenesis like pax7, cmet, cav1, cav3, and wnt5b, indicating the beginning of muscle remodeling. In fact, vegfa also showed a clear peak of expression at 16 days post-lesion, which parallels the Mrfs and the fusogens myomaker and myomixer patterns of expression previously described (65), thus suggesting that progression of myogenesis and angiogenesis is taking place resulting in new fibers formation.
Furthermore, it could be hypothesized that some bone-derived genes such as bmp2, ogn1, ocn and mstn2 could be participating by sending information to the muscle mainly at the early stages of the regenerative process to control harmonic musculoskeletal growth, thus potentially presenting a role as osteokines in fish.
Overall, this study shows how different mechanisms are orchestrated to recover an injured muscle and allows knowing for the first time in gilthead sea bream the timely contribution of different molecules along the regeneration process of damaged skeletal muscle fibers.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Materials. Further inquiries can be directed to the corresponding author.
Ethics statement
The study was conducted according to the guidelines of the European Union Council directive (EU 2010/63) and approved by the Ethics Committee of the University of Barcelona (protocol code CEEA 37/20).
Author contributions
JG conceptualized the study. AO-T, MP-A, VJ-P, IG-P, AS-M and JG performed the sampling. AO-T, MP-A, VJ-P, FM and IG-P performed the laboratory analyses. AO-T, MP-A, VJ-P and JG analyzed and interpreted the data. IN, EC, JB and JG acquired funding. AO-T, MP-A, IN, EC, JF-B, DG, JB and JG drafted and critically reviewed the manuscript. All authors read and approved the final paper. The authors have declared no conflict of interests. All authors contributed to the article and approved the submitted version.
Funding
This publication is part of the R+D+i projects AGL2015-70679-R and RTI2018-100757-B-I00 to JG and JB, and AGL2017-89436-R and PID2020-116172RB-I00 to IN and EC funded by the Spanish “Ministerio de Ciencia e Innovación” (MCIN/AEI/10.13039/501100011033/) and by “ERDF A way of making Europe”, by the “European Union”, and the “Xarxa de Refèrencia d’R+D+I en Aqüicultura” and the 2017-SGR1574 from the “Generalitat de Catalunya”. MP-A and IG-P were supported by predoctoral fellowships (BES-2016-078697 and PRE2019-089578 respectively) from the MCIN/AEI/10.13039/501100011033 and by “ESF Investing in your future”. FM was supported by a grant from the Department of Scholarship and Students’ Affairs Abroad from the Ministry of Science, Research and Technology of I.R. Iran.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2023.1101356/full#supplementary-material
References
1. FAO. El Estado mundial de la pesca y la acuicultura 2020. la sostenibilidad en acción. Roma: Marine Pollution Bulletin (2020). doi: 10.4060/ca9229es
2. Talbot C. Some aspects of the biology of feeding and growth in fish. Proc Nutr Soc (1993) 52:403–16. doi: 10.1079/PNS19930081
3. Johnston IA. Muscle development and growth: Potential implications for flesh quality in fish. Aquaculture (1999) 177(1–4):99–115. doi: 10.1016/S0044-8486(99)00072-1
4. Mommsen TP. Paradigms of growth in fish. Comp Biochem Physiol - B Biochem Mol Biol (2001) 129(2–3):207–19. doi: 10.1016/S1096-4959(01)00312-8
5. Mommsen TP, Moon TW. Hormonal regulation of muscle growth. Fish. Physiol (2001) 18(C):251–308. doi: 10.1016/S1546-5098(01)18009-X
6. Fuentes EN, Valdés JA, Molina A, Björnsson BT. Regulation of skeletal muscle growth in fish by the growth hormone - insulin-like growth factor system. Gen Comp Endocrinol (2013) 192(June):136–48. doi: 10.1016/j.ygcen.2013.06.009
7. Reinecke M, Björnsson BT, Dickhoff WW, McCormick SD, Navarro I, Power DM, et al. Growth hormone and insulin-like growth factors in fish: Where we are and where to go. Gen Comp Endocrinol (2005) 142(1-2 SPEC. ISS.):20–4. doi: 10.1016/j.ygcen.2005.01.016
8. Vélez EJ, Lutfi E, Azizi S, Montserrat N, Riera-Codina M, Capilla E, et al. Contribution of in vitro myocytes studies to understanding fish muscle physiology. Comp Biochem Physiol Part - B: Biochem Mol Biol (2016) 199:67–73. doi: 10.1016/j.cbpb.2015.12.003
9. Rius-Francino M, Acerete L, Jiménez-Amilburu V, Capilla E, Navarro I, Gutiérrez J. Differential effects on proliferation of GH and IGFs in sea bream (Sparus aurata) cultured myocytes. Gen Comp Endocrinol (2011) 172(1):44–49(1). doi: 10.1016/j.ygcen.2011.03.024
10. Tiago DM, Laizé V, Cancela ML. Alternatively spliced transcripts of sparus aurata insulin-like growth factor 1 are differentially expressed in adult tissues and during early development. Gen Comp Endocrinol (2008) 157(2):107–15. doi: 10.1016/j.ygcen.2008.04.006
11. Duan C, Ren H, Gao S. Insulin-like growth factors (IGFs), IGF receptors, and IGF-binding proteins: Roles in skeletal muscle growth and differentiation. Gen Comp Endocrinol (2010) 167(3):344–51. doi: 10.1016/j.ygcen.2010.04.009
12. Froehlich JM, Galt NJ, Charging MJ, Meyer BM, Biga PR. In vitro indeterminate teleost myogenesis appears to be dependent on Pax3. In Vitro Cell Dev Biol Anim (2013) 49(5):371–85. doi: 10.1007/s11626-013-9616-2
13. Bentzinger CF, Wang YX, Rudnicki MA. Building muscle: molecular regulation of myogenesis. Cold Spring Harbor Perspect Biol (2012) 4(2):a008342. doi: 10.1101/cshperspect.a008342
14. Baldwin KM, Haddad F. Skeletal muscle plasticity. Am J Phys Med Rehabil (2002) 81(Supplement):S40–51. doi: 10.1097/00002060-200211001-00006
15. Salmerón C, García de la serrana D, Jiménez-Amilburu V, Fontanillas R, Navarro I, Johnston IA, et al. Characterisation and expression of calpain family members in relation to nutritional status, diet composition and flesh texture in gilthead Sea bream (Sparus aurata). PLos One (2013) 8(9):e75349. doi: 10.1371/journal.pone.0075349
16. Salmerón C, Navarro I, Johnston IA, Gutiérrez J, Capilla E. Characterisation and expression analysis of cathepsins and ubiquitin-proteasome genes in gilthead sea bream (Sparus aurata) skeletal muscle. BMC Res Notes (2015) 8(1):1–15. doi: 10.1186/s13104-015-1121-0
17. Musarò A. The basis of muscle regeneration. Adv Biol (2014) 2014(Table 1):1–16. doi: 10.1155/2014/612471
18. Panci G, Chazaud B. Inflammation during post-injury skeletal muscle regeneration. Semin Cell Dev Biol (2021) 119(February):32–8. doi: 10.1016/j.semcdb.2021.05.031
19. Yang W, Hu P. Skeletal muscle regeneration is modulated by inflammation. In J Orthopaedic Translation (2018) 13:25–32. doi: 10.1016/j.jot.2018.01.002
20. Ciciliot S, Schiaffino S. Regeneration of mammalian skeletal muscle: Basic mechanisms and clinical implications. Curr Pharm Design (2010) 16(8):906–14. doi: 10.2174/138161210790883453
21. Johnson SE, Allen RE. Proliferating cell nuclear antigen (PCNA) is expressed in activated rat skeletal muscle satellite cells. J Cell Physiol (1993) 154(1):39–43. doi: 10.1002/jcp.1041540106
22. Anderson JE. Hepatocyte growth factor and satellite cell activation - growth factors and cytokines in skeletal muscle development, growth, regeneration and disease. In: White J, Smythe G, editors. Growth factors and cytokines in skeletal muscle development, growth, regeneration and disease. (Switzerland: Springer International Publishing) (2016). p. 1–25. doi: 10.1007/978-3-319-27511-6_1
23. Gal-Levi R, Leshem Y, Aoki S, Nakamura T, Halevy O. Hepatocyte growth factor plays a dual role in regulating skeletal muscle satellite cell proliferation and differentiation. Biochim Biophys Acta - Mol Cell Res (1998) 1402(1):39–51. doi: 10.1016/S0167-4889(97)00124-9
24. Millay DP. Regulation of the myoblast fusion reaction for muscle development, regeneration, and adaptations. Exp Cell Res (2022) 415(2):113134. doi: 10.1016/j.yexcr.2022.113134
25. Thomas CM, Smart EJ. Caveolae structure and function: Caveolae review series. J Cell Mol Med (2008) 12(3):796–809. doi: 10.1111/j.1582-4934.2008.00295.x
26. Volonte D, Liu Y, Galbiati F. The modulation of caveolin‐1 expression controls satellite cell activation during muscle repair. FASEB J (2005) 19(2):1–36. doi: 10.1096/fj.04-2215fje
27. Shang L, Chen T, Deng Y, Huang Y, Huang Y, Xian J, et al. Caveolin-3 promotes glycometabolism, growth and proliferation in muscle cells. PLos One (2017) 12(12):e0189004. doi: 10.1371/journal.pone.0189004
28. Girardi F, le Grand F. Wnt signaling in skeletal muscle development and regeneration. In: Progress in molecular biology and translational science. (United States: Elsevier Inc) (2018) 153. doi: 10.1016/bs.pmbts.2017.11.026
29. von Maltzahn J, Chang NC, Bentzinger CF, Rudnick MA. Wnt signaling in myogenesis. Trends Cell Biol (2012) 22(11):602–9. doi: 10.1016/j.tcb.2012.07.008
30. Shahar R, Dean MN. The enigmas of bone without osteocytes. BoneKEy Rep (2013) 2(March):1–9. doi: 10.1038/bonekey.2013.77
31. Witten PE, Huysseune A. A comparative view on mechanisms and functions of skeletal remodelling in teleost fish, with special emphasis on osteoclasts and their function. Biol Rev (2009) 84(2):315–46. doi: 10.1111/j.1469-185X.2009.00077.x
32. Ytteborg E, Torgersen J, Baeverfjord G, Takle H. Four stages characterizing vertebral fusions in Atlantic salmon. J Appl Ichthyology (2012) 28(3):453–9. doi: 10.1111/j.1439-0426.2012.01984.x
33. Riera-Heredia N, Martins R, Mateus AP, Costa RA, Gisbert E, Navarro I, et al. Temperature responsiveness of gilthead sea bream bone; an in vitro and in vivo approach. Sci Rep (2018) 8(1):1–14. doi: 10.1038/s41598-018-29570-9
34. Vieira FA, Thorne MAS, Stueber K, Darias M, Reinhardt R, Clark MS, et al. Comparative analysis of a teleost skeleton transcriptome provides insight into its regulation. Gen Comp Endocrinol (2013) 191:45–58. doi: 10.1016/j.ygcen.2013.05.025
35. Teitelbaum SL. Bone resorption by osteoclasts. Science (2000) 289(September):1504–7. doi: 10.1126/science.289.5484.1504
36. To TT, Witten PE, Huysseune A, Winkler C. An adult osteopetrosis model in medaka reveals the importance of osteoclast function for bone remodeling in teleost fish. Comp Biochem Physiol Part - C: Toxicol Pharmacol (2015) 178:68–75. doi: 10.1016/j.cbpc.2015.08.007
37. Tagliaferri C, Wittrant Y, Davicco MJ, Walrand S, Coxam V. Muscle and bone, two interconnected tissues. Ageing Res Rev (2015) 21:55–70. doi: 10.1016/j.arr.2015.03.002
38. Kaji H. Effects of myokines on bone. BoneKEy Rep (2016) 5(July):1–6. doi: 10.1038/bonekey.2016.48
39. Kirk B, Feehan J, Lombardi G, Duque G. Muscle, bone, and fat crosstalk: the biological role of myokines, osteokines, and adipokines. Curr Osteoporos. Rep (2020) 18(4):388–400. doi: 10.1007/s11914-020-00599-y
40. Brotto M, Bonewald L. Bone and muscle: Interactions beyond mechanical. Bone (2015) 80:109–14. doi: 10.1016/j.bone.2015.02.010
41. Lara-Castillo N, Johnson ML. Bone-muscle mutual interactions. Curr Osteoporos. Rep (2020) 18(4):408–21. doi: 10.1007/s11914-020-00602-6
43. Pfaffl MW. A new mathematical model for relative quantification in real-time RT–PCR. Nuleic. Acids Res (2001) 29(9):e45. doi: 10.1093/nar/29.9.e45.
44. Rowlerson A, Radaelli G, Mascarello F, Veggetti A. Regeneration of skeletal muscle in two teleost fish: Sparus aurata and brachydanio rerio. Cell Tissue Res (1997) 289(2):311–22. doi: 10.1007/s004410050878
45. le Roith D, Scavo L, Butler A. What is the role of circulating IGF-I? Trends Endocrinol Metab (2001) 12(2):48–52. doi: 10.1016/s1043-2760(00)00349-0
46. Louie KW, Saera-Vila A, Kish PE, Colacino JA, Kahana A. Temporally distinct transcriptional regulation of myocyte dedifferentiation and myofiber growth during muscle regeneration. BMC Genomics (2017) 18(1):1–17. doi: 10.1186/s12864-017-4236-y
47. Florini JR, Ewton DZ, Coolican SA. Growth hormone and the insulin-like growth factor system in myogenesis. Endocrine. Rev (1996) 17(5):423–80. doi: 10.1210/er.17.5.423
48. Wilson EM, Rotwein P. Control of MyoD function during initiation of muscle differentiation by an autocrine signaling pathway activated by insulin-like growth factor-II. J Biol Chem (2006) 281(40):29962–71. doi: 10.1074/jbc.M605445200
49. Wilson EM, Rotwein P. Selective control of skeletal muscle differentiation by Akt1. J Biol Chem (2007) 282(8):5106–10. doi: 10.1074/jbc.C600315200
50. Ge Y, Wu AL, Warnes C, Liu J, Zhang C, Kawasome H, et al. mTOR regulates skeletal muscle regeneration in vivo through kinase-dependent and kinase-independent mechanisms. Am J Physiol - Cell Physiol (2009) 297(6):1434–44. doi: 10.1152/ajpcell.00248.2009
51. Jiménez-Amilburu V, Salmerón C, Codina M, Navarro I, Capilla E, Gutiérrez J. Insulin-like growth factors effects on the expression of myogenic regulatory factors in gilthead sea bream muscle cells. Gen Comp Endocrinol (2013) 188(1):151–8. doi: 10.1016/j.ygcen.2013.02.033
52. Vélez EJ, Lutfi E, Jiménez-Amilburu V, Riera-Codina M, Capilla E, Navarro I, et al. IGF-I and amino acids effects through TOR signaling on proliferation and differentiation of gilthead sea bream cultured myocytes. Gen Comp Endocrinol (2014) 205:296–304. doi: 10.1016/j.ygcen.2014.05.024
53. Pérez-Sánchez J. The involvement of growth hormone in growth regulation, energy homeostasis and immune function in the gilthead sea bream (Sparus aurata): A short review. Fish. Physiol Biochem (2000) 22(2):135–44. doi: 10.1023/A:1007816015345
54. Duguez S, le Bihan MC, Gouttefangeas D, Féasson L, Freyssenet D. Myogenic and nonmyogenic cells differentially express proteinases, Hsc/Hsp70, and BAG-1 during skeletal muscle regeneration. Am J Physiol - Endocrinol Metab (2003) 285(1 48-1):206–15. doi: 10.1152/ajpendo.00331.2002
55. Perelló-Amorós M, García-Pérez I, Sánchez-Moya A, Innamorati A, Vélez EJ, Achaerandio I, et al. Diet and exercise modulate GH-IGFs axis, proteolytic markers and myogenic regulatory factors in juveniles of gilthead Sea bream (Sparus aurata). Animals (2021) 11(8):2182. doi: 10.3390/ani11082182
56. Vélez EJ, Azizi S, Lutfi E, Capilla E, Moya A, Navarro I, et al. Moderate and sustained exercise modulates muscle proteolytic and myogenic markers in gilthead sea bream (Sparus aurata). Am J Physiol - Regul Integr Comp Physiol (2017) 312(5):R643–53. doi: 10.1152/ajpregu.00308.2016
57. Lavajoo F, Perelló-Amorós M, Vélez EJ, Sánchez-Moya A, Balbuena-Pecino S, Riera-Heredia N, et al. Regulatory mechanisms involved in muscle and bone remodeling during refeeding in gilthead sea bream. Sci Rep (2020) 10(1):1–14. doi: 10.1038/s41598-019-57013-6
58. Vélez EJ, Azizi S, Verheyden D, Salmerón C, Lutfi E, Sánchez-Moya A, et al. Proteolytic systems’ expression during myogenesis and transcriptional regulation by amino acids in gilthead sea bream cultured muscle cells. PLos One (2017) 12(12):1–21. doi: 10.1371/journal.pone.0187339
59. Tintignac LA, Lagirand J, Batonnet S, Sirri V, Leibovitch MP, Leibovitch SA. Degradation of MyoD mediated by the SCF (MAFbx) ubiquitin ligase. J Biol Chem (2005) 280(4):2847–56. doi: 10.1074/jbc.M411346200
60. Taillandier D, Combaret L, Pouch M-N, Samuels SE, Béchet D, Attaix D. The role of ubiquitin–proteasome-dependent proteolysis in the remodelling of skeletal muscle. Proc Nutr Soc (2004) 63(2):357–61. doi: 10.1079/par2004358
61. Huang J, Forsberg NE. Role of calpain in skeletal-muscle protein degradation. Proc Natl Acad Sci United States America (1998) 95(21):12100–5. doi: 10.1073/pnas.95.21.12100
62. Saido TC, Sorimachi H, Suzuki K. Calpain : new perspectives in molecular involvement diversity. FASEB J (1994) 8(11):814–22. doi: 10.1096/fasebj.8.11.8070630
63. Balcerzak D, Poussard S, Brustis JJ, Elamrani N, Soriano M, Cottin P, et al. An antisense oligodeoxyribonucleotide to m-calpain mRNA inhibits myoblast fusion. J Cell Sci (1995) 108(5):2077–82. doi: 10.1242/jcs.108.5.2077
64. Balcerzak D, Cottin P, Poussard S, Cucuron A, Brustis JJ, Ducastaing A. Calpastatin-modulation of m-calpain activity is required for myoblast fusion. Eur J Cell Biol (1998) 75(3):247–53. doi: 10.1016/S0171-9335(98)80120-9
65. Perelló-Amorós M, Otero-Tarrazón A, Jorge-Pedraza V, García-Pérez I, Sánchez-Moya A, Gabillard J-C, et al. Myomaker and myomixer characterization in gilthead sea bream under different myogenesis conditions. International Journal of Molecular Sciences (2022) 23(23):14639. doi: 10.21203/rs.3.rs-1123861/v1
66. Manneken JD, Dauer MVP, Currie PD. Dynamics of muscle growth and regeneration: Lessons from the teleost. Exp Cell Res (2022) 411(2):112991. doi: 10.1016/j.yexcr.2021.112991
67. Forcina L, Cosentino M, Musarò A. Mechanisms regulating muscle regeneration: Insights into the interrelated and time-dependent phases of tissue healing. Cells (2020) 9(5):1297. doi: 10.3390/cells9051297
68. Perello-Amorós M, Rallière C, Gutiérrez J, Gabillard JC. Myomixer is expressed during embryonic and post-larval hyperplasia, muscle regeneration and differentiation of myoblats in rainbow trout (Oncorhynchus mykiss). Gene (2021) 790(April):145688. doi: 10.1016/j.gene.2021.145688
69. Williams TM, Lisanti MP. The caveolin proteins. Genome Biol (2004) 5(3):1–8. doi: 10.1186/gb-2004-5-3-214
70. Costa RA, Martins RST, Capilla E, Anjos L, Power DM. Vertebrate SLRP family evolution and the subfunctionalization of osteoglycin gene duplicates in teleost fish. BMC Evolutionary Biol (2018) 18(1):1–17. doi: 10.1186/s12862-018-1310-2
71. Vélez EJ, Lutfi E, Azizi S, Perelló M, Salmerón C, Riera-Codina M, et al. Understanding fish muscle growth regulation to optimize aquaculture production. Aquaculture (2017) 467:28–40. doi: 10.1016/j.aquaculture.2016.07.004
72. Gabillard JC, Sabin N, Paboeuf G. In vitro characterization of proliferation and differentiation of trout satellite cells. Cell Tissue Res (2010) 342(3):471–7. doi: 10.1007/s00441-010-1071-8
73. Schmidt JG, Andersen EW, Ersbøll BK, Nielsen ME. Muscle wound healing in rainbow trout (Oncorhynchus mykiss). Fish. Shellfish Immunol (2016) 48:273–84. doi: 10.1016/j.fsi.2015.12.010
74. Chazaud B. Inflammation and skeletal muscle regeneration: Leave it to the macrophages! Trends Immunol (2020) 41(6):481–92. doi: 10.1016/j.it.2020.04.006
75. Chazaud B, Sonnet C, Lafuste P, Bassez G, Rimaniol AC, Poron F, et al. Satellite cells attract monocytes and use macrophages as a support to escape apoptosis and enhance muscle growth. J Cell Biol (2003) 163(5):1133–43. doi: 10.1083/jcb.200212046
76. Zhang C, Li Y, Wu Y, Wang L, Wang X, Du J. Interleukin-6/signal transducer and activator of transcription 3 (STAT3) pathway is essential for macrophage infiltration and myoblast proliferation during muscle regeneration. J Biol Chem (2013) 288(3):1489–99. doi: 10.1074/jbc.M112.419788
77. Ebisui C, Tsujinaka T, Morimoto T, Kan K, Iijima S, Yano M, et al. Interleukin-6 induces proteolysis by activating intracellular proteases (cathepsins b and l, proteasome) in C2C12 myotubes. Clin Sci (1995) 89(4):431–9. doi: 10.1042/cs0890431
78. Tsujinaka T, Fujita J, Ebisui C, Yano M, Kominami E, Suzuki K, et al. Interleukin 6 receptor antibody inhibits muscle atrophy and modulates proteolytic systems in interleukin 6 transgenic mice. J Clin Invest (1996) 97(1):244–9. doi: 10.1172/JCI118398
79. Kang X, Yang MY, Shi YX, Xie MM, Zhu M, Zheng XL, et al. Interleukin-15 facilitates muscle regeneration through modulation of fibro/adipogenic progenitors. Cell Communication Signaling (2018) 16(1):1–11. doi: 10.1186/s12964-018-0251-0
80. O’Leary MF, Wallace GR, Bennett AJ, Tsintzas K, Jones SW. IL-15 promotes human myogenesis and mitigates the detrimental effects of TNFα on myotube development. Sci Rep (2017) 7(1):1–11. doi: 10.1038/s41598-017-13479-w
81. Chen S, Gerken E, Zhang Y, Zhan M, Mohan RK, Li AS, et al. Role of TNF-{alpha} signaling in regeneration of cardiotoxin-injured muscle. Am J Physiol Cell Physiol (2005) 289(5):C1179–87. doi: 10.1152/ajpcell.00062.2005
82. Warren GL, Hulderman T, Jensen N, McKinstry M, Mishra M, Luster MI, et al. Physiological role of tumor necrosis factor alpha in traumatic muscle injury. FASEB Journal : Off Publ Fed Am Societies Exp Biol (2002) 16(12):1630–2. doi: 10.1096/fj.02-0187fje
83. Chen SE, Jin B, Li YP. TNF-α regulates myogenesis and muscle regeneration by activating p38 MAPK. Am J Physiol - Cell Physiol (2007) 292(5):C1660–71. doi: 10.1152/ajpcell.00486.2006
84. Arnold L, Henry A, Poron F, Baba-Amer Y, van Rooijen N, Plonquet A, et al. Inflammatory monocytes recruited after skeletal muscle injury switch into antiinflammatory macrophages to support myogenesis. J Exp Med (2007) 204(5):1057–69. doi: 10.1084/jem.20070075
85. Tidball JG. Inflammatory processes in muscle injury and repair. Am J Physiol - Regul Integr Comp Physiol (2005) 288(2 57-2):345–53. doi: 10.1152/ajpregu.00454.2004
86. Tidball JG. Regulation of muscle growth and regeneration by the immune system. Nat Rev Immunol (2017) 17(3):165–78. doi: 10.1038/nri.2016.150.Regulation
87. Hamilton JA. Colony-stimulating factors in inflammation and autoimmunity. Nat Rev Immunol (2008) 8(7):533–44. doi: 10.1038/nri2356
88. Tanaka KI, Matsumoto E, Higashimaki Y, Katagiri T, Sugimoto T, Seino S, et al. Role of osteoglycin in the linkage between muscle and bone. J Biol Chem (2012) 287(15):11616–28. doi: 10.1074/jbc.M111.292193
89. Karsenty G, Mera P. Molecular bases of the crosstalk between bone and muscle. Bone (2018) 115(1):43–9. doi: 10.1016/j.bone.2017.04.006
90. Battafarano G, Rossi M, Marampon F, Minisola S, del Fattore A. Bone control of muscle function. Int J Mol Sci (2020) 21(4):1–14. doi: 10.3390/ijms21041178
91. Mera P, Laue K, Ferron M, Confavreux C, Wei J, Lacampagne A, et al. (2016). Osteocalcin Signaling in Myofibers Is Necessary and Sufficient for Optimum Adaptation to Exercise. Cell Metabolism 23(6), 1078–1092. doi: 10.1016/j.cmet.2016.05.004
92. Riera-Heredia N, Vélez EJ, Gutiérrez J, Navarro I, Capilla E. Gene expression analyses in malformed skeletal structures of gilthead sea bream (Sparus aurata). J Fish. Dis (2019) 42(8):1169–80. doi: 10.1111/jfd.13019
93. Ytteborg E, Baeverfjord G, Torgersen J, Hjelde K, Takle H. Molecular pathology of vertebral deformities in hyperthermic atlantic salmon (Salmo salar). BMC Physiol (2010) 10(1):2. doi: 10.1186/1472-6793-10-12
Keywords: skeletal muscle, injury, regeneration, GH-IGFs axis, proteolytic systems, myogenesis, crosstalk, bone
Citation: Otero-Tarrazón A, Perelló-Amorós M, Jorge-Pedraza V, Moshayedi F, Sánchez-Moya A, García-Pérez I, Fernández-Borràs J, García de la serrana D, Navarro I, Blasco J, Capilla E and Gutierrez J (2023) Muscle regeneration in gilthead sea bream: Implications of endocrine and local regulatory factors and the crosstalk with bone. Front. Endocrinol. 14:1101356. doi: 10.3389/fendo.2023.1101356
Received: 17 November 2022; Accepted: 03 January 2023;
Published: 23 January 2023.
Edited by:
Hamid R. Habibi, University of Calgary, CanadaReviewed by:
Peggy Biga, University of Alabama at Birmingham, United StatesCunming Duan, University of Michigan, United States
Copyright © 2023 Otero-Tarrazón, Perelló-Amorós, Jorge-Pedraza, Moshayedi, Sánchez-Moya, García-Pérez, Fernández-Borràs, García de la serrana, Navarro, Blasco, Capilla and Gutierrez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Joaquin Gutierrez, jgutierrez@ub.edu
†These authors have contributed equally to this work
‡ORCID: Daniel García de la serrana, orcid.org/0000-0001-6630-3379