Next Article in Journal
Simple Does Not Mean Trivial: Behavior of Phosphatidic Acid in Lipid Mono- and Bilayers
Next Article in Special Issue
Neuroprotective Properties of Quinone Reductase 2 Inhibitor M-11, a 2-Mercaptobenzimidazole Derivative
Previous Article in Journal
Prediction of Anti-Glioblastoma Drug-Decorated Nanoparticle Delivery Systems Using Molecular Descriptors and Machine Learning
Previous Article in Special Issue
A Transcriptomic Meta-Analysis Shows Lipid Metabolism Dysregulation as an Early Pathological Mechanism in the Spinal Cord of SOD1 Mice
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Subnormal GM1 in PBMCs: Promise for Early Diagnosis of Parkinson’s Disease?

1
Department of Pharmacology, Physiology, and Neuroscience, Rutgers, The State University of NJ, Newark, NJ 07103, USA
2
Department of Neurology, Columbia University Irving Medical Center, New York, NY 10032, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally.
Int. J. Mol. Sci. 2021, 22(21), 11522; https://doi.org/10.3390/ijms222111522
Submission received: 25 June 2021 / Revised: 7 October 2021 / Accepted: 15 October 2021 / Published: 26 October 2021

Abstract

:
The fact that Parkinson’s disease (PD) pathologies are well advanced in most PD patients by the time of clinical elucidation attests to the importance of early diagnosis. Our attempt to achieve this has capitalized on our previous finding that GM1 ganglioside is expressed at subnormal levels in virtually all tissues of sporadic PD (sPD) patients including blood cells. GM1 is present in most vertebrate cells, is especially abundant in neurons where it was shown essential for their effective functioning and long term viability. We have utilized peripheral blood mononuclear cells (PBMCs) which, despite their low GM1, we found to be significantly lower in sPD patients compared to age-matched healthy controls. To quantify GM1 (and GD1a) we used high performance thin-layer chromatography combined with cholera toxin B linked to horseradish peroxidase, followed by densitometric quantification. GM1 was also deficient in PBMCs from PD patients with mutations in the glucocerebrosidase gene (PD-GBA), apparently even lower than in sPD. Reasons are given why we believe these results obtained with patients manifesting fully developed PD will apply as well to PD patients in preclinical stages—a topic for future study. We also suggest that these findings point to a potential disease altering therapy for PD once the early diagnosis is established.

1. Introduction

Parkinson’s disease (PD) has been termed the fastest growing neurological disorder in the world [1], progressing in synchrony with the increasing age span of individuals worldwide. The absence of treatment that would affect disease progression attests to the importance of early diagnosis. It has become evident that the low efficacy of treating PD patients is due at least in part to the late diagnosis and start of therapy. By the time of clear clinical diagnosis of PD half or more of the dopaminergic (DA) neurons of the nigrostriatal pathway have suffered viability decline or complete loss. Those DA neurons have received primary attention owing to their direct involvement in the movement disorders that constitute the clinical hallmarks of sporadic PD (sPD). However, the neuropathology of sPD is now recognized as much broader in scope, major cell loss also occurring in other brain regions and cell types. A prime example is the basal forebrain where the cholinergic neurons start to degenerate early in sPD and gradually worsen with the appearance of dementia [2]. The noradrenergic neurons in the locus coeruleus are among the first to suffer degenerative loss in sPD [3] while serotonergic neurons are also involved [4]. Moreover, neurons throughout the body are also affected before, during, and after the motor dysfunction sets in. Treatment of these diverse populations with their associated symptoms would clearly be more effective if commenced early enough to prevent the extensive neuronal losses destined to occur, especially if accompanied by disease-altering therapy that would terminate progressive pathological changes.
Numerous efforts have been made and are currently in progress to develop such early diagnosis based on a variety of behavioral and imaging methodologies as well as biochemical markers [5]. Some of these have proved to be reliable indicators of preclinical sPD symptoms but what is additionally needed is the connection of the indicator to appropriate therapy to effectively terminate disease progression at the early stage. We present here our findings on peripheral blood mononuclear cells (PBMCs) and their depressed level of GM1 ganglioside that shows promise as a reliable biochemical marker while also pointing to potential disease-altering therapy. GM1 along with its oligosialo-analogs occurs in numerous mammalian cell types and is especially abundant in neurons where it mediates a wide variety of regulatory and neuroprotective functions essential for effective signaling and long-term viability [6,7] (see Discussion). Our findings suggesting a systemic deficiency of GM1 in the brain and peripheral tissues of sPD patients point to this as a major contributor to PD pathogenesis and GM1 replacement as a promising therapy for disease alteration; some evidence that this is the case is already at hand (see Discussion) [8]. Here we show that PBMCs despite their low level of GM1 expression, show GM1 deficiency compared to healthy controls and also the subnormal levels characteristic of sPD tissues in general. The PBMCs employed in this study were obtained from PD patients well advanced in the disease but reasons are given why we consider it likely the same will be found in patients with pre-clinical PD (see Discussion).
In diagnosis and treatment of PD, it is beneficial to distinguish those patients with mutations in the glucocerebrosidase (GBA) gene from those with sporadic PD lacking such mutations; such mutations are found in 5–20% of PD cases and their symptoms and treatment options differ somewhat [9]. GBA encodes the lysosomal enzyme glucosylceramide hydrolase which showed lower activity in GBA homozygotes/compound heterozygotes and also GBA heterozygotes compared to non-carriers [10]. Mutations in this gene are the most significant genetic risk factors for sporadic PD. The clinical course of PD-GBA has been shown to be more aggressive with an earlier age of onset and a lower median survival time than sporadic PD [9]. The present study using a revised and simplified procedure shows GM1 levels in PBMCs from such PD-GBA patients to be subnormal like PD (including genetic cases) in general and even lower than in sPD.

2. Results

2.1. Clinical Characteristics of sPD and PD-GBA Patients and Healthy Control Subjects

Lipid analysis was completed for a total of 23 sPD patients, 13 PD-GBA patients, and 21 age-matched healthy controls, with males and females in each group. The demographical information and clinical assessment of all subjects are shown in Table 1 and Table 2. PBMCs from sPD and PD-GBA patients and controls were analyzed in two separate cohorts (as identified in Table 1 and Table 2).

2.2. Original Method: Analysis of GM1 and GD1a in PBMCs

In our first procedure GM1 ganglioside from PBMCs was separated via high performance thin-layer chromatography (HPTLC) and detected with cholera toxin B (CtxB) subunit linked to horseradish peroxidase (HRP) (Figure 1).
The high sensitivity of CtxB-HRP permitted the detection of the limited levels of GM1 in these cells. GM1 and other ganglio-series gangliosides were previously reported to be absent (as major gangliosides) in human monocytes/macrophages; that study employed resorcinol spray for detection on thin-layer chromatography [11]. The latter report mentions the conflicting reports on the types of gangliosides present in such cells. Typical PBMCs preparations are known to include a number of cell types including T cells, and our previous studies demonstrated these to clearly express GM1 and GD1a [12].
Because of its high selectivity for GM1, the CtxB reagent could not detect other gangliosides initially, but three other members of the ganglio-series were revealed following desialylation by neuraminidase (N’ase) on the HPTLC plate (see Methods). Neuraminidase catalyzes the removal of all sialic acid attached to ganglio-series ganglioside except GM1, which is why it is the only ganglioside left to react with CtxB. In this process neuraminidase converts GD1a, GD1b, GT1a, and GT1b to GM1 and the GM1 produced by this procedure is detected by CtxB. We have confirmed that overnight treatment with neuraminidase to complete the process, which can be visualized on a HPTLC plate with CtxB reagent [13]. These too are shown in Figure 1. GD1a is one of the primary interests because of its function as a metabolic precursor and reserve pool for GM1. Quantification of GM1 and GD1a with the aid of ImageJ densitometric analysis (see Methods) revealed significantly less of both gangliosides in PBMCs from sporadic PD and PD-GBA patients, compared to age-matched controls. Importantly, there was minimal overlap between controls and either form of PD (Figure 1). On the other hand, while PBMCs from PD-GBA patients expressed lower GM1 on average than PBMCs from sPD patients this difference did not prove significant. From this, we conclude that the GM1 level in PBMCs constitutes a reliable biomarker for both forms of PD but can only suggestively distinguish PD-GBA from sPD—at least by this method (see below for revised procedure). We were able to achieve this distinction with the revised method described below.

2.3. Modified Method: Analysis of GM1 in PBMCs

Because this was a rather lengthy and involved procedure, we next attempted to modify/shorten it by omitting the N’ase application. In that case, GM1 was the only ganglioside revealed by CtxB-HRP (Figure 2). In addition to omitting N’ase, this modified procedure contained a few other relatively minor changes (see Materials and Methods).
Quantification of the available samples (six each of sPD and PD-GBA) also revealed significantly lower GM1 in both sPD and PD-GBA, and fortuitously, a significantly lower level in PD-GBA than sPD (Figure 2). Accordingly, we conclude that this modified/shortened method provides a differentiated diagnosis of sPD and PD-GBA.

2.4. ROC Analysis of GM1 and GD1a in the PMBCs in sPD and PD-GBA Patients and Healthy Controls

2.4.1. Original Analytical Method

Receiver operating characteristics (ROC) curves were performed to analyze the diagnostic accuracy of using ganglioside levels detectable in PBMCs of healthy controls, sPD and PD-GBA patients. The area under the curve (AUC) was statistically significant (p < 0.05) for GM1 when comparing healthy controls with sPD (AUC = 0.97) and PD-GBA (AUC = 1.0) patients; and GD1a when compared with healthy controls (sPD: AUC = 0.98; PD-GBA: AUC = 1.0) (Figure 3). The determination of the sensitivity and specificity of GM1 as a biomarker revealed both high sensitivity (100%) and specificity (>85.7%) in distinguishing sPD and PD-GBA patients from controls with an optimal cutoff value (calculated using Youden’s index) of 397.5 and 324.7 µg/mg protein respectively. These cutoff values are below the lowest GM1 amount from healthy controls, indicating that this assay can distinguish between samples taken from sPD and PD-GBA patients and healthy controls with high probability. The optimal cutoff value for differentiating sPD patients from PD-GBA patients using GM1 levels was 284.5 µg/mg protein with high sensitivity (100%). The sensitivity correctly identifies the sample as originating from a patient with the disease. However, given that the specificity rate was 55.6%, this indicates that the probability of incorrectly identifying a sPD sample as a PD-GBA sample is near chance using this analytical method.

2.4.2. Modified Analytical Method

The accuracy of the modified analytical method was evident when comparing both patient groups with healthy controls (sPD: AUC = 0.98, p < 0.0001; PD-GBA: AUC = 1.0; p = 0.0005) (Figure 4). GM1 levels in PBMCs can also be used to distinguish sPD from PD-GBA patients using the modified analytical method with better accuracy (original method: AUC = 0.79, n.s.; modified method: AUC = 0.93, p = 0.003). Additionally, at the optimal cutoff value of <27.62 ng, the sensitivity (>83.3%) and specificity (92.9%) of the modified analytical method can distinguish PD-GBA patients from sPD patients’ above random chance.
The specificity of the modified analytical method (92.9%) was higher than the original analytical method (55.6%) indicating a smaller chance of incorrectly identifying a sPD patient as a PD-GBA patient, further supporting the advantage of using this modified method of GM1 measurement.

3. Discussion

This study revealed subnormal expression of GM1 ganglioside in PBMCs to constitute a reliable indicator of both sPD and PD-GBA. Ganglioside GD1a, revealed by N’ase application to the HPTLC plate, was also significantly deficient in sPD and PD-GBA. This was validated by assessing the diagnostic accuracy of GM1 and GD1a measurements in PBMCs using ROC curves. The high AUC values demonstrate the accuracy of these molecules for the differentiation of sPD and PD-GBA patients from controls. This was similarly robust when using our modified analytical method with simplified procedure that omitted N’ase. Additionally, the latter method proved able to distinguish between sPD and PD-GBA patients with high sensitivity and specificity. The modified analytical method had higher specificity in differentiating between sPD and PD-GBA patients than the original analytical method.
Because this study employed PBMCs from PD patients well advanced in clinical symptoms, additional work is required to substantiate this as a reliable method for diagnosis in the preclinical stage. If successful it would offer the additional advantage of also pointing to a potential disease altering treatment for PD to follow such diagnosis; this refers to GM1 replacement therapy to alleviate the neuronal deficiency of GM1 which our studies point to as the underlying cause of sPD [14].
GM1 ganglioside has been identified as a key player in the pathogenesis of PD due to its multiple regulatory and modulatory functions within the neuron [6,7]. These include—of special relevance to PD—GM1 facilitation of neuroprotection by glial cell line derived neurotrophic factor (GDNF) whose receptor complex includes GM1 as essential constituent [15]. GDNF was shown essential for the survival of adult catecholaminergic neurons, such as those of the substantia nigra [16]. It is noteworthy that NGF and BDNF, two other essential neurotrophic factors, also have receptors (TrkA, Trk B) with strongly associated GM1 [17,18]. This points to the additional neuronal types dependent on an adequate level of neuronal GM1 and is consistent with the broad array of neurons that dysfunction and gradually die, albeit at different rates, during PD progression [8]. Another GM1 function of major relevance to PD is its association with alpha-synuclein (aSyn), whose aggregation (e.g., Lewy bodies/Lewy neurites) is generally associated with and widely considered a major contributor to neuronal loss in PD. This protein was shown to bind GM1 with high affinity and specificity in a manner that maintains it in a non-aggregating conformation [19,20]. Application of GM1 to mice expressing such pathological aSyn was shown an effective means of dispersing such aggregates [21,22].
In view of these and additional vital functions of GM1 referred to above it is perhaps not surprising that even a partial deficiency of GM1 leads to gradual neuronal dysfunction and neuron death throughout the brain and body. That this would apply to virtually all neurons, albeit at different rates, would explain the enormous complexity and systemic nature of PD pathologies. It is becoming well recognized that PD does not involve selective damage to the dopaminergic neurons of the mesencephalon but is rather a severe multisystemic disorder [23]; this is well characterized, for example, in the numerous non-dopamine lesions that occur in PD [24]. Nigrostriatal neurons and those mediating the various prodromal symptoms would appear to be the most vulnerable and hence the earliest to succumb. That this is the underlying cause of sPD is supported by the mouse PD model based on mono-allelic disruption of the B4galnt1 gene which, in the process of diminishing GM2, also reduces GM1 to subnormal levels. Such mice manifest virtually all the neuropathology of PD in addition to motor dysfunctions [15,21]. The partial GM1 deficiency in such mice, caused without the use of a neurotoxin, approximates that subsequently found in human PD tissues [8]. What is particularly noteworthy is that such GM1 deficiency includes non-nervous tissues as well as neurons of the CNS and PNS, as demonstrated in the present study. PD is thus a truly systemic disorder.
In view of this model, it was speculated that human PD may have a similar genetic cause. In fact, mutated B4galnt1 has been demonstrated in patients with hereditary spastic paraplegia [25]. However, genome-wide association studies showed no evidence that this or other genes involved in GM1 synthesis are affected in sPD [26]. On the other hand, in situ hybridization histochemistry demonstrated a significant decrease in B3galt4 (GM1 synthase) and ST3gal2 (GD1a synthase) gene expression in residual neuromelanin-containing cells in the substantia nigra of PD patients; these decreases in GM1-synthesizing enzymes were cell-type specific and were not observed in non-neuromelanin containing neurons outside the substantia nigra in sPD brain [27]. Epigenetic effects, of course, are possible, such factors having been shown capable of influencing ganglioside biosynthesis [28,29]. We consider it likely a major cause of GM1 deficiency in PD derives from the aging process itself, as suggested in two key discoveries by Svennerholm’s group: (a) brain GM1 and GD1a levels vary significantly among individuals of the same age, and (b) gangliosides GM1 and GD1a both decline progressively with age [30]. This decrease was most striking for GD1a [31], the above mentioned metabolic precursor of GM1 which functions as a reserve pool for the latter through co-association with NEU3 neuraminidase [32]. We have speculated [14] that individuals whose endowed levels of GM1/GD1a fall in the very low normal range could start life with adequate levels of these two gangliosides but cross a threshold in later years when both GM1 and GD1a fall below the level necessary to maintain normal function of GDNF, aSyn and the other GM1-dependent processes. Such gradual age-related decline makes it likely that the subthreshold deficiency of GM1 characteristic of sPD would also be prevalent in the preclinical stages. An effort is underway to test this.
The promising preclinical experiments with GM1 applied to animal PD models [14] led to equally promising clinical trials with PD patients by Schneider and coworkers. This included a five-year open-label study [33] in which patients were administered an initial loading dose of 1000 mg GM1 (intravenous) followed by two subcutaneous injections of 100 mg GM1 per day. This resulted in reduced, slow linear projection of motor dysfunction based on UPDRS motor scores, and improved activities of daily living scores. It was significant that after five years the motor disability scores were lower than at baseline. This was followed by a randomized, controlled, delayed-start phase II trial that demonstrated GM1 superiority to placebo in reducing motor impairment and slowing symptom progression over a two-year period; UPDRS scores were improved at the trial end of 120 weeks compared to baseline [34]. This suggested to the investigators that “GM1 may have symptomatic and potentially disease modifying effects”. The above studies employed GM1 isolated from bovine brain and such studies were regrettably terminated due to FDA-perceived danger of contamination by prion-like proteins. Therapeutic use of E. coli derived GM1 proved equally effective as brain derived GM1 in a mouse PD model [35] but that too had to be discontinued due to manufacturing problems.
These examples suggest that PD origin and pathologies are causally related to the GM1 deficiencies in PD patient tissues. The present study indicates that GM1 in PBMCs serves well as a vehicle for the detection of such deficiency. PBMCs from PD patients were previously used to show a reduction of intracellular DA [36] and DA transporter [37]. What remains to be seen is whether such deficiency also occurs in preclinical PD. We have indicated our reasons for believing that such is the case, and that belief is strengthened by a recent report showing that GM1 is reduced in the serum of patients with REM sleep behavior disorder, a strong predictor of PD [38]. The latter study also demonstrated subnormal GM1 in substantia nigra, CSF, and serum. Beyond that is support by these findings that use of GM1-replacement therapy would prove truly disease modifying in a manner that would stop pathological progression. That the simplified procedure described herein also serves to distinguish between sPD and PD-GBA was an interesting additional discovery of this study.

4. Materials and Methods

4.1. Sample Collection and Patient Assessment

The PBMCs samples were provided for biochemical analysis by Roy N. Alcalay, these samples were collected during a visit to the Columbia University Irving Medical Center (NY, USA. Patients and healthy controls were assessed the same day for the presence and severity of Parkinsonian symptoms using The Unified Parkinson’s Disease Rating Scale (UPDRS). The following UPDRS exams were conducted: (1) UPDRS Part I for behavior and mood; (2) UPDRS Part II for activities of daily living; (3) UPDRS Part III motor examination.

4.2. PBMCs Isolation

PBMCs were isolated from the three populations mentioned above: (a) sPD, (b) PD-GBA, (c) healthy, age-matched controls—as described [39]. All PD-GBA samples were heterozygotes except one which was a homozygotic carrier. Whole blood samples were collected in PBS-diluted heparin and transferred into a Ficoll-filled Leucosep tube. This was centrifuged and the upper phase extracted and centrifuged again. The supernatant was decanted, and the cells washed. Cell counts were measured prior to resuspension and final centrifugation. The cell pellet was stored at −80 °C until shipment.

4.3. GM1 and GD1a Analysis

To extract total lipids, mixtures of 3.2 × 106–3.8 × 106 PBMC were treated with 2 mL of chloroform/methanol/H2O (C/M/W; 5/5/1) in a glass tube and homogenized 0.5 h with cell homogenizer and water bath sonicator. The mixture was then centrifuged at 5000 RPM for 15 min and the upper phase containing the lipids was removed from the protein-containing residue and transferred to another glass tube. The lipid solution was evaporated to dryness with a nitrogen stream and the protein quantified by a modified method of Lowry [40]; the latter results were used to determine the amount of lipid to apply to the HPTLC plate.
For HPTLC the lipid mixtures were dissolved in 55–65 µL C/M (1/2), adjusted to a volume that would allow application of 20 µL to the HPTLC plate corresponding to 100 µg protein. The Silica Gel 60 aluminum plates (Millipore Sigma, Burlington, USA; 20 × 20 cm) were placed in the paper-lined glass tank containing 150 mL C/M/0.25% aq. KCl (5/4/1) and the solvent allowed to ascend the plate for 1.5 h. The plate was dried for at least one hour and immersed in 1% polyisobutyl methacrylate (Sigma Aldrich, St. Louis, MO, USA; ave. MW 70,000—crystalline) for one minute. The plate was dried, then immersed in acetic acid buffer, pH 5.4 for one hour. The plate was immersed overnight with mild shaking in 20 mL of the same buffer containing 15 units of neuraminidase (N’ase; Clostridium—Sigma Aldrich). The plate was then washed twice for two min in phosphate-buffered saline (PBS). For blocking, the plate was immersed one hour in 15 mL PBS containing 6.7% powdered milk with mild shaking. For GM1 detection the plate was immersed one hour in 20 mL PBS containing 10 µg CtxB-HRP (Thermofisher, NY, USA), then washed 3× for 2–3 min in PBS. The plate was moved to a dark room and covered with five mL ECL Western Blotting Detection Reagent for one minute. The reagent was removed, and the plate was placed in a Hypercassette with HyBlot ES autoradiographic film and developed for five minutes. Quantification of the visible GM1 bands was carried out with ImageJ [41]. The quantities of GM1 and GD1a were determined using a ratiometric method in which the percentage change in optical density between the standards and samples was used to calculate the amount of GM1 and GD1a.
In the above procedure it was not necessary to first remove the other lipids from the gangliosides prior to HPTLC because they did not interfere due to non-reaction with CtxB; also, they generally moved well ahead of the gangliosides on the HPTLC plate.

4.4. Modified Analytical Method for Analysis of GM1 Alone

In this procedure we normalized the lipid quantity on the basis of cell number rather than protein content to determine the amount of lipid extract to apply to the HPTLC plate. As before, phospho- and other lipids ran well ahead of gangliosides and did not interfere. Following extraction of the lipids as described above and evaporation of the extract to dryness, C/M (1/2) was added (50–70 µL) in a calculated amount corresponding to 1.2 × 106 PBMCs. In addition to omission of N’ase, other changes from the above procedure included CtxB-HRP treatment for 1.5 h and development time with HyBlot ES autoradiographic film for six minutes.

4.5. Statistical Analysis

The UPDRS Part I and Part II scores were combined and analyzed together. The Total UPDRS scores were also reported and is the sum of the scores from Part I, II, and III exams. For the original analytical procedure—data are shown as the average amount and range of GM1 and GD1a in PBMCs (µg ganglioside/mg protein). For the modified analytical procedure—data are shown as the average amount and range of GM1 in PBMCs (ng per 1.2 × 106 PBMCs). One-way ANOVA with Tukey’s HSD post-hoc testing was used to analyze differences in ganglioside levels between controls and patient groups. The diagnostic accuracy of using ganglioside levels in PBMCs as a biomarker in sPD and PD-GBA was assessed using the receiver operating characteristic (ROC). This was completed using individual ROC curves for 3 comparisons for each biomarker (GM1 and GD1a): (1) control and sPD; (2) control and PD-GBA; and (3) sPD and PD-GBA. The area under the curve (AUC) and the optimal cutoff values for sensitivity and specificity (calculated using the Youden’s index) were further used to determine the diagnostic accuracy of these biomarkers [42,43]. All statistical analyses were completed using GraphPad Prism 9.0 (GraphPad Software, San Diego, CA, USA).

5. Conclusions

In this study we demonstrate that PBMCs from both sPD and PD-GBA patients express subnormal levels of GM1 and GD1a, while a modified/simplified procedure accomplishes that and also differentiates PD-GBA from sPD by demonstrating the lowest levels of GM1 in the former. Accordingly, this constitutes a promising method for early diagnosis of both forms of PD. It is clear that the benefits of early diagnosis would be significantly enhanced if accompanied by disease-altering therapy—a treatment that stops the progression of the pathology and even revives those neurons that have only become dormant/quiescent prior to final demise. We discuss the evidence suggesting that our procedure for early diagnosis points to such disease-altering therapy.

Author Contributions

S.K.A. and M.C., analytical procedures, data analysis, manuscript preparation; S.C., data analysis, manuscript preparation; R.N.A. and M.S., funding for PBMC acquisition (from the M.J. Fox Foundation), provided patient clinical data, manuscript preparation; R.L., project conceptualization and organization, work supervision, manuscript preparation. All authors have read and agreed to the published version of the manuscript.

Funding

Sample collections were funded to Roy N. Alcalay by the Michael J Fox Foundation and the Parkinson’s Foundation. Robert Ledeen received intra-departmental funding.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki, and approved by the Institutional Review Board of Columbia University protocol code IRB-AAAP9604.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to ethical reasons.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

aSynAlpha-synuclein
AUCArea under the ROC Curve
B3galt4Beta-1,3-Galatosyl transferase 4; GM1 Synthase
B4galnt1Beta-1,4-N-Acetyl-Galactosaminyltransferase 1; GM2 Synthase
BDNFBrain-derived neurotrophic factor
CNSCentral nervous system
CSFCerebrospinal Fluid
CtxBCholera toxin B
DADopamine; Dopaminergic neurons
GBAGlucocerebrosidase
GDNFGlial cell line derived neurotrophic factor
HPTLCHigh performance thin-layer chromatography
HRPHorseradish peroxidase
HSDHonestly significant difference
N’aseNeuraminidase
NEU3Neuraminidase 3
NGFNerve growth factor
PBMCPeripheral blood mononuclear cells
PBSPhosphate-buffered saline
PDParkinson’s disease
PNSPeripheral nervous system
REMRapid eye movement
ROCReceiver operating characteristic
RPMRevolutions per minute
sPDSporadic Parkinson’s disease
ST3gal2Beta-galactoside alpha-2,3-sialyltransferase 2; GD1a Synthase
TrkATropomyosin receptor kinase A
TrkBTropomyosin receptor kinase B
UPDRSUnified Parkinson’s disease rating scale

References

  1. Dorsey, E.R.; Sherer, T.; Okun, M.S.; Bloem, B.R. The Emerging Evidence of the Parkinson Pandemic. J. Parkinsons Dis. 2018, 8, S3–S8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Bohnen, N.I.; Albin, R.L. The cholinergic system and Parkinson disease. Behav. Brain Res. 2011, 221, 564–573. [Google Scholar] [CrossRef] [Green Version]
  3. Paredes-Rodriguez, E.; Vegas-Suarez, S.; Morera-Herreras, T.; De Deurwaerdere, P.; Miguelez, C. The noradrenergic system in Parkinson’s disease. Front. Pharmacol. 2020, 11, 435. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Politis, M.; Niccolini, F. Serotonin in Parkinson’s disease. Behav. Brain Res. 2015, 277, 136–145. [Google Scholar] [CrossRef] [PubMed]
  5. Ugrumov, M. Development of early diagnosis of Parkinson’s disease: Illusion or reality? CNS Neurosci. Ther. 2020, 26, 997–1009. [Google Scholar] [CrossRef] [PubMed]
  6. Ledeen, R.W.; Wu, G. The multi-tasked life of GM1 ganglioside, a true factotum of nature. Trends Biochem. Sci. 2015, 40, 407–418. [Google Scholar] [CrossRef] [PubMed]
  7. Chiricozzi, E.; Lunghi, G.; Di Biase, E.; Fazzari, M.; Sonnino, S.; Mauri, L. GM1 ganglioside is a key factor in maintaining the mammalian neuronal functions avoiding neurodegeneration. Int. J. Mol. Sci. 2020, 21, 868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Ledeen, R.W.; Chowdhury, S.; Lu, Z.H.; Chakraborty, M.; Wu, G. Systemic deficiency of GM1 ganglioside in Parkinson’s disease tissues and its relation to the disease etiology. Glycoconj. J. 2021, in press. [Google Scholar]
  9. Mullin, S.; Hughes, D.; Mehta, A.; Schapira, A.H.V. Neurological effects of glucocerebrosidase gene mutations. Eur. J. Neurol. 2019, 26, 388–393. [Google Scholar] [CrossRef]
  10. Alcalay, R.N.; Levy, O.A.; Waters, C.C.; Fahn, S.; Ford, B. Glucocerebrosidase activity in Parkinson’s disease with and without GBA mutation. Brain 2015, 135, 2648–2658. [Google Scholar]
  11. Yohe, H.C.; Wallace, P.K.; Berenson, C.S.; Ye, S.; Reinhold, B.B.; Reinhold, V.N. The major gangliosides of human peripheral blood monocytes/macrophages: Absence of ganglio series structures. Glycobiology 2001, 11, 831–841. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wu, G.; Lu, Z.-H.; Gabius, H.-J.; Ledeen, R.W.; Bleich, D. Ganglioside GM1 deficiency in effector T cells from NOD mice induces resistance to regulatory T-cell suppression. Diabetes 2011, 60, 2341–2349. [Google Scholar] [CrossRef] [Green Version]
  13. Ledeen, R.W.; Wu, G. Gangliosides of the Nervous System. In Methods in Molecular Biology; Sonnino, S., Prinetti, A., Eds.; Elsevier: New York, NY, USA, 2018; Volume 1804, pp. 19–25. [Google Scholar]
  14. Ledeen, R.W.; Wu, G. Gangliosides, α-Synuclein, and Parkinson’s disease. In Progress in Molecular Biology and Translational Science: Gangliosides in Health and Disease; Schnaar, R.L., Lopez, P.H.H., Eds.; Academic Press: London, UK, 2018; Volume 156, pp. 435–454. [Google Scholar]
  15. Hadaczek, P.; Wu, G.; Sharma, N.; Ciesielska, A.; Bankiewicz, K.; Davidow, A.L.; Lu, Z.-H.; Forsayeth, J.; Ledeen, R.W. GDNF signaling implemented by GM1 ganglioside; failure in Parkinson’s disease and GM1-deficient murine model. Exp. Neurol. 2015, 263, 177–189. [Google Scholar] [CrossRef] [PubMed]
  16. Pascual, A.; Hidalgo-Figueroa, M.; Piruat, J.L.; Pintado, C.O.; Gomez-Diaz, R.; Lopez-Barneo, J. Absolute requirement of GDNF for adult catecholaminergic neuron survival. Nat. Neurosci. 2008, 11, 755–761. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Mutoh, T.; Tokuda, A.; Miyadai, T.; Hamaguchi, M.; Fujiki, N. Ganglioside GM1 binds to the Trk protein and regulates receptor function. Proc. Natl. Acad. Sci. USA 1995, 92, 5087–5091. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Pitto, J.; Mutoh, T.; Kuriyama, M.; Ferraretto, A.; Palestini, P.; Masserini, M. Influence of endogenous GM1 ganglioside on TrkB activity in cultured neurons. FEBS Lett. 1998, 439, 93–96. [Google Scholar] [CrossRef] [Green Version]
  19. Martinez, Z.; Zhu, M.; Han, S.; Fink, A.L. GM1 specifically interacts with α-synuclein and inhibits fibrillation. Biochemistry 2007, 46, 1868–1877. [Google Scholar] [PubMed]
  20. Bartels, T.; Kim, N.C.; Luth, E.S.; Selkoe, D.J. N-Alpha-acetylation of α-synuclein increases its helical folding propensity, GM1 binding specificity and resistance to aggregation. PLoS ONE 2014, 9, 103727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Wu, G.; Lu, Z.-H.; Kulkarni, N.; Ledeen, R.W. Deficiency of ganglioside GM1 correlates with Parkinson’s disease in mice and humans. J. Neurosci. Res. 2012, 90, 1997–2008. [Google Scholar] [PubMed]
  22. Schneider, J.S.; Aras, R.; Williams, C.K.; Koprich, J.B.; Brotchie, J.M.; Singh, V. GM1 ganglioside modifies α-synuclein toxicity and is neuroprotective in a rat α-synuclein model of Parkinson’s disease. Sci. Rep. 2019, 9, 8362. [Google Scholar]
  23. Pellicano, C.; Benincasa, D.; Pisani, V.; Buttarelli, F.R.; Giovannelli, M.; Pontieri, F.E. Prodromal non-motor symptoms of Parkinson’s disease. Neuropsych. Dis. Treat. 2007, 3, 145–151. [Google Scholar] [CrossRef] [PubMed]
  24. Haliday, G.M.; Barker, R.A.; Rowe, D.B. Non-Dopamine Lesions in Parkinson’s Disease; Oxford University Press: Oxford, UK, 2011. [Google Scholar]
  25. Bhuiyan, R.H.; Ohmi, Y.; Ohkawa, Y.; Zhang, P.; Takano, M.; Hashimoto, N.; Okajima, T.; Furukawa, K.; Furukawa, K. Loss of enzyme activity in mutated B4GALNT1 gene products in patients with Hereditary Spastic Paraplegia results in relatively mild neurological disorders: Similarity with phenotypes of B4galnt1 knockout mice. Neuroscience 2019, 397, 94–106. [Google Scholar] [CrossRef]
  26. Nalls, M.A.; Blauwendraat, C.; Vallerga, C.L.; Heilbron, K.; Bandres-Ciga, S.; Chang, D.; Tan, M.; Kia, D.A.; Noyce, A.J.; Xue, A.; et al. Identification of novel risk loci, causal insights, and heritable risk for Parkinson’s disease: A meta-analysis of genome-wide association studies. Lancet Neurol. 2019, 18, 1091–1102. [Google Scholar] [CrossRef]
  27. Schneider, J. Altered expression of genes involved in ganglioside biosynthesis in substantia nigra neurons in Parkinson’s disease. PLoS ONE 2018, 13, e0199189. [Google Scholar] [CrossRef] [PubMed]
  28. Tsai, Y.T.; Yu, R.K. Epigenetic activation of mouse ganglioside synthase genes: Implications for neurogenesis. J. Neurochem. 2014, 128, 101–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Itokazu, Y.; Tsai, Y.T.; Yu, R.K. Epigenetic regulation of ganglioside expression in neural stem cells and neuronal cells. Glycoconj. J. 2016, 34, 749–756. [Google Scholar] [CrossRef] [PubMed]
  30. Svennerholm, L.; Bostrom, K.; Jungbjer, B.; Olsson, L. Membrane lipids of adult human brain: Lipid composition of frontal and temporal lobe in subjects of age 20 to 100 years. J. Neurochem. 1994, 63, 1802–1811. [Google Scholar] [CrossRef]
  31. Svennerholm, L.; Bostrom, K.; Fredman, P.; Mansson, J.E.; Rosengren, D.; Rynmark, B.M. Human brain gangliosides: Developmental changes from early fetal stage to advanced age. Biochim. Biophys. Acta 1989, 1005, 109–117. [Google Scholar] [CrossRef]
  32. Miyagi, T.; Yamaguchi, K. Mammalian sialidases: Physiological and pathological roles in cellular functions. Glycobiology 2012, 22, 880–896. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Schneider, J.S.; Sendek, S.; Daskalakis, C.; Cambi, F. GM1 ganglioside in Parkinson’s disease: Results of a five-year open study. J. Neuro. Sci. 2010, 292, 45–51. [Google Scholar] [CrossRef] [PubMed]
  34. Schneider, J.S.; Gollomp, S.M.; Sendek, S.; Colcher, A.; Cambi, F.; Du, W. A randomized, controlled, delayed start trial of GM1 ganglioside in treated Parkinson’s disease patients. J. Neurol. Sci. 2013, 324, 140–148. [Google Scholar] [CrossRef] [Green Version]
  35. Wu, G.; Lu, Z.-H.; Seo, J.H.; Alselehdar, S.K.; DeFrees, S.; Ledeen, R.W. Mice deficient in GM1 manifest both motor and non-motor symptoms of Parkinson’s disease; successful treatment with synthetic GM1 ganglioside. Exp. Neurol. 2020, 329, 113284. [Google Scholar] [CrossRef]
  36. Caronti, B.; Tanda, G.; Colosimo, C.; Ruggieri, S.; Calderaro, C.; Palladini, G.; Pontieri, F.E.; Di Chiara, G. Reduced dopamine in peripheral blood lymphocytes in Parkinson’s disease. Neuroreport 1999, 10, 2907–2910. [Google Scholar] [CrossRef] [PubMed]
  37. Caronti, B.; Antonini, G.; Caslderaro, C.; Ruggieri, S.; Palladini, G.; Pontieri, F.E.; Colosimo, C. Dopamine transporter immunoreactivity in peripheral blood lymphocytes in Parkinson’s disease. J. Neural Trans. 2001, 108, 803–807. [Google Scholar] [CrossRef] [PubMed]
  38. Huebecker, M.; Moloney, E.B.; van der Spoel, A.C.; Priestman, D.A.; Isacson, O.; Hallett, P.J.; Platt, F.M. Reduced sphingolipid hydrolase activities, substrate accumulation and ganglioside decline in Parkinson’s disease. Mol. Neurodegener. 2019, 14, 1–21. [Google Scholar] [CrossRef] [Green Version]
  39. Padmanabhan, S.; Lanz, T.A.; Gorman, D.; Wolfe, M.; Joyce, A.; Cabrera, C.; Lawrence-Henderson, R.; Levers, N.; Joshi, N.; Ma, T.C.; et al. An assessment of LRRK2 serine 935 phosphorylation in human peripheral blood mononuclear cells in idiopathic Parkinson’s disease and G2019S LRRK2 cohorts. J. Parkinson’s Dis. 2020, 10, 623–629. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Lees, M.B.; Paxman, S. Modification of the Lowry procedure for the analysis of proteolipid protein. Anal. Biochem. 1972, 47, 184–192. [Google Scholar] [CrossRef]
  41. Skowron, M.; Zakrzewski, R.; Ciesielski, W. Application of thin-layer chromatography image analysis technique in quantitative determination of sphingomyelin. J. Anal. Chem. 2016, 71, 808–813. [Google Scholar] [CrossRef]
  42. Youden, W.J. Index for Rating Diagnostic Tests. Cancer 1950, 3, 32–35. [Google Scholar] [CrossRef]
  43. Schisterman, E.F.; Perkins, N.J.; Liu, A.; Bondell, H. Optimal cut-point and its corresponding Youden index to discriminate individuals using pooled blood samples. Epidemiology 2005, 16, 73–81. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Detection and quantification of GM1 and GD1a in PBMCs. (a) Image of autoradiographic film depicting HPTLC separation and CtxB-HRP detection of GM1, GD1a and other ganglio-series following N’ase treatment. (b) GM1 and GD1a are shown as the mean and range of levels by upper and lower lines for PBMCs from seven healthy controls, nine sPD, and seven PD-GBA patients. One-way ANOVA revealed a significant difference in GM1 between controls and the other two groups: F(2, 20) = 19.45, p < 0.0001. Tukey’s HSD post-hoc testing revealed that sPD and PD-GBA patients both had significantly less GM1 than healthy controls (p = 0.0008 and p < 0.0001 resp.). One-way ANOVA revealed a significant difference in GD1a between controls and the other two groups: F(2, 20) = 33.81, p < 0.0001. Tukey’s HSD post-hoc testing revealed that sPD and PD-GBA patients both had significantly less GD1a than healthy controls (p < 0.0001 for both). Additionally, PD-GBA patients had lower amount of GM1 (p = 0.13) and GD1a (p = 0.74) compared to sPD patients, though these were not significant.
Figure 1. Detection and quantification of GM1 and GD1a in PBMCs. (a) Image of autoradiographic film depicting HPTLC separation and CtxB-HRP detection of GM1, GD1a and other ganglio-series following N’ase treatment. (b) GM1 and GD1a are shown as the mean and range of levels by upper and lower lines for PBMCs from seven healthy controls, nine sPD, and seven PD-GBA patients. One-way ANOVA revealed a significant difference in GM1 between controls and the other two groups: F(2, 20) = 19.45, p < 0.0001. Tukey’s HSD post-hoc testing revealed that sPD and PD-GBA patients both had significantly less GM1 than healthy controls (p = 0.0008 and p < 0.0001 resp.). One-way ANOVA revealed a significant difference in GD1a between controls and the other two groups: F(2, 20) = 33.81, p < 0.0001. Tukey’s HSD post-hoc testing revealed that sPD and PD-GBA patients both had significantly less GD1a than healthy controls (p < 0.0001 for both). Additionally, PD-GBA patients had lower amount of GM1 (p = 0.13) and GD1a (p = 0.74) compared to sPD patients, though these were not significant.
Ijms 22 11522 g001
Figure 2. Quantification of GM1 in PBMCs. (a) Image of autoradiographic film depicting HPTLC separation and CtxB-HRP detection of GM1 derived from PBMCs from healthy controls, sPD and PD-GBA patients. (b) GM1 levels are shown as the mean represented by a solid line and range by upper and lower lines for PBMCs from healthy controls (n = 14), sPD (n = 14), and PD-GBA patients (n = 6). The data were normalized to the amount of lipid fraction applied for 1.2 × 106 PBMCs. One-way ANOVA showed a significant difference in GM1 between groups; F(2,31) = 59.62, p < 0.0001. Tukey’s HSD post-hoc testing revealed that sPD and PD-GBA patients both had a significantly lower amount of GM1 compared to healthy controls (p < 0.0001 respectively). Additionally, PD-GBA patients had a significantly lower amount of GM1 compared to sPD (p = 0.0005).
Figure 2. Quantification of GM1 in PBMCs. (a) Image of autoradiographic film depicting HPTLC separation and CtxB-HRP detection of GM1 derived from PBMCs from healthy controls, sPD and PD-GBA patients. (b) GM1 levels are shown as the mean represented by a solid line and range by upper and lower lines for PBMCs from healthy controls (n = 14), sPD (n = 14), and PD-GBA patients (n = 6). The data were normalized to the amount of lipid fraction applied for 1.2 × 106 PBMCs. One-way ANOVA showed a significant difference in GM1 between groups; F(2,31) = 59.62, p < 0.0001. Tukey’s HSD post-hoc testing revealed that sPD and PD-GBA patients both had a significantly lower amount of GM1 compared to healthy controls (p < 0.0001 respectively). Additionally, PD-GBA patients had a significantly lower amount of GM1 compared to sPD (p = 0.0005).
Ijms 22 11522 g002
Figure 3. ROC curves of GM1 and GD1a in PBMCs measured using the original analytical method. The ROC curves were used to determine the probability of using GM1 expression in PBMCs to differentiate sPD (a) and PD-GBA (b) patients from controls; and PD-GBA patients from sPD patients (c). The same was analyzed using GD1a expression in PBMCs comparing healthy controls and sPD patients (d); healthy controls and PD-GBA patients; (e) sPD and PD-GBA patients (f). The area under the curve (AUC) and the sensitivity (%) and specificity (%) for cutoff levels are indicated on the plots.
Figure 3. ROC curves of GM1 and GD1a in PBMCs measured using the original analytical method. The ROC curves were used to determine the probability of using GM1 expression in PBMCs to differentiate sPD (a) and PD-GBA (b) patients from controls; and PD-GBA patients from sPD patients (c). The same was analyzed using GD1a expression in PBMCs comparing healthy controls and sPD patients (d); healthy controls and PD-GBA patients; (e) sPD and PD-GBA patients (f). The area under the curve (AUC) and the sensitivity (%) and specificity (%) for cutoff levels are indicated on the plots.
Ijms 22 11522 g003
Figure 4. ROC curves of GM1 in PBMCs measured using the modified analytical method. The ROC curves were used to determine the probability of differentiating sPD (a) and PD-GBA (b) patients from controls; and PD-GBA patients from sPD patients (c) using GM1 expression levels in PBMCs. The area under the curve (AUC) and the sensitivity (%) and specificity (%) for cutoff levels are indicated on the plots.
Figure 4. ROC curves of GM1 in PBMCs measured using the modified analytical method. The ROC curves were used to determine the probability of differentiating sPD (a) and PD-GBA (b) patients from controls; and PD-GBA patients from sPD patients (c) using GM1 expression levels in PBMCs. The area under the curve (AUC) and the sensitivity (%) and specificity (%) for cutoff levels are indicated on the plots.
Ijms 22 11522 g004
Table 1. Patient demographics and clinical assessment of sPD and PD-GBA patients in Cohort 1. Data are mean ± SEM. GM1 and GD1a levels in PBMCs were measured using our original analytical method. UPDRS I and II were calculated for seven healthy controls, nine sPD, and five PD-GBA. UPDRS III and total UPDRS were calculated for seven healthy controls, nine sPD, and four PD-GBA.
Table 1. Patient demographics and clinical assessment of sPD and PD-GBA patients in Cohort 1. Data are mean ± SEM. GM1 and GD1a levels in PBMCs were measured using our original analytical method. UPDRS I and II were calculated for seven healthy controls, nine sPD, and five PD-GBA. UPDRS III and total UPDRS were calculated for seven healthy controls, nine sPD, and four PD-GBA.
Cohort I
ParameterControlsPDPD-GBA
N797
Sex (male/female)2/56/32/5
Age62.39 ± 3.5866.28 ± 3.1565.96 ± 2.94
Age at disease onsetN/A57.00 ± 2.8962.57 ± 2.81
Disease durationN/A9.28 ± 2.113.39 ± 0.78
UPDRS I and II0.14 ± 0.148.56 ± 1.437.60 ± 3.5
UPDRS III0.29 ± 0.1816.89 ± 2.5415.75 ± 3.94
UPDRS Total0.43 ± 0.2025.44 ± 3.1120.25 ± 5.36
Table 2. Patient demographics and clinical assessment of sPD and PD-GBA patients in Cohort II. Data are mean ± SEM. GM1 in PBMCs were measured using our modified analytical method. UPDRS I and II were calculated for 12 healthy controls, 11 sPD, and four PD-GBA. UPDRS III and total UPDRS were calculated for 11 healthy controls, 11 sPD, and five PD-GBA.
Table 2. Patient demographics and clinical assessment of sPD and PD-GBA patients in Cohort II. Data are mean ± SEM. GM1 in PBMCs were measured using our modified analytical method. UPDRS I and II were calculated for 12 healthy controls, 11 sPD, and four PD-GBA. UPDRS III and total UPDRS were calculated for 11 healthy controls, 11 sPD, and five PD-GBA.
Cohort II
ParameterControlsPDPD-GBA
N14146
Sex (male/female)6/89/53/3
Age68.86 ± 2.2564.38 ± 2.8465.12 ± 2.79
Age at disease onsetN/A57.07 ± 3.0857.00 ± 4.28
Disease durationN/A7.76 ± 1.668.12 ± 1.93
UPDRS I and II Score0.50 ± 0.267.64 ± 1.448.00 ± 1.29
UPDRS III Score1.83 ± 0.7219.45 ± 2.4118.60 ± 5.47
UPDRS Total2.55 ± 0.7727.09 ± 3.3122.00 ± 4.80
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alselehdar, S.K.; Chakraborty, M.; Chowdhury, S.; Alcalay, R.N.; Surface, M.; Ledeen, R. Subnormal GM1 in PBMCs: Promise for Early Diagnosis of Parkinson’s Disease? Int. J. Mol. Sci. 2021, 22, 11522. https://doi.org/10.3390/ijms222111522

AMA Style

Alselehdar SK, Chakraborty M, Chowdhury S, Alcalay RN, Surface M, Ledeen R. Subnormal GM1 in PBMCs: Promise for Early Diagnosis of Parkinson’s Disease? International Journal of Molecular Sciences. 2021; 22(21):11522. https://doi.org/10.3390/ijms222111522

Chicago/Turabian Style

Alselehdar, Samar K., Monami Chakraborty, Suman Chowdhury, Roy N. Alcalay, Matthew Surface, and Robert Ledeen. 2021. "Subnormal GM1 in PBMCs: Promise for Early Diagnosis of Parkinson’s Disease?" International Journal of Molecular Sciences 22, no. 21: 11522. https://doi.org/10.3390/ijms222111522

APA Style

Alselehdar, S. K., Chakraborty, M., Chowdhury, S., Alcalay, R. N., Surface, M., & Ledeen, R. (2021). Subnormal GM1 in PBMCs: Promise for Early Diagnosis of Parkinson’s Disease? International Journal of Molecular Sciences, 22(21), 11522. https://doi.org/10.3390/ijms222111522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop