Response of plasma microRNAs to nusinersen treatment in patients with SMA

Abstract Objective Spinal muscular atrophy (SMA) is a common genetic cause of infant mortality. Nusinersen treatment ameliorates the clinical outcome of SMA, however, some patients respond well, while others have limited response. We investigated microRNAs in blood samples from SMA patients and their response to nusinersen treatment evaluating the potential of circulating microRNAs as biomarkers for SMA. Methods In a discovery cohort study, microRNA next‐generation sequencing was performed in blood samples from SMA patients (SMA type 2, n = 10; SMA type 3, n = 10) and controls (n = 7). The dysregulated microRNAs were further analysed in the therapeutic response cohort comprised of SMA type 1 patients (n = 22) who had received nusinersen treatment, at three time points along the treatment course (baseline, 2 and 6 months of treatment). The levels of the studied microRNAs were correlated to the SMA clinical outcome measures. Results In the discovery cohort, 69 microRNAs were dysregulated between SMA patients and controls. In the therapeutic response cohort, the baseline plasma levels of miR‐107, miR‐142‐5p, miR‐335‐5p, miR‐423‐3p, miR‐660‐5p, miR‐378a‐3p and miR‐23a‐3p were associated with the 2 and 6 months response to nusinersen treatment. Furthermore, the levels of miR‐107, miR‐142‐5p, miR‐335‐5p, miR‐423‐3p, miR‐660‐5p and miR‐378‐3p at 2 months of treatment were associated with the response after 6 months of nusinersen treatment. Interpretation Blood microRNAs could be used as biomarkers to indicate SMA patients’ response to nusinersen and to monitor the efficacy of the therapeutic intervention. In addition, some of these microRNAs provide insight into processes involved in SMA that could be exploited as novel therapeutic targets.


Introduction
Spinal muscular atrophy (SMA) is an autosomal recessive neuromuscular disorder, characterised by progressive muscle atrophy and paralysis, resulting from motor neuron degeneration in the spinal cord. 1 SMA is caused by bi-allelic loss of function of the survival motor neuron 1 gene (SMN1). 2 Its paralogous neighbouring centromeric gene, SMN2, is intact in all patients but contains a single nucleotide variation in exon 7. The variation affects a splice enhancer and determines the exclusion of exon 7 in the majority of its transcripts, leading to an unstable SMN protein that cannot substitute the function of fulllength SMN1. 3,4 Nusinersen is an intrathecally administered antisense oligonucleotide drug and is the first approved diseasemodifying treatment for SMA. It increases full-length SMN protein by augmenting SMN2 exon 7 splicing in the mature mRNA transcripts. [5][6][7] While the extent of clinical response observed in SMA patients is unequivocal and life changing compared to the natural history of the disease, it is increasingly evident that patients respond differently to nusinersen treatment. Some affected children respond well to the acquisition of unexpected motor abilities, while others have a more limited clinical response. 5 Variables, including disease severity and the interval between disease onset and the start of treatment, are reported to play a role, but they do not account for the entire spectrum of clinical responses. 8,9 It is therefore important to identify informative biomarkers which may be used to better understand the variables that determine the response to nusinersen treatment.
MicroRNAs are a class of small (~22 nt) endogenous non-protein-coding RNA molecules that posttranscriptionally regulate gene expression. They are emerging biomarkers for diagnosis and prognosis of diseases, as well as treatment response, in a number of disease fields. 10,11 MicroRNAs are known to play important roles in the regulation of muscle and central nervous system development. [12][13][14] Dysregulation of microRNAs has been extensively investigated in animal models and patients with different neuromuscular disorders such as Duchenne muscular dystrophy and SMA. [15][16][17][18] Furthermore, disrupting microRNA biogenesis pathways by deleting Dicer 1 from spinal motor neurons in mice caused an SMA-like neurodegenerative phenotype with motor neuron dysfunction and death, indicating the close connection between microRNAs and motor neuron health. 19 Specific microRNA signatures in biofluids can provide useful non-invasive tools to monitor disease progression and response to therapeutic intervention. We have previously reported the potential of serum microRNAs as noninvasive and informative biomarkers of disease progression and response to experimental antisense therapy in SMA transgenic mice. 20 MiR-9, miR-206 and miR-132 were differentially expressed in the spinal cord, skeletal muscle and serum samples in SMA mice. After a singledose therapeutic morpholino antisense oligomer treatment in SMA mice, the pathological miR-132 levels in the spinal cord, muscle and serum reversed to the normal levels. In addition, altered levels of miR-9 and miR-132 were detected in serum samples from SMA patients.
To further assess the role of microRNAs as biomarkers in SMA, we here conducted a deep microRNA discovery profiling study using next-generation sequencing (NGS) in blood samples from SMA patients. We identified a list of microRNAs with significantly differential expression between the SMA patient cohort and healthy controls. The microRNA signatures were further studied in nusinersen-treated SMA patients to assess the effect that therapeutic intervention had on their profile. The molecular data were correlated with the longitudinal clinical response to therapy measured by motor function scores.
Our data reveals the differential expression of microRNAs in response to nusinersen treatment in individuals with SMA and underlines the potential application of selected circulating microRNAs as biomarkers in this condition.

Patient cohort
For the discovery cohort study, serum samples from SMA patients and healthy controls were obtained from the MRC Center for Neuromuscular Diseases Biobank London (REC reference number 06/Q0406/33, http://www. cnmd.ac.uk). A summary of the patients is presented in Table 1. All patients had confirmed genetic diagnosis of SMA with a homozygous genomic deletion in the SMN1 gene and two or three copies of the SMN2 gene. For the therapeutic response cohort study, blood samples for plasma extraction were collected from SMA patients (n = 22) treated with nusinersen at the Institute I-Motion, Hôpital Armand Trousseau, Paris, France and the CHR Citadelle Hospital, Li ege, Belgium. 21,22 All patients had a diagnosis of SMA type 1 and two or three SMN2 copies. The mean age of diagnosis of the patients was 3.6 months and the mean age of start of treatment was 65 months (median = 19 months; age range 1-420 months) ( Table 2). The samples were collected before the first intrathecal injection of nusinersen on day 1 as a baseline assessment, at 2 and 6 months of treatment. At each visit, the patients completed physiotherapy assessment using Hammersmith Infant Neurological Examination Section 2 (HINE-2) and Children's Hospital of Philadelphia Infant Test of Neuromuscular Disorders (CHOP-INTEND) clinical outcome assessment scales.

Serum microRNA next-generation sequencing
MicroRNAs were extracted from 500 lL of serum, and after microRNA library preparation, the libraries were sequenced on Illumina Nextseq500 with average of 10 million reads per sample, 50 nt singe-end reads.
The counts of microRNA were normalised by applying the trimmed mean of M-values method. 23 The analysis between the different groups was performed using a pairwise comparison of the normalised microRNA counts.
The p values were derived by an exact test on the negative binomial distribution and were corrected for multiple testing using the Benjamini-Hochberg method. 24 The Benjamini-Hochberg method, or false discovery rate (FDR), is the expected ratio of false-positive results to the total number of positive results. FDR is less stringent than Bonferroni correction and has a greater power of detecting positive results; however, it might lead to increased Type I errors. 25

miRCURY LNA PCR
Custom miRCURY microRNA panels were designed for the validation of the significant microRNAs from the NGS screening. MicroRNAs were extracted from 200 lL of plasma using miRNeasy Serum/Plasma Kit (cat no. 217184) and then converted to cDNA using an miRCURY LNA RT Kit (cat. no. 339340). After diluting the cDNA

Assessment of patients' motor function
The CHOP-INTEND 26 and the HINE 27 total scores were collected at baseline, 2 and 6 months after starting nusinersen treatment.

Statistical analysis
Univariate linear regression analysis was carried out to explore the relationship between levels of microRNAs and response to nusinersen treatment. Mann-Whitney test for independent samples or Wilcoxon-signed rank tests for the analysis of paired samples were used to analyse the change in the levels of microRNAs during treatment with nusinersen. The statistical analyses were performed using IBM SPSS Statistics V25 software and the p value of <0.05 was considered statistically significant. Graphs were created using IBM SPSS Statistics V25 software and Prism v8 GraphPad.

Results
Circulating microRNAs profiling in SMA patients by next-generation sequencing To identify circulating microRNAs differentially expressed between SMA patients and healthy controls, we carried out microRNA NGS in serum samples of SMA patients (type 2, n = 10; type 3, n = 10) and healthy controls (n = 7). On average 16.2 million reads were obtained from each sample. In the SMA patients, of the reads mapped to the human genome (hg19), on average 7% were mapped to microRNAs (mirBase release 20), and 36% were mapped to other small RNA sequences such as Y RNA, snRNA, snoRNA and tRNA. Forty-four percent of the reads were mapped to locations in the genome where no microRNAs or small RNAs were located, 5% of the reads mapped to abundant sequences, for example polyA and polyC homopolymers or ribosomal RNA, and 12% of the reads failed to map to the reference genome. In the control samples, a higher number of reads mapped to micro-RNA sequences (23%) and а comparable number of reads mapped to small RNA sequences (36%), other genomic regions (28%), abundant sequences (5%) and 12% did not align to the reference genome (Fig. 1A).
We then set to explore the expression pattern of micro-RNAs between the SMA and the control samples by performing hierarchical clustering. To visualise the microRNAs with the highest variation between SMA and control samples, we normalised the data by converting the expression values to a log2 scale and the microRNAs with the highest coefficient of variation were selected for the analysis. Three main clades of the dendrogram were observed with the largest being formed of SMA samples, followed by a clade of control samples and a clade of SMA patients' samples ( Fig. 1B). Next, we analysed the separation of SMA and control samples and the multidimensional scaling plot illustrated the subtle difference in the expression profiles of SMA and control samples as they did not cluster separately (Fig. 1C).
The differentially enriched serum microRNAs were first compared between SMA patients and controls and then between SMA type 2 and type 3 patients. A total of 42 microRNAs were differentially expressed between SMA and controls (FDR <0.05). Of these, 14 microRNAs were upregulated and 28 were downregulated in SMA patients (Table S1). When the serum microRNAs detected in SMA type 2 were compared to the microRNAs in type 3 SMA patients, a total of 27 microRNAs were significantly different (uncorrected p < 0.05), however, they did not pass the threshold for multiple testing of FDR <0.05. Of these micro-RNAs, 16 were upregulated and 11 were downregulated in SMA type 2 compared to type 3 patients (Table S2).
For the subsequent analysis in the therapeutic response cohort study, we selected the 42 microRNAs differentially expressed between SMA patients and healthy controls, and in addition, the 27 microRNAs with different levels between SMA type 2 and type 3, even though they did not pass the significance threshold for multiple testing.

Analysis of selected microRNAs in SMA patients treated with nusinersen
We next analysed the selected microRNAs in nusinersentreated SMA patients to explore the potential of these microRNAs to serve as biomarkers to monitor response to treatment.
Custom miRCURY LNA miRNA qPCR panels were designed for the analysis of the 69 microRNAs significant in the discovery study. Due to the limited availability of serum samples collected from nusinersentreated SMA patients, we carried out the screening in plasma samples. The plasma samples from nusinersentreated SMA patients were collected at three time points: pre-treatment (baseline), 2 and 6 months of treatment.
MicroRNAs with a cycle threshold (Ct) value above 36 in more than 25% of the samples was considered as not reliably detected and were excluded from the analysis. This left us with 41 microRNAs for further statistical analysis. Data were normalised to the expression of hsa-miR-191-5p which was detected as stably expressed across the different groups in the discovery study. 28

Prognostic microRNAs biomarkers in SMA type 1 patients treated with nusinersen
To understand the potential of these microRNAs as prognostic biomarkers for the response to therapeutic intervention, first, we analysed if the levels of microRNAs at baseline were associated with the response to nusinersen treatment assessed at 2 or 6 months of treatment. As a measure of the therapeutic response, we calculated the difference in CHOP-INTEND or HINE scores at 2 or 6 months from the baseline score for each patient. A larger difference in the scores indicates a better response to treatment with a higher number of test items achieved in CHOP-INTEND or HINE scales. Linear regression analysis between the baseline level of microRNAs and the change of scores in motor scales was then carried out. The baseline levels of six microRNAs (miR-107, miR-142-5p, miR-328-3p, miR-335-5p, miR-423-3p and miR-660-5p) showed significant association with the response at 2 months of treatment measured using CHOP-INTEND and four microRNAs (miR-181b-5p, miR-378a-3p, miR-125a-5p and miR-23a-3p) using HINE ( Fig. 2 and Table 3).
Further, we analysed if the relative expression of the studied microRNAs in plasma at 2 months of treatment correlated with the functional improvement of the patients after 6 months of treatment. When the relative expression of the microRNAs at 2 months and the change of CHOP-INTEND scores at 6 months to baseline were analysed, nine microRNAs were correlated with the functional improvement of patients at 6 months of nusinersen treatment. MiR-142-5p (R 2 = 0.51, p = 0.003) and miR-378a-3p (R 2 = 0.46, p = 0.006) showed the highest prediction potential ( Table 5). The linear regression plots are presented in Figure 4.
The relative expression of miR-340-5p (R 2 = 0.255, p = 0.04) at 2 months was the only significant microRNA associated with the motor function improvement after 6 months of treatment assessed using HINE (Fig. 4 and Table 5).
In addition, we carried out correlation analysis between the baseline levels of the analysed microRNAs and the functional improvement of the patients after  (Table S3).
We then set to examine if other factors such as the age of start of treatment and SMN2 copy number could have any effect on the plasma levels of the microRNAs. At baseline, the age of the SMA patients ranged between 1 month and 35 years (mean age = 5 years; median age = 3.5 years) with the majority of patients being between 1 month and 1.6 years (12 patients), followed by six patients with age between 2.3 and 4.3 years and four patients above 10 years. Linear regression analysis between the age at the start of nusinersen treatment and the baseline levels of the most promising microRNAs (miR-107, miR-142-5p, miR-335-5p, miR-423-3p, miR-660-5p, miR-378a-3p and miR-23a-3p) did not detect any significant positive association between the levels of microRNAs and the age of SMA patients except for miR-107, which was higher in older patients (R 2 = 0.27, p = 0.01) (Fig. S1). This result indicates that while most of the microRNAs identified in this study are stable with age, the impact of age on the levels of miR-107 should be considered when using miR-107 as a biomarker in SMA. No significant difference was detected between the levels of the microRNAs and SMN2 copy numbers.
The response of plasma microRNA levels to nusinersen treatment To determine whether the nusinersen treatment affected the plasma levels of the 41 microRNAs identified above, we carried out Wilcoxon matched-pairs t test comparing the microRNAs plasma levels at baseline, 2 and 6 months of treatment. Increased levels at 2 months compared to the baseline levels were detected in miR-335-5p (p = 0.02) and miR-328-3p (p = 0.04) (Fig. 5). Increased plasma levels at 6 months of treatment compared to the baseline levels were detected, in addition to miR-335-5p (p = 0.02), also for miR-423-3p (p = 0.002) and miR-142-5p (p = 0.03) (Fig. 5). When the levels after 6 months of treatment were compared to the levels at 2 months, a significant difference was detected only for miR-26b-5p (p = 0.04) which showed a decrease of miR-26b-5p at 6 months of nusinersen treatment (Fig. 5).
The differences between the baseline, 2 and 6 months levels of miR-335-5p, miR-328-3p, miR-423-3p, miR-142-5p and miR-26b-5p were also analysed in the nusinersentreated SMA patients, grouped based on their clinical response into two categories: low responders and good responders. As a measure of response, we used the difference in the assessment score at 2 or 6 months of treatment relative to the baseline score. A difference in the CHOP-INTEND score of at least 2 points or HINE score of at least 1 point was considered as a clinically significant improvement for the patients (good responders). A nonparametric test (Mann-Whitney test) was carried out. We did not detect any significant difference in the levels of the five microRNAs between patients with the low response and good response.
Further, to determine if the identified changes in the levels of the five microRNAs were due to the age of the patients rather than an effect of the nusinersen treatment, we conducted a linear regression analysis between the age of the patients and the baseline plasma microRNA levels. No significant correlation was detected between the age of SMA patients and the baseline levels of miR-335-5p, miR-328-3p, miR-423-3p, miR-142-5p and miR-26b-5p.

Discussion
Large heterogeneity in terms of clinical response to currently available treatments, ranging from the absence of response to impressive acquisition of motor milestones and survival, has been reported for all disease-modifying therapies in SMA, with disease duration negatively correlated with patient response. 29 Nevertheless, as other determinants of response are likely to exist, it is important to understand both the causes leading to the variable response as well as to identify biomarkers that could predict or monitor the response to treatment at an early stage.
Plasma levels of phosphorylated neurofilaments (NFs) have been previously used in SMA, 30,31 where elevated NF levels were detected in SMA type 1 patients compared to aged-matched controls. 30 The levels of NFs negatively correlated with the age of start of treatment and the severity of SMA patients. Moreover, nusinersen treatment was able to reduce NF levels, indicating that plasma NF levels could serve as biomarkers for disease activity and response to nusinersen treatment. 30 Investigating the blood NF levels in older SMA patients (above 1 year old), however, demonstrated that NF declined with age to a level comparable to healthy controls, 30 indicating the limited application of NFs as biomarkers in older SMA patients and pointing to the need of identifying additional biomarkers.
Muscle-specific microRNAs (myomiRs) in the blood and CSF from SMA patients have recently been investigated as biomarkers to monitor patients' response to nusinersen. 32,33 In the blood, miR-133a reduction predicted patients' response to nusinersen. 32 In the CSF, lower miR-206 and miR-133 predicted a more robust clinical response to nusinersen in later-onset SMA patients. 33 Here we have conducted a high-throughput microRNA NGS in blood samples of SMA patients and healthy controls to identify dysregulated microRNAs. Notably, a lower number of microRNAs were detected in SMA patients compared to controls (Fig. 1A). A hypothesis is   that the SMN complex affects the microRNA biogenesis by interacting with a specific RNA-binding protein 18,34,35 and in the SMA context the deficiency of SMN may cause defects in microRNA biogenesis leading to a lower number of microRNAs in SMA patients.
Although at the time of conducting the study in nusinersen-treated SMA 1 patients, a report demonstrated comparable levels of microRNAs in serum and plasma, 36 it is possible that some microRNAs might have been missed by analysing them only in plasma. Indeed, one of the microRNAs that we found significantly enriched in SMA serum compared to controls was miR-206 (>40 fold) (Table S1), which has been previously identified as dysregulated in SMA animal models. 20,37,38 However, in the therapeutic cohort study conducted in plasma, miR-206 was barely detectable and hence was not included in the following analysis. Discrepant microRNA levels between serum and plasma have also been reported previously. 39,40 Another factor that could have influenced the findings could be due to the different methods usedwhile NGS was used in the discovery cohort study, a Taq-Man qPCR assay was used in the therapeutic cohort study, and we cannot exclude different sensitivities in detecting miR-206 between these two techniques.
Different sets of microRNAs were identified in this study when the two motor function assessment scales, CHOP-INTEND and HINE, were analysed. The CHOP-INTEND scale has been designed to assess motor function in weak infants and the rapid progression of the disease. HINE motor milestones scale could longitudinally monitor the development of motor milestones, such as sitting, head control, crawling or standing, which are milestones that some nusinersen-treated SMA type 1 children can now acquire. 41 Thus, the different aspects of these two assessment tools could lead to capturing different microRNAs that monitor the distinct functional abilities of the patients. Indeed, fewer microRNAs were found significant in the HINE scale analysis compared to the CHOP-INTEND analysis. This is likely related to the fact that the CHOP-INTEND measures more granular aspects of function, while the HINE concentrates on more gross functional improvements. Future studies on more longitudinal time points beyond 10 months of treatment may reveal more microRNAs able to indicate the response measured using the HINE scale.
Among the seven microRNAs whose baseline levels showed potential in predicting patients' response to nusinersen treatment, miR-378a-3p was in common between both CHOP-INTEND and HINE scales. Elevated levels of miR-378a-3p were associated with better motor function in nusinersen-treated patients.
Pathway analysis of miR-107, miR-142-5p, miR-335-5p, miR-423-3p, miR-660-5p, miR-378a-3p and miR-23a-3p    identified three pathways involved in neurotrophin signalling (p = 0.0005,) mTOR signalling (p = 0.005) and FoxO signalling (p = 0.003) as enriched with genes targeted by the analysed microRNAs. These three pathways are all involved in the regulation of cell survival, autophagy and apoptosis. miR-107 has been reported to inhibit the apoptosis of myoblasts but also to negatively regulate myoblast differentiation by suppressing the expression of MyoD, MyoG and MyHC. 42 A counteracting role was reported from miR-378a-3p that promotes skeletal muscle differentiation. 43,44 As a dynamic system, some myogenic cells will go through self-renewal and return to quiescence, while others will proliferate and differentiate into myofibers, thus the two microRNAs likely reflect the different processes of skeletal muscle development.
A direct target of miR-142-5p is FOXO3, a transcriptional activator that induces the expression of autophagy genes, triggering proteolysis of skeletal muscles and activating apoptosis leading to neuronal cell death upon oxidative stress. 45,46 Overexpression of miR-142-5p significantly reduces the mRNA levels of FOXO3 and enhances the expression of genes involved in skeletal muscle growth. 47 In our study, we observed that higher baseline levels of miR-142-5p were associated with a better response to therapeutic intervention. We hypothesise that the increased levels of miR-142-5p, through suppression of FOXO3, may lead to reduced muscular atrophy and improved motor function acquirement.
Previously, miR-335-5p had been found to control selfrenewal or differentiation of mouse embryonic stem cells. 48 Decreased level of miR-335-5p was detected in SMAinduced pluripotent stem cells (iPSCs) compared to wildtype iPSCs during the early steps of differentiation. 49 In this study, we demonstrated that higher levels of miR-335-5p were associated with better motor function improvement in SMA patients. We also detected increased levels of miR-335-5p at 2 and 6 months after treatment. Our data may suggest a potential protective effect of miR-335 in SMA.
Worth noting is miR-23a-3p, which we found increased in nusinersen patients with the higher motor ability both at baseline and following therapeutic intervention. MiR-23a has been previously reported as having neuroprotective and skeletal muscle atrophy inhibiting properties 50 and has been identified as a protective modifier of SMA. 51 In the central nervous system, miR-23a inhibits the expression of the apoptotic protease activating factor-1 (Apaf-1) and thus regulates the sensitivity of neurons to apoptosis. 50 Overexpression of miR-23a has also been shown to promote the expression of myelin genes through directly regulating the PTEN gene. 52 In SMA iPSC-derived motor neurons miR-23a was found to be significantly reduced compared to iPSC-derived motor neurons from healthy individuals. 51 Transfection of synthetic miR-23a into the SMA iPSC-derived motor neurons led to reduced motor neuronal death, demonstrating that miR-23a deficiency contributes to the pathology observed in SMA motor neurons. 51 Furthermore, AAV9mediated expression of miR-23a in the Smn 2B /À SMA mouse model 53 resulted in reduced neuromuscular junction defects, improved motor neuronal soma size, increased muscle fibre area and extended the survival of the Smn 2B /À mice. 51 In skeletal muscle, overexpression of miR-23a downregulated MAFbx/atrogin-1 and MuRF1 genes and provided resistance against atrophy exemplified by preserved muscle weight and cross-sectional area of the muscles. 54 These data indicate miR-23a as an important microRNA to be characterised further in longer duration studies not only in patients receiving nusinersen but also in other SMA therapeutic interventions, including risdiplam and onasemnogene abeparvovec. In addition to representing a promising biomarker, this microRNA may also serve as an independent target for therapeutic intervention.
It is acknowledged that while the discovery cohort study included SMA type 2 and 3 patients, the therapeutic cohort only consisted of SMA type 1 patients who are more severe and younger than type 2 and type 3. In addition, while the median age in the discovery cohort was 10.2 years, the median age in the therapeutic cohort was only 19 months. The discrepancy in age and clinical severity between the two study cohorts might have   resulted in some microRNAs specific to the severe type 1 patients being missed as the selection of microRNAs was based on data obtained in milder patients in the discovery cohort. Therefore, a further study in young and severe type 1 patients might be needed in the future to compliment our findings. Another limitation of this study is the relatively small number of patients and the short duration of follow-up. While this is not uncommon for rare genetic diseases, our data need to be further validated in a larger cohort, and the actual added value of treatment response predictive value should be quantified after adjusting for age, disease duration, baseline CHOP-INTEND and weight/height. Additionally, it will be interesting to investigate these microRNAs in patients treated by other disease-modifying treatments such as risdiplam or onasemnogene abeparvovec.
In conclusion, we have identified a set of microRNAs in the blood of nusinersen-treated SMA patients that could be considered as potential biomarkers to predict or monitor patients' response to treatment. To further develop microRNAs in blood as biomarkers for SMA, more longitudinal studies in larger SMA cohorts will be needed to better understand the correlation of the microRNAs to disease progression and the improvement in motor function acquisition following different therapeutic interventions.

Supporting Information
Additional supporting information may be found online in the Supporting Information section at the end of the article. Figure S1. Linear regression analysis between the age of start of nusinersen treatment and the relative expression at baseline of miR-107. Table S1. Differentially expressed microRNA in serum of SMA patients compared to controls. Table S2. Differentially expressed microRNA (p < 0.05) in serum of SMA type 2 to SMA type 3 patients. Table S3. Correlation between the baseline levels of microRNAs and the functional improvement of the patients at 2 and 6 months after treatment assessed using CHOP-INTEND and HINE scales.