Preprint
Article

This version is not peer-reviewed.

Neuronal Pentraxin 2 as a Potential Biomarker for Nusinersen Therapy Response in Adults with Spinal Muscular Atrophy: A Pilot Study

A peer-reviewed article of this preprint also exists.

Submitted:

26 June 2025

Posted:

30 June 2025

You are already at the latest version

Abstract
Background: The treatment landscape for spinal muscular atrophy (SMA) has changed significantly with the approval of gene-based therapies such as nusinersen for adults with SMA (pwSMA). Despite their efficacy, high costs and treatment burden highlight the need for biomarkers to objectify or predict treatment response. This study aimed to identify such biomarkers. Methods: A proteomic analysis of cerebrospinal fluid (CSF) from pwSMA (n = 7), who either significantly improved (SMA Improvers) or did not improve in motor function (SMA Non-Improvers) under nusinersen therapy, was performed. Data are available via ProteomeXchange with identifier PXD065345. Candidate biomarkers—Neuronal Pentraxin 2 (NPTX2), Contactin 5 (CNTN5), and Anthrax Toxin Receptor 1 (ANTXR1)—were investigated by ELISA in serum and CSF from an independent pwSMA cohort (n = 14) at baseline, 2 and 14 months after therapy initiation. Biomarker concentrations were correlated with clinical outcomes. Additionally, NPTX2 was stained in spinal cord sections from a mild SMA mouse model (FVB.Cg-Smn1tm1Hung Tg(SMN2)2Hung/J). Results: CSF NPTX2 levels decreased in pwSMA after 14 months of nusinersen therapy, independent of clinical response. In serum, NPTX2 levels were lower in SMA Improvers compared to Non-Improvers at 14 months. CNTN5 and ANTXR1 showed no significant changes. In the SMA mouse model, NPTX2 immunoreactivity increased at motoneuron loss onset. Conclusion: NPTX2 emerges as a potential biomarker of treatment response to nusinersen in pwSMA suggesting its significant pathophysiological role in late-onset SMA, warranting further investigation.
Keywords: 
;  ;  ;  

1. Introduction

5q-associated Spinal Muscular Atrophy (SMA) is a hereditary motor neuron disease leading to progressive weakness and atrophy of limb, bulbar and respiratory muscles caused by a homozygous deletion or compound heterozygosity with deletion and point mutation in the survival motor neuron 1 (SMN1) gene [1]. Survival motor neuron (SMN) protein deficiency causes progressive motoneuronal degeneration in the anterior horn of the spinal cord and can only be insufficiently compensated by truncated SMN protein deriving from the survival motor neuron 2 (SMN2) gene of variable copy number [2,3]. The disease classification ranges from type 1, characterized by a low SMN2 copy number and a naturally early-onset, severe disease course with a shortened life expectancy, to type 4, featuring a higher SMN2 copy number, later disease onset in adulthood, and normal life expectancy [4,5].
Recently, disease-modifying therapies (nusinersen, risdiplam and onasemnogene abeparvovec) have been approved for treatment of SMA, with nusinersen and risdiplam showing improvement of motor function even in adults [6,7,8,9]. While these treatments are proven effective, therapy costs and efforts are high. Treatment decisions are currently not based on clear evidence-based criteria but on patient preference, unless, for example, technical conditions such as severe scoliosis or spondylodesis prohibit the intrathecal administration of nusinersen.
Considering expanding therapeutic options for SMA, it is essential to establish objective decision criteria for the selection of specific therapies. In this regard, the utilization of biomarkers could serve to objectively measure treatment success and predict treatment response. Although clinical, structural-morphological, and functional assessments are valuable, they are often challenging to assess, rater-dependent, and limited in their sensitivity [10,11]. Blood and cerebrospinal fluid (CSF) biomarkers are easier to obtain, quantitative, and can be evaluated objectively. Neurofilaments are becoming established blood and CSF biomarkers for Amyotrophic Lateral Sclerosis and Multiple Sclerosis [12,13]. They have also been investigated in SMA [14,15,16]. A recent large-scale study investigating Neurofilament light chain (NfL) levels longitudinally in 113 adults with SMA under nusinersen therapy found higher NfL levels in the SMA than in the control group at baseline, but no predictive value of NfL regarding clinical parameters [15]. A moderate decrease in NfL concentration in CSF was observed under nusinersen treatment, but no correlation with clinical parameters was found [15].
Proteomics enables the comprehensive analysis of proteins in biological samples, allowing for the detection of changes in protein expression or modification [17]. These changes can serve as potential biomarkers indicating treatment response at the molecular or even clinical level. To identify useful biomarkers for adult SMA, we quantified candidate biomarkers identified through CSF-based proteome analysis in adults with SMA who showed different responses to nusinersen therapy.

2. Materials and Methods

2.1. Study Design

This study was conducted in the Department of Neurology, University Hospital, Essen, Germany. Participants provided written informed consent prior to their inclusion in the study. The study was approved by the local ethics committee of the University of Duisburg-Essen, Germany (approval number 18-8071-BO).
We included adult individuals with molecularly confirmed (homozygous deletion of SMN1) 5q-associated SMA, referred to as persons with SMA (pwSMA), under treatment with nusinersen as well as healthy non-SMA individuals as controls (lumbar puncture for exclusion diagnostics and without inflammatory, degenerative, ischemic, neoplastic or metabolic diseases of the central nervous system). Assessments and samples of pwSMA (serum and CSF) were collected within the regular treatment schedule for administration of nusinersen between 2017 and 2021.
In a first approach, an unbiased proteomic CSF analysis from 7 pwSMA with different responses to nusinersen therapy, as well as from 2 controls, was conducted. In a second, hypothesis-driven approach, single proteins were analyzed through quantification by ELISA in the CSF and serum in an independent sample of 14 pwSMA and 4 controls. PwSMA who improved by at least 3 points in the Hammersmith Functional Motor Scale Expanded (HFMSE) Score after 6 months of nusinersen therapy are classified as “SMA Improvers”, while those who showed less or no improvement in the HFMSE are classified as “SMA Non-Improvers”. To further validate the biological significance of the CSF and serum findings in pwSMA, a small-scale analysis of spinal cord tissue from a late-onset SMA mouse model was conducted. An overview of the analyses performed and the samples included is provided in Figure 1.

2.2. Assessment of Motor Function

We used the HFMSE and Revised Upper Limb Module (RULM) to assess motor function. Tests were conducted within regular visits at baseline, 2 months, 6 months and 14 months after treatment initiation with nusinersen. HFMSE and RULM are both established and validated tools for evaluating motor function and disease progression in individuals with SMA types 2 and 3 [18,19]. The HFMSE assesses motor function of the extremities, trunk, and head with 33 items scored from 0–2, resulting in a maximum score of 66. Higher scores indicate better motor function. The RULM evaluates upper extremity function using 19 items. One item is scored 0 or 1, and the remaining 18 are scored 0–2 points resulting in a maximum score of 37 with higher scores indicating better motor function.

2.3. Proteome Analysis from CSF

Unbiased proteomic profiling was performed on CSF samples from 7 pwSMA and 2 age-matched controls. Demographic and clinical characteristics of pwSMA are summarized in Table 1. Sample composition and measurement time points (baseline and 6 months) are depicted in Figure 1.
For each sample, 50 µL of CSF were reduced, alkylated, and digested overnight with trypsin. The resulting peptides were desalted using C18 spin columns, dried, and reconstituted in LC solvent containing indexed retention time (iRT) peptides for calibration. Peptide concentrations were quantified via BCA assay. Equal peptide volumes per group and time point were pooled and fractionated by high-pH reversed-phase chromatography (HPRP) into 6 fractions. LC-MS/MS analysis was performed on a Thermo Scientific™ Q Exactive™ HF mass spectrometer coupled to an Easy-nLC 1200 nano-LC system. For spectral library generation, pooled study samples and commercial CSF were analyzed using data-dependent acquisition (DDA). Data-independent acquisition (DIA, HRM) with 22 windows per sample was used for quantification. All raw data were processed using Spectronaut Pulsar (Biognosys), employing a hybrid spectral library generated from DDA and DIA data. Protein and peptide identifications were filtered at a 1% false discovery rate (FDR). Data normalization was performed using local regression.
Differentially expressed proteins were defined by a p-value < 0.05 and an average fold change > 1.5. Further statistical analysis and data visualization were conducted in R. Hierarchical clustering was based on Manhattan distance and Ward’s linkage. Principal component analysis (PCA) was performed using prcomp, and partial least squares discriminant analysis (PLS-DA) using the mixOmics package.
The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD065345. Data can be assessed by logging in to the PRIDE website using the following account details: Username: reviewer_pxd065345@ebi.ac.uk; Password: mXK9pLXENO4D.

2.4. Sample Selection for Enzyme-Linked Immunosorbent Assay (ELISA)

To reduce the risk of sampling bias and to ensure a broad spectrum of treatment responses, an independent sample was created. From a total of 28 pwSMA with available baseline CSF or serum samples, 14 were selected to represent distinct ranges of change in HFMSE score (Δ6 HFMSE) after 6 months of nusinersen treatment, covering the full observed response range (Δ6 HFMSE -6 to +20). One individual per response level was randomly selected to ensure maximal variability in treatment outcomes, rather than representativeness of the cohort. HFMSE and RULM scores were available for all pwSMA at all time points. Demographic and clinical characteristics of all included pwSMA are given in Table 1. Four non-SMA individuals served as controls.
Table 1. Demographic and clinical characteristics of all pwSMA included in ELISA (n=14) and proteome (n=7) analysis.
Table 1. Demographic and clinical characteristics of all pwSMA included in ELISA (n=14) and proteome (n=7) analysis.
Sex Male (%) Female (%)
Proteome analysis 5 (71) 2 (29)
ELISA Analysis 10 (71) 4 (29)
SMA type 2 3
Proteome analysis 2 (29) 5 (71)
ELISA Analysis 1 (7) 13 (93)
SMN2 copy number 2 (%) 3 (%) 4 (%)
Proteome analysis 0 (0) 3 (43) 4 (57)
ELISA Analysis 1 (7) 6 (43) 7 (50)
Clinical classification Non-Sitter Sitter Walker
Proteome analysis 0 (0) 4 (57) 3 (43)
ELISA Analysis 4 (29) 1 (7) 9 (64)
Spondylodesis Yes (%) No (%)
Proteome analysis 1 (14) 6 (86)
ELISA Analysis 2 (14) 12 (86)
Non-invasive ventilation Yes (%) No (%)
Proteome analysis 2 (29) 5 (71)
ELISA Analysis 3 (21) 11 (79)
Age Mean (sda) Minimum Maximum
Proteome analysis 40.00 (14.28) 26 65
ELISA Analysis 37.71 (12.86) 22 62
Disease duration in years Mean (sd) Minimum Maximum
Proteome analysis 35.01 (16.22) 14 63
ELISA Analysis 29.08 (14.34) 12 61
Baseline HFMSE score Mean (sd) Minimum Maximum
Proteome analysis 28.00 (22.06) 0 57
ELISA Analysis 32.50 (21.38) 1 62
Baseline RULM score Mean (sd) Minimum Maximum
Proteome analysis 25.00 (13.90) 0 37
ELISA Analysis 29.64 (10.49) 7 37
a standard deviation.

2.5. ELISA from CSF and Serum

In a hypothesis-driven approach, single proteins were selected for quantification by ELISA in the cohort described in section 2.4. Based on the results of the proteome analysis and the function of these proteins, Contactin 5 (CNTN5), mediating cell surface interactions during the development of the nervous system, and Anthrax Toxin Receptor 1 (ANTXR1), a transmembrane protein regulating cell attachment and migration, were further investigated. As the proteome analysis also revealed regulatory differences of proteins of the neuronal pentraxin family, Neuronal Pentraxin 2 (NPTX2) was additionally selected for further investigation. This decision was supported by previous literature reports describing NPTX2 as a biomarker candidate for synaptic impairment in neurodegenerative diseases with altered expression reported in conditions such as Alzheimer's disease and amyotrophic lateral sclerosis [20,21].
The levels of NPTX2, CNTN5 and ANTXR1 were measured in CSF and serum at 3 time points: baseline (T0), 2 months (T1) and 14 months (T2) after starting treatment with nusinersen. Measuring time points and types and numbers of samples that were analyzed are shown in Figure 1. To determine the levels of NPTX2, CNTN5 and ANTXR1 in CSF and serum samples, ELISA assays were performed according to the manufacturer’s protocol (#MBS8802191, #MBS1606034, #MBS941904; MyBioSource, USA).

2.6. Animals

The SMN-deficient mouse model FVB.Cg-Smn1tm1Hung Tg(SMN2)2Hung/J (Jackson #005058) reflects a milder SMA phenotype. The animals exhibit reduced weight after 20 days and decreased grip strength after 33 days, while maintaining a normal life expectancy. The model was chosen to reflect the disease course of milder affected pwSMA [22]. This mouse model was obtained from Jackson Laboratory and bred at the University Hospital Essen. It is homozygous for the murine SMN1 knockout and carries 4 copies of human SMN2. Mice were maintained on a 12/12-hour light/dark cycle with free access to water and standard food pellets. Animals were monitored weekly to examine body condition, weight, and general health. The experiments were conducted under the animal welfare guidelines of the University Duisburg Essen. The use of the SMA mouse model was approved by the State Agency for Nature, Environment and Consumer Protection (LANUV) in North Rhine-Westphalia (reference number 81–02.04.2020.A335).

2.7. NPTX2 Immunostaining of Murine Spinal Cord Tissue Slices

Spinal cord tissue was collected from three mice of the SMA mouse model, referred to as “SMA mice”, each at P20 (postnatal day 20), P42, P70 and P>400. The same number of age-matched FVB/N wild type (wt) mice served as controls.
Lumbar spinal cord tissue was frozen in liquid nitrogen and stored at -80 °C before cryosections of 20 µm were prepared, with every fifth section placed on a separate microscopy slide. The tissue was fixed in 4% paraformaldehyde (PFA in PBS, 15 min), washed, permeabilized (PBS, 0.1 v/w Triton X-100, 15 min) and blocked (PBS, 5% bovine serum albumin, 1h). Primary antibodies for motor neurons (anti-SMI-32, mouse, 1:500, #801701, BioLegends, San Diego, CA, USA) and NPTX2 (anti-NPTX2, rabbit, 1:250, #NBP2-19572, Novus Biologicals, Centennial, CO 80112, USA) were diluted in blocking solution and incubated at 4 °C overnight. Sections were washed and secondary antibodies (goat anti-rabbit, goat anti-mouse, 1:300, Dianova, Hamburg, Germany) and DAPI (1:1000, Sigma-Aldrich, Taufkirchen, Germany) were diluted in blocking solution and incubated at room temperature.
Three slices per animal were analyzed. Images were captured with a Zeiss Axio Observer.Z1 Apotome fluorescence microscope and Zeiss Zen software. Immunoreactivity was quantified using ImageJ software (NIH). Immunoreactive cells were selected using the free-hand tool and the fluorescence intensity of NPTX2 was measured and normalized to the background in each image. Fluorescence intensity in SMA mice tissue was normalized to that of wt tissue. Following immunostaining with the motoneuronal marker SMI-32, motoneurons were manually counted.

2.8. Statistical Analysis

Statistical analyses were performed on the results of protein quantification by ELISA and clinical scores. Additionally, NPTX2 immunostaining and motoneuron count in spinal cord tissue of the SMA mice were analyzed.
A Wilcoxon signed-rank test was performed for pre-post comparisons of the biomarker levels between T0, T1 and T2. The Mann-Whitney U test was used to compare SMA Improvers and SMA Non-Improvers. Correlations between the protein levels and motor scores HFMSE and RULM were calculated using Spearman’s rank correlation coefficient. Calculations were performed using blank-subtracted values, which were adjusted by subtracting the absorbance value of the blank (a sample without the target analyte) from the actual sample value. This helps eliminate background noise and nonspecific signals, ensuring more accurate results. The absorbance values of pwSMA were normalized to those of the controls.
To compare the NPTX2 immunostaining fluorescence intensity and number of motoneurons in the spinal cord slices between SMA and wt mice at P20, P42, P70, and P>400, the one-sample Wilcoxon signed-rank test was used. For each time point, 6 animals (three SMA mice and three wt mice) were examined, with up to 3 spinal cord sections per animal. The statistical analysis was performed based on the values per spinal cord section, which were normalized to the wt values on the respective slide. If a section or the staining was not fully intact, it was excluded from the analysis.

3. Results

3.1. Proteomic Profiling of CSF

A total of 117 proteins were found to be differentially regulated (p-value <0.05 and average fold change >1.5) between pwSMA and controls. Figure 2 shows the differences in protein expression between pwSMA and controls. Seventeen proteins (ANTXR1, CNTN5, CSF2RA, C1QTNF5, FUCA2, HAMP, IGHG3, IGHV3-43D, IGHV3-74, IGHV6-1, IGLC7, LARGE1, MT3, MYOC, PLOD3, RTN1, SERPINE2) were differentially regulated between SMA Improvers and SMA Non-Improvers at baseline. Between baseline and 6 months after treatment initiation, 4 proteins (DNS2A_HUMAN, FSTL4, NPTN, PLA2G15) were differentially regulated in SMA Improvers and 3 (CACHD1, DDAH1, HEXB) in SMA Non-Improvers. Among the proteins differentially regulated between SMA Improvers and SMA Non-Improvers at baseline were CNTN5 (ratio 0.21; log ratio 2.24; p-ANOVA = 0.01) and ANTXR1 (ratio 1.91; log ratio 0.93; p-ANOVA = 0.05) (Figure 2). NPTX2 levels were not different between SMA Improvers and SMA Non-Improvers at baseline (ratio 1.21; log ratio 0.28; p-ANOVA = 0.48). Regarding other proteins of the NPTX family, NPTX1 was lower at 6 months than at baseline, measured across all pwSMA (ratio 1.3; log ratio 0.38; p-ANOVA = 0.03). All protein ratios between the groups are available as a supplementary online resource (Online Resource 1).

3.2. ELISA-Based Quantification of NPTX2, CNTN5 and ANXTR1 in CSF and Serum

In CSF, the NPTX2 concentration showed no differences between pwSMA and controls (Figure 3a). Across all pwSMA, the NPTX2 concentration in CSF was decreased at T2 compared to T0 (z=2.666, p=0.008, n=9) and to T1 (z=2.366, p=0.018, n=7) (Figure 3a). There were no differences between SMA Improvers and SMA Non-Improvers at T0, T1 or T2 (Figure 3b) and no correlations of change in the NPTX2 CSF level and HFMSE or RULM scores from T0 to T2.
In serum, the NPTX2 concentration showed no differences between pwSMA and controls (Figure 3c), nor between SMA Improvers and SMA Non-Improvers at T0 or T1 (Figure 3d). NPTX2 serum levels did not change over time across all pwSMA (Figure 3c). The NPTX2 serum level was higher at T2 in SMA Non-Improvers than in SMA Improvers (z=-2.030, p=0.48, n=12) (Figure 3d). It increased in all 5 SMA Non-Improvers from T0 to T2 and decreased in 5 of 6 SMA Improvers from T0 to T2. Changes in the NPTX2 serum level from T0 to T2 correlated negatively with changes in the HFMSE from T0 to T2 (r=-0.688, p=0.019, n=11) (Figure 4).
The CNTN5 serum but not CSF level was higher in controls than in pwSMA at baseline (z= 2.219, p=0.026, n=15). There were no changes in any other CNTN5 and ANTXR1 levels over time or differences between controls and SMA Improvers and SMA Non-Improvers. Levels of CNTN5 and ANTXR1 did not correlate with motor scores.

3.3. NPTX2 Immunostaining of Murine Spinal Cord Tissue Slices

The NPTX2 immunostaining fluorescence intensity was higher in SMA mice compared to wt mice at P42 (z=2.366, p=0.018, n=7; 3 animals per condition with 2 or 3 slices per animal; wt mice not included in the analysis as they serve as reference values). NPTX2 intensities did not differ between SMA and wt mice at P20, P70 or P>400. The motoneuron count was significantly lower in SMA mice compared to wt mice at P42 (z=-2.375, p=0.018, n=7), P70 (z=-2.366, p=0.018, n=7) and P>400 (z=-2.677, p=0.007, n=9) but not at P20, indicating motoneuronal loss had begun by P42 (Figure 5).

4. Discussion

The level of NPTX2, a protein of the neuronal pentraxin family involved in synaptic function and plasticity, decreased in CSF in adult pwSMA in the first 14 months of treatment with nusinersen independent of motor improvement. The change in NPTX2 serum levels over 14 months of nusinersen treatment correlated with the change in HFMSE during this period. This demonstrates the potential of NPTX2 as a CSF biomarker indicating a biological response to nusinersen therapy and a serum biomarker displaying clinical treatment response. The upregulation of NPTX2 in the spinal cord with the onset of motor neuron loss in the SMA mouse model supports the biological significance of the CSF and serum findings.
NPTX2 is expressed in the brain and spinal cord but also in peripheral tissues, such as the Langerhans islets and the adrenal medulla, and is involved in the assembly of excitatory synapses [23]. Notably, abnormalities in the Langerhans islets have been reported in SMA mouse models and in children with severe infantile form of the disease, linking this peripheral expression site of NPTX2 to pathological processes in SMA [24,25].
The presynaptic expression and release of NPTX2 is upregulated by synaptic activity and also by brain-derived neurotrophic factor (BDNF) [26]. After being released, NPTX2 can accumulate at excitatory synapses on interneurons within perineuronal nets [27]. Its function there has not yet been fully clarified. However, building a complex with Neuronal Pentraxin 1 (NPTX1), it binds and aggregates α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) receptors [28,29]. It thus can modulate synaptic plasticity and in particular the transmission of excitatory signals via glutamatergic synapses and appears to stabilize AMPA receptors and thus increase their activity [30]. Under certain stress conditions, e.g. with upregulation of tumor necrosis factor-α converting enzyme, Neuronal Pentraxin Receptor (NPTXR) clusters with NPTX2 and AMPA receptors, leading to the internalization of this complex [31]. Investigations in an ischemic stroke rat model suggested that the induced clustering of AMPA receptors by NPTX2 leads to their internalization, resulting in reduced response to glutamate, reflecting a rescue-mechanism. Thus, upregulation of NPTX2 may protect neurons from glutamate excitotoxicity in environments with excess of glutamate in the synaptic cleft as in ischemic stroke [32], but also in adult SMA. Glutamate excitotoxicity has been shown to play a critical role in motor neuron loss in adult SMA. One of our previous studies demonstrated that the excitatory amino acid transporter 1 (EAAT1), a glutamate transporter expressed in astrocytes, was downregulated in the spinal cord of a late-onset SMA mouse model and in cultured SMN-deficient astrocytes [33]. Decreased EAAT1 expression was associated with increased glutamate levels in the spinal cord [33]. Another group validated the reduction of EAAT1 in lumbar spinal cord tissue of individuals with SMA and in human iPSC-derived astrocytes [34]. The activation of spinal astrocytes and the decreased expression of EAAT1 with increased glutamate levels even preceded the loss of spinal motor neurons in the late-onset SMA mouse model [22]. This early astrocytic pathomechanism of SMA may involve NPTX2 via its influence on glutamatergic synapses.
Considering the biological function of NPTX2, the decrease in NPTX2 CSF level during nusinersen therapy might reflect a downregulated rescue mechanism as a consequence of SMN restoring therapy. However, an upregulation of NPTX2 at baseline could not be demonstrated in pwSMA compared to controls. The divergent dynamics observed in CSF and serum remain unclear, and may, at least in part, be attributable to the limited sample size as well as to the cross-reactivity and limited sensitivity of the ELISA, particularly in the context of low protein concentrations in CSF. Furthermore, the dynamics of NPTX2 should be interpreted in light of the overall low disease dynamics in our cohort, predominantly consisting of SMA type 3, characterized by long disease duration and a slow rate of clinical progression.
Serum NPTX2 levels have been analyzed in previous studies and have been found to be altered in individuals with acute psychotic episodes in schizophrenia [35] and to be associated with cognitive function in patients with dementia [36]. Also, CSF NPTX2 levels as measured by ELISA, have been found to be reduced in individuals with Alzheimer´s Dementia compared to age-matched controls. NPTX2 concentration was even shown to correlate with cognitive performance and hippocampal volume in this patient group [20]. Postmortem analyses of human cortex tissue from this study revealed a decreased NPTX2 level, supporting the biological significance of the CSF findings [20].
Our small-scale analyses of spinal cord sections in the established late-onset SMA mouse model [22] revealed an upregulation of NPTX2 in the ventral horn of the spinal cord on day 42, when motor neuron loss had already begun, aligning with our human CSF and serum findings. While these analyses do not aim to elucidate pathomechanisms, they underscore the relevance of NPTX2 concentrations in human CSF and serum as potential biomarkers.
Our findings suggest that NPTX2 plays a role in the pathomechanism of adult SMA. As a biomarker in CSF, it may provide biological-mechanistic insights into the therapeutic response to SMN-restoring therapy with nusinersen. Measured in serum, it could be useful as a biomarker to objectify the clinical-functional treatment response, though our results do not support that it allows for its prediction.
Multi-omics analyses, performed by mass spectrometric analyses, have become an essential tool in understanding underlying molecular mechanisms and identifying biomarkers especially of rare diseases like SMA [37,38]. Proteomic studies, for example, have demonstrated that CSF proteins involved in axonogenesis, as well as bioenergetic and inflammatory pathways, exhibit differential expression in pwSMA compared to controls, and also before and after nusinersen treatment [39,40]. Notably, these effects were evident even in small sample sizes, with as few as 10 to 13 pwSMA, and [38,41]. However, the technique faces challenges, including difficulties in detecting low-abundance proteins, the complexity of data analysis, and technical limitations such as sample degradation and incomplete coverage of the proteome. Furthermore, techniques are not standardized, making the results of different studies often poorly comparable.
CNTN5 is a protein of the contactin family, playing a role in cell adhesion and neuronal development. It is primarily expressed in the central nervous system and is involved in synapse formation and function. It has been studied as a biomarker in autoimmune nodopathies [42]. ANTXR1 is a membrane protein which is primarily expressed on the surface of cells, particularly endothelial cells and macrophages. Anthrax toxins receptors were shown to be expressed in the nervous system, especially in sensory neurons [43]. While significant differences in CNTN5 and ANTXR1 levels between improvers and non-improvers were observed in the proteomic analysis, we could not verify this in the quantitative measurement using ELISA in serum or CSF. The lack of concordance between the results of the proteomic analysis and the quantification via ELISA may be due to the different sensitivity and specificity of the tests, a different type of sample preparation and processing for the proteome analysis, different calibration strategies and reference standards and lower resolving power with e.g. interferences in the proteome analysis. In view of the increasing use of new explorative methods such as proteome analysis, our data also demonstrate that the results should be interpreted with caution and should be confirmed by quantitative analyses.
Limitations of our study included the sample size with some missing data for later time points. The sample selection, characterized by very heterogeneous treatment responses as indicated by changes in HFMSE scores, may potentially confound the correlational statistics with the biomarkers and is not optimally representative of a real-world SMA cohort. Limitations of proteomic analysis primarily include the complexity and dynamics of the proteome, with proteins existing in various isoforms, complicating identification and quantification, as well as the general detection threshold and quantification accuracy. The limitations of the ELISA assay include potential cross-reactivity and limited sensitivity.

5. Conclusions

In this study, we were able to identify NPTX2 through a simple ELISA analysis in a small cohort as a promising biomarker to objectify the treatment response of adults with SMA to nusinersen. These findings should be confirmed in a larger cohort to validate the results within a more robust large-scale framework. Furthermore, our results imply a significant pathophysiological role of NPTX2 in late-onset SMA, demanding further investigation. It would be particularly interesting to closer investigate the interactions between NPTX2, NPTX1 and NPTXR, considering their formation of complexes and interactions in their expression patterns [30].

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org, Table S1: Protein ratios between groups.

Author Contributions

"Conceptualization, S.N., L.-I.S., M.L. and T.H.; Data collection and analysis, S.N., L.-I.S., N.I.T., S.H. and K.C.L.; Writing – Original Draft Preparation, S.N.; Writing – Review & Editing, L.-I.S., K.C.L., N.I.T., S.H., C.K., M.L. and T.H.; Funding Acquisition, T.H.”. All authors read and approved the final manuscript.

Funding

This work was financially supported by the Myelin Repair Foundation. The authors also want to thank the West German Biobank Essen (WBE) for sample processing and storage.

Institutional Review Board Statement

This study was performed in line with the principles of the Declaration of Helsinki. It was approved by the local ethics committee of the University of Duisburg Essen, Germany (approval number 18-8071-BO). The experiments with animals were conducted under the animal welfare guidelines of the University Duisburg Essen. The use of the SMA mouse model was approved by the State Agency for Nature, Environment and Consumer Protection (LANUV) in North Rhine-Westphalia (reference number 81–02.04.2020.A335).

Informed Consent Statement

Informed consent was obtained from all individual participants included in the study.

Data Availability Statement

The protein ratios between the groups obtained from the proteomic analysis are available as a supplementary online resource (Table S1). The mass spectrometry proteomics data have been deposited to the Proteo-meXchange Consortium via the PRIDE partner repository with the dataset identifier PXD065345. The other datasets generated and analysed during the current study are not publicly available due to the inclusion of clinical patient data that cannot be disclosed for data protection and privacy reasons but are available from the corresponding author on reasonable request.

Acknowledgments

We thank Biognosys AG for performing the proteomic analysis and providing all associated data.

Conflicts of Interest

Financial interests: Authors L.-I.S., K.C.L, S.H., N.I.T. and M.L. declare that they have no financial interests. Author S.N. has received speaker honoraria from Roche and Biogen and travel support from Johnson & Johnson. C.K. received lecture and consultancy fees from Biogen, Roche, and Novartis. T.H. received lecture and consultancy fees from Biogen, Roche and Novartis, as well as research support from Biogen, Roche, and Novartis Gene Therapies.

References

  1. Lefebvre S, Bürglen L, Reboullet S, et al (1995) Identification and characterization of a spinal muscular atrophy-determining gene. Cell 80:155–165. [CrossRef]
  2. Lefebvre S, Burlet P, Liu Q, et al (1997) Correlation between severity and SMN protein level in spinal muscular atrophy. Nat Genet 16:265–269. [CrossRef]
  3. Kolb SJ, Battle DJ, Dreyfuss G (2007) Molecular functions of the SMN complex. J Child Neurol 22:990–994. [CrossRef]
  4. Calucho M, Bernal S, Alías L, et al (2018) Correlation between SMA type and SMN2 copy number revisited: An analysis of 625 unrelated Spanish patients and a compilation of 2834 reported cases. Neuromuscul Disord NMD 28:208–215. [CrossRef]
  5. Wadman RI, Wijngaarde CA, Stam M, et al (2018) Muscle strength and motor function throughout life in a cross-sectional cohort of 180 patients with spinal muscular atrophy types 1c-4. Eur J Neurol 25:512–518. [CrossRef]
  6. Hagenacker T, Wurster CD, Günther R, et al (2020) Nusinersen in adults with 5q spinal muscular atrophy: a non-interventional, multicentre, observational cohort study. Lancet Neurol 19:317–325. [CrossRef]
  7. Günther R, Wurster CD, Brakemeier S, et al (2024) Long-term efficacy and safety of nusinersen in adults with 5q spinal muscular atrophy: a prospective European multinational observational study. Lancet Reg Health Eur 39:100862. [CrossRef]
  8. Mercuri E, Deconinck N, Mazzone ES, et al (2022) Safety and efficacy of once-daily risdiplam in type 2 and non-ambulant type 3 spinal muscular atrophy (SUNFISH part 2): a phase 3, double-blind, randomised, placebo-controlled trial. Lancet Neurol 21:42–52. [CrossRef]
  9. Brakemeier S, Lipka J, Schlag M, et al (2024) Risdiplam improves subjective swallowing quality in non-ambulatory adult patients with 5q-spinal muscular atrophy despite advanced motor impairment. J Neurol. [CrossRef]
  10. Maggi L, Bello L, Bonanno S, et al (2020) Nusinersen safety and effects on motor function in adult spinal muscular atrophy type 2 and 3. J Neurol Neurosurg Psychiatry 91:1166–1174. [CrossRef]
  11. Coratti G, Pera MC, Montes J, et al (2021) Different trajectories in upper limb and gross motor function in spinal muscular atrophy. Muscle Nerve 64:552–559. [CrossRef]
  12. Gaiani A, Martinelli I, Bello L, et al (2017) Diagnostic and Prognostic Biomarkers in Amyotrophic Lateral Sclerosis: Neurofilament Light Chain Levels in Definite Subtypes of Disease. JAMA Neurol 74:525–532. [CrossRef]
  13. Martin S-J, McGlasson S, Hunt D, Overell J (2019) Cerebrospinal fluid neurofilament light chain in multiple sclerosis and its subtypes: a meta-analysis of case-control studies. J Neurol Neurosurg Psychiatry 90:1059–1067. [CrossRef]
  14. Bayoumy S, Verberk IMW, Vermunt L, et al (2024) Neurofilament light protein as a biomarker for spinal muscular atrophy: a review and reference ranges. Clin Chem Lab Med 62:1252–1265. [CrossRef]
  15. Cordts I, Fuetterer C, Wachinger A, et al (2025) Long-Term Dynamics of CSF and Serum Neurofilament Light Chain in Adult Patients With 5q Spinal Muscular Atrophy Treated With Nusinersen. Neurology 104:e213371. [CrossRef]
  16. Totzeck A, Stolte B, Kizina K, et al (2019) Neurofilament Heavy Chain and Tau Protein Are Not Elevated in Cerebrospinal Fluid of Adult Patients with Spinal Muscular Atrophy during Loading with Nusinersen. Int J Mol Sci 20:5397. [CrossRef]
  17. Birhanu AG (2023) Mass spectrometry-based proteomics as an emerging tool in clinical laboratories. Clin Proteomics 20:32. [CrossRef]
  18. Mc P, G C, N F, et al (2017) Content validity and clinical meaningfulness of the HFMSE in spinal muscular atrophy. BMC Neurol 17:. [CrossRef]
  19. Mazzone ES, Mayhew A, Montes J, et al (2017) Revised upper limb module for spinal muscular atrophy: Development of a new module. Muscle Nerve 55:869–874. [CrossRef]
  20. Xiao M-F, Xu D, Craig MT, et al NPTX2 and cognitive dysfunction in Alzheimer’s Disease. eLife 6:e23798. [CrossRef]
  21. Oh S, Jang Y, Na CH (2023) Discovery of Biomarkers for Amyotrophic Lateral Sclerosis from Human Cerebrospinal Fluid Using Mass-Spectrometry-Based Proteomics. Biomedicines 11:1250. [CrossRef]
  22. Schmitt LI, David C, Steffen R et al. (2023) Spinal astrocyte dysfunction drives motor neuron loss in late-onset spinal muscular atrophy. Acta Neuropathol 145:611–635. [CrossRef]
  23. Chapman G, Shanmugalingam U, Smith PD (2020) The Role of Neuronal Pentraxin 2 (NP2) in Regulating Glutamatergic Signaling and Neuropathology. Front Cell Neurosci 13:575. [CrossRef]
  24. Bowerman M, Swoboda KJ, Michalski J-P, et al (2012) Glucose Metabolism and Pancreatic Defects in Spinal Muscular Atrophy. Ann Neurol 72:256–268. [CrossRef]
  25. Reilly A, Beauvais A, Al-Aarg M, et al (2024) Peripheral defects precede neuromuscular pathology in the Smn2B/- mouse model of spinal muscular atrophy. J Neuromuscul Dis 11:1200–1210. [CrossRef]
  26. Tsui CC, Copeland NG, Gilbert DJ, et al (1996) Narp, a novel member of the pentraxin family, promotes neurite outgrowth and is dynamically regulated by neuronal activity. J Neurosci Off J Soc Neurosci 16:2463–2478. [CrossRef]
  27. Chang MC, Park JM, Pelkey KA, et al (2010) Narp regulates homeostatic scaling of excitatory synapses on Parvalbumin interneurons. Nat Neurosci 13:1090–1097. [CrossRef]
  28. O’Brien RJ, Xu D, Petralia RS, et al (1999) Synaptic clustering of AMPA receptors by the extracellular immediate-early gene product Narp. Neuron 23:309–323. [CrossRef]
  29. O’Brien R, Xu D, Mi R, et al (2002) Synaptically targeted narp plays an essential role in the aggregation of AMPA receptors at excitatory synapses in cultured spinal neurons. J Neurosci Off J Soc Neurosci 22:4487–4498. [CrossRef]
  30. Lee S-J, Wei M, Zhang C, et al (2017) Presynaptic Neuronal Pentraxin Receptor Organizes Excitatory and Inhibitory Synapses. J Neurosci Off J Soc Neurosci 37:1062–1080. [CrossRef]
  31. Cho RW, Park JM, Wolff SBE, et al (2008) mGluR1/5-dependent long-term depression requires the regulated ectodomain cleavage of neuronal pentraxin NPR by TACE. Neuron 57:858–871. [CrossRef]
  32. Schwarz DA, Barry G, Mackay KB, et al (2002) Identification of differentially expressed genes induced by transient ischemic stroke. Brain Res Mol Brain Res 101:12–22. [CrossRef]
  33. Leo M, Schmitt L-I, Fleischer M, et al (2022) Induction of Survival of Motor Neuron (SMN) Protein Deficiency in Spinal Astrocytes by Small Interfering RNA as an In Vitro Model of Spinal Muscular Atrophy. Cells 11:558. [CrossRef]
  34. Welby E, Ebert AD (2023) Diminished motor neuron activity driven by abnormal astrocytic EAAT1 glutamate transporter activity in spinal muscular atrophy is not fully restored after lentiviral SMN delivery. Glia 71:1311–1332. [CrossRef]
  35. Göverti D, Büyüklüoğlu N, Kaya H, et al (2022) Neuronal pentraxin-2 (NPTX2) serum levels during an acute psychotic episode in patients with schizophrenia. Psychopharmacology (Berl) 239:2585–2591. [CrossRef]
  36. Shao K, Shan S, Ru W, Ma C (2020) Association between serum NPTX2 and cognitive function in patients with vascular dementia. Brain Behav 10:e01779. [CrossRef]
  37. Lunke S, Bouffler SE, Patel CV, et al (2023) Integrated multi-omics for rapid rare disease diagnosis on a national scale. Nat Med 29:1681–1691. [CrossRef]
  38. Kessler T, Latzer P, Schmid D, et al (2020) Cerebrospinal fluid proteomic profiling in nusinersen-treated patients with spinal muscular atrophy. J Neurochem 153:650–661. [CrossRef]
  39. Panicucci C, Sahin E, Bartolucci M, et al (2024) Proteomics profiling and machine learning in nusinersen-treated patients with spinal muscular atrophy. Cell Mol Life Sci 81:393. [CrossRef]
  40. Meneri M, Abati E, Gagliardi D, et al (2023) Identification of Novel Biomarkers of Spinal Muscular Atrophy and Therapeutic Response by Proteomic and Metabolomic Profiling of Human Biological Fluid Samples. Biomedicines 11:1254. [CrossRef]
  41. Zandl-Lang M, Züllig T, Holzer M, et al (2025) Multi-omics profiling in spinal muscular atrophy (SMA): investigating lipid and metabolic alterations through longitudinal CSF analysis of Nusinersen-treated patients. J Neurol 272:183. [CrossRef]
  42. Wieske L, Martín-Aguilar L, Fehmi J, et al (2021) Serum Contactin-1 in CIDP: A Cross-Sectional Study. Neurol Neuroimmunol Neuroinflammation 8:e1040. [CrossRef]
  43. Yang NJ, Isensee J, Neel DV, et al (2022) Anthrax toxins regulate pain signaling and can deliver molecular cargoes into ANTXR2+ DRG sensory neurons. Nat Neurosci 25:168–179. [CrossRef]
Figure 1. Overview of the study design. Control samples were only available from a single time point. pwSMA = adult individuals with 5q-SMA; SMA Improvers = increase ≥ 3 points in HMFSE after 6 months of treatment with nusinersen; SMA Non-Improvers = increase of < 3 points in HFMSE after 6 months of treatment with nusinersen; HFMSE = Hammersmith Functional Motor Scale Expanded; RULM = Revised Upper Limb Module; mo. = months; wt = wild type; P20, P40, P72 and P>400 refer to postnatal days, with P20 indicating postnatal day 20.
Figure 1. Overview of the study design. Control samples were only available from a single time point. pwSMA = adult individuals with 5q-SMA; SMA Improvers = increase ≥ 3 points in HMFSE after 6 months of treatment with nusinersen; SMA Non-Improvers = increase of < 3 points in HFMSE after 6 months of treatment with nusinersen; HFMSE = Hammersmith Functional Motor Scale Expanded; RULM = Revised Upper Limb Module; mo. = months; wt = wild type; P20, P40, P72 and P>400 refer to postnatal days, with P20 indicating postnatal day 20.
Preprints 165388 g001
Figure 2. Results of the proteomic analysis with a) Heatmap of differentially expressed proteins between pwSMA and controls (Ctrl). The color scale on the right represents Z-score values, with yellow shades indicating increased expression levels and purple shades indicating decreased expression levels. SMA BS (Baseline) I to V represent SMA Improvers (HFMSE improvement ≥3 points after 6 months of nusinersen therapy), while SMA BS VI and VII represent SMA Non-Improvers (HFMSE improvement <3 points). Each row corresponds to a protein and each column represents a sample. The visual cluster structure of the heatmap reveals differences in protein expression between individuals with SMA and controls; b) Volcano Plot of differentially expressed proteins between pwSMA and controls (ctrl). The logarithmic change in expression levels (log2 fold change) is represented on the X-axis, while the statistical significance is shown as the negative decadic logarithm of the p-value (-log10(p-value)) on the Y-axis. Each point represents a different protein, with color-coded points indicating statistically significant differences (p < 0.05) between pwSMA and controls. Red points indicate a positive change in protein expression compared to controls, while blue points represent a negative change; c) Results of proteomic profiling from CSF with baseline levels of CNTN5, ANTXR1 and NPTX2 in controls and SMA Improvers and Non-Improvers. Significant differences (p <0.05) are marked with an asterix (*). CNTN5 levels are decreased in SMA Improvers compared to Non-Improvers, with no significant difference compared to controls. ANTXR1 is elevated in SMA Improvers compared to Non-Improvers, again without a significant difference compared to controls. For NPTX2, no significant group differences were observed.
Figure 2. Results of the proteomic analysis with a) Heatmap of differentially expressed proteins between pwSMA and controls (Ctrl). The color scale on the right represents Z-score values, with yellow shades indicating increased expression levels and purple shades indicating decreased expression levels. SMA BS (Baseline) I to V represent SMA Improvers (HFMSE improvement ≥3 points after 6 months of nusinersen therapy), while SMA BS VI and VII represent SMA Non-Improvers (HFMSE improvement <3 points). Each row corresponds to a protein and each column represents a sample. The visual cluster structure of the heatmap reveals differences in protein expression between individuals with SMA and controls; b) Volcano Plot of differentially expressed proteins between pwSMA and controls (ctrl). The logarithmic change in expression levels (log2 fold change) is represented on the X-axis, while the statistical significance is shown as the negative decadic logarithm of the p-value (-log10(p-value)) on the Y-axis. Each point represents a different protein, with color-coded points indicating statistically significant differences (p < 0.05) between pwSMA and controls. Red points indicate a positive change in protein expression compared to controls, while blue points represent a negative change; c) Results of proteomic profiling from CSF with baseline levels of CNTN5, ANTXR1 and NPTX2 in controls and SMA Improvers and Non-Improvers. Significant differences (p <0.05) are marked with an asterix (*). CNTN5 levels are decreased in SMA Improvers compared to Non-Improvers, with no significant difference compared to controls. ANTXR1 is elevated in SMA Improvers compared to Non-Improvers, again without a significant difference compared to controls. For NPTX2, no significant group differences were observed.
Preprints 165388 g002
Figure 3. a) CSF and c) Serum levels of NPTX2 in controls and pwSMA at T0 (baseline), T1 (2 months after start of treatment nusinersen) and T2 (14 months after start of treatment with nusinersen) across all pwSMA shown in a bar plot with individual data points and error bars representing the standard deviation. The bar height represents the mean value. Across the entire pwSMA cohort, CSF NPTX2 levels were lower at T2 compared to T0, and lower at T2 compared to T1. No significant differences were observed in serum levels across time points. b) CSF and d) Serum levels of NPTX2 at T0, T1 and T2 in SMA Improvers (improvement in HFMSE ≥3 points after 6 months of nusinersen therapy) vs. SMA Non-improvers (improvement in HFMSE <3 points after 6 months of nusinersen therapy). Significant differences (p<0.05) are marked by an asterisk *. Note that the Y-axis scaling differs between serum and CSF values due to a different range of protein concentration. At T2, serum NPTX2 levels were lower in SMA Improvers compared to Non-Improvers. No significant differences in CSF NPTX2 levels were observed between SMA Improvers and Non-Improvers at the respective time points.
Figure 3. a) CSF and c) Serum levels of NPTX2 in controls and pwSMA at T0 (baseline), T1 (2 months after start of treatment nusinersen) and T2 (14 months after start of treatment with nusinersen) across all pwSMA shown in a bar plot with individual data points and error bars representing the standard deviation. The bar height represents the mean value. Across the entire pwSMA cohort, CSF NPTX2 levels were lower at T2 compared to T0, and lower at T2 compared to T1. No significant differences were observed in serum levels across time points. b) CSF and d) Serum levels of NPTX2 at T0, T1 and T2 in SMA Improvers (improvement in HFMSE ≥3 points after 6 months of nusinersen therapy) vs. SMA Non-improvers (improvement in HFMSE <3 points after 6 months of nusinersen therapy). Significant differences (p<0.05) are marked by an asterisk *. Note that the Y-axis scaling differs between serum and CSF values due to a different range of protein concentration. At T2, serum NPTX2 levels were lower in SMA Improvers compared to Non-Improvers. No significant differences in CSF NPTX2 levels were observed between SMA Improvers and Non-Improvers at the respective time points.
Preprints 165388 g003
Figure 4. Correlation of the change in NPTX2 serum concentration and HFMSE score from T0 (baseline) to T2 (14 months after treatment initiation with nusinersen) (r=-0.688, p=0.019, n=11). The points represent individual pwSMA. The line of best fit was calculated using simple linear regression, R² = 0.52, p = 0.01.
Figure 4. Correlation of the change in NPTX2 serum concentration and HFMSE score from T0 (baseline) to T2 (14 months after treatment initiation with nusinersen) (r=-0.688, p=0.019, n=11). The points represent individual pwSMA. The line of best fit was calculated using simple linear regression, R² = 0.52, p = 0.01.
Preprints 165388 g004
Figure 5. Immunostaining of a) spinal motoneurons (SMI-32, red) and b) NPTX2 (green) in the ventral horn of lumbar spinal cord slices from wild type (wt) and SMA mice (SMA) at P20 (postnatal day 20), P42 and P70 and P>400. The SMA values normalized to the wt values for the immunoreactivity of NPTX2 and the absolute number of motoneurons at the different time points are shown on the right (n=9 for P20 and P>400; n=7 for P42 and P70; 3 SMA mice per condition with 2 or 3 slices per animal). Scale bar: 20 µm. From P42 onwards, the number of motor neurons is reduced in SMA mice compared to wt mice. At P42 - but not at any other time point - NPTX2 immunofluorescence is increased in SMA compared to wt mice. c) Dynamics of wt-normalized values for motoneuron count and NPTX2 immunoreactivity in SMA mice across the different time points. In SMA mice, the number of motor neurons continues to decline progressively from P42 onwards, while NPTX2 immunofluorescence peaks at P42.
Figure 5. Immunostaining of a) spinal motoneurons (SMI-32, red) and b) NPTX2 (green) in the ventral horn of lumbar spinal cord slices from wild type (wt) and SMA mice (SMA) at P20 (postnatal day 20), P42 and P70 and P>400. The SMA values normalized to the wt values for the immunoreactivity of NPTX2 and the absolute number of motoneurons at the different time points are shown on the right (n=9 for P20 and P>400; n=7 for P42 and P70; 3 SMA mice per condition with 2 or 3 slices per animal). Scale bar: 20 µm. From P42 onwards, the number of motor neurons is reduced in SMA mice compared to wt mice. At P42 - but not at any other time point - NPTX2 immunofluorescence is increased in SMA compared to wt mice. c) Dynamics of wt-normalized values for motoneuron count and NPTX2 immunoreactivity in SMA mice across the different time points. In SMA mice, the number of motor neurons continues to decline progressively from P42 onwards, while NPTX2 immunofluorescence peaks at P42.
Preprints 165388 g005
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

Disclaimer

Terms of Use

Privacy Policy

Privacy Settings

© 2025 MDPI (Basel, Switzerland) unless otherwise stated