Abstract
In muscular dystrophies, muscle fibers loose integrity and die, causing significant suffering and premature death. Strikingly, the extraocular muscles (EOMs) are spared, functioning well despite the disease progression. Although EOMs have been shown to differ from body musculature, the mechanisms underlying this inherent resistance to muscle dystrophies remain unknown. Here, we demonstrate important differences in gene expression as a response to muscle dystrophies between the EOMs and trunk muscles in zebrafish via transcriptomic profiling. We show that the LIM-protein Fhl2 is increased in response to the knockout of desmin, plectin and obscurin, cytoskeletal proteins whose knockout causes different muscle dystrophies, and contributes to disease protection of the EOMs. Moreover, we show that ectopic expression of fhl2b can partially rescue the muscle phenotype in the zebrafish Duchenne muscular dystrophy model sapje, significantly improving their survival. Therefore, Fhl2 is a protective agent and a candidate target gene for therapy of muscular dystrophies.
Similar content being viewed by others
Introduction
Muscular dystrophies, caused by mutations in more than 40 genes1 share similar features including muscle fiber disruption leading to muscle weakness, loss of ambulation, and often premature death predominantly due to respiratory failure. Although muscular dystrophies affect 1:4–5000 births worldwide2 and induce pronounced suffering, there is currently no cure, and treatment options are sparse. Thus, there is a need for new effective treatment options to prolong and improve the quality of life of patients suffering from these detrimental diseases.
The most common form of muscular dystrophy is Duchenne muscular dystrophy (DMD). In DMD, the giant protein dystrophin is lost or truncated3,4. Dystrophin is a crucial member of the dystrophin-glycoprotein-complex (DGC). The DGC links the extracellular matrix, across the cell membrane, to F-actin in the cytoskeleton within the myofibrils, which is fundamental for myofiber integrity5. The DGC contributes to multiple functions in the myofiber such as force transmission across the sarcolemma, but also acts as a signaling hub6,7,8. Dystrophin is therefore a key element for myofiber integrity. The most frequently used zebrafish DMD model is the sapje line9,10. Sapje zebrafish carry an A-to-T transversion in exon 4 of the dystrophin gene, resulting in a premature stop codon9. The lack of dystrophin in zebrafish subsequently leads to detachment of trunk myofibers from the myosepta and failure of the contractile apparatus which ultimately results in the premature death of up to 50% of zebrafish larvae at the age of 5 days post fertilization (dpf)11. Hence, the sapje zebrafish are severely affected by the lack of dystrophin, but generally mimic the human DMD condition and constitute a good model for experimental treatment studies.
The EOMs have shown an innate resistance towards muscular dystrophies12,13,14,15. Even though EOMs share attributes with other striated skeletal muscles they differ in terms of gene expression and protein content16,17 as well as neuromuscular and myotendinous junction composition18,19. We have previously shown that zebrafish EOMs are a good model to study the cytoskeleton as well as myofiber and neuromuscular junction cytoarchitecture20. Recently, the muscle-specific intermediate filament protein desmin was shown to be naturally lacking in a subset of EOM myofibers in humans and zebrafish20,21. Additionally, other intermediate filament proteins such as nestin and keratin-19 show a complex pattern in human EOM myotendinous junctions19. Altogether these findings suggest that the cytoskeletal composition of the EOMs differs from that of other skeletal muscles in the body.
Desmin is the most abundant cytoskeletal protein in skeletal muscle fibers22 and, besides its interaction with the DGC, desmin has been suggested to contribute to gene regulation23 and also as a mechanosensor, utilizing mechanical stretch to trigger intracellular signaling24. Additionally, desmin anchors myonuclei and mitochondria to the sarcolemma and myofibrils25, reviewed in26,27. Despite this important role in the myofiber, studies of Des knockout mice have shown that the EOMs are relatively unaffected28. Patients with desminopathy have near normal life expectancy, with minor complications compared to DMD patients29,30. Desmin mutant animal models therefore offer a route to study the EOMs in a dystrophic background over extended periods of time. We hypothesized that the transcriptome of the EOMs in a dystrophic setting would provide information regarding their innate resistance towards muscular dystrophies, and that this could be utilized to rescue other dystrophic skeletal muscle.
In the current study, we used a desmin knockout zebrafish line to identify a subset of genes specifically upregulated in EOMs, via transcriptome analysis. We further investigated one of these genes, fhl2b, and demonstrate that fhl2b has an important role in the inherent resistance of EOMs towards muscular dystrophies. Additionally, we show that fhl2b significantly improves survival, muscle integrity, and function of sapje zebrafish making it a novel target in the treatment of muscular dystrophies.
Results
A zebrafish model of desminopathy displays skeletal muscle defects
To study the EOMs in a non-lethal muscular dystrophy context, we generated desma+/-;desmb+/- zebrafish with premature stop codons in exon 1 of both genes (Fig. S1a, b). These mutations lead to truncated desmin lacking α-helix-rod domains, essential for coil formation and tertiary structure formation. We confirmed that the desma−/−;desmb−/− double mutants, unlike desma+/-;desmb+/- controls, lacked desmin immunolabeling at 3 dpf (Fig. S1c). Next, we performed functional experiments to investigate the impact of lack of desmin in zebrafish. Larvae were reared in 1% methyl cellulose, to generate swimming resistance, between days 4 and 5. This triggered significant myofiber detachment and breaks, both in desma−/−;desmb−/− mutants and desma−/− single mutants, at the 10–12th somite level (Fig. 1a, b), whereas desmb−/− and desma+/-;desmb+/- displayed no myofiber damage (Fig. 1b), indicating that desma is the main contributor to muscle tensile strength. Previous studies have shown both desma and desmb expression in developing zebrafish somite musculature31. Therefore, to investigate the complete knockout of desmin genes and avoid potential confounding factors, we decided to continue our studies using only desma−/−;desmb−/− zebrafish. Force/tension relationship measurements revealed that 5-6 dpf desma−/−;desmb−/− zebrafish larvae display a significant decrease in trunk myofiber maximal force generation when compared to desma+/-;desmb+/- (Fig. 1c, d). Furthermore, desma−/−;desmb−/− mutants showed a significant decrease in spontaneous movement in comparison to desma+/-;desmb+/- larvae (Fig. 1e, f). This was not due to myofiber loss, as both genotypes had equal numbers of Tg(mylz2:EGFP) fast and Tg(smyhc1:tdTomato) slow twitch positive myofibers in cross sections of the trunk (Fig. S1d, e). Altogether, these results show that desmin contributes to myofiber integrity, and is needed to maintain proper function of muscle tissue in the embryonic zebrafish trunk.
Extraocular muscle integrity is preserved in the zebrafish desminopathy model
Patients with desminopathy are often asymptomatic until mid-30s29. Therefore, to characterize our desminopathy model at later stages, we reared desma−/−;desmb−/− mutants and controls under equal conditions until 20–24 months of age, approximately two thirds of domesticated zebrafish life span, at which point they were histologically analyzed. desma−/−;desmb−/− mutants were fertile and survived until adulthood when reared separately from siblings (Fig. S1f). A consistent feature of muscular dystrophies is a glycolytic shift in myofiber identity from fast to slow32,33. To assess whether the lack of desmin had an impact on myofiber identity, trunk and EOMs of 24 months old zebrafish were immunolabeled using antibodies against fast myosin light chain (F310) and against slow myosin heavy chain 1, 2 and 3 (S58)34. The zebrafish trunk consists of myofibers subdivided into compartments separated by thin layers of connective tissue, myosepta, easily identified on cross-sections (Fig. 1, dashed line in g, i, k, m). Slow myofibers are generally found in the most laterally positioned compartments, separated from fast myofibers by a hyperplastic growth zone positioned near the myosepta35. Adult zebrafish EOMs also display a distinct location of slow and fast myofibers, although not separated by a myoseptum20 (Fig. 1o, q). Cross-sections of trunk muscle of desma−/−;desmb−/− mutants revealed a decrease of fast F310 positive myofibers inside the slow domain (Fig. 1g, h, open arrowheads), a decrease in the slow myofiber diameter together with a significant increase in myofiber number in the slow domain, compared to wild type (WT) controls (desma+/+;desmb+/+) (Fig. 1I, j). Additionally, the number of slow myofibers inside the fast-medial compartments was significantly increased (Fig. 1I, j, arrowheads). Cross-sections of EOMs showed that they remained unaffected in terms of proportion of fast and slow myofibers (Fig. 1o–r). In summary, adult desma−/−;desmb−/− mutant trunk muscles show a glycolytic-fast to oxidative-slow metabolic shift, consistent with a muscular dystrophy phenotype whereas the EOMs remain unaffected in this regard.
To further define the desma−/−;desmb−/− muscle phenotype, we analyzed zebrafish trunk muscles and EOMs for signs of regeneration, proliferation and cell death. desma−/−;desmb−/− mutant trunk muscle showed a significantly increased proportion of both Pax7 positive nuclei (Fig. 1k, l) and PCNA positive nuclei (Fig. 1m, n) in the slow myofiber domains compared to controls. No change in TUNEL positive nuclei in the slow domains was observed, however, the fast myofiber domains contained clusters of myofibers where most nuclei were TUNEL positive (Fig. S1g, arrows), along with a significant increase in centrally positioned nuclei, a hallmark of muscular dystrophy (Fig. S1g, h, asterisk). In contrast, controls only appeared to have sporadic signs of DNA fragmentation dispersed among nuclei in the entire myofiber population (Fig. S1g). Interestingly, in EOM cross sections of desma−/−;desmb−/− mutants, we instead noted a significant decrease of PCNA positive myonuclei (Fig. 1s, t) and TUNEL positive myonuclei (Fig. 1u, v) compared to controls, suggesting a decreased cellular turn over in response to the lack of desmin in the EOMs. Overall, EOM myofiber composition remains unchanged despite the lack of desmin, whereas trunk muscle is significantly affected and shows clear signs of muscular dystrophy.
Upregulation of fhl2b is an EOM response to muscular dystrophy
To investigate the adaptations observed in desmin deficient EOMs, we performed RNA-sequencing of EOMs and trunk muscles at 5 and 20 months of age, representing pre- and symptomatic stages, respectively. For each stage, we obtained transcriptome profiles from both trunk muscles and EOMs from desma−/−;desmb−/− mutant and WT controls (Fig. S2a, schematic illustration). To identify genes involved specifically in the EOMs resistance towards muscular dystrophy, we performed differential expression analysis in desma−/−;desmb−/− compared to WT controls in EOMs (Fig. 2a, comparison I-II) and trunk (Fig. 2a, comparison III-IV) at both time points (Fig. S2a–e, Supplementary Data 1), additionally, we also compared EOMs to trunk within the same genetic background (Fig. 2, comparisons V-VIII, Fig. S2f–i, Supplementary Data 1). Subsequently, we performed Gene Ontology (GO) analysis and retrieved differentially expressed genes (DEGs) contained within the cellular compartment GO terms related to myofiber function/structure (Fig. S3a–f, Supplementary Data 2), which we further cross-compared (Fig. 2a, comparisons I-VIII). Interestingly, when analyzing 20-month-old EOMs, we found upregulation of several DEGs related to cytoskeletal rearrangement (Fig. 2a, comparison III). One of these genes, fhl2b, was consistently identified across comparisons III-VIII highlighting it as a potential candidate gene (Fig. 2a). Notably, four members of the fhl family, including fhl2b, were identified in all EOM vs trunk comparisons (Fig. 2a). We did not find fhl2b to be differentially expressed in 5 months old desma−/−;desmb−/− EOMs which can likely be attributed to the early stage of disease progression. Next, we analyzed fhl2b gene expression across the different comparisons (Fig. 2a, comparisons I-VIII) and found that fhl2b consistently was more highly expressed in EOMs compared to trunk muscle (Fig. 2b). Previous studies have shown that Fhl2 is mainly expressed in cardiac muscle36, localized to Z-discs, but has to our knowledge not been studied in the EOMs prior to this study. EOM longitudinal sections immunolabeled for Fhl2 showed localization to the Z-disc (Fig. S3g). The number of Fhl2-immunolabelled myofibers in EOMs of desma-/-;desmb-/- mutants was significantly increased relative to controls at both 5 and 20 months of age and increased significantly with age (Fig. 2c–e) confirming our transcriptomic analysis. Next, we asked whether an increase in Fhl2 positive myofibers of EOMs could be found across muscular dystrophies and therefore immunolabeled plecb−/− (plectin) and obscnb−/− (obscurin) mutant zebrafish EOMs (Fig. 2f), both resulting in muscular dystrophy when mutated in other models37,38. Interestingly, plecb−/− and obscnb−/− mutant zebrafish displayed increased numbers of Fhl2 positive myofibers compared to controls, similar to desma-/-;desmb-/- EOMs (Fig. 2c, d, f). This shows that Fhl2 is more widespread in several different cytoskeletal gene mutations and models for muscular dystrophy. Importantly, healthy adult human and mouse EOMs were also found to contain FHL2-positive myofibers (Fig. 2g, i), confirmed by western blots (Fig. 2h, j), indicating a putatively conserved role for Fhl2 in EOMs across species. In summary, Fhl2 is present in EOMs across different muscular dystrophies and species and is increased with the progression of the disease in the EOMs of desma-/-;desmb-/- mutant zebrafish.
Knockout of fhl2 leads to EOM myofiber hypertrophy
Previous studies have found Fhl2 to be mainly expressed in cardiac muscle, however, Fhl2 knockout mice were shown to be viable and no phenotype was observed unless challenged with isoproterenol, triggering adrenergic stimuli which resulted in cardiac hypertrophy39. We hypothesized that knockout of fhl2 in the background of desma−/−;desmb−/− mutations could cause similar phenotypes in EOMs of adult zebrafish. We generated knockout lines of both zebrafish fhl2 genes, fhl2a, and fhl2b, to avoid possible confounding effects of redundancy (Fig. S4a, b). In situ hybridization of fhl2a showed low EOM expression whereas fhl2b was found to be distinctly expressed in the EOMs at 5 dpf (Fig. S4c, d, arrowheads) suggesting a larger role for fhl2b compared to fhl2a in the EOMs. Immunolabeling using Fhl2 antibodies confirmed lack of Fhl2 in desma−/−;desmb−/−;fhl2a−/−;fhl2b−/− compared to desma−/−;desmb−/−;fhl2a+/-;fhl2b+/- sibling controls (Fig. S4e). To further elucidate the effect of lack of fhl2 genes on other fhl family members, we analyzed the gene expression of fhl1a, fhl1b, fhl3a and fhl3b in the different mutant larvae. There were no significant differences between WT and desma−/−;desmb−/− double mutants for these genes, however, we found that knockout of fhl2a, fhl2b or both genes rendered a general downregulation of fhl1a and fhl1b (Fig. S4f, g). fhl3a expression remained unchanged across all comparisons, however, a similar, less pronounced, downregulation trend was observed for fhl3b (Fig. S4h, i). These results are in line with previous findings40, and indicate that knockout of fhl2 genes is not compensated by other fhl family members. Instead, our results imply that the expression of fhl1 and fhl3 may be partly dependent on Fhl2. In addition, previous studies have linked Fhl2 to WNT signaling via β-catenin41,42,43,44 and we therefore investigated if knockout of fhl2 genes would affect this signaling pathway. Similar to the previous analysis, we found no significant differences between WT and desma−/−;desmb−/− double mutants. Our results showed that knockout of either fhl2a, fhl2b or both fhl2a and fhl2b significantly altered the levels of wnt5a, wnt5b, wnt11, and ctnnb2 (Fig. S4j–n). These data further support a role for Fhl2 in the WNT/ β-catenin signaling pathway.
Quantification of EOM myofiber size at 12 months of age in single fhl2a−/− and fhl2b−/− mutants and fhl2a−/−;fhl2b−/− double mutants, all in a desma-/-;desmb-/- background, alongside the corresponding controls revealed that the myofiber areas were significantly increased in desma−/−;desmb−/−;fhl2b−/− and desma−/−;desmb−/−;fhl2a−/−;fhl2b−/− mutants compared to WT controls (Fig. 3a, b). This indicates that fhl2b has a role in hypertrophic protection, in line with results observed for Fhl2 in mouse cardiac muscle39. Given the reduction on cell death and proliferation observed in desma−/−;desmb−/− mutant EOMs (Fig. 1s–v), we wondered whether the knockout of fhl2a and fhl2b would revert these effects. For this purpose, we addressed cell death and proliferation on EOM cross sections. TUNEL (Fig. 3c–e) and PCNA (Fig. 3f–h) labeling of myonuclei showed a redundant relationship between fhl2a and fhl2b. Both desma−/−;desmb−/−;fhl2a−/− and desma−/−;desmb−/−;fhl2b−/− triple mutants displayed no change in cell death whereas a moderate increase in proliferation was found. However, in desma−/−;desmb−/−;fhl2a−/−;fhl2b−/− quadruple mutant zebrafish EOMs, a highly significant increase in both cell death and proliferation was observed (Fig. 3e, h). Given that hypertrophic myofibers were observed in desma−/−;desmb−/−;fhl2b−/− and desma−/−;desmb−/−;fhl2a−/−;fhl2b−/− EOMs at 12 months, we further analyzed 12-month-old cardiac and skeletal muscle tissues for hypertrophy in our mutants, however, no significant changes were observed at this stage (Fig. S5a, h). These results indicate that both fhl2a and fhl2b maintain myonuclei integrity in the EOMs under muscle dystrophy conditions, however, fhl2b likely has an additional role in hypertrophic protection given the increase in myofiber area observed in the examined mutants lacking fhl2b (Fig. 3a, b). Collectively, we show that Fhl2 is needed to maintain myofiber homeostasis in the EOMs in desmin knockout background and can therefore be considered as an EOM protective gene.
Muscle specific overexpression of fhl2b significantly improves survival and myofiber integrity of the Duchenne muscular dystrophy zebrafish model
Since fhl2b was the only Fhl2 paralog differentially expressed in our transcriptomic data of EOMs lacking desmin, we chose to investigate whether fhl2b also could protect muscles other than EOMs in muscular dystrophy. To test this in a more severe context, we utilized the lethal zebrafish sapje (dmdta222a) line9, a model of Duchenne muscular dystrophy, hereby referred to as dmd−/−. We overexpressed fhl2b under the muscle-specific 503unc promoter45 coupled to EGFP via a T2A linker and analyzed trunk muscle in the dmd−/− background (Fig. 4a). We initially examined mosaic overexpression of fhl2b, which showed an overlap between EGFP positive myofibers and Fhl2 antibody positive myofibers (Fig. S6a). dmd−/− larvae with a high number of EGFP positive myofibers generally showed less myofiber disruption than seen in dmd−/− larvae with low number EGFP positive fibers (Fig. S6b). We therefore reared mosaic larvae to adulthood and generated stable overexpression transgenic lines from three different founders with different EGFP intensity levels, one low and two high (Fig. S6c). To functionally test the effects of fhl2b overexpression on dmd−/− larvae, we analyzed survival among low and high level fhl2b expressing lines and found a significant correlation between fhl2b expression (Fig. S4c) and survival of dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae (Fig. 4b). The median survival time of dmd−/− larvae was found to be 18 dpf, whereas dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) and sibling control median survival time was longer than 30 days (Fig. 4b). Notably, even low levels of fhl2b overexpression increased survival beyond 30 days in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP). Additionally, spontaneous swimming distance was also significantly increased in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae compared to dmd−/− larvae at 5 dpf (Fig. 4c, d). Collectively, these results show that muscle specific fhl2b overexpression improves motor function and significantly prolongs the lifespan in this model of Duchenne muscular dystrophy.
To examine if fhl2b overexpression improved myofiber integrity, we labeled 5 dpf dmd−/−;Tg(503unc:fhl2b-T2A-EGFP), dmd−/− and sibling control (dmd+/+, dmd+/-) larvae using phalloidin and DAPI. Additionally, we generated a Tg(503unc:EGFP) line which was crossed into the dmd+/- line, used as a negative control. dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae showed detached myofibers similar to dmd−/− and dmd−/−;Tg(503unc:EGFP) larvae (Fig. 4e, f). However, dmd−/− and dmd−/−;Tg(503unc:EGFP) detached myofibers left gaps in the muscle devoid of F-actin (Fig. 4f, closed arrowheads). Interestingly, detached myofibers in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) myotomes seldom had these gaps and thin F-actin positive myofibers were present (Fig. 4f, open arrowheads) which exhibited a high EGFP intensity, suggesting that they are newly formed as the 503unc promoter is more active in early stages of myogenesis46 (Fig. 4g). We quantified myofiber detachment in dmd−/− and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae using phalloidin. We found that dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) consistently showed reduced number of myofiber detachments per somite compared to dmd−/− controls (Fig. 4h, i). Overall, these data indicate that dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) myofibers have a lower tendency to detach compared to dmd−/− myofibers.
Muscle specific fhl2b overexpression improves motor axon integrity and neuromuscular junctions in dmd −/− larvae
To further understand how overexpression of fhl2b improves the dmd−/− phenotype, we analyzed the transcriptomes of trunk muscle tissue from sibling controls (dmd+/+, dmd+/-), sibling controls with fhl2b overexpression (Tg(503unc:fhl2b-T2A-EGFP)), dmd−/− and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae using RNA-sequencing at 5 dpf (Fig. S6d–j, Supplementary Data 3). We then intersected the DEGs from the pairwise comparisons dmd−/− vs sibling controls and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) vs sibling controls and identified 1054 DEGs unique to dmd−/− larvae (Fig. 5a, gene set C). These 1054 DEGs correspond to dmd disease-related genes, whose expression is partially rescued by fhl2b overexpression (Fig. 5b). Gene ontology analysis of all intersections (Fig. 5c, Fig. S7a, b) revealed a unique enrichment for axon and neuron guidance-related terms in dmd−/− larvae (Fig. 5c), which notably included several semaphorin genes, crucial for axon guidance and motor neuron survival (Fig. 5d). In addition, mir-206, a regulatory microRNA suggested to mediate communication between myofibers and neurons, and histone deacetylase 4 (hdac4) are both upregulated in Duchenne muscular dystrophy patients and suggested as biomarkers47,48. We determined the expression of these genes in our system and observed that they were significantly upregulated in dmd−/− larvae (Fig. 5e, f). However, in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae, the expression levels of both genes were significantly closer to that of sibling controls, which comparatively might be indicative of healthier tissue (Fig. 5e, f). To test whether axon and neuromuscular junction (NMJ) integrity were improved in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae, we immunolabeled 5 dpf larvae using acetylated tubulin antibodies and α-bungarotoxin, labeling axons and post-synaptic NMJs, respectively (Fig. 5g–i, Fig. S8). dmd−/− larvae displayed thin axons with reduced branching compared to sibling controls (Fig. 5g, i, arrowhead), additionally, NMJs were disorganized and fragmented (Fig. 5h, i arrow) and an almost complete loss of contact was registered between axons and NMJs. In contrast, dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae displayed near normal axonal branching (Fig. 5g) and NMJ organization (Fig. 5i, open arrowhead). Because the dmd−/− NMJs were extremely fragmented, we defined healthy NMJs as α-bungarotoxin positive puncta larger than the smallest WT puncta. The healthy NMJs were significantly reduced in dmd−/−, but not in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae (Fig. 5k). Additionally, double positive Synaptic vesicle 2 (SV2)/α-bungarotoxin NMJs were significantly fewer in dmd−/− as compared to dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) and sibling control larvae (Fig. 5j arrowhead, Fig. 5l), indicating decreased innervation of myofibers. To quantify axon integrity, we measured the length (Fig. 5m, o) and volume (Fig. 5n, p) of axons using 3D-renderings of acetylated tubulin immunolabeled larvae and found significant improvements in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) compared to dmd−/− larvae. Together, our data indicate that axon retraction or degradation is a commonly occurring phenomenon in dmd−/− larvae. Furthermore, our data suggest that in myofibers rescued by fhl2b expression present a significant improvement in axons, NMJs, innervation and functionality as evidenced by swimming capacity (Fig. 4d). However, we recognize that this may be the consequence of improved muscle tissue integrity.
Muscle specific overexpression of fhl2b accelerates macrophage activity and muscle regeneration
To further define mechanisms underlying survival of dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae, we intersected DEGs from dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) vs dmd−/− and Tg(503unc:fhl2b-T2A-EGFP) vs sibling controls (dmd+/+, dmd+/-) (Fig. 6a, Fig. S7c). We identified 78 genes differentially expressed in response to fhl2b overexpression in the dmd−/− disease background (Fig. 6a, gene set F). Among these 78 DEGs, we found several genes involved in muscle regeneration (Fig. 6b, flnca, rbfoxl1, ptgs2b, kdm6ba, mmp14b), suggesting that this is a possible mechanism for improved muscle integrity in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae. In addition, we analyzed 5 dpf dmd−/−;Tg(503unc:fhl2b-T2A-EGFP), dmd−/−, dmd−/−;Tg(503unc:EGFP) and sibling controls for muscle regeneration, proliferation and cell death (Fig. S9a–g). Interestingly, we noted a significant decrease in Pax7 (Fig. S9a, b), BrdU (24 h pulse, Fig. S9c–e) and TUNEL (Fig. S9f–g) positive nuclei in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) compared to dmd−/− larvae, indicating overall healthier myofibers. These data confirm that dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) myofibers are partially rescued by fhl2b overexpression and support our findings that dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) myofibers are less prone to detachment. Interestingly, our study indicates that fhl2b induces faster formation of new myofibers and we therefore hypothesized that the muscle regeneration would be quicker in fhl2b overexpressing larvae. To address this, we performed timelapse studies on dmd−/−;Tg(503unc:EGFP) and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae during 24 h to observe muscle regeneration in real time (Fig. 6c). We observed that detached myofiber debris was quickly removed and myofibers were apparently replaced in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae whereas practically no change was observed in dmd−/−;Tg(503unc:EGFP) controls during the 24 h window (Fig. 6c). To quantify this, we measured muscle integrity using birefringence. We found a significant improvement in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) over time as compared to dmd−/−;Tg(503unc:EGFP) controls (Fig. 6d, e). Collectively, these results suggest that fhl2b overexpression accelerates muscle regeneration.
Given our results showing that detached myofibers quickly disappeared in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae (Fig. 6c), we hypothesized that increased leukocyte activity could be involved in the muscle regeneration process. Leukocytes have previously been shown to be critical in zebrafish muscle regeneration49. Additionally, macrophage specific genes mmp13a and mpeg1.2 as well as il-8 were upregulated in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) and uninjured Tg(503unc:fhl2b-T2A-EGFP) as compared to sibling controls (Fig. 7a), supporting this hypothesis. To address this, we analyzed fhl2b´s effect on various aspects of muscle regeneration in a controlled setting using Tg(503unc:fhl2b-T2A-EGFP) larvae compared to sibling controls. We performed needle-stick injuries and analyzed the muscle regeneration process by immunolabeling for neutrophils (Mpx), macrophages (Mfap4) and satellite cells (Pax7) coupled with DAPI and phalloidin, over time (Fig. 7b). All three antibodies were found to label an equal number of cells in 3 dpf uninjured controls (Fig. S.9h-m), and the neutrophil population labeled by Mpx antibodies was unchanged across all stages tested, from uninjured to 24 h post injury (hpi) (Fig. 7c, d). Noticeably, we found a significant increase in macrophage count at the wound site in Tg(503unc:fhl2b-T2A-EGFP) larvae already at 6 hpi (Fig. 7e, f). This was followed by a significantly reduced number of macrophages at 24 hpi in Tg(503unc:fhl2b-T2A-EGFP) wounds compared to sibling control wounds (Fig. 7e, f). Pax7 positive satellite cells were present at wound sites faster in Tg(503unc:fhl2b-T2A-EGFP) compared to sibling controls and peaked at 48 hpi whereas sibling control satellite cell numbers peaked at 72 hpi (Fig. 7g, h). Additionally, F-actin intensity was found to be significantly higher at 72 hpi as compared to controls (Fig. 7I, j). Taken together, these results suggest that Tg(503unc:fhl2b-T2A-EGFP) larvae have faster myofiber regeneration via enhanced macrophage recruitment.
Discussion
In the current study we show that induced expression of fhl2b can provide resistance to muscular dystrophy and that the EOMs offer a novel approach regarding protective cellular strategies. Here we propose a model (Fig. 8) in which dmd−/− larvae exhibit prolonged survival due to protection by fhl2b, muscle injury does not occur as frequently or to the same extent and axons and neuromuscular junctions are less affected. Additionally, fhl2b improves macrophage response to myofiber injury, thereby accelerating recovery from myofiber damage. We therefore propose that Fhl2 based therapy has the potential to be utilized to alleviate muscular dystrophy in body musculature.
Fhl2 is a co-transcription factor and can translocate to the nucleus to co-activate gene expression50 and has been shown to be expressed in myogenic precursor cells to promote differentiation44 and in Pax7 positive satellite cells during muscle regeneration51,52. However, we did not detect any Fhl2 positive myonuclei in the EOMs, instead Fhl2 was present mainly on the Z-discs of myofibers and likely executes its function there, similar to what has been observed in rat myocardiocytes53. fhl2 is not widely expressed in muscle derived from the paraxial mesoderm in zebrafish between mid somitogenesis to 5 dpf46 and its knockout did not have obvious consequences on muscle development in our study, and others39,40, arguing against a primary role in skeletal muscle development. It has been shown that Fhl2 inhibits β-adrenergic stimulated calcineurin/NFAT signaling53. However, we did not find any DEGs known to be involved in this signaling cascade in desma−/−;desmb−/− EOMs In rodent cardiomyocytes, Fhl2 overexpression has also been shown to inhibit SRF/ERK pathway activity54 and block ERK2 translocation to the nucleus and subsequent transcriptional activity, offering protection from hypertrophy55. The WNT/β-catenin signaling pathway has previously also been linked to both Fhl2 and hypertrophy43,44,56,57. Fhl2 can act both to repress β-catenin mediated transcription in muscle as well as to co-activate gene transcription in kidney and colon cell lines, showcasing Fhl2’s diverse functions depending on cellular context. In this study, knockout of fhl2 genes led to increased expression of ctnnb2 (β-catenin), indicating a repressive function of fhl2 on ctnnb2 in our mutants. Furthermore, knockout of fhl2 genes also affected the expression of wnt5a, wnt5b and wnt11, in line with previous studies58, suggesting a role for Fhl2 in both canonical and non-canonical WNT-signaling in zebrafish, which may explain the hypertrophic changes observed in EOMs lacking Fhl2. As zebrafish hearts lacking Fhl2 were not hypertrophic, and hypertrophy in mice hearts lacking Fhl2 show increased calcineurin/NFAT signaling53,54, it is therefore possible that the function of Fhl2 and the process of hypertrophy differs in EOMs compared to cardiac muscles.
Fhl2 has been linked to metabolic enzyme targeting to the N2B region of titin, suggesting that it acts as an adaptor protein for high energy demanding regions within the sarcomere59. The N2B region has also been proposed to be an anchoring hub for mechanotransduction60 and Fhl2 has previously been shown to be important for mechanotransduction in several scenarios where it is situated on F-actin and can translocate to the nucleus upon a change of tension to regulate gene expression50,61. Interestingly, zyx, pdlim2 and ldb3a, identified in our transcriptional profiling of dystrophic EOMs, are also linked to mechanotransduction62,63,64, introducing mechanosensing as a potentially important mechanism for EOM specific resistance to muscular dystrophy. This notion is strengthened by the fact that we found Fhl2 to be more abundant in a number of muscular dystrophy models, suggesting it is part of a general mechanism of defense in the EOMs. Knockout of fhl2 genes leads to increased levels of cellular turn around and hypertrophy, emphasizing the importance of Fhl2 in EOM myofiber homeostasis, similar to what has been shown in rodent cardiomyocytes39,55. There was no redundancy between Fhl family members, however, we cannot exclude that the downregulation of fhl1 genes in the quadruple mutant larvae might contribute to the cellular turnaround phenotype in the EOMs. We show that ectopic fhl2b expression in dmd−/− larvae leads to protection from myofiber integrity loss. However, whether ectopic fhl2b has the same function in trunk muscle, EOMs and cardiomyocytes remains to be determined.
Duchenne muscular dystrophy models have shown alterations in NMJ patterning65, non-surprising considering that the dystrophin-glycoprotein-complex (DGC) show accumulation at the NMJs in healthy muscle. Interestingly, our RNA-sequencing and in vivo data clearly demonstrate preserved expression of several genes related to axon and motor neuron guidance and a clear improvement in axon and NMJ structure in dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) compared to dmd−/− larvae. fhl2b overexpressing larvae notably showed higher levels of axon guidance cues such as semaphorins, one of the major families of axon guidance molecules66. While muscle-secreted molecules have been shown to affect motor neurons directly67, it is difficult to assess whether our results are a direct effect of fhl2b overexpression or a consequence of preserved myofiber integrity. Denervated myofibers undergo degeneration68, likely adding to the overall muscle phenotype observed in dmd−/− larvae. Previous studies have shown the beneficial effects of electrical muscle stimulation on DMD-deficient myofibers69. In our study, fhl2b overexpressing dmd−/− larvae showed improved motor function, which is likely a consequence of both improved muscle and axon/NMJ structure. The DGC is abundant at cell adhesion sites, including the NMJ70. In the absence of a functioning DGC in our DMD models, the second major costameric stabilization unit, the laminin-integrin-talin complex, likely becomes more important. Fhl2 has been shown to bind several integrins in yeast and human HEK293 cells (integrin-α3A, α3B and α7A) and colocalize with integrins in myoblasts and mouse cardiac muscle at focal adhesion sites and the sarcolemma (integrin-α7β1)71,72. Furthermore, lack of integrin-α3 causes nerve terminal detachment in mouse NMJs73. We speculate that Fhl2 may contribute towards stabilizing the sarcolemma, and in extension the NMJs, potentially through interaction with integrins.
Fhl2 is linked to wound healing in several tissues, including muscle51,74,75, and lack of Fhl2 or its suppression in these processes is coupled with chronical inflammation or deteriorating wound healing76,77. In muscle, Fhl2 has been suggested to regulate il-6 and il-8 production via MAPK signaling, playing a role in post-injury inflammation78,79. Our data show increased il-8 levels under uninjured and dmd−/− fhl2b overexpression conditions, supporting these claims. Recently, macrophages have been shown to be critical in muscle injury regeneration by inducing satellite cell proliferation49. Interestingly, our RNA-sequencing data show upregulation of macrophage-specific genes (mmp13a, mpeg1.2) in Tg(503unc:fhl2b-T2A-EGFP and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP larvae. Furthermore, in needle-stick injured Tg(503unc:fhl2b-T2A-EGFP larvae, macrophages are significantly more abundant in the early stage of myofiber repair and regeneration compared to controls. Subsequently, the number of Pax7 positive satellite cells is significantly higher at an early stage of regeneration. As a result, myofibers are restored more rapidly at the wound site. In summary, we propose that overexpression of fhl2b results in accelerated muscle regeneration via enhanced macrophage recruitment.
Taking advantage of the EOMs innate resistance to muscular dystrophies combined with genetic models constitutes a novel approach to identify treatment strategies for these devastating conditions. Our study shows Fhl2 upregulation in EOMs of several disease models, suggesting that Fhl2 is a strong candidate for the development of future therapeutic strategies for a common management of muscular dystrophies.
Methods
Animal husbandry and ethical approval
All experiments were performed in compliance with national and institutional laws and guidelines and the study is reported in accordance with ARRIVE guidelines. All zebrafish animal experiments were ethically approved by the Regional Ethics Committee at the court of Appeal of Northern Norrlands Umeå djurförsöksetiska nämnd, Dnr: A6 2020. The Mouse experiments were approved by the Animal Review Board at the Court of Appeal of Northern Norrland in Umeå Dnr: A22–2023. Human samples were collected with the approval of the Regional Ethical Review Board in Umeå (Dnr: 2010-373-31 M), in accordance with the principles of the Declaration of Helsinki.
Zebrafish larvae (Danio rerio) and adult fish were maintained from AB WT. Mutant lines used were desmaumu10, desmbumu11, obscnbumu16, plecbumu25, fhl2aumu32, fhl2bumu33 and sapjet222a (referred to as dmd+/-). Transgenic lines used were Tg(mylz2:EGFP)i135, Tg(smyhc1:tdTomato)i261, Tg(503unc:EGFP)umu37 and Tg(503unc:fhl2b-T2A-EGFP)umu34. All zebrafish of the same genotype were reared from the same parental couple, to minimize genetic background bias across age. Additionally, WT zebrafish used originated from the same line utilized when generating the desma−/−;desmb−/− mutant. Zebrafish were maintained by standard procedures on a 10/14 h dark/light cycle at 28 °C at the Umeå University Zebrafish Facility and fed twice daily (Live artemia cysts 10309, ZEBCARE B.V, Nederweert, The Netherlands). Leftover mouse tissue from 11 mice that were used in terminal experiments performed by other researchers was kindly donated (Leif Carlsson, Umeå University) and were of mixed WT background (western blot, eight weeks old: 129/Sv:CBA/J:C57BL/6 J:DBA2/J), immunohistochemistry, four weeks old: 129/Sv:CBA/J:C57BL/6 J:DBA2/J). Mice were maintained under 12:12 h light/dark cycles at constant temperature and humidity (22 °C and 50% humidity, fed formula 1310 breeding diet ad libitum (Altromin Spezialfutter GmbH, Lage, Germany)) and were euthanized by cervical dislocation. A total of ten EOM muscle samples were obtained at autopsy from five human donors (four men and one woman, ages 47–80) who, when alive, had consented to donate their eyes and other tissues post-mortem for transplantation and research purposes, according to Swedish law. There was no previously known neuromuscular disease among the donors. Schematic images were adapted from https://www.biorender.com under a subscription license.
Generation of desma, desmb, fhl2a and fhl2b mutant zebrafish using CRISPR/Cas9
desma and desmb zebrafish mutants were generated using methods previously described80. For desma, and desmb a guide RNA (gRNA) targeting exon 1 was synthesized using the sequences ATTCAGCCTCCGCCGAGTCGG and GGTGGGTCGGGCAGCTCTCGG respectively, and was coupled with a gRNA scaffold80. The gRNA was then transcribed using the MegaShortScript T7 (Invitrogen) kit and co-injected with Cas9 protein (New England Biolabs) into one-cell stage zebrafish eggs. Injected larvae were grown to adulthood, outcrossed into WT zebrafish, and screened to identify founders containing germline mutations. Mutant zebrafish larvae carrying a 5 bp and a 20 bp deletion in the desma and desmb genes, respectively, were chosen for further examination. Genotyping of desma mutant larvae was accomplished by standard PCR protocols using forward 5’-ATAGAAGTGGGCGCCAATG-3’ and reverse 5’-GTCTTGAGGAGCCAGAGGAA-3’ primers and desmb mutant larvae were genotyped using forward 5’-AGCCACTCTTATGCCACCTC-3’ and reverse 5’GCGGTCATTTAGATGCTGAAG-3’ primers. The PCR products were then digested overnight using HinfI and AluI for desma and desmb, respectively, and analyzed on a 2% agarose gel. fhl2a and fhl2b mutants were generated as described above, using gRNA sequences AAGAAGTATGTCCTGCGTGAGG and CGGGAAGAAGTACGTCCTGCGG respectively. To genotype fhl2a mutants forward 5’ATAGAAGTGGGCGCCAATG-3’ and reverse 5’TGGGTTTCTTGCATTCCTCG-3’ primers were used and the resulting PCR product was digested with EcoNI to screen for mutations. fhl2b mutants were genotyped using 5’CCCTTTCAACCTGTCTGCAC-3’ and reverse 5’GCAGGTATTGGAGTAGAGGCT-3’ primers and the PCR product was digested using HpyCH4IV. The resulting fhl2a and fhl2b mutants chosen for investigation both carried 7 bp deletions in exon 1.
Transgenesis
An overexpression vector (503unc:fhl2b-T2A-EGFP) containing the muscle-specific promoter 503unc driving expression of fhl2b coupled to EGFP via T2A was acquired from Vectorbuilder (Neu-Isenburg, Germany). Subsequently, AgeI and NheI were used to remove the fhl2b cassette and fused using ligase, to generate a 503unc:EGFP vector. Both vectors were co-injected in one cell-stage AB WT zebrafish eggs with Tol2-transposase RNA, at 30 ng/ul, respectively. Founders were screened using EGFP expression detected under a fluorescent dissecting microscopy and out crossed with wild type fish to generate stable lines used in our experiments.
Immunohistochemistry and TUNEL assay
The muscle samples from both humans and mouse were immediately mounted on cardboard after collection and rapidly frozen in propane chilled with liquid nitrogen and then stored at −80 °C until sectioned. Serial transverse sections, 5–8 µm thick, were cut using a cryostat (Reichert Jung; Leica, Heidelberg, Germany) at a temperature of −23 °C and collected on glass slides which were processed for immunostaining as described previously18. Adult zebrafish at the appropriate age were deeply anesthetized with 0.01% ethyl 3-aminobenzoate methane sulfonate (Tricaine, MS-222, Sigma Aldrich) and sacrificed by decapitation. The heads and posterior trunks were fixed in 2% paraformaldehyde (PFA) for 1 h at RT and stepwise incubated in 10, 20 and 30% sucrose for 12 h each, at 4 °C. The heads and trunks were then mounted separately on cardboard using OCT cryomount and snap frozen in liquid nitrogen chilled propane, and finally serially cut into 12 µm (trunk tissue) or 14 µm (head tissue) thick sections using a cryostat (Reichert Jung; Leica, Heidelberg, Germany). Trunk preparations were cut from adult zebrafish at the anal opening and cut again 8–10 mm caudally to the initial cut. Sections were always made from the proximal end of the specimen to assure the highest possible section similarity between fish. Sectioned zebrafish muscle tissue was rinsed in phosphate buffered saline (PBS) for 15 min before the addition of blocking solution (1% blocking reagent (Roche Diagnostics GmBH, Manheim, Germany) with 0.4% Triton X, 5% dimethyl sulfoxide (DMSO) and 0,1 % TWEEN20) for 1 h, at RT. Primary antibodies were diluted in blocking solution and were applied for 48 h, at 4 °C. Sections were then washed 3 times in PBS and incubated with secondary antibodies for 24 h at 4 °C, washed 3 times in PBS and mounted using 80% glycerol. Adult 12 months old zebrafish hearts were dissected as previously described81 and fixed for 1 h at RT in 2% PFA before stepwise incubation in sucrose as described above. Next, the hearts were orientated in the same position in OCT cryomount before they were carefully placed in −80 °C for storage until sectioned at 20 µm.
For whole-mount immunohistochemistry, larvae were fixed in 2% PFA for 1 h at RT, washed 3 times in PBS, acetone cracked for 1 h at −20 °C, washed 3 times in PBS and incubated for 1 h with blocking solution at RT. The blocking solution was replaced by a blocking solution containing the primary antibody and incubated 48–96 h depending on the stage. Larvae were washed 3 times 15 min in PBS before the addition of secondary antibodies diluted in blocking solution and incubated 24–48 h depending on the stage before being washed 3 times 15 min in PBS again. Lastly, larvae were allowed to equilibrate to 80% glycerol for 1 h before being mounted on glass slides for imaging. DAPI and phalloidin were added with the secondary antibodies. All primary and secondary antibodies used are presented in Supplementary Data 4.
Additionally, a TUNEL assay (Click-iT, Alexa Flour, Invitrogen, Thermo Fisher Scientific) was performed to label nuclei containing degraded DNA molecules, following the manufacturer’s description.
Western blot
Protein extraction was performed in Pierce™ RIPA buffer (Thermo Fisher Scientific, 89901), supplemented with a protease and phosphatase inhibitor cocktail (Thermo Fisher Scientific, 1861282), employing a handheld tissue ruptor. Protein concentration was quantified using Pierce™ BCA Protein Assay Kit (Thermo Fisher Scientific, 23225). An equal amount of proteins were run in Any KD precast polyacrylamide gel (Bio-Rad Laboratories) and transferred onto polyvinylidene fluoride (PVDF) membrane. Primary antibodies used were as follows: FHL2 (Medical and Biological laboratory, K0055-3, 1:1000) for mouse muscle lysate, FHL2 (Atlas Antibodies, HPA005922, 1:1000) for human muscle lysate, and GAPDH (Abcam, AB8245, 1:2000). Anti-rabbit IgG, HRP-linked secondary antibody (Cell signaling, 7074, 1:2000) and anti-mouse IgG HRP-linked secondary antibody (Cell signaling, 7076, 1:2000) were used before membrane underwent incubation with SuperSignal West Pico PLUS Chemiluminescent Substrate (Thermo Fisher Scientific, 34580) and analyzed using an Odyssey Fc Dual-Mode Imaging System (LI-COR Biotechnology). Uncropped blots, including all samples, are displayed in the Source data file for Fig. 3.
Whole-mount in situ hybridization
RNA probes for in situ hybridization were synthesized using a PCR method for RNA probes as previously described82. Zebrafish larvae were fixed in 4% paraformaldehyde (PFA) overnight at 4 oC and dehydrated in 30%, 50%, 70% and 100% methanol and stored in 100% methanol at −20 oC until use. Whole-mount in situ hybridization was performed as described previously83. The RNA probes were based on the sequences of fhl2a (NM_001003732.1) and fhl2b (NM-001006028.2). Primers used to generate probes were: fhl2a forward 5’-CCTGCGTGAGGACAACCCATAC-3’ and reverse 5’-GGTCTCATGCCAGCTGTTTCC-3’. fhl2b forward 5’-GCAAAAAGCCCATTGGCTGC-3’ and reverse 5’-CAGGTCTCATGCCAGCTGTTC-3’.
BrdU treatment
To label the proliferating cell population in our model, zebrafish larvae were collected at 4 dpf treated with BrdU at a concentration of 10 mM and incubated overnight at 28.5 oC. 5 dpf larvae were fixed in 4% PFA for 2 h at room temperature. BrdU was counter-stained with anti BrdU-555 (1:500).
Wound assay
3 dpf larvae were anesthetized in tricane in embryo medium and placed dorsal side up on a 2% agarose gel in a petri dish. Mechanical injuries were targeted to the somites on the dorsal side of embryos, above the cloaca, using a single stab with a sharpened glass capillary. This generated extensive muscle tissue damage localized to approximately one somite. The larvae were then washed in fresh embryo medium and reared until the appropriate stage. More than 95% of larvae survived this procedure.
Imaging and automated measurements
Imaging of sectioned tissue and whole-mount larvae was performed using a Nikon A1 confocal microscope (Nikon, Tokyo, Japan). Images are presented as five merged Z-stacks of representative areas or Z-depth in all figures. Automated measurements of EOM myofiber areas was performed using CellProfiler84. Pre-photographed confocal images of sibling control, dmd−/− and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) 5 dpf larvae immunolabeled for acetylated tubulin and α-bungarotoxin were uploaded to the Imaris Arena (10.1, Bitplane, UK) for quantitative measurements of axons and NMJs. Each image was quantified on four somites directly dorsal to the larval cloaca. To assess acetylated tubulin lengths and volume, images were first individually segmented in the Labkit extension of ImageJ85 to remove background signal. A newly generated layer containing a robust signal was passed back to Imaris Arena. Next, the FilamentTracer tool was utilized to generate 3D filament structures, to quantify axon length (µm). Settings utilized were: Algorithm – Autopath (no loops), shortest distance from distance map, Threshold – Auto. To quantify axon volume (µm3), the Surface tool was applied to the same layer using the following settings: Threshold – Absolute intensity, filter – auto. To quantify normal sized NMJs, the surface tool was applied to the α-bungarotoxin layer using the following settings: Threshold – Absolute intensity. A voxel filter excluding surfaces <37 voxels was applied to all images, determined by the smallest sibling control NMJs detected in 5 samples.
RNA-sequencing and analysis
For each group, 30 five months old and size matched (24 mm ± 3 mm) adult desma−/−;desmb−/−;mylz2:EGFP mutants and WTAB;mylz2:EGFP were euthanized by Tricaine. These were then swiftly subjected to EOM dissection, essentially as described20 but without paraformaldehyde fixation. All six EOMs still attached to the sclera were transferred into RNA-Later solution kept on ice and the EOMs were further cleaned, pooled and frozen in RNA-later. A piece of the lateral trunk including the slow domain myofibers from each fish was also excised and treated in the same manner. For EOM and trunk sequencing of 20 months old zebrafish, the same process was applied except the size-matched animals were 31 mm ± 3 mm and each group contained 12 animals. In total, 254 fish were harvested. For 5 dpf zebrafish larvae, four groups of 10 larvae each containing sibling controls, Tg(503unc:fhl2b-T2A-EGFP), dmd−/− and dmd−/−;Tg(503unc:fhl2b-T2A-EGFP) larvae were cut diagonally distal to the swim bladder to exclude the majority of the gastrointestinal region and the head from the analysis. These larvae were then stored in RNA-later at −80 °C until further use.
Pooled tissue was thawed and briefly rinsed in PBS before RNA extraction was performed using the TRIzol (Invitrogen, Thermo Fisher Scientific) reagent standard procedures and isolated in 15 µl of RNAse free H2O. Quality controls were performed using a bioanalyzer and a RIN value greater than 8 was considered acceptable for analysis. Library preparation was performed using the TruSeq Stranded mRNA Library Prep kit (cat#20020595, Illumina, Inc.), including poly-A selection, according to the manufacturer’s instructions. Unique index adapters were used (cat#20022371, Illumina, Inc. 15 cycles of amplification). The RNA-seq libraries were sequenced on a Novaseq 6000 Sequencing system (Illumina, Inc.) obtaining in average ~78.5 million 150 paired end reads per library. Library preparation and sequencing was performed at SciLifeLab Stockholm.
Fastq files were quality controlled using FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and raw reads were mapped to the zebrafish genome (GRCz11) using STAR (2.7.6a)86. Normalization and differential expression analysis were performed using DESeq287, only genes with 10 or more reads were processed. Differentially expressed genes (DEGs, padj = <0.05) from different comparisons were intersected using the UpSet plot function from the ComplexHeatmap package88 and Gene Ontology (GO) analysis was performed using the clusterProfiler package89 (padj = <0.05) and plotted using the enrichplot package. Expression of single genes was plotted using ggplot. All statistical analysis related to RNA-sequencing data was performed within the R environment (version 4.2.3) using basic built-in functions and publicly available packages listed above. These are open-source tools and can be accessed via Bioconductor (hppts://Bioconductor.org/).
qPCR
Quantitative PCR was performed on using the same RNA-extraction method as described above on whole 5 dpf larvae. cDNA was synthesized using SuperScript IV (Invitrogen, Thermo Fisher Scientific). Primers used for all genes are presented in Supplementary Data 5. β-actin was used as a reference gene. The samples were run using an Applied Biosystems VIIA-7 Real Time PCR system (Thermo Fisher Scientific) using FastStart universal SYBR green master mix (Roche).
Physiological properties of zebrafish muscle force/tension relationship
Zebrafish larvae were examined with length- force experiments as previously described by Dou et al90. and Li et. al.91. In brief, 5 and 6 dpf WT and desma−/−;desmb−/− larvae were euthanized and mounted with aluminum clips between a force transducer and a puller for length adjustment. The bath was perfused at 22oC in a MOPS buffered physiological solution with a pH of 7.4. The preparations were allowed to acclimatize in the solution for at least 10 min before initiating contractions using single twitch stimuli with 0.5 ms pulses at 2-min intervals and supramaximal voltage, via platinum electrodes placed on each side of the larvae. Measurements were initiated at slack length (Ls) and length was gradually increased by 10% steps every 4 min, between stimuli, until a decline in active force was observed. At each length, both active and passive force were recorded twice, and an average was calculated. The values were plotted against relative stretch (lambda=length/Ls).
Spontaneous swimming and resistance swimming
5 dpf WT and desma−/−;desmb−/− zebrafish larvae were placed in a 48 well-plate inside the Viewpoint ZebraBox system (Viewpoint Behavior Technology) to determine spontaneous movement patterns. Larvae were carefully touched by the tail to ensure mobility and were subsequently allowed to acclimatize to the environment for 15 min before a 60 min protocol with a steady light was initiated. The movement thresholds were set to inactivity = 0 and large movements = 1. Inactivity, small movement and large movement counts, swimming distance and swimming duration were recorded. For resistance swimming, 4 dpf larvae were placed in 1% methyl cellulose in 1x E3 medium and reared at 28.5 °C overnight. The following day, larvae were analyzed using fluorescent microscopy.
Statistical analysis
Myofiber and myonuclei counts were performed in the whole slow domain of adult zebrafish trunks except when F310 positive myofibers were counted outside of the slow domain. F310 positive myofibers were then counted in 2 × 0.6 mm squares (20x magnification) just medial to the intermediate fast domain. The myofibers of each fish was counted twice, once for each side of the fish for all quantifications except for F310, where a total of four squares were counted, two for each side of the fish. For cardiac dimension measurements, the major and minor axis of the heart were added together and divided by the size of individual fish before comparisons between genotypes. The myofibers and myonuclei in the EOMs were counted separately in a total of two per fish. The medial rectus muscle was used. In embryonic experiments including Pax7, BrdU or TUNEL positive cells in the trunk muscle, somite numbers 13–22 was counted and the total number of positive cells were divided by the numbers of somites. For quantifications of α-bungarotoxin positive NMJs somite number 17–20 were counted. For quantifications of SV2/ α-bungarotoxin double positive NMJs, somites 18-19 were counted. Acetylated tubulin volume and length were measured for somites 17–20. For broken myofiber counts, somites 8–24 were counted. The number of larvae in the experiments are presented in the graphs. Measurements were taken from distinct samples and not from repeated measurements of the same samples. Normality tests were performed prior to the statistical analysis for each experiment.
All data was collected in Microsoft Excel and plotted in GraphPad Prism 10.0. Statistical analysis was performed using two-sided t-tests with Welch correction, p = <0.05 was considered significant (*p < 0.05, **p < 0.005, ***p < 0.0005, ****p < 0.0001). All data are presented as mean ± standard error of mean (SEM). Kaplan-Meier log rank test was used to determine the difference between genotypes in the survival analysis. p = <0.05 was considered significant.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The RNA-sequencing data generated in this study have been deposited in the Gene Expression Omnibus (GEO) database under accession code GSE242137. The processed RNA-sequencing data generated in this study are provided in the Supplementary Data 1, 2, and 3. Source data are provided with this paper.
References
Benarroch, L., Bonne, G., Rivier, F. & Hamroun, D. The 2021 version of the gene table of neuromuscular disorders (nuclear genome). Neuromuscul. Disord. 30, 1008–1048 (2020).
Theadom, A. et al. Prevalence of muscular dystrophies: a systematic literature review. Neuroepidemiology 43, 259–268 (2014).
Ankala, A. et al. Aberrant firing of replication origins potentially explains intragenic nonrecurrent rearrangements within genes, including the human DMD gene. Genome Res 22, 25–34 (2012).
Emery, A. E. The muscular dystrophies. Lancet 359, 687–695 (2002).
Campbell, K. P. & Kahl, S. D. Association of dystrophin and an integral membrane glycoprotein. Nature 338, 259–262 (1989).
Eid Mutlak, Y. et al. A signaling hub of insulin receptor, dystrophin glycoprotein complex and plakoglobin regulates muscle size. Nat. Commun. 11, 1381 (2020).
Hynes, R. O. Integrins: bidirectional, allosteric signaling machines. Cell 110, 673–687 (2002).
Pasternak, C., Wong, S. & Elson, E. L. Mechanical function of dystrophin in muscle cells. J. Cell Biol. 128, 355–361 (1995).
Bassett, D. I. et al. Dystrophin is required for the formation of stable muscle attachments in the zebrafish embryo. Development 130, 5851–5860 (2003).
Granato, M. et al. Genes controlling and mediating locomotion behavior of the zebrafish embryo and larva. Development 123, 399–413 (1996).
Stocco, A. et al. Treatment with a triazole inhibitor of the mitochondrial permeability transition pore fully corrects the pathology of sapje zebrafish lacking dystrophin. Pharm. Res 165, 105421 (2021).
Kaminski, H. J., al-Hakim, M., Leigh, R. J., Katirji, M. B. & Ruff, R. L. Extraocular muscles are spared in advanced Duchenne dystrophy. Ann. Neurol. 32, 586–588 (1992).
Formicola, L., Marazzi, G. & Sassoon, D. A. The extraocular muscle stem cell niche is resistant to ageing and disease. Front Aging Neurosci. 6, 328 (2014).
Li, A. et al. Distinct transcriptomic profile of satellite cells contributes to preservation of neuromuscular junctions in extraocular muscles of ALS mice. Preprint at bioRxiv, https://doi.org/10.1101/2023.02.12.528218 (2023).
Khurana, T. S. et al. Absence of extraocular muscle pathology in Duchenne’s muscular dystrophy: role for calcium homeostasis in extraocular muscle sparing. J. Exp. Med 182, 467–475 (1995).
Fischer, M. D. et al. Definition of the unique human extraocular muscle allotype by expression profiling. Physiol. Genomics 22, 283–291 (2005).
Gargan, S. et al. Mass Spectrometric Profiling of Extraocular Muscle and Proteomic Adaptations in the mdx-4cv Model of Duchenne Muscular Dystrophy. Life (Basel) 11, https://doi.org/10.3390/life11070595 (2021).
Liu, J. X. & Pedrosa Domellof, F. Complex Correlations Between Desmin Content, Myofiber Types, and Innervation Patterns in the Human Extraocular Muscles. Invest Ophthalmol. Vis. Sci. 61, 15 (2020).
Liu, J. X. & Pedrosa Domellof, F. Cytoskeletal Proteins in Myotendinous Junctions of Human Extraocular Muscles. Invest Ophthalmol. Vis. Sci. 62, 19 (2021).
Dennhag, N., Liu, J. X., Nord, H., von Hofsten, J. & Pedrosa Domellof, F. Absence of Desmin in Myofibers of the Zebrafish Extraocular Muscles. Transl. Vis. Sci. Technol. 9, 1 (2020).
Janbaz, A. H., Lindstrom, M., Liu, J. X. & Pedrosa Domellof, F. Intermediate filaments in the human extraocular muscles. Invest Ophthalmol. Vis. Sci. 55, 5151–5159, (2014).
Price, M. G. Molecular analysis of intermediate filament cytoskeleton–a putative load-bearing structure. Am. J. Physiol. 246, H566–572, (1984).
Fuchs, C. et al. Desmin enters the nucleus of cardiac stem cells and modulates Nkx2.5 expression by participating in transcription factor complexes that interact with the nkx2.5 gene. Biol. Open 5, 140–153 (2016).
Lockard, V. G. & Bloom, S. Trans-cellular desmin-lamin B intermediate filament network in cardiac myocytes. J. Mol. Cell Cardiol. 25, 303–309 (1993).
Chapman, M. A. et al. Disruption of both nesprin 1 and desmin results in nuclear anchorage defects and fibrosis in skeletal muscle. Hum. Mol. Genet 23, 5879–5892 (2014).
Carlsson, L. & Thornell, L. E. Desmin-related myopathies in mice and man. Acta Physiol. Scand. 171, 341–348 (2001).
Henderson, C. A., Gomez, C. G., Novak, S. M., Mi-Mi, L. & Gregorio, C. C. Overview of the Muscle Cytoskeleton. Compr. Physiol. 7, 891–944 (2017).
Rodriguez, M. A., Liu, J. X., Parkkonen, K., Li, Z. & Pedrosa Domellof, F. The Cytoskeleton in the Extraocular Muscles of Desmin Knockout Mice. Invest Ophthalmol. Vis. Sci. 59, 4847–4855, (2018).
Goldfarb, L. G., Olive, M., Vicart, P. & Goebel, H. H. Intermediate filament diseases: desminopathy. Adv. Exp. Med Biol. 642, 131–164 (2008).
Maddison, P. et al. Clinical and myopathological characteristics of desminopathy caused by a mutation in desmin tail domain. Eur. Neurol. 68, 279–286 (2012).
Kayman Kurekci, G. et al. Knockout of zebrafish desmin genes does not cause skeletal muscle degeneration but alters calcium flux. Sci. Rep. 11, 7505 (2021).
Hardee, J. P. et al. Metabolic remodeling of dystrophic skeletal muscle reveals biological roles for dystrophin and utrophin in adaptation and plasticity. Mol. Metab. 45, 101157 (2021).
Molotsky, E., Liu, Y., Lieberman, A. P. & Merry, D. E. Neuromuscular junction pathology is correlated with differential motor unit vulnerability in spinal and bulbar muscular atrophy. Acta Neuropathol. Commun. 10, 97 (2022).
Elworthy, S., Hargrave, M., Knight, R., Mebus, K. & Ingham, P. W. Expression of multiple slow myosin heavy chain genes reveals a diversity of zebrafish slow twitch muscle fibres with differing requirements for Hedgehog and Prdm1 activity. Development 135, 2115–2126 (2008).
Nguyen, P. D. et al. Muscle Stem Cells Undergo Extensive Clonal Drift during Tissue Growth via Meox1-Mediated Induction of G2 Cell-Cycle Arrest. Cell Stem Cell 21, 107–119 e106 (2017).
Genini, M. et al. Subtractive cloning and characterization of DRAL, a novel LIM-domain protein down-regulated in rhabdomyosarcoma. DNA Cell Biol. 16, 433–442 (1997).
Zrelski, M. M., Kustermann, M. & Winter, L. Muscle-Related Plectinopathies. Cells 10, https://doi.org/10.3390/cells10092480 (2021).
Randazzo, D. et al. Exercise-induced alterations and loss of sarcomeric M-line organization in the diaphragm muscle of obscurin knockout mice. Am. J. Physiol. Cell Physiol. 312, C16–C28 (2017).
Kong, Y. et al. Cardiac-specific LIM protein FHL2 modifies the hypertrophic response to beta-adrenergic stimulation. Circulation 103, 2731–2738 (2001).
Chu, P. H., Bardwell, W. M., Gu, Y., Ross, J. Jr. & Chen, J. FHL2 (SLIM3) is not essential for cardiac development and function. Mol. Cell Biol. 20, 7460–7462 (2000).
Duan, Y. et al. Deletion of FHL2 in fibroblasts attenuates fibroblasts activation and kidney fibrosis via restraining TGF-beta1-induced Wnt/beta-catenin signaling. J. Mol. Med (Berl.) 98, 291–307 (2020).
Hamidouche, Z. et al. FHL2 mediates dexamethasone-induced mesenchymal cell differentiation into osteoblasts by activating Wnt/beta-catenin signaling-dependent Runx2 expression. FASEB J. 22, 3813–3822 (2008).
Wei, Y. et al. Identification of the LIM protein FHL2 as a coactivator of beta-catenin. J. Biol. Chem. 278, 5188–5194 (2003).
Martin, B. et al. The LIM-only protein FHL2 interacts with beta-catenin and promotes differentiation of mouse myoblasts. J. Cell Biol. 159, 113–122 (2002).
Berger, J. & Currie, P. D. 503unc, a small and muscle-specific zebrafish promoter. Genesis 51, 443–447 (2013).
Bradford, Y. M. et al. Zebrafish Information Network, the knowledgebase for Danio rerio research. Genetics 220, https://doi.org/10.1093/genetics/iyac016 (2022).
Ma, G. et al. MiR-206, a key modulator of skeletal muscle development and disease. Int J. Biol. Sci. 11, 345–352 (2015).
Sandonà, M. et al. Histone Deacetylases: Molecular Mechanisms and Therapeutic Implications for Muscular Dystrophies. Int. J. Mol. Sci. 24, https://doi.org/10.3390/ijms24054306 (2023).
Ratnayake, D. et al. Macrophages provide a transient muscle stem cell niche via NAMPT secretion. Nature 591, 281–287 (2021).
Nakazawa, N., Sathe, A. R., Shivashankar, G. V. & Sheetz, M. P. Matrix mechanics controls FHL2 movement to the nucleus to activate p21 expression. Proc. Natl Acad. Sci. USA 113, E6813–E6822 (2016).
Shi, X., Bowlin, K. M. & Garry, D. J. Fhl2 interacts with Foxk1 and corepresses Foxo4 activity in myogenic progenitors. Stem Cells 28, 462–469 (2010).
Zhu, Y. et al. miR-377 Inhibits Proliferation and Differentiation of Bovine Skeletal Muscle Satellite Cells by Targeting FHL2. Genes (Basel) 13, https://doi.org/10.3390/genes13060947 (2022).
Hojayev, B., Rothermel, B. A., Gillette, T. G. & Hill, J. A. FHL2 binds calcineurin and represses pathological cardiac growth. Mol. Cell Biol. 32, 4025–4034 (2012).
Okamoto, R. et al. FHL2 prevents cardiac hypertrophy in mice with cardiac-specific deletion of ROCK2. FASEB J. 27, 1439–1449 (2013).
Purcell, N. H. et al. Extracellular signal-regulated kinase 2 interacts with and is negatively regulated by the LIM-only protein FHL2 in cardiomyocytes. Mol. Cell Biol. 24, 1081–1095 (2004).
Methatham, T., Tomida, S., Kimura, N., Imai, Y. & Aizawa, K. Inhibition of the canonical Wnt signaling pathway by a β-catenin/CBP inhibitor prevents heart failure by ameliorating cardiac hypertrophy and fibrosis. Sci. Rep. 11, 14886 (2021).
Armstrong, D. D. & Esser, K. A. Wnt/beta-catenin signaling activates growth-control genes during overload-induced skeletal muscle hypertrophy. Am. J. Physiol. Cell Physiol. 289, C853–859, (2005).
Brun, J. et al. The LIM-only protein FHL2 controls mesenchymal cell osteogenic differentiation and bone formation through Wnt5a and Wnt10b. Bone 53, 6–12 (2013).
Lange, S. et al. Subcellular targeting of metabolic enzymes to titin in heart muscle may be mediated by DRAL/FHL-2. J. Cell Sci. 115, 4925–4936 (2002).
van der Pijl, R. J. et al. The titin N2B and N2A regions: biomechanical and metabolic signaling hubs in cross-striated muscles. Biophys. Rev. 13, 653–677 (2021).
Sun, X. et al. Mechanosensing through Direct Binding of Tensed F-Actin by LIM Domains. Dev. Cell 55, 468–482 e467 (2020).
Wang, Y. X., Wang, D. Y., Guo, Y. C. & Guo, J. Zyxin: a mechanotransductor to regulate gene expression. Eur. Rev. Med Pharm. Sci. 23, 413–425 (2019).
Fisher, L. A. B. & Schock, F. The unexpected versatility of ALP/Enigma family proteins. Front Cell Dev. Biol. 10, 963608 (2022).
Pathak, P. et al. Myopathy associated LDB3 mutation causes Z-disc disassembly and protein aggregation through PKCalpha and TSC2-mTOR downregulation. Commun. Biol. 4, 355 (2021).
Shiao, T. et al. Defects in neuromuscular junction structure in dystrophic muscle are corrected by expression of a NOS transgene in dystrophin-deficient muscles, but not in muscles lacking alpha- and beta1-syntrophins. Hum. Mol. Genet 13, 1873–1884 (2004).
Dickson, B. J. Molecular mechanisms of axon guidance. Science 298, 1959–1964 (2002).
Correia, J. C. et al. Muscle-secreted neurturin couples myofiber oxidative metabolism and slow motor neuron identity. Cell Metab. 33, 2215–2230 e2218 (2021).
Borisov, A. B., Dedkov, E. I. & Carlson, B. M. Interrelations of myogenic response, progressive atrophy of muscle fibers, and cell death in denervated skeletal muscle. Anat. Rec. 264, 203–218 (2001).
Yamauchi, N. et al. High-intensity interval training in the form of isometric contraction improves fatigue resistance in dystrophin-deficient muscle. J. Physiol., https://doi.org/10.1113/JP284532 (2023).
Belhasan, D. C. & Akaaboune, M. The role of the dystrophin glycoprotein complex on the neuromuscular system. Neurosci. Lett. 722, 134833 (2020).
Wixler, V. et al. The LIM-only protein DRAL/FHL2 binds to the cytoplasmic domain of several alpha and beta integrin chains and is recruited to adhesion complexes. J. Biol. Chem. 275, 33669–33678 (2000).
Samson, T. et al. The LIM-only proteins FHL2 and FHL3 interact with alpha- and beta-subunits of the muscle alpha7beta1 integrin receptor. J. Biol. Chem. 279, 28641–28652 (2004).
Ross, J. A. et al. Multiple roles of integrin-alpha3 at the neuromuscular junction. J. Cell Sci. 130, 1772–1784 (2017).
Wixler, V. et al. Deficiency in the LIM-only protein Fhl2 impairs skin wound healing. J. Cell Biol. 177, 163–172 (2007).
Alnajar, A. et al. The LIM-only protein FHL2 attenuates lung inflammation during bleomycin-induced fibrosis. PLoS One 8, e81356 (2013).
Wixler, V. et al. FHL2 regulates the resolution of tissue damage in chronic inflammatory arthritis. Ann. Rheum. Dis. 74, 2216–2223 (2015).
van de Pol, V. et al. LIM-only protein FHL2 attenuates inflammation in vascular smooth muscle cells through inhibition of the NFkappaB pathway. Vasc. Pharm. 125-126, 106634 (2020).
Wong, C. H., Mak, G. W., Li, M. S. & Tsui, S. K. The LIM-only protein FHL2 regulates interleukin-6 expression through p38 MAPK mediated NF-kappaB pathway in muscle cells. Cytokine 59, 286–293 (2012).
Wan, S. et al. Expression of FHL2 and cytokine messenger RNAs in human myocardium after cardiopulmonary bypass. Int J. Cardiol. 86, 265–272 (2002).
Varshney, G. K. et al. A high-throughput functional genomics workflow based on CRISPR/Cas9-mediated targeted mutagenesis in zebrafish. Nat. Protoc. 11, 2357–2375 (2016).
Arnaout, R., Reischauer, S. & Stainier, D. Y. Recovery of adult zebrafish hearts for high-throughput applications. J. Vis. Exp., https://doi.org/10.3791/52248 (2014).
Hua, R., Yu, S., Liu, M. & Li, H. A PCR-based method for RNA probes and applications in neuroscience. Front. Neurosci. 12, 266 (2018).
Thisse, C. & Thisse, B. High-resolution in situ hybridization to whole-mount zebrafish embryos. Nat. Protoc. 3, 59–69 (2008).
Stirling, D. R. et al. CellProfiler 4: improvements in speed, utility and usability. BMC Bioinforma. 22, 433 (2021).
Arzt, M. et al. LABKIT: Labeling and Segmentation Toolkit for Big Image Data. Front. Comp Sci-Switz 4, ARTN 777728 https://doi.org/10.3389/fcomp.2022.777728 (2022).
Dobin, A. et al. STAR: ultrafast universal RNA-seq aligner. Bioinformatics 29, 15–21 (2013).
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 550 (2014).
Gu, Z., Eils, R. & Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 32, 2847–2849 (2016).
Wu, T. et al. clusterProfiler 4.0: a universal enrichment tool for interpreting omics data. Innov. (Camb.) 2, 100141 (2021).
Dou, Y., Andersson-Lendahl, M. & Arner, A. Structure and function of skeletal muscle in zebrafish early larvae. J. Gen. Physiol. 131, 445–453 (2008).
Li, M., Andersson-Lendahl, M., Sejersen, T. & Arner, A. Knockdown of desmin in zebrafish larvae affects interfilament spacing and mechanical properties of skeletal muscle. J. Gen. Physiol. 141, 335–345 (2013).
Acknowledgements
The authors acknowledge support from the National Genomics Infrastructure in Stockholm and Uppsala funded by Science for Life Laboratory, the Knut and Alice Wallenberg Foundation and the Swedish Research Council, and SNIC/Uppsala Multidisciplinary Center for Advanced Computational Science for assistance with massively parallel sequencing and access to the UPPMAX computational infrastructure. The computations were enabled by resources in projects SNIC 2022/22-442 and NAISS 2023/22-287 provided by the Swedish National Infrastructure for Computing (SNIC) and the National Academic Infrastructure for Supercomputing in Sweden (NAISS) at UPPMAX. The authors acknowledge Philip W. Ingham for generously sharing the Tg(smyhc1:tdTomato)i261 and Tg(mylz2:EGFP)i135 zebrafish lines. The authors thank Leif Carlsson for the donation of mouse tissue. AA was funded by grants from Hans-Gabriel and Alice Trolle-Wachtmeister’s Foundation for Medical Research. JvH was funded by biotechnology grant for basic science FS 2.1.6-1911-22, the Medical Faculty, Umeå University. FPD was supported by research grants from the Swedish Research Council (Dnr 2018-02401), Västerbotten County Council (Central ALF and Spjutspetsmedel), Kronprinsessan Margaretas Arbetsnämnd för synskadade (Stiftelsen KMA), Ögonfonden. ND was supported by Kempestiftelserna, Ögonfonden and Arnerska forskningsfonden.
Funding
Open access funding provided by Umea University..
Author information
Authors and Affiliations
Contributions
N.D., Jv.H., and F.P.D. designed the study. Experiments were performed, analyzed, and interpreted by N.D., A.K., I.N., H.N., M.C., J.-X.L., J.L., A.A., and L.J.B. Transcriptional data was analyzed and interpreted by ND and IN, with supervision from S.R. N.D. wrote the manuscript with supervision and input from S.R., J.v.H. and F.P.D. N.D., Jv.H. and F.P.D. secured funding for the project.
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Dennhag, N., Kahsay, A., Nissen, I. et al. fhl2b mediates extraocular muscle protection in zebrafish models of muscular dystrophies and its ectopic expression ameliorates affected body muscles. Nat Commun 15, 1950 (2024). https://doi.org/10.1038/s41467-024-46187-x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-46187-x
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.