# Exercise, disease state and sex influence the beneficial effects of Fn14-depletion on survival and muscle pathology in the SOD1G93A amyotrophic lateral sclerosis (ALS) mouse model

Gareth Hazell, Eve Mccallion, Nina Ahlskog, Emma R Sutton, Magnus Okoh, Emad I H Shaqoura, Joseph M Hoolachan, Taylor Scaife, Sara Iqbal, Amarjit Bhomra, et al.

► To cite this version:

Gareth Hazell, Eve Mccallion, Nina Ahlskog, Emma R Sutton, Magnus Okoh, et al. Exercise, disease state and sex influence the beneficial effects of Fn14-depletion on survival and muscle pathology in the SOD1G93A amyotrophic lateral sclerosis (ALS) mouse model. Skeletal Muscle, 2024, 14 (1), pp.23. 10.1186/s13395-024-00356-0 . hal-04767133

HAL Id: hal-04767133

https://hal.science/hal-04767133v1

Submitted on 5 Nov 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

Check for
updates

# Exercise, disease state and sex influence the beneficial effects of Fn14-depletion on survival and muscle pathology in the SOD1<sup>G93A</sup> amyotrophic lateral sclerosis (ALS) mouse model

Gareth Hazell¹, Eve McCallion², Nina Ahlskog¹,³, Emma R. Sutton², Magnus Okoh², Emad I. H. Shaqoura², Joseph M. Hoolachan², Taylor Scaife⁴, Sara Iqbal⁴, Amarjit Bhomra¹,³, Anna J. Kordala¹,³, Frederique Scamps⁵, Cedric Raoul⁵,⁶, Matthew J. A. Wood¹,³ and Melissa Bowerman¹,²,⁷*

**Abstract**

**Background** Amyotrophic lateral sclerosis (ALS) is a devastating and incurable neurodegenerative disease. Accumulating evidence strongly suggests that intrinsic muscle defects exist and contribute to disease progression, including imbalances in whole-body metabolic homeostasis. We have previously reported that tumour necrosis factor (TNF)-like weak inducer of apoptosis (TWEAK) and fibroblast growth factor inducible 14 (Fn14) are significantly upregulated in skeletal muscle of the SOD1<sup>G93A</sup> ALS mouse model. While antagonising TWEAK did not impact survival, we did observe positive effects in skeletal muscle. Given that Fn14 has been proposed as the main effector of the TWEAK/Fn14 activity and that Fn14 can act independently from TWEAK in muscle, we suggest that manipulating Fn14 instead of TWEAK in the SOD1<sup>G93A</sup> ALS mice could lead to differential and potentially improved benefits.

**Methods** We thus investigated the contribution of Fn14 to disease phenotypes in the SOD1<sup>G93A</sup> ALS mice. To do so, Fn14 knockout mice (Fn14<sup>-/-</sup>) were crossed onto the SOD1<sup>G93A</sup> background to generate SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> mice. Investigations were performed on both unexercised and exercised (rotarod and/or grid test) animals (wild type (WT), Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup>).

**Results** Here, we firstly confirm that the TWEAK/Fn14 pathway is dysregulated in skeletal muscle of SOD1<sup>G93A</sup> mice. We then show that Fn14-depleted SOD1<sup>G93A</sup> mice display increased lifespan, myofiber size, neuromuscular junction endplate area as well as altered expression of known molecular effectors of the TWEAK/Fn14 pathway, without an impact on motor function. Importantly, we also observe a complex interaction between exercise (rotarod and grid test), genotype, disease state and sex that influences the overall effects of Fn14 deletion on survival, expression of known molecular effectors of the TWEAK/Fn14 pathway, expression of myosin heavy chain isoforms and myofiber size.

<footnote>*Correspondence:
Melissa Bowerman
m.bowerman@keele.ac.uk
Full list of author information is available at the end of the article</footnote>
**Conclusions** Our study provides further insights on the different roles of the TWEAK/Fn14 pathway in pathologi-cal skeletal muscle and how they can be influenced by age, disease, sex and exercise. This is particularly relevant in the ALS field, where combinatorial therapies that include exercise regimens are currently being explored. As such, a better understanding and consideration of the interactions between treatments, muscle metabolism, sex and exer-cise will be of importance in future studies.

**Keywords** Amyotrophic lateral sclerosis, Skeletal muscle, TWEAK, Fn14, Exercise, Sex, Metabolism

**Background**

Amyotrophic lateral sclerosis (ALS) is a devastating and currently incurable neurodegenerative disease. Once symptomatic, the median survival of patients is usually between 3 and 5 years. Clinical manifestations typically occur in mid-life, followed by the rapid and progressive wasting of muscles and subsequent paralysis [1]. ALS can be sporadic (~80%) or familial (~20%) [2], and in the lat-ter case can be caused by numerous genetic mutations with the most common being in *chromosome 9 open reading frame 72 (C9ORF72)* [3, 4], *superoxide dismutase 1 (SOD1)* [5], *Fused in Sarcoma (FUS)* [6, 7] and *TAR DNA-binding protein 43 (TDP-43)* [8–10]. Both sporadic and familial ALS patients present similar symptoms and pathophysiology. While the primary pathological tar-get of ALS is undeniably the motor neurons (both upper and lower), accumulating evidence strongly suggests that intrinsic muscle defects exist and contribute to disease progression and presentation [11]. Indeed, the muscle-restricted expression of mutant SOD1 results in a canon-ical ALS pathophysiology [12, 13]. Furthermore, aberrant genetic, biochemical, developmental, regulatory and physiological changes prior to, or accompanying, motor neuron loss are observed in ALS muscle and progenitor cells [11]. As muscle plays an important role in maintain-ing systemic energy homeostasis [14], intrinsic muscle defects can have severe consequences on whole-body metabolic homeostasis. Interestingly, instances of insulin resistance [15], hyperlipidemia [16], hyperglycemia [17], aberrant fatty acid metabolism [18], hyperglucagonemia [19], glucose intolerance [18] and development of diabe-tes [20] have all been reported in ALS patients and animal models. Furthermore both dietary and exercise interven-tions, which are direct modulators of muscle metabolism [21], have been demonstrated to impact disease progres-sion in ALS patients and animal models [22–24]. Thus, uncovering and targeting pathological molecular effec-tors in ALS muscle may lead to tissue-specific and whole-body improvements [11, 25, 26].

One important pathway that contributes to skeletal muscle health, function and metabolism is controlled by the binding of the tumour necrosis factor (TNF)-like weak inducer of apoptosis (TWEAK) ligand to the TNF fibroblast growth factor inducible 14 (Fn14) receptor [27,

28]. Interestingly, the TWEAK/Fn14 pathway can impact muscle positively or negatively depending on the levels of TWEAK present. High levels are typically associated with detrimental effects while low levels have a benefi-cial impact [27, 28]. Similarly, Fn14 expression is typically very low in healthy muscle and becomes upregulated in muscle atrophy conditions, which can lead to sustained muscle pathology if not restored to normal levels [27, 28]. Furthermore, TWEAK and Fn14 have both been impli-cated in the regulation of key muscle metabolic effec-tors such as peroxisome proliferative activated receptor, gamma, coactivator 1 alpha (PGC-1α), Slc2a4 solute car-rier family 2, member 4 (GLUT4), hexokinase 2 (HKII) and Krüppel-like transcription factor 15 (KLF15) [29].

What still remains unclear however, is the potential role of the TWEAK/Fn14 pathway in neuromuscular conditions, where chronic muscle wasting occurs due to motor neuron loss and muscle denervation [30]. In an attempt to explore this further, we have previously inves-tigated the TWEAK/Fn14 signalling cascade in mouse models of ALS and spinal muscular atrophy (SMA), a childhood neuromuscular disease [31]. In pre-weaned SMA mice, we observed a significant downregulation of *Tweak* and *Fn14* in various skeletal muscles dur-ing disease progression, accompanied by the expected dysregulation of *PGC-1α*, *Glut4*, *HKII* and *Klf15* [32]. Interestingly, administering Fc-TWEAK, an agonist of the pathway, to SMA mice, improved several canonical disease phenotypes [32]. Conversely, we have previously observed that *Tweak* and *Fn14* are significantly upregu-lated in skeletal muscle of *SOD1*<sup>G93A</sup> ALS mice during disease progression [33]. While antagonising TWEAK, either genetically or pharmacologically, did not impact survival, we did observe positive effects in skeletal muscle [33]. Since the receptor has been proposed as the main effector of the TWEAK/Fn14 pathway activity [34] and that Fn14 can act independently from TWEAK in mus-cle [35], it is possible that manipulating Fn14 instead of TWEAK in the *SOD1*<sup>G93A</sup> ALS mice could lead to differ-ential and/or improved benefits.

In this study, we investigated the effect of Fn14 depletion on disease progression and muscle pathol-ogy in *SOD1*<sup>G93A</sup> ALS mice by crossing Fn14 knock-out mice (*Fn14*<sup>-/-</sup>) with the *SOD1*<sup>G93A</sup> mouse model.
We confirmed that the TWEAK/Fn14 pathway is dysregulated in the skeletal muscle of *SOD1*<sup>G93A</sup> mice. We then showed that Fn14-depleted *SOD1*<sup>G93A</sup> mice had an increased lifespan and decreased muscle pathology, which was dependent on exposure to exercise and sex. Our study provides further insights into the different roles of the TWEAK/Fn14 pathway in skeletal muscle and how they may be influenced by age, disease, sex and exercise.

# Methods

## Animals and animal procedures

*SOD1*<sup>G93A</sup> mice (B6.Cg-Tg(SOD1\*G93A)1Gur/J) were obtained from Jackson Laboratories (Strain #: 004435). The *Fn14*<sup>-/-</sup> mouse model [36] was provided by Biogen. Both strains were on a C57BL/6 genetic background.

Experimental procedures with live animals were authorized and approved by the University of Oxford ethics committee and UK Home Office (Project licenses PDFEDC6F0 and 30/2907) in accordance with the Animals (Scientific Procedures) Act 1986.

For survival studies, mice were weighed and monitored daily and culled upon reaching their defined humane endpoint as specified in the project license.

For all experiments, litters were randomly assigned treatment at birth. Sample sizes were determined based on similar studies with *SOD1*<sup>G93A</sup> mice.

For the grid test, we used our previously described protocol [33], whereby starting with a 40 g metal grid (followed by 30, 20 and 10 g grids), we measured the time (maximum 30 s) the animal held on to the grid before dropping it. The experiment was repeated three times with each grid. Muscle strength (arbitrary units) was quantified with the following formula: $(40 \text{ g} \times \text{best time}) + (30 \text{ g} \times \text{best time}) + (20 \text{ g} \times \text{best time}) + (10 \text{ g} \times \text{best time})$.

For the rotarod test, we followed the previously described protocol [37], whereby mice were placed on the rotarod (opposite orientation to rotation) with an acceleration protocol of 4 to 40 rpm in 300 s. The latency to fall (s) and highest rpm reached was recorded.

To reduce the total number of mice used, the fasttwitch *tibialis anterior* (TA) and gastrocnemius muscles from the same mice were used for molecular and histological analyses, respectively.

## qPCRs

RNA was extracted from tissues with the RNeasy kit (Qiagen) or with a Isolate II RNA Mini Kit (Bioline) as per the manufacturers' instructions. The same RNA extraction method was employed for similar experiments and equal RNA amounts were used between samples within the same experiments. cDNA was prepared

with the High-capacity cDNA Kit (Life Technologies) or qPCRBIO cDNA Synthesis Kit (PBCR Biosystems) according to the manufacturers' instructions. The same reverse transcription method was employed for similar experiments. The cDNA template was amplified on a StepOnePlus Real-Time PCR Thermocycler (Life Technologies) with SYBR Green Mastermix (Applied Biosystems) or with qPCRBIO SyGreen Blue Mix Hi-ROX (PCR Biosystems). The same amplification method was used for similar experiments. qPCR data was analysed using the StepOne Software v2.3 (Applied Biosystems). Primers used for qPCR were obtained from IDT and sequences for primers were self-designed (Supplementary Table 1). Relative gene expression was quantified using the Pfaffl method [38] and primer efficiencies were calculated with the LinRegPCR software. The relative expression of all genes of interest was normalised to the expression of *RNA polymerase II polypeptide J (PolJ)* [39].

## Immunoblots

Freshly prepared RIPA buffer (50 mM Tris pH 8.8, 150 mM NaCl, 1% NP-40, 0.5% sodium dodecycholate, 0.1% sodium dodecyl-sulphate (SDS) and complete mini-proteinase inhibitors (Roche)) was used to homogenize tissue. Equal amounts of total protein were loaded in the wells, as measured by Bradford Assay. Protein samples were first diluted 1:1 with Laemmli sample buffer (Bio-Rad, Hemel Hempstead, UK) containing 5% β-mercaptoethanol (Sigma) and heated at 100 °C for 10 min. Next, samples were loaded on freshly made 1.5 mm 12% polyacrylamide separating and 5% stacking gel and electrophoresis was performed at 120 V for ~ 1.5 h in running buffer. Proteins were then transferred from the gel onto a polyvinylidene fluoride membrane (Merck Millipore) via electroblotting at 120 V for 60 min in transfer buffer containing 20% methanol. Membranes were then incubated for 2 h in Odyssey Blocking Buffer (Licor). The membrane was then probed overnight at 4 °C with the primary antibodies (p105/p50, Abcam ab32360, 1:1000; Actin, Abcam ab3280, 1:1000) in Odyssey Blocking Buffer and 0.1% Tween-20. The next day, after three 10-min washes in phosphate buffered saline (PBS), the membrane was incubated for 1 h at room temperature with secondary antibodies (goat anti-rabbit IgG 680RD, LI-COR 926–68,071, 1:1000; goat anti-mouse IgG 800CW, LI-COR, 926–32,210, 1:1000). Lastly, the membrane was washed three times for 10 min in PBS and visualized by scanning the 700 nm and 800 nm channels on the LI-COR Odyssey CLx infrared imaging system (LI-COR) for 2.5 min per channel. The background was subtracted and signal of protein of interest was divided by signal of the housekeeping protein (actin).
## Laminin staining of skeletal muscles

*Tibialis anterior* (TA) muscles were fixed in 4% para-
formaldehyde (PFA) overnight. Tissues were sectioned
(13 µm) and incubated in blocking buffer (0.3% Triton-
X, 20% foetal bovine serum (FBS) and 20% normal
goat serum in PBS) for 2 h. After blocking, tissues were
stained overnight at 4 °C with rat anti-laminin (Sigma
L0663, 1:1000) in blocking buffer. The next day, tissues
were washed in PBS and probed using a goat-anti-rat IgG
488 secondary antibody (Invitrogen A-11006, 1:500) for
one hour. PBS-washed tissues were mounted in Fluoro-
mount-G (Southern Biotech). Images were taken with a
microscope equipped with a 20X objective. Quantita-
tive assays were performed blinded on 3–5 mice for each
group and five sections per mouse. Myofiber area was
measured using Fiji (ImageJ) [40].

## Endplate staining of skeletal muscles

Endplates were stained as previously described [41]. Briefly, whole TA muscle was harvested and fixed in 4% PFA for 15 min. Muscles were incubated with α-bungarotoxin (α-BTX) conjugated to tetramethylrhodamine (BT00012, Biotium, 1:100) at RT for 30 min with ensuing PBS washes. Finally, 2–3 thin filets per muscle were sliced and mounted in Fluoromount-G (Southern Biotech). Images were taken with a confocal microscope, with a 20X objective. The experimenter quantifying end-plate size was blinded to the genotype of the animals until all measurements were finalized.

## Statistical analyses

All statistical analyses were done with the most up to
date GraphPad Prism software at time of writing. When
appropriate, a Student's unpaired two-tail *t*-test, a one-
way analysis of variance (ANOVA) or a two-way ANOVA
was used. Post-hoc analyses used are specified in Figure
Legends. Outliers were identified via the Grubbs' test.
For the Kaplan-Meier survival analysis, the log-rank test
was used and survival curves were considered signifi-
cantly different at *p* < 0.05.

## Results

### The Fn14 signalling cascade is dysregulated in skeletal muscle of *SOD1*<sup>G93A</sup> mice during disease progression

We firstly compared Fn14 expression in the skeletal
muscle of 20-week-old symptomatic *SOD1*<sup>G93A</sup> and
*SOD1*<sup>G93A</sup>;*Tweak<sup>-/-</sup>* mice, which showed no significant
difference in *Fn14* mRNA expression (Fig. 1A). This sug-
gests that genetically depleting the ligand (TWEAK)
was not sufficient to reduce the expression of the recep-
tor (Fn14). Since Fn14 is a key factor in modulating the
activity of the TWEAK/Fn14 pathway [34], its persistent
expression despite Tweak depletion may have limited the
benefits on muscle pathology and disease progression.

We thus set out to further characterize the Fn14 signal-
ling cascade in skeletal muscle of *SOD1*<sup>G93A</sup> males. We
started by reproducing our previously published data
[33] and demonstrated that *Fn14* mRNA levels in the *tib-
ialis anterior* (TA) of *SOD1*<sup>G93A</sup> and wild type (WT) mice
are similar in 4- (pre-symptomatic) and 12-week-old
(early symptomatic) animals, while there is a significant
increase in 20-week-old (late symptomatic) *SOD1*<sup>G93A</sup>
mice (Fig. 1B). We next assessed the expression of
nuclear factor kappa-light-chain-enhancer of activated
B cells (NF-κB) subunit p50, a direct downstream effec-
tor of TWEAK/Fn14 signalling in skeletal muscle [28,
42] that mediates pathological events in muscle when
chronically activated [43]. We found that the expres-
sion of NF-κB subunit p50 was significantly upregulated
in the TAs of *SOD1*<sup>G93A</sup> mice at both early symptomatic
(Fig. 1C) and late symptomatic (Fig. 1D) time-points
compared to WT animals, supporting an increased activ-
ity of TWEAK/Fn14 activity in skeletal muscle of ALS
mice. Next, we evaluated the gene expression of *PGC-1α*,
*Klf15*, *HKII* and *Glut4*. Interestingly, we observed a sig-
nificant decrease in the expression of *PGC-1α* (Fig. 1E),
*Klf15* (Fig. 1F), *HKII* (Fig. 1G) and *Glut4* (Fig. 1H) in
the TA muscles of 12- and 20-week-old *SOD1*<sup>G93A</sup> mice
compared to WT animals, providing further support for
increased Fn14 expression in *SOD1*<sup>G93A</sup> mice.

Together, our results demonstrate an aberrant hyperac-
tivity of TWEAK/Fn14 signalling in the skeletal muscle

(See figure on next page.)

**Fig. 1** Aberrant expression of the TWEAK/Fn14 signaling pathway in skeletal muscle of *SOD1*<sup>G93A</sup> ALS mice. A qPCR analysis of *Fn14* mRNA expression in gastrocnemius muscle from 20-week-old *SOD1*<sup>G93A</sup>;*Tweak<sup>+/+</sup>* and *SOD1*<sup>G93A</sup>;*Tweak<sup>-/-</sup>* males. Normalized relative expressions are compared to *SOD1*<sup>G93A</sup>;*Tweak<sup>+/+</sup>*. Data are scatter dot plot mean ± SEM, *n* = 3 animals per genotype, unpaired *t* test; ns, not significant. B qPCR analysis of *Fn14* mRNA expression in the *tibialis anterior* (TA) of *SOD1*<sup>G93A</sup> and wild type (WT) mice at 4 (pre-symptomatic), 12 (early symptomatic) and 20 (late symptomatic) weeks. Normalized relative expressions are compared to WT 4 weeks. Data are mean ± SEM, *n* = 3–4 animals per experimental group, two-way ANOVA, *****p* < 0.0001. C–D Quantification of NF-κB p50/actin protein levels in the TAs of 12- (C) and 20-week-old (D) *SOD1*<sup>G93A</sup> and WT. Images are representative immunoblots. Data are scatter dot plot mean ± SEM, *n* = 3–4 animals per experimental group, unpaired *t* test, *p* = 0.0302 (12 weeks), *p* = 0.0088 (20 weeks). E–F qPCR analysis of *PGC-1α* (E), *Klf15* (F), *HKII* (G) and *Glut4* (H) mRNAs in TAs of 4-, 12- and 20-week-old *SOD1*<sup>G93A</sup> and WT. Normalized relative expressions are compared to WT 4 weeks. Data are mean ± SEM, *n* = 3–4 animals per experimental group, two-way ANOVA, \* *p* < 0.05, *****p* < 0.001
<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 1** (See legend on previous page.)</figcaption></figure>
of *SOD1*<sup>G93A</sup> mice, impacting key regulatory downstream effectors known to influence overall skeletal muscle health and metabolic homeostasis.

## Genetic deletion of Fn14 increases survival of SOD1<sup>G93A</sup> mice

We sought to determine if decreasing TWEAK/Fn14 activity in *SOD1*<sup>G93A</sup> mice would improve muscle health and slow disease progression. As described above, we have previously modulated TWEAK activity both genetically and pharmacologically [33]. We thus decided to investigate the impact of depleting the activity of the receptor to abolish downstream signalling effector of the TWEAK/Fn14 pathway [34]. We crossed *SOD1*<sup>G93A</sup> mice with *Fn14*<sup>-/-</sup> mice [36], to generate ALS mice with a homozygous deletion of Fn14. Interestingly, we found that *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice had a significantly increased lifespan compared to *SOD1*<sup>G93A</sup> mice (females and males combined) (Fig. 2A) without any substantial improvements in weight (Fig. 2B, C). In

fact, *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> females tended to weigh less than *SOD1*<sup>G93A</sup> females (Fig. 2B), while there were no significant differences between *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> and *SOD1*<sup>G93A</sup> males (Fig. 2C). Nevertheless, Fn14 depletion appears to have an overall positive impact on disease progression in *SOD1*<sup>G93A</sup> mice.

## Genetic deletion of Fn14 improves muscle pathology in SOD1<sup>G93A</sup> mice

We next determined the impact of Fn14 depletion on previously characterised skeletal muscle pathologies in 20-week-old *SOD1*<sup>G93A</sup> males. We firstly measured the myofiber area in the gastrocnemius muscle of WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice as muscle wasting is evident in these ALS mice at that symptomatic time-point [33]. We found that while Fn14 depletion in WT animals had no impact on myofiber size, there was a significant increase in myofiber size in

<figure><img src="image_11.png" /><figcaption>**Fig. 2** Genetic deletion of *Fn14* increases survival of *SOD1*<sup>G93A</sup> ALS mice. **A** Survival curves of untreated *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice (males and females combined). Data are represented as Kaplan-Meier survival curves, *n* = 20 animals per experimental group, Log-rank (Mantel-Cox), *p* = 0.0009. **B–C** Weekly weights of *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> females (B) and males (C) from 6 weeks to humane endpoint. Data are mean ± SEM, *n* = 9–11 animals per experimental group, two-way ANOVA, \* *p* < 0.05, \*\* *p* < 0.01</figcaption></figure>

(See figure on next page.)

**Fig. 3** Genetic deletion of *Fn14* improves muscle phenotypes in *SOD1*<sup>G93A</sup> ALS mice. **A** Quantification of myofiber area of laminin-stained cross-sections of gastrocnemius muscles from 20-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> males. Data are dot plot and mean, *n* = 3–4 animals per experimental group (> 800 myofibers per experimental group), one-way ANOVA, *ns* = not significant, \*\*\*\* *p* < 0.0001. **B** Relative frequency distribution of myofiber size in gastrocnemius muscles from 20-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. **C** Representative images of laminin-stained cross-sections of gastrocnemius muscles from 20-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. **D** Quantification of neuromuscular junction endplate (EP) area of alpha-bungarotoxin-stained TA muscles from 20-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Data are dot plot and mean, *n* = 3–4 animals per experimental group (> 80 myofibers per experimental group), one-way ANOVA, *ns* = not significant, \*\* *p* < 0.01. **E** Representative images of alpha-bungarotoxin-stained TA muscles from 20-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. **F** qPCR analysis of *Fn14*, *Tweak*, *Klf15*, *Glut4*, *HKII*, *PGC-1α*, *Atrogin-1* and *MuRF-1* mRNA in TA muscles from *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Normalized relative expressions are compared to *SOD1*<sup>G93A</sup> for each gene. Data are scatter dot plot mean ± SEM, *n* = 3–8 animals per experimental group, two-way ANOVA, \* *p* < 0.05. **G** qPCR analysis of *osteopontin*, *ChAT* and *NeuN* in spinal cord from *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Normalized relative expressions are compared to *SOD1*<sup>G93A</sup> for each gene. Data are scatter dot plot mean ± SEM, *n* = 4–5 animals per experimental group, two-way ANOVA, \* *p* < 0.05
<figure><img src="image_12.png" /><figcaption>**Fig. 3** (See legend on previous page.)</figcaption></figure>
SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice compared to SOD1<sup>G93A</sup> animals (Fig. 3A-C).

We also investigated the impact of Fn14 deletion on post-synaptic neuromuscular junction (NMJ) pathologies by evaluating endplate size, which is typically reduced in ALS mice [44] and associated with muscle size [45]. Similar to myofiber size, we observed that Fn14 depletion did not influence the NMJ endplate size in the TA muscles of WT animals, but it significantly increased endplate size in SOD1<sup>G93A</sup> mice (Fig. 3D, E).

To determine the impact of Fn14 deletion on skeletal muscle at a molecular level, we investigated the gene expression of molecular effectors associated with the TWEAK/Fn14 signalling cascade (Fn14, Tweak, Klf15, Glut4, HKII and PGC-1α) [29] and muscle atrophy markers (Atrogin-1 and MuRF-1) [46]. We found that the complete elimination of Fn14 in TA muscles of SOD1<sup>G93A</sup> mice did not influence the expression of Tweak, Glut4, HKII, PGC-1α and MuRF-1 (Fig. 3F). However, we observed a significant decrease in the expression of Klf15 and, importantly, the atrogene Atrogin-1 (Fig. 3F).

Finally, given the effects observed in skeletal muscle and at the NMJ, we evaluated the effect of systemic Fn14 deletion on the expression of motor neuron markers (osteopontin [47, 48], choline acetyltransferase (ChAT) [49], and neuronal nuclear antigen (NeuN) [50]) in the spinal cord. We observed that the levels of osteopontin, ChAT and NeuN were not significantly different between SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice (Fig. 3G), implying that the origin of the benefits of Fn14 depletion on the NMJ in SOD1<sup>G93A</sup> animals is most likely the muscle and not the spinal cord.

Combined, our analyses of symptomatic mice reveal that deletion of *Fn14* in *SOD1*<sup>G93A</sup> mice improves several muscle wasting phenotypes without impacting the expression of spinal cord motor neuron markers. This suggests that the aberrant increased expression of Fn14 in skeletal muscle of *SOD1*<sup>G93A</sup> animals may contribute to the muscle pathologies that characterise the disease.

## Enhanced physical activity and Fn14 depletion both have positive effects on survival of SOD1<sup>G93A</sup> mice

We assessed if the observed molecular and histological benefits in the muscles of Fn14-depleted ALS mice translated into improved motor performance. SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice therefore performed a weekly rotarod [51] and grid test [33, 52], starting at 8 weeks and ending when the animals reached their defined humane endpoint. Both tests have previously been used in SOD1<sup>G93A</sup> mice [33, 53] and are aimed at evaluating motor balance and coordination (rotarod) and strength (grid test). We found that there were no significant differences in the time spent on the rotarod between SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> female and male mice (Fig. 4A, B). With the grid test, no significant difference in muscle strength was observed between SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> females (Fig. 4C), while SOD1<sup>G93A</sup> males were significantly stronger than SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> males at the very early pre-symptomatic time-points (Fig. 4D). Although these results suggest that Fn14 depletion does not enhance muscle strength and/or performance in SOD1<sup>G93A</sup> mice, this might be due to the independent benefits provided by the weekly rotarod and grid tests exercises. Indeed, exercised SOD1<sup>G93A</sup> animals had a significantly greater lifespan than unexercised SOD1<sup>G93A</sup> mice (Fig. 4E). As such, exercised SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice had similar survivals, suggesting that both exercise and Fn14 depletion can improve survival in SOD1<sup>G93A</sup> mice. While the median lifespan of exercised SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice were not significantly different, there did appear to be a delay in the early deaths in the exercised SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> group, pointing towards a potential combination of independent and dependent mechanisms. Of note, there was also no significant difference between the survival of heterozygous SOD1<sup>G93A</sup>; Fn14<sup>+/-</sup> and homozygous SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice that underwent weekly rotarod and grid test assessments (Supplementary Fig. 1).

(See figure on next page.)

**Fig. 4** Weekly exercise tests do not reveal any improvements in motor function of Fn14-depleted SOD1<sup>G93A</sup> mice and induce benefits on lifespan independent of Fn14 depletion. SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice performed both the rotarod and grid test weekly from 8 weeks to humane endpoint. A-B) Time in seconds (s) spent on rotarod before falling (maximum 300 s) for SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> female (A) and male (B) mice. Data are mean ± SEM, n = 5-7 animals per experimental group, two-way ANOVA, ns = not significant. C-D) Muscle strength (arbitrary units (a.u.) for SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> female (C) and male (D) mice. Data are mean ± SEM, n = 5-7 animals per experimental group, two-way ANOVA, ns = not significant, \*p < 0.05. E) Survival curves of SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup> mice that performed both the rotarod and grid test weekly from 8 weeks to humane endpoint (males and females combined). Data are represented as Kaplan-Meier survival curves, n = 11-12 animals per experimental group, Log-rank (Mantel-Cox), ns = not significant, \*\*\*p < 0.001, \*\*\*\*p < 0.0001
A Females

<figure><img src="image_13.png" /></figure>

B

<figure><img src="image_14.png" /></figure>

C Females

<figure><img src="image_15.png" /></figure>

D

<figure><img src="image_16.png" /></figure>

E

<figure><img src="image_17.png" /><figcaption>**Fig. 4** (See legend on previous page.)</figcaption></figure>
## Fn14 depletion changes molecular response of SOD1<sup>G93A</sup> muscle to exercise

To further elucidate the potential complex interactions between exercise, disease state and Fn14 depletion, 12-week-old mice underwent the rotarod or grid test for 5 consecutive days. The 12-week time point was chosen as it is an early symptomatic age for SOD1<sup>G93A</sup> mice that still allows them to complete both exercise regimens to the same extent as WT and Fn14<sup>-/-</sup> animals. The TAs were harvested 2 h after the last test and compared to those of unexercised sex- and age-matched mice for the expression of the known TWEAK/Fn14 molecular effectors and atrogenes investigated above.

First, we assessed and compared TAs from unexercised and rotarod-trained males. Interestingly, we observed that Fn14 expression was significantly upregulated in rotarod-trained SOD1<sup>G93A</sup> mice compared to unexercised SOD1<sup>G93A</sup> animals, while Fn14 levels remained unchanged in WT animals (Fig. 5A), suggesting a yet to be determined role for Fn14 in exercised SOD1<sup>G93A</sup> muscle. Next, we compared the expression of *Tweak*, *MuRF-1*, *Atrogin-1*, *Glut4*, *Klf15*, *HKII* and *PGC-1α* in unexercised and rotarod-trained WT, Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> mice. We found that *Tweak* was significantly increased only in rotarod-trained SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> animals compared to unexercised mice (Fig. 5B), suggesting a compensatory mechanism that is plausibly due to reduced levels of its ligand and exercise. The atrogene *MuRF-1* was significantly increased only in the muscles of rotarod-trained SOD1<sup>G93A</sup> mice compared to unexercised animals (Fig. 5C), indicating that depletion of Fn14 prevents exercise-induced *MuRF-1* upregulation. However, this effect appears to be specific to *MuRF-1* as *Atrogin-1*, which is significantly upregulated in rotarod-trained SOD1<sup>G93A</sup> mice, was also increased in rotarod-trained Fn14<sup>-/-</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> mice compared to unexercised cohorts (Fig. 5D). Similarly, the expression of *Glut4* was significantly increased only in rotarod-trained SOD1<sup>G93A</sup> mice compared to unexercised animals and remained unchanged in Fn14-depleted

groups (Fig. 5E). Interestingly, the expression of *Klf15* was significantly upregulated only in rotarod-trained Fn14<sup>-/-</sup> animals compared to unexercised mice (Fig. 5F). As for the expression of *HKII*, it was significantly increased in rotarod-trained SOD1<sup>G93A</sup> and Fn14<sup>-/-</sup> mice, while it remained unchanged in rotarod-trained WT mice and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> compared to unexercised groups (Fig. 5G). Finally, the expression of *PGC-1α* was significantly upregulated only in rotarod-trained SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> animals compared to unexercised mice (Fig. 5H).

Next, we performed the same investigations in TAs from unexercised and grid test-trained males. Contrary to what was observed in rotarod-trained SOD1<sup>G93A</sup> males (Fig. 5A), we found that grid test-trained SOD1<sup>G93A</sup> mice expressed significantly less Fn14 than unexercised SOD1<sup>G93A</sup> animals (Fig. 5I), suggesting a distinct response between rotarod and grid test activities. On the other hand, *Tweak* expression was significantly increased only in grid test-trained WT animals compared to unexercised mice, while it remained unchanged in grid test-trained animals of the same genotype (Fig. 5J). The expression of both atrogenes, *MuRF-1* and *Atrogin-1*, was significantly upregulated only in grid-test trained SOD1<sup>G93A</sup> mice compared to unexercised animals and restored to low levels when Fn14 was depleted (Fig. 5K, L). *Glut4* levels were unchanged in all experimental groups when comparing unexercised animals to grid test-trained mice (Fig. 5M). Interestingly, *Klf15* expression was significantly upregulated only in SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> animals compared to unexercised mice, as it remained unchanged in all other groups (Fig. 5N). Similar to *Glut4*, *HKII* levels were also unchanged in all experimental groups when comparing unexercised animals to grid test-trained mice (Fig. 5O). Finally, *PGC-1α* expression was significantly increased only in grid test-trained SOD1<sup>G93A</sup> mice compared to unexercised animals and returned to lower levels in Fn14-depleted animals (Fig. 5P).

The same analyses were then performed in females to see if sex was an additional factor influencing the

(See figure on next page.)

**Fig. 5** Type of exercise and genotype impact the expression of *Tweak*, *Fn14* and their downstream effectors in male mice. 12-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> males were either placed on the rotarod (A-H) or performed the grid test (I-P) daily for 5 consecutive days. *Tibialis anterior* (TA) muscles were harvested approximately 2 h after the last bout of exercise. A qPCR analysis of *Fn14* mRNA expression in unexercised and rotarod-exercised WT and *SOD1*<sup>G93A</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 3-4 animals per experimental group, two-way ANOVA, ns = not significant, \*\*p < 0.01. B-H qPCR analysis of *Tweak* (*B*), *MuRF-1* (*C*), *Atrogin-1* (*D*), *Glut4* (*E*), *Klf15* (*F*), *HKII* (*G*), *PGC-1α* (*H*) mRNA expression in unexercised and rotarod-exercised WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 4-9 animals per experimental group, two-way ANOVA, ns = not significant, \*p < 0.05, \*\*\*p < 0.001, \*\*\*\*p < 0.0001. I qPCR analysis of *Fn14* mRNA expression in unexercised and grid test-exercised WT and *SOD1*<sup>G93A</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 3-4 animals per experimental group, two-way ANOVA, ns = not significant, \*\*p < 0.01. J-P qPCR analysis of *Tweak* (*J*), *MuRF-1* (*K*), *Atrogin-1* (*L*), *Glut4* (*M*), *Klf15* (*N*), *HKII* (*O*), *PGC-1α* (*P*) mRNA expression in unexercised and grid test-exercised WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 4-9 animals per experimental group, two-way ANOVA, ns = not significant, \*p < 0.05, \*\*\*p < 0.001, \*\*\*\*p < 0.0001
<figure><img src="image_18.png" /><figcaption>**Fig. 5** (See legend on previous page.)</figcaption></figure>
interactions between exercise, disease state and Fn14 depletion. Indeed, in the rotarod female experimental groups, we found that *Fn14* expression was significantly elevated in rotarod-trained WT mice compared to unexercised WT animals, while remaining unchanged in *SOD1*<sup>G93A</sup> mice (Fig. 6A). *Tweak* levels were significantly upregulated only in rotarod-trained *Fn14*<sup>-/-</sup> animals compared to unexercised mice (Fig. 6B). As for the atrogenes, we observed a significant increased expression of *MuRF-1* only in rotarod-trained WT animals compared to unexercised mice (Fig. 6C) and no significant changes in any experimental groups for *Atrogin-1* (Fig. 6D). *Glut4* levels were significantly elevated only in rotarod-trained *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> animals compared to unexercised mice (Fig. 6E) while the expression of *Klf15*, *HKII* and *PGC-1α* expression were significantly increased only in rotarod-trained WT mice compared to unexercised animals (Fig. 6F-H).

When performing the same comparisons in grid test female experimental groups, we observed a significant increased expression of *Fn14* only in grid test-trained WT animals compared to unexercised mice (Fig. 6I). *Tweak* levels were significantly elevated in both grid test-trained WT and *Fn14*<sup>-/-</sup> mice compared to unexercised animals (Fig. 6J). Interestingly, *MuRF-1* expression was significantly increased in all grid test-trained groups compared to untrained animals (Fig. 6K) while *Atrogin-1* levels were significantly elevated only in grid test-trained *Fn14*<sup>-/-</sup> mice compared to unexercised animals (Fig. 6L). *Glut4* expression was significantly increased in both grid test-trained WT and *Fn14*<sup>-/-</sup> mice compared to unexercised animals (Fig. 6M). *Klf15* levels were significantly elevated only in grid test-trained *Fn14*<sup>-/-</sup> animals compared to unexercised mice (Fig. 6N) while *HKII* expression was significantly increased in all grid test-trained experimental groups compared to unexercised animals except for *SOD1*<sup>G93A</sup> mice where the levels remained unchanged (Fig. 6O). Finally, similar to *Glut4*, *PGC-1α* levels were significantly elevated in both grid test-trained WT

and *Fn14*<sup>-/-</sup> mice compared to unexercised animals (Fig. 6P).

Our overall results suggest that exercise regimens have a differential impact on the skeletal muscle of our 12-week-old experimental cohorts, pointing towards specific interactions between sex, genotype, Fn14 depletion and exercise (summarised in Table 1).

## Fn14 depletion influences the impact of exercise on expression of myosin heavy chain isoforms

The physiological and biochemical properties of skeletal muscle are, in part, conferred by myosin heavy chain (MyHC) isoforms. Of those, MyHC isoforms IIa, IIx and IIb are specific to fast twitch type of muscles, have been demonstrated to co-exist in the same muscle and their expression can be influenced by exercise [54, 55]. Given that the TA is considered a fast twitch skeletal muscle, we investigated the impact of rotarod and grid test exercise regimens on the gene expression of MyHC IIa, IIx and IIb isoforms in 12-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> males and females.

In rotarod-trained males, we found that the expression of all three *MyHC* isoforms remained unchanged in WT and *SOD1*<sup>G93A</sup> animals (Fig. 7A, B). Interestingly, we observed a significant increased expression of *MyHC* IIb in rotarod-trained *Fn14*<sup>-/-</sup> mice compared to unexercised animals (Fig. 7C) while the levels of all three *MyHC* isoforms were similar between rotarod-trained and unexercised *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice (Fig. 7D).

In grid test-trained males, only *MyHC* IIa expression was significantly increased in WT animals (Fig. 7E). Similar to the rotarod, the levels of all three *MyHC* isoforms remained unchanged in grid test-trained *SOD1*<sup>G93A</sup> mice (Fig. 7F). Similar to WT animals, only *MyHC* IIa expression was significantly upregulated in grid test-trained *Fn14*<sup>-/-</sup> mice compared to unexercised animals (Fig. 7G). In contrast, the expression of all three *MyHC* isoforms was significantly downregulated in grid test-trained *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> animals compared to unexercised mice (Fig. 7H).

(See figure on next page.)

**Fig. 6** Type of exercise and genotype impact the expression of *Tweak*, *Fn14* and their downstream effectors in female mice. 12-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> females were either placed on the rotarod (A-H) or performed the grid test (I-P) daily for 5 consecutive days. *Tibialis anterior* (TA) muscles were harvested approximately 2 h after the last bout of exercise. **A** qPCR analysis of *Fn14* mRNA expression in unexercised and rotarod-exercised WT and *SOD1*<sup>G93A</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 4–6 animals per experimental group, two-way ANOVA, ns = not significant, ***p < 0.001. **B-H** qPCR analysis of *Tweak* (**B**), *MuRF-1* (**C**), *Atrogin-1* (**D**), *Glut4* (**E**), *Klf15* (**F**), *HKII* (**G**), *PGC-1α* (**H**) mRNA expression in unexercised and rotarod-exercised WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 3–8 animals per experimental group, two-way ANOVA, ns = not significant, *p* < 0.05, **p* < 0.01, ***p* < 0.001. **I** qPCR analysis of *Fn14* mRNA expression in unexercised and grid test-exercised WT and *SOD1*<sup>G93A</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 3–4 animals per experimental group, two-way ANOVA, ns = not significant, **p* < 0.01. **J-P** qPCR analysis of *Tweak* (**J**), *MuRF-1* (**K**), *Atrogin-1* (**L**), *Glut4* (**M**), *Klf15* (**N**), *HKII* (**O**), *PGC-1α* (**P**) mRNA expression in unexercised and grid test-exercised WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> mice. Data are scatter dot plot mean ± SEM, *n* = 3–5 animals per experimental group, two-way ANOVA, ns = not significant, *p* < 0.05, **p* < 0.01, ***p* < 0.001, ****p* < 0.0001
<figure><img src="image_19.png" /><figcaption>**Fig. 6** (See legend on previous page.)</figcaption></figure>
<table><caption>**Table 1** Effect of exercise, sex and genotype on expression of the TWEAK/Fn14 signalling pathway and atrogenes</caption><thead><tr><th rowspan="3"></th><th colspan="4">WT</th><th colspan="4">SOD1<sup>G93A</sup></th><th colspan="4">Fn14<sup>-/-</sup></th><th colspan="4">SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup></th></tr><tr><th colspan="2">Rotarod</th><th colspan="2">Grid test</th><th colspan="2">Rotarod</th><th colspan="2">Grid test</th><th colspan="2">Rotarod</th><th colspan="2">Grid test</th><th colspan="2">Rotarod</th><th colspan="2">Grid test</th></tr><tr><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th></tr></thead><tbody><tr><td>Fn14</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>Tweak</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>MuRF-1</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>Atrogin-1</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>Glut4</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>Klf15</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>HKII</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr><tr><td>PGC-1α</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td><td>□</td></tr></tbody><tfoot><tr><td colspan="18">Red = upregulated<br>Blue = downregulated<br>Grey = no change<br>Black = not applicable<br>M = Males<br>F = Females</td></tr></tfoot></table>

In rotarod-trained females, we observed a significant increased expression of only *MyHC IIb* in both rotarod-trained WT and *SOD1*<sup>G93A</sup> mice compared to unexercised animals (Fig. 8A, B). In contrast, the expression of all three *MyHC* isoforms remained unchanged in rotarod-trained *Fn14*<sup>-/-</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> animals compared to unexercised mice (Fig. 8C, D).

Finally, in grid test-trained females, only the expression of *MyHC IIa* was significantly increased in grid test-trained WT animals compared to unexercised mice (Fig. 8E) while the levels of both *MyHC IIa* and *IIb* were significantly elevated in grid test-trained *SOD1*<sup>G93A</sup> mice compared to unexercised animals (Fig. 8F). The

expression of all three *MyHC* isoforms remained unchanged in grid test-trained *Fn14*<sup>-/-</sup> animals compared to unexercised mice (Fig. 8G). In *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> animals, we observed a significant increase only in *MyHC IIx* expression compared to unexercised mice (Fig. 8H).

Combined, our results suggest that the rotarod and grid test exercise regimens did not lead to changes in *MyHC* isoform expression in all cases. Indeed, similar to our analyses of molecular effectors of the TWEAK/Fn14 signalling pathway and atrogenes, our data points to the influence of sex, genotype, exercise type and disease state on the expression of *MyHC IIa*, *IIx* and *IIb* isoforms in the TA muscle (summarised in Table 2).

(See figure on next page.)

**Fig. 7** Type of exercise and genotype impact the expression of *MyHC* isoforms *IIa*, *IIx* and *IIb* in male mice. 12-week-old WT, *Fn14*<sup>-/-</sup>, *SOD1*<sup>G93A</sup> and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> males were either placed on the rotarod (A-D) or performed the grid test (E-H) daily for 5 consecutive days. *Tibialis anterior* (TA) muscles were harvested approximately 2 h after the last bout of exercise. **A-D** qPCR analysis of *MyHC IIa*, *IIx* and *IIb* mRNA expression in unexercised and rotarod-exercised WT (**A**), *SOD1*<sup>G93A</sup> (**B**), *Fn14*<sup>-/-</sup> (**C**) and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> (**D**) males. Data are scatter dot plot mean ± SEM, *n* = 3–4 animals per experimental group, two-way ANOVA, ns = not significant, ***p < 0.001. **E-H** qPCR analysis of *MyHC IIa*, *IIx* and *IIb* mRNA expression in unexercised and grid test-exercised WT (**E**), *SOD1*<sup>G93A</sup> (**F**), *Fn14*<sup>-/-</sup> (**G**) and *SOD1*<sup>G93A</sup>;*Fn14*<sup>-/-</sup> (**H**) males. Data are scatter dot plot mean ± SEM, *n* = 3–5 animals per experimental group, two-way ANOVA, ns = not significant, *p* < 0.05
Rotarod (Males)

A

WT

<figure><img src="image_20.png" /></figure>

E

<figure><img src="image_21.png" /></figure>

SOD1<sup>G93A</sup>

B

<figure><img src="image_22.png" /></figure>

F

<figure><img src="image_23.png" /></figure>

Fn14<sup>-/-</sup>

C

<figure><img src="image_24.png" /></figure>

G

<figure><img src="image_25.png" /></figure>

SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup>

D

<figure><img src="image_26.png" /></figure>

H

<figure><img src="image_27.png" /><figcaption>**Fig. 7** (See legend on previous page.)</figcaption></figure>
## Fn14 depletion and exercise influence myofiber size

We next determined if Fn14 depletion impacted myofiber size in the gastrocnemius muscle of rotarod- or grid test-trained 12-week-old WT, Fn14⁻/⁻, SOD1G93A and SOD1G93A; Fn14⁻/⁻ males and females.

Following rotarod training in males, we observed a significant increase in myofiber size of rotarod-trained WT animals compared to unexercised WT mice while this type of exercise did not impact myofiber size in SOD1G93A animals (Fig. 9A). Interestingly, the myofiber size of rotarod-trained Fn14⁻/⁻ and SOD1G93A; Fn14⁻/⁻ mice was significantly smaller than in unexercised control animals (Fig. 9A), indicating that the combination of Fn14 depletion and rotarod exercise can reduce muscle size.

Similar to rotarod-trained WT males, grid test-trained WT males displayed a significant increase in myofiber size compared to unexercised animals (Fig. 9B). However, unlike rotarod-trained SOD1G93A males, grid test-trained SOD1G93A animals had a significant decrease in myofiber size compared to unexercised SOD1G93A mice (Fig. 9B). A significant decrease in myofiber area was also observed in grid test-trained Fn14⁻/⁻ and SOD1G93A; Fn14⁻/⁻ mice compared to unexercised control animals (Fig. 9B), suggesting that both the SOD1G93A genotype and Fn14 depletion negatively impact muscle size following a resistance exercise regimen.

Interestingly, rotarod training in females resulted in a significant decrease in myofiber size in all rotarod-trained experimental groups when compared to unexercised animals (Fig. 9C).

In grid test-females, a significant decrease in myofiber size was also observed in grid test-trained WT, SOD1G93A and SOD1G93A; Fn14⁻/⁻ mice compared to unexercised animals (Fig. 9D). However, a significant increase in myofiber size was found in grid test-trained Fn14⁻/⁻ animals compared to unexercised mice (Fig. 9D).

Together, our data points to independent influences of exercise type, sex and genetics on muscle fiber size.

## Discussion

In this study, we aimed to better understand how increased Fn14 expression in an ALS mouse model with chronic denervation and muscle wasting could contribute to muscle pathology and disease progression. To

achieve this, we genetically deleted Fn14 in both WT and SOD1G93A mice and observed behavioural, molecular and histological changes that were dependent on exercise, sex and disease progression.

In the first instance, we not only confirmed our previous observation of increased expression of Fn14 in the skeletal muscle of SOD1G93A ALS mice during disease progression [33], but we also validated the previously reported negative correlation between the activity of the TWEAK/Fn14 pathway and the expression of the known molecular and metabolic effectors Glut4, Klf15, HKII and PGC-1α [29]. Interestingly, we recently demonstrated a similar but inverse negative correlation in the skeletal muscle of another neuromuscular mouse model, SMA mice, whereby the expression levels of Tweak and Fn14 decreased during disease progression while those of Glut4, Klf15, HKII and PGC-1α increased [32]. One important distinction between these two studies is the developmental stage investigated. Indeed, SMA mice were of pre-weaned age [32] while the SOD1G93A ALS mice were at adult stages (current study), suggesting that the TWEAK/Fn14 signalling pathway is differentially regulated at different stages of muscle development. This differential regulation might have an impact on downstream metabolic requirements and regulation as well as therapeutic interventions in cases of dysregulation.

One of our key findings is the extended lifespan of Fn14-depleted SOD1G93A ALS mice, which is contrary to the absence of impact following genetic Tweak deletion in the same mouse model, which we have previously reported [33]. This suggest that the detrimental effect of the aberrant activity of the TWEAK/Fn14 pathway in skeletal muscle of SOD1G93A ALS mice is driven by the receptor (Fn14) and not the ligand (TWEAK). This aligns with previous work that points to a greater role for Fn14 than TWEAK in enabling pathway activity [34]. It is also possible that the differential impacts observed in TWEAK- and Fn14-depleted SOD1G93A ALS mice are due to Fn14-independent TWEAK signalling [56] and/or TWEAK-independent Fn14 signalling [57]. Furthermore, the distinct effects of TWEAK and Fn14 depletion in SOD1G93A ALS mice could further be caused by their known roles in other tissues such as the heart, gastrointestinal tract, kidney, liver, central nervous system and

(See figure on next page.)

**Fig. 8** Type of exercise and genotype impact the expression of *MyHC* isoforms *IIa*, *IIx* and *IIb* in female mice. 12-week-old WT, *Fn14*⁻/⁻, *SOD1*G93A and *SOD1*G93A;*Fn14*⁻/⁻ females were either placed on the rotarod (A-D) or performed the grid test (E-H) daily for 5 consecutive days. *Tibialis anterior* (TA) muscles were harvested approximately 2 h after the last bout of exercise. A-D qPCR analysis of *MyHC* *IIa*, *IIx* and *IIb* mRNA expression in unexercised and rotarod-exercised WT (**A**), *SOD1*G93A (**B**), *Fn14*⁻/⁻ (**C**) and *SOD1*G93A;*Fn14*⁻/⁻ (**D**) females. Data are scatter dot plot mean ± SEM, *n* = 4–8 animals per experimental group, two-way ANOVA, *ns* = not significant, ***p* < 0.001. ***p* < 0.0001. E-H qPCR analysis of *MyHC* *IIa*, *IIx* and *IIb* mRNA expression in unexercised and grid test-exercised WT (**E**), *SOD1*G93A (**F**), *Fn14*⁻/⁻ (**G**) and *SOD1*G93A;*Fn14*⁻/⁻ (**H**) females. Data are scatter dot plot mean ± SEM, *n* = 3–8 animals per experimental group, two-way ANOVA, *ns* = not significant, ***p* < 0.001. ***p* < 0.0001
Rotarod (Females)

A

WT

<figure><img src="image_28.png" /></figure>

E

<figure><img src="image_29.png" /></figure>

Grid test (Females)

B

SOD1<sup>G93A</sup>

<figure><img src="image_30.png" /></figure>

F

<figure><img src="image_31.png" /></figure>

C

Fn14<sup>-/-</sup>

<figure><img src="image_32.png" /></figure>

G

<figure><img src="image_33.png" /></figure>

D

SOD1<sup>G93A</sup>; Fn14<sup>-/-</sup>

<figure><img src="image_34.png" /></figure>

H

<figure><img src="image_35.png" /><figcaption>**Fig. 8** (See legend on previous page.)</figcaption></figure>
<table><caption>**Table 2** Effect of exercise, sex and genotype on expression of the myosin heavy chain isoforms</caption><thead><tr><th rowspan="3"></th><th colspan="4">WT</th><th colspan="4">SOD1<sup>G93A</sup></th><th colspan="4">Fn14<sup>-/-</sup></th><th colspan="4">SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup></th></tr><tr><th colspan="2">Rotarod</th><th colspan="2">Grid test</th><th colspan="2">Rotarod</th><th colspan="2">Grid test</th><th colspan="2">Rotarod</th><th colspan="2">Grid test</th><th colspan="2">Rotarod</th><th colspan="2">Grid test</th></tr><tr><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th><th>M</th><th>F</th></tr></thead><tbody><tr><td>MyHC IIa</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MyHC IIx</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>MyHC IIb</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Red = upregulated

Blue = downregulated

Grey = no change

M = Males

F = Females

epithelium [58–60]. As the genetic knock-out of TWEAK and Fn14 was systemic in both cases, we cannot exclude additional benefits or detrimental effects stemming from altered function in other cells and tissues. Regardless of the reasons, our combined studies point to a greater therapeutic value in modulating Fn14 over TWEAK.

In addition to lifespan, we also observed improvements in skeletal muscle pathology at molecular and histological levels in Fn14-depleted SOD1<sup>G93A</sup> mice. These changes did not occur in Fn14<sup>-/-</sup> mice when compared to WT animals, suggesting that the effects were dependent on disease stage. Interestingly, we have previously shown increased muscle fibre and NMJ endplate sizes in TWEAK-depleted SOD1<sup>G93A</sup> ALS mice [33], further supporting a role for the TWEAK/Fn14 pathway in muscle pathology in this mouse model and in more general adult denervation-induced muscle atrophy [34]. Of note is that in both TWEAK- and Fn14-depleted SOD1<sup>G93A</sup>

ALS mice, there were no significant improvements in motor function [33], suggesting that simply targeting the TWEAK/Fn14 pathway is not sufficient for the recovery of the neuromuscular unit.

Surprisingly, the beneficial impact of Fn14 depletion on the survival of SOD1<sup>G93A</sup> ALS mice was almost masked when the mice underwent weekly rotarod and grid test assessments for approximately 16 weeks as the enhanced physical activity itself had a positive impact on survival of the SOD1<sup>G93A</sup> ALS mice. Interestingly, there are reports of both beneficial and detrimental effects of exercise in ALS mouse models and patients that suggest that exercise regimen (type and length) and sex are important factors that contribute to the observed outcomes [61–65]. In our study, short weekly bouts of grid test and rotarod were sufficient to improve survival. Further investigations showed that the combination of 5 consecutive days of exercise (rotarod or grid test) and Fn14 depletion was

(See figure on next page.)

**Fig. 9** Type of exercise, genotype and sex impact muscle fibre size. 12-week-old WT, Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> males and females either performed the grid test or were placed on the rotarod daily for 5 consecutive days. Gastrocnemius muscles were harvested approximately 2 h after the last bout of exercise. **A** Quantification of myofiber area of laminin-stained cross-sections of gastrocnemius muscles from 12-week-old unexercised and rotarod-exercised WT, Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> males. Data are dot plot and mean, *n* = 3–4 animals per experimental group (> 1100 myofibers per experimental group), two-way ANOVA, ns = not significant, ****p < 0.0001. **B** Quantification of myofiber area of laminin-stained cross-sections of gastrocnemius muscles from 12-week-old unexercised and grid test-exercised WT, Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> males. Data are dot plot and mean, *n* = 3–4 animals per experimental group (> 1100 myofibers per experimental group), two-way ANOVA, *p* < 0.05, ****p < 0.0001. **C** Quantification of myofiber area of laminin-stained cross-sections of gastrocnemius muscles from 12-week-old unexercised and rotarod-exercised WT, Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> females. Data are dot plot and mean, *n* = 3–4 animals per experimental group (> 1400 myofibers per experimental group), two-way ANOVA, ns = not significant, ****p < 0.0001. **D** Quantification of myofiber area of laminin-stained cross-sections of gastrocnemius muscles from 12-week-old unexercised and grid test-exercised WT, Fn14<sup>-/-</sup>, SOD1<sup>G93A</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> females. Data are dot plot and mean, *n* = 3–4 animals per experimental group (> 1100 myofibers per experimental group), two-way ANOVA, **p* < 0.01, ****p < 0.0001
A

Rotarod_Males

<figure><img src="image_36.png" /></figure>

C

Rotarod_Females

<figure><img src="image_37.png" /></figure>

B

Grid_test_Males

<figure><img src="image_38.png" /></figure>

D

Grid_test_Females

<figure><img src="image_39.png" /><figcaption>**Fig. 9** (See legend on previous page.)</figcaption></figure>

sufficient to induce changes at molecular and histological levels in the skeletal muscle of 12-week-old animals. These changes were dependent on disease stage, exercise and sex. Combining exercise and Fn14 depletion may therefore lead to potentially, additive, synergistic and/or

antagonistic interactions that may be dependent on the exercise regimen itself and individual characteristics. However, in our study, we did not identify clear commonalities between Fn14 depletion and exercises that would point to shared mechanisms. Furthermore, the results
of some analyses such as myofiber size and expression of atrogenes did not necessarily align with the increased survival in Fn14-depleted and exercised $SOD1^{G93A}$ mice, suggesting more complex and possibly multi-systemic mechanisms that influence overall disease progression in these mice.

One key observation was that changes in *Tweak* and *Fn14* expression appeared dependent on the type of exercise, sex and genotype of the animal. For example, *Fn14* levels displayed a differential expression in $SOD1^{G93A}$ males only, whereby it was increased following rotarod and decreased following the grid test. Conversely, in females, *Fn14* levels specifically increased in WT mice following both rotarod and grid test, when compared to unexercised animals. These diverse patterns may reflect the complex metabolic adaptations impacted by disease, Fn14 presence/absence, sex and type of exercise. Typically, endurance exercises promote the use of aerobic/ oxidative metabolic pathways in skeletal muscle while resistance exercises favour anaerobic/glycolytic metabolic pathways [66]. In ALS, skeletal muscle metabolism during rest and exercise is altered in both pre-clinical models and patients [67–70], which could alter how ALS muscle adapts to different types of exercises and the overall beneficial vs detrimental outcomes [62, 63]. As for Fn14, it is typically increased in skeletal muscle of healthy individuals and adult mice following exercise, irrespective of type (endurance vs resistance) [71–74]. Conversely, the muscle-specific deletion of Fn14 and the ubiquitous TWEAK deletion in mice both improved exercise capacity and oxidative metabolism [75, 76], suggesting that sustained and/or aberrant increase in TWEAK/Fn14 activity expression during exercise may be detrimental. It is therefore unclear why the expression of both the ligand and effector are commonly reported as being elevated following exercise. Of note, we only observed changes in *Fn14* expression in exercised WT female mice in our study, which may be due to our selected exercise regimens (length and type of exercise). Nevertheless, our results, combined with previous studies, suggest and support a complex interaction between Fn14 regulation, disease state, exercise and the metabolic status of muscle.

Another noticeable result is the influence of genotype, sex and exercise on the expression of the atrogenes *MuRF-1* and *Atrogin-1*. For example, we found that the expression of *MuRF-1* is significantly elevated in $SOD1^{G93A}$ males following both the rotarod and grid test, supporting previous studies on the negative impact of exercise in ALS patients [63, 77]. In $SOD1^{G93A}$;*Fn14*⁻/⁻ males, *MuRF-1* levels remained low in both rotarod and grid test groups, aligning with the previous report of reduced neurogenic muscle atrophy in muscle-specific Fn14-depleted animals [75]. However, in females,

*MuRF-1* expression was significantly increased in WT mice only following the rotarod and all experimental groups after the grid test. While the differential expression patterns of both atrogenes might appear contradictory, previous studies have demonstrated that their regulation can be controlled by distinct pathways and in a sex-dependent fashion [78–82]. Of note, our analysis of muscle fibres shows an absence of perfect correlation between the expression of atrogenes and myofiber size, suggesting that changes in *MuRF-1* and *Atrogin-1* are not sufficient to improve muscle size and that other molecular effectors and regulatory pathways may be responsible for modulating muscle mass [83]. Indeed, Fn14 has previously been demonstrated to modulate myoblast fusion [84, 85], a process that contributes to muscle size growth during regeneration, which typically occurs following bouts of exercise. As such, Fn14 depletion may have affected myofiber size in a subset of our rotarod- and grid test-trained animals.

The expression of molecular and metabolic effectors previously shown to be regulated by TWEAK/Fn14 signalling also appear to be dependent on genotype, sex and type of exercise. For example, *PGC-1α* expression was upregulated in both the rotarod- and grid test-trained $SOD1^{G93A}$ males and Fn14 depletion restored the levels to normal only in the grid-test exercised $SOD1^{G93A}$;*Fn14*⁻/⁻ mice. In females however, *PGC-1α* levels were significantly elevated in rotarod- and grid-test trained WT animals and Fn14 depletion restored the levels to normal only in the rotarod-trained WT females. These differential patterns and relationships between genotype, exercise, sex and metabolic effectors most likely result from the combination of different metabolic pathways favoured by different types of exercise [66], as well as the impact of ALS-causing mutations and sex on muscle metabolism [67–70, 86, 87]. Indeed, the previously reported role of the TWEAK/Fn14 pathway in the regulation of *PGC-1α* and mitochondrial content in skeletal muscle [88] may have contributed to our *PGC-1α* observations with added influences from sex, genotype and exercise that still need to be explored.

MyHC isoforms IIa, IIx and IIb are specific to fast twitch skeletal muscles such as the TA and they each confer distinct metabolic and functional properties to skeletal muscle fibres [89]. Furthermore, skeletal muscle can adapt rapidly to metabolic changes induced by exercise, which can be reflected by altered MyHC expression [90]. In our study, we indeed observed some changes in MyHC isoform expression that were dependent on exercise type, genotype and sex, suggesting that the different combinations of these factors had distinct effects on the properties of skeletal muscle. In some instances (e.g.
grid test-trained SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> males), changes in the same direction of more than one MyHC isoform was observed, a phenomenon previously reported following a short bout of exercise [91].

Finally, our analyses of myofiber size further emphasized the differential responses of skeletal muscle from females and males to types of exercise, disease state and Fn14 depletion. In fact, throughout our molecular and histological analyses, sex-dependent differences were observed. These align with previous studies showing the effect of sex on responses to exercise [80, 81, 87], in Fn14<sup>-/-</sup> mice exposed to neonatal hypoxia-ischemia [92], on general muscle properties [93], on expression of MyHC isoforms [94] and on disease onset and response to exercise in SOD1<sup>G93A</sup> mice [95]. All of these extrinsic and intrinsic factors are therefore important to consider when assessing skeletal muscle adaptations in ALS. While our work provides some interesting insights, it is important to note its key limitations. Firstly, the impact of exercise on SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> mice was observed in animals that performed both types of exercise weekly from 8 weeks of age to humane endpoint. However, the rotarod and grid test experiments were done separately on 12-week-old animals for 1 week only. Furthermore, our study focused on the known metabolic effectors downstream of TWEAK and Fn14, which means that additional genes and signalling cascades could be impacted by exercise, sex and/or genotype and contribute to our observed results. Finally, our research was aimed at investigating skeletal muscle but as Fn14 depletion is systemic, some of the beneficial and detrimental effects reported may be due to other cells and tissues.

## Conclusions

Our study provides additional insights on the role of the TWEAK/Fn14 pathway in a denervation-induced muscle pathology as modelled in the SOD1<sup>G93A</sup> ALS mice. Importantly, we demonstrate that the benefits of Fn14 depletion are impacted by exercise and sex. This is particularly relevant in the context of the current therapeutic landscape of the ALS field, where combinatorial therapies that include exercise regimens are being explored by many research and clinical teams. As such, a better understanding and consideration of the interactions between treatments, muscle metabolism, exercise and sex will be of importance in future studies.

### Abbreviations

α-BTX α-Bungarotoxin

ALS Amyotrophic lateral sclerosis

C9ORF72 Chromosome 9 open reading frame 72

ChAT Choline acetyltransferase

FBS Fetal bovine serum

Fn14 Fibroblast growth factor inducible 14

FUS Fused in Sarcoma

Glut4 Glucose transporter 4

HKII Hexokinase II

Klf15 Krüppel-like factor 15

MyHC Myosin heavy chain

NeuN Neuronal nuclear antigen

NF-κB Nuclear factor kappa-light-chain-enhancer of activated B cells

PBS Phosphate buffered saline

PFA Paraformaldehyde

PGC-1α Peroxisome proliferator-activated receptor-gamma coactivator 1α

PolJ RNA polymerase II polypeptide J

SDS Sodium dodecyl-sulfate

SMA Spinal muscular atrophy

SOD1 Superoxide dismutase 1

TA Tibialis anterior

TDP-43 TAR DNA-binding protein 43

TWEAK Tumor necrosis factor-like weak inducer of apoptosis

## Supplementary Information

The online version contains supplementary material available at https://doi.org/10.1186/s13395-024-00356-0.

Additional file 1: Supplementary Figure 1. Survival of SOD1<sup>G93A</sup>;Fn14<sup>+/-</sup> and are not significantly different. Survival curves of SOD1<sup>G93A</sup>;Fn14<sup>+/-</sup> and SOD1<sup>G93A</sup>;Fn14<sup>-/-</sup> mice that performed both the rotarod and grid test weekly from 8 weeks to humane endpoint (males and females combined). Data are represented as Kaplan-Meier survival curves, n= 12-13 animals per experimental group, Log-rank (Mantel-Cox), ns = not significant.

Additional file 2: Supplementary Table 1. Mouse primers used for quantitative real-time PCR.

### Acknowledgements

We would like to thank the staff at the BMS facility at the University of Oxford, Dr. Linda Burkly and Biogen.

### Authors' contributions

Conceptualisation: G.H., C.R., M.B.; Methodology: G.H., M.B.; Validation: G.H., M.B.; Formal analysis: G.H., E.M., N.A., E.R.S., M.O., E.I.H.S., J.M.H., T.S., S.I., A.B., A.K., F.S., M.B.; Investigation: G.H., E.M., N.A., E.R.S., M.O., E.I.H.S., J.M.H., T.S., S.I., A.B., A.K., F.S., M.B.; Writing-original draft preparation: M.B.; Writing-review and editing: G.H., E.M., N.A., E.R.S., M.O., E.I.H.S., J.M.H., T.S., S.I., A.B., A.K., F.S., C.R., M.J.A.W., M.B.; Visualisation: M.B.; Supervision: C.R., M.J.A.W., M.B.; Project administration: G.H., M.B.; Funding acquisition: C.R., M.J.A.W., M.B.

### Funding

J.M.H. was funded by a Ph.D. studentship from the Keele University of School of Medicine. E.M. was supported by an Academy of Medical Sciences grant (SBF006/1162) and is currently supported by a Medical Research Council grant (MR/Y003640/1). E.R.S. was funded by a MDUK Ph.D. studentship (18GRO-PS48-0114). E.I.H.S. is funded by a studentship from the Islamic University of Gaza.

### Data availability

All data generated or analysed during this study are either included in this published article [and its supplementary information files] or are available from the corresponding author on reasonable request.

### Declarations

#### Ethics approval and consent to participate

Not applicable.

#### Consent for publication

Not applicable.

#### Competing interests

The authors declare no competing interests.

#### Author details

¹Department of Physiology, Anatomy and Genetics, University of Oxford, Oxford, UK. ²School of Medicine, Keele University, Staffordshire, UK.
³Department of Paediatrics, University of Oxford, Oxford, UK. ⁴School of Life Sciences, Keele University, Staffordshire, UK. ⁵INM, Univ Montpellier, INSERM, Montpellier, France. ⁶ALS Reference Center, Univ Montpellier, CHU Montpellier, Montpellier, France. ⁷Wolfson Centre for Inherited Neuromuscular Disease, RJAH Orthopaedic Hospital, Oswestry, UK.

Received: 5 July 2024 Accepted: 4 October 2024

Published online: 14 October 2024

## References

1. Al-Chalabi A, Jones A, Troakes C, King A, Al-Sarraj S, van den Berg LH. The genetics and neuropathology of amyotrophic lateral sclerosis. Acta Neuropathol. 2012;124:339–52.

2. Andersen PM, Al-Chalabi A. Clinical genetics of amyotrophic lateral sclerosis: what do we really know? Nat Rev Neurol. 2011;7:603–15.

3. DeJesus-Hernandez M, Mackenzie IR, Boeve BF, Boxer AL, Baker M, Rutherford NJ, et al. Expanded GGGGCC hexanucleotide repeat in non-coding region of C9ORF72 causes chromosome 9p-linked FTD and ALS. Neuron. 2011;72:245–56.

4. Renton AE, Majounie E, Waite A, Simón-Sánchez J, Rollinson S, Gibbs JR, et al. A hexanucleotide repeat expansion in C9ORF72 is the cause of chromosome 9p21-linked ALS-FTD. Neuron. 2011;72:257–68.

5. Rosen DR, Siddique T, Patterson D, Figlewicz DA, Sapp P, Hentati A, et al. Mutations in Cu/Zn superoxide dismutase gene are associated with familial amyotrophic lateral sclerosis. Nature. 1993;362:59–62.

6. Kwiatkowski TJ, Bosco DA, Leclerc AL, Tamrazian E, Vanderburg CR, Russ C, et al. Mutations in the FUS/TLS gene on chromosome 16 cause familial amyotrophic lateral sclerosis. Science. 2009;323:1205–8.

7. Vance C, Rogelj B, Hortobágyi T, De Vos KJ, Nishimura AL, Sreedharan J, et al. Mutations in FUS, an RNA processing protein, cause familial amyotrophic lateral sclerosis type 6. Science. 2009;323:1208–11.

8. Gitcho MA, Baloh RH, Chakraverty S, Mayo K, Norton JB, Levitch D, et al. TDP-43 A315T mutation in familial motor neuron disease. Ann Neurol. 2008;63:535–8.

9. Kabashi E, Valdmanis PN, Dion P, Spiegelman D, McConkey BJ, Vande Velde C, et al. TARDBP mutations in individuals with sporadic and familial amyotrophic lateral sclerosis. Nat Genet. 2008;40:572–4.

10. Sreedharan J, Blair IP, Tripathi VB, Hu X, Vance C, Rogelj B, et al. TDP-43 mutations in familial and sporadic amyotrophic lateral sclerosis. Science. 2008;319:1668–72.

11. Manzano R, Toivonen JM, Moreno-Martínez L, de la Torre M, Moreno-García L, López-Royo T, et al. What skeletal muscle has to say in amyotrophic lateral sclerosis: Implications for therapy. Br J Pharmacol. 2021;178:1279–97.

12. Dobrowolny G, Aucello M, Rizzuto E, Beccafico S, Mammucari C, Boncompagni S, et al. Skeletal muscle is a primary target of SOD1G93A-mediated toxicity. Cell Metab. 2008;8:425–36.

13. Wong M, Martin LJ. Skeletal muscle-restricted expression of human SOD1 causes motor neuron degeneration in transgenic mice. Hum Mol Genet. 2010;19:2284–302.

14. Baskin KK, Winders BR, Olson EN. Muscle as a and quot;Mediator and quot; of systemic metabolism. Cell Metab. 2015;21:237–48.

15. Reyes ET, Perurena OH, Festoff BW, Jorgensen R, Moore WV. Insulin resistance in amyotrophic lateral sclerosis. J Neurol Sci. 1984;63:317–24.

16. Dedic SIK, Stevic Z, Dedic V, Stojanovic VR, Milicev M, Lavrnic D. Is hyperlipidemia correlated with longer survival in patients with amyotrophic lateral sclerosis? Neurol Res. 2012;34:576–80.

17. Shimizu T, Honda M, Ohashi T, Tsujino M, Nagaoka U, Kawata A, et al. Hyperosmolar hyperglycemic state in advanced amyotrophic lateral sclerosis. Amyotroph Lateral Scler. 2011;12:379–81.

18. Pradat P-F, Bruneteau G, Gordon PH, Dupuis L, Bonnefont-Rousselot D, Simon D, et al. Impaired glucose tolerance in patients with amyotrophic lateral sclerosis. Amyotroph Lateral Scler. 2010;11:166–71.

19. Hubbard RW, Will AD, Peterson GW, Sanchez A, Gillan WW, Tan SA. Elevated plasma glucagon in amyotrophic lateral sclerosis. Neurology. 1992;42:1532–4.

20. Hamasaki H, Takeuchi Y, Masui Y, Ohta Y, Abe K, Yoshino H, et al. Development of diabetes in a familial amyotrophic lateral sclerosis patient

carrying the I113T SOD1 mutation. Case Report Neuro Endocrinol Lett. 2015;36:414–6.

21. López-Otín C, Galluzzi L, Freije JMP, Madeo F, Kroemer G. Metabolic Control of Longevity. Cell. 2016;166:802–21.

22. Drory VE, Goltsman E, Reznik JG, Mosek A, Korczyn AD. The value of muscle exercise in patients with amyotrophic lateral sclerosis. J Neurol Sci. 2001;191:133–7.

23. Nieves JW, Gennings C, Factor-Litvak P, Hupf J, Singleton J, Sharf V, et al. Association Between Dietary Intake and Function in Amyotrophic Lateral Sclerosis. JAMA Neurol. 2016;73:1425–32.

24. Zhao Z, Sui Y, Gao W, Cai B, Fan D. Effects of diet on adenosine monophosphate-activated protein kinase activity and disease progression in an amyotrophic lateral sclerosis model. J Int Med Res. 2015;43:67–79.

25. Pikatza-Menoio O, Elicegui A, Bengoetxea X, Naldiz-Gastesi N, López de Munain A, Gerenu G, et al. The Skeletal Muscle Emerges as a New Disease Target in Amyotrophic Lateral Sclerosis. J Pers Med. 2021;11:671.

26. Scaricamazza S, Salvatori I, Ferri A, Valle C. Skeletal Muscle in ALS: An Unappreciated Therapeutic Opportunity? Cells. 2021;10:525.

27. Pascoe AL, Johnston AJ, Murphy RM. Controversies in TWEAK-Fn14 signaling in skeletal muscle atrophy and regeneration. Cell Mol Life Sci. 2020;

28. Enwere EK, Lacasse EC, Adam NJ, Korneluk RG. Role of the TWEAK-Fn14-cIAP1-NF-κB Signaling Axis in the Regulation of Myogenesis and Muscle Homeostasis. Front Immunol. 2014;5:34.

29. Sato S, Ogura Y, Tajrishi MM, Kumar A. Elevated levels of TWEAK in skeletal muscle promote visceral obesity, insulin resistance, and metabolic dysfunction. FASEB J. 2015;29:988–1002.

30. Boyer JG, Ferrier A, Kothary R. More than a bystander: the contributions of intrinsic skeletal muscle defects in motor neuron diseases. Front Physiol. 2013;4:356.

31. Crawford TO, Pardo CA. The neurobiology of childhood spinal muscular atrophy. Neurobiol Dis. 1996;3:97–110.

32. Meijboom KE, Sutton ER, McCallion E, McFall E, Anthony D, Edwards B, et al. Dysregulation of Tweak and Fn14 in skeletal muscle of spinal muscular atrophy mice. Skelet Muscle. 2022;12:18.

33. Bowerman M, Salsac C, Coque E, Eiselt É, Deschaumes RG, Brodovitch A, et al. Tweak regulates astrogliosis, microgliosis and skeletal muscle atrophy in a mouse model of amyotrophic lateral sclerosis. Hum Mol Genet. 2015;24:3440–56.

34. Mittal A, Bhatnagar S, Kumar A, Lach-Triflieff E, Wauters S, Li H, et al. The TWEAK-Fn14 system is a critical regulator of denervation-induced skeletal muscle atrophy in mice. J Cell Biol. 2010;188:833–49.

35. Dogra C, Hall SL, Wedhas N, Linkhart TA, Kumar A. Fibroblast growth factor inducible 14 (Fn14) is required for the expression of myogenic regulatory factors and differentiation of myoblasts into myotubes. Evidence for TWEAK-independent functions of Fn14 during myogenesis. J Biol Chem. 2007;282:15000–10.

36. Jakubowski A, Ambrose C, Parr M, Lincecum JM, Wang MZ, Zheng TS, et al. TWEAK induces liver progenitor cell proliferation. J Clin Invest. 2005;115:2330–40.

37. Rotarod Protocol. IMPReSS | International Mouse Phenotyping Resource of Standardised Screens. [cited 2013 Jul 26]. Available from: https://www mousephenotype.org/impress/ProcedureInfo?action=list&proclD=168.

38. Pfaffl MW. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001;29: e45.

39. Radonić A, Thulke S, Mackay IM, Landt O, Siegert W, Nitsche A. Guideline to reference gene selection for quantitative real-time PCR. Biochem Biophys Res Commun. 2004;313:856–62.

40. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012;9:676–82.

41. Bowerman M, Beauvais A, Anderson CL, Kothary R. Rho-kinase inactivation prolongs survival of an intermediate SMA mouse model. Hum Mol Genet. 2010;19:1468–78.

42. Sato S, Ogura Y, Kumar A. TWEAK/Fn14 Signaling Axis Mediates Skeletal Muscle Atrophy and Metabolic Dysfunction. Front Immunol. 2014;5:18.

43. Peterson JM, Bakkar N, Guttridge DC. NF-κB signaling in skeletal muscle health and disease. Curr Top Dev Biol. 2011;96:85–119.
44. Alhindi A, Boehm I, Chaytow H. Small junction, big problems: Neuro-muscular junction pathology in mouse models of amyotrophic lateral sclerosis (ALS). J Anat. 2022;241(5):1089-107.

45. Tintignac LA, Brenner H-R, Rüegg MA. Mechanisms Regulating Neuro-muscular Junction Development and Function and Causes of Muscle Wasting. Physiol Rev. 2015;95:809-52.

46. Bodine SC, Latres E, Baumhueter S, Lai VK, Nunez L, Clarke BA, et al. Identification of ubiquitin ligases required for skeletal muscle atrophy. Science. 2001;294:1704-8.

47. Misawa H, Hara M, Tanabe S, Niikura M, Moriwaki Y, Okuda T. Osteopontin is an alpha motor neuron marker in the mouse spinal cord. J Neurosci Res. 2012;90:732-42.

48. Morisaki Y, Niikura M, Watanabe M, Onishi K, Tanabe S, Moriwaki Y, et al. Selective Expression of Osteopontin in ALS-resistant Motor Neurons is a Critical Determinant of Late Phase Neurodegeneration Mediated by Matrix Metalloproteinase-9. Sci Rep. 2016;6:27354.

49. Oda Y, Nakanishi I. The distribution of cholinergic neurons in the human central nervous system. Histol Histopathol. 2000;15:825-34.

50. Kang Y, Saito M, Toyoda H. Molecular, Morphological and Electrophysiological Differences between Alpha and Gamma Motoneurons with Special Reference to the Trigeminal Motor Nucleus of Rat. Int J Mol Sci. 2024;25:5266.

51. Assessment of Motor Coordination and Balance in Mice Using the Rotarod, Elevated Bridge, and Footprint Tests - Brooks - 2012 - Current Protocols in Mouse Biology - Wiley Online Library. [cited 2021 Sep 21]. Available from: https://currentprotocols.onlinelibrary.wiley.com/doi/10.1002/9780470942390.mo110165.

52. Saxena S, Cabuy E, Caroni P. A role for motoneuron subtype-selective ER stress in disease manifestations of FALS mice. Nat Neurosci. 2009;12:627-36.

53. Liu KX, Edwards B, Lee S, Finelli MJ, Davies B, Davies KE, et al. Neuron-specific antioxidant OXR1 extends survival of a mouse model of amyotrophic lateral sclerosis. Brain. 2015;138:1167-81.

54. Pette D, Staron RS. Myosin isoforms, muscle fiber types, and transitions. Microsc Res Tech. 2000;50:500-9.

55. Plotkin DL, Roberts MD, Haun CT, Schoenfeld BJ. Muscle Fiber Type Transitions with Exercise Training: Shifting Perspectives. Sports (Basel). 2021;9:127.

56. Ratajczak W, Atkinson SD, Kelly C. The TWEAK/Fn14/CD163 axis-implications for metabolic disease. Rev Endocr Metab Disord. 2022;23:449-62.

57. Brown SAN, Cheng E, Williams MS, Winkles JA. TWEAK-independent Fn14 self-association and NF-kB activation is mediated by the C-terminal region of the Fn14 cytoplasmic domain. PLoS ONE. 2013;8: e65248.

58. Burkly LC, Michaelson JS, Hahm K, Jakubowski A, Zheng TS. TWEAKing tissue remodeling by a multifunctional cytokine: role of TWEAK/Fn14 pathway in health and disease. Cytokine. 2007;40:1-16.

59. Winkles JA. The TWEAK-Fn14 cytokine-receptor axis: discovery, biology and therapeutic targeting. Nat Rev Drug Discov. 2008;7:411-25.

60. Burkly LC, Dohi T. The TWEAK/Fn14 pathway in tissue remodeling: for better or for worse. Adv Exp Med Biol. 2011;691:305-22.

61. Scaricamazza S, Nesci V, Salvatori I, Fenili G, Rosina M, Gloriani M, et al. Endurance exercise has a negative impact on the onset of SOD1-G93A ALS in female mice and affects the entire skeletal muscle-motor neuron axis. Front Pharmacol. 2024;15:1360099.

62. Ortega-Hombrados L, Molina-Torres G, Galán-Mercant A, Sánchez-Guerrero E, González-Sánchez M, Ruiz-Muñoz M. Systematic Review of Therapeutic Physical Exercise in Patients with Amyotrophic Lateral Sclerosis over Time. Int J Environ Res Public Health. 2021;18:1074.

63. Tsitkanou S, Della Gatta P, Foletta V, Russell A. The Role of Exercise as a Non-pharmacological Therapeutic Approach for Amyotrophic Lateral Sclerosis: Beneficial or Detrimental? Front Neurol. 2019;10:783.

64. Carreras I, Yuruker S, Aytan N, Hossain L, Choi J-K, Jenkins BG, et al. Moderate exercise delays the motor performance decline in a transgenic model of ALS. Brain Res. 2010;1313:192-201.

65. Kirkinezos IG, Hernandez D, Bradley WG, Moraes CT. Regular exercise is beneficial to a mouse model of amyotrophic lateral sclerosis. Ann Neurol. 2003;53:804-7.

66. Hargreaves M, Spriell LL. Skeletal muscle energy metabolism during exercise. Nat Metab. 2020;2:817-28.

67. Quessada C, Bouscary A, René F, Valle C, Ferri A, Ngo ST, et al. Skeletal Muscle Metabolism: Origin or Prognostic Factor for Amyotrophic Lateral Sclerosis (ALS) Development? Cells. 2021;10:1449.

68. Desseille C, Deforges S, Biondi O, Houdebine L, D'amico D, Lamazière A, et al. Specific Physical Exercise Improves Energetic Metabolism in the Skeletal Muscle of Amyotrophic-Lateral-Sclerosis Mice. Front Mol Neurosci. 2017;10:332.

69. Lanfranconi F, Ferri A, Corna G, Bonazzi R, Lunetta C, Silani V, et al. Inefficient skeletal muscle oxidative function flanks impaired motor neuron recruitment in Amyotrophic Lateral Sclerosis during exercise. Sci Rep. 2017;7:2951.

70. Mezzani A, Pisano F, Cavalli A, Tommasi MA, Corrà U, Colombo S, et al. Reduced exercise capacity in early-stage amyotrophic lateral sclerosis: Role of skeletal muscle. Amyotroph Lateral Scler. 2012;13:87-94.

71. Trappe S, Luden N, Minchev K, Raue U, Jemiolo B, Trappe TA. Skeletal muscle signature of a champion sprint runner. J Appl Physiol. 1985;15(118):1460-6.

72. Cui D, Drake JC, Wilson RJ, Shute RJ, Lewellen B, Zhang M, et al. A novel voluntary weightlifting model in mice promotes muscle adaptation and insulin sensitivity with simultaneous enhancement of autophagy and mTOR pathway. FASEB J. 2020;34:7330-44.

73. Raue U, Jemiolo B, Yang Y, Trappe S. TWEAK-Fn14 pathway activation after exercise in human skeletal muscle: insights from two exercise modes and a time course investigation. J Appl Physiol. 1985;2015(118):569-78.

74. Ribeiro MBT, Guzzoni V, Hord JM, Lopes GN, Marqueti R de C, de Andrade RV, et al. Resistance training regulates gene expression of molecules associated with intramyocellular lipids, glucose signaling and fiber size in old rats. Sci Rep. 2017;7:8593.

75. Tomaz da Silva M, Joshi AS, Koike TE, Roy A, Mathukumalli K, Sopariwala DH, et al. Targeted ablation of Fn14 receptor improves exercise capacity and inhibits neurogenic muscle atrophy. FASEB J. 2022;36:e22666.

76. Sato S, Ogura Y, Mishra V, Shin J, Bhatnagar S, Hill BG, et al. TWEAK promotes exercise intolerance by decreasing skeletal muscle oxidative phosphorylation capacity. Skelet Muscle. 2013;3:18.

77. Rahmati M, Malakoutinia F. Aerobic, resistance and combined exercise training for patients with amyotrophic lateral sclerosis: a systematic review and meta-analysis. Physiotherapy. 2021;113:12-28.

78. Okamoto T, Torii S, Machida S. Differential gene expression of muscle-specific ubiquitin ligase MAFbx/Atrogin-1 and MuRF1 in response to immobilization-induced atrophy of slow-twitch and fast-twitch muscles. J Physiol. 2011;61:537-46.

79. Nishimura Y, Chunthorng-Om J, Lord S, Musa I, Dawson P, Holm L, et al. Ubiquitin E3 ligase Atrogin-1 protein is regulated via the rapamycin-sensitive mTOR-S6K1 signaling pathway in C2C12 muscle cells. Am J Physiol Cell Physiol. 2022;323:C215-25.

80. West DWD, Burd NA, Churchward-Venne TA, Camera DM, Mitchell CJ, Baker SK, et al. Sex-based comparisons of myofibrillar protein synthesis after resistance exercise in the fed state. J Appl Physiol. 1985;2012(112):1805-13.

81. Skelly LE, Gillen JB, MacInnis MJ, Martin BJ, Safdar A, Akhtar M, et al. Effect of sex on the acute skeletal muscle response to sprint interval exercise. Exp Physiol. 2017;102:354-65.

82. Taillandier D, Polge C. Skeletal muscle atrogenes: From rodent models to human pathologies. Biochimie. 2019;166:251-69.

83. Vainshtein A, Sandri M. Signaling Pathways That Control Muscle Mass. Int J Mol Sci. 2020;21:4759.

84. Tomaz da Silva M, Joshi AS, Castillo MB, Koike TE, Roy A, Gunaratne PH, et al. Fn14 promotes myoblast fusion during regenerative myogenesis. Life Sci Alliance. 2023;6:e202302312.

85. Enwere EK, Holbrook J, Lejmi-Mrad R, Vineham J, Timusk K, Sivaraj B, et al. TWEAK and cIAP1 regulate myoblast fusion through the noncanonical NF-kB signaling pathway. Sci Signal. 2012;5:ra75.

86. Rosa-Caldwell ME, Greene NP. Muscle metabolism and atrophy: let's talk about sex. Biol Sex Differ. 2019;10:43.

87. Pataky MW, Dasari S, Michie KL, Sevits KJ, Kumar AA, Klaus KA, et al. Impact of biological sex and sex hormones on molecular signatures of skeletal muscle at rest and in response to distinct exercise training modes. Cell Metab. 2023;35:1996-2010.e6.

88. Hindi SM, Mishra V, Bhatnagar S, Tajrishi MM, Ogura Y, Yan Z, et al. Regulatory circuitry of TWEAK-Fn14 system and PGC-1α in skeletal muscle atrophy program. FASEB J. 2014;28:1398-411.
89. Schiaffino S, Reggiani C. Fiber Types in Mammalian Skeletal Muscles. Physiol Rev. 2011;91:1447–531.

90. Goldspink G. Gene expression in muscle in response to exercise. J Muscle Res Cell Motil. 2003;24:121–6.

91. Willoughby DS, Nelson MJ. Myosin heavy-chain mRNA expression after a single session of heavy-resistance exercise. Med Sci Sports Exerc. 2002;34:1262–9.

92. Kichev A, Baburamani AA, Vontell R, Gressens P, Burkly L, Thornton C, et al. TWEAK Receptor Deficiency Has Opposite Effects on Female and Male Mice Subjected to Neonatal Hypoxia-Ischemia. Front Neurol. 2018;9:230.

93. Nuzzo JL. Narrative Review of Sex Differences in Muscle Strength, Endurance, Activation, Size, Fiber Type, and Strength Training Participation Rates, Preferences, Motivations, Injuries, and Neuromuscular Adaptations. The Journal of Strength & Conditioning Research. 2023;37:494.

94. Nuzzo JL. Sex differences in skeletal muscle fiber types: A meta-analysis. Clin Anat. 2024;37:81–91.

95. Jh V, Pr B, Ea J, M O, Jh W, Lh van den B. Sexual differences in onset of disease and response to exercise in a transgenic model of ALS. Neuro-muscular disorders : NMD. 2003 [cited 2024 Sep 26];13. Available from: https://pubmed.ncbi.nlm.nih.gov/14561497/

**Publisher's Note**

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.