PLD3 affects axonal spheroids and network defects in Alzheimer’s disease

The precise mechanisms that lead to cognitive decline in Alzheimer’s disease are unknown. Here we identify amyloid-plaque-associated axonal spheroids as prominent contributors to neural network dysfunction. Using intravital calcium and voltage imaging, we show that a mouse model of Alzheimer’s disease demonstrates severe disruption in long-range axonal connectivity. This disruption is caused by action-potential conduction blockades due to enlarging spheroids acting as electric current sinks in a size-dependent manner. Spheroid growth was associated with an age-dependent accumulation of large endolysosomal vesicles and was mechanistically linked with Pld3—a potential Alzheimer’s-disease-associated risk gene1 that encodes a lysosomal protein2,3 that is highly enriched in axonal spheroids. Neuronal overexpression of Pld3 led to endolysosomal vesicle accumulation and spheroid enlargement, which worsened axonal conduction blockades. By contrast, Pld3 deletion reduced endolysosomal vesicle and spheroid size, leading to improved electrical conduction and neural network function. Thus, targeted modulation of endolysosomal biogenesis in neurons could potentially reverse axonal spheroid-induced neural circuit abnormalities in Alzheimer’s disease, independent of amyloid removal.

The precise mechanisms that lead to cognitive decline in Alzheimer's disease are unknown. Here we identify amyloid-plaque-associated axonal spheroids as prominent contributors to neural network dysfunction. Using intravital calcium and voltage imaging, we show that a mouse model of Alzheimer's disease demonstrates severe disruption in long-range axonal connectivity. This disruption is caused by action-potential conduction blockades due to enlarging spheroids acting as electric current sinks in a size-dependent manner. Spheroid growth was associated with an age-dependent accumulation of large endolysosomal vesicles and was mechanistically linked with Pld3-a potential Alzheimer's-disease-associated risk gene 1 that encodes a lysosomal protein 2,3 that is highly enriched in axonal spheroids. Neuronal overexpression of Pld3 led to endolysosomal vesicle accumulation and spheroid enlargement, which worsened axonal conduction blockades. By contrast, Pld3 deletion reduced endolysosomal vesicle and spheroid size, leading to improved electrical conduction and neural network function. Thus, targeted modulation of endolysosomal biogenesis in neurons could potentially reverse axonal spheroidinduced neural circuit abnormalities in Alzheimer's disease, independent of amyloid removal.
Alzh eimer's disease (AD) is a neurodegenerative condition that is characterized by widespread disruption in neural circuits and network connectivity 4 . The extracellular deposition of the β-amyloid (Aβ) peptide is thought to trigger a cascade of events, eventually leading to cognitive decline 5 . However, the cellular underpinnings linking Aβ deposition and neural network disruption are not well understood 6 , limiting the rational design of new therapies. Extensive previous research has focused on mechanisms such as synapse loss and cell death as potential causes of neural dysfunction 7,8 , and therapeutic efforts have mainly focused on strategies for extracellular amyloid removal 9 . However, an important and understudied pathological hallmark of AD is the markedly enlarged neuronal processes, traditionally termed dystrophic neurites, that are found around Aβ deposits [10][11][12] . These processes were previously shown to be all axonal rather than dendritic in origin 13,14 (Fig. 1a,b), so we therefore name them plaque-associated axonal spheroids (PAASs). Although various hypotheses regarding their development have been proposed over the years [15][16][17][18][19] , these structures have not been a major focus of mechanistic investigation and their pathophysiological importance remains uncertain.
Here we undertook a multipronged approach to investigate PAASs using high-resolution structural imaging, time-lapse intravital calcium (Ca 2+ ) and voltage imaging, as well as single-axon molecular manipulations and computational modelling. We found that hundreds of axons develop PAASs around each amyloid deposit and, rather than being associated with degenerative retracting neuronal processes, they are persistent and undergo dynamic changes in size over extended imaging intervals. PAASs markedly disrupt the propagation of action potentials (AP) by acting as electric current sinks in a spheroid-size-dependent manner, leading to abnormal long-range axonal connectivity and neuronal network dysfunction. We uncovered neuronal endolysosomal/ multivesicular body (MVB) biogenesis as a critical determinant of PAAS size, and identified the neuronal lysosomal protein PLD3-a possible genetic risk factor for AD 1,20,21 -as a key mediator of endolysosomal abnormalities and PAAS enlargement. We further show that modulation of PLD3 levels can reverse axonal conduction abnormalities and restore neuronal network function in an AD mouse model. Together, our findings suggest a paradigm in which modulation of neuronal endolysosomal biogenesis can have a substantial effect on axonal electrical conduction and neural circuit function. Thus, targeting PAAS formation could be a strategy for ameliorating neural network abnormalities in AD.  Article and 2i,l). On the basis of the average volume of individual PAASs and the total volume of the PAAS halo around plaques, we estimated that individual plaques can affect hundreds of axons on average (Extended Data Fig. 1a-c and Supplementary Discussion 1). Given the abundance of amyloid plaques in the AD brain, this suggests that substantial numbers of axons and their downstream interconnected neurons can be affected, highlighting the potential importance of PAASs as a mechanism of neural network dysfunction. In 5xFAD mice, time-lapse imaging of virally labelled axons around plaques revealed that PAASs can be very stable over intervals of up to months ( Fig. 1c and Extended Data Fig. 1d-f).
Although most PAASs increased in size over time, a substantial number decreased in volume or disappeared during this interval without the loss of the parent axon ( Fig. 1c and Extended Data Fig. 1d-f), consistent with previous reports 11,15 . This supports the idea that PAASs are not a feature of degenerating axons but, rather, are stable structures that may affect neuronal circuits for extended intervals, while at the same time having the potential for reversibility. To examine how PAASs might disrupt neuronal circuits, we computationally modelled how these spheroids could affect axonal electrical conduction. We found that the likelihood of disruption in the conduction of APs, for a particular axonal segment, markedly increased as a function of the total PAAS surface area (Extended Data Fig. 2). This in silico model predicted that PAASs behave as capacitors that function as electric current sinks for incoming APs and can therefore cause conduction blocks or prolonged delays (Extended Data Fig. 2, Supplementary Discussion 2 and Supplementary Video 1). To experimentally examine the electrical conduction properties of axons, we developed a strategy for measuring the propagation of APs in individual axons through Ca 2+ imaging in the live-mouse brain. We virally expressed the calcium sensor GCaMP6f through delivery of adeno-associated viral (AAV) vectors to one brain hemisphere (as in Fig. 1a) and performed Ca 2+ imaging of individual projection axons on the contralateral cortex (Fig. 1d). We measured AP propagation after electrically stimulating the ipsilateral hemisphere with trains of electrical pulses and compared the rise times of Ca 2+ transients at two regions of interest (ROIs) located on both sides of individual PAASs, along selected axonal segments (Fig. 1d). We found that the onset of the rise times was consistently delayed over intervals of tens to hundreds of milliseconds ( Fig. 1e-g and Supplementary Video 2). Given that we used a series of pulses of electrical stimulation to induce trains of AP spikes, we concluded that the unusually long delays in the Ca 2+ rise times observed were due to conduction blocks of a substantial proportion of individual AP spikes, once they reached individual spheroids (Extended Data Fig. 2b). By contrast, a comparison of the rise times at two ROIs in axon segments without spheroids or at ROIs located on the same side of axonal segments adjacent to spheroids ( Fig. 1e (right)) demonstrated no difference in the onset of rise times, strongly suggesting a deleterious effect of PAASs, rather than a broad spheroid-independent disruption in axonal electrical conduction. Furthermore, in addition to electrically evoked responses, we also imaged spontaneous Ca 2+ transients in individual axons and observed similar conduction abnormalities in segments with PAASs (Extended Data Fig. 3 and Supplementary Video 2). Regardless of whether the APs were the result of spontaneous or electrically induced neuronal activity, we observed that larger spheroids caused more severe conduction blocks ( Fig. 1g and Extended Data Fig. 3d). These experimental observations were consistent with our computational modelling demonstrating that the size of individual PAASs is a critical determinant of the degree of axonal conduction defects (Extended Data Fig. 2f,f′).

Disruption of long-range connectivity
Given that we found marked abnormalities in local axonal conduction around plaques, we further investigated whether this was associated with more widespread defects in long-range cortical connectivity. Thus, we developed a strategy for measuring interhemispheric conduction velocity through calcium imaging in live 5xFAD mice. To achieve this, we stereotaxically injected AAV9-Syn-GCaMP6f to label a homogeneous population of closely located cortical neurons in the somatosensory cortex, which assured comparable axonal distances to the imaged regions on the contralateral hemisphere across different mice. We then imaged contralateral projecting axons while electrically stimulating the ipsilateral GCaMP6f-labelled neurons (Fig. 1h). We found that the interhemispheric conduction velocities in wild-type (WT) mice were of a similar magnitude to those previously reported using slice electrophysiology recordings 22 . However, in 5xFAD mice, we found that Ca 2+ rise times in projecting axons were markedly delayed ( Fig. 1i-j). This suggests that local AP conduction abnormalities caused by PAASs lead to a disruption in long-range axonal conduction. To further validate that our observations, as measured using calcium imaging, reflected actual AP delays or blockades, we implemented in vivo voltage imaging using the genetically encoded voltage sensor ASAP3 23 . We intracortically injected AAV2-Syn-ASAP3 in the same manner as we did for GCaMP6f (Fig. 1k). Given the relatively low signal-to-noise ratio when imaging with genetically encoded voltage sensors 23 , it was not feasible to perform single-trial experiments to visualize AP propagation in individual axons. Instead, we stimulated axons in one hemisphere and recorded antidromic APs at neuronal cell bodies on the contralateral hemisphere ( Fig. 1k,l), which markedly improved the signal-to-noise ratio. Using this strategy, we found that, in 5xFAD mice, there was a marked increase in the electric current required to induce the interhemispheric propagation of APs through single-trial stimulations ( Fig. 1m-o), consistent with PAASs acting as current sinks. Furthermore, we also observed frequent delays in AP propagation when comparing 5xFAD and WT mice (Fig. 1p,q and Supplementary Discussion 3) in agreement with our Ca 2+ imaging experiments. Together, our imaging data as well as our computational modelling highlight the prevalence of AP conduction blocks resulting from spheroid pathology in AD-like mice.
Given the long interhemispheric distances, the probability of axons encountering amyloid plaques and developing spheroids is high. Although the density of amyloid plaques in humans is lower than that in mice, we hypothesize that the much greater axonal lengths in humans markedly increase the probability of adjacency to amyloid plaques and the likelihood of disruption in axonal connectivity. In support of this idea, through a quantitative analysis of post-mortem brains, we found a greater average number of axonal spheroids per amyloid plaque and larger PAAS size in individuals with moderate to severe AD, compared with those with mild cognitive impairment (MCI) (Extended Data Fig. 4a-c), consistent with previous literature 12,24 . Although this clinical-pathological correlation has various limitations, it suggests that both PAAS number and size could be important factors that determine the degree of neural circuit disruption and cognitive deficits in AD.

Aberrant endolysosomes drive PAAS growth
High-resolution confocal imaging revealed that, as mice aged, there was a progressive accumulation of aberrantly enlarged lysosome-associated membrane protein 1 (LAMP1)-positive vesicles (ELPVs) within axonal spheroids (Fig. 2a,b). Transmission electron microscopy imaging revealed the accumulation of a heterogeneous population of vesicular organelles, including endosomes/MVBs, endolysosomes, amphisomes and autolysosomes ( Fig. 2c and Extended Data Fig. 5), some of which could overlap with the ELPVs that we observed by optical imaging. This diversity potentially reflects distinct stages of organelle maturation through the lysosomal biogenesis and autophagic pathways 25,26 . There was a notable correlation between the presence of ELPVs and the overall size of individual spheroids (Fig. 2d). Furthermore, we also found that small PAASs were predominantly filled with vesicles that contained higher levels of the protease cathepsin D and were acidiccharacteristic of more mature lysosomes (Fig. 2e-h and Extended Data Fig. 6f,g). By contrast, as PAASs increased in size, their overall acidification and cathepsin D levels declined (Fig. 2e-h and Extended Data Fig. 6f,g), consistent with the accumulation of ELPVs, which have not acquired sufficient lysosomal proteases and acidic pH 27 . Overall, this suggests that spheroid enlargement could be mechanistically linked to the accumulation of ELPVs.
Similar to mice, in post-mortem human brains, spheroids containing enlarged vesicles with low levels of cathepsin D (Fig. 2i,j and Extended Data Fig. 6a,b) were observed by immunolabelling of the lysosomal proton pump v-ATPase subunit V0A1 (Extended Data Fig. 6e,e′), and their presence was also associated with an overall larger spheroid size (Fig. 2k). Consistent with our observation that PAAS size is inversely correlated with premortem cognitive function (Extended Data Fig. 4b), we found a similar correlation between premortem cognition, the abundance of large ELPVs and low levels of cathepsin D within PAASs ( Fig. 2o and Extended Data Fig. 6c,d). Together, our mouse and human data support a hypothesis in which the accumulation of ELPVs may drive the enlargement of axonal spheroids, leading to the disruption of axonal conduction and ultimately cognitive dysfunction.

Role of PLD3 in spheroid enlargement
We next investigated the potential mechanisms of ELPV accumulation within axonal spheroids. PLD3 is a lysosomal protein 2,3,28 that is (left) and 2 μm (right). i, ELPV size in spheroids of 10-month-old 5xFAD mice with PLD3 or GFP overexpression. n = 3 and n = 4 mice for the GFP and PLD3 groups, respectively. Each dot represents the average measurement from 500-1,000 ELPVs. j, Confocal (left two images) and expansion microscopy (right two images) images of Aβ42 immunofluorescence (red) and LAMP1-GFP (green) within spheroids. The arrows indicate Aβ42 puncta contained within ELPVs. Scale bars, 2 μm (left two images) and 5 μm (right two images). k, Confocal images showing incorporation of FM 1-43 dye into spheroids in cultured brain slices after treatment with vehicle or PitStop2. Scale bars, 10 μm. l, Quantification of FM 1-43 incorporation into PAASs after PitStop or Dynasore treatment. n = 20 spheroids for each group. The red dashed lines show regression to a sigmoid inhibition curve. m, In vivo assay to measure Aβ endocytosis into spheroids after intraparenchymal brain microinjections of fluorescently labelled Aβ42 peptide. n, Confocal images of fluorescently tagged Aβ42 (magenta) incorporated into spheroids (white). The arrows indicate Aβ42 puncta. Scale bars, 10 μm (left) and 5 μm (middle and right). o, Quantification of Aβ42 incorporation into spheroids. n = 3 mice, each with average measurements from n = 10 fields of view. Statistical analysis was performed using two-tailed Mann-Whitney U-tests (e and g) and two-tailed Welch's t-tests (i and o), and F-tests were used to compare the fitted top and bottom parameters for each group in l. Data are mean ± s.e.m.
of potential interest because it strongly accumulates in axonal spheroids in both humans and mice 28,29 (Figs. 2l and 3a) and its expression is not detectable in other cell types such as microglia and astrocytes (Extended Data Fig. 7), despite mRNA presence in glial cells 30 . Moreover, PLD3 genetic variants may increase the risk of AD, although this remains a controversial topic 1,20,21 . We found that there was an overall increased abundance of ELPVs in axonal spheroids of patients with AD who have the PLD3 variant V232M 1 (Fig. 2m-o), suggesting a role of this variant in aberrant axonal endolysosomal function. During endolysosomal maturation, the limiting membrane of late endosomes invaginates to form intraluminal vesicles (ILVs) and become MVBs that fuse with lysosomes 25 . Furthermore, autophagic degradation also requires the fusion of autophagosomes and MVBs, forming amphisomes that later fuse with lysosomes 31 . Thus, MVBs are crucial intermediate organelles that connect with various components of endocytosis, autophagy and lysosomal degradation 26,31 (Extended Data Fig. 5). PLD3 is unique among lysosomal-resident proteins in that it is sorted to the ILVs of MVBs 3 , in contrast to most lysosomal-resident proteins, which are sorted to the limiting membranes of MVBs 32 . Indeed, immunofluorescence confocal and expansion microscopy imaging of axonal spheroids showed an accumulation of a punctate PLD3 signal within the lumen of LAMP1-positive vesicular structures (Fig. 3b), similar to previous immunogold electron microscopy analysis of cultured cells 3 . Together, these observations raise the possibility that PLD3 may have a role in MVB biogenesis, thereby affecting various organelles that interact with MVBs. This may lead to the accumulation of ELPVs, driving the expansion of axonal spheroids.
To further understand the role of PLD3 in the evolution of axonal spheroid pathology, we implemented in vivo AAV2-mediated overexpression of PLD3 in neurons of 5xFAD mice (Fig. 3c). We observed that spheroids in PLD3-overexpressing axons were markedly larger than those expressing only green fluorescent protein (GFP) (Fig. 3d,e and Extended Data Fig. 8a,b′). Notably, confocal microscopy of individual spheroids revealed an increase in the number of large ELPVs, even beyond what is seen in older 5xFAD mice (Fig. 3f,g and Extended Data Fig. 8c,d). Furthermore, there was a marked increase in ELPV size within PLD3-overexpressing spheroids compared with the GFP-expressing controls (Fig. 3h,i). This manipulation was not associated with changes in amyloid plaque number or size (Extended Data Fig. 8e,f), suggesting that PLD3 overexpression does not affect the processing of amyloid precursor protein (APP), in agreement with previous in vivo and in vitro experiments 2 .
The effect of PLD3 overexpression on the accumulation of large LAMP1-positive vesicles was predominantly seen in axonal spheroids around plaques rather than in neuronal cell bodies (Extended Data Fig. 8g,h). This suggests that extracellular Aβ deposits are critical for PLD3-induced endolysosomal abnormalities in axonal spheroids. Indeed, using an antibody that specifically recognizes the Aβ42 peptide, we observed Aβ accumulation within large ELPVs in spheroids (Fig. 3j). The source of Aβ42 within ELPVs is likely to be the endocytosis of oligomeric peptides from adjacent amyloid plaques. This is supported by our data showing that spheroids are sites of active endocytosis (Fig. 3k,l), and that administration of fluorescently labelled Aβ42 to 5xFAD mice leads to robust uptake into vesicular structures within spheroids (Fig. 3m-o). Previous reports have demonstrated the formation of ELPVs in vitro after Aβ administration 33 . Furthermore, our data with PLD3 overexpression in WT mice also led to occasional formation of small axonal swellings with ELPVs (Extended Data Fig. 8i), suggesting that excessive PLD3 by itself can have detrimental effects, independent of amyloidosis. Together, these data suggest that the accumulation of PLD3, observed in both mice and humans, is mechanistically linked with endolysosomal abnormalities and the subsequent enlargement of axonal spheroids. This may be compounded by the concurrent effect of PLD3 and Aβ accumulation within the same subcellular compartments.

Pld3 deletion restores axonal conduction
To assess whether reducing PLD3 levels would ameliorate axonal spheroid pathology, we deleted Pld3 in neurons using AAV2-mediated CRISPR-Cas9 gene editing in 5xFAD mice, using either of two single guide RNAs (sgRNAs), targeting different Pld3 exons ( Fig. 4a and Extended Data Fig. 9). We found that treatment with both sgRNAs led to a marked decrease in the abundance of large ELPVs (Fig. 4b,c and Extended Data Fig. 10d), which was associated with an overall reduction in PAAS size, regardless of whether the treatment was initiated at 3 or 7 months of age in 5xFAD mice (Fig. 4d-e and Extended Data Fig. 10a-c′). These data demonstrate that Pld3 deletion at early or later stages of amyloid deposition can decrease ELPV accumulation in axons, leading to a marked reduction in spheroid growth, without any changes in amyloid plaque number or size (Extended Data Fig. 10e,f).
To test whether deletion of Pld3 in neurons and the consequent reduction in spheroid size had a beneficial effect on axonal conduction, we co-infected neurons with AAV2-U6-sgRNA(Pld3)-CAG-Tomato-P2A-Cre to delete Pld3 and AAV9-Syn-GCamP6f to implement our Ca 2+ -imaging approach for measuring interhemispheric axonal conduction (Fig. 4f). This coinfection strategy enabled us to image adjacent axons with or without Pld3 deletion within the same mouse and compare the Ca 2+ rise times in the contralateral cortex after electrical stimulation of the ipsilateral hemisphere (Fig. 4f). Using the two sgRNAs, we found that axons with Pld3 deletion had a marked improvement in the propagation of APs (Fig. 4g,h and Extended Data Fig. 11a,a′,c,d′) that approached what we observed in control non-AD mice (Fig. 1j). By contrast, PLD3 overexpression had the opposite effect, with worsening of AP propagation observed in axons with increased PLD3 compared with adjacent controls (Fig. 4i-j and Extended Data Fig. 11b,b′). Together, these data demonstrate that PLD3 reduction can ameliorate endolysosomal abnormalities in axons near plaques, leading to reduced spheroid size and restored axonal conduction properties.

PAAS reduction improves network function
To examine the effect of reversing spheroid-associated axonal conduction defects on neural circuit function, we focused on the basal forebrain nucleus of Meynert, which is a major source of cholinergic neurotransmission with widespread axonal projections that exert complex neuromodulatory effects on cortical neurons 34 . Furthermore, cholinergic networks are critical for normal cognitive function and are affected in the early stages of AD 35 . We therefore investigated the potential network effects of spheroid-induced conduction deficits and their potential reversibility by PLD3 modulation.
To achieve this, we injected AAV2-U6-sgRNA(Pld3)-CAG-Tomato-P2A-Cre to delete Pld3 in neurons of the basal forebrain in 7-month-old 5xFAD mice. To examine the effect of improved basal forebrain neurotransmission, we imaged spontaneous Ca 2+ transients during awake resting sessions in neurons of layer 2/3 of the somatosensory cortex, previously infected with AAV9-Syn-GCamP6f (Fig. 5a-c and Supplementary  Video 3). These neurons were within the immediate vicinity of projecting basal forebrain axons (Fig. 5a (right)). We observed a higher proportion of hyperactive neurons in 5xFAD mice (Fig. 5d), consistent with previous reports 36 . Moreover, we also found that 5xFAD mice showed increased correlated activity in neurons that were in proximity to each other (Fig. 5e) and displayed activity patterns with higher spatiotemporal similarities (Fig. 5f,g). These aberrant patterns of activity are relevant because they are predicted to markedly disrupt efficient information encoding 37 . Pld3 deletion in basal forebrain neurons led to a reduction in such aberrant cortical neuronal activity patterns to levels that were similar to WT controls ( Fig. 5d-g). These data demonstrate that reversing spheroid-induced AP blockades in basal forebrain projection axons, which potentially improves neuromodulatory neurotransmission, can ameliorate the aberrant patterns of activity of downstream cortical neurons.

Discussion
Here we show that hundreds of axons around each amyloid plaque develop spheroids and, rather than being retraction bulbs from degenerating axons, these structures are stable for extended periods of time and could therefore have an ongoing detrimental effect on neuronal connectivity. Given the similarity in the morphology, organelle and biochemical content of PAASs in mice and humans, it is probable that, in humans, these are also stable structures that could disrupt neural circuits for extended intervals. To better understand the effect of PAASs on axonal function, we implemented in vivo Ca 2+ and voltage imaging in individual cortical axons and cell bodies. Both Ca 2+ and voltage imaging revealed that a substantial proportion of axons in 5xFAD mice had disrupted AP conduction and an overall increase in the threshold for AP propagation manifested by conduction blockades. This was due to the presence of axonal spheroids and was shown to be correlated with their size. The finding that larger PAASs caused more severe conduction blocks was consistent with computational modelling (Supplementary Discussion 2) showing that PAASs resemble electrical capacitors that function as current sinks, and that PAAS size is a major determinant of the degree of conduction defects. Together, our data suggest that the large number of amyloid deposits present in the AD brain have the potential to substantially affect neural networks by widespread disruption of axonal connectivity. Using in vivo Ca 2+ imaging, we found that neurons in the cortex of 5xFAD mice showed a pattern of aberrant activity, in agreement with previous reports 36 . We found that these abnormal activity patterns can be rectified by deleting Pld3 from basal forebrain neurons, which provide cholinergic inputs to cortical neurons. These results raise the possibility that the aberrant neuronal cortical activity could be partly due to conduction defects in long-range cholinergic projections caused by PAASs. The potential relevance of these findings is further highlighted by the observation of neuronal hyperactivity in humans in early stages of AD 38 and the susceptibility of cholinergic neurons to AD neuropathology 35 .
In addition to cholinergic projections, strategically located amyloid plaques could have deleterious effects in other brain regions, such as the hippocampus, in which parallel compact axonal bundles follow a stereotyped projection path along a trisynaptic loop. Memory formation may be particularly vulnerable to axonal conduction delays and blockades, given that memory consolidation depends on the precise timing of hippocampal replay and sharp-wave ripples 39 (Supplementary Discussion 4). Consistent with this view, our quantitative histopathology analysis in a limited number of human post-mortem brain samples from AD or MCI showed that PAAS size and number correlate well with the degree of premortem cognitive decline.
Mechanistically, we found that ELPVs-which probably include MVBs, endolysosomes and autolysosomes-accumulate within axonal spheroids and that their presence is correlated with spheroid size. Moreover, we found an increased presence of ELPVs within spheroids in older 5xFAD mice and in more severely impaired human patients with AD, indicating that ELPV accumulation may be a key feature of disease progression. MVBs are crucial intermediate organelles that evolve through the maturation of endosomes and fuse with autophagosomes and lysosomes 26,31 . Thus, dysregulation in MVB biogenesis has the potential to affect the normal generation of fusion vesicles.
The endosomal sorting complex required for transport (ESCRT) machinery has a major role in MVB biogenesis by regulating the formation of ILVs within MVBs and the sorting of ubiquitinated proteins into ILVs destined for degradation 25 . In contrast to lysosomal-resident  Article proteins sorted to the limiting membrane of MVBs 32 , PLD3 is unique because it is sorted into ILVs through the ESCRT pathway in mammals 3 . Indeed, we found that PLD3 was present within ELPVs and was highly enriched in PAASs, suggesting a potential role of PLD3 in MVB biogenesis. Consistent with this, we found that overexpression of PLD3 in neurons led to a marked enlargement and accumulation of ELPVs and resulted in an overall increase in PAAS size. Furthermore, the human PLD3(V232M) variant was associated with an increased abundance of ELPVs within PAASs, similar to PLD3 overexpression in mice, suggesting that this variant may exert a gain-of-function effect. A study in young Pld3-knockout mice showed that classical lysosomes in neuronal cell bodies were enlarged 2 , whereas a separate study in HeLa cells showed no lysosomal changes with PLD3 deletion 3 . Together with our data, this suggests that PLD3 has complex effects on endolysosomes that are age-and context-dependent and that axons are particularly susceptible to endolysosomal abnormalities, especially in the presence of amyloidosis.
Dysfunction of ESCRT components can lead to the enlargement of endolysosomal compartments 3,40 . Given the accumulation of PLD3 in axonal spheroids, it is possible that PLD3 accumulation leads to endolysosomal enlargement by interfering with ESCRT machinery. Consistent with this, we found that PLD3 overexpression led to the enlargement of LAMP1-positive vesicles and the formation of small axonal swellings even in WT mice. However, the substantial enlargement of ELPVs after PLD3 overexpression predominantly occurred in spheroids in the vicinity of amyloid plaques, suggesting that this process is amyloid-dependent. We observed robust endocytic activity at axonal spheroids that was associated with the uptake of Aβ42 into endolysosomal compartments. Moreover, we found that Aβ42 was present within ELPVs at axonal spheroids, consistent with previous immunogold electron microscopy in patients with AD, showing that the most prominent subcellular localization of Aβ42 is within MVBs 41 . Consistent with this, administration of Aβ42 to cultured neurons has been shown to result in MVB enlargement, possibly through interference with ESCRT proteins 33 . Thus, internalization of Aβ from extracellular deposits may be critical for PLD3-induced ELPV accumulation.
Aβ42 accumulation within MVBs may also impair the ubiquitin-proteasome system, leading to defects in protein sorting and ILV invagination 42 . Protein-sorting abnormalities could in turn exacerbate the accumulation of PLD3. Thus, PLD3 could work synergistically with Aβ42 in the same subcellular compartment, leading to greater endolysosomal abnormalities. On the other hand, given that APP and β-site APP-cleaving enzyme (BACE1) also accumulate within PAASs 43 , we do not exclude the possibility that intracellularly produced Aβ could also contribute to abnormalities in MVB biogenesis. APP is also sorted into ILVs of MVBs through the ESCRT machinery 40 , and deletion of ESCRT components promotes APP processing and increased intracellular Aβ accumulation 40 . Thus, PLD3 and Aβ may constitute a cycle in which axonal endocytosis and/or intracellularly produced Aβ facilitate the generation and accumulation of ELPVs (Extended Data Fig. 12). However, the precise molecular mechanisms explaining the potential interactions between PLD3 and Aβ remain to be elucidated.
Although the focus of this investigation was on PLD3, it is probable that manipulation of other proteins in the endolysosomal pathway could lead to changes in PAAS size. However, given the negligible expression of PLD3 in non-neuronal cells 29 (Extended Data Fig. 7), this molecule could be a promising therapeutic target because global modulation of the endolysosomal pathway may negatively affect glial cells and their roles in controlling protein aggregation and amyloid brain accumulation 44 . Thus, modulation of neuronal MVB biogenesis through PLD3 or other endolysosomal molecules in neurons could constitute strategies for ameliorating PAAS pathology, independent of amyloid plaque removal. Furthermore, although cytoskeletal abnormalities in axons were not the focus of this study, they probably have important roles in the accumulation of vesicles and spheroid enlargement.
Moreover, the presence of hyperphosphorylated tau protein within PAASs as a potential source of tau propagation throughout the neuronal soma 45 suggests important links between plaques, axonal spheroids and neurofibrillary tangles that need further investigation. Finally, in addition to neuronal intrinsic mechanisms, the glial microenvironment around plaques has been shown to have a key role in preventing PAAS formation 46,47 . Thus, the interplay between intrinsic neuronal and extrinsic glial mechanisms may contribute to the formation and enlargement of PAASs and should be considered when designing therapies.
Our findings reveal a cell-intrinsic neuronal mechanism that modulates the size of axonal spheroids and the consequent axonal conduction defects, with important implications for neuronal network dysfunction. Given that axonal spheroids are a prominent feature in various neurological disorders in addition to AD [48][49][50] , it remains to be studied whether spheroids in these conditions share mechanistic properties with those in AD. Thus, our study opens a theoretical and experimental framework for systematically investigating axonal spheroid pathology in a variety of neurological conditions.

Online content
Any methods, additional references, Nature Portfolio reporting summaries, source data, extended data, supplementary information, acknowledgements, peer review information; details of author contributions and competing interests; and statements of data and code availability are available at https://doi.org/10.1038/s41586-022-05491-6.

AAV production and delivery
The GCaMP6f and GCaMP6s viruses were purchased (UPenn Virus Core, AV-9-PV2822 and AV-9-PV2824; Addgene, 100837 and 100843). The ASAP3 construct was purchased from Addgene (Addgene, 132331). Customized AAV vectors for overexpression were constructed based on plasmid 28014 from Addgene, in which the GFP sequence was deleted and replaced by the customized sequences described below. In many cases in which the virus transduced both a target protein and a fluorescent protein reporter, a GFP without the stop codon and P2A sequence was placed in front of the target protein sequence in the same open reading frame, as described previously 51 . The target proteins used in this study were as follows: tdTomato: sequence is available online (http://www.tsienlab.ucsd.edu/ Samples/PDF/tdTomato-map%20&%20sequence.pdf), synthesized by Integrated DNA Technologies; mCherry-SEpHluorin: the sequence was cut from Addgene 32001; Lamp1-GFP: the Lamp1 sequence was amplified from mouse brain mRNA, using the 5′ primer TGCGTCGCGCCATG-GCGGCC and 3′ primer GATGGTCTGATAGCCGGCGT; GFP-P2A-Pld3: the Pld3 sequence was amplified from mouse Pld3 cDNA (GE open biosystem), using the 5′ primer ATGAAGCCCAAACTGATGTACCAGG and the 3′ primer TCAAAGCAGGCGGCAGGC. The sgRNA constructs for Pld3 deletion were cloned using plasmid 60229 from Addgene. The sequences of the sgRNAs were: Pld3 sgRNA 1: GTCCTGATCCTGG CGGTAGT; Pld3 sgRNA 2: GCTAGTGGAGGGGTTGCTCG and control scrambled sgRNA: GGAAGAGCGAGCTCTTCT. All of the constructs were verified by DNA sequencing, and the expression or deletion of the target proteins was tested using immunohistochemistry.
AAV2 viruses were produced in HEK293T cells (American Type Culture Collection (ATCC)) and purified according to previously described procedures 52 using a two-plasmid helper-free system (PlasmidFactory). The virus titre was determined by counting the infection in HEK293 cells. AAV vectors were injected into the subarachnoid space in one hemisphere as previously described 51 . The total viral particles injected per mouse was approximately 10 7 .

Cranial window implant
Eight-month-old 5xFAD mice were anaesthetized with ketamine/xylazine solution (100 mg per kg and 10 mg per kg, respectively) and hair was removed from the skull area. Buprenex (0.1 mg per kg), dexamethasone (2 mg per kg) and carprofen (5 mg per kg) were given subcutaneously at this point. The mouse was placed on a heating pad during the surgery and anaesthesia was checked periodically. Povidone-iodine solution was applied to the skin and cleaned with ethanol, and eye ointment was applied to the eyes. A small piece of skin was removed to expose the skull, and the membranous layer on the skull surface was removed by forceps. A 4 mm diameter circle was drilled on the contralateral hemisphere of virus infusion (approximate location of the centre is −2.5 mm from bregma and 2.5 mm from the midline). The skull was rinsed with sterile PBS periodically to avoid excessive heating. The skull was thinned in a circumferential area and then lifted with fine forceps without causing injury to the underlying pila surface. Gelfoam sponge (Pfizer) was used to absorb blood after lifting the skull. Using a pair of very fine forceps, the dura was removed within the circle area and a 4 mm cover glass was gently pressed onto the brain surface and glued to the skull. A customized head bar was glued (for acute imaging) or chronically implanted (with dental cement, for chronic imaging) onto the skull. For chronic imaging, mice were placed onto a heating pad to recover after the surgery and given buprenex (0.1 mg per kg) and carprofen (5 mg per kg) for 3 days. Imaging procedures started 1 month after the surgery.

In vivo two-photon imaging
In vivo imaging was performed using a two-photon microscope equipped with a Ti-sapphire tuneable laser (Spectra Physics), a gallium arsenide phosphide (GaAsP) detector (Prairie technology) and a ×20/1.0 NA water-immersion objective (Leica), or the Ultima Investigator multiphoton microscope (Bruker) with the Insight X3 tuneable ultrafast laser (Spectra Physics) and a ×20/1.0 NA water-immersion objective (Olympus), using the associated Prairie View software. GFP was excited at 920 nm; dTomato and tdTomato were excited at 920 nm/1,045 nm; and FSB was excited at 850 nm. For chronic imaging, a location close to the centre of the cranial window was selected as the starting point and the blood vessel pattern was recorded. The coordinates of each ROI were recorded as well. To relocate in the next imaging session, the starting point was relocated on the basis of the recorded coordinates and the field of view was adjusted to match the recorded blood vessel pattern.

Aβ(1-42) preparation and injection
Fluorescently tagged Aβ(1-42) peptide (AnaSpec, 60480-1) was reconstituted in DMSO to a final concentration of 1 mg ml −1 . The solution was 1:10 (v/v) diluted in fresh artificial cerebrospinal fluid before being injected into the subarachnoid space as previously described 53 . We injected 10 μl of Aβ(1-42) solution into each mouse and collected the brain the next day. Brain tissue was prepared for immunohistochemistry analysis.

Calcium imaging of cortical axons and related analysis
Six-to-eight-month-old 5xFAD mice were injected with GCaMP6 virus through the subarachnoid space on one hemisphere to label cortical neurons and measure local axonal conduction properties. For measurements of interhemispheric axonal conduction, GCaMP6f virus was injected stereotaxically with the following coordinates: bregma (AP, −0.34; ML, 1.65; DV, 0.45; angle, 0°) 54 (Allen Mouse Brain Connectivity Atlas (2011)). After more than 2 weeks of the injection, an acute cranial imaging window was implanted onto the contralateral hemisphere as described above. For stimulated calcium imaging, an additional opening on the skull was made on the ipsilateral side of the virus infusion. A glass electrode was inserted through this opening using a motorized micromanipulator and used for electrical stimulation.
The ROI was identified under a two-photon microscope. GCaMP6-labelled neurons were imaged through excitation at a wavelength of 920 nm. A limited field of view was used to improve the sampling rate. GCaMP6s was imaged at 2 Hz and GCaMP6f was imaged at 10 to 20 Hz. Only axons that displayed spontaneous calcium transients at least once per minute were used for analysis. For stimulated calcium events, mice were anaesthetized using 0.5% isoflurane. Stimulation trains of 2 ms pulses were delivered to the glass electrode at 50 Hz (18 ms interval) with 10 to 60 μA currents for 500 ms. The calcium responses within the imaging window were monitored after stimulation. The stimulating electrode was adjusted to different depths within the cortex to trigger responses in contralateral hemisphere axons. Three consecutive trials were acquired for each axon and the responses were averaged.
The raw GCaMP6 fluorescence intensity was normalized to ΔF/F for analysis. Images were then spatially smoothed with a 3 × 3 window. Several ROIs were selected on each axon. The average ΔF before stimulation was used as the baseline measurement. To estimate the calcium rise time (a surrogate for AP spike time 55 ), the calcium trace from the event-specific peak to the first data point exceeding the baseline was used as the rising phase of the event. This rising phase trace was fitted to an exponential equation: Y = 1 − exp(−k × (x − t)). The spike timing estimation (t 0 ) was then calculated by extrapolating the x-intercept. For analysis of the spontaneous Ca 2+ transients, we calculated the correlation coefficients between two ROIs chosen on each axonal side of a particular spheroid using the Pearson correlation coefficient.
Calcium imaging of cortical neuronal networks and related analysis AAV2 viruses encoding either Pld3 sgRNA or control sgRNA were injected stereotaxically into the basal forebrain of 5xFAD mice (aged 6 to 8 months) with the following coordinates: bregma (AP, 0.62; ML, 1.2; DV, 4.85; angle, 0°) 54 (Allen Mouse Brain Connectivity Atlas (2011)). GCaMP6 virus was injected through the subarachnoid space on the ipsilateral hemisphere of basal forebrain to label cortical neurons.
The ROI with intermingled projecting axons from the basal forebrain and GCaMP6f-labelled cortical neurons was identified under a two-photon microscope. Calcium imaging was performed in cortical neurons of awake mice in the same region as those projecting forebrain axons. Spontaneous calcium activities of neurons were recorded at 25 Hz for 30 min.
To analyse the cortical neuron activities, we applied a rigid motion correction 56 to the raw time-lapse data and segmented the neuronal signal using a previously established algorithm based on robust estimation 57 . Cell locations were defined as the centroid points from the spatial masks of each cell. We further estimated the neuronal spikes on the basis of calcium traces as previously shown 58 . To calculate the pairwise mutual information, we used previously described methodologies 59 . Mutual information from each cell pair in the same imaging session was grouped by the distances between neurons. To calculate neuron clusters with similar activity patterns, we used previously described graph-theory-based community detection algorithms 60 , with forced deterministic behaviour for better reproducibility, and determined the sizes of each cluster of cells identified in a single imaging session. Calculation of population entropy as performed according to a previously reported method 61 . To account for different numbers of cells imaged, we drew a random subsample of 100 neurons from each mouse and calculated the population entropy. This process was repeated 100 times and the average results were recorded.
Voltage imaging and analysis AAV2-Syn-ASAP3 virus was injected stereotaxically into 5xFAD mice (aged 6 to 8 months) with the following coordinates: bregma (AP, −0.34; ML, 1.65; DV, 0.45; angle, 0°) 54 (Allen Mouse Brain Connectivity Atlas (2011)). After more than 2 weeks of the injection, an acute cranial imaging window was implanted onto the ipsilateral hemisphere as described above. An additional opening on the skull was made on the contralateral side of the virus infusion. A glass electrode was inserted through this opening using a motorized micromanipulator and used for electrical stimulation.
The ROI was identified under a two-photon microscope. Line scan imaging of ASAP3-labelled neuronal soma was performed through excitation at a wavelength of 920 nm. The scanning speed was set to 1 kHz. Stimulation trains of 5 ms pulses were delivered at 10 Hz (95 ms intervals) with 10 to 100 μA currents for 1 s.
Line scan images were used to extract the voltage response of ASAP3. On each line of the kymograph, the intensity of the pixels covering the ASAP3-labelled neuronal soma was averaged to represent the membrane potential of the neuron at given time, which was used to generate the voltage response trace later. The response trace showed a large dip after each electrical stimulation pulse. The time of the dip indicated the antidromic AP time 62 on the neuronal soma. To quantify the AP conduction failure, we performed FFT of the voltage response curve and took the 10 Hz (the same frequency as the electrical stimulation) component power as an indicator. The amplitude of the electrical current pulse was gradually increased; and the FFT power-current amplitude could be plotted. We fit the plot with a logistic function: Y = a/ (1 + exp(−k × (x − b))). The current amplitude at the half-height of the logistic curve was defined as the threshold of the electrical stimulation.

Human post-mortem brain tissues
Formalin-fixed human post-mortem brain tissue blocks were acquired from brain banks, with ethics approval from each providing organization. The middle frontal gyrus, a cortical region affected in early stages of Alzheimer's disease 63 , was used for this study. Detailed information is provided in Extended Data Fig. 4, including 12 cases of AD, four cases of AD with PLD3 variants and six cases of MCI. Cases were matched for age, sex and APOE genotype for most quantifications except for PLD3(V232M)-variant samples. For immunohistochemistry analysis of human tissue, 30-μm-thick slices were prepared and treated with sodium citrate solution at 95 °C for 45 min, before staining with primary antibodies for 3 days.
Acute organotypic brain-slice culture Brain-slice cultures were prepared from 5xFAD mice (aged 8 months) according to a previously established protocol 64 . In brief, we dissected the hippocampal region in a sterile hood from anaesthetized mice. We then manually cut coronal sections approximately 300 μm thick and transferred the slices onto a Millicell culture membrane (Thermo Fisher Scientific, PICM03050). The culture membrane was placed in a six-well plate filled with 1 ml of the culture medium, and was placed in an incubator at 37 °C under 5% CO 2 . We examined the slice condition after 7 days with a light microscope. Healthy slices were then used for experiments. To measure endocytosis, we added endocytosis marker FM 1-43 dye (Thermo Fisher Scientific, T35356) to the culture medium to a final concentration of 1 μM. We incubated the slices with the dye for 1.5 h and then washed the slices with fresh medium twice before fixing the slice in 4% paraformaldehyde. For experiments blocking endocytosis, we dissolved Dynasore (Tocris Bioscience) or the endocytosis inhibitor PitStop2 (Abcam, ab120687) in dimethyl sulfoxide (DMSO) and added it to the culture medium with different final concentrations of FM 1-43 dye. The control groups used DMSO with no drug following the same dilutions.

Plaque-associated axonal spheroid imaging and quantification
Fixed-tissue imaging was performed using a confocal microscope (Leica SP5 or Lecia SP8), and the images were taken using a ×63/1.4 NA oil-immersion objective (Leica). Individual PAAS sizes were measured on the basis of LAMP1, APP or V0A1 immunohistochemistry. Specifically, we examined the cross-sections of the target PAASs through the high-resolution z-stack and selected the outlines of individual cross-sections using NIH/Fiji software. The largest cross-sectional area was recorded as the PAAS size. For quantification of the proportion of PAASs with ELPVs, we defined PAASs in a binary manner as either having or not having at least one obvious LAMP1-positive enlarged ring, consisting of a clear LAMP1-labelled circular outline and a darker lumen with an area greater than 0.25 μm 2 , as measured using NIH/ Fiji. To classify PAASs into neutral versus acidic, using the genetically encoded pH sensor SEpHluorin 65 , we established an arbitrary threshold of green/red fluorescence intensity ratio of 0.75. To classify PAASs into cathepsin D low versus high groups using immunofluorescence, we used an arbitrary threshold of 100 fluorescence intensity as measured using NIH ImageJ/Fiji.

Expansion microscopy
Expansion of brain sections was performed according to conventional immunostaining. Brain sections were treated with glutaraldehyde (TCI Chemicals, G0068) and then processed for gelation, digestion and expansion, as described previously 66,67 . In brief, brain sections were first incubated with monomer solution (1× PBS, 2 M NaCl, 8.625% (w/w) sodium acrylate, 2.5% (w/w) acrylamide, 0.15% (w/w) N,N′-methylenebisacrylamide) at 4 °C for 45 min. The sections were then transferred into a gel chamber and incubated in gelling solution (concentrated stocks (10%, w/w) of ammonium persulfate initiator and tetramethyl-ethylenediamine accelerator added to the monomer solution for up to 0.2% (w/w) each and the inhibitor 4-hydroxy-2,2,6, 6-tetramethylpiperidin-1-oxyl added up to 0.01% (w/w) from a 0.5% (w/w) stock) at 37 °C for 1.5-2 h for gelation. The gels were then fully immersed in proteinase solution (proteinase K (New England Biolabs, P8107S) diluted 1:100 to 8 U ml −1 in digestion buffer (50 mM Tris (pH 8), 1 mM EDTA, 0.5% Triton X-100, 1 M NaCl)) at 37 °C overnight. Digested gels were next placed in excess volumes of double deionized water (double-distilled H 2 O) for 25 min to expand. This step was repeated 3−5 times in double-distilled H 2 O, until the size of the expanding sample plateaued.

Transmission electron microscopy
Twelve-month-old 5xFAD mice were perfused with 4% PFA and the brain tissues were sectioned into 50-μm-thick slices using a vibratome (VT1000S, Leica). The slices were refixed in 2% glutaraldehyde in 0.1 M cacodylate buffer (pH 7.4) for 1 h, then post-fixed in 1% OsO 4 in the same buffer at room temperature for 1 h. After en bloc staining with 2% aqueous uranyl acetate for 30 min, tissues were dehydrated in a graded series of ethanol to 100%, followed by propylene oxide and finally embedded in EMBed 812 resin. Tissue blocks were polymerized overnight in an oven at 60 °C. Thin sections (60 nm) were cut by a Leica ultramicrotome (UC7) and post-stained with 2% uranyl acetate and lead citrate. Sample grids were examined on the FEI Tecnai transmission electron microscope with an accelerating voltage of 80 kV, and digital electron micrographs were recorded with an Olympus Morada CCD camera and iTEM imaging software.

AAV-mediated molecular manipulations
AAV vectors were injected into approximately 3-month-old and approximately 7-month-old 5xFAD mice. AAVs were infused through the subarachnoid space as previously shown 51 . Brain tissues were collected around 1.5 months and 3 months after virus injection for treatments initiated at 3 and 7 months of age, respectively, and fixed with 4% paraformaldehyde. Brain slices of 50 μm thickness were prepared and stained with anti-LAMP1 antibodies (DSHB, 1D4B) and thioflavin S.
Imaging and quantification of PAASs were performed as described previously 51,66 . In brief, tiled z-stack images of the infected cortical regions were taken at zoom 1 and 1 μm z-steps. Individual plaques were segmented from the tiled images and blinded to the mouse and treatment information. For measurement of PAAS bulb size, the z plane with the largest cross-sectional area of each individual spheroid was selected, and the cross-sectional area was measured using NIH ImageJ/ Fiji by manually selecting the outlines of that cross-section based on the outline of virally expressed cytoplasmic GFP fluorescence.

Computational modelling
All modelling experiments of axon spheroids used methodologies that were previously validated 68,69 . The chosen morphology parameters and ion-channel distributions were held constant within each compartment for each simulation. The axon length was 566 μm and diameters were set from 0.1 to 0.9 μm. The spheroids were modelled as a 'cylinder and stick'. The stick (5 μm length and diameter varying from 0.3 to 8.3 μm) was connected to the middle of the axon and on the other end to the cylinder (a bulb head) of which the surface area varied from 0 to 7,100 μm 2 (equivalent spherical diameter varied from 0 to 150 μm). In most simulations, the axon contained voltage-gated sodium and potassium channels and a leak current. The densities of the voltage-gated channels were varied from 0 to an amount 1.1× as strong as needed to produce regenerative (sustained) APs in the axon. Channel conductances were set in the spheroids with the same (varying) strength as was present in the axon. One or more strong-current injections (0.2 ms, 2 nA) were applied to one end of the axon to reliably evoke a single or multiple input AP(s) that subsequently propagated to the spheroids. Current injections just over the threshold generated qualitatively similar results (larger diameter axons required more current to evoke an AP). We checked for the presence of output APs on the other side of the spheroids by testing for voltages exceeding 10 mV. We counted the number of output APs that were present in a 10 s simulation where either the single input AP occurred at, or the train of 20 Hz input APs began at 100 ms.

Statistics and reproducibility
The methods to analyse ELPVs, cathepsin D levels or the pH in individual axon spheroids, as well as spheroid size after treatment, the specific numbers of plaques or spheroids measured from an individual human or mouse sample for each experiment are provided in the figure legend of the particular experiment, and the average results were used as the representative outcome for that individual. The number of post-mortem human tissues or mice was used as the sample size in these cases. To analyse differences in the calcium rise time between ROIs along axon segments, the number of axons was used as the sample size. This is because individual in vivo experiments had very few axons measured due to the necessary sparse labelling method used in these experiments. To analyse the interhemispheric calcium rise time delay, both the number of axons and the number of mice were used as sample sizes. For analysis of interhemispheric conduction using a voltage sensor, the number of cell bodies were used as the sample size. For analysis of neural circuit function, the number of mice was used as the sample size.
In all graphs, individual data points are shown. In all statistical comparisons, non-parametric tests were used unless otherwise justified. Specific tests used for each graph can be found in the corresponding figure legend. When more than two groups were considered and compared, corrections for multiple comparisons were performed as part of the post hoc analysis. All statistical analysis was performed using GraphPad Prism.
Representative confocal images shown in the figures were repeated independently by at least two experimenters in multiple mice or human samples.

Inclusion and ethics statement
We balanced both sexes in all of our experiments to the best of our ability. No sex-specific effects of axonal spheroids on conduction, or modulation of PLD3 on the degree of axonal spheroids were observed in our study. The brain tissues were obtained from the Sun Health Research Institute Brain and Body Donation Program, the Mayo Clinic, the University of Washington Alzheimer's Disease Research Center Neuropathology Core and the Alzheimer's Disease Research Center at Washington University in St Louis, in complete compliance with the human tissue research use regulation in each organization.

Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability
The raw datasets generated and analysed during this study are not deposited to a publicly accessible repository due to the large amount of storage required but are available from the corresponding author on reasonable request.

Code availability
All of the custom codes for FIJI and MATLAB used in this study have been deposited at GitHub (https://github.com/PaulYJ/Axon-spheroid). The scripts are free to the public to view, distribute and use. Our NEURON simulation environment computer code is available at ModelDB (http:// modeldb.yale.edu/187612).  Fig. 3 | Axonal spheroids markedly disrupt spontaneous action potential conduction. a, Two-photon in vivo calcium imaging of spontaneous activity in axons near amyloid plaques (blue) with and without PAAS. Given the relatively low frequency of spontaneously active neurons that can be captured in the vicinity of plaques, we used lower frame rates to image larger fields of view and were thus unable to measure precisely the Ca 2+ rise times like in Fig. 1. b, Example traces of GCaMP6s fluorescence signal were obtained from ROIs (green and magenta circles) at the two sides of the plaques indicated in a (2Hz imaging frame rate). Mismatched Ca 2+ transients are indicated with orange arrows. c, Correlation maps were calculated using the average fluorescence intensity within ROI1 (green circle in a, as reference, and colour-coded for correlation coefficient to every other pixel within the field of view. d, Decorrelation of GCaMP6s fluorescence in ROIs at the two axonal sides with respect to the plaque, during spontaneous Ca 2+ transients (n = 12 axons without PAAS; and n = 10 axons with PAAS). Two-tailed Mann-Whitney test was used for comparison.