Patent Application
20230393152
The invention concerns the identification of novel combinations of metabolites specific for neurodegenerative diseases, in particular Parkinson's disease (PD), for a use as biomarkers in the early diagnosis of neurodegenerative diseases, in particular Parkinson's disease at prodromal or de novo stage.
In the description below, references in square brackets ([ ]) refer to the list of references at the end of the text.
Parkinson's disease (PD) is one of the most common neurodegenerative diseases worldwide, affecting approximately 1% of adults over the age of 60 [1]. As soon as it appears, the pathology evolves gradually and continuously according to three stages, i.e. preclinical, prodromal, and clinical stage, corresponding respectively to the appearance of neurodegeneration without symptoms (preclinical phase), then the progression of the neurodegenerative process leads to the appearance of non-motor symptoms (prodromal phase), and finally to the appearance of motor symptoms (clinical phase) [2]. It is these symptoms, present only in the last stage, which are sufficiently characteristic of PD to allow its diagnosis [3]. However, they appear when 70-80% of the dopaminergic neurons of the nigrostriatal pathway have already degenerated [1], therefore generally too late to consider neuroprotective intervention [4]. Furthermore, post-mortem histopathological analyzes of brains show that currently the diagnostic accuracy of PD in patients with motor symptoms for less than five years [5] is only 53%. Thus, too late diagnosis and its inaccuracy [3, 6, 7] strongly contribute to limited and only symptomatic patient management and therapies which, moreover, become ineffective after several years.
Moreover, although many molecules have already been tested and have shown promising theoretical or preclinical evidence of efficacy, there is still no agent or neuroprotective therapy allowing slowing down, stopping or effectively reversing neurodegeneration in patients with PD, in part because these therapeutic strategies are applied too late [8-10].
Thus, being able to detect PD with a high degree of specificity, during the early phases, especially in patients in the prodromal phase or newly and recently diagnosed (de novo early clinical phase), represents a crucial challenge to enable the success of curative or preventive strategies of PD, and the development of new therapies capable of delaying or slowing down its progression [11, 12].
The «prodromal» phase is characterized by a set of non-motor symptoms including neuropsychiatric disorders (e.g. apathy, depression, anxiety), Although they are expressed in 40 to 50% of Parkinson's patients, they are not exclusive to the disease and can be encountered in other pathologies. Patients newly and recently diagnosed with Parkinson's disease correspond to an early clinical phase called «de novo» where neuroprotective intervention can still be considered to slow the progression of the disease.
Thus, finding reliable and specific biomarkers of the prodromal or de novo phase, suitable and highly predictive for a clinical use, is an important challenge, which motivated the study of the inventors. Their starting hypothesis is that the metabolome—i.e. the global pool of metabolites—which results from the expression of genome but also reflects the impact of many environmental factors [15], could contain potential biomarkers of the disease [16-19], now considered as a multifactorial disease resulting from complex interactions between genetic and environmental factors [14]. Moreover, metabolic alterations have already been described in patients expressing neuropsychiatric symptoms similar to those observed in the prodromal phase of PD [20]. In addition, it has been shown that the analysis of metabolic deregulations observed in the serum of patients newly diagnosed with PD made it possible to precisely distinguish them from healthy subjects [21]. Finally, alterations in the plasma metabolome have been correlated with disease progression in Parkinsonian patients [22]. However, the study of the metabolome as a tool for predicting the prodromal phase of PD has never been carried out. To access this prodromal phase, the use of animal models is essential. Nevertheless, although many animal models of PD have been developed, none recapitulates all the phenotypic and etiological characteristics of the disease [23].
Thus, although the early detection of PD represents a key step for the development of curative therapies or capable of delaying or slowing down its progression, there is currently no specific tool to allow the diagnosis of patients with PD before the appearance of motor symptoms, in particular during the prodromal phase. The discovery of easily measurable and highly predictive biomarkers of the prodromal phase of the disease appears as a milestone for accelerating the development of curative therapies and improving patient care.
Also, the inventors have developed an original strategy making it possible to propose new biomarkers capable of specifically and early diagnosing a neurodegenerative disease, for example Parkinson's disease (PD), in particular Parkinson's disease at the prodromal or de nova stage. This strategy, based on the use of different animal models and different patient cohorts, overcame the intrinsic limitations and increased the incomplete predictive value of each model taken individually, revealing a characteristic set of deregulated metabolites usable as a potential biomarker. Finally, the study was extended to human serum samples from two biobanks in order to validate the clinical relevance of this new PD biomarker.
Thus, via behavioral and histological assessment and using metabolomic approaches based on nuclear magnetic resonance (NMR) spectroscopy followed by advanced multivariate statistics, the inventors studied the metabolic changes (i) in serum and brain tissue of three different animal models of PD, two rodent models and a non-human primate model, mimicking the different phases of the disease [24-26], including the prodromal phase, and (ii) in the serum of two cohorts of patients with de novo PD (n=95) (United States, Italy). Deregulations common to serum and tissue samples have been observed in animals, showing that blood might reflect central (brain) dysfunction [27], at least partially.
The 6-hydroxydopamine (6-OHDA) rat model was first used because it stably mimics each stage of PD, with no evolution between stages, and it has already been validated for the study of prodromal phases [24]. Moreover, it represents a phenotypically correct model of PD, with good isomorphism and good predictive value [28], although it is etiologically limited. The α-synuclein rat model was then used because it has the advantage of being progressive in the establishment of the phenotype, which makes longitudinal studies possible. In addition, it is closer to human pathology, especially with regard to the anatomo-pathological characteristics and the physiopathological mechanisms of the disease [30].
Finally, non-human MPTP (1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine) primates were used because they share greater genome sequence identity with humans and important neuroanatomical similarities.
Metabolic changes were analyzed at the blood (serum) level in order to find metabolites that are easy to be transposed into the clinic, and in specific regions of the brain differently affected by the neurodegenerative process (the dorsal striatum (DS) and the nucleus accumbens (NAcc)) in order to potentially associate changes in these metabolites to the pathophysiological processes involved in PD.
From a translational perspective, the metabolomic study was extended to serums from newly diagnosed parkinsonian patients (de novo phase) from two different biobanks (United States (NIH)+/−2 patients, Italy). The translational strategy revealed common metabolic deregulations in the serum of different animal models and patients with PD. Some of them have been found in tissue samples (animal models), possibly reflecting pathophysiological mechanisms associated with development of the PD.
Thus, this inter-species study has made it possible to reinforce the relevance of circulating metabolites for the diagnosis of PD. In addition, it proposes a composite biomarker constructed with a combination of 6 metabolites to distinguish animals mimicking PD from control animals (sham), regardless of the stage of the disease, and in particular from the prodromal or de novo phase. It is important to note that the combination of the same metabolites made it possible to distinguish patients with PD from healthy subjects.
It is interesting to note that in animals, certain metabolic dysregulations appeared before the motor symptoms and could thus represent early biomarkers of the disease.
Thus, the inventors have demonstrated that the relative concentrations of 6 metabolites: acetoacetate, betaine, β-hydroxybutyrate (BHB), creatine, pyruvate, and valine, combined together according to a logistic regression method, could constitute a precise composite biomarker of PD, in animals and in humans.
The use of this biomarker made it possible to distinguish American parkinsonian patients from healthy controls (comparison group) with a high sensitivity (0.95), and Italian parkinsonian patients with a slightly lower sensitivity (0.83). The divergent pre-analytical procedure applied for serum collection in the two cohorts (that is to say used tubes, sample processing, etc.) could be a major source of experimental variability that could explain this difference [31, 32]. Furthermore, it should be noted that the human samples were collected according to a very strict and standardized experimental protocol, but which was not designed for metabolomics. This reflects the power of metabolomics and the robustness of the results, and supports the view that serum biomarkers, easier to use and less invasive than imaging methods, could show similar efficacy [33].
Interestingly, of the six metabolites used to construct the composite PD biomarker in animal or human models, betaine, BHB, and pyruvate dysregulation were corrected by chronic treatment with an agonist dopaminergic (i.e. Pramipexole or Pra), which corrects also motivational deficits in 6-OHDA rats mimicking the prodromal phase. This pharmacological validation highlights their potential involvement in cellular dysregulations involved in neuropsychiatric symptoms associated with the prodromal or de novo stage. The metabolic deregulations reported in the present study could therefore also provide some clues about the mechanisms involved in the progression of PD.
In summary, the present study reported, for the first time, common serum metabolic alterations in three animal models mimicking PD, although from different cohorts and species. Consequently, this common signature is likely to reflect specific alterations related to the pathophysiology of PD. In particular, specific deregulations characteristic of the prodromal or de novo stage have been highlighted. In addition to providing clues for further study of metabolic alterations in animal models, metabolic dysregulations have also been found in the serums of newly diagnosed Parkinsonian patients (that is to say in the de novo clinical phase) from two biobanks, not specifically designed for metabolomics. Logistic regression, using the levels of the same metabolites used for animals, allowed a very adequate discrimination between healthy subjects and patients with PD in two cohorts of different origins (at least 0.83).
Overall, the present cross-species study demonstrated the utility and reproducibility of untargeted metabolomics as a non-invasive approach for PD biomarker research. This approach seems particularly promising for patients with PD, as the serum is easily accessible, and the biomarker may be relevant even in the early stages of disease when other methods fail.
The subject of the present invention is therefore an in vitro method for the early diagnosis of a neurodegenerative disease from a biological sample of a patient comprising:
According to a particular mode of the method of the present invention, step a) further comprises the in vitro measurement of the level of expression of at least one metabolite selected from the group consisting of alanine, lactate, DMSO2, glycine, serine, threonine, myo-inositol, and leucine.
According to a particular mode of the method of the present invention, step a) comprises or consists of measuring the levels of expression of acetoacetate, betaine, β-hydroxybutyrate (BHB), creatine, pyruvate, and the valine.
The term «early diagnosis», within the meaning of the present invention, means a diagnosis made at a stage of the neurodegenerative disease where the discriminating symptoms have not yet appeared or in newly and recently diagnosed patients. For example, a diagnosis would be considered early if performed at the prodromal or de novo stage of PD.
The term «neurodegenerative disease» within the meaning of the present invention means a progressive pathology which affects the brain or more generally the central nervous system, resulting in the death of nerve cells. The most famous and most frequent are Alzheimer's disease and Parkinson's disease, but there are others (multiple sclerosis, Huntington's disease, Creutzfeldt-Jakob disease, amyotrophic lateral sclerosis (also called Charcot's disease), . . . ). In particular, it is about Parkinson's disease.
The term «biological sample» within the meaning of the present invention means a biological sample (e.g. of bodily fluid (serum, blood, urine, saliva, etc.), or tissue (brain, marrow, liver, skin, etc.) taken after a 12-hour fast in a patient.
The term «biological reference sample» within the meaning of the present invention means a biological sample (e.g. of bodily fluid (serum, blood, urine, saliva), or tissue (brain, marrow, liver, skin, etc.)) taken after a 12-hour fast in a healthy individual (control/comparison group/sham), that is to say presenting no neurological disease, in particular no symptoms of PD or who may be a forerunner of PD, no family history of movement disorder and/or PD, and no specific medical conditions.
The term «reference value» within the meaning of the present invention, the value of the level of expression of a metabolite in a healthy individual, that is to say not presenting any neurological disease, in particular no symptoms of PM or may be a precursor to PD, no family history of movement disorder and/or PD, and no specific medical conditions.
The term «difference in expression levels» means, within the meaning of the present invention, the variation in the expression levels of the metabolites measured in a patient compared to the expression levels of the same metabolites in a healthy individual (reference values), that is to say having no neurological disease; in particular no symptoms of PD or which may be a precursor to PD, no family history of movement disorder and/or PD, and no specific medical condition. For example, in the case of Parkinson's disease.
The subject of the present invention is an in vitro method for the early diagnosis of Parkinson's disease from a biological sample from a patient comprising:
−45.60+0.28EBHE−4.54EPyruvate−11.23Evaline+7.05EAcetoacetate+2.85Ecreatine−7.79EBetaine; or
−42.78+0.16EBHB−3.31EPyruvate−12.01Evaline+7.0G6EAcetoacetate+3.14ECreatine−8.27EBetain.
The term «calculated value of Logit(P)» within the meaning of the present invention means a standard value obtained from the application of a logistic regression using the measurement of the levels of expression of the metabolites mentioned above. In the case of Parkinson's disease, according to the two logistic regressions above using the following 6 metabolites: acetoacetate, betaine, p-hydroxybutyrate (BHB), creatine, pyruvate, and valine, any individual obtaining a value greater than 0.316 (arbitrary unit, a.u.) will be considered to have PD and any individual scoring less than 0.279 (arbitrary unit, a.u.) will be considered not to have PD.
According to a particular mode of the method of the present invention, the biological sample of the patient comes from a patient at the prodromal or de novo stage of PD.
The term «prodromal stage of PD» means, within the meaning of the present invention, a stage where the diagnosis of PD could not be established definitively and/or with certainty, on the basis of the presented symptoms. For example, this concerns any patient presenting with non-motor symptoms frequently observed in patients who eventually develop the disease, and which could justly represent prodromal symptoms (neuropsychiatric symptoms (apathy, depression, anxiety), loss of sense of smell, constipation, sleep disorders, digestive disorders, cognitive deficits) and/or any patient who presents in his genealogical line people already suffering from the disease.
The term «de nova stage of PD» or «de novo patients» within the meaning of the present invention means newly and recently diagnosed patients who are at an early clinical stage where a neuroprotective intervention can still be envisaged to slow the progression of the disease.
According to a particular embodiment of the method of the present invention, the biological sample is for example a serum sample.
According to a particular embodiment of the present invention, the step of measuring expression levels is carried out, for example, by NMR spectroscopy measurements.
The present invention also relates to a biomarker intended for the early diagnosis of a neurodegenerative disease comprising the combination of the following at least six metabolites: acetoacetate, betaine, β-hydroxybutyrate (BHB), creatine, pyruvate, and valine.
According to a particular mode of the present invention, the biomarker further comprises the combination of at least one metabolite selected from the group consisting of alanine, lactate, DMSO2, glycine, serine, threonine, myo-inositol, and leucine.
According to a particular embodiment of the present invention, the biomarker comprises or consists of the combination of the following six metabolites: acetoacetate, betaine, β-hydroxybutyrate (BHB), creatine, pyruvate, and valine.
A subject of the present invention is also the in vitro use of a combination of metabolites according to the present invention, as a biomarker of a neurodegenerative disease, in particular of Parkinson's disease.
Mean±SEM, one-way ANOVA, followed by Tukey post-hoc test and correction for multiple comparisons. *: p≤0.05; **: p≤0.01; ***: p≤0.001.
a.u.: arbitrary units.
Scheme
The second cohort of rats (
Monkeys (n=8) received intramuscular injections of MPTP (0.3 to 0.5 mg/kg) every 4-5 days for 3 weeks [77]. The repeated administration of low doses of MPTP has been used to mimic a moderate stage, with slow disease development, and the onset of moderate dopaminergic injury. Serum samples were collected at the start of the experiment and after the monkeys reached a stable parkinsonian state (
All serum samples were analyzed by 1H NMR at 950 MHz and brain samples were subjected to HRMAS 1H NMR at 500 MHz.
Animals
Rats: the experiments were carried out on adult male Sprague-Dawley rats (Janvier, Le Genest-Saint-Isle, France), weighing approximately 300 g (7 weeks) at the start of the experiment, They were housed under standard laboratory and ethical conditions with an inverted light-dark cycle (12 p.m./light-dark cycle, with lights on at 7 p.m.) and with food and water available ad libitum.
Monkeys: the experiments were carried out on adult male Macaca fascicularis, weighing between 5 and 8 kg, aged 3 to 5 years and housed under standard conditions (light-dark cycle of 12:12 hours; 23° C.; 50% humidity).
The used protocols complied with the European Union Animal Protection Act 2010 and French Directive 2010/63, and were approved by the French National Ethics Committee (2013/113) No. 004 and by the local ethics committee C2EA84 and CELYNE C2EA.
Surgery
Bilateral injection of 6-OHDA: As described previously [18, 23], the rats received a subcutaneous injection of desipramine (15 mg/kg) 30 minutes before the injection of 6-OHDA in order to protect the noradrenergic neurons. All animals were then anesthetized by intraperitoneal injection of ketamine-xylazine (100-7 mg/kg), and placed in a stereotaxic frame (Kopf Instruments). The stereotaxic coordinates of the injection sites were as follows, according to the stereotaxic atlas of Paxinos and Watson [78] and in relation to the bregma: incisor bar placed at −3.2 mm: anteroposterior (AP)=−5.4 mm/lateral (L)=+/−1.8 mm/dorsoventral (V)=−8.1 mm. The animals received a bilateral injection of 2.3 μl of 6-OHDA (3 mg/ml) or physiological saline (NaCl 0.9%; control/sham), at a flow rate of 0.5 μl/min. For animals injected with 6-OHDA and presenting with transient starvation 2-3 days after the operation, supplementation with a high-calorie liquid diet and appetizing foods was started for 1-2 weeks, and stopped 10 days before serum collection. Animals that did not recover were excluded.
Bilateral injection of α-synuclein: As previously described [19], rats were anesthetized with 2.5% isoflurane and placed in a stereotactic frame (Kopf Instruments) in order to receive two bilateral injections of AAV-hA53Tα-syn (1 μl-7.0×1012 vg/ml) or AAV-GFP (7.0×1012 vg/ml, control/sham animals) in the SNc (AP=−5.1 and −5.6 mm/L=+1-2.2 mm/V=8 mm from the bregma). After being awakened, the animals were returned to the animal facility.
Behavioral Assessment
6-OHDA rats: The rats were subjected to 2 behavioral tests, 3 weeks after the operation, while the 6-OHDA lesion was stabilized [79].
Operant self-administration of sucrose (motivational component): Rats were trained to self-administer a 2.5% sucrose solution in operant chambers (Med Associates, St Albans, VT, USA) containing an active lever, reinforced, for which pressing causes the delivery of 0.2 ml of a sucrose solution associated with a low-intensity light stimulus, and an inactive, unreinforced lever, for which pressing causes neither delivery of sucrose nor light stimulus. This task was performed with a fixed ratio of 1, with a single press of the reinforced lever resulting in the delivery of a reward. Each session ended when 100 rewards were obtained, or at the end of the allotted time (1 hour). The number of obtained rewards was counted for each session by the MED-PC IV software.
Pitch readjustment test (motor component): The animals, held by the posterior third of their body, were moved over a length of 90 cm by a straight and regular movement from left to right, and vice versa, along a smooth surface table. The number of forelimb adjustments made during movement was counted [80]. The test was carried out three times, by two different experimenters, blind (ignorance of the nature of the experimental groups).
MPTP monkeys: The severity of the parkinsonian state was assessed using a rating scale taking into account classic motor symptoms (bradykinesia, rigidity, tremors, freezing, posture and position of the arms), spontaneous activities (arm movements, spontaneous eye movements and activity in the housing cage) and other activities (vocalization, triggered eye movements and feeding) [81].
Histology: Immunolabeling of TH and Quantification of Denervation in 6-OHDA Rats.
Brain harvest: The rats were sacrificed by decapitation 1 day after the last session of sucrose self-administration and 1 hour after the last injection of Pra or Veh (that is to say 6 weeks after the injection of 6-OHDA or saline solution). The brains were immediately frozen in liquid nitrogen and stored at −80° C. They were then sliced using a cryostat thermostated at −20° C. Sections of the regions of interest were taken to carry out both the immunostaining of TH and to carry out the metabolomic study. For 1HHRMAS NMR, thick sections of DS and Nacc tissue were pooled in inserts, and kept at −80° C. until NMR analysis.
For TH immunostaining, 14 μm coronal sections of the striatal levels of interest were taken and stored at −20° C. Immunostaining was then performed as previously described [18,23]. Briefly, after fixation with 4% paraformaldehyde, striatal sections were incubated with an anti-TH antibody (mouse monoclonal MAB5280, Millipore, France, 1:2500) overnight at 4° C. Then, the sections were incubated with a biotinylated goat anti-mouse IgG antibody (BA-9200, Vector Laboratories, Burlingame, CA, USA; 1:500) and the immunoreactivity revealed by the avidin-peroxidase conjugate (Vectastain ABC Elite, Vector Laboratories Burlingame, CA, USA).
Quantification of the extent of striatal dopaminergic denervation was determined with the computerized image analysis system ICS FrameWork (Calopix, version 2.9.2, TRIBVN, Châtillon, France) coupled to an optical microscope (Nikon, Eclipse 80i). After drawing the masks of the three striatal levels of interest, the optical densities (OD) were measured for each striatal subregion (DS and Nacc). The ODs were expressed as percentages relative to the mean optical density obtained from the homologous regions of the operated control animals (intracerebral injection of physiological saline).
Progression of Parkinson's Disease (PD) Score in the 6-OHDA Model
In order to assess the progression of PD (PDP) in 6-OHDA rats, a scale called «PDP score», and combining 3 criteria, was developed: 1) criterion 1: assessment of neuropsychiatric symptoms (ie. behavior apathetic type), 2) criterion 2: assessment of motor symptoms (i.e. fine motor deficit), these two criteria also being used clinically in the scales allowing the diagnosis of PD, and 3) criterion 3: the level of dopaminergic denervation in the DS, available only after postmortem histological analysis in the patient. Each criterion was quantified for each animal, with values ranging from 0 to 4, as shown in
For the assessment of behavioral deficits in the self-administration and pitch readjustment tests (criteria 1 and 2, respectively), the performance of each animal was compared before and after the injection of 6-OHDA (or physiological saline (sham). A value of 0 for each of these two criteria corresponded to an absence of deficit, or to a deficit of less than 30%, i.e. the normal daily variability of the used tests. The following subcategories have been set as follows: 1=low deficit (30%-50%); 2=mean symptoms (50%-70%); 3=strong symptoms (70%-90%); 4=significant or total deficit (>90%).
Considering that motor symptoms appear after 70% loss of striatal dopaminergic neurons in patients with PD, the value of dopaminergic denervation in DS (criterion 3) was used to define the boundary between the prodromal phase and the clinical phase.
With regard to the PDP score obtained by adding the values of the 3 criteria, the animal mimicking PD was classified as: 1) asymptomatic rat with weak lesions but no behavioral symptoms (PDP score=1 to 2), 2) prodromal type animal, i.e. presenting only neuropsychiatric disorders and a small DS lesion (PDP score=3 to 7) and 3) clinical type animal, i.e. presenting neuropsychiatric disorders; motor symptoms and an extensive DS lesion (PDP score=8 at 12). Control (sham) animals were assigned a score of 0.
Given the low number of asymptomatic animals and the fact that they are of no clinical relevance, they were excluded from the study.
Human Cohorts
Blood samples from PD patients and corresponding comparison groups subjects were obtained through 1) the Parkinson's Disease Biomarker Program (PDBP) Consortium, supported by the National Institute of Neurological Disorders and Vascular Accidents National Institutes of Health (NIH-USA), and 2) Santa Lucia Foundation (Rome, Italy). All subjects (patients and comparison groups) gave written informed consent. Inclusion criteria for patients from the United States included recent diagnosis (<1 year) of PD according to the new Movement Disorders Society (MDS) criteria. For the Italian cohort, the diagnostic criteria were less restrictive and patients diagnosed for longer than those of the US cohort were included (0 to 3 years after diagnosis), In both cohorts, the patients had not received any antiparkinsonian treatment at the time of inclusion, did not have dementia or psychiatric disorders, or any other medical condition that could compromise the study. Inclusion criteria for comparison groups were: absence of neurological disease, absence of family history of movement disorders, and absence of specific medical conditions. The US cohort included 21 PD patients without antiparkinsonian therapy, 9 PD patients who received Pra treatment, and 30 age- and gender-matched comparison groups. The Italian cohort included 21 patients with PD and 23 age- and gender-matched comparison groups. All clinical and demographic information collected from the two human cohorts has been summarized in Table 1 below (data are represented by mean value±SEM, M: Male; F: Female; *2 patients of unknown sex),
NMR Experiments
HRMAS 1H NMR of 6-OHDA rat brain samples: Just before HRMAS analysis, 10 μl of D2O was added to the inserts containing the brain biopsies, which were then sealed and placed in a 4 mm zirconium rotor.
Data acquisition: All HRMAS NMR spectra were acquired, as previously described [82], on a Bruker Advance Ill spectrometer (Billerica, USA) (IRMaGe, CEA Grenoble, France) at 500 MHz. The samples were rotated at 4000 Hz at the magic angle (54.7°) and the temperature was maintained at 4° C. for all experiments. One-dimensional spectra were acquired using a CPMG pulse sequence (TE=30 ms, 256 repetitions, 17 min). The residual water signal was pre-saturated for 1.7 s of relaxation.
Data processing: Quantification was performed with JMRUI software using quantification based on a quantum estimation procedure (QUEST) [83]. This procedure requires the use of a database of metabolites and a complete assignment of spectra. Nineteen metabolites were attributed and quantified: acetate, alanine, ascorbate, choline, creatine, gamma-aminobutyrate (GABA), glutamate, glutamine, glycine, glycerophosphocholine, glutathione, lactate, myo-inositol, N-acetylaspartate, phosphocholine, phosphocreatine, phosphoethanolamine, scyllo-inositol and taurine. The amplitude of the metabolite calculated by QUEST was normalized with respect to the total signal of the spectrum. CRLBs were calculated for each metabolite as an 2.5 estimate of the standard deviation of the fit.
1H NMR of Serum
Serum collection: For the rats, the blood was collected from the tail vein after 2 hours of fasting and under gaseous anesthesia with isoflurane (2%) and was stored in ice before being quickly centrifuged at 1600 g for 15 min at 4° C. The supernatant serum was collected and stored at −80° C. until the day of the NMR. The time before freezing never exceeded 30 min [84].
For monkeys, blood was collected from the saphenous vein after fasting for 12 hours, under anesthesia (atropine 0.05 mg/kg i.m. followed by Zoletil 15 mg/kg i.m.) and subjected to the same protocol as rat blood.
For humans, blood was collected after overnight fasting, stored at room temperature for 15-60 min before being centrifuged at 1200-1500 g at 4° C. for 10-15 min. The supernatant serum was then collected and stored at −80° C.
On the day of the NMR experiments, the samples were slowly thawed on ice and then quickly centrifuged to remove any cryo-precipitates. NMR tubes were filled with 60 μl of serum sample and 120 μl of 0.1 M phosphate buffered saline (PBS) in D2O (50% D2O, pH=7.4), and stored at 4° C. until acquisition by NMR.
Data acquisition: All animal and human serum samples were subjected to the same NMR protocol. 1H NMR experiments were performed on a Bruker Advance III NMR spectrometer at 950 MHz (IBS, Grenoble, France) by using a 3 mm cryoprobe. Spectra were recorded using a CPMG pulse sequence. The residual water signal was pre-saturated for 2 seconds of relaxation.
Peak assignment was performed using the 2-dimensional homonuclear 1H-1H (TOCSY) and heteronuclear 1H-13C (HSQC) experiments, and databases. If necessary, addition of selected metabolites was used to clarify ambiguous identifications.
Data processing: The time signals were transformed by the Fourier method and phased manually with Topspin software from Bruker, version 3.6.2. Then, other pre-processing steps (baseline correction, alignment, segmentation) were performed using the online NMRProcFlow v1.2 software (http://nmrprocflow.org). Segments of ppm were made over the 0-10 ppm spectral range, excluding the residual water peak, macromolecule/lipid signals and other regions containing pollutions (at δ 5.0-4.7; δ 5.25-5.38; δ 2.2-2.0; δ 1.4-1.25; δ 1.15-1; δ 0.9-0.5). Each segment was normalized to the total sum of the segments.
Statistical Analysis
Multivariate analysis: Data from liquid NMR or HRMAS were imported into SIMCA v.14 (Malmo, Sweden) for multivariate statistics. Unsupervised principal component analysis (PCA) was first used for the overall visualization of the distribution of all samples, followed by OPLS to find the discriminating metabolites associated with a specific disease stage. For the latter, either a continuous variable (that is to say the PDP score), either a discrete variable (that is to say group membership) was used to label each sample, leading to OPLS-DA in the latter case. The total number of components was determined using the cross-validation procedure, which produces the R2Y and Q2 factors that indicate the goodness of fit and the predictability of the model, respectively. A model is considered robust and predictive when both factors are ≥0.5. The scores were plotted in 2D against the first two components of the OPLS models, while the weights were plotted in 1D to mimic an NMR spectrum but with positive/negative peaks indicating up-regulated/down-regulated metabolites, respectively. Indeed, in this «statistical spectrum», the correlation of each NMR variable with the continuous ranking variable (score), ranging from 0 to 1, was color coded, ranging from dark gray (high correlation) to light gray (low correlation). Metabolites with a correlation ≥0.5, considered the most discriminating, were subjected to univariate analysis. Analysis of variance of cross-validated predictive residuals (CV-ANOVA) was used to assess the significance of the model.
Metabolites significantly varying in at least three of the four studied models (6-OHDA, α-synuclein, MPTP and human) were subjected to multiple logistic regression to assess their relevance as a potential predictive biomarker of PD at using the R software (version 3.6.1—R core team). The regressions, including cross-validation, were performed for prodromal type animals (6-OHDA and α-synuclein prodromal type rats). Clinical-type animals were used for external validation (clinical-type 6-OHDA and α-synuclein rats, and all primates). The same analyzes were performed for the US patient cohort. The Italian patient cohort was used for external validation. All combinations of retained metabolites were tested, and ROC curves were then generated for each of them to assess the goodness of fit using the area under the curve (AUC) and the optimal threshold that maximizes sensitivity and specificity.
Univariate Analysis:
All univariate analyzes were performed using Graphpad Prism 8 software (San Diego—USA). All results were expressed as mean values±SEM with a significance level set at 0.05.
6-OHDA models: For the self-administration tests and the step readjustment tests, a univariate analysis of the effect of Pra was performed by applying a one-way repeated measures (MR) ANOVA followed by a Sidak post-hoc test.
For histological and metabolomic data, a one-way ANOVA followed by Tukey post hoc test with correction for multiple comparisons was performed.
α-synuclein models: For the longitudinal study performed in α-synuclein animals, some values were missing, due to artifacts during NMR measurement, Data were therefore analyzed by fitting a proposed mixed model using a compound symmetry covariance matrix, and were fitted using the restricted maximum likelihood (REML) method. This model was followed by a Tukey post hoc test.
MPTP monkeys: Given the small number of animals, the non-parametric Mann-Whitney test was used.
Humans; Data were log-transformed given the non-Gaussian distribution of metabolite levels. First, t-tests were performed individually in the two different cohorts. Then a two-way ANOVA was applied with i) origin of cohort and ii) group (comparison group or PD) as factors, followed by Sidak post hoc test.
The 6-OHDA Rat Model Makes it Possible to Study the Different Stages of Parkinson's Disease
Although, in humans, the diagnosis of PD relies on a scale that incorporates assessments of neuropsychiatric and motor symptoms [22], its definitive confirmation is based on postmortem histological assessments, primarily of neuronal loss in the nigrostriatal pathway leading to DS denervation. To characterize 6-OHDA rats for clinical disease progression, a score, called the Parkinson's disease progression score (PDP), based on the same type of criteria as those used in humans, ie. the neuropsychiatric component—assessed by self-administration performance (motivation)—, fine motor skills—assessed by the performance of the step readjustment test—, and the extent of the DS lesion—assessed by the post-mortem histological analysis was assigned to them. Based on this score, animals with PD were divided into two categories, namely prodromal type animals and clinical type animals.
Tyrosine hydroxylase immunoreactivity (TH-IR) revealed that bilateral injection of 6-OHDA into the SNc led to partial nigrostriatal dopaminergic injury resulting in dopaminergic denervation in the DS and to a lesser extent in the Nacc (
With respect to fine motor skills, all animals showed a similar number of paw readjustments before 6-OHDA/physiological saline injection, which was not changed after 6-OHDA injection, except in the clinical type group where it was significantly reduced (
Co-Evolve Serum Metabolic Signatures with PD Progression in 6-OHDA Rats
1H
13C
and C27
.3
.9
6
4
indicates data missing or illegible when filed
The graph of the OPLS scores, made with the NMR data and the PDP scores for each animal (i.e. the «6-OHDA serum OPLS model») shows a clear gradation of the dot colors, from «cold» (white) on the left to «warm» (dark gray) on the right. This indicates that the metabolic profile is changing along with disease progression (
Metabolic Profiles of Brain Tissue in the 6-OHDA Rat Model Strengthen Blood Findings
The «High Resolution Magic Angle Spining» (HRMAS) 1H NMR” spectroscopy was used to study the metabolic profiles of intact brain biopsies (DS and Nacc) of 6-OHDA animals, sampled at the same level as sections used for histology to quantify dopaminergic denervation. A representative spectrum of a DS biopsy from a 6-OHDA rat is shown in
The OPLS analysis constructed with the Nacc data (i.e. the «OPLS 6-OHDA Nacc model») showed the same trend, although less clear (
Although some metabolites such as pyruvate and citrate are not detected in tissue samples, due to their immediate transformation upon death (post-mortem effect), these experiments show dysregulations of alanine and lactate levels common to serum and tissue samples.
Pramipexole Modifies Serum and Tissue Metabolic Dysregulations Induced by 6-OHDA
It was checked whether the serum metabolic profile was altered by a classically used clinical dopamine agonist, Pramipexole (Pra), known to improve only neuropsychiatric symptoms when used at low dose in animal models [23,24]. As expected, after 15 days of subchronic Pra administration, deficits in the self-administration task measured in 6-OHDA rats were fully or partially ameliorated in prodromal-type and clinical-type animals, respectively. On the other hand, the performance of Veh-treated rats was not improved (
With respect to serum metabolic profiles, to visualize simultaneously the effect of 6-OHDA lesion and treatment, metabolite levels of prodromal-like animals were normalized to the level of control (sham) animals (lesion effect, solid bar) and metabolite levels of Pra-treated prodromal-type animals were normalized to Veh-treated prodromal-type animals (Pra effect, hatched bar). In prodromal animals, alanine BHB and pyruvate levels increased after injury (
Unlike serum, for which multiple longitudinal samples could be collected from the same animal, brain samples could only be collected at the end of the experiment, when all animals had received a Pra or Veh chronic treatment. Therefore, control (sham) animals treated with Veh were used as a reference to normalize metabolite levels for prodromal type animals treated with Veh (Lesion effect, solid bar). The metabolite levels of Pra-treated prodromal-type animals were normalized relative to prodromal-type animals given Veh (Pra effect, hatched bar). As in serum, administration of Pra caused a decrease in some metabolites previously increased by 6-OHDA lesion. That was true for alanine in DS and Nacc, and for lactate in Nacc in prodromal type animals (
Overall, significant alterations in serum and tissue metabolites were measured in the 6-OHDA rat model and associated with the different phases of the disease and, interestingly, were detected as early as the prodromal stage. Moreover, some of these disturbances were partially reversed by chronic administration of Pra, which also corrected the neuropsychiatric deficits.
In order to test the specificity of the biomarkers found in the 6-OHDA model regarding disease progression and pathophysiological mechanisms, the serum metabolomics study was extended to two other animal models of PD.
6-OHDA, α-Synuclein Rats and MPTP Monkeys Share Serum Metabolic Disturbances
First, using the same strain of rat (Sprague Dawley) as for the 6-OHDA study, an α-synuclein rat model was selected, in which overexpression of human A53T alpha-synuclein was induced in SNc using adeno-associated viral (AAV) vectors. This targeted expression leads to progressive neurodegeneration in the targeted region, which makes it possible to follow the different stages of PD in the same animal [19], namely the prodromal-type phase around 3 weeks and the clinical-type phase around 10 weeks after the injection.
An OPLS-DA model was used to assess whether a specific metabolic signature was associated with these different phases (i.e. the «serum α-synuclein OPLS-DA model»). Sham animals (control injected with a comparison group virus encoding GFP) showed no difference between the three different time points 0, 3 and 16 weeks after injection (
These results revealed similar metabolic disturbances between the two rat models of PD, which are associated with the disease progression. To go further, a metabolic analysis was carried out in a non-human primate model of PD, e.g. the MPTP model, which mimics the clinical phase of PD and shows greater homology with humans [25].
The OPLS-DA model shown in
Taken together, these results reveal common metabolic alterations associated with PD progression in three different animal models of PD, performed in two different species. Thus, to validate the clinical relevance of these markers, these preclinical results were compared with those obtained in samples from patients with PD.
Serum Metabolic Signatures of Newly Diagnosed Parkinson's Patients and Animal Models of PD Show Similarities
Each cohort was first analyzed individually, ie. newly diagnosed parkinsonian patients (n=21 in each cohort) compared to age- and sex-matched comparison groups (n=30/23 NIH/Italy).
With the exception of some metabolites, e.g. alanine (data not shown), the variations observed were consistent between the 2 cohorts, as shown in Table 3. The latter presents the main individual and collective changes in metabolites between patients with PD and controls in the Italian and United States (NIH) cohorts, analyzed individually by t-test and collectively by two-way ANOVA, followed by a Sidak post hoc test. *: p≤0.05; **: p≤0.01; ***: p≤0.001; ns: not significant.
For example, a marginal increase in BHB and acetoacetate was observed in Italian PD patients (+28.3%, p=0.24 and +17.2%, p=0.05 respectively), and was even higher in patients from the United States (+103.2%, p=0.09 and +28.1%, p=0.12).
In order to increase the statistical power, the two cohorts were pooled (PD n=42, comparison group n=53). Nevertheless, since the variance of the data was mainly affected by the origin of the sample, that is to say Italy or the United States, independently of the case (PD or comparison group), a two-way ANOVA was performed with i) cohort origin and ii) case (comparison group or PD) as factors to extract PD-only information. Table 3 illustrates the significant metabolites for the case factor, with no interaction between the origin cohort and the case. The increase in BHB and acetoacetate, observed in each cohort, appeared significant when pooled, A significant decrease in betaine in the serum of patients with PD was also detected, as in all animal models. In addition, a significant decrease in valine was observed as in the 6-OHDA rat model. Finally, a decrease in leucine was observed. Additionally, the analysis of a small number of serum samples from other newly diagnosed and treated with Pra 0.5-100 mg/day (n=9) revealed a trend towards normalization of BHB and betaine levels compared to untreated PD patients, as observed in rats 6-OHDA (
Serum Metabolites as Potential Markers for the Early Diagnosis of Parkinson's Disease
To assess the potential of the identified metabolites to distinguish prodromal type animals from control (sham) animals, multiple logistic regression models were developed, based on all possible combinations of levels of acetoacetate, betaine, creatine, BHB, pyruvate and valine, grouping prodromal-type 6-OHDA and α-synuclein rats.
The best regression model, evaluated by the ROC curves (see materials and methods), retained the six metabolites, using a regression based on the following predictive algorithm (p: 1.107×10−15):
Logit(P)=log(P/(1−P))=−54.72+0.49BHB+0.39Pyruvate−0.42Valine+0.06Acetoacetate+0.52Creatine−0.10Betaine
The ROC curve had an area under the curve (AUC) of 0.936. By applying the optimal threshold of 0.618 (sensitivity 0.85, specificity 0.88) (
Regarding human samples, the United States (NIH) cohort was used to build a model because i) it is the largest cohort and ii) it includes only recently diagnosed patients, unlike the Italian cohort, which includes patients with longer PD duration. The latter was used as an external validation cohort. Here again, the best model retained the six metabolites combined with the following regression algorithm (p: 1.95×10−4): Logit(P1)=log(P1/(1−P1))=−45.60+0.28BHB−4.54Pyruvate−11.23Valine+7.05Acetoacetate+2.85Creatine−8.27Betaine.
The following ROC curve had an AUC of 0.88. The optimal threshold was 0.316, and the corresponding sensitivity and specificity were 0.957 and 0.714, respectively (
The United States (NIH) cohort minus 2 unretained patients was also used to build a model. Here again, the best model retained the six metabolites combined with the following regression algorithm (p: 1.95×10−4):
Logit(P2)=log(P2/(1−P2))=−42.78+0.16BHB−3.31Pyruvate−12.01Valine+7.06Acetoacetate+3.14Creatine−7.79Betaine.
The following ROC curve had an AUC of 0.87. The optimal threshold was and the corresponding sensitivity and specificity were 0.952 and 0.714, respectively (
Finally, multiple logistic regression was performed using serum data from all MP animals (prodromal and clinical type) from all models using the same six metabolites and was evaluated with the following ROC curve showing an AUC of (
In total, the precision of the biomarker developed from the same six serum metabolites was at least 82.6% regardless of the species.
| Number | Date | Country | Kind |
|---|---|---|---|
| FR2011817 | Nov 2020 | FR | national |
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/FR2021/052035 | 11/18/2021 | WO |
