Neurofilament light chain is increased in the parahippocampal cortex and associates with pathological hallmarks in Parkinson’s disease dementia
Translational Neurodegeneration volume 12, Article number: 3 (2023)
Increased neurofilament levels in biofluids are commonly used as a proxy for neurodegeneration in several neurodegenerative disorders. In this study, we aimed to investigate the distribution of neurofilaments in the cerebral cortex of Parkinson’s disease (PD), PD with dementia (PDD) and dementia with Lewy bodies (DLB) donors, and its association with pathology load and MRI measures of atrophy and diffusivity.
Using a within-subject post-mortem MRI-pathology approach, we included 9 PD, 12 PDD/DLB and 18 age-matched control donors. Cortical thickness and mean diffusivity (MD) metrics were extracted respectively from 3DT1 and DTI at 3T in-situ MRI. After autopsy, pathological hallmarks (pSer129-αSyn, p-tau and amyloid-β load) together with neurofilament light-chain (NfL) and phosphorylated-neurofilament medium- and heavy-chain (p-NfM/H) immunoreactivity were quantified in seven cortical regions, and studied in detail with confocal-laser scanning microscopy. The correlations between MRI and pathological measures were studied using linear mixed models.
Compared to controls, p-NfM/H immunoreactivity was increased in all cortical regions in PD and PDD/DLB, whereas NfL immunoreactivity was increased in the parahippocampal and entorhinal cortex in PDD/DLB. NfL-positive neurons showed degenerative morphological features and axonal fragmentation. The increased p-NfM/H correlated with p-tau load, and NfL correlated with pSer129-αSyn but more strongly with p-tau load in PDD/DLB. Lastly, neurofilament immunoreactivity correlated with cortical thinning in PD and with increased cortical MD in PDD/DLB.
Taken together, increased neurofilament immunoreactivity suggests underlying axonal injury and neurofilament accumulation in morphologically altered neurons with increased pathological burden. Importantly, we demonstrate that such neurofilament markers at least partly explain MRI measures that are associated with the neurodegenerative process.
Parkinson’s disease (PD) is a common neurodegenerative disease with a heterogeneous clinical presentation, which is diagnosed when bradykinesia and tremor and/or rigidity are present . Cognitive impairment is common in PD, leading to Parkinson’s disease dementia (PDD) in up to 80% of PD patients [2,3,4,5]. However, when patients develop dementia before or within one year of motor symptom onset, the disease is diagnosed as dementia with Lewy bodies (DLB) [6, 7]. Pathologically, PD, PDD and DLB are defined as synucleinopathies , since they are characterized by the accumulation of alpha-synuclein (αSyn) in the form of Lewy bodies (LBs) and Lewy neurites (LNs), which are abundant in the cortex of demented patients . In addition, PDD and DLB patients frequently show Alzheimer’s disease (AD) co-pathology, such as amyloid-beta (Aβ) plaques and phosphorylated-tau (p-tau) neurofibrillary tangles (NFT) [10, 11]. The load of AD co-pathology varies greatly between cases, ranging from cases with no to high AD neuropathological changes [6, 7, 12], but generally DLB patients are likely to show high prevalence of AD neuropathological changes . However, the pathological hallmarks of PDD and DLB largely overlap, and only few neuropathological differences between the two have been described [6, 13, 14].
In addition to pathology load, there is growing evidence for the presence of axonal degeneration in PD, PDD and DLB [15, 16]. Neurofilaments in cerebrospinal fluid (CSF) and blood are commonly used as a proxy of axonal degeneration in many neurodegenerative disorders [15, 16]. Neurofilaments are neuronal-specific proteins which consist of neurofilament light (NfL), medium (NfM) and heavy chains (NfH) . CSF and blood NfL levels have been shown to be elevated in PD, more abundantly in PDD and DLB [17,18,19,20,21], and associated with higher levels of CSF αSyn and p-tau biomarkers  as well as with measures of cognitive decline [16,17,18, 20, 23, 24]. However, the distribution pattern of NfL and phosphorylated-neurofilament medium and heavy chain (p-NfM/H) in cortical brain regions, and to which extent these neurofilaments are related to pathological accumulation in PD, PDD and DLB are yet to be determined.
Atrophy and altered microstructural integrity of the cortex can be captured with MRI outcome measures such as cortical thickness and mean diffusivity (MD), and may be a proxy of underlying neurodegeneration . While PDD and DLB patients show more pronounced cortical atrophy and increased cortical MD compared to cognitively unimpaired elderly [26,27,28,29], PD patients show only subtle changes [30,31,32]. Cortical atrophy and increased cortical MD have been shown to correlate with plasma NfL in several cortical regions in de-novo PD . However, little is known about the relation between regional distribution of neurofilaments and MRI-derived microstructural measures in the brains of PD, PDD and DLB.
The aim of the current study was to unravel the regional distribution of cortical neurofilaments, and its association with pathology load and MRI measures of cortical thickness and diffusivity in the brains of clinically-defined and pathologically-confirmed PD, PDD/DLB and control brain donors, using a within-subject post-mortem MRI-pathology approach . Results of this study will increase knowledge on the regional distribution of neurofilaments in the cortex in PD, PDD and DLB, and to what extent this is related to pathological accumulation and reflected by cortical thickness and diffusivity measures, thereby contributing to the understanding of the pathological underpinnings of MRI.
Materials and methods
In collaboration with the Netherlands Brain Bank (NBB; http://brainbank.nl) we included 21 Lewy Body disease donors who could be further subdivided into 9 PD, 7 PDD and 5 DLB based on clinical presentation [1, 5, 6]. Age at diagnosis (at symptom onset) and disease duration (age at death minus age at diagnosis) were extracted from the clinical files of all donors. Neuropathological diagnosis was confirmed by an expert neuropathologist (A.R.) and performed according to the international guidelines of the Brain Net Europe II (BNE) consortium (https://www.brainbank.nl/about-us/brain-net-europe/) [35,36,37]. Additionally, 18 age- and sex-matched and pathologically confirmed controls with no records of neurological symptoms at diagnosis, were selected from the Normal Aging Brain Collection Amsterdam (NABCA; http://nabca.eu) . All donors signed an informed consent for brain donation and the use of material and clinical information for research purposes. The procedures for brain tissue collection of NBB and NABCA have been approved by the Medical Ethical Committee of Amsterdam UMC, Vrije Universiteit Amsterdam. For donor characteristics, see Additional file 1: Table S1.
Post-mortem in-situ MRI acquisition
Post-mortem in-situ (brain in cranium) MRI scans were acquired according to a previously described pipeline  (Fig. 1), with a post-mortem delay (interval between death and MRI) of maximum 13 h for all brain donors. Briefly, scans were acquired on a 3T scanner (Signa-MR750, General Electric Medical Systems, Milwaukee, WI) with an eight-channel phased-array head-coil. The following pulse-sequences were performed for all subjects: (i) sagittal 3D T1-weighted fast spoiled gradient echo (repetition time [TR] = 7 ms, echo time [TE] = 3 ms, flip angle = 15°, 1-mm-thick axial slices, in-plane resolution = 1.0 × 1.0 mm2); (ii) sagittal 3D fluid attenuation inversion recovery (FLAIR; TR = 8000 ms, TE = 130 ms, inversion time [TI] = 2000–2500 ms, 1.2-mm-thick axial slices, in-plane resolution = 1.11 × 1.11 mm2), with TI corrected for post-mortem delay; (iii) diffusion-weighted imaging (DWI) axial 2D echo-planar imaging with diffusion gradients applied in 30 non-collinear directions, TR/TE = 7400/92 ms, slices thickness of 2.0 mm, in-plane resolution = 2.0 × 2.0 mm2, b = 1000 s/mm2, and 5 b0 scans. To allow for offline distortion correction of the images, b0 images with reversed phase-encoding direction were acquired using the same sequence parameters.
Cortical thickness and brain volume assessment
To minimize the impact of age-related white matter abnormalities (e.g. vascular change) on automated segmentations, the 3D-T1 images were lesion-filled  as previously described [39, 40]. Image processing was performed using Freesurfer, version 6.0 (http://surfer.nmr.mgh.harvard.edu) . For each subject, nine cortical regions of interest (ROIs) in the right hemisphere were selected from the Lausanne atlas parcellation [42, 43], to closely match the location of the cortical regions dissected at autopsy. Details on the Lausanne atlas labels corresponding to the ROIs in the present study can be found in Additional file 1: Table S2, and included the superior frontal gyrus, anterior and posterior cingulate gyrus, anterior insula, middle temporal gyrus, superior parietal gyrus, entorhinal cortex, parahippocampal and fusiform gyrus. Cortical thickness was measured as the perpendicular distance from the grey/white matter boundary to the corresponding pial surface. Quality of parcellations was assessed using the ENIGMA guidelines for quality control (see http://enigma.ini.usc.edu/protocols/imaging-protocols/). Furthermore, post-mortem normalized brain volume, and normalized grey and white matter volumes were obtained from the T1-weighted images using Structural Image Evaluation, using Normalisation, of Atrophy (SIENAX) (part of FSL 5.0.9; http://fsl.fmrib.ox.ac.uk/), which estimates brain tissue volume normalized for skull size .
Diffusion tensor imaging (DTI) pre-processing and free water correction
DTI and reference scans of opposite phase encodings (acquired along anterior–posterior and posterior-anterior directions) were collected for all 9 PD donors, 11 out of 12 PDD/DLB donors, and 14 out of 18 control donors (see Additional file 1: Table S1). DTI was first corrected for susceptibility distribution and eddy current-induced geometric distortion using FSL (Eddy and topup)  and fitted for single tensors diffusion maps, deriving MD , which is a good (bio)marker for grey matter microstructure, where diffusion does not conform to a specific direction . To avoid underestimation by partial volume effects of the CSF, a bi-tensor model for free water correction was performed using an open-resource, python-based script of DIPY, which has been shown to be plausible for using single shell DWI acquisitions to derive free-water corrected MD (FWC-MD) maps [47, 48]. Using Freesurfer version 7.1.1 , we performed a within-subject registration of T1-weighted to FWC-MD images. With this registration, the cortical thickness of ROIs was registered to the FWC-MD map and corrected for overestimation into the white matter, by applying a threshold at 50% probability and a limit to the cortical ribbon mask, after which regional cortical FWC-MD values were obtained. Quality of parcellations was manually assessed (I.F.) and, when necessary, poorly parcellated regions were manually corrected or excluded from statistical analysis (i.e. entorhinal cortex of 1 PD case and 3 PDD/DLB cases was excluded). For clarity, FWC-MD is referred to as MD in the manuscript.
Subsequent to MRI acquisition, brain tissue was collected at autopsy, resulting in a total post-mortem delay (interval between death and autopsy) of maximum 13 h for all brain donors. Formalin-fixed paraffin-embedded tissue blocks (4% formalin, four weeks of fixation) from the following seven regions within the right hemisphere were obtained and processed for immunohistochemistry (IHC): superior frontal gyrus, anterior and posterior cingulate gyri, anterior (dysgranular) insula, middle temporal gyrus, superior parietal gyrus, and hippocampus (including the entorhinal cortex, and parahippocampal and fusiform gyri, as described before) [39, 49].
IHC for quantification
Sections (6-µm thick) from tissue blocks of the above-mentioned regions were cut and mounted onto superfrost + glass slides (Thermo Scientific, Waltham, MA). All sections were stained for pSer129-αSyn (clone EP1536Y), Aβ (clone 4G8) and p-tau (clone AT8). Additionally, sections from 7 out of the 9 cortical brain regions (showing significantly higher LB count across groups, i.e., all regions except the posterior cingulate cortex and the superior parietal gyrus) were stained for NfL (amino acid sequence 1–376) and p-NfM-H (clone SMI312) (for information of primary antibodies see Additional file 1: Table S3). Briefly, the sections were deparaffinised, immersed in buffer, and heated to 100 °C in a steam cooker for 30 min for antigen retrieval. The sections were blocked for endogenous peroxidase using 1% hydrogen peroxide and in Tris-buffered saline (TBS; pH 7.4), and consequently in 3% normal donkey serum in TBS (Triton 0.5%). Primary antibodies were diluted in 1% normal donkey serum in TBS (Triton 0.1%) and incubated overnight at 4 °C. Primary antibodies were detected using EnVision (Dako, Glostrup, Denmark), and visualized using 3.3’-Diaminobenzidine (DAB, Dako) with Imidazole (50 mg DAB, 350 mg Imidazole and 30 μl of H2O2 per 100 ml of Tris–HCl 30 mM, pH 7.6). Between steps, TBS was used to wash the sections. After counterstaining with haematoxylin, the sections were dehydrated and mounted with Entellan (Merck, Darmstadt, Germany).
Pathology image processing
Using a whole-slide scanner (Vectra Polaris, 20 × objective), images of immunostained sections were taken and quantified using QuPath 0.2.3 (https://qupath.readthedocs.io/en/0.2/) . ROIs containing all cortical layers were delineated in straight areas of the cortex to avoid over- or underestimation of pathology in sulci and gyri, respectively . Hippocampal sections were segmented according to the method described by Adler et al. , where the entorhinal cortex, parahippocampal gyrus and fusiform gyrus were delineated to match the MRI-derived ROIs. Briefly, the entorhinal cortex was delineated from the end of the parasubiculum until layer IV started to be visible ; the parahippocampal gyrus started from this point and ended at the collateral sulcus; the fusiform gyrus started at this point and ended at the inferior temporal sulcus (Additional file 2: Fig. S1). For the pSer129-αSyn, NfL and p-NfM/H stainings, the cortex was further segmented into superficial (layers I-III) and deep (layers IV-VI) cortical layers based on the haematoxylin counterstaining. For Aβ and p-tau, the same was done in the entorhinal cortex, parahippocampal and fusiform gyrus, as these are among the first regions to be affected by AD co-pathology [53, 54]. Note that the entorhinal cortex was subdivided into layers I-III (superficial) and lamina dissecans plus layers V-VI (deep) . DAB immunoreactivity was quantified with in-house QuPath scripts, using pixel and object classifiers. The outcome measures for pSer129-αSyn were both LB count per mm2 (LB/mm2) and pSer129-αSyn %area load excluding LBs (for additional information, see Additional file 3: Supplementary Material), while outcome measures for Aβ and p-tau were %area load (Additional file 2: Fig. S2). The outcome measures for neurofilaments were %area load, expressed in the text as %immunoreactivity (Additional file 2: Fig. S3).
Multi-label immunofluorescence and confocal microscopy for morphological characterisation of neurofilaments
To study the morphological features of neurofilaments, we performed a multi-label immunofluorescence staining in combination with 3D confocal laser scanning microscopy (CSLM) of p-NfM/H and NfL on a representative case of each group. This was done on the parahippocampal gyrus, as this was the region showing the largest neurofilament alterations between groups. Briefly, sections were deparaffinised, immersed in Tris EDTA pH 9.0, and heated to 100 °C in a steam cooker for 30 min for antigen retrieval. The sections were blocked for endogenous peroxidase using 1% hydrogen peroxide in TBS, and consequently in 3% normal donkey serum in TBS (Triton 0.5%; pH 7.4). Primary antibodies were diluted in 1% normal donkey serum in TBS (Triton 0.1%; pH 7.4) and incubated overnight at 4 °C. Primary antibodies were detected and visualized with donkey anti-mouse Alexa 488 (ThermoFisher, Pittsburgh, PA) and donkey anti-rabbit Alexa 594 (ThermoFisher) targeting p-NfM/H and NfL, respectively. For all stainings, TBS was used to wash the sections between steps. After counterstaining with DAPI (Sigma-Aldrich, St. Louis, MO), the sections were mounted with Mowiol (Sigma-Aldrich) plus anti-fading agent DABCO. Confocal imaging was performed with a Leica TCS SP8 (Leica Microsystems, Germany) using a HC PL PAO CS2 60 × oil objective lens, NA 1.40 and a zoom factor of 2.0. Sections were sequentially scanned for each fluorochrome with a pulsed white light laser at different wavelengths (excitation wavelengths: DAPI at 405 nm; Alexa 488 at 499 nm; Alexa 594 at 598 nm). All signals were detected using gated hybrid detectors in counting mode. Z-stacks (Z = 6 μm; 1024 × 1024 pixels) were taken in the parahippocampal gyrus of a representative control, PD and DLB. After scanning, the images were deconvoluted using CMLE algorithms in Huygens Professional (Scientific Volume imaging; Huygens, The Netherlands; https://svi.nl/Huygens-Professional), and their maximum projections (ImageJ Fiji, National Institute of Health, USA; https://imagej.nih.gov/ij/) or 3D surface rendering reconstructions (Imaris, Oxford instruments 2022: https://imaris.oxinst.com/) were used to represent graphically the structures of interests and their morphologies. In some cases, signal brightness was increased in ImageJ for clarity. Figures were created using Adobe Illustrator (CS6, Adobe Systems incorporated).
Statistical analyses were performed in SPSS 26.0 (Chicago, IL). Normality was tested, and demographics of PD, PDD/DLB and control groups were compared using parametric or non-parametric tests for continuous data, and Fisher exact test for categorical data. The associations between neurofilament immunoreactivity and pathological staging were calculated with Spearman’s correlation. To account for multiple (i.e., 9) brain regions within cases (i.e., nested data), the group differences across all regions and the associations between variables were assessed with multiple linear mixed models. In the case of pathological group differences, the pathological variable of interest was the dependent variable, and age and gender the main effects (i.e., covariates). In the case of MRI group differences, cortical thickness or MD was the dependent variable, and age, sex and post-mortem delay the main effects (i.e., covariates). In the case of global associations between neurofilament immunoreactivity and pathology load (for example NfL and p-tau), the neurofilament immunoreactivity was the dependent variable and the pathological load was the main effect, together with age and gender (i.e., covariates). And lastly, in the case of global associations between MRI measures and pathological markers (for example cortical thickness and NfL), cortical thickness or MD was the dependent variable and neurofilament immunoreactivity or pathology load was the main effect, together with age, gender and post-mortem delay (i.e. covariates). The intercept was included in all analyses as random effect. Post-mortem delay was added as a covariate in MRI-related analyses since it might influence post-mortem MRI-derived biomarkers [55, 56]. For linear mixed model associations, cortical thickness metrics were transformed into z-scores for each brain region as different cortical regions inherently have different thicknesses, which may drive the association. In all analyses, multiple comparisons between groups were corrected with post-hoc Bonferroni, and group comparisons were expressed as percentage differences, as in the formula: %difference = [(absolute difference)/mean control]*100. Statistics at the brain region level were corrected for multiple comparisons using the false discovery rate approach (FDR), and FDR-corrected P-values are expressed as P-FDR . After Bonferroni and FDR correction, P-values less than 0.05 were considered significant. Correlation coefficients (r) of linear mixed model associations were calculated as in the formula: r = (estimate fixed effect * standard deviation-fixed effect) / standard deviation-dependent variable.
Demographic, clinical, radiological, and pathological data of PD, PDD/DLB and non-neurological control donors are summarized in Table 1 (and per donor in Additional file 1: Table S1). Sex, age at diagnosis, age at death, and post-mortem delay did not differ between groups, whereas disease duration was shorter in PDD/DLB compared to PD donors (P = 0.006), due to the inclusion of DLB donors with shorter disease duration (n = 5, mean ± SD, 5 ± 1.6 years). On MRI, normalized brain volume (P = 0.011) and normalized grey matter volume (P = 0.019), but not normalized white matter volume (P = 0.095) were lower in PDD/DLB cases compared to controls. In terms of pathology load, the PDD/DLB group showed more abundant LB count, p-tau and Aβ load compared to controls (number of LB/mm2: P < 0.001; p-tau: P = 0.015; Aβ: P = 0.018) and PD donors (number of LB/mm2: P = 0.033; p-tau: P = 0.035; Aβ: P = 0.033). Moreover, LBs were more abundant in deep than superficial cortical layers (PD: P = 0.004; PDD/DLB: P = 0.035), while p-tau and Aβ were more abundant in superficial than deep cortical layers in PDD/DLB (p-tau: P = 0.046; Aβ: P < 0.001). The p-tau load strongly correlated with both LB count (r = 0.68, P < 0.001) and pSer129-αSyn load (r = 0.60, P < 0.001) in the PDD/DLB group. More details on regional and layer-specific distribution of pSer129-αSyn, Aβ, and p-tau load, are shown in Additional file 2: Fig. S4. For more details on correlations between pathological hallmarks, see Additional file 1: Table S4.
Increase in neurofilament immunoreactivity in the cortex of PD and PDD/DLB
P-NfM/H and NfL immunoreactivity were investigated in seven cortical areas. P-NfM/H staining was observed in axons, while NfL staining was seen in the neuronal soma and its processes (Fig. 2a, b). Both p-NfM/H and NfL immunoreactivity positively correlated with age (p-NfM/H: r = 0.33, R2 = 11%, P = 0.011; NfL in superficial cortical layers: r = 0.23, R2 = 5%, P = 0.020), and p-NfM/H was significantly higher in males compared to females (P = 0.007) (Additional file 2: Fig. S5). P-NfM/H immunoreactivity was higher in deep compared to superficial cortical layers in all groups (controls: + 38% compared to superficial layers, P < 0.001; PD: + 22%, P < 0.001; PDD/DLB: + 23%, P < 0.001), while NfL immunoreactivity showed the opposite, being more abundant in superficial than in deep cortical layers in all groups (controls: + 26% compared to deep layers, P < 0.001; PD: + 27%, P < 0.001; PDD/DLB: + 21%, P = 0.003). For more details on layer-specific distribution of neurofilaments, see Additional file 2: Fig. S6.
Overall, p-NfM/H immunoreactivity was significantly increased in both PD (P < 0.001) and PDD/DLB (P < 0.001) compared to controls by 40% and 32% respectively, and did not differ between PD and PDD/DLB (P = 1.000) (Fig. 2a, and Additional file 2: Fig. S3a for representative images). To investigate whether the increase in p-NfM/H immunoreactivity in the PDD/DLB group was driven by a specific clinical phenotype (PDD or DLB), we explored the difference in p-NfM/H immunoreactivity between these groups, but observed no differences. Regionally, p-NfM/H immunoreactivity was significantly increased in both PD and PDD/DLB groups compared to controls in the superior frontal gyrus (PD: + 36%, P-FDR = 0.009; PDD/DLB: + 27%, P-FDR = 0.036), entorhinal cortex (PD: + 52%, P-FDR = 0.001; PDD/DLB: + 38%, P-FDR = 0.003), parahippocampal gyrus (PD: + 54%, P-FDR < 0.001; PDD/DLB: + 46%, P-FDR < 0.001), fusiform gyrus (PD: + 42%, P-FDR = 0.019; PDD/DLB: + 39%, P-FDR = 0.005) and anterior insula (PD: + 46%, P-FDR = 0.003; PDD/DLB: + 29%, P-FDR = 0.049) (Fig. 2c).
Overall, NfL immunoreactivity was significantly increased in the PDD/DLB group compared to controls (+ 30%, P = 0.037), and tended to be increased in PDD/DLB compared to PD (P = 0.058), while it did not differ in PD compared to controls (P = 1.000) (Fig. 2b, and Additional file 2: Fig. S3b for representative images). To investigate whether the increase in NfL immunoreactivity in the PDD/DLB group was driven by a specific clinical phenotype (PDD or DLB), we explored the difference in NfL immunoreactivity between these groups, but observed no differences. Regionally, NfL immunoreactivity was increased specifically in the entorhinal cortex (+ 56%, P-FDR = 0.028) and the parahippocampal gyrus (+ 71%, P-FDR = 0.016) in PDD/DLB compared to controls (Fig. 2d).
Taken together, p-NfM/H immunoreactivity was uniformly increased across all cortical areas analysed in all patient groups, while increased NfL immunoreactivity was observed specifically in the entorhinal cortex and the parahippocampal gyrus in PDD/DLB.
Accumulation and fragmentation of neurofilaments in diseased neurons in PDD/DLB
To better understand the observed increase in p-NfM/H and NfL immunoreactivity in PD and PDD/DLB, we examined the morphological differences of both neurofilaments in the region that was most affected, the parahippocampal gyrus, using multi-label immunofluorescence and 3D CSLM.
First, NfL was much more present in neuronal somas and apical axons of PDD/DLB compared to PD and controls (Fig. 3a–c). Specifically, NfL seemed to accumulate in neurons showing neurodegenerative morphologies, such as nuclear fragmentation and swelling (Fig. 3c), and ballooning and corkscrew deformation of the soma (Fig. 3d). Second, accumulation and fragmentation of NfL staining patterns were observed in several axons, where axonal swellings intermitted axonal fragmentations in PD and, more abundantly, in PDD/DLB (Fig. 3e, g, h, and i), which were not observed in controls. Third, we also observed glial cells close to an apparently fragmented axon in PDD/DLB (Fig. 3c, e). Lastly, while p-NfM/H and NfL co-localized in most of the control axons, axons in the deep cortical layers of PD and PDD/DLB cases were mostly positive for p-NfM/H or NfL, rarely showing co-localization of the two neurofilaments (Fig. 3f–h). This was also supported by our quantitative data, showing an absence of significant correlation between the two neurofilaments in PD and PDD/DLB groups, which was present in controls (Additional file 2: Fig. S7).
Taken together, NfL seemed to accumulate in diseased neurons in PDD/DLB and to fragment axons in PD and, more abundantly, in PDD/DLB. Moreover, both neurofilaments showed less overlap in immunoreactivity profile in PD and PDD/DLB groups.
Cortical neurofilament immunoreactivity correlates with pathological staging
To better understand the relationship between the observed increases in p-NfM/H and NfL immunoreactivity in the parahippocampal gyrus (the most affected region) and pathological disease staging, we tested the correlation of neurofilament immunoreactivity with Braak LB stage, Braak NFT stage and Thal phase.
P-NfM/H immunoreactivity in the parahippocampal gyrus positively and strongly correlated with Braak LB stage (rs = 0.75, P < 0.001) and Braak NFT stage (rs = 0.59, P < 0.001), and moderately with Thal phase (rs = 0.44, P = 0.007) in the whole cohort (Additional file 2: Fig. S8). Similarly, NfL immunoreactivity in the parahippocampal gyrus positively and moderately correlated with Braak LB stage (rs = 0.44, P = 0.007), and strongly with Braak NFT stage (rs = 0.55, P < 0.001) and Thal phase (rs = 0.54, P = 0.001) (Additional file 2: Fig. S8).
Cortical neurofilament immunoreactivity correlates with pSer129-αSyn load, and more strongly with p-tau load in PDD/DLB
No correlations were found between p-NfM/H immunoreactivity and pSer129-αSyn pathology in any group across all regions (P > 0.05, Fig. 4a). On the other hand, p-NfM/H positively correlated with p-tau load in the PDD/DLB group across all regions (r = 0.29, R2 = 8%, P = 0.043), which was a trend in PD (r = 0.32, P = 0.066), and absent in controls (P = 0.149) (Fig. 4b). No correlation was found with Aβ load (P > 0.05). The associations were not specific to the superficial or deep cortical layers (for all correlations, see Additional file 1: Table S5).
NfL positively correlated with pSer129-αSyn pathology in PDD/DLB across all regions (LB count: r = 0.40, R2 = 16%, P = 0.002; pSer129-αSyn load: r = 0.36, R2 = 13%, P = 0.005), but not in the PD group (LB count: P = 0.392; pSer129-αSyn load: P = 0.910) (Fig. 4c). Moreover, NfL also positively correlated with p-tau load in PDD/DLB (r = 0.42, R2 = 18%, P < 0.001) and PD (r = 0.31, R2 = 10%, P = 0.014), but not in controls across all regions (P = 0.579) (Fig. 4d). When both pSer129-αSyn and p-tau load were considered, NfL correlated with p-tau load (r = 0.24, R2 = 6%, P = 0.030) but not pSer129-αSyn load (P = 0.195). Taken together, it seems that NfL correlated more strongly with p-tau than with pSer129-αSyn load. No correlation was found between NfL and Aβ load (P > 0.05). Additionally, the associations were not specific to superficial or deep cortical layers (for all correlations, see Additional file 1: Table S6).
Taken together, in PDD/DLB, higher p-NfM/H immunoreactivity correlated with higher p-tau load, while NfL immunoreactivity correlated with higher pSer129-αSyn load, but more strongly with higher p-tau load. None of the neurofilaments correlated with Aβ load.
Cortical neurofilament immunoreactivity is reflected by MRI cortical thickness in PD
To investigate whether MRI-derived cortical thickness is sensitive to regional neurofilament changes in the cortex of PD and PDD/DLB donors, we assessed whether regional neurofilament immunoreactivity correlated with cortical thickness across the included brain regions.
On MRI, neither PD (P = 1.000) nor PDD/DLB cases (P = 0.756) showed cortical thickness differences compared to controls across and within all cortical brain regions (all P-FDR > 0.05) (Additional file 1: Fig. S9a and b). Moreover, cortical thickness did not correlate with any pathological hallmark (P > 0.05) (Additional file 2: Fig. S10a and c, and Additional file 1: Table S7).
P-NfM/H immunoreactivity negatively correlated with cortical thickness in the full PD + PDD/DLB cohort across all cortical regions (r = − 0.23, R2 = 5%, P = 0.020). Particularly, the association was driven by the PD group (r = − 0.32, R2 = 10%, P = 0.047), and not significant in PDD/DLB (r = − 0.21, P = 0.106) (Fig. 5a). No correlations were found in the control group (P > 0.05), or for individual brain regions (P-FDR > 0.05).
Additionally, NfL immunoreactivity negatively correlated with cortical thickness in the full PD + PDD/DLB cohort across all cortical regions (r = − 0.22, R2 = 5%, P = 0.012). Again, the association was driven by PD (r = − 0.27, R2 = 7%, P = 0.025), and absent in PDD/DLB (r = − 0.15, P = 0.220) (Fig. 5c). No correlations were found in the control group (P > 0.05), or for individual brain regions (P-FDR > 0.05).
Taken together, cortical thickness was not sensitive to pathological load. However, both p-NfM/H and NfL immunoreactivity levels correlated with cortical thickness, especially in the PD group.
Cortical neurofilament immunoreactivity is reflected by cortical diffusivity in PDD/DLB
To investigate whether cortical MD is sensitive to regional neurofilament variation in the cortex of PD and PDD/DLB donors, we assessed whether regional neurofilament immunoreactivity correlated with MD across the cortical regions of interest.
Neither PD (P = 1.000) nor PDD/DLB donors (P = 0.487) showed overall cortical MD differences compared to controls across and within the regions analysed (P-FDR > 0.05) (Additional file 2: Fig. S9c and d). Overall, cortical MD positively correlated with pSer129-αSyn (r = 0.32, R2 = 10%, P = 0.002) and p-tau pathology (r = 0.36, R2 = 13%, P < 0.001) in PDD/DLB across all cortical regions, but not with Aβ load (P = 0.339) (Additional file 2: Fig. S10b and d, and Additional file 1: Table S8).
P-NfM/H immunoreactivity did not correlate with cortical MD in the full PD + PDD/DLB cohort across all cortical regions (P = 0.182). However, p-NfM/H immunoreactivity positively correlated with cortical MD in the PDD/DLB group (r = 0.28, R2 = 8%, P = 0.003) (Fig. 5b). No correlations were found in the control group (P = 0.430), or the PD group (P = 0.586), or for individual brain regions (P-FDR > 0.05).
Additionally, NfL immunoreactivity positively correlated with cortical MD in the full PD + PDD/DLB cohort across all cortical regions (r = 0.21, R2 = 4%, P = 0.024). Particularly, the association was present in the PDD/DLB group (r = 0.41, R2 = 17%, P = 0.003), but not in the PD (P = 0.184) or control group (P = 0.738) (Fig. 5d). No correlations were found for individual brain regions (P-FDR > 0.05).
Taken together, cortical MD is a sensitive marker for pSer129-αSyn and p-tau load. In addition, both p-NfM/H and NfL immunoreactivity correlated with increased cortical MD, especially in the PDD/DLB group.
Using post-mortem within-subject in-situ MRI and histopathology approaches , we investigated cortical regional neurofilament immunoreactivity and its relation to pathological hallmarks and MRI outcome measures in PD and PDD/DLB from a microscale to a mesoscale level. We found that p-NfM/H was increased in PD and PDD/DLB across almost all cortical regions, while NfL was increased specifically in the parahippocampal gyrus and entorhinal cortex of PDD/DLB donors. NfL accumulated in somas of diseased neurons showing neurodegenerative morphologies, and in axons showing fragmentation. Furthermore, we found that in PDD/DLB, p-NfM/H immunoreactivity levels positively correlated with p-tau load, whereas NfL positively correlated with pSer129-αSyn pathology and more strongly with p-tau load. Lastly, we showed that neurofilament immunoreactivity correlated with cortical thinning in PD and with higher cortical MD in PDD/DLB.
For p-NfM/H, we found that both PD and PDD/DLB donors had higher p-NfM/H immunoreactivity levels compared to controls, and that this increase was uniform across the analysed cortical areas. An increase in phosphorylation of neurofilament medium and heavy chains (NfM/H) has been described in several neurological disorders, such as AD, multiple sclerosis, amyotrophic lateral sclerosis and stroke , but it has never been investigated in the cerebral cortex of PD patients. Phosphorylation is assumed to be the dominant phosoform of both NfM and NfH, and under normal conditions axons show extensive neurofilament phosphorylation, while there is little or no phosphorylation in the neuronal soma and its dendrites [58, 59], which we also show in our study. Neurofilament phosphorylation is a highly regulated and complex process, and increased phosphorylation of subunits M (medium) and H (heavy) has been linked to several pathological processes [58, 59], even if its function in disease is still under debate. The increase in p-NfM/H immunoreactivity we found may indicate two processes: (i) an increase in the amount of axons or (ii) an increase in phosphorylation of NfM/H. While the first option is unlikely in neurodegenerative diseases, the second has been widely described and suggests that neurofilament hyper-phosphorylation might be linked to axonal transport impairment [58, 59]. Furthermore, a putative neurofilament hyper-phosphorylation might be a by-product of dysregulated kinase activity due to stress factors, such as chronic oxidant stress in neurodegenerative diseases . Cellular stress factors are associated with extensive phosphorylation of cytoskeletal elements through several proline-directed kinases, which are capable of phosphorylating neurofilaments and tau [59, 60], and have been shown to be dysregulated in AD and PD [60, 61]. Moreover, we found that both p-NfM/H and NfL correlated to p-tau load in PDD/DLB, suggestive of cytoskeletal alterations in axons in cortical brain regions [59, 60]. However, more data should be gathered by future research to investigate neurofilament hyper-phosphorylation in PD. Interestingly, we found that p-NfM/H immunoreactivity was higher in males compared to female donors. Males also show a higher prevalence of PD, especially between 50 and 59 years old . Since we show that PD(D) and DLB donors had higher p-NfM/H immunoreactivity, our data suggests that an innate increased p-NfM/H immunoreactivity in males might be related to the higher predominance of PD in this sex group. However, these results need to be taken with caution as the study was not designed to answer this specific research question. Taken together (Fig. 6), our data suggests that the NfM/H subunits might be hyper-phosphorylated across the cortex not only in PDD/DLB, but also in PD donors with long disease duration, suggesting that NfM/H hyper-phosphorylation might be a marker of axonal stress. However, this should be confirmed in future studies.
Furthermore, we found that PDD/DLB cases had higher cortical NfL immunoreactivity levels compared to age-matched controls, and that this increase was specific to the entorhinal cortex and parahippocampal gyrus. An increase in NfL immunoreactivity seems counterintuitive, as it is generally believed that the breakdown of axons leads to an increase of this marker in the extracellular fluid, and therefore a decrease of it in brain tissue [15, 63]. However, higher NfL immunoreactivity has been described before in an animal model and in the brain tissues of stroke patients [64, 65]. Specifically, it has been shown that the same NfL antibody used in this study, the polyclonal rabbit NfL (AA 1–284), not only targets NfL, but most importantly also targets its degradation products . Translating this information to our study, we recorded an increase of NfL immunoreactivity in PDD/DLB, which may result from increased NfL degradation products. In support of this, we also found that NfL accumulated in neuronal somas that showed patterns of neurodegeneration (e.g., nuclear material fragmentation, ballooned and corkscrew cell bodies), and NfL-positive axons often showed fragmentation, as described before in a mouse model of stroke [64, 65]. Taken together, increased NfL in parahippocampal regions seems to point towards evidence of neurofilament accumulation and fragmentation in PDD and DLB. These results are consistent with several CSF and plasma studies, which report that higher CSF or plasma NfL levels correlate with cognitive decline in PD [16,17,18, 20, 23]. The information that we add to these studies is regional, by showing that NfL immunoreactivity is specifically increased in the entorhinal and parahippocampal cortex, regions that are strongly involved in cognitive processes , and are the first to show p-tau accumulation in the adult and aged human brain . This process might therefore hint towards why NfL CSF and plasma levels are increased in several neurological disorders which are per definition p-tau positive, such as AD, frontotemporal dementia and multiple system atrophy [15, 16, 63]. Moreover, besides NfL correlation with p-tau load, we show that NfL immunoreactivity levels were also associated with pSer129-αSyn load. CSF NfL levels have been shown to positively correlate with CSF αSyn  and with decline in memory, attentional and executive functioning in PD . Additionally, we recently showed that neurofilaments accumulate and cluster around the core of Lewy bodies in post-mortem brain tissue of PD and PDD/DLB donors , illustrating a close relationship between NfL and pathological αSyn, suggesting a role for neurofilament in encapsulating toxic proteins. Taken together (Fig. 6), aggregation and fragmentation of NfL, and its association to p-tau and pSer129-αSyn pathology, suggest that NfL undergoes structural changes with increasing pathological burden, confirming its important role in the context of cognitive decline.
To investigate the possibility of measuring the above-mentioned pathological changes in clinical practice, we evaluated whether cortical pathology and neurofilament immunoreactivity would be reflected by MRI biomarkers of neurodegeneration, including cortical thickness and cortical microstructural integrity. We found that increased neurofilament immunoreactivity correlated with lower cortical thickness in PD and with higher cortical MD in PDD/DLB. This discrepant association in non-demented and demented cases may be explained by the fact that we found cortical MD to be a more sensitive marker than cortical thickness in picking up pSer129-αSyn and p-tau pathology, which was more abundantly present in PDD/DLB than in non-demented PD. Moreover, we found an association between neurofilament immunoreactivity and reduced cortical thickness in PD patients, in the absence of obvious cortical atrophy. This was to be expected, since cortical degeneration on MRI has been shown to be extensive but very subtle in PD , and only MRI studies with large sample sizes would be able to pick up such changes.
Our results partially go in line with previous literature, where serum NfL was shown to correlate with cortical thinning of several brain areas in de-novo PD, but also with increased cortical MD . Unfortunately, no similar studies have been carried out in PDD and DLB, making this the first study to investigate this relationship. Taken together (Fig. 6), we found that cortical thinning in PD, and lower microstructural integrity in PDD/DLB seem to reflect not only accumulation of neurofilaments (i.e., increased NfL), but also axonal stress (i.e., increased p-NfH/M).
The main strength of this study is the combination of MRI and gold-standard histological data of multiple cortical regions, which were collected from well-characterized donors. All donors had pathological confirmation of clinical diagnosis, as clinical-pathological discrepancies occur in about 15% of cases and might obscure in-vivo studies [68, 69]. Moreover, we both quantitatively and qualitatively studied regional cortical neurofilament changes in PD. NfL is widely used as a CSF or blood biomarker to detect axonal degeneration in patient cohorts, and with this study we hope to shed light on the regional and morphological changes in neurofilaments, linking neurodegenerative features from micro- to mesoscale and to the clinic. Therefore, results of this study provide evidence for the important roles of neurofilaments in neurodegeneration and their MRI signatures.
There are also some limitations. Neurofilament changes are not a process specific to PD(D) or DLB. In fact, increased CSF and plasma NfL levels have been described in several other neurological disorders compared to age-matched controls, and are considered a measure of neuroaxonal damage independent of casual pathways [15, 16]. In addition, this study focused on neurofilament changes of only 7 cortical regions and their correlations to cortical thickness, rather than subcortical neurofilament changes and their correlations with volume, even if differential neurofilament immunoreactivity might be found at this level. In fact, PD patients might show higher NfL accumulation in the substantia nigra rather than the cortex, and DLB patients might have differential NfL accumulation in the basal ganglia and the occipital cortex. Moreover, neurofilament changes only account for part of the variance of the damage seen on MRI and DTI (in our study up to 17%), and it is likely that other cellular and molecular components contribute to neurodegeneration, such as synaptic loss or glial changes. In addition, we did not find associations between MRI outcome measures and neurofilament immunoreactivity within the included brain regions (only across). This lack of regional sensitivity to neurofilament changes might be due to our small sample size, since large PD(D) cohorts are needed to pick up subtle changes in cortical thinning . Of note, even though post-mortem diffusivity measures were acquired in situ (brain in cranium), which is closer to the in vivo brain scanning than the ex vivo setting (brain extracted from the cranium), several factors (i.e., lower body temperature, tissue decomposition and partial volume effects of decreased CSF diffusivity, swelling, hypoxia) may influence the diffusivity properties [55, 56]. Therefore, even though we took into account post-mortem delay as a covariate in our analyses, this should be kept in mind in the comparisons between post-mortem and in vivo studies. Taken together, future research should therefore be performed to investigate the mechanisms behind the increase in neurofilaments and their association with MRI biomarkers in larger cohorts and in other neurological disorders showing high NfL levels and more pronounced regional atrophy, such as AD, frontotemporal dementia and multiple system atrophy [15, 16].
Taken together, we provide evidence for increased phosphorylation of NfM/H, suggesting axonal stress across the cortex of PD, PDD and DLB donors. Furthermore, we show that NfL immunoreactivity is increased in the parahippocampal gyrus and entorhinal cortex of PDD and DLB donors, reflecting aggregation and fragmentation of NfL in diseased neurons and axons, suggesting that NfL undergoes structural changes with increasing pathological burden in the context of cognitive decline. Importantly, we provide evidence for the relation between pathological burden and neurofilament levels, and demonstrate that such microscopic markers at least partly explain MRI markers that are associated with the neurodegenerative process.
Availability of data and materials
The data that support the findings of this study are available from the corresponding author, upon reasonable request.
Dementia with Lewy bodies
Diffusion tensor imaging
False discovery rate
Free-water corrected mean diffusivity
Neurofilament light chain
Parkinson’s disease dementia
Phosphorylated neurofilament medium and heavy chain
Phosphorylated Ser129 alpha synuclein
Postuma RB, Berg D, Stern M, Poewe W, Olanow CW, Oertel W, et al. MDS clinical diagnostic criteria for Parkinson’s disease. Mov Disord. 2015;30(12):1591–601.
Aarsland D, Zaccai J, Brayne C. A systematic review of prevalence studies of dementia in Parkinson’s disease. Mov Disord. 2005;20(10):1255–63.
Jellinger KA. Morphological basis of Parkinson disease-associated cognitive impairment: an update. J Neural Transm (Vienna). 2022;129(8):977–99.
E Sousa CS, Alarcão J, Martins IP, Ferreira JJ. Frequency of dementia in Parkinson's disease: a systematic review and meta-analysis. J Neurol Sci. 2022;432:120077.
Emre M, Aarsland D, Brown R, Burn DJ, Duyckaerts C, Mizuno Y, et al. Clinical diagnostic criteria for dementia associated with Parkinson’s disease. Mov Disord. 2007;22(12):1689–707.
McKeith IG, Boeve BF, Dickson DW, Halliday G, Taylor JP, Weintraub D, et al. Diagnosis and management of dementia with Lewy bodies: Fourth consensus report of the DLB Consortium. Neurology. 2017;89(1):88–100.
Yamada M, Komatsu J, Nakamura K, Sakai K, Samuraki-Yokohama M, Nakajima K, et al. Diagnostic criteria for dementia with lewy bodies: updates and future directions. J Mov Disord. 2020;13(1):1–10.
Jellinger KA. Neuropathological spectrum of synucleinopathies. Mov Disord. 2003;18(S6):2–12.
Irwin DJ, White MT, Toledo JB, Xie SX, Robinson JL, Van Deerlin V, et al. Neuropathologic substrates of Parkinson disease dementia. Ann Neurol. 2012;72(4):587–98.
Lippa C, Duda J, Grossman M, Hurtig H, Aarsland D, Boeve B, et al. DLB and PDD boundary issues: diagnosis, treatment, molecular pathology, and biomarkers. Neurology. 2007;68(11):812–9.
Hepp DH, Vergoossen DL, Huisman E, Lemstra AW, Bank NB, Berendse HW, et al. Distribution and load of amyloid-β pathology in Parkinson disease and dementia with Lewy bodies. J Neuropathol Exp Neurol. 2016;75(10):936–45.
McKeith IG, Dickson DW, Lowe J, Emre M, O’Brien JT, Feldman H, et al. Diagnosis and management of dementia with Lewy bodies: third report of the DLB Consortium. Neurology. 2005;65(12):1863–72.
Tsuboi Y, Uchikado H, Dickson DW. Neuropathology of Parkinson’s disease dementia and dementia with Lewy bodies with reference to striatal pathology. Parkinsonism Relat Disord. 2007;13:S221–4.
Halliday GM, Song YJC, Harding AJ. Striatal β-amyloid in dementia with Lewy bodies but not Parkinson’s disease. J Neural Transm (Vienna). 2011;118(5):713–9.
Khalil M, Teunissen CE, Otto M, Piehl F, Sormani MP, Gattringer T, et al. Neurofilaments as biomarkers in neurological disorders. Nat Rev Neurol. 2018;14(10):577–89.
Ashton NJ, Janelidze S, Al Khleifat A, Leuzy A, van der Ende EL, Karikari TK, et al. A multicentre validation study of the diagnostic value of plasma neurofilament light. Nat Commun. 2021;12(1):1–12.
Yong AC, Tan YJ, Ng EY, Lu Z, Ng SY, Chia NS, et al. Association between plasma neurofilament light chain levels and cognition in early Parkinson’s disease: Biomarkers (non-neuroimaging)/Plasma/Serum/Urine biomarkers. Alzheimers Dement. 2020;16: e040206.
Lerche S, Wurster I, Röben B, Zimmermann M, Machetanz G, Wiethoff S, et al. CSF NFL in a longitudinally assessed PD cohort: age effects and cognitive trajectories. Mov Disord. 2020;35(7):1138–44.
Mollenhauer B, Dakna M, Kruse N, Galasko D, Foroud T, Zetterberg H, et al. Validation of serum neurofilament light chain as a biomarker of Parkinson’s disease progression. Mov Disord. 2020;35(11):1999–2008.
Lin CH, Li CH, Yang KC, Lin FJ, Wu CC, Chieh JJ, et al. Blood NfL: a biomarker for disease severity and progression in Parkinson disease. Neurology. 2019;93(11):e1104–11.
Oosterveld LP, Verberk IM, Majbour NK, El-Agnaf OM, Weinstein HC, Berendse HW, et al. CSF or serum neurofilament light added to α-synuclein panel discriminates Parkinson’s from controls. Mov Disord. 2020;35(2):288–95.
Hall S, Surova Y, Öhrfelt A, Zetterberg H, Lindqvist D, Hansson O. CSF biomarkers and clinical progression of Parkinson disease. Neurology. 2015;84(1):57–63.
Bäckström D, Linder J, Mo SJ, Riklund K, Zetterberg H, Blennow K, et al. NfL as a biomarker for neurodegeneration and survival in Parkinson disease. Neurology. 2020;95(7):e827–38.
Oosterveld LP, Kuiper TI, Majbour NK, Verberk IM, van Dijk KD, Twisk JW, et al. CSF biomarkers reflecting protein pathology and axonal degeneration are associated with memory, attentional, and executive functioning in early-stage parkinson’ s disease. Int J Mol Sci. 2020;21(22):8519.
Sampedro F, Kulisevsky J. Intracortical surface-based MR diffusivity to investigate neurologic and psychiatric disorders: a review. J Neuroimaging. 2022;32(1):28–35.
Burton EJ, McKeith IG, Burn DJ, Williams ED, O’Brien JT. Cerebral atrophy in Parkinson’s disease with and without dementia: a comparison with Alzheimer’s disease, dementia with Lewy bodies and controls. Brain. 2004;127(4):791–800.
Beyer MK, Larsen JP, Aarsland D. Gray matter atrophy in Parkinson disease with dementia and dementia with Lewy bodies. Neurology. 2007;69(8):747–54.
Watson R, Blamire AM, Colloby SJ, Wood JS, Barber R, He J, et al. Characterizing dementia with Lewy bodies by means of diffusion tensor imaging. Neurology. 2012;79(9):906–14.
Atkinson-Clement C, Pinto S, Eusebio A, Coulon O. Diffusion tensor imaging in Parkinson’s disease: review and meta-analysis. Neuroimage. 2017;16:98–110.
Laansma MA, Bright JK, Al-Bachari S, Anderson TJ, Ard T, Assogna F, et al. International multicenter analysis of brain structure across clinical stages of Parkinson’s disease. Mov Disord. 2021;36(11):2583–94.
Sampedro F, Martínez-Horta S, Marín-Lahoz J, Pagonabarraga J, Kulisevsky J. Longitudinal intracortical diffusivity changes in de-novo Parkinson’s disease: a promising imaging biomarker. Parkinsonism Relat Disord. 2019;68:22–5.
Taylor KI, Sambataro F, Boess F, Bertolino A, Dukart J. Progressive decline in gray and white matter integrity in de novo Parkinson’s disease: an analysis of longitudinal Parkinson progression markers initiative diffusion tensor imaging data. Front Aging Neurosci. 2018;10:318.
Sampedro F, Pérez-González R, Martínez-Horta S, Marín-Lahoz J, Pagonabarraga J, Kulisevsky J. Serum neurofilament light chain levels reflect cortical neurodegeneration in de novo Parkinson’s disease. Parkinsonism Relat Disord. 2020;74:43–9.
Jonkman LE, Galis-de Graaf Y, Bulk M, Kaaij E, Pouwels PJ, Barkhof F, et al. Normal Aging Brain Collection Amsterdam (NABCA): A comprehensive collection of postmortem high-field imaging, neuropathological and morphometric datasets of non-neurological controls. Neuroimage. 2019;22:101698.
Alafuzoff I, Thal DR, Arzberger T, Bogdanovic N, Al-Sarraj S, Bodi I, et al. Assessment of β-amyloid deposits in human brain: a study of the BrainNet Europe Consortium. Acta Neuropathol. 2009;117(3):309–20.
Alafuzoff I, Arzberger T, Al-Sarraj S, Bodi I, Bogdanovic N, Braak H, et al. Staging of neurofibrillary pathology in Alzheimer’s disease: a study of the BrainNet Europe Consortium. Brain Pathol. 2008;18(4):484–96.
Alafuzoff I, Ince PG, Arzberger T, Al-Sarraj S, Bell J, Bodi I, et al. Staging/typing of Lewy body related α-synuclein pathology: a study of the BrainNet Europe Consortium. Acta Neuropathol. 2009;117(6):635–52.
Steenwijk MD, Pouwels PJ, Daams M, van Dalen JW, Caan MW, Richard E, et al. Accurate white matter lesion segmentation by k nearest neighbor classification with tissue type priors (kNN-TTPs). Neuroimage. 2013;3:462–9.
Frigerio I, Boon BD, Lin C-P, Galis-de Graaf Y, Bol J, Preziosa P, et al. Amyloid-β, p-tau and reactive microglia are pathological correlates of MRI cortical atrophy in Alzheimer’s disease. Brain Commun. 2021;3(4):fcab281.
Lin C-P, Frigerio I, Boon BD, Zhou Z, Rozemuller AJ, Bouwman F, et al. Structural (dys) connectivity associates with cholinergic cell density of the nucleus basalis of Meynert in Alzheimer’s disease. Brain. 2022;145(8):2869–81.
Dale AM, Fischl B, Sereno MI. Cortical surface-based analysis: I. Segmentation and surface reconstruction. Neuroimage. 1999;9(2):179–94.
Daducci A, Gerhard S, Griffa A, Lemkaddem A, Cammoun L, Gigandet X, et al. The connectome mapper: an open-source processing pipeline to map connectomes with MRI. PLoS ONE. 2012;7(12):e48121.
Hagmann P, Cammoun L, Gigandet X, Meuli R, Honey CJ, Wedeen VJ, et al. Mapping the structural core of human cerebral cortex. PLoS Biol. 2008;6(7):e159.
Jenkinson M, Beckmann CF, Behrens TE, Woolrich MW, Smith SM. Fsl. Neuroimage. 2012;62(2):782–90.
Cicchetti DV. Guidelines, criteria, and rules of thumb for evaluating normed and standardized assessment instruments in psychology. Psychol Assess. 1994;6(4):284.
Chiapponi C, Piras F, Piras F, Fagioli S, Caltagirone C, Spalletta G. Cortical grey matter and subcortical white matter brain microstructural changes in schizophrenia are localised and age independent: a case-control diffusion tensor imaging study. PLoS ONE. 2013;8(10):e75115.
Garyfallidis E, Brett M, Amirbekian B, Rokem A, Van Der Walt S, Descoteaux M, et al. Dipy, a library for the analysis of diffusion MRI data. Frontiers Neuroinform. 2014;8:8.
Hoy AR, Koay CG, Kecskemeti SR, Alexander AL. Optimization of a free water elimination two-compartment model for diffusion tensor imaging. Neuroimage. 2014;103:323–33.
Adler DH, Pluta J, Kadivar S, Craige C, Gee JC, Avants BB, et al. Histology-derived volumetric annotation of the human hippocampal subfields in postmortem MRI. Neuroimage. 2014;84:505–23.
Bankhead P, Loughrey MB, Fernandez JA, Dombrowski Y, McArt DG, Dunne PD, et al. QuPath: Open source software for digital pathology image analysis. Sci Rep. 2017;7(1):16878.
Arendt T, Morawski M, Gartner U, Frohlich N, Schulze F, Wohmann N, et al. Inhomogeneous distribution of Alzheimer pathology along the isocortical relief. Are cortical convolutions an Achilles heel of evolution? Brain Pathol. 2017;27(5):603–11.
Insausti R, Munoz-Lopez M, Insausti AM, Artacho-Perula E. The human periallocortex: layer pattern in presubiculum, parasubiculum and entorhinal cortex. A review. Front Neuroanat. 2017;11:84.
Braak H, Alafuzoff I, Arzberger T, Kretzschmar H, Del Tredici K. Staging of Alzheimer disease-associated neurofibrillary pathology using paraffin sections and immunocytochemistry. Acta Neuropathol. 2006;112(4):389–404.
Thal DR, Rüb U, Orantes M, Braak H. Phases of Aβ-deposition in the human brain and its relevance for the development of AD. Neurology. 2002;58(12):1791–800.
Braak H, Del Tredici K, Rub U, de Vos RA, Jansen Steur EN, Braak E. Staging of brain pathology related to sporadic Parkinson's disease. Neurobiol Aging. 2003;24:197-211 https://doi.org/10.1016/s0197-4580(02)00065-9.
Montine TJ, Phelps CH, Beach TG, Bigio EH, Cairns NJ, Dickson DW, Duyckaerts C, Frosch MP, Masliah E, Mirra SS. National Institute on Aging–Alzheimer’s Association guidelines for the neuropathologic assessment of Alzheimer’s disease: a practical approach. Acta neuropathologica. 2012;123:1–11.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57(1):289–300.
Petzold A. Neurofilament phosphoforms: surrogate markers for axonal injury, degeneration and loss. J Neurol Sci. 2005;233(1–2):183–98.
Holmgren A, Bouhy D, Timmerman V. Neurofilament phosphorylation and their proline-directed kinases in health and disease. J Peripher Nerv Syst. 2012;17(4):365–76.
Lee H-g, Perry G, Moreira PI, Garrett MR, Liu Q, Zhu X, et al. Tau phosphorylation in Alzheimer’s disease: pathogen or protector? Trends Mol Med. 2005;11(4):164–9.
Cheung ZH, Ip NY. Cdk5: a multifaceted kinase in neurodegenerative diseases. Trends Cell Biol. 2012;22(3):169–75.
Pringsheim T, Jette N, Frolkis A, Steeves TD. The prevalence of Parkinson’s disease: a systematic review and meta-analysis. Mov Disord. 2014;29(13):1583–90.
Ashton NJ, Leuzy A, Lim YM, Troakes C, Hortobagyi T, Hoglund K, et al. Increased plasma neurofilament light chain concentration correlates with severity of post-mortem neurofibrillary tangle pathology and neurodegeneration. Acta Neuropathol Commun. 2019;7(1):5.
Härtig W, Krueger M, Hofmann S, Preißler H, Märkel M, Frydrychowicz C, et al. Up-regulation of neurofilament light chains is associated with diminished immunoreactivities for MAP2 and tau after ischemic stroke in rodents and in a human case. J Chem Neuroanat. 2016;78:140–8.
Mages B, Aleithe S, Altmann S, Blietz A, Nitzsche B, Barthel H, et al. Impaired neurofilament integrity and neuronal morphology in different models of focal cerebral ischemia and human stroke tissue. Front Cell Neurosci. 2018;12:161.
Weintraub D, Doshi J, Koka D, Davatzikos C, Siderowf AD, Duda JE, et al. Neurodegeneration across stages of cognitive decline in Parkinson disease. Arch Neurol. 2011;68(12):1562–8.
Moors TE, Maat CA, Niedieker D, Mona D, Petersen D, Timmermans-Huisman E, et al. The subcellular arrangement of alpha-synuclein proteoforms in the Parkinson’s disease brain as revealed by multicolor STED microscopy. Acta Neuropathol. 2021;142(3):423–48.
Geut H, Hepp D, Foncke E, Berendse H, Rozemuller J, Huitinga I, et al. Neuropathological correlates of parkinsonian disorders in a large Dutch autopsy series. Acta Neuropathol Commun. 2020;8(1):1–14.
Hughes AJ, Daniel SE, Ben-Shlomo Y, Lees AJ. The accuracy of diagnosis of parkinsonian syndromes in a specialist movement disorder service. Brain. 2002;125(4):861–70.
We would like to thank all brain donors and their next of kin for brain donation. We would also like to thank the autopsy teams of the Netherlands Brain Bank (NBB) and Normal Aging Brain collection Amsterdam (NABCA).
This study was funded by The Michael J. Fox Foundation (grant #17253) and Stichting ParkinsonFonds (grant #1881). F.B. is supported by the NIHR biomedical research centre at the University College Hospital of London (UCLH). The authors have no relevant financial or non-financial interests to disclose.
Ethical approval and consent to participate.
All donors signed an informed consent for brain donation and the use of material and clinical information for research purposes. The procedures for brain tissue collection of NBB and NABCA have been approved by the Medical Ethical Committee of Amsterdam UMC, Vrije Universiteit Amsterdam.
Consent for publication
The authors declare that they have no competing interests.
. Supplementary tables. Table S1. Donor characteristics. Table S2. Concatenated labels from Lausanne atlas for the regions of interest of this study. Table S3. Information on primary and secondary antibodies. Table S4. Correlations between pathological hallmarks. Table S5. Correlations of p-NfM/H immunoreactivity with pathology load in superficial (layer I-III) and deep cortical layers (layer IV-VI). Table S6. Correlations of NfL immunoreactivity with pathology load in superficial (layer I-III) and deep cortical layers (layer IV-VI). Table S7. Correlations of MRI-derived cortical thickness with neurofilaments immunoreactivity and pathology load. Table S8. Correlations of DTI-derived cortical MD with neurofilaments immunoreactivity and pathology load.
. Supplementary figures. Fig. S1. Segmentation of entorhinal cortex, parahippocampal and fusiform gyrus in hippocampal section. Fig. S2. Quantification of pathology load using QuPath algorithms. Fig. S3. Photographs of representative sections of controls, PD and PDD/DLB donors stained for NfL and p-NfM/H. Fig. S4. Cortical pathology load in controls, PD and PDD/DLB donors. Fig. S5. Differences in cortical neurofilaments immunoreactivity in sex, and neurofilaments association with age. Fig. S6. Cortical neurofilaments immunoreactivity in superficial (layer I-III) and deep cortical layers (layer IV-VI). Fig. S7. Cortical neurofilaments correlation in superficial (layer I-III) and deep cortical layers (layer IV-VI). Fig. S8. Neurofilaments immunoreactivity in the parahippocampal gyrus correlates with pathological staging. Fig. S9. Cortical thickness and MD do not differ between controls, PD and PDD/DLB donors. Fig. S10. Cortical MD, but not cortical thickness, correlates with pSer129-αSyn and p-tau pathology.
. Supplementary materials. Pathological pSer129-αSyn quantification.
About this article
Cite this article
Frigerio, I., Laansma, M.A., Lin, CP. et al. Neurofilament light chain is increased in the parahippocampal cortex and associates with pathological hallmarks in Parkinson’s disease dementia. Transl Neurodegener 12, 3 (2023). https://doi.org/10.1186/s40035-022-00328-8