Prerpints.org logo
Preprint
Article

Evaluation of the Effect of β-Wrapin AS69 in a Mouse Model Based on Alpha-Synuclein Overexpression

Altmetrics

Downloads

111

Views

59

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

02 May 2024

Posted:

03 May 2024

You are already at the latest version

Alerts
Abstract
Aggregation of the protein α-Synuclein (αSyn) is a hallmark of Parkinson's disease (PD), dementia with Lewy bodies (DLB) and multiple systems atrophy, and alleviating the extent of αSyn pathology is an attractive strategy against neurodegeneration. The engineered binding protein β-wrapin AS69 binds monomeric αSyn. AS69 reduces primary and secondary nucleation as well as fibril elongation in vitro. It also mitigates aSyn pathology in a mouse model based on intrastriatal injection of aSyn pre-formed fibrils (PFF). Since the PFF-based model does not represent all aspects of PD, we tested here whether AS69 can reduce neurodegeneration resulting from αSyn overexpression. Human A53T-αSyn was overexpressed in the mouse Substantia nigra (SN) by using recombinant adeno-associated viral vector (rAAV). AS69 was also expressed by rAAV transduction. Behavioral tests and immunofluorescence stainings were used as outcomes.Transduction with rAAV-αSyn resulted in αSyn pathology as reported by phospho-αSyn staining and caused degeneration of dopaminergic neurons in the SN. The co-expression of rAAV-AS69 did not reduce αSyn pathology or degeneration of dopaminergic neurons. We conclude that αSyn monomer binding by rAAV-AS69 was insufficient to protect from aSyn pathology resulting from αSyn overexpression.
Keywords: 
Subject: Medicine and Pharmacology  -   Neuroscience and Neurology

1. Introduction

Parkinson’s disease (PD), multiple systems atrophy and dementia with Lewy bodies (LBD) are characterized by the aggregation of the small synaptic protein α-Synuclein (αSyn). The αSyn proteinopathy spreads between neurons and glial cells as well as interconnected brain regions. Alleviating the extent of αSyn pathology, e.g., by preventing the formation of synuclein aggregates, is therefore considered an attractive strategy against neurodegeneration in PD and related disorders [1,2].
AS69 is an engineered dimeric protein that induces a local folding of αSyn residues 37-54 into a β-hairpin [3]. It is thereby able to prevent primary and secondary nucleation as well as fibril elongation by sequestering monomeric αSyn [3,4]. AS69′s ability to inhibit the formation of αSyn aggregates has been demonstrated in HEK293T-cells and a Drosophila melanogaster model based on neuronal αSyn overexpression [4]. Furthermore, AS69 alleviates αSyn pathology in primary mouse neurons and in a mouse model based on intrastriatal injection of preformed αSyn fibrils [5]. Small molecules that inhibit the aggregation of αSyn, such as UCB0599 and Anle138b, are being studied in clinical trials for Parkinson’s disease (NCT04685265, NCT04875962).
Injecting αSyn fibrils into mouse brain can induce αSyn pathology and dopaminergic neuron death [6]. Yet, fibril-based models only reproduce a part of αSyn pathology [7,8]. In patients, PD can be triggered by increased expression of αSyn - resulting from either triplication of the αSyn locus or from genetic variants in the αSyn promoter region [9,10,11]. Overexpression of aSyn therefore constitutes another realistic model of PD pathogenesis. In order to determine the utility of monomer stabilization, we tested AS69 - expressed by transduction with a recombinant adeno-associated vector (rAAV) - in a mouse model that is based on rAAV-mediated overexpression of human A53T-αSyn (H-αSyn) in the Substantia nigra (SN).

2. Materials and Methods

2.1. Animals

All animal experiments were carried out in accordance with the European Communities Council Directive of November 24, 1986 (86/609/EEC) and approved by the Landesdirektion Dresden, Germany (25-5131/49615). Twenty 12-week-old C57BLl/6 male mice from Charles River Laboratories (Wilmington, Massachusetts, United States) were housed under a 12 h light and dark cycle with free access to pellet food and water in the Experimental Center, TU Dresden, Dresden, Germany.

2.2. Recombinant rAAV Production and Purification

The production of rAAVs was performed at the Leuven Viral Vector Core as described previously [12]. Pseudotype rAAV2/7 was used due to its efficient uptake in neurons. The ssDNA virus encodes for the following sequences: ITR - CMVie-enhanced synapsin1 promoter - transgene - WPRE - bovine growth hormone polyadenylation sequence - ITR. Either green fluorescent protein (GFP), human A53T-αSyn or AS69 were inserted as a transgene. As technical titers, genomic copies (GC) were determined by real-time PCR analysis. All viral vectors were diluted to a concentration of 7*1011 GC/mL. Viral vector solutions were aliquoted and freshly thawed before stereotactic injection. The GFP control cohort received 2 µl of rAAV-GFP, while the other groups received a 2 µl of a 1:1 mix of the respective viral vectors.

2.3. Stereotactic Surgery & Tissue Processing

Stereotactic surgeries and injection of viral vector solutions were performed unilaterally into the SN. The following stereotactic coordinates relative to Bregma were used: Anterior-posterior -3 mm, lateral -1.2 mm, dorso-ventral -4.1 mm. All procedures were performed using aseptic techniques. Mice were anesthetized using Ketamine (100 mg/kg i.p) and Xylazine (10 mg/kg i.p.). Animals were placed in a flat skull position in a stereotactic head frame. A bore hole craniotomy was performed and was followed by a dural incision. A 30 G needle and a 10 µl Hamilton syringe (Hamilton Bonaduz AG, Bonaduz, Switzerland) were used for the injection of 2 µl vector solution. The needle was advanced at a speed of 0.3 mm/min and was left in place at the final coordinates for 5 min before and after the injection of the vector. The vector dilution was injected at a rate of 200 nl/min. After a slow withdrawal, the skin was sutured and lidocaine (Aspen Pharmacare, Durban, South Africa) was topically applied. Eight weeks later, mice were sacrificed with an overdose of isoflurane (Baxter, Lessines, Belgium) and transcardially perfused with 4 % paraformaldehyde (PFA) diluted in Tris-buffered saline (TBS, pH 7.6). The brains were kept in PFA 4 % for another 48 h at 4 °C and were afterwards transferred to 30 % sucrose in TBS for cryoprotection. The brains were frozen at -55 °C in isopentane and stored at -80 °C until 30 µm thick coronal brain sections were obtained using a Cryostat (Leica Biosystems, Nussloch, Germany).

2.4. Cylinder Test

The cylinder test was performed at three time points: one day before stereotactic injection, four weeks and eight weeks after surgery. One animal at a time was transferred into a 15-cm-wide glass cylinder and up to 25 contacts with the cylinder wall were videotaped. The time until the dominant forepaw performed 25 contacts was scored based on the videos.

2.5. Immunofluorescence Stainings

The brain sections were stained for the following antigens: tyrosine hydroxylase (TH), human αSyn (H-αSyn), phosphorylated αSyn at Serine 129 (P-αSyn), glial fibrillary acidic protein (GFAP) and ionized calcium-binding adapter molecule 1 (Iba1). First, the sections were rinsed three times for 10 minutes in TBS and afterwards incubated in 10 % donkey serum (BIOZOL Diagnostica Vertrieb GmbH, Eching, Germany), 0.2 % Triton X-100 (Thermo Scientific, Waltham, United States) diluted in TBS for 1 h at room temperature. Next, the sections were incubated for 18 h in sheep anti-TH (1:2000, P60101, Pel-Freez), rabbit anti-P-αSyn (1:1000, ab51253, Abcam), rat anti-H-αSyn (1:1000, ALX-804-258-L001, ENZO Life Sciences), chicken anti-GFAP (1:500, ab4674, Abcam) or guinea pig anti-Iba1 (1:2000, HS-234308, Histo Sure) at 4°C. After three additional 10 min rinses, fluorophore-conjugated secondary antibodies were added for 1 h at 21°C: Alexa 647-conjugated donkey anti-sheep (1:1000, A21448, Molecular Probes), donkey anti-rat 647 (1:1000, ab150155, Abcam), Alexa 647-conjugated donkey anti-chicken (1:500, 703-605-155, Jackson ImmunoResearch), Alexa 555-conjugated donkey anti-rabbit (1:1000, A31572, Molecular Probes) and Alexa 555-conjugated donkey anti-guinea pig (1:1000, A-21435, Molecular Probes). Nuclei were counterstained with Hoechst (1:2000 for 5 min; Invitrogen, Waltham, United States) and the sections were mounted in Fluoromount-G (Invitrogen).

2.6. Quantification of TH-Positive Neurons, Analysis of TH-Positive Dendrites and Striatal TH-Positive Axon Terminals

For the quantification of dopaminergic neurons in the SN, every fourth section (about 10 sections per animal) of the SN was stained for TH. Slide scans were acquired by the Light Microscopy Facility (DZNE Bonn, Germany) using a slide scanning microscope - AxioScan.Z1 (Zeiss, Jena, Germany) equipped with a 20x/0.8NA objective. Five Z-levels were imaged at intervals of one micrometer. TH-positive neurons in the Substantia nigra pars compacta (SNc) were counted manually and investigator-blinded in each hemisphere using Zeiss Zen 3.1 software (Zeiss). To obtain an estimate of the total number of TH-positive neurons per hemisphere, the number of TH-positive SNc neurons of all slices was summed; the sum of neurons in the injected hemispheres was expressed relative to the uninjected hemisphere.
Dopaminergic neurons extend long apical dendrites into the Substantia nigra pars reticulata (SNr). The density of these dendrites was assessed by measuring the average intensity of TH staining within each hemisphere’s SNr in every slide scan mentioned above. The signal intensity was averaged across each hemisphere for each animal. Finally, mean intensity in the injected hemispheres was expressed relative to the uninjected hemisphere.
Dopaminergic neurons project their axons into the striatum. In order to analyze their integrity, we collected ten to fifteen high magnification images from two to three sections from the dorsal striatum using a spinning disk confocal microscope. The set up uses a Zeiss Axio Observer.Z1 Inverted Microscope (Zeiss) and a Yokogawa CSU-X1 unit (Yokogawa Life Science, Musashino-shi, Tokyo) running a 40x/0.95 objective. Image analysis was performed in CellProfiler using the unbuilt “Enhance Neurites” function followed by thresholding the image using the “robust background” algorithm [13]. Next we calculated for each animal the relative mean count of positive pixels per hemisphere and presented the values from the injected hemisphere normalized to the contralateral hemisphere.

2.7. Evaluation of αSyn Pathology

The expression of H-αSyn and P-αSyn within the SNc was analyzed by the following protocol: Scans from every fourth section were acquired using a spinning disk confocal microscope. The set up uses a Zeiss Axio Observer.Z1 Inverted Microscope (Zeiss) and a Yokogawa CSU-X1 unit (Yokogawa Life Science, Musashino-shi, Tokyo) running a 20x/0.8NA objective. Five Z-levels at intervals of one µm were scanned per field of view. Maximum intensity projections were created. In order to quantify the number of H-αSyn and P-αSyn positive cells we worked with an open source software, namely QuPath [14]. First, the SNc within the treated hemisphere was carefully encircled. Next, the cell detection function was applied using its default parameters, which detects nuclei based on their Hoechst staining. We manually annotated examples for H-αSyn positive cells and P-αSyn positive cells as well as a background signal in every SNc which allowed us to train the in-build object classifier algorithm which analyzed the whole SNc based on our manually annotated examples. The count of positive cells (H-αSyn, P-αSyn) per animal was summed and the sum was multiplied by four since every fourth section was analyzed.

2.8. Analysis of Astrogliosis & Microgliosis

For quantification of astrogliosis and microgliosis, fifteen to twenty images within the SNc of both hemispheres’ per animal spanning four sections were acquired using the spinning disc confocal microscope described above featured by a 40x/0.95NA objective. Seven z-levels at a step size of 0.48 µm were covered per image. Images were processed as maximum intensity projections. The GFAP channel was thresholded via “robust background method” in CellProfiler and the resulting binary images were opened once to reduce salt-and-pepper noise. The count of positive pixels was divided by the total count of pixels per image. Furthermore, the binary image was used to create a masked image of the original GFAP channel and the signal intensity in the GFAP positive area was determined. We calculated the mean GFAP positive area and the mean GFAP intensity per hemisphere and per animal. The results of the injected hemisphere were expressed relative to the uninjected hemisphere.
In order to evaluate the degree of microglia activation, we analyzed the Iba1 channel from the same images as described for GFAP (15-20 images per animal from four sections and both hemispheres). The analysis was performed manually by classifying each Iba1 positive cell into one of the three following categories based on its morphology: I) ramified, II) hypertrophic or dystrophic or III) amoeboid or rod-shaped microglia [15]. Next, we determined the mean count of microglia per image for each animal and hemisphere. The average number of microglia per image in the injected hemisphere is shown relative to the average per image number of microglia in the uninjected hemisphere. We also show, based on data from the injected hemisphere, the mean count of each microglia subtype relative to the mean count of all microglia categories.

2.9. Statistical Analysis

Data analysis and visualization were performed in R (version 4.2.0) using RStudio (“Spotted Wakerobin” Release, 7872775e, 2022-07-22) and the following packages: Tidyverse for data wrangling and visualization; Ggpubr, FSA and Corrplot for statistical testing [16,17,18,19,20,21]. The Shapiro-Wilk normality test was used to determine normal distribution of the data and the Bartlett’s test to assess equal population variances. For normally distributed data we applied Welch’s two sample t-Test, Welch’s two sample t-Tests with Benjamini-Hochberg p-value correction for multiple testing or a two-way ANOVA followed by Tukey’s HSD depending on the number of comparisons. Data that did not meet above described standards was compared using the Wilcoxon rank sum exact test for pairwise comparisons or Kruskal-Wallis test with Holm’s p-value correction for multiple testing followed by Dunn-test as post-hoc-test. In order to assess correlations we calculated Pearson’s product-moment correlation. Error bars represent mean ± SEM. Exact p-values are listed in the figure legends.

3. Results

3.1. Overexpression of αSyn in Mouse SN Induced a Subtle Motor Deficit

In this study, we induced αSyn pathology and neurodegeneration in mice by transducing SN neurons with human A53T-αSyn using rAAV vectors [12,22]. Transduction with rAAV-GFP was used as negative control. Cylinder tests were performed at three time points: one day before stereotactic injections and 4 weeks as well as 8 weeks post-injection (Figure 1A). We measured the time it took each animal to touch the cylinder wall 25 times (Figure 1B). At baseline, it took all groups equally long. At 4 weeks, it took animals transduced with rAAV-αSyn+rAAV-GFP (blue in Figure 1B) somewhat longer to complete 25 contacts than observed in the other groups, but the difference was not statistically significant. At eight weeks, mice transduced with rAAV-αSyn+ rAAV-GFP required on average almost twice as long to complete 25 contacts than controls did. We conclude that αSyn overexpression induced deficits in animals’ motor abilities.
Furthermore we studied the consequences of transducing SN neurons with AS69, also using rAAV vectors. Again, transduction with rAAV-GFP was used as negative control, resulting in four treatment groups (Figure 1A). Cylinder test performance of mice transduced with rAAV-AS69+rAAV-GFP did not differ from mice transduced with rAAV-GFP only. In summary, mice transduced with rAAV-αSyn+rAAV-GFP resulted in a slower phenotype, which was somewhat restored in mice transduced with rAAV-αSyn+rAAV-AS69 (Figure 1B). Yet, this difference was not statistically significant.

3.2. rAAV-AS69 did Not Prevent the Degeneration of TH-positive Neurons

Next, we quantified the number of dopaminergic neurons in the SNc, using tyrosine hydroxylase (TH) as a marker (Figure 2A). Hemispheres transduced with rAAV-GFP showed a similar number of TH-positive neurons as uninjected hemispheres (Figure 2B), suggesting that transduction with rAAV-GFP did not cause degeneration of dopaminergic neurons in the SN. This also holds true for mice transduced with rAAV-GFP+rAAV-AS69. Mice transduced with rAAV-αSyn+rAAV-GFP showed about 20 % fewer TH-positive neurons than rAAV-GFP expressing animals. The reduction of dopaminergic neurons by αSyn overexpression is therefore modest, but it can explain the subtle motor deficit (Figure 1B). In mice transduced with rAAV-αSyn+rAAV-AS69, the loss of TH-positive neurons was similar as in mice with rAAV-αSyn+rAAV-GFP (Figure 2B). This suggests that transduction with rAAV-AS69 did not reduce the αSyn-induced loss of dopaminergic neurons in the SN.
We also analyzed the immunoreactivity of the long apical dendrites of dopaminergic neurons in the SNr (Figure 2C). Specifically, the mean signal intensity in the TH channel within the SNr was normalized to the contralateral hemisphere. Normalized TH immunoreactivity in the SNr was significantly reduced in mice transduced with rAAV-αSyn as compared to rAAV-GFP, but not significantly different between mice transduced with rAAV-αSyn+rAAV-AS69 compared to rAAV-αSyn+rAAV-GFP (Figure 2C), consistent with the observations in the SNc.
Moreover, we analyzed the density of TH-positive axon terminals in the striatum (Figure 2D & E). The density of TH-positive axon terminals in the striatum was expressed relative to the intact contralateral hemisphere. With this method, we did not observe degeneration of axon terminals with αSyn overexpression (Figure 2E).

3.3. Transduction with rAAV-AS69 Is Associated with an Increase in αSyn Pathology

In order to better understand the consequences of transducing SN neurons with human αSyn (H-αSyn), we stained for H-αSyn and quantified the extent of αSyn pathology using p-S129-αSyn (P-αSyn) as an established epitope for pathologically aggregated αSyn [23,24,25]. We acquired slide scans of the SN (Figure 3A) for quantification of positive cells (Figure 3C) as well as high magnification images of dopaminergic neurons co-stained for H-αSyn and P-αSyn (Figure 3B). The count of H-αSyn-positive cells seemed to be increased in animals transduced with rAAV-αSyn+rAAV-AS69 compared to animals transduced with rAAV-αSyn+rAAV-GFP (Figure 3C), albeit not reaching statistical significance. The average count of P-αSyn positive cells within the SNc was significantly higher in rAAV-αSyn+rAAV-AS69 expressing mice than in rAAV-αSyn+rAAV-GFP expressing mice (Figure 3C). The counts for H-αSyn positive cells and P-αSyn positive cells were highly correlated in mice transduced with rAAV-αSyn+rAAV-AS69 (Figure 3E). Mice transduced with rAAV-αSyn+rAAV-GFP did not show a significant correlation (Figure 3D).

3.4. Microglia Are Activated in Mice Transduced with rAAV-αSyn+rAAV-AS69

Microglia react to αSyn pathology and change morphology upon activation. Residential microglia exhibit long filigree dendrites (Figure 4A, green arrow). Hypertrophic and dystrophic microglia (Figure 4A, yellow arrow) are considered as an activated subtype [26]. Amoeboid and rod-shaped microglia are also associated with neurodegenerative processes and can be recognized by an enlarged soma with ruffled extensions (Figure 4A, cluster of microglia: red arrow). In order to investigate microglia activation in our model, we quantified and categorized Iba1-positive cells in the SNc.
In mice transduced with rAAV-GFP, rAAV-αSyn+rAAV-GFP or rAAV-αSyn+rAAV-AS69, microglia numbers were 50 % higher in the injected hemisphere as compared to the non-injected hemisphere (Figure 4B), suggesting that rAAV transduction caused microglia activation. With respect to microglia subtypes, about 75 % of Iba1-positive microglia in mice transduced with rAAV-GFP, rAAV-GFP+rAAV-AS69 or rAAV-αSyn+rAAV-GFP were in a residential state (Figure 4C). This fraction was lower in mice transduced with rAAV-αSyn+rAAV-AS69. Conversely, activated microglia subtypes were significantly more common in mice transduced with rAAV-αSyn+rAAV-AS69 than in the other groups (Figure 4C). The counts for hypertrophic microglia correlated with the count of P-αSyn positive cells (Figure 4D), consistent with the hypothesis that both effects might be causally linked.

3.4.1. Astrogliosis Is Not Attenuated by rAAV-AS69

In order to quantify astrogliosis, we determined the GFAP-positive area in the SNc (Figure 4A). In mice transduced with rAAV-GFP or rAAV-GFP+rAAV-AS69, the size of the GFAP-positive area in the injected hemispheres was similar to the uninjected hemisphere (Figure 4E). In mice transduced with rAAV-αSyn+rAAV-GFP or rAAV-αSyn+rAAV-AS69, the GFAP-positive area was 15% larger in the injected hemisphere than in the uninjected hemisphere (Figure 4E), suggesting that αSyn expression could contribute to astroglia activation whereas AS69 transduction had no effect.

4. Discussion

In this study, transduction of SN neurons with rAAV-AS69 did not reduce the extent of neurodegeneration resulting from overexpression of αSyn.
Transduction of SN neurons with rAAV-αSyn caused αSyn pathology as evidenced by cells positive for H-αSyn and P-αSyn (Figure 3C), degeneration of dopaminergic neurons (Figure 2B, C) and a motor phenotype (Figure 1B). This finding is consistent with previous studies using the same paradigm [22].
Even though aSyn transduction affected mouse behavior in the cylinder test (Figure 1B), we did not find statistically significant degeneration of TH-positive axon terminals in the striatum (Figure 2E). We attribute the discrepancy to the modest degree of degeneration of TH-positive neurons in the SNc compared to previous experiments [22]. Furthermore, due to the viral injection, we observed non-homogenous degeneration in the SNc. Transduction with rAAV-AS69 did not reduce αSyn pathology as quantified by H-αSyn and P-αSyn staining (Figure 3C); rAAV-AS69 also did not reduce degeneration of dopaminergic neurons (Figure 2B, C). These findings were unexpected given the reduction in αSyn pathology we observed with AS69 in cultured neurons and mouse striatum [5]. This discrepancy could be explained by the fact that in our previous work, αSyn pathology was induced by adding PFF, whereas αSyn overexpression was used in the present study. Collectively, these findings would be consistent with the hypothesis that PFF-based models do not capture the entire spectrum of αSyn pathology and toxicity. Along similar lines, we and others previously found that inhibiting αSyn fibril formation by synthetic mutations did not block αSyn toxicity [27,28,29].
Activation of microglia was most pronounced in mice transduced with rAAV-αSyn+ rAAV-AS69 (Figure 4C). Because constitutive activation of inflammatory changes is sufficient to induce αSyn pathology [30], this finding would be consistent with the hypothesis that any reduction in αSyn pathology resulting from AS69 co-expression was consumed by the consequences of increased neuroinflammation. rAAV-AS69 alone did not activate microglia (Figure 4C), suggesting that the microglial activation observed with rAAV-αSyn and rAAV-AS69 could result from an altered processing of αSyn with AS69. Indeed, neuronal expression of AS69 did rescue the behavioral phenotype of neuronal αSyn expression in drosophila [4], which lack neuroinflammation by microglia.
In addition, aSyn pathology spreading to non-neuronal cells might be rescued insufficiently by expression of rAAV-AS69 with the neuronal promoter we used. Administering AS69-protein showed beneficial effects in our previous experiment. Indeed, both small molecules UCB0599 and Anle138b that are currently tested in clinical trials have been administered systemically in various mouse models of synucleinopathies [31,32]. Oral Anle138b can prevent neurodegeneration and synuclein pathology in toxin-based models (rotenone, MPTP) as well as in transgenic mice that overexpress human A30P-α-Syn [33,34]. UCB0599 has been studied in transgenic mice that overexpress human wild type αSyn using intraperitoneal injections five times a week for 3 months [35]. Neither Anle138b nor UCB0599 have been tested in a model that is based on viral overexpression of αSyn.
AS69 is a dimer of two identical subunits that are covalently linked by a disulfide bond between the Cys-28 residues of both subunits. Intracellular disulfide bonds can be reduced in a disulfide exchange reaction with excess glutathione [36]. Potentially, the dimeric complex can therefore not be formed. On the other hand we previously observed protective effects of protein AS69 in PFF-mice [5]. Since protein AS69 was also in the extracellular space after intrastriatal injection, it also acted in an additional environment compared to rAAV-AS69 while being exposed to less redox reactions extracellularly.
In summary, AS69′s inhibitory activities did not reduce neurodegeneration in the SN resulting from αSyn overexpression, highlighting the importance to test potential therapeutic strategies in a variety of animal models and to study multiple routes of drug administration.

Author Contributions

Conceptualization, LH, ES, KW, VB, WH, BF. Investigation, LH, ES, WK, CvdH, DG. Writing – first draft, LH, BF. Writing - Review & Editing, all authors.

Funding

This research received no external funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Ethics Committee of Technische Universität Dresden (protocol code 25-5131/496/15).

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We thank Andrea Kempe, Anett Böhme and Sylvia Kanzler for their excellent technical assistance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Obeso, J.A.; Stamelou, M.; Goetz, C.G.; Poewe, W.; Lang, A.E.; Weintraub, D.; Burn, D.; Halliday, G.M.; Bezard, E.; Przedborski, S.; et al. Past, Present, and Future of Parkinson’s Disease: A Special Essay on the 200th Anniversary of the Shaking Palsy. Mov. Disord. 2017, 32, 1264–1310. [CrossRef]
  2. Rodger, A.T.; ALNasser, M.; Carter, W.G. Are Therapies That Target α-Synuclein Effective at Halting Parkinson’s Disease Progression? A Systematic Review. Int. J. Mol. Sci. 2023, 24, 11022. [CrossRef]
  3. Mirecka, E.A.; Shaykhalishahi, H.; Gauhar, A.; Akgül, Ş.; Lecher, J.; Willbold, D.; Stoldt, M.; Hoyer, W. Sequestration of a β-Hairpin for Control of α-Synuclein Aggregation. Angew. Chem. Int. Ed. 2014, 53, 4227–4230. [CrossRef]
  4. Agerschou, E.D.; Flagmeier, P.; Saridaki, T.; Galvagnion, C.; Komnig, D.; Heid, L.; Prasad, V.; Shaykhalishahi, H.; Willbold, D.; Dobson, C.M.; et al. An Engineered Monomer Binding-Protein for α-Synuclein Efficiently Inhibits the Proliferation of Amyloid Fibrils. eLife 2019, 8, e46112. [CrossRef]
  5. Szegő, É.M.; Boß, F.; Komnig, D.; Gärtner, C.; Höfs, L.; Shaykhalishahi, H.; Wördehoff, M.M.; Saridaki, T.; Schulz, J.B.; Hoyer, W.; et al. A β-Wrapin Targeting the N-Terminus of α-Synuclein Monomers Reduces Fibril-Induced Aggregation in Neurons. Front. Neurosci. 2021, 15.
  6. Luk, K.C.; Kehm, V.; Carroll, J.; Zhang, B.; O’Brien, P.; Trojanowski, J.Q.; Lee, V.M.-Y. Pathological α-Synuclein Transmission Initiates Parkinson-like Neurodegeneration in Nontransgenic Mice. Science 2012, 338, 949–953. [CrossRef]
  7. Brundin, P.; Melki, R. Prying into the Prion Hypothesis for Parkinson’s Disease. J. Neurosci. Off. J. Soc. Neurosci. 2017, 37, 9808–9818. [CrossRef]
  8. Stopschinski, B.E.; Diamond, M.I. The Prion Model for Progression and Diversity of Neurodegenerative Diseases. Lancet Neurol. 2017, 16, 323–332. [CrossRef]
  9. Singleton, A.B.; Farrer, M.; Johnson, J.; Singleton, A.; Hague, S.; Kachergus, J.; Hulihan, M.; Peuralinna, T.; Dutra, A.; Nussbaum, R.; et al. α-Synuclein Locus Triplication Causes Parkinson’s Disease. Science 2003, 302, 841–841. [CrossRef]
  10. Maraganore, D.M. Collaborative Analysis of α-Synuclein Gene Promoter Variability and Parkinson Disease. JAMA 2006, 296, 661. [CrossRef]
  11. Nalls, W. Imputation of Sequence Variants for Identification of Genetic Risks for Parkinson’s Disease: A Meta-Analysis of Genome-Wide Association Studies. The Lancet 2011, 377, 641–649. [CrossRef]
  12. Van der Perren, A.; Toelen, J.; Carlon, M.; Van den Haute, C.; Coun, F.; Heeman, B.; Reumers, V.; Vandenberghe, L.H.; Wilson, J.M.; Debyser, Z.; et al. Efficient and Stable Transduction of Dopaminergic Neurons in Rat Substantia Nigra by rAAV 2/1, 2/2, 2/5, 2/6.2, 2/7, 2/8 and 2/9. Gene Ther. 2011, 18, 517–527. [CrossRef]
  13. Stirling, D.R.; Swain-Bowden, M.J.; Lucas, A.M.; Carpenter, A.E.; Cimini, B.A.; Goodman, A. CellProfiler 4: Improvements in Speed, Utility and Usability. BMC Bioinformatics 2021, 22, 433. [CrossRef]
  14. Bankhead, P.; Loughrey, M.B.; Fernández, J.A.; Dombrowski, Y.; McArt, D.G.; Dunne, P.D.; McQuaid, S.; Gray, R.T.; Murray, L.J.; Coleman, H.G.; et al. QuPath: Open Source Software for Digital Pathology Image Analysis. Sci. Rep. 2017, 7, 16878. [CrossRef]
  15. Lier, J.; Streit, W.J.; Bechmann, I. Beyond Activation: Characterizing Microglial Functional Phenotypes. Cells 2021, 10, 2236. [CrossRef]
  16. Wickham, H.; Averick, M.; Bryan, J.; Chang, W.; McGowan, L.; François, R.; Grolemund, G.; Hayes, A.; Henry, L.; Hester, J.; et al. Welcome to the Tidyverse. J. Open Source Softw. 2019, 4, 1686. [CrossRef]
  17. Kassambara, A. Ggpubr: “ggplot2” Based Publication Ready Plots; 2023;
  18. Ogle, D.H.; Doll, J.C.; Wheeler, A.P.; Dinno, A. FSA: Simple Fisheries Stock Assessment Methods; 2023;
  19. Wilke, C.; Fox, S.J.; Bates, T.; Manalo, K.; Lang, B.; Barrett, M.; Stoiber, M.; Philipp, A.; Denney, B.; Hesselberth, J.; et al. Wilkelab/Cowplot: 1.1.1 2021.
  20. Wei, T.; Simko, V. R Package “Corrplot”: Visualization of a Correlation Matrix; 2021;
  21. Kolde, R.; others Pheatmap: Pretty Heatmaps. R Package Version 2012, 1, 726.
  22. Oliveras-Salvá, M.; Van der Perren, A.; Casadei, N.; Stroobants, S.; Nuber, S.; D’Hooge, R.; Van den Haute, C.; Baekelandt, V. rAAV2/7 Vector-Mediated Overexpression of Alpha-Synuclein in Mouse Substantia Nigra Induces Protein Aggregation and Progressive Dose-Dependent Neurodegeneration. Mol. Neurodegener. 2013, 8, 44. [CrossRef]
  23. Mahul-Mellier, A.-L.; Burtscher, J.; Maharjan, N.; Weerens, L.; Croisier, M.; Kuttler, F.; Leleu, M.; Knott, G.W.; Lashuel, H.A. The Process of Lewy Body Formation, Rather than Simply α-Synuclein Fibrillization, Is One of the Major Drivers of Neurodegeneration. Proc. Natl. Acad. Sci. 2020, 117, 4971–4982. [CrossRef]
  24. Weston, L.J.; Stackhouse, T.L.; Spinelli, K.J.; Boutros, S.W.; Rose, E.P.; Osterberg, V.R.; Luk, K.C.; Raber, J.; Weissman, T.A.; Unni, V.K. Genetic Deletion of Polo-like Kinase 2 Reduces Alpha-Synuclein Serine-129 Phosphorylation in Presynaptic Terminals but Not Lewy Bodies. J. Biol. Chem. 2021, 296, 100273. [CrossRef]
  25. Fujiwara, H.; Hasegawa, M.; Dohmae, N.; Kawashima, A.; Masliah, E.; Goldberg, M.S.; Shen, J.; Takio, K.; Iwatsubo, T. α-Synuclein Is Phosphorylated in Synucleinopathy Lesions. Nat. Cell Biol. 2002, 4, 160–164. [CrossRef]
  26. Savage, J.C.; Carrier, M.; Tremblay, M.-È. Morphology of Microglia Across Contexts of Health and Disease. In Microglia; Garaschuk, O., Verkhratsky, A., Eds.; Methods in Molecular Biology; Springer New York: New York, NY, 2019; Vol. 2034, pp. 13–26 ISBN 978-1-4939-9657-5.
  27. Karpinar, D.P.; Balija, M.B.G.; Kügler, S.; Opazo, F.; Rezaei-Ghaleh, N.; Wender, N.; Kim, H.-Y.; Taschenberger, G.; Falkenburger, B.H.; Heise, H.; et al. Pre-Fibrillar Alpha-Synuclein Variants with Impaired Beta-Structure Increase Neurotoxicity in Parkinson’s Disease Models. EMBO J. 2009, 28, 3256–3268. [CrossRef]
  28. Prots, I.; Grosch, J.; Brazdis, R.-M.; Simmnacher, K.; Veber, V.; Havlicek, S.; Hannappel, C.; Krach, F.; Krumbiegel, M.; Schütz, O.; et al. α-Synuclein Oligomers Induce Early Axonal Dysfunction in Human iPSC-Based Models of Synucleinopathies. Proc. Natl. Acad. Sci. U. S. A. 2018, 115, 7813–7818. [CrossRef]
  29. Rockenstein, E.; Nuber, S.; Overk, C.R.; Ubhi, K.; Mante, M.; Patrick, C.; Adame, A.; Trejo-Morales, M.; Gerez, J.; Picotti, P.; et al. Accumulation of Oligomer-Prone α-Synuclein Exacerbates Synaptic and Neuronal Degeneration in Vivo. Brain J. Neurol. 2014, 137, 1496–1513. [CrossRef]
  30. Szego, E.M.; Malz, L.; Bernhardt, N.; Rösen-Wolff, A.; Falkenburger, B.H.; Luksch, H. Constitutively Active STING Causes Neuroinflammation and Degeneration of Dopaminergic Neurons in Mice. eLife 2022, 11, e81943. [CrossRef]
  31. Smit, J.W.; Basile, P.; Prato, M.K.; Detalle, L.; Mathy, F.; Schmidt, A.; Lalla, M.; Germani, M.; Domange, C.; Biere, A.; et al. Phase 1/1b Studies of UCB0599 , an Oral Inhibitor of α-Synuclein Misfolding, Including a Randomized Study in Parkinson’s Disease. Mov. Disord. 2022, 37, 2045–2056. [CrossRef]
  32. Levin, J.; Sing, N.; Melbourne, S.; Morgan, A.; Mariner, C.; Spillantini, M.G.; Wegrzynowicz, M.; Dalley, J.W.; Langer, S.; Ryazanov, S.; et al. Safety, Tolerability and Pharmacokinetics of the Oligomer Modulator Anle138b with Exposure Levels Sufficient for Therapeutic Efficacy in a Murine Parkinson Model: A Randomised, Double-Blind, Placebo-Controlled Phase 1a Trial. eBioMedicine 2022, 80, 104021. [CrossRef]
  33. Wagner, J.; Ryazanov, S.; Leonov, A.; Levin, J.; Shi, S.; Schmidt, F.; Prix, C.; Pan-Montojo, F.; Bertsch, U.; Mitteregger-Kretzschmar, G.; et al. Anle138b: A Novel Oligomer Modulator for Disease-Modifying Therapy of Neurodegenerative Diseases Such as Prion and Parkinson’s Disease. Acta Neuropathol. (Berl.) 2013, 125, 795–813. [CrossRef]
  34. Levin, J.; Schmidt, F.; Boehm, C.; Prix, C.; Bötzel, K.; Ryazanov, S.; Leonov, A.; Griesinger, C.; Giese, A. The Oligomer Modulator Anle138b Inhibits Disease Progression in a Parkinson Mouse Model Even with Treatment Started after Disease Onset. Acta Neuropathol. (Berl.) 2014, 127, 779–780. [CrossRef]
  35. Price, D.L.; Khan, A.; Angers, R.; Cardenas, A.; Prato, M.K.; Bani, M.; Bonhaus, D.W.; Citron, M.; Biere, A.-L. In Vivo Effects of the Alpha-Synuclein Misfolding Inhibitor Minzasolmin Supports Clinical Development in Parkinson’s Disease. Npj Park. Dis. 2023, 9, 114. [CrossRef]
  36. Bechtel, T.J.; Weerapana, E. From Structure to Redox: The Diverse Functional Roles of Disulfides and Implications in Disease. PROTEOMICS 2017, 17, 1600391. [CrossRef]
Preprints 105521 g001
Figure 1. rAAV-AS69 Does Not Attenuate αSyn-Induced Motor Deficits. (A) Illustration of the study design. Five mice per group received unilateral stereotactic injections of rAAV into the substantia nigra pars compacta. Cylinder tests were performed one day before injection, four weeks and eight weeks after injection of rAAV. (B) Time required for 25 forepaw contacts with the glass cylinder wall at the indicated time points. Markers represent individual animals (n = 5). Kruskal-Wallis multiple comparisons followed by Dunn’s test as post-hoc-analysis, p-values were adjusted with the Holm method (p-value = 0.016 for GFP vs. αSyn+GFP).
Figure 1. rAAV-AS69 Does Not Attenuate αSyn-Induced Motor Deficits. (A) Illustration of the study design. Five mice per group received unilateral stereotactic injections of rAAV into the substantia nigra pars compacta. Cylinder tests were performed one day before injection, four weeks and eight weeks after injection of rAAV. (B) Time required for 25 forepaw contacts with the glass cylinder wall at the indicated time points. Markers represent individual animals (n = 5). Kruskal-Wallis multiple comparisons followed by Dunn’s test as post-hoc-analysis, p-values were adjusted with the Holm method (p-value = 0.016 for GFP vs. αSyn+GFP).
Preprints 105521 g001
Preprints 105521 g002
Figure 2. rAAV-AS69 Does Not Prevent the Degeneration of TH-positive Neurons and Dendrites. (A) Representative images for each cohort of the SN in the mouse midbrain stained for TH. Cohorts as indicated above each image. SNc is highlighted in blue and SNr is encircled in red. Scale bar: 200 μm. (B) Number of TH-positive neurons in the substantia nigra pars compacta relative to the contralateral hemisphere. Markers represent individual animals (n = 5). Comparison by two-way ANOVA (p-value = 0.0001 for factor GFP vs. αSyn; p-value = 0.942 for factor GFP vs. AS69). (C) Signal intensity of TH-positive dendrites in the SNr relative to the contralateral hemisphere. Labels represent individual animals (n = 5 animals per group). Comparison by two-way ANOVA (p-value = 0.0002 for GFP vs. αSyn; p-value = 0.531 for GFP vs. AS69). (D) Representative images of dopaminergic axon terminals, labeled by TH, in the striatum. Cohorts as indicated above each image. Scale bar: 30 μm. (E) TH-positive area in the striatum relative to the contralateral hemisphere. Markers represent individual animals (n = 5 for αSyn+GFP and n = 3 for αSyn+AS69).
Figure 2. rAAV-AS69 Does Not Prevent the Degeneration of TH-positive Neurons and Dendrites. (A) Representative images for each cohort of the SN in the mouse midbrain stained for TH. Cohorts as indicated above each image. SNc is highlighted in blue and SNr is encircled in red. Scale bar: 200 μm. (B) Number of TH-positive neurons in the substantia nigra pars compacta relative to the contralateral hemisphere. Markers represent individual animals (n = 5). Comparison by two-way ANOVA (p-value = 0.0001 for factor GFP vs. αSyn; p-value = 0.942 for factor GFP vs. AS69). (C) Signal intensity of TH-positive dendrites in the SNr relative to the contralateral hemisphere. Labels represent individual animals (n = 5 animals per group). Comparison by two-way ANOVA (p-value = 0.0002 for GFP vs. αSyn; p-value = 0.531 for GFP vs. AS69). (D) Representative images of dopaminergic axon terminals, labeled by TH, in the striatum. Cohorts as indicated above each image. Scale bar: 30 μm. (E) TH-positive area in the striatum relative to the contralateral hemisphere. Markers represent individual animals (n = 5 for αSyn+GFP and n = 3 for αSyn+AS69).
Preprints 105521 g002
Preprints 105521 g003
Figure 3. Transduction with rAAV-AS69 is associated with an increase in P-αSyn pathology. (A) Representative images of the SN showing GFP (yellow), H-αSyn (gray) and P-αSyn (magenta). Cohorts as indicated by row labels. Scale bar: 200 μm. (B) Representative high magnification images of TH-positive neurons (yellow), H-αSyn (gray) and P-αSyn (magenta) in the SNc. Scale bar: 30 μm. (C) Count of positive cells for H-αSyn and P-αSyn in the SNc of indicated cohorts. Markers represent individual animals (n = 5). Comparisons by Wilcoxon rank sum exact test (p-value = 0.22 for H-αSyn; p-value = 0.032 for P-αSyn). (D) Correlation of H-αSyn and P-αSyn in animals transduced with rAAV-αSyn and rAAV-GFP. Pearson’s product-moment correlation at r = 0.63 (p-value = 0.26). Dotted line represents the corresponding generalized linear mixed model. (D) Correlation of H-αSyn and P-αSyn in animals transduced with rAAV-αSyn and rAAV-AS69. Pearson’s product-moment correlation at r = 0.93 (p-value = 0.021). Dotted line represents the corresponding generalized linear mixed model.
Figure 3. Transduction with rAAV-AS69 is associated with an increase in P-αSyn pathology. (A) Representative images of the SN showing GFP (yellow), H-αSyn (gray) and P-αSyn (magenta). Cohorts as indicated by row labels. Scale bar: 200 μm. (B) Representative high magnification images of TH-positive neurons (yellow), H-αSyn (gray) and P-αSyn (magenta) in the SNc. Scale bar: 30 μm. (C) Count of positive cells for H-αSyn and P-αSyn in the SNc of indicated cohorts. Markers represent individual animals (n = 5). Comparisons by Wilcoxon rank sum exact test (p-value = 0.22 for H-αSyn; p-value = 0.032 for P-αSyn). (D) Correlation of H-αSyn and P-αSyn in animals transduced with rAAV-αSyn and rAAV-GFP. Pearson’s product-moment correlation at r = 0.63 (p-value = 0.26). Dotted line represents the corresponding generalized linear mixed model. (D) Correlation of H-αSyn and P-αSyn in animals transduced with rAAV-αSyn and rAAV-AS69. Pearson’s product-moment correlation at r = 0.93 (p-value = 0.021). Dotted line represents the corresponding generalized linear mixed model.
Preprints 105521 g003
Preprints 105521 g004
Figure 4. Microglia Are Particularly Activated in αSyn+AS69 Expressing Mice and Astrogliosis Is Not Attenuated by rAAV-AS69. (A) Representative images of Iba1-positive microglia and GFAP-positive astrocytes in the SNc. Ramified microglia are highlighted by the green arrows. The yellow arrow indicates a dystrophic microglial cell and the red arrow shows an amoeboid microglial cell. Cohorts as indicated by the column names. Scale bar: 40 μm. Individual cells are shown in higher magnification in the inset. Scale bar: 20 μm. (B) Count of Iba1-positive cells in the SNc relative to the contralateral hemisphere. Labels represent individual animals (n = 5). Comparisons by multiple t-tests, p-values were adjusted for multiple testing with the Benjamini-Hochberg method (p-value= 0.04 for GFP vs. GFP+AS69, GFP+AS69 vs. αSyn+GFP and GFP+AS69 vs. αSyn+AS69 ). (C) Occurrence of each morphological microglia phenotype relative to the total count of microglia in the injected SNc. Labels represent individual animals (n = 5). Comparisons by two-way ANOVA followed by Tukey’s HSD (Dendrite-like: p-value = 0.0008 for GFP vs. αSyn, p-value = 0.001 for GFP vs. AS69, p-value = 0.004 for the interaction; post-hoc-analysis: p-value = 0.0001 for GFP vs. αSyn+AS69, p-value = 0.0004 for GFP+AS69 vs. αSyn+AS69, p-value = 0.0005 for αSyn+GFP vs. αSyn+AS69, Hypertrophic or dystrophic: p-value = 0.001 for GFP vs. αSyn, p-value = 0.001 for GFP vs. AS69, p-value = 0.007 for the interaction; post-hoc-analysis: p-value = 0.0003 for GFP vs. αSyn+AS69, p-value = 0.0007 for GFP+AS69 vs. αSyn+AS69, p-value = 0.0009 for αSyn+GFP vs. αSyn+AS69, Amoeboid or Rod-shaped: p-value = 0.029 for GFP vs. αSyn, p-value = 0.027 for GFP vs. AS69, p-value = 0.0256 for the interaction; post-hoc-analysis: p-value = 0.017 for GFP vs. αSyn+AS69, p-value = 0.016 for GFP+AS69 vs. αSyn+AS69, p-value = 0.015 for αSyn+GFP vs. αSyn+AS69). (D) Correlation of the count of P-αSyn positive cells in the SN and microglia classified as ”hypertrophic or dystrophic microglia” in the SN. Pearson’s product-moment correlation r = 0.71 (p-value = 0.022). Dotted line represents the corresponding generalized linear mixed model. (E) GFAP positive area in the SNc relative to the contralateral hemisphere. Labels represent individual animals (n = 5). Kruskal-Wallis multiple comparisons followed by Dunn’s test as post-hoc-analysis, p-values were adjusted with the Holm method (p-value = 0.009 for GFP vs. αSyn+GFP, p-value = 0.023 for GFP vs. αSyn+AS69).
Figure 4. Microglia Are Particularly Activated in αSyn+AS69 Expressing Mice and Astrogliosis Is Not Attenuated by rAAV-AS69. (A) Representative images of Iba1-positive microglia and GFAP-positive astrocytes in the SNc. Ramified microglia are highlighted by the green arrows. The yellow arrow indicates a dystrophic microglial cell and the red arrow shows an amoeboid microglial cell. Cohorts as indicated by the column names. Scale bar: 40 μm. Individual cells are shown in higher magnification in the inset. Scale bar: 20 μm. (B) Count of Iba1-positive cells in the SNc relative to the contralateral hemisphere. Labels represent individual animals (n = 5). Comparisons by multiple t-tests, p-values were adjusted for multiple testing with the Benjamini-Hochberg method (p-value= 0.04 for GFP vs. GFP+AS69, GFP+AS69 vs. αSyn+GFP and GFP+AS69 vs. αSyn+AS69 ). (C) Occurrence of each morphological microglia phenotype relative to the total count of microglia in the injected SNc. Labels represent individual animals (n = 5). Comparisons by two-way ANOVA followed by Tukey’s HSD (Dendrite-like: p-value = 0.0008 for GFP vs. αSyn, p-value = 0.001 for GFP vs. AS69, p-value = 0.004 for the interaction; post-hoc-analysis: p-value = 0.0001 for GFP vs. αSyn+AS69, p-value = 0.0004 for GFP+AS69 vs. αSyn+AS69, p-value = 0.0005 for αSyn+GFP vs. αSyn+AS69, Hypertrophic or dystrophic: p-value = 0.001 for GFP vs. αSyn, p-value = 0.001 for GFP vs. AS69, p-value = 0.007 for the interaction; post-hoc-analysis: p-value = 0.0003 for GFP vs. αSyn+AS69, p-value = 0.0007 for GFP+AS69 vs. αSyn+AS69, p-value = 0.0009 for αSyn+GFP vs. αSyn+AS69, Amoeboid or Rod-shaped: p-value = 0.029 for GFP vs. αSyn, p-value = 0.027 for GFP vs. AS69, p-value = 0.0256 for the interaction; post-hoc-analysis: p-value = 0.017 for GFP vs. αSyn+AS69, p-value = 0.016 for GFP+AS69 vs. αSyn+AS69, p-value = 0.015 for αSyn+GFP vs. αSyn+AS69). (D) Correlation of the count of P-αSyn positive cells in the SN and microglia classified as ”hypertrophic or dystrophic microglia” in the SN. Pearson’s product-moment correlation r = 0.71 (p-value = 0.022). Dotted line represents the corresponding generalized linear mixed model. (E) GFAP positive area in the SNc relative to the contralateral hemisphere. Labels represent individual animals (n = 5). Kruskal-Wallis multiple comparisons followed by Dunn’s test as post-hoc-analysis, p-values were adjusted with the Holm method (p-value = 0.009 for GFP vs. αSyn+GFP, p-value = 0.023 for GFP vs. αSyn+AS69).
Preprints 105521 g004
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.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated