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

# Nigral ATP13A2 depletion induces Parkinson's disease-related neurodegeneration in a pilot study in non-human primates

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

Check for updates

Joanna Sikora<fnref n="1" />,²,⁴, Sandra Dovero<fnref n="1" />,⁴, Rémi Kinet❶, Marie-Laure Arotcarena❶, Sylvain Bohic❶,⁴, Erwan Bezard❶, Pierre-Olivier Fernagut❷ & Benjamin Dehay❶,⁴

Lysosomal impairment is strongly implicated in Parkinson's disease (PD). Among the several PD-linked genes, the *ATP13A2* gene, associated with the PARK9 locus, encodes a transmembrane lysosomal P5-type ATPase. Mutations in the *ATP13A2* gene were primarily identified as the cause of Kufor-Rakeb syndrome (KRS), a juvenile-onset form of PD. Subsequently, an increasing list of several mutations has been described. These mutations result in truncation of the ATP13A2 protein, leading to a loss of function but surprisingly causing heterogeneity and variability in the clinical symptoms associated with different brain pathologies. In vitro studies show that its loss compromises lysosomal function, contributing to cell death. To understand the role of ATP13A2 dysfunction in disease, we disrupted its expression through a viral vector-based approach in nonhuman primates. Here, in this pilot study, we injected bilaterally into the substantia nigra of macaques, a lentiviral vector expressing an ATP13A2 small hairpin RNA. Animals were terminated five months later, and brains were harvested and compared with historical non-injected control brains to evaluate cerebral pathological markers known to be affected in KRS and PD. We characterised the pattern of dopaminergic loss in the striatum and the substantia nigra, the regional distribution of α-synuclein immunoreactivity in several brain structures, and its pathological status (i.e., S129 phosphorylation), the accumulation of heavy metals in nigral sections and occurrence of lysosomal dysfunction. This proof-of-concept experiment highlights the potential value of lentivirus-mediated ATP13A2 silencing to induce significant and ongoing degeneration in the nigrostriatal pathway, α-synuclein pathology, and iron accumulation in nonhuman primates.

The *ATP13A2* gene (OMIM#610513) encodes a transmembrane lysosomal P5-type ATPase. The first mutations described in *ATP13A2* cause an autosomal recessive form of early-onset parkinsonism called Kufor-Rakeb Syndrome (KRS, OMIM#606693)<fnref n="1" /> and neuronal ceroid lipofuscinosis (NCL, OMIM#256730), later assigned to PARK9 locus. The clinical spectrum of KRS comprises parkinsonism, dystonia, supranuclear gaze palsy, and severe cognitive decline². Brain MRI of KRS patients revealed generalised atrophy and putaminal and caudate iron accumulation, classifying KRS amongst neurodegeneration with brain iron accumulation (NBIA)³.

NCL is a lysosomal storage disorder (LSD) that belongs to a genetically heterogeneous group of rare neurodegenerative diseases resulting from an accumulation of lipopigments, predominantly in the brain and retina. Typical clinical manifestations of these LSDs include progressive motor and cognitive impairment, seizures, and visual impairment or blindness. The outcome is invariably fatal, usually during the second or third decade.

Today, there is a growing list of homozygous and compound-heterozygous mutations⁴,⁵ linked to the truncation of the ATP13A2 protein, resulting in a loss of its functional capabilities. Remarkably, all mutations

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

e-mail: benjamin.dehay@u-bordeaux.fr

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

Parkinson's
Foundation

<footnote n="1">¹Univ. Bordeaux, CNRS, IMN, Bordeaux, France. ²Univ. De Poitiers, INSERM, LNEC, Poitiers, France. ³Univ. Grenoble Alpes, Synchrotron Radiation for Biomedicine (STROBE), Grenoble, France. ⁴These authors contributed equally: Joanna Sikora, Sandra Dovero.</footnote>
associated with KRS that have undergone experimental scrutiny exhibit varying degrees of interference with ATP13A2 function. These disruptions are brought about through diverse mechanisms, including nonsense-mediated messenger RNA decay, protein mislocalization, and premature degradation of protein products by the proteasomal system. However, the phenotypic spectrum associated with genetic variants in *ATP13A2* now becomes more and more complex since the (partial) loss of function leads to a whole range of clinical phenotypes, including previously comprised of KRS, NCL, multiple system atrophy (MSA)⁶, hereditary spastic paraplegia (SPG78)⁷⁻¹⁰, and juvenile-onset amyotrophic lateral sclerosis (ALS)¹¹, as supported by both genetic and functional data, which makes the *ATP13A2* gene a puzzling gene¹². To date, neuroprotective strategies for all those diseases that would stop or slow down the yet unrelenting ATP13A2-related degenerative process are eagerly awaited.

Mechanistic studies have previously reported that KRS-linked mutations in *ATP13A2* lead to several lysosomal alterations in *ATP13A2* KRS patient-derived fibroblasts, including impaired lysosomal acidification, decreased proteolytic processing of lysosomal enzymes, reduced degradation of lysosomal substrates and diminished lysosomal-mediated clearance of autophagosomes¹³⁻¹⁶. In addition, ATP13A2 levels are decreased in nigral dopaminergic neurons from sporadic Parkinson's disease (PD) patients¹³, a result later corroborated by others¹⁷. Overall, those results unravel an instrumental role for ATP13A2 on lysosomal function and cell viability. Recently, the cellular function of this P5-type ATPase was unmasked¹⁶ as a lysosomal H⁺, K⁺-ATPase, and polyamine exporter, confirming its link to protein degradation via the lysosomal machinery¹³⁻¹⁵,¹⁸ and metal homeostasis (as reviewed in¹⁹), which cryo-electron microscopy has further structurally characterised²⁰⁻²⁴, altogether improving our understanding behind ATP13A2-associated functions.

To date, only one human brain neuropathological study from one ATP13A2-mutant KRS patient is available², reporting widespread neuronal and glial lipofuscin accumulation, no Lewy-body pathology and α-synuclein pathology and sparse iron deposits in several brain areas, with skin and muscle containing electron-dense lamellated structures in this patient²⁵. To decipher the pathological consequences of ATP13A2 deficiency, a vast spectrum of model organisms has been used, from yeast²⁶ to different animal species (fly²⁷, worm²⁸, zebrafish²⁹, and mouse³⁰,³¹). ATP13a2 null mice exhibit increased autofluorescence in multiple brain regions, indicative of lipofuscin accumulation and modest age-related motor dysfunction preceded by neuropathological changes, including gliosis, accumulation of ubiquitinated protein aggregates, and endolysosomal abnormalities. Still, no overt neurodegeneration and iron accumulation have been observed. Additionally, a recessive mutation has been described in canine ATP13A2, causing adult-onset NCL in Tibetan terriers where autofluorescent cytoplasmic inclusions were detected in brain tissue consistent with NCL³². Altogether, model organisms used until now with the (partial) loss of ATP13A2 exhibit subtle age-dependent motor dysfunction, which at the cellular level correlates to lysosomal impairment. In mammalian model organisms, gliosis and lipofuscinosis were additional from ATP13A2 deficiency. These findings demonstrate the significance of lysosomal functionality at the core of ATP13A2 depletion. However, to truly decipher the clinical spectrum observed in humans, the role of ATP13A2 in neuronal survival, and how this dysfunction leads to neurodegeneration, we hypothesised that, in nonhuman primates, a species phylogenetically closer to humans, targeted silencing of ATP13A2 in dopaminergic neurons might result in neuronal death, α-synuclein pathology, metal dyshomeostasis, and autophagy dysfunction.

To this end, in a proof-of-concept pilot study, we injected bilaterally into the substantia nigra (SN) of macaques a lentiviral vector-mediated short-hairpin RNA (LV-shRNA) to downregulate the expression of ATP13A2 under the control of a ubiquitous promoter. Five months after administration, we observed nigrostriatal neurodegeneration associated with α-synuclein pathology, nigral iron accumulation, and autophagy lysosomal pathway (ALP) disturbances upon LV-mediated ATP13A2 silencing.

## Results

### Nigral ATP13A2 depletion induces nigrostriatal neurodegeneration

The lentiviral vector (LV) system is highly effective for delivering short hairpin RNA (shRNA) expression cassettes³³ associated with no indications of increased non-human primate (NHP) tissue damage or inflammatory responses³⁴. We used a functional and efficient sequence of shRNA targeting human ATP13A2 cloned into a LV plasmid, which displayed a 95% reduction in ATP13A2 levels by immunoblotting, as previously reported¹³. In addition, we showed that LV-mediated ATP13A2 knockdown in primary mesencephalic dopaminergic neurons resulted in neurotoxicity¹³. After confirming that lentiviral injection of non-targeting (Scramble) shRNA did not induce dopaminergic neurodegeneration in SN of macaques five months after injection (Supplementary Fig. 1), we decided to give priority to the LV encoding shATP13A2 to be injected in both hemispheres to reduce animal use. To further determine the pathological consequences of silencing ATP13A2 in a species closer to humans, adult female macaques received bilateral stereotactic injections (10 µl per hemisphere) of LV encoding shATP13A2 into the SN (Supplementary Fig. 2a). Tissues from 4 to 6 additional age- and sex-matched non-injected historical control animals were processed, stored and used under the same conditions for post-mortem experiments³⁵,³⁶. Five months after administration, macaques were tested to evaluate the integrity of the nigrostriatal pathway. ATP13A2 immunofluorescence revealed a decrease in ATP13A2 levels in Tyrosine Hydroxylase (TH)-positive dopaminergic neurons of the shATP13A2 versus the control group (Supplementary Fig. 2b). We observed a downward trend in protein levels, with an estimated reduction of around ~60% in endogenous ATP13A2 five months after shATP13A2 LV injection obtained by immunoblot of homogenised formalin-fixed tissues collected on sections from the SNpc (Supplementary Fig. 2c). Stereological counts showed that shATP13A2-injected animals exhibited a ~28% TH-positive cell loss in the SNpc (Fig. 1a) with a typical morphology of dying dopaminergic neurons displaying shrunken eosinophilic cytoplasm and cytoplasmic vacuolation, as signs of oxidative stress and autophagic degeneration³⁷,³⁸ (Supplementary Fig. 3). However, no overt parkinsonism was observed because the extent of the lesion remained below the threshold for the onset of motor symptoms, i.e., 45% cell loss³⁹. At the striatal level, LV-shATP13A2 injection induced a significant 30% loss of TH immunoreactivity in the caudate nucleus (Fig. 1b) and 45% in the putamen (Fig. 1b). AADC, the enzyme that decarboxylates L-DOPA into dopamine (Fig. 1b), and striatal DAT immunostaining (Fig. 1b) were also both reduced, with a mean difference between groups of 60-70% (AADC) and 35% (DAT) in the caudate and putamen respectively.

### ATP13A2 depletion alters α-synuclein homeostasis in substantia nigra

We then investigated the consequences of ATP13A2 depletion-induced nigrostriatal degeneration on nigral α-synuclein (α-syn) protein expression. We performed a quantitative post-mortem analysis of total α-syn (syn211 clone antibody) and phosphorylated α-syn at S129 (EP1536Y clone antibody) levels as a surrogate marker of the presence of α-syn pathology⁴⁰,⁴¹ within the SN. In contrast to predictions from in vitro studies, we observed no changes in total α-syn-positive immunostaining in the SN (Fig. 2a), while the α-syn pathogenic form, i.e., S129 phosphorylated α-syn, increased in DA neurons in shATP13A2-injected animals (Fig. 2b). Even though we could not obtain fresh tissues from the substantia nigra of macaques for Western blotting analysis, we performed protein extraction from formalin-fixed substantia nigra to further investigate α-syn pathology. Notably, our findings indicated that there were no discernible changes in total α-syn levels when using the syn211 antibody on mesencephalic total protein extracts, aligning with the results obtained from the immunohistochemical analysis, which also showed a significant accumulation of phosphorylated α-synuclein (pSyn) (Fig. 2c). The absence of modifications in total α-syn levels could potentially be attributed to (i) an increase in its phosphorylation state and/or (ii) a conformational alteration in the α-syn protein, causing it to adopt a pathological form that conceals the syn211 epitope (located within
<figure><img src="image_5.png" /><figcaption>**Fig. 1** | Decreased ATP13A2 levels in midbrain dopamine neurons lead to a PD-like pattern of nigrostriatal degeneration. **a** Representative images (top) and quantification (bottom) of Tyrosine Hydroxylase (TH)-positive neurons by stereological counting in the substantia nigra (SN) of control and shATP13A2-injected macaques five months after injection (SNpc: *p* = 0.0025, *t* = 3.814). Scale bar: 500 µm (top) 200 µm (inset). **b** Illustrative images (left) and scatter plots (right) of striatal TH, AADC and DAT immunostaining in control and shATP13A2-injected macaques, expressed as a percentage of controls (TH put: *p* = 0.0038, *t* = 3.947; TH cd: *p* = 0.00605, *t* = 3.548; AADC put: *p* = 0.0001, *t* = 10.45; AADC cd: *p* = 0.0001,</figcaption></figure>

*t* = 10.04; DAT put: *p* = 0.0557, *t* = 1.866; DAT cd: *p* = 0.0129, *t* = 2.944). A green fire blue LUT (lookup table) was used to enhance contrast and highlight the differences between non-injected and shATP13A2-injected macaques at the striatum level. Scale bars = 5 mm. Each dot represents one hemisphere of the control (black) and shATP13A2-injected NHPs (red) (light red is one macaque/deep red is the second macaque). The horizontal line indicates the average value per group ± SD. The bootstrapped mean difference with 95% CI (error bar) is shown on the right side of each graph. Comparisons were made using unpaired *t*-tests, \* *p* < 0.05.

amino acids 121–125)⁴². Consequently, this post-translational modification might alter its conformation and render it undetectable by the syn211 antibody while remaining accessible to the anti-p-S129 α-synuclein antibody. This phenomenon could account for the observed increase in one species coinciding with no changes in the other. To get insight on the ratio of soluble versus aggregated α-syn, mesencephalic total protein extracts were subjected to Proteinase-K (PK) digestion and to dot-blotting with α-syn antibody (syn211 antibody) (Fig. 2d). No clear accumulation of α-syn PK-resistant aggregates was observed in the shATP13A2 group compared to control.

To further examine the distribution of α-syn, we performed triple immunofluorescence at the nigral level with the well-established dopaminergic marker TH, the lysosomal marker LAMP1 and the Syn211 antibody (Supplementary Fig. 4). We observed clustered compartments stained for LAMP1 with some clustering with α-syn in dopaminergic neurons of shATP13A2-injected animals compared to controls. Altogether, these results indicate no discernible changes in total α-syn or PK-resistant α-syn levels but increased pSyn protein expression in SN five months after nigral LV-shATP13A2 injection.

## Accumulation of lysosomal-related vesicles in dopaminergic neurons of shATP13A2-injected animals

In vitro studies suggest that loss of ATP13A2 function leads to a rise in the number and size of lysosomes¹³,¹⁴, potentially as compensation for poor lysosomal function. Thus, we evaluated lysosomal-related vesicle number, surface and intracellular distribution in surviving nigral TH-positive neurons in vivo by segmentation analysis of confocal images. Following cell segmentation into whole cell area, perinuclear and cytosolic regions of

interest (ROI) and detection of LAMP2- and LC3-positive vesicles, we calculated the average number and area of LAMP2- and LC3-positive puncta (Fig. 3). Consistent with in vitro reports, we found an increased number of LAMP2-positive lysosomes (Fig. 3a) and an increased number and surface of LC3-positive autophagosomes (Fig. 3b) in surviving nigral TH-positive neurons, consistent with a perturbation of these organelles. Strikingly, we observed an accumulation of both LAMP2- and LC3-positive puncta in the perinuclear area, which might indicate an accretion of immature and undegraded autolysosomes at this location as ATP13A2 inactivation is associated with impairment of autophagosome-lysosome fusion²⁷. These findings are consistent with previous in vitro studies and may reflect compensatory changes in the ALP following improper clearance of lysosomal substrates. Lipofuscin is a pigmented, heterogeneous byproduct of failed intracellular catabolism conventionally found within lysosomes or the cytosol in the brains of patients with ATP13A2 mutations²⁵. We next attempted to examine the presence of lipofuscin in the nigral sections. Compared with control animals, LV-shATP13A2 macaques showed an increased deposition of autofluorescent storage material in the substantia nigra five months after nigral LV-shATP13A2 injection (Supplementary Fig. 5). Overall, these data indicate a dysfunction of the ALP associated with an accumulation of lipofuscin deposits characteristic of NCL.

## ATP13A2 depletion alters metal homeostasis in substantia nigra

Because the loss of ATP13A2 function is associated with dysregulated metal homeostasis in different experimental models²⁶,²⁸,⁴³,⁴⁴ and KRS patients⁴⁵,⁴⁶, we measured nigral heavy metal levels by synchrotron radiation x-ray fluorescence (SR-XRF) (Fig. 4a) as previously described⁴⁷. Iron (Fig. 4b),
<figure><img src="image_6.png" /><figcaption>**Fig. 2** | Decreased ATP13A2 levels induce α-syn pathology in substantia nigra.</figcaption></figure>

a Representative nigral sections of endogenous α-syn immunostaining (top) and scatter plot (bottom) of the mean grey value of syn211 immunostaining in the substantia nigra of control and shATP13A2-injected NHPs by quantification of surface staining, expressed as a percentage of controls macaques (p = 0.16355, t = 1.067). Scale bar: 200 µm. b Illustrative photomicrographs (top) and quantification (bottom) of phosphorylated S129 α-synuclein in the SNpc of control and shATP13A2-injected NHPs (p = 0.0109, t = 3.077). Scale bars = 100 µm (sections) and 10 µm (insets). c Midbrain protein extracts were analysed by dot-blot assay on

nitrocellulose membrane and revealed by immunoblot with syn (p = 0.0528, t = 1.972) and pSyn (p = 0.0001, t = 13.39) antibodies. d Midbrain protein extracts were treated with proteinase K (PK) (0.1 µg/ml) for 0, 15, 30 and 60 min and analyzed by dot-blotting with syn211 antibody. Each dot represents one hemisphere of the control (black) and shATP13A2-injected NHPs (red) (light red is one macaque/deep red is the second macaque). Bootstrapped mean difference with 95% CI (error bar) is shown on the right side of each graph. Comparisons were made using an unpaired t-test. \*p < 0.05 compared to control animals.
<figure><img src="image_7.png" /><figcaption>**Fig. 3** | Accumulation of lysosomal-related structures after reduction of ATP13A2 levels. a Representative nigral sections of endogenous LAMP2 immunostaining in TH-positive cells (left) and scatter plot (right) depicting the average number of LAMP2-positive puncta in the whole cell area ($p = 0.0099, t = 3.004$), the perinuclear ($p = 0.0042, t = 3.631$) and cytosolic ($p = 0.0139, t = 2.768$) area as well as the LAMP2-positive puncta area in the whole cell ($p = 0.2547, t = 0.6951$), the perinuclear ($p = 0.0256, t = 0.6907$) and cytosolic ($p = 0.25605, t = 0.6905$) area. b Representative nigral sections of endogenous LC3 immunostaining in TH-positive cells (left) and scatter plot (right) depicting the average number of LC3-positive</figcaption></figure>

puncta in the whole cell area ($p = 0.02435, t = 2.322$), the perinuclear ($p = 0.0172, t = 2.546$) and cytosolic ($p = 0.0382, t = 2.034$) area as well as the LC3-positive puncta area in the whole cell ($p = 0.01025, t = 2.880$), the perinuclear ($p = 0.01025, t = 2.882$) and cytosolic ($p = 0.00995, t = 2.900$) area. Scale bar: 10 µm. Each dot represents one hemisphere of the control (black) and shATP13A2-injected NHPs (red) (light red is one macaque/deep red is the second macaque). Bootstrapped mean difference with 95% CI (error bar) is shown on the right side of each graph. Comparisons were made using an unpaired *t*-test. \* *$p < 0.05$ compared to control animals.

sulfur (Fig. 4c), and manganese (Fig. 4d) concentrations were significantly increased in LV-shATP13A2 macaques compared to control animals. In the case of copper, concentrations were statistically decreased in LV-shATP13A2 macaques compared to control macaques (Fig. 4e), whereas no statistical differences were observed for calcium (Fig. 4f) and zinc (Fig. 4g). Interestingly, iron concentrations were a three-fold increase in shATP13A2 macaques, reminiscent of human pathology⁴⁵,⁴⁶,⁴⁸⁻⁵⁰. Next, we sought to determine the cell-type specificity and distribution of iron accumulation in the SN. We performed a triple immunofluorescence with the dopaminergic neuronal marker TH, the microglial marker Iba1, and Ferritin, the main iron-storage protein. We detected an increase of ferritin-positive cells in LV-shATP13A2 macaques compared to controls (Supplementary Fig. 6), and we observed increased expression in iron-storage protein ferritin in both Iba1-positive microglia (Supplementary Fig. 6a, d) and TH-positive neurons (Supplementary Fig. 6b, c). Upon closer examination, the iron-positive cells seemed to display distinctive morphological features consistent with activated microglia, including a compact soma and numerous thin processes⁵¹, indicative of an ongoing inflammatory response. This immunofluorescence analysis suggests that iron accumulates in dopaminergic neurons and microglial cells, with more abundant smaller puncta distributed in the microglia (Supplementary Fig. 6c, d).

## Discussion

Decreased ATP13A2 expression in the substantia nigra through an LV-shRNA strategy leads to major pathological features in nonhuman primates reminiscent of PD pathology. Knockdown of ATP13A2 by ~60% in nigral dopaminergic neurons affected their integrity, as evidenced here by (i) changes of α-synuclein immunoreactivity associated with an increase of one of its pathological forms (i.e., S129 phosphorylation); (ii) a dyshomeostasis in heavy metals and (iii) occurrence of ALP alterations. All these features were prominent in the SN of shATP13A2-injected NHP, resulting in a loss of striatal dopaminergic terminals and SNpc dopaminergic neurons. Relevant to KRS, we also show that silencing of ATP13A2 induces iron accumulation in the substantia nigra. These observations indicate an instrumental role for ATP13A2 deficiency in lysosomal function, α-synuclein and heavy metal homeostasis, and eventually cell viability. This new PD model provides a mechanistic framework that recapitulates some key pathological features of ATP13A2-related diseases (Fig. 5).

However, it is essential to acknowledge the limitations inherent in this study, foremost among them being the relatively small cohort sizes, typically consisting of only 2 to 3 animals per group. Nevertheless, this size is commonly recognised as a standard size for pilot studies in non-human primate (NHP) research, grounded in scientific considerations
<figure><img src="image_8.png" /><figcaption>**Fig. 4** | Decreased ATP13A2 levels alter SN heavy metal concentrations.
a Representative elemental mapping in nigral slices of the macaque groups. Levels of iron (b), sulfur (c), manganese (d), copper (e), calcium (f), and zinc (g) were measured using synchrotron X-ray fluorescence in the SN of control and shATP13A2-injected macaques (Iron: *p* = 0.0022, *t* = 4.908; Zinc: *p* = 0.1292, *t* = 1.275; Copper: *p* = 0.0028, *t* = 4.644; Calcium: *p* = 0.15715, *t* = 1.118; Manganese: *p* = 0.0011,</figcaption></figure>

*t* = 5.785; Sulfur: *p* = 0.0038, *t* = 4.322). Each dot represents one hemisphere of the control (black) and shATP13A2-injected NHPs (red) (light red is one macaque/deep red is the second macaque). Bootstrapped mean difference with 95% CI (error bar) is shown on the right side of each graph. Comparisons were made using an unpaired Student's *t*-test. \* *p* < 0.05 compared to control animals.

and ethical constraints³². Notably, our study's statistical power is diminished, a common characteristic in NHP research, where a significant proportion often involves minimal animal numbers per group, typically ranging from 2 to 4 individuals, due to cost and difficulty accessing macaques. With two groups of two to three macaques, either LV-shATP13A2-injected or non-injected historical controls, it is challenging to obtain for each variable homogenous results that allow giving definitive answers. Some endpoints do not reach statistical significance because of the small sample size. Hence, the importance of the Gardner–Altman plots that focus upon the effect size of the difference, when any. Many studies conducted by our team and other colleagues working on NHPs have generally used these types of low group sizes³³⁻³⁶. Our team is, however, known to involve a much larger number of NHPs when moving to definitive studies¹⁷,⁶⁴⁻⁶⁷ that have always confirmed experiments with lower n.

To investigate the pathological consequences of ATP13A2 deficiency, a vast panel of models of ascending complexity, i.e., in human dopaminergic neuroblastoma-derived cells, mesencephalic primary cultures, mice, and dogs, have been developed. However, in ATP13A2 knockout mice, behavioural phenotype impairment and pathophysiological abnormalities are not detectable until old age (20- to 29-month-old), without overt neurodegeneration³⁰,³¹, suggesting the involvement of strong compensatory mechanisms in transgenic mice. Dogs with ATP13A2 mutations show clinical signs of NCL similar to those of KRS patients, but without parkinsonism, associated with astrogliosis and autofluorescent storage material in nervous and cardiac tissues at an average age of 6 years, progressively worsening over time³²,⁶⁸,⁶⁹. A striking result of the current study is the efficacy of the induction of dopaminergic neurodegeneration over five months, corresponding to a short duration for a study in nonhuman primates.
<figure><img src="image_9.png" /><figcaption>**Fig. 5** | Direct comparison between experimental groups to highlight the pathological profile following the decrease ATP13A2 expression in the substantia nigra via a lentiviral shRNA approach. Heat map representing the normalized expression (Z scoring from control values) on the different variables used in this study, for each hemisphere of the two experimental NHPs. Controls are pulled as mean. From left to right: stereological counting Tyrosine Hydroxylase (TH)-positive neurons in the substantia nigra (SN) (stereo TH SN), TH, AADC and DAT immunostaining levels measured in the putamen and in the caudate nucleus (TH putamen; TH caudate; AADC putamen; AADC caudate; DAT putamen; DAT caudate), α-syn and pSyn signals in the SN (syn211 SN, pSyn SN, DB syn211 SN, DB pSyn SN, ratio pSyn/Syn211), the average number of LAMP2-positive puncta and LC3-positive puncta in the whole cell, perinuclear and cytosolic areas (LAMP2 nb whole cell, LAMP2 nb perinuclear, LAMP2 nb cytosolic, LC3 nb whole cell, LC3 nb perinuclear, LC3 nb cytosolic), the heavy metals content in SN (Iron, Zinc, Copper, Calcium, Manganese, Sulfur). The colour bars represent the z-score value of the ratio of each brain region.</figcaption></figure>

We hypothesise that follow-up long-term studies are required to observe the occurrence of a full spectrum of behavioural deficits, including progressive motor and cognitive impairment, i.e., when the threshold of the nigrostriatal lesion for manifest PD will be reached³⁹. Additionally, we should also target other brain regions to extend the disease modelling panel associated with ATP13A2 loss of function. A second aspect to be envisioned for future studies is to compare the effects obtained with an adeno-associated virus (AAV)-based strategy. Several advantages can be emphasised by moving to AAV due to its unique biological and biophysical properties. First, AAV is compatible with the viral capacity to carry recombinant DNA such as

shRNA. Second, AAV has a diffusion advantage over LV due to its particle size, which is 25 nm in diameter, whereas LV particles are no less than 100 nm in diameter. Third, specific serotypes of AAV share the ability to cross the blood-brain barrier. They may be delivered through various administration routes, making them a powerful tool for genetic material delivery to the CNS⁷⁰⁻⁷². Fourth, AAV has proved safe and effective in preclinical and clinical settings and has been the virus of choice for developing gene therapies for PD⁷³.

While this proof-of-concept study is encouraging, pending questions await future investigation and should be conducted with a larger sample size and a parallel control group to increase statistical power. To date, it is evident that numerous neurological disorders have been linked to ATP13A2 with diverse and complex neurological deficits, complicating our understanding of ATP13A2. However, the fact that mutations in a single lysosomal gene result in such a diversity of disease phenotypes suggests that lysosomal burden is a common thread running amongst them. Therefore, deciphering the ATP13A2-dependent cellular roles and pathways would enhance our understanding of brain physiology and the broader mechanisms underpinning neurodegeneration. First, one aspect to investigate in future studies will be the analysis of polyamine levels. Recent studies demonstrated the polyamine transport function of ATP13A2¹⁶, which may support heavy metal chelation⁷⁴ and provide antioxidative effects⁷⁵. Polyamine levels decline with ageing and are altered in PD patients. Analysing polyamine levels in this new experimental animal model should provide valuable information regarding the precise role of polyamines in PD, KRS, and beyond. Second, in humans, *ATP13A2* mutations cause neurodegeneration with brain iron accumulation¹⁴; in particular, KRS patients present iron deposits in the brain⁴⁵,⁴⁶ whereas, in mammalian cells, ATP13A2 expression and transport activity protect against several heavy metals, such as manganese, zinc, and iron toxicity⁴³,⁷⁶. Our study measured an iron increase of over 300%, demonstrating a close relationship between iron accumulation and ATP13A2 depletion. Several experiments proposed that the ferric iron load accelerates α-synuclein aggregation⁷⁷, and an increase in iron levels plus a decrease in ferritin can contribute to nigral degeneration⁷⁸. An increased iron burden has been found in early PD stages in the SN⁴⁸⁻⁵⁰,⁷⁸. The pathological consequences of the increased iron SN burden could contribute to α-syn pathology, as we observed an increased nigral pSyn expression.

Noteworthy, we also observed a significant increase in manganese in the SN, consistent with the fact that ATP13A2 has been reported to be involved in manganese homeostasis²⁶,²⁸,⁴⁴. This result is reminiscent of previous reports showing that ATP13A2 protects cells against manganese toxicity by transporting manganese into lysosomes; however, this protection is lost in the case of mutations in this gene⁴³. A recent study showed that ATP13A2 polymorphism could be a risk marker for the neurotoxic effects of manganese in humans⁷⁹. Both of our findings are of interest regarding the work of Fleming and colleagues that show increased manganese and iron in ATP13A2 knockout mice administered with manganese chloride compared to treated wild-type mice but also increased α-synuclein aggregation and higher lipofuscin deposits in the SN⁸⁰. Of note, the copper reduction is also of interest and deserves further studies as it has been reported to be reduced in human SN dopaminergic neurons of PD and incidental Lewy body disease cases⁸¹.

Overall, we show that ATP13A2 depletion triggers PD/KRS-related features in nonhuman primates, with stimulating perspectives for understanding human ATP13A2 pathology per se. Further studies are needed to extend our observations and appreciate the full spectrum of behavioural and pathological outcomes of ATP13A2 deficiency. Understanding the complex contribution of the dysfunctional or absent lysosomal ATP13A2 protein to ALP dysfunction, heavy metal dyshomeostasis, and neurodegeneration would likely open new therapeutic opportunities for slowing down the degenerative process in PD/KRS patients.
# Methods

## Ethics statement

All experimental protocols comply with the Council Directive of 2010 (2010/63/EU) of the European Community and the National Institute of Health Guide for the Care and Use of Laboratory Animals. The proposed research has received approbation from the French Ethical Committee for Animal Research CE50 (agreement number: 50120102-s).

## Recombinant LV-shRNA plasmid and virus production

shRNA lentiviral plasmids were purchased from Sigma-Aldrich. shRNAs were validated based on the Genbank human ATP13A2 sequence NM_022089. All sequences were BLAST confirmed for specificity. Synthetic DNA encoding shRNA sequences targeting human ATP13A2 were cloned into shRNA lentiviral plasmids (pLK0.1-puro) under the control of the U6 promoter, which is a type III RNA polymerase III promoter commonly used for driving small hairpin RNA (shRNA) expression in vector-based RNAi. shRNA lentiviral plasmids were transfected into BE(2)-M17 cells for 48 h to determine the level of knockdown of gene expression by immunoblotting¹³. The ATP13A2 shRNA-expressing viruses have been tested as described above. The virus titers were identical and corresponds to 6.73 × 10⁸ pI/ml for ATP13A2 shRNA and Scramble shRNA sequence. The ATP13A2 shRNA sequence used in the study is ccggGCCCCATCAACTTCAAGTTC-TACTCGAGTAGAACTTGAAGTTGATGGGCttttt and the Scramble shRNA sequence used in the study is ccggCAACAAGATGAA-GAGCACCAACTCGAGTTGGTGCTCTTCATCTTGTTGttttt. Lentiviral vector production was done by the service platform for lentiviral vector production "Vect'UB' of the TMB-Core of Bordeaux University. A lentiviral vector was produced by transient transfection of 293 T cells according to standard protocols. In brief, subconfluent 293 T cells were co-transfected with the lentiviral genome (psPAX2)⁸, with an envelope coding plasmid (pMD2G-VSVG) and with vector constructs by calcium phosphate precipitation. LVs were harvested 48 h post-transfection and concentrated by ultracentrifugation. Viral titers of pLV lentivectors were determined by transducing 293 T cells with serial dilutions of viral supernatant, and a control pLV-EGFP expression was quantified 5 days later by flow cytometry analysis or provirus copies number by qPCR method.

## Animals and stereotactic injections

Two female adult rhesus macaques (Macaca mulatta) weighing ~6-7 kg (9-10 years old) were used in this study. The macaques were housed in the animal facility of the Institute of Neurodegenerative Diseases (UMR CNRS 5293) under standard conditions (12/12 h day/night cycle, ~60% humidity, and ~22 °C). Tissues from 4 to 6 additional age- and sex-matched historical control animals were used for post-mortem experiments when needed and available³,⁵,³⁶. Bilateral intranigral LV injections of 10 µl were performed with electrophysiological guidance considering the line passing through the anterior (AC)-posterior (PC) commissures as a reference: -7mm posterior to AC, -4 mm below the AC-PC line, and 3 mm lateral from midline. At the end of the experiment (5 months post-injection), as previously described⁴⁷,⁶²,⁶⁶, all macaques were euthanised by sodium pentobarbital overdose (150 mg/kg iv), followed by perfusion with saline solution (containing 1% heparin) and 4% paraformaldehyde, and brains were quickly removed and processed for histological studies. Brains were sectioned in 50 µm thick serial free-floating coronal sections in a Leica CM3050S cryostat (Leica Microsystems, Wetzlar, Germany) at -20 °C, collected in PBS Azide 0.2% and stored at 4 °C until they were processed for analysis.

## Histological analysis

**Dopaminergic neurodegeneration assessments.** To assess the effect of ATP13A2 silencing on dopaminergic neurones and fibres, TH (tyrosine hydroxylase), AADC (Aromatic L-amino acid decarboxylase), and DAT (Dopamine Transporter) immunohistochemistry and quantifications were performed on striatal and SN sections as previously described⁴⁷,⁶²,⁶⁶. Briefly, 50-µm free-floating sections from one representative level of the striatum (post-anterior commissure) and serial

sections (1 every 12) corresponding to the whole SNpc were selected. For TH staining, anterior striatum (post-anterior commissure) sections and a whole rostrocaudal series (1 to 12) encompassing the SNpc were selected. For AADC and DAT staining, anterior striatum (post-anterior commissure) sections were used. Sections were incubated with mouse monoclonal TH antibody (Millipore, MAB318, 1:5000), rabbit polyclonal AADC antibody (Merck AB136, 1:1000), or rat monoclonal DAT antibody (Merck MAB369, 1:500) for one night at room temperature (RT) and revealed the next day with the corresponding peroxidase EnVision secondary antibody, followed by DAB visualisation. SNpc sections were mounted on gelatine-coated slides, counterstained with 0.1% cresyl violet solution, dehydrated and cover-slipped. Striatal sections were mounted on gelatine-coated slides, dehydrated and cover-slipped. Striatal sections were analysed by optical density (OD) in the caudate nucleus and putamen. The slides were scanned using Epson Expression 10000XL high-resolution scanner. Images were analysed using ImageJ open-source software (version 1.53) to compare mean grey levels in the caudate nucleus and putamen. TH-positive neurons of the SNpc were counted by stereology blind about the experimental condition using a Leica DM6000B microscope coupled with the Mercator software (IMASCOPE, France). The SNpc was delineated for each slide, and disector probes for stereological counting (100 × 80 µm spaced by 600 × 400 µm) were applied to the map obtained. Each TH-positive cell with the nucleus included in the probe was counted following the stereological rules of exclusion⁴⁷,⁶²,⁶⁶. The total number of TH-positive neurons in the whole SNpc was then assessed per hemisphere using the optical fractionator method.

**α-Synuclein pathology assessment.** α-synuclein status was assessed with a mouse monoclonal antibody raised against human α-syn (ThermoFisher [Syn211], 32-8100, 1:1000) and another against phosphorylated α-syn (Abcam, [EP1536Y], ab51253, 1:5000) as previously described⁴⁷,⁶²,⁶⁶. Briefly, selected sections of all animals were explicitly identified and incubated in the same well to allow direct comparison of immunostaining intensity. Sections were incubated overnight at room temperature with the antibodies. The following day, the revelation was performed with the corresponding peroxidase EnVision secondary antibody followed by 3,3'-diaminobenzidine (DAB, DAKO, K346811-2) incubation. Sections were then mounted on gelatinized slides, counterstained with 0.1% cresyl violet solution, if necessary, dehydrated, and cover-slipped until further analysis. Slides were scanned with a high-resolution scanner (Panoramic Scan II, 3D HISTECH Ltd, France) at x20 magnification and on five layers spaced by 1.4 µm each. Quantifications were estimated by immunostaining-positive surface quantification at regional levels with the Mercator software (IMASCOPE, France), as previously described⁴⁷,⁶²,⁶⁶.

The “surface” is an additional quantification method used in tissue sections. For these analyses, a specific staining process was used to keep all tissues together in the same solution during the staining process. Then, high-resolution whole colour slide images were acquired with the 3D Histech Panoramic Scanner at 20X magnification, with 5 layers in extended mode. Each image was opened in the off-line MERCATOR PRO 7.12.3 software (Explora Nova, France), and all regions of interest were mapped. The same brightness and contrast rules were systemically applied to the RGB pictures to optimise the staining’s details without any image saturation. The colour thresholding tool was then used to select the threshold corresponding to the brown colour revealed by the DAB staining. The threshold was established based on the intensity of the staining to detect the maximum DAB staining. The file of the threshold parameters was saved and applied to all measurements for each animal/staining. Before performing the quantification, the threshold was randomly applied to some images of different treatment groups to verify the accuracy of the settings. In each region, the software extracted the surface corresponding to the threshold defined. The “surface” parameter was finally expressed as a ratio of the total surface of each area of interest. This parameter has been previously used and published⁸³.
**Immunofluorescence and image analysis.** To assess ATP13A2 expression and autophagy-related markers (i.e., autophagosomes and lysosomes), nigral macaque sections were incubated simultaneously with two antibodies: a mouse monoclonal antibody against TH (Millipore, MAB318, 1:5000), and, either a rabbit monoclonal antibody against ATP13A2 (Novus Biological, NB110-41486, 1:250) or a rabbit monoclonal against human LAMP1 (Genetex, GTX62434, 1:1000) or LAMP2 (Santa Cruz Biotechnology, sc-18822, 1:1000) or LC3 (Novus Biological, NB100-2220, 1:1000) or human α-syn (ThermoFisher [Syn211], 32-8100, 1:1000) for one night at room temperature. Incubation with secondary antibodies was done sequentially with a goat anti-mouse Alexa-Fluor 488 (Invitrogen, A28751, 1:400) or a goat anti-rabbit AlexaFluor 568 (Invitrogen, A11036, 1:400) 1h30 at room temperature. Sections were mounted on non-gelatinised slides, and cover slipped using fluorescent mounting media without DAPI (Vector Labs). Images were acquired using a Zeiss SP5 confocal microscope.

For LAMP2 and LC3 count and area analysis in TH-positive cells at 5 months post-injection, image analysis was performed in Fiji/ImageJ. All image acquisitions and analyses were performed blinded to the researcher. Image analysis was performed in Fiji/ImageJ using custom scripts (available at https://github.com/Soria FN/Lysosome_analysis). Cells were segmented by manual cell selection using the TH channel. The mask was then applied to the LAMP2 and LC3 channels, where LAMP2 and LC3 staining was quantified by automatic thresholding followed by binarisation. Total staining was then normalised per cell. For subcellular localisation analysis, cytoplasmic and nuclear ROIs were manually segmented using the TH channel of single z-planes where the nucleus was clearly visible. To prevent a bias due to low cytoplasmic area and to ensure that the comparisons were made between cells with similar cytoplasmic ROIs, the analysis was subsequently performed only in cells where the ratio of cytoplasmic to nuclear ROI area was 2:1. Lysosomal-related puncta were segmented from the LAMP2 or LC3 channel, and the average number and area from the nucleus centre to the centre of mass of each lysosomal ROI was calculated. Finally, to estimate the distance to the nucleus border, the average radius of an ellipse fitted to the nuclear ROI was subtracted from the distance to the nucleus centre.

To assess ferritin pattern expression, each nigral macaque section punches was made from free-floating sections and mounted on the same slide to compare the pattern between TH, Iba1 and ferritin-expressing cells in SNpc. The resulting glass slide was incubated simultaneously with a mix of antibodies, including chicken anti-TH (Abcam, ref ab76442), goat anti-Iba1 (Abcam, ref ab5076) and rabbit anti-ferritin [EPR3004Y] (Abcam, ref ab75973), all diluted at 1:1000 and incubated for one night at room temperature. Incubation with secondary antibodies was done after thorough rinses, with a mix of donkey fluorescent secondary antibodies raised against the different species (respectively anti-chicken AlexaFluor 568 (Invitrogen, A11041), anti-goat AlexaFluor 647 (Invitrogen, A21447) and anti-rabbit AlexaFluor488 (Invitrogen, A21206) diluted at 1:500 and incubated 1 h at room temperature. After several rinses with PBS, a nuclear stain was made with 1:5000 diluted Hoechst solution for 30 s, and sections were then covered slipped using ProlongGlass fluorescent mounting media (Thermo Fischer, ref P36984). Images were acquired using a 3D-HISTECH Panoramic Scan II with 20X objective and extended focus mode (5 layers and 7 steps).

**Lipofuscin imaging.** For imaging of lipofuscin, untreated free-floating nigral slices were mounted with a Vectashield Antifade Mounting Medium to observe lipofuscin autofluorescence. Z-stack images were acquired using a Leica DM6 CFS TCS SP8 confocal microscope at x63 magnification, with an excitation of 500 nm and an emission of 600-650 nm. Analysis was performed through ImageJ Fiji.

**Biochemical analysis**

**Protein Extraction and immunoblot analysis.** We used 2 mm³ formalin-fixed tissue per animal. Tissue patches were extracted using the

Qproteome FFPE TissueKit (Qiagen)⁸⁴⁻⁸⁶ and then quantified using the Lowry method (Biorad RC DC Protein Assay Kit) following the manufacturer's instructions. Based on total protein concentrations calculated from the assay, aliquots of tissue lysates corresponding to 20 µg were prepared for each structure in Laemmli buffer (Tris-HCl 25 mM pH 6.8, glycerol 7.5%, SDS 1%, DTT 250 mM and bromophenol blue 0.05%).

For immunoblot analysis, samples were then boiled at 100 °C for 5 min and loaded into a 12% SDS-PAGE gel, followed by electrophoretic transfer onto 0.2 µm nitrocellulose membranes (Bio-Rad, Hercules, CA, United States). Membranes were blocked in phosphate buffer saline (PBS) with 5% dry skimmed milk and probed with an antibody against ATP13A2 (Novus Bio, NB110-41486, 1:1000) overnight. Anti-β-actin (1:10 000, Sigma, A5441) was used to control equal loading. Appropriate secondary antibodies coupled to peroxidase were revealed using a Super Signal West Pico Chemiluminescent kit (Immobilon Western, Chemiluminescent HRP substrate, Millipore) and analysed with Image Lab software 5.0 (Bio-Rad, Hercules, United States).

**Dot blotting analysis.** The dot blot analysis was performed as previously described⁸⁷. After heating at 100 °C for 5 min, 20 µg of protein extract was diluted in buffer (25 mM Tris-HCl, 200 mM Glycine, 1% SDS) and filtered through a nitrocellulose membrane (Bio-Rad, 0.2µm pore size). Membranes were then saturated in 5% dry milk before incubation with total α-syn (1:1000, ThermoFisher [Syn211], 32-8100) and anti-phosphorylated-α-syn at Ser129 (1:5000, Abcam [EP1536Y], ab51253). To evaluate proteinase K-resistant α-syn contained in midbrain protein extracts from control and shATP13A2-injected animals, we subjected 50 µg of protein extract to digestion with proteinase K (0.1 µg/ml) for 0, 15, 30, and 60 min. The reaction was stopped by boiling for 5 min at 100 °C before dot blotting with syn211 antibody. Revelation was done as described above.

**SR-XRF microscopy elemental mapping of brain tissue cryosections**

The synchrotron experiments were carried out at the Diamond Light Source, Harwell Science and Innovation Campus (Didcot, UK) with 3-GeV energy of the storage ring and 300-mA currents with top-up injection mode. All SR-XRF microscopy investigations reported herein were carried out on the microfocus spectroscopy beamline (I18). The micro-XRF elemental mapping was acquired at room temperature with an incident x-ray energy set to 12 keV using a Si(111) monochromator. This resulted in an x-ray photon flux of $2 \times 10^{11}$ photons/s. The SN cryostat section (50 µm in thickness) of each animal was collected from free-floating sections and mounted onto an x-ray transparent metal-free 4-µm-thick Ultralene foil (SPEX CertiPrep, Metuchen, NJ, USA) secured to a customised polyetheretherketone (PEEK) holder ensuring contamination-free samples and reduced x-ray scattering contribution. The samples were then left and kept at room temperature in a moisture-free box until used for analysis. The samples were affixed to a magnetic plate connected to the sample stage of the microfocus beamline. The four-element Si drift Vortex ME4 energy-dispersive detector (Hitachi High-Technologies Science America), with Xspress 3 processing electronics, was operated in the 90° geometry; hence, it minimizes the background signal. The sample-detector distance was fixed (75 mm). The sample was held at 45° to the incident X-ray beam and raster in front of the beam while the XRF spectra were collected. An area of 500 µm by 500 µm within the SNpc was mapped for each sample with a step size that matches the beam size (5 µm) and a dwell time of 1 s per pixel due to the low concentration of the element. A thin (100 µm) pellet of the National Institute of Standards and Technology (NIST) standard reference materials SRM1577c (bovine liver material, NIST, Gaithersburg, MD, USA) was measured to calibrate experimental parameters as well as a thin-film XRF reference material (AXO DRESDEN GmbH) with nanoscale uniform mass depositions in the range of ng/mm². This was followed by elemental quantification calculated using the Fundamental Parameters (FP) approach through the open-source software PyMCA⁸⁸, in which the reference
material and the sample are modelled in terms of main composition, density, and thickness. The fluorescence spectrum obtained from each pixel was fitted, the elemental concentration (micrograms per gram of dry weight or parts per million) maps were generated, and an average elemental concentration of the SNpc regions was obtained.

## Statistical analysis

For all experiments, comparisons among means were performed using raw data. Student's two-tailed unpaired *t*-test (GraphPad Prism 10.0, San Diego, CA). Statistical significance was set at *p* < 0.05. The debate about the need to move beyond *p*-value is raging. Data must be further analysed with estimation graphics<sup>89</sup> emphasising the effect size. Therefore, all data appear as estimation graphics called 'Gardner-Altman plots': on the left of each graph, data of controls and sh-ATP13A2 groups are presented as scatter plots showing the observed values along with the defined descriptive statistics (mean ± standard deviation). Each dot represents one hemisphere of the control (shade of black) and shATP13A2-injected NHPs (shade of red). One circle corresponds to one hemisphere. The dark circles correspond to the control group, and the red circles correspond to the shATP13A2-injected groups (Light red is one macaque/deep red is the second macaque). On the right of each graph is a contrast graph using the different axes to display an effect size, which means difference. Horizontally aligned with the mean of the test group, the mean difference is indicated by the black circle. The black vertical line illustrates the 95% confidence interval (CI) of the mean difference. Given the observed data, the curve shows the resampled distribution of the effect size.

## Data availability

The data supporting the findings of this study are provided in Supplementary Table 1. The material supporting the findings of this study is available from the corresponding author on request.

Received: 8 January 2024; Accepted: 24 July 2024;

Published online: 01 August 2024

## References

1. Ramirez, A. et al. Hereditary parkinsonism with dementia is caused by mutations in ATP13A2, encoding a lysosomal type 5 P-type ATPase. *Nat. Genet* **38**, 1184-1191 (2006).

2. Chien, H. F. et al. Neuropathologic findings in a patient with juvenile-onset levodopa-responsive parkinsonism due to ATP13A2 mutation. *Neurology* **97**, 763-766 (2021).

3. Kruer, M. C. & Boddaert, N. Neurodegeneration with brain iron accumulation: a diagnostic algorithm. *Semin Pediatr. Neurol.* **19**, 67-74 (2012).

4. Yang, X. & Xu, Y. Mutations in the ATP13A2 gene and Parkinsonism: a preliminary review. *Biomed. Res. Int.* **2014**, 371256 (2014).

5. Park, J. S., Blair, N. F. & Sue, C. M. The role of ATP13A2 in Parkinson's disease: clinical phenotypes and molecular mechanisms. *Mov. Disord.* **30**, 770-779 (2015).

6. Gagliardi, M. et al. Mutation analysis of the ATP13A2 gene in patients with PD and MSA from Italy. *J. Neurol. Sci.* **430**, 120031 (2021).

7. Estrada-Cuzcano, A. et al. Loss-of-function mutations in the ATP13A2/PARK9 gene cause complicated hereditary spastic paraplegia (SPG78). *Brain* **140**, 287-305 (2017).

8. Estiar, M. A. et al. Clinical and genetic analysis of ATP13A2 in hereditary spastic paraplegia expands the phenotype. *Mol. Genet Genom. Med.* **8**, e1052 (2020).

9. Algahtani, H. et al. Autosomal recessive spastic paraplegia type 78 associated with a homozygous variant in the ATP13A2 gene. *Mov. Disord. Clin. Pract.* **9**, 997-1002 (2022).

10. Zhang, F., Liu, P., Li, J., Cen, Z. & Luo, W. A novel ATP13A2 variant causing complicated hereditary spastic paraplegia. *Neurol. Sci.* **45**, 1749-1753 (2024).

11. Spataro, R. et al. Mutations in ATP13A2 (PARK9) are associated with an amyotrophic lateral sclerosis-like phenotype, implicating this locus in further phenotypic expansion. *Hum. Genom.* **13**, 19 (2019).

12. Croucher, K. M. & Fleming, S. M. ATP13A2 (PARK9) and basal ganglia function. *Front. Neurol.* **14**, 1252400 (2023).

13. Dehay, B. et al. Loss of P-type ATPase ATP13A2/PARK9 function induces general lysosomal deficiency and leads to Parkinson disease neurodegeneration. *Proc. Natl. Acad. Sci. USA* **109**, 9611-9616 (2012).

14. Usenovic, M., Tresse, E., Mazzulli, J. R., Taylor, J. P. & Krainc, D. Deficiency of ATP13A2 leads to lysosomal dysfunction, alpha-synuclein accumulation, and neurotoxicity. *J. Neurosci.* **32**, 4240-4246 (2012).

15. Bento, C. F., Ashkenazi, A., Jimenez-Sanchez, M. & Rubinsztein, D. C. The Parkinson's disease-associated genes ATP13A2 and SYT11 regulate autophagy via a common pathway. *Nat. Commun.* **7**, 11803 (2016).

16. van Veen, S. et al. ATP13A2 deficiency disrupts lysosomal polyamine export. *Nature* **578**, 419-424 (2020).

17. Murphy, K. E., Cottle, L., Gysbers, A. M., Cooper, A. A. & Halliday, G. M. ATP13A2 (PARK9) protein levels are reduced in brain tissue of cases with Lewy bodies. *Acta Neuropathol. Commun.* **1**, 11 (2013).

18. Fujii, T. et al. Parkinson's disease-associated ATP13A2/PARK9 functions as a lysosomal H(+),K(+)-ATPase. *Nat. Commun.* **14**, 2174 (2023).

19. van Veen, S. et al. Cellular function and pathological role of ATP13A2 and related P-type transport ATPases in Parkinson's disease and other neurological disorders. *Front Mol. Neurosci.* **7**, 48 (2014).

20. Chen, X. et al. Cryo-EM structures and transport mechanism of human P5B type ATPase ATP13A2. *Cell Discov.* **7**, 106 (2021).

21. Tomita, A. et al. Cryo-EM reveals mechanistic insights into lipid-facilitated polyamine export by human ATP13A2. *Mol. Cell* **81**, 4799-4809.e4795 (2021).

22. Tillinghast, J., Drury, S., Bowser, D., Benn, A. & Lee, K. P. K. Structural mechanisms for gating and ion selectivity of the human polyamine transporter ATP13A2. *Mol. Cell* **81**, 4650-4662.e4654 (2021).

23. Sim, S. I., von Bulow, S., Hummer, G. & Park, E. Structural basis of polyamine transport by human ATP13A2 (PARK9). *Mol. Cell* **81**, 4635-4649.e4638 (2021).

24. Mu, J. et al. Conformational cycle of human polyamine transporter ATP13A2. *Nat. Commun.* **14**, 1978 (2023).

25. Malandrini, A., Rubegni, A., Battisti, C., Berti, G. & Federico, A. Electron-dense lamellated inclusions in 2 siblings with Kufor-Rakeb syndrome. *Mov. Disord.* **28**, 1751-1752 (2013).

26. Schmidt, K., Wolfe, D. M., Stiller, B. & Pearce, D. A. Cd2+, Mn2+, Ni2+ and Se2+ toxicity to Saccharomyces cerevisiae lacking YPK9p the orthologue of human ATP13A2. *Biochem. Biophys. Res. Commun.* **383**, 198-202 (2009).

27. Wang, R. et al. ATP13A2 facilitates HDAC6 recruitment to lysosome to promote autophagosome-lysosome fusion. *J. Cell Biol.* **218**, 267-284 (2019).

28. Anand, N. et al. Dysregulated iron metabolism in C. elegans catp-6/ATP13A2 mutant impairs mitochondrial function. *Neurobiol. Dis.* **139**, 104786 (2020).

29. Heins-Marroquin, U., Jung, P. P., Cordero-Maldonado, M. L., Crawford, A. D. & Linster, C. L. Phenotypic assays in yeast and zebrafish reveal drugs that rescue ATP13A2 deficiency. *Brain Commun.* **1**, fcz019 (2019).

30. Schultheis, P. J. et al. Atp13a2-deficient mice exhibit neuronal ceroid lipofuscinosis, limited alpha-synuclein accumulation and age-dependent sensorimotor deficits. *Hum. Mol. Genet* **22**, 2067-2082 (2013).

31. Kett, L. R. et al. alpha-Synuclein-independent histopathological and motor deficits in mice lacking the endolysosomal Parkinsonism protein Atp13a2. *J. Neurosci.* **35**, 5724-5742 (2015).
32. Farias, F. H. et al. A truncating mutation in ATP13A2 is responsible for adult-onset neuronal ceroid lipofuscinosis in Tibetan terriers. *Neurobiol. Dis.* **42**, 468–474 (2011).

33. Van den Haute, C., Eggermont, K., Nuttin, B., Debyser, Z. & Baekelandt, V. Lentiviral vector-mediated delivery of short hairpin RNA results in persistent knockdown of gene expression in mouse brain. *Hum. Gene Ther.* **14**, 1799–1807 (2003).

34. Barnes, L., Widdowson, P., Bezard, E., Kelleher, M. & Mitrophanous, K. Optimising lentiviral vector delivery to the brain using infusion techniques. *Hum. Gene Ther.* **22**, A80–A81 (2011).

35. Deffains, M. et al. L-DOPA regulates alpha-synuclein accumulation in experimental parkinsonism. *Neuropathol. Appl. Neurobiol.* **47**, 532–543 (2021).

36. Recasens, A. et al. Lewy body extracts from Parkinson disease brains trigger alpha-synuclein pathology and neurodegeneration in mice and monkeys. *Ann. Neurol.* **75**, 351–362 (2014).

37. Anglade, P., Vyas, S., Hirsch, E. C. & Agid, Y. Apoptosis in dopaminergic neurons of the human substantia nigra during normal aging. *Histol. Histopathol.* **12**, 603–610 (1997).

38. Jackson-Lewis, V., Jakowec, M., Burke, R. E. & Przedborski, S. Time course and morphology of dopaminergic neuronal death caused by the neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine. *Neurodegeneration* **4**, 257–269 (1995).

39. Bezard, E. et al. Relationship between the appearance of symptoms and the level of nigrostriatal degeneration in a progressive 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine-lesioned macaque model of Parkinson's disease. *J. Neurosci.* **21**, 6853–6861 (2001).

40. Anderson, J. P. et al. Phosphorylation of Ser-129 is the dominant pathological modification of alpha-synuclein in familial and sporadic Lewy body disease. *J. Biol. Chem.* **281**, 29739–29752 (2006).

41. Ghanem, S. S. et al. alpha-Synuclein phosphorylation at serine 129 occurs after initial protein deposition and inhibits seeded fibril formation and toxicity. *Proc. Natl. Acad. Sci. USA* **119**, e2109617119 (2022).

42. Giasson, B. I. et al. A panel of epitope-specific antibodies detects protein domains distributed throughout human alpha-synuclein in Lewy bodies of Parkinson's disease. *J. Neurosci.* **Res. 59**, 528–533 (2000).

43. Gitler, A. D. et al. Alpha-synuclein is part of a diverse and highly conserved interaction network that includes PARK9 and manganese toxicity. *Nat. Genet* **41**, 308–315 (2009).

44. Tan, J. et al. Regulation of intracellular manganese homeostasis by Kufor-Rakeb syndrome-associated ATP13A2 protein. *J. Biol. Chem.* **286**, 29654–29662 (2011).

45. Schneider, S. A. et al. ATP13A2 mutations (PARK9) cause neurodegeneration with brain iron accumulation. *Mov. Disord.* **25**, 979–984 (2010).

46. Bruggemann, N. et al. Recessively inherited parkinsonism: effect of ATP13A2 mutations on the clinical and neuroimaging phenotype. *Arch. Neurol.* **67**, 1357–1363 (2010).

47. Bourdenx, M. et al. Identification of distinct pathological signatures induced by patient-derived alpha-synuclein structures in nonhuman primates. *Sci. Adv.* **6**, eaaz9165 (2020).

48. Dexter, D. T. et al. Increased nigral iron content in postmortem parkinsonian brain. *Lancet* **2**, 1219–1220 (1987).

49. Groger, A. & Berg, D. Does structural neuroimaging reveal a disturbance of iron metabolism in Parkinson's disease? Implications from MRI and TCS studies. *J. Neural Transm. (Vienna)* **119**, 1523–1528 (2012).

50. Wallis, L. I. et al. MRI assessment of basal ganglia iron deposition in Parkinson's disease. *J. Magn. Reson Imaging* **28**, 1061–1067 (2008).

51. Paolicelli, R. C. et al. Microglia states and nomenclature: a field at its crossroads. *Neuron* **110**, 3458–3483 (2022).

52. Jean, L. The statistical power of three monkeys. *bioRxiv*, 2022.005.2010.491373, https://doi.org/10.1101/2022.05.10.491373 (2022).

53. Koprich, J. B., Johnston, T. H., Reyes, G., Omana, V. & Brotchie, J. M. Towards a non-human primate model of alpha-synucleinopathy for

development of therapeutics for Parkinson's disease: optimization of AAV1/2 delivery parameters to drive sustained expression of alpha synuclein and dopaminergic degeneration in macaque. *PLoS One* **11**, e0167235 (2016).

54. McCormack, A. L. et al. Alpha-synuclein suppression by targeted small interfering RNA in the primate substantia nigra. *PLoS One* **5**, e12122 (2010).

55. Bassil, F. et al. Viral-mediated oligodendroglial alpha-synuclein expression models multiple system atrophy. *Mov. Disord.* **32**, 1230–1239 (2017).

56. Mandel, R. J. et al. Novel oligodendroglial alpha synuclein viral vector models of multiple system atrophy: studies in rodents and nonhuman primates. *Acta Neuropathol. Commun.* **5**, 47 (2017).

57. Marmion, D. J. et al. Viral-based rodent and nonhuman primate models of multiple system atrophy: fidelity to the human disease. *Neurobiol. Dis.* **148**, 105184 (2021).

58. Jarraya, B. et al. Dopamine gene therapy for Parkinson's disease in a nonhuman primate without associated dyskinesia. *Sci. Transl. Med* **1**, 2ra4 (2009).

59. Palfi, S. et al. Lentivirally delivered glial cell line-derived neurotrophic factor increases the number of striatal dopaminergic neurons in primate models of nigrostriatal degeneration. *J. Neurosci.* **22**, 4942–4954 (2002).

60. Galvan, A. et al. Intracerebroventricular administration of AAV9- PHP.B SYN1-EmGFP induces widespread transgene expression in the mouse and monkey central nervous system. *Hum. Gene Ther.* **32**, 599–615 (2021).

61. Darricau, M. et al. Tau seeds from patients induce progressive supranuclear palsy pathology and symptoms in primates. *Brain.* **146**, 2524–2534 (2023).

62. Teil, M. et al. Brain injections of glial cytoplasmic inclusions induce a multiple system atrophy-like pathology. *Brain* **145**, 1001–1017 (2022).

63. Teil, M. et al. Cortical Lewy body injections induce long-distance pathogenic alterations in the non-human primate brain. *NPJ Parkinsons Dis.* **9**, 135 (2023).

64. Capogrosso, M. et al. A brain-spine interface alleviating gait deficits after spinal cord injury in primates. *Nature* **539**, 284–288 (2016).

65. Rosenblad, C. et al. Vector-mediated I-3,4-dihydroxyphenylalanine delivery reverses motor impairments in a primate model of Parkinson's disease. *Brain* **142**, 2402–2416 (2019).

66. Arotcarena, M. L. et al. Bidirectional gut-to-brain and brain-to-gut propagation of synucleinopathy in non-human primates. *Brain* **143**, 1462–1475 (2020).

67. Fridjonsdottir, E. et al. Mass spectrometry imaging identifies abnormally elevated brain I-DOPA levels and extrastriatal monoaminergic dysregulation in I-DOPA-induced dyskinesia. *Sci. Adv.* **7**, eabe5948 (2021).

68. Wohlke, A. et al. A one base pair deletion in the canine ATP13A2 gene causes exon skipping and late-onset neuronal ceroid lipofuscinosis in the Tibetan terrier. *PLoS Genet* **7**, e1002304 (2011).

69. Schmutz, I. et al. ATP13A2 missense variant in Australian Cattle Dogs with late onset neuronal ceroid lipofuscinosis. *Mol. Genet Metab.* **127**, 95–106 (2019).

70. Arotcarena, M. L. et al. Pilot study assessing the impact of intrathecal administration of variants AAV-PHP.B and AAV-PHP.eB on brain transduction in adult rhesus macaques. *Front Bioeng. Biotechnol.* **9**, 762209 (2021).

71. Chuapoco, M. R. et al. Adeno-associated viral vectors for functional intravenous gene transfer throughout the non-human primate brain. *Nat. Nanotechnol.* **18**, 1241–1251 (2023).

72. Dehay, B., Dalkara, D., Dovero, S., Li, Q. & Bezard, E. Systemic scAAV9 variant mediates brain transduction in newborn rhesus macaques. *Sci. Rep.* **2**, 253 (2012).
73. Bjorklund, T. & Davidsson, M. Next-generation gene therapy for Parkinson's disease using engineered viral vectors. *J. Parkinsons Dis.* **11**, S209-S217 (2021).

74. Lovaas, E. Antioxidative and metal-chelating effects of polyamines. *Adv. Pharm.* **38**, 119-149 (1997).

75. Ha, H. C. et al. The natural polyamine spermine functions directly as a free radical scavenger. *Proc. Natl Acad. Sci. USA* **95**, 11140-11145 (1998).

76. Martin, S. et al. Protection against Mitochondrial and Metal Toxicity Depends on Functional Lipid Binding Sites in ATP13A2. *Parkinsons Dis.* **2016**, 9531917 (2016).

77. Levin, J. et al. Generation of ferric iron links oxidative stress to alpha-synuclein oligomer formation. *J. Parkinsons Dis.* **1**, 205-216 (2011).

78. Dexter, D. T. et al. Alterations in the levels of iron, ferritin and other trace metals in Parkinson's disease and other neurodegenerative diseases affecting the basal ganglia. *Brain* **114**, 1953-1975 (1991).

79. Rentschler, G. et al. ATP13A2 (PARK9) polymorphisms influence the neurotoxic effects of manganese. *Neurotoxicology* **33**, 697-702 (2012).

80. Fleming, S. M. et al. The effect of manganese exposure in Atp13a2-deficient mice. *Neurotoxicology* **64**, 256-266 (2018).

81. Davies, K. M. et al. Copper pathology in vulnerable brain regions in Parkinson's disease. *Neurobiol. Aging* **35**, 858-866 (2014).

82. Dull, T. et al. A third-generation lentivirus vector with a conditional packaging system. *J. Virol.* **72**, 8463-8471 (1998).

83. Rusanescu, G. & Mao, J. Immature spinal cord neurons are dynamic regulators of adult nociceptive sensitivity. *J. Cell Mol. Med* **19**, 2352-2364 (2015).

84. Bellet, V. et al. Proteomic analysis of RCL2 paraffin-embedded tissues. *J. Cell Mol. Med* **12**, 2027-2036 (2008).

85. Blaszczyk, L. et al. Sequential alteration of microglia and astrocytes in the rat thalamus following spinal nerve ligation. *J. Neuroinflamm.* **15**, 349 (2018).

86. Gustafsson, O. J., Arentz, G. & Hoffmann, P. Proteomic developments in the analysis of formalin-fixed tissue. *Biochim. Biophys. Acta* **1854**, 559-580 (2015).

87. Rousseau, E. et al. Targeting expression of expanded polyglutamine proteins to the endoplasmic reticulum or mitochondria prevents their aggregation. *Proc. Natl. Acad. Sci. USA* **101**, 9648-9653 (2004).

88. Solé, V. A., Papillon, E., Cotte, M., Walter, P. & Susini, J. A multiplatform code for the analysis of energy-dispersive X-ray fluorescence spectra. *Spectrochim. Acta B At. Spectrosc.* **62**, 63-68 (2007).

89. Ho, J., Tumkaya, T., Aryal, S., Choi, H. & Claridge-Chang, A. Moving beyond *P* values: data analysis with estimation graphics. *Nat. Methods* **16**, 565-566 (2019).

## Acknowledgements

We thank Dr. Gregory Porras, Marie-Laure Thiolat, Tho Hai Nguyen, and Hugues Orignac for their invaluable help conducting the reported experiments. Histological analysis was performed by the HistoCARE facility of IMN (UMR5293, unit of the CNRS and Bordeaux University). The authors thank Nathalie Biendon, for the help with sample preparation. Experiments were performed in the Service des Animaleries of the University of Bordeaux, supported by the Région Aquitaine. The synchrotron Diamond Light Source is acknowledged for providing 118 beam time (exp. SP29838). We would also

like to thank Tina Geraki for her scientific support. This study received financial support from the French government in the framework of the University of Bordeaux's IdEx "Investments for the Future" program/GPR BRAIN_2030. This work was supported by grants from ANR-21-CE18-0025-01 NOVEL, France Parkinson (B.D.) and the Parkinson's Foundation 2023 Bill and Amy Gurley Impact Award PF-IMP-941400. J.S. was supported by a Région Nouvelle-Aquitaine fellowship (France). R.K. is a recipient of a Clément Fayat Foundation fellowship (France). INSERM, CNRS, and the University of Bordeaux provided financial and infrastructural support.

## Author contributions

P.O.F., E.B. and B.D. conceived and designed the study. E.B. and B.D. performed surgeries. J.S. and S.D. performed histological and immunohistochemical analyses of the data. R.K. and M.L.A. performed imaging experiments. J.S., S.B. and B.D. performed SR-XRF microscopy investigations. M.L.A. and B.D. performed biochemical experiments. J.S., S.D., R.K., M.L.A., S.B., P.O.F., E.B. and B.D. analyzed the data. P.O.F., E.B. and B.D. wrote the paper. All authors discussed the results, assisted in the preparation, and contributed to the manuscript. All authors approved the final version of the document.

## Competing interests

E.B. is a director and shareholder of Motac Neuroscience Ltd. The other authors declare no conflicts of interest.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41531-024-00757-4.

**Correspondence** and requests for materials should be addressed to Benjamin Dehay.

**Reprints and permissions information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024