The present invention relates to the use of molecular markers for detection, prognosis and follow up of Parkinson's disease.
Abbreviations: ACTB: beta actin; AD: Alzheimer's disease; ALAS1: aminolevulinate delta synthase 1; ALDH1A1: aldehyde dehydrogenase 1 family, member A1; ALS: Amyotrophic Lateral Sclerosis; ARPP: 21-cyclic AMP-regulated phosphoprotein; DA: dopamine; DEPC: diethyl pyrocarbonate; EGLN1: eg1 nine homolog 1; EIF4G1: eukaryotic translation initiation factor 4 gamma, 1; GAPDH: glyceraldehyde-3-phosphate dehydrogenase; HSPA8/HSC70/HSC54: chaperone heat shock 70 kDa protein 8; L13a: RPL13A; LGALS9: lectin, galactoside binding, soluble, 9; LOC56920: semaphorin sem2; LRP6: low density lipoprotein receptor-related protein 6; MAN2B1: mannosidase, alpha, class 2B, member1; MPTP: N-methyl-4-phenyl-1,2,3,6-tetrahydropyridine; OS: Oxidative stress; PARVA: parvin, alpha; pc: pars compacta; PD: Parkinson's disease; PENK: proenkephalin; PET: positron emission tomography; PMD: post-mortem delay; ROS: Reactive oxygen species; SELPLG: selectin P ligand; SKP1A: S-phase kinaseassociated protein 1A; SN: Substantia nigra; SPECT: single photon emission tomography; SPHK1: sphingosine kinase 1; SRPK2: SFRS protein kinase 2; SRRM2: serine/arginine repetitive matrix 2; TMEFF1: transmembrane protein with EGF-like and two follistatin-like domains 1; TRIM36: tripartite motif-containing 36; UCHL: 1-ubiquitin C-terminal hydrolase-L1; UPS-ubiquitin-proteasome system; UPDRS: Unified Parkinson's Disease Rating Scale; VMAT: vesicular monoamine member; ZSIG11: putative secreted protein ZSIG11.
Parkinson's disease (PD) is a progressive neurodegenerative disorder, with a prevalence of 1% in the population above 65 years of age, that results in degeneration of dopamine neurons in the substantia nigra (SN), and a consequent striatal dopamine deficiency (Bernheimer et al., 1973). The causes and mechanism for the degeneration of dopaminergic neurons is still elusive. Sporadic (non-hereditary) Parkinson's disease constitutes the most common form of the disease.
Both hereditary and sporadic forms of PD primarily converge in impairment of protein handling, catabolism and in oxidative-stress damage.
Potential contributing factors to the etiology of PD include ongoing selective oxidative stress (OS) resulting from mitochondrial dysfunction, auto-oxidation or enzymatic (monoamine-oxidase) oxidation of dopamine (DA) and excessive iron accumulation in the SN pars compacta (pc) (Riederer et al., 1989; Youdim et al., 1993; Gotz et al., 1994; Jenner and Olanow, 1996; Olanow and Youdim, 1996; Youdim and Riederer, 1997; Jenner, 1998). Specifically, redox-active iron has been observed within the melanin containing neurons that selectively die and in the rim of Lewy body, the morphological hallmark of PD.
Lewy body is composed of lipids, aggregated alpha synuclein (concentrating in its peripheral halo) and ubiquitinated, hyperphosphorylated neurofilament proteins (Jellinger, 2003). A number of studies (Ostrerova-Golts et al., 2000; Ebadi et al., 2001; Turnbull et al., 2001) have shown that alpha synuclein forms toxic aggregates in the presence of iron and this is considered to contribute to the formation of Lewy body via OS. Recently, Lewy bodies with aggregated ubiquitinated proteins are beginning to be viewed as a defensive measure aimed at removing toxic misfolded damaged proteins, rather than a mere general site for precipitated misfolded proteins (Hashimoto et al., 2004; Tanaka et al., 2004).
Misregulation of brain iron metabolism has taken center stage in neurodegenerative diseases since a significant number of mutated iron metabolism genes have now been shown to be directly involved in neurodegeneration (Felletschin et al., 2003; Youdim and Riederer, 2004). Thus, iron redox status constitutes a pivotal factor contributing to the extent of protein misfolding and aggregation in the ageing and disease affected brain.
Defects in ubiquitination and proteasomal protein handling are common features in PD and other chronic neurodegenerative diseases such as Alzheimer's disease (AD), Amyotrophic Lateral Sclerosis (ALS) and Huntington disease, and in ageing (Ciechanover and Brundin, 2003; Dawson and Dawson, 2003). This may, in turn, lead to impairment in several cellular processes linked to ubiquitination such as cell cycle, processing and regulation of transcription, intracellular trafficking, signalling pathways and degradation of normal and damaged intracellular proteins (Ciechanover and Brundin, 2003).
Accumulating evidence points out to a crucial role for protein misfolding and aggregation into protein inclusions in sporadic PD, which constitutes the most common form of the disease. For example, a loss of 20S proteasome alpha-subunits (McNaught et al., 2002; McNaught et al., 2003) and reduction in the activity of the 26/20S proteasome system in SNpc of sporadic PD (McNaught et al., 2003) has also been reported.
Three apparently independent gene mutations in alpha synuclein, parkin and ubiquitin C-terminal hydrolase-L1 (UCHL-1), which are capable of impairing the activity of the ubiquitin-proteasome system (UPS), have been described in rare forms of hereditary PD (Dauer and Przedborski, 2003). More recently, recessive mutations in the DJ-1 gene were proposed to play a role in cellular response to oxidative stress (Bonifati et al., 2003). However, none of these genes have been demonstrated to be mutated in sporadic PD, which constitutes more than 95% of total PD cases.
An approach to explore the etiology of PD is based on the study of changes in the up or down regulation of gene expression, which might increase the vulnerability of the neurons to cell death or even cause it. A number of studies have reported alterations in the expressions of various genes, such as decrease in calcium-binding protein (28 kDa calbindin-D) in the SN (Iacopino and Christakos, 1990) and D3 receptor mRNA in lymphocytes (Nagai et al., 1996) from PD patients.
Among the neurotoxins used to induce parkinsonism in animals, N-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) best replicates many of the neurochemical and anatomical characteristics of the Parkinsonian syndrome in rodents, primates and other species (Dauer and Przedborski, 2003). In spite of certain limitations intrinsic to this model, namely lack of progressive nature and Lewy bodies (Inclusion bodies), much has been learned about the molecular events that lead to dopaminergic neurodegeneration.
Using the MPTP animal model alterations in the expression of several genes were observed. For example, an increase of glutamate decarboxylase mRNA was observed in a subpopulation of neurons in the putamen of parkinsonian monkeys, which provides further evidence that striato-pallidal GABAergic neurons are hyperactive in MPTP-treated parkinsonian monkeys (Soghomonian and Laprade, 1997). An increase in Bax (a cell death effector) mRNA expression in SN, with a concomitant increase in Bax immunoreactivity was also detected in the same model (Hassouna at al., 1996).
More recently, the MPTP model has been employed to asses differential gene expression changes in the midbrain of mice (Grunblatt et al., 2001), to obtain a more global picture of the series of events occurring during degeneration. The median density microarray employed in the study contained around 1200 genes. Thus, a limited picture of the differential gene expressions that may be altered was observed. Nevertheless, the study demonstrated alterations in genes related to oxidative-stress, inflammation, nitric oxide, glutamate excitotoxicity and neurotrophic factors pathways, which were up or down regulated. Additional gene pathways, not described previously, including, cell cycle and iron metabolism regulation, apoptosis, intermediary metabolism and signal transduction have also been observed.
Furthermore, an increase in interleukin-1 beta (IL-1β) mRNA, in the pro-death genes BAD, BAX, and BID was reported in methamphetamine treated rats (Yamaguchi et al., 1991). Concomitantly, there were significant decreases in the anti-death genes Bcl-2 and Bcl-XL (Cadet et al., 2001; Jayanthi et al., 2001).
In another animal model for PD, caused by depletion of dopamine via 6-hydroxydopamine(6-OHDA), an increase in the glutamate decarboxylase mRNA expression was found in adult rats striatum as well as in neonates. A parallel increase in preproenlcephalin and a decrease in preprodynorphin mRNA levels was observed (Laprade and Soghomonian, 1999).
Current accepted clinical criteria for the diagnosis of PD, such as Unified Parkinson's Disease Rating Scale (UPDRS) (Fahn, 1987), provide high sensitivity for detecting parkinsonism (Brooks, 1998). However, there are no sensitive and specific biochemical markers that can be used to reliably diagnose clinical and especially preclinical PD for sporadic and most familial PD cases.
One approach in the diagnosis of PD is functional imaging, which provides a means of discriminating typical from atypical PD, revealing characteristic patterns of loss of dopaminergic function. In addition, positron emission tomography (PET) and single photon emission tomography (SPECT) show preserved levels of striatal metabolism and dopamine receptor binding in PD, whereas levels are reduced in the atypical variants (Brooks, 1998). Still these tools do not give an exact diagnosis of PD, and often, experts in PD changed their diagnoses infrequently during the 7.6-year follow-up (Jankovic et al., 2000).
Moreover, all these diagnostic methods are able to detect subjects with PD only after nearly 70% of the neurons have been degenerated, as only at this point symptoms appears (Zigmond et al., 1989). This, of course, makes treatment and maybe even rescue of the neurons nearly impossible. In view of this, it is desirable to diagnose the disease at an early stage and, for this reason, it is important to develop an early diagnostic method.
There are already many drugs that show neuroprotective effects in vitro and in vivo (Drukarch and van Muiswinkel, 2001; Grunblatt et al.; 2001 and 2003; Koller, 2002; Olanow et al., 1996 and 1998; Riederer et al., 2000 and 2002; Soto-Otero et al., 2002). Additionally, gene therapy was shown to be most successful in delaying the neurodegenerative process (Azzouz et al., 2002; Le and Frim, 2002; Luo et al., 2002; McBride and Kordower, 2002; Monville, 2002; Olanow, 2002; Senior, 2002; Tenenbaum et al., 2002). The problem is that these drugs and methods did not show any success in patients with PD, because the beginning of therapy is too late and the number of surviving neurons is too small. Therefore, an early diagnosis may provide a better time point for the submission of therapeutic strategies that can protect against the cell death occurring in PD and to prevent the progress of the disease.
Therefore, it is crucial to treat PD patients before substantial damage of neurons occurs, and typical symptoms appear. Thus, there is a need for an early diagnostic method for detection and monitoring of PD at an early stage.
The present invention relates, in one aspect, to the use of molecular markers for Parkinson's disease for the purpose of diagnosis, monitoring disease progression and effect of treatment.
In another aspect, the present invention relates to methods for diagnosis, prognosis and follow up of Parkinson's disease comprising the use of measurable biological markers, wherein said markers are one or more genes exhibiting an altered expression pattern, or gene products encoded thereby.
The genes for use as molecular markers according to the invention are genes which expression is upregulated or downregulated in PD patients, as detailed hereinafter.
The present invention relates to the use of specific molecular or biological markers for detection, prognosis and follow up of treatment of Parkinson's disease.
According to the present invention, genes with altered expression pattern in PD were found, such as for example the genes ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, PSMC4, HIP2, PACE4, COX6A1, PFKP, OXCT, GBE1, UQCRC2, LANCL1, TRIP15, PIK3CA, PLCL1, GNG5, GNAI1, VEGF, RHOB, NR4A2, SCL31A2, SCP2, PIGH, ARIH2, GMPR2, PP, IKBKAP, PRKACB, PTPRN2, BCAS2, IARS, PPP1R8, SEP15, TAF9, ZFP103, WRB, TMEM4, SMARCA3, FMR1, PDE6D, SGCE, AUH, SLC16A7, ATP6V1E1, UGTREL1, SEC22L1, CD9, CDH19, DUSP1, HSA6591, ACTR3, KIF2, TUBB2, ASPA, HELO1, C3orf4, CBR1, XPOT, LOC51142, NY-REN-45, SET0-2, EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, ZSIG11, ITGB3BP, ITGAM, COL18A1, TM4SF9, LAMB2, HS3ST2, TSTA3, COL5A3, PALM, MYOM1, FLNB, HMBS, KRT2A, CSK, NUDC, HYPE, GAK, SIAT1, CSF1R, ICSBP1, CD22, ERCC1, DNAJB5, TRAF3, MMP9, EIF4G1, RPL36, SRPK1, CSNK1G2, RPS6KA1, JIK, LNK, INPP5D, TCOF1, NAPG, SLC19A1, ITSN1, LOC51035, PMVK, C21orf2, EFEMP2, TBL1X, APRT, SPUF, GLTSCR2, ADIR, PSCD4, CBFA2T1, CUGBP1, ING4, STAT6, ZNF239, TAL1, TAF11, MXD4, RDHL, LOC51157, LRP6, MBD3, and C9orf7.
Thus, in one embodiment, the invention relates to the use of molecular markers for the detection, prognosis and follow up of treatment of Parkinson's disease, wherein the molecular markers are one or more of the above genes or their gene products.
It was further found by the inventors that some of the above mentioned genes show an increased expression in PD patients while other genes show a decreased expression, as compared to healthy individuals. Tables 3 and 4 hereinafter present the downregulated and upregulated genes, respectively, and their characterization including, for each gene, the Gene Bank number, the full gene name and its function.
In accordance to the present invention, the following genes were found to be upregulated, namely, to exhibit an increased level of expression in PD, and will be referred to hereinafter as “the upregulated genes”: EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, ZSIG11, ITGB3BP, ITGAM, COL18A1, TM4SF9, LAMB2, HS3ST2, TSTA3, COL5A3, PALM, MYOM1, FLNB, HMBS, KRT2A, CSK, NUDC, HYPE, GAK, SIAT1, CSF1R, ICSBP1, CD22, ERCC1, DNAJB5, TRAF3, MMP9, EIF4G1, RPL36, SRPK1, CSNK1G2, RPS6KA1, JIK, LNK, INPP5D, TCOF1, NAPG, SLC19A1, ITSN1, LOC51035, PMVK, C21orf2, EFEMP2, TBL1X, APRT, SPUF, GLTSCR2, ADIR, PSCD4, CBFA2T1, CUGBP1, ING4, STAT6, ZNF239, TALI, TAF11, MXD4, RDHL, LOC51157, LRP6, MBD3, and C9orf7.
Further in accordance with the present invention, the following genes were found to be downregulated, namely, to exhibit a decreased level of expression in PD, and will be referred to hereinafter as “the downregulated genes”: ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, PSMC4, HIP2, PACE4, COX6A1, PFKP, OXCT, GBE1, UQCRC2, LANCL1, TRIP15, PIK3CA, PLCL1, GNG5, GNAI1, VEGF, RHOB, NR4A2, SCL31A2, SCP2, PIGH, ARIH2, GMPR2, PP, IKBKAP, PRKACB, PTPRN2, BCAS2, IARS, PPP1R8, SEP15, TAF9, ZFP103, WRB, TMEM4, SMARCA3, FMR1, PDE6D, SGCE, AUH, SLC16A7, ATP6V1E1, UGTREL1, SEC22L1, CD9, CDH19, DUSP1, HSA6591, ACTR3, KIF2, TUBB2, ASPA, HELO1, C3orf4, CBR1, XPOT, LOC51142, NY-REN-45, SET0-2.
As used herein the terms “molecular markers” and “biological markers” are used interchangeably to include gene products that can be identified and measured for diagnosing PD in samples obtained from individuals.
The term “gene products” as used herein refers to the expression product of a gene by conversion of the information contained in a gene into a gene product. A gene product can be the direct transcript of the gene, i.e., an RNA such as mRNA, tRNA, or any other type of RNA, or a protein produced by translation of a mRNA.
As used herein, the term “genes with altered expression pattern in PD” refers to genes that are upregulated or downregulated in PD patients in comparison, e.g., to healthy individuals. The term “upregulation” as used herein in the context of upregulation of a gene refers to any process which results in an increase in production of a gene product, e.g. RNA or protein expressed by the gene. The term “downregulation” as used herein in the context of downregulation of a gene refers to any process which results in a decrease in production of a gene product, e.g. RNA or protein expressed by the gene.
The level of gene expression, including the level of gene upregulation or downregulation, can be measured by known procedures as described hereinafter. In general, gene upregulation or downregulation comprises any detectable increase or decrease in the production of a gene product. The terms upregulated and gene activation can also mean that the observed activity relative to a baseline level is a statistically significant difference (i.e., increase or decrease). A difference is typically considered to be “statistically significant” if the probability of the observed difference occurring by chance (the p-value) is less than some predetermined level. As used herein a “statistically significant difference” refers to a p-value that is <0.05, preferably <0.02 for the downregulated genes (Table 3), and more preferably <0.005 for the upregulated genes (Table 4). In accordance with the present invention, the increase of expression in the upregulated genes is at least 1.5-fold, preferably 1-3 fold (Table 4) and the decrease of expression in the down-regulated genes is at least 0.66-fold or below, preferably 0.3-0.66 (Table 3).
In one preferred embodiment of the invention, the biological markers for PD are one or more downregulated genes selected from the group consisting of ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, PSMC4, and gene products thereof, and/or one or more upregulated genes selected from the group consisting of EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, ZSIG11, and gene products thereof.
In one embodiment, the invention relates to a method for diagnosis, prognosis and/or follow up of treatment of Parkinson's disease comprising the use of molecular markers, wherein the molecular markers are one or more of the upregulated and/or downregulated genes of the invention. Preferably, the molecular markers are one or more downregulated genes selected from the group consisting of ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, PSMC4, and gene products thereof, and/or one or more upregulated genes selected from the group consisting of EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, ZSIG11, and gene products thereof.
The methods of the invention comprise obtaining a sample from the individual to be tested and detecting in the sample decreased or increased levels of one or more gene products specific to the upregulated and/or downregulated genes of the invention.
The tested individual may be an individual suspected of having PD, an individual exhibiting Parkinsonian-like symptoms, a de novo PD patient, a PD patient under treatment with a neuroprotective or another suitable PD drug, to follow up progress of the disease and/or to monitor effectiveness of the treatment. Thus, the invention provides methods for diagnosing PD at an early stage and for monitoring treatment of PD patients, both familial and sporadic PD patients.
Any suitable sample can be used for detection, but preferably blood, serum or biopsy samples of skin obtained from the tested individual are used.
A variety of different assays can be utilized to detect alteration in gene pattern expression, including methods that detect gene transcript levels and gene encoded protein levels. More specifically, the diagnostic and prognostic methods disclosed herein involve obtaining a sample from an individual and determining the expression level of one or more of the upregulated and/or downregulated genes in the sample. Usually this determined value or test value is compared against some type of reference such as the expression level in control (healthy individuals, non-PD patients) or a baseline value.
In accordance with the invention, also the use of the molecular markers in vivo is contemplated.
When the gene product to be tested is RNA extracted from the sample, the level of gene expression, i.e the level of gene upregulation or downregulation, can be measured by known procedures including, but not limited to, nucleic acid probe arrays, Northern blots, RNase protection assays (RPA), quantitative reverse-transcription PCR(RT-PCR), dot blot assays and in-situ hybridization.
In one preferred embodiment, nucleic acid probe arrays are used to detect and quantitate gene transcript (as described in the Examples below). The arrays can be of different types and may include probes of varying types such as, for example, short-length synthetic probes (20-mer or 25-mer), full length cDNA or fragments of gene), amplified DNA, fragments of DNA (generated by restriction enzymes, for example) and reverse transcribed DNA. The array may be a custom array, including probes that hybridize to particular preselected subsequences of mRNA gene sequences of the genes or amplification products thereof, or a generic array, designed to analyze mRNAs irrespective of sequence.
In the method using probe array, nucleic acids obtained from a test sample are usually reversed transcribed into labeled cDNA, although labeled mRNA can be used directly. The test sample containing the labeled nucleic acids is then contacted with the probes of the array and, after any labeled nucleic acid related to the tested genes present in the sample hybridizes to the probes, the array is typically subjected to one or more high stringency washes to remove unbound nucleic acids and to minimize nonspecific binding to the nucleic acid probes of the arrays. Binding of labeled nucleic acid is detected using any of a variety of commercially available scanners and accompanying software programs. For example, if the nucleic acids from the sample are labeled with fluorescent labels, hybridization intensity can be determined by, for example, a scanning confocal microscope in photon counting mode. The label can provide a signal that can be amplified by enzymatic methods or other labels can be used including, for example, radioisotopes, chromophores, magnetic particles and electron dense particles.
Those locations on the probe array that are hybridized to labeled nucleic acid are detected using a reader as commercially available. For customized arrays, the hybridization pattern can then be analyzed to determine the presence and/or relative amounts or absolute amounts of known mRNA species in samples being analyzed.
In another preferred embodiment, real time reverse-transcription PCR (real time RT-PCR) methods can be used to determine the quantity of the gene mRNA present in the sample (see Example below). These methods involve measuring the amount of amplification product formed during an amplification process, for example, by fluorogenic nuclease assays, to detect and quantitate specific transcripts of the genes of interest. These assays continuously measure PCR product accumulation using a dual-labeled fluorogenic oligonucleotide probe as in the approach frequently referred to in the literature simply as the TaqMan® method.
The probe used in real time PCR assays is typically a short (ca. 20-25 bases) polynucleotide that is labeled with two different fluorescent dyes: a reporter dye at the 5′ terminus of the probe and a quenching dye at the 3′ terminus, although the dyes can be attached at other locations on the probe as well. For measuring a specific transcript, the probe is designed to have at least substantial sequence complementarity with a probe binding site on the specific transcript. Upstream and downstream PCR primers that bind to regions that flank the specific transcript are also added to the reaction mixture for use in amplifying the nucleic acid.
When the probe is intact, energy transfer between the two fluorophores occurs and the quencher quenches emission from the reporter. During the extension phase of PCR, the probe is cleaved by the 5′ nuclease activity of a nucleic acid polymerase such as Taq polymerase, thereby releasing the reporter dye from the polynucleotide-quencher complex and resulting in an increase of reporter emission intensity that can be measured by an appropriate detection system. The fluorescence emissions created during the fluorogenic assay is measured by commercially available detectors that comprise computer software capable of recording the fluorescence intensity of reporter and quencher over the course of the amplification. These recorded values can then be used to calculate the increase in normalized reporter emission intensity on a continuous basis and ultimately quantify the amount of the mRNA being amplified.
In the assay based on dot blots and in-situ hybridization, a sample from an individual being tested for PD is spotted on a support, e.g., a filter, and then probed with labeled nucleic acid probes that specifically hybridize with nucleic acids derived from one or more of the upregulated or downregulated genes of the invention. After hybridization of the probes with the immobilized nucleic acids on the filter, unbound nucleic acids are rinsed away and the presence of hybridization complexes detected and quantitated on the basis of the amount of labeled probe bound to the filter.
Thus, according to the invention, gene expression pattern is determined by extraction of RNA from biological material, preferably blood or biopsy samples of skin. The RNA is isolated rapidly by a commercially available kit. The RNA is then tested through hybridization to a customized GeneChip array containing the selected genes and relevant house-keeping genes serving for normalization, or by means of real time PCR for each of the selected genes. The gene expression pattern is determined via comparison to the expression of positive and negative control RNA (with de-novo PD and healthy subjects, respectively). The pattern of the gene expression obtained via one of the techniques should be similar to the pattern described in Table 3 or 4 in order to define the subject as PD patient.
In another embodiment, the gene product obtained by expression of the upregulated or downregulated gene is a protein, that can be detected by an antibody, or a fragment thereof, able to bind to that protein. The antibody, or a fragment derivative, may be detectably labeled with any appropriate marker, for example, a radioisotope, an enzyme, a fluorescent label, a paramagnetic label, or a free radical.
In a preferred embodiment, the invention relates to a method for diagnosing occurrence of Parkinson's disease in an individual exhibiting Parkinsonian-like symptoms, comprising detecting in a sample obtained from said individual a decreased level of expression of one or more of the downregulated genes, more preferably one or more of the genes ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, and PSMC4, and/or an increased level of expression of one or more of the upregulated genes, more preferably one or more of the genes EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, and ZSIG11. An increase of the expression of one or more of the upregulated genes by at least about 1.5-fold and/or a decrease of the expression of one or more of the downregulated genes to about 0.66-fold or below will indicate that the individual suffers from PD.
In comparison with the current methods used for diagnosis of PD, this method has the advantage not to use rare genetically mutations or familial history of the disease, but rather use general gene expression changes which occur also in sporadic PD. These gene expression alterations may be caused not only as a consequence of specific genetic background, but also as a consequence of environmental background. Therefore, this method will detect both PD patients carrying gene mutations and sporadic PD patients very early in the development of the disease, before extensive cell death and irreversible substantial damage of neurons occurs.
When untreated, PD progresses to total disability, often accompanied by general deterioration of all brain functions, and may lead to an early death. Treated, the disorder still impairs people in varying ways. Most people respond to some extent to medications. The extent of symptom relief, and how long this control of symptoms lasts, is highly variable.
Thus, in another embodiment, the invention provides a method for prognosticating or monitoring the treatment of a PD patient, which comprises detecting in a sample obtained from said patient the level of expression of one or more of the downregulated or upregulated genes of the invention, whereby an increased level of expression of one or more of the downregulated genes, more preferably of one or more of the genes ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, and PSMC4, and/or a decreased level of expression of one or more of the upregulated genes, more preferably one or more of the genes EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, and ZSIG11, indicate the effectiveness of the treatment.
As shown in the Examples below, the inventors identified global major differentially expressed genes in the most affected brain area in PD, the substantia nigra pars compacta (See Tables 3 and 4). Employing Affymetrix high density DNA microarray, differentially expressed genes compared to aged-matched controls were identified. The cerebellum of PD, an unaffected brain region, served as control for tissue specificity. Confirmation of gene expressions was achieved by analysis with quantitative real-time PCR.
It was found that the downregulated genes belong to signal transduction, protein degradation (e.g. ubiquitin-proteasome subunits), dopaminergic transmission/metabolism, ion transport, protein modification/phosphorylation and energy pathways/glycolysis functional classes. The upregulated genes clustered mainly in biological processes involving cell adhesion/cytoskeleton, extracellular matrix components, cell cycle, protein modification/phosphorylation, protein metabolism, transcription and inflammation/stress (e.g. key iron and oxygen sensor EGLN1).
One major finding was the particular decreased expression of SKP1A, a member of the SCF (E3) ligase complex specifically in the substantia nigra (SN) of sporadic parkinsonian patients. SKP1A decline was accompanied by decreased expression in various subunits of the 26S proteasome, in energy pathways and signal transduction, in parallel with a prominent increase in a number of genes with functional activities related to cell adhesion/cytoskeleton, extracellular matrix components and inflammation/stress.
SKP1 is part of the Rbx family of RING proteins (Kamura et al., 1999) functioning within modular multiprotein Skp1, Cullin and a substrate-recognizing F-box protein (SCF). This unit allows the formation of multiple E3 complexes, which, in turn, are able to recognize a wide spectrum of different protein substrates.
SCF complexes are modular: SKP1 can interact with several F-box proteins, which are responsible for specific target recognition, thereby providing functional diversity and increasing the repertoire of proteins processed by this complex. Humans express only one functional SKP1 isoform (Semple, 2003). Thus, the decrease in its expression observed, may constitute a rate limiting factor and may account for the accumulation of a wide spectrum of ubiquitinated protein aggregates in brains of PD patients such as tyrosine hydroxylase, synphilin-1, alpha-synuclein, phosphorylated tau (Liani et al., 2004; Meredith et al., 2004; Zhang and Goodlett, 2004).
Quantitative real time PCR analysis revealed that the SN expresses extremely low mRNA levels of SKP1A in both normal and PD subjects. For example, control samples display very low expression values (0.013±0.011), differing by more than two orders of magnitude from Rab3B (1.80±0.70), whose expression was not affected in SNpc of PD. A clear trend of decreased expression of SKPA1A in the SNpc of PD was observed. The expression of SKP1A was not affected in the SNr or in cerebellum of PD.
In addition to this selective decrement in SKP1A, a concomitant decrease in the 20S proteasome subunits alpha-5 (PSMA5), alpha-3 (PSMA3) and alpha-2 (PSMA2) mRNAs and in two subunits of the 19S regulatory complex of the 26S proteasome, the non-ATPase subunit 8 (PSMD8/Rpn12) and the ATPase subunit 4 (PSMC4/TBP7/Rpt3), was observed in the SN of PD brains, which may further contribute to dopaminergic neuronal damage in PD.
These findings reveal novel players in the neurodegenerative scenario and provide potential targets for the development of novel drug compounds.
The present invention thus further provides a method for screening for an agent useful for treating Parkinson's disease, which comprises identifying an agent that upregulates the expression of one or more of the downregulated genes, preferably the downregulated genes ALDH1A1, ARPP-21, HSPA8, SKP1A, SLC18A2, SRPK2, TMEFF1, TRIM36, ADH5, PSMA3, PSMA2, PSMA5, PSMC4, and/or downregulates the level of expression of one or more of the upregulated genes, preferably EGLN1, EIF4EBP2, LGALS9, LOC56920, LRP6, MAN2B1, PARVA, PENK, SELPLG, SPHK1, SRRM2, ZSIG11.
High throughput screening methods as known in the art can be used to identify compounds that will be useful for treatment of PD.
The invention will now be illustrated by the following non-limiting examples.
Post-mortem human brains were from the National German Brain-Net Center (Project no GA10) and the Netherlands Brain-Bank (Project no. 350). PD tissues were obtained from moderately to severe parkinsonism individuals based on the Hoehn & Yahr criteria in accordance with established guidelines (Hoehn and Yahr, 1967). All the subjects were negative for AD pathology according to Braak & Braak (Braak and Braak, 1997) (for details see Table 1). The average age for PD and control is 76.6 and 77.8 years, respectively. The average post-mortem delay (PMD, time interval from death to brain-freezing at −80° C.) for PD and control is 26.2 and 19.8 hours, respectively. At autopsy, brains were dissected; one half was snap frozen in liquid nitrogen as slices and the other half was stored in neutral buffered formalin for histopathology. Stringent criteria were used in the case selection of human brain tissues employed in these studies. Human SN and cerebellum regions were obtained from control and PD donors. All procedures were in accordance with the with consensus criteria established by the German and Netherlands Brain-Bank systems and in accordance NIH Guide for the Care and Use of Laboratory Human Tissue and approved by the University of Würzhurg Ethics Committee (Würzburg, Germany).
Total RNA was prepared with a phenol-guanidine isothiocyanate reagent (peqGOLD TriFast; PeQLabGmbH, Erlangen, Germany) (Lukiw and Bazan, 1997). RNA isolation reagents were prepared with 0.2 μM filtered diethyl pyrocarbonate (DEPC)-treated water (Fermentas Inc., Hanover, Md., USA) and used throughout the isolation procedure. Total RNA samples were spectrophotometrically scanned from 220 to 320 nm; the A260/A280 of total RNA was typically >1.9. In addition, formaldehyde agarose gel electrophoresis was conducted for quality control of total RNA. For all total RNAs extracted, the 28S/18S ratio was >1.5. Importantly, no significant differences in the spectral purity, rate of degradation, molecular size, or yield of SN and cerebellum total RNA between PD and control, were noted. Total RNAs was subjected to DNase-I digestion to get rid of genomic DNA residues and subsequently cleaned by the RNeasy Mini Kit (RNeasy Mini Kit; Qiagen Inc., Valencia, Calif., USA).
(iii) Array Processing
All experiments were performed using Affymetrix HG-FOCUS oligonucleotide arrays, as described at //www.affymetrix.com/support/technical/datasheets/human_datasheet.pdf. Total RNA from each sample was used to prepare biotinylated target RNA, with minor modifications from the manufacturer's recommendations (http://www.affymetrix.com/support/technical/manual/expression_manual.affx). Briefly, 10 μg of mRNA was used to generate first-strand cDNA by using a T7-linked oligo(dT) primer. After second-strand synthesis, in vitro transcription was performed with biotinylated UTP and CTP (Enzo Diagnostics), resulting in approximately 100-fold amplification of RNA. A complete description of procedures is available at: http://bioinf.picr.man.ac.uk/mbcf/downloads/GeneChip_Target_Prep Protocol_CR UK_v—2.pdf. The target cDNA generated from each sample was processed as per manufacturer's recommendation using an Affymetrix GeneChip Instrument System (http://www.affymetrix.com/support/technical/manual/expression_manual.affx). Briefly, spike controls were added to 15 μg fragmented cDNA before overnight hybridisation. Arrays were then washed and stained with streptavidin-phycoerythrin, before being scanned on an Affymetrix GeneChip scanner. A complete description of these procedures is available at: http://bioinf.picr.man.ac.uk/mbcf/downloads/GeneChip_Hyb_Wash_Scan_Protocol_v—2_web.pdf. Additionally, quality and amount of starting RNA was confirmed using an agarose gel. After scanning, array images were assessed by eye to confirm scanner alignment and the absence of significant bubbles or scratches on the chip surface. BioB spike controls were found to be present on all the arrays, with BioC, BioD and CreX also present in increasing intensity. A total of 17 gene chips were employed (12 for the SNpc and 5 for the cerebellum of PD patients and age-matched controls).
Genes were filtered using Mas 5 algorithm results. One expected difficulty when working with post mortem samples is the various degrees of degradation in the RNA preparations. Therefore, cRNA originating from these samples contain more 3′ ends than 5′ prime ends. Since the Affymetrix array probe sets are designed so that the 3′ ends are selected whenever possible (http://www.affymetrix.com/support/technical/technotes/hgu133_design_technote.p df), the analysis of relatively high-3′ content samples becomes feasible. Probe sets detected by MAS 5 as Present (P) have a signal with a p-value lower than 0.04 (determined from the probes contained in the probe set). This allows us to use present call with confidence.
Samples that are partially degraded may hybridize to part of the probes in the probe sets leading to a high p-value (and detected as absent). If the probe set in at least four out of six samples was detected as present and all of the signals were greater than 20, we determined this probe set to be present. A list of 3517 probe sets representing genes with signals higher than 20 and detected as present (P) in all control samples, or with signals higher than 20 in all PD samples and detected as present in four out of the six samples was generated from the 8763 probe sets contained on the array (supplement http://eng.sheba.co.il/genomics).
262 probe sets differentiated between patient and control samples as determined by the Wilcoxon rank sum test (p-value<0.05) (supplement http://eng.sheba.co.il/genomics). Probe sets were further filtered selecting those where the ratio of the average (geomean) signals of the PD samples and the average signals of the control samples exceeded 1.5 or was lower than 0.66. The Bonferroni's correction has not been applied because it would result in loss of most of the valid hits. However the use of the aforementioned cut-offs together with the statistical test provide a high stringency analysis and the possibility to focus on defined subset of genes engaged in a common biological process(s). 69 probe sets were up regulated and 68 were down regulated.
In order to validate the microarray results, we performed quantitative real-time RT-PCR for mRNA samples from PD and controls in the SNpc, pars reticulate (SNr) and cerebellum. Total RNA (1-0.4 mg) from each sample was reverse transcribed with random hexamer and oligodT primers using iScript™ cDNA Synthesis Kit (BioRad Laboratories, Hercules, Calif., USA, 170-8890). The genes were normalized to the house-keeping genes: beta actin (ACTB), ribosomal protein L13a (RPL13A), aminolevulinate delta synthase 1 (ALAS1) and glyceraldehyde-3-phosphate dehydrogenase (GAPDH) (QuantiTect™ Gene Expression Assay, Qiagen Inc., Valencia, Calif., USA, Hs ACTB Assay 241013 & Hs GAPD Assay 241011), chosen after analysis according to the program geNorm (vs. 3.3; to be download from: Intp://medgen31.ugent.be/jvdesomp/genorm/) (Vandesompele et al., 2002; Schulz et al., 2004). The geNorm program determined the most stable house-keeping genes from a set of tested genes in a given cDNA sample panel and calculates a gene expression normalization factor for each tissue sample based on the geometric mean of a user-defined number of house keeping genes. Originally, we tested a total of six house-keeping genes, including ribosomal 18S and cyclophilin A, but the latter two were found not suitable.
These house-keeping genes were tested for their stability and found be the most accurate for normalization. Absence of DNA contamination was verified by amplifying the house-keeping gene ribosomal 18S, and running the reaction solution on agarose gel to confirm the absence of product. Minus RT samples tested simultaneously with experimental samples by quantitative RT-PCR with FAM consistently yielded no amplification below 35 cycles using the above protocol. Real-time PCR was performed in the iCycler iQ system (BioRad Co., Hercules, Calif., USA) as described previously (Svaren et al., 2000; Ugozzoli et al., 2002). Briefly, 30-100 ng of cDNA and gene specific primers and probes produced by QuantiTect™ Custome Assay (Qiagen Inc., Valencia, Calif., USA,) (Table 2) were added to QuantiTect Probe PCR Master Mix (Qiagen Inc., Valencia, Calif., USA, 204343). Real-time PCR was subjected to PCR amplification (1 cycle at 95° C. for 15 min, 30-45 cycles at 94° C. C for 15 s, annealing and detecting with FAM at 56° C. for 30 s and extension at 76° C. for 30 s). All PCR reactions were run in duplicate. The amplified transcripts were quantified using the comparative CT method analyzed with the BioRad iCycler iQ system program. Standard curves for each amplification product were generated from 10-fold dilutions of pooled cDNA amplicons, isolated from agarose gel using MinElute™ Gel Extraction Kit (Qiagen Inc., Valencia, Calif., USA), to determine primer efficiency and quantization. Data was analyzed with Microsoft Excel 2000 to generate raw expression values. The differences in gene expression in the different brain regions were compared using the analysis of variance (ANOVA), the StatView software program (Stat View 5.0. software, SAS Institute Inc. Cary, N.C., USA) on a PC computer.
Paraffin-embedded substantia nigra from parkinsonian patients and age-matched controls underwent serial coronal sections (6 micrometer thickness) by microtome. Deparaffinized sections were unmasked in digital decloacking chamber (Biockare Medical, Walnut Creek, Calif., USA) using citrate buffer (pH 6.8), blocked with 10% nalinal goat serum for 2 h at RT and incubated overnight at 4° C., with anti-SKP1. Detection was obtained by appropriate biotinylated second antibody with streptavidin-peroxidase conjugate and S-(2-aminoethyl)-L-cysteine (AEC) as substrate, using Histostain-Plus kit (Zymed, South San Francisco Calif., USA) according to manufacture instructions, Counterstaining was performed with hematoxylin.
DNA microarray analysis was performed on postmortem SNpc from 6 parkinsonian patients and 6 aged-matched controls to identify alteration in gene transcription associated with the disease.
We confined our analysis to genes with at least a 1.5 fold differential expression and significance level of p<0.05. Of a total of 3517 valid probe sets analyzed (see Materials and Methods), 137 met the criteria. A heat map was performed depicting the relative expression levels of these 137 genes in PD and control samples (
Genes are clustered by their relative expression levels over the 12 samples. For negative control, five oligonucleotide hybridizations were performed in a non-related brain area, the cerebellum of both PD (2 samples) and control (3 samples), to assess tissue specificity of gene changes. Hierarchical clustering showed that both PD and normal cerebellum display a complete different pattern of expression (data not shown).
Afterwards, the genes were classified into functional groups according to GeneOntology annotation tools (Dennis et al., 2003) (http://apps 1 niaid nih.gov/David/upload.asp). The Program allows finding regulation trends in group of genes organized according to biological processes, molecular function or cellular components, as defined by the Gene Ontology Consortium (hap://www geneontology.org).
It is recognized that a given gene may be assigned to more than one function or biological pathway. This may result in overestimation of the true size of each functional group. We overcame this by assignment to each gene a single functional class. As shown in
With respect to differentially upregulated genes in PD, they clustered mainly in biological processes involving cell adhesion/cytoskeleton, extracellular matrix components, cell cycle, protein modification/phosphorylation, protein metabolism, transcription and inflammation/stress (
Statistical analysis of gene categories, to identify genes that are the most over-represented or enriched, was accomplished with the EASE program (Hosack et al., 2003). Such groupings increase confidence in the results when the proportion of genes that change within a functional group is significantly greater than the fraction of such genes in the whole chip.
Functional classification with “Ease score” lower than 0.05 were marked as over represented. This analysis revealed that the most prominent changes in the differentially downregulated genes occurred in proteolysis and peptidolysis (10%) biological process, with an EASE score of p<0.02 (Table 3, marked by an asterisk) and in molecular functions related to ion transporter and hydrolase activities (p<0.05, data not presented). Within this proteolytic group we observed decreased expression of both catalytic and regulatory subunits of the UPS such as 20S proteasome subunits alpha-5 (PSMA5), alpha-3 (PSMA3) and alpha-2 (PSMA2), two subunits of the 19S regulatory complex of the 26S proteasome, the non-ATPase subunit 8 (PSMD8/Rpn12) and the ATPase subunit 4 (PSMC4/TBP7/Rpt3). PSMD8 slightly exceeded the upper limit of significance, but a clear trend of reduction was seen in the six patients (three of them >1.5 fold of control). Since PSMD8 clusters together with the other UPS components, we considered it as functionally affected in PD.
On the other hand, statistical significant trends in biological processes and molecular functions with increased gene expression in PD were found for cell adhesion/cellular matrix (˜15%) with an EASE score of p<0.005 (Table 4, marked by an asterisk) and for structural molecule and protein kinase activities (p<0.05, data not shown).
In parallel, an exceptionally stringent analysis was done selecting those genes that display a 1.5 fold differential expression in at least five out of the six parkinsonian SN. Each patient sample signals were compared to the mean of the control sample signals (geomean). This evaluation is of major importance, since it allows identifying variations in gene expression within each of the post-mortem tissue samples that may derive from a combination of factors such as age, cause of death, gender, post-mortem interval, severity of the disease, etc.
This approach restricted the total differentially expressed genes in the PD cases to only 20, of which 8 were decreased and 12 were increased respective to the control group (
Importantly, of the 8 downregulated genes, 3 (37.5%) belong to DA neurotransmission and metabolism (
The 12 upregulated genes (
To confirm our findings of decreased expression of one of the major genes SKP1A, playing a key role in the UPS, we conducted real-time quantitative PCR for it and for two additional genes, HSC70 and VMAT2, all of them affected in at least five out of six parkinsonian patients by a factor <0.67.
We added two more samples to the experimental groups (Control 7 and Patient 7) in cases where the RNA used for the array hybridization was not adequate for the PCR. We analyzed three separate areas: the SNpc, SNr and the cerebellum. Cerebellum was chosen as an area not related to PD, for assessing specificity of gene changes.
Real-time PCR analysis revealed that the SN expresses extremely low mRNA levels of SKP1A in both normal and PD subjects (
We confirmed significant decreased gene expression of HSP8, coding for HSC70/HSC54, both in the SNpc and SNr of parkinsonian patients (FIGS. 4D,E), whereas no significant alterations were observed in cerebellum from PD and controls (
In addition to these novel findings, our array analysis also confirmed evidence for decreased expression of the vesicular DA transporter (VMAT2). Real-time PCR revealed a specific reduction of VMAT2 mRNA in the SNpc while the changes in the SNr were significantly less pronounced (FIGS. 4G,H). These alterations were specific for the SN since no significant difference was seen in the cerebellum between the PD and control samples (
To determine whether the reduced expressions found in both the arrays and in real-time PCR analysis was merely a reflection of a general reduction in mRNA levels in PD SN, or loss of synapses, we investigated the expression of three genes: Rab3b GTPase (FIGS. 4J,K,L), syntaxin 6 and the coatomer protein complex subunit zeta 2 (COP zeta 2) (data not presented), which are regulators of vesicle trafficking from endoplasmic reticulum to Golgi, vesicle membrane docking and neuronal exocytosis of hormones and neurotransmitters (Gonzalez and Scheller, 1999; Futatsumori et al., 2000; Wendler and Tooze, 2001). Their expression did not differ between PD and control samples in the different brain areas.
The stringent analysis performed in the experiments above, focusing on selected biological meaningful pathways in PD and the substantial agreement of gene changes reported in literature and confirmed by real-time PCR analysis. make the analysis of the present results valid.
The immunohistology analysis of SKP1 protein was carried out in samples of melanized dopaminergic neurons of the SNpc from parkinsonian brain and healthy control (as described in Materials and Methods). The results summarized in
The results of the present invention show, for the first time, the global gene expression changes in the SNpc of post-mortem parkinsonian brains as compared to aged-matched controls. One major novel finding is the negative regulation of an essential component of protein catabolism, the SKP1A gene. Its decline was accompanied by decreased expression in various subunits of the 26S proteasome, in energy pathways and signal transduction, in parallel with a prominent increase in a number of genes with functional activities related to cell adhesion/cytoskeleton, extracellular matrix components and inflammation/stress.
Reduction in SKP1A, a Component of the SCF Ubiquitin Ligase Multiprotein Complex, in Parkinsonian SNpc
SKP1 is part of the Rbx family of RING proteins (Kamura et al., 1999) functioning within modular multiprotein Skp1, Cullin and a substrate-recognizing F-box protein (SCF). This unit allows the formation of multiple E3 complexes, which, in turn, are able to recognize a wide spectrum of different protein substrates. SCF complexes are modular: SKP1 can interact with several F-box proteins, which are responsible for specific target recognition, thereby providing functional diversity and increasing the repertoire of proteins processed by this complex. Humans express only one functional SKP1 isoform (Semple, 2003). Thus, the decrease in its expression observed herein, may constitute a rate-limiting factor and may account for the accumulation of a wide spectrum of ubiquitinated protein aggregates in brains of PD patients such as tyrosine hydroxylase, synphilin-1, synuclein, phosphorylated tau (Liani et al., 2004; Meredith et al., 2004; Zhang and Goodlett, 2004).
Quantitative real time PCR analysis revealed that the SN expresses extremely low mRNA levels of SKP1A in both normal and PD subjects. For example, control samples display very low expression values (0.013±0.011), differing by more than two orders of magnitude from Rab3B (1.80±0.70), whose expression was not affected in SNpc of PD. It is also probable that the mRNA of SKP1A is particularly unstable and it degrades fast in post-mortem tissue. This may explain the great variation in the expression levels of the samples, although a clear trend of decreased expression in the SNpc of PD was observed. Because of the special design of the Affymetrix chip, even samples that are partially degraded, such as in post-mortem, may in part hybridize to the probes in the probe sets leading to a high p-value and considered as present. Thus, the chip strategy together with the stringent analysis (p<0.05, decrease in five patients by a factor of at least 1.5), make this result confident. The low SKP1A expression levels may explain, at least in part, the particular susceptibility of the SNpc to different types of stresses, such as increased iron concentration, enzymatic (monoamine oxidase) and non-enzymatic (auto-oxidation) DA metabolism, abnormal protein accumulation, proteasomal inhibition and neurotoxin-induced OS. In this context, the PA28 multisubunit proteasome activator (a component of the 26S proteasome) protein levels were shown to be very low in the SNpc of both normal and idiopatic PD subjects compared to other brain areas (McNaught et al., 2003), possibly exacerbating the already compromised DA-containing neurons to the various stress insults. The expression of SKP1A was not affected in the SNr or in cerebellum of PD.
A recent study has demonstrated that Parkin can also function within a novel SCF-like complex, along with the F-box/WD repeat protein hSel-10, responsible for substrate recognition, and Cullin1, which however does not include SKP1 (Staropoli et al., 2003). According to this finding, it is possible that Parkin associates with adaptor proteins other than hSel-10. Such complexes would likely display diverse substrate specificity and may explain the several targets that have been reported for Parkin (Dev et al., 2003).
Impairment in Proteasome Subunits
In addition to this selective decrement in SKP1A, a concomitant decrease in the 20S proteasome subunits alpha-5 (PSMA5), alpha-3 (PSMA3) and alpha-2 (PSMA2) mRNAs and in two subunits of the 19S regulatory complex of the 26S proteasome, the non-ATPase subunit 8 (PSMD8/Rpn12) and the ATPase subunit 4 (PSMC4/TBP7/Rpt3), was observed in the SN of PD brains, which may further contribute to dopaminergic neuronal damage in PD. The 19S complexes are located at one or both extremities of the 20S proteolytic core of the proteasome (Coux et al., 1996; Voges et al., 1999) and comprise at least 18 subunits (Glickman and Ciechanover, 2002). They are subdivided into two subcomplexes, the “base” and the “lid” which form the portions proximal and distal to the 20S proteolytic core respectively, and all together form the 26S complex. The lid consists of eight regulatory particle non-ATPase (Rpn) subunits, including PSMD8/Rpn12. One important function of the lid is to recognize multiubiquitinated proteins and other potential substrates of the proteasome (Glickman et al., 1998). The base contains three non-ATPases subunits (Rpn1, Rpn2 and Rpn10, this last common to both base and lid) and six putative ATPases subunits (Rpt1-6), including PSMC4/TBP7. These interact with non-native conformations of nonubiquitinated target proteins (Glickman et al., 1998; Voges et al., 1999; Strickland et al., 2000) and one of them S6′/TBP1, was found to bind aggregated or monomeric α-synuclein (Ghee et al., 2000; Snyder et al., 2003). In the same context, it has been recently reported that PSMC4, also called regulatory proteasomal protein S6, specifically interacts with both wild type and mutant synphilin-1 (Duke et al., 2004), another presynaptic protein, associated with synaptic vesicles (Ribeiro et al., 2002). Synphilin-1 was found associated with α-synuclein and, similarly, it accumulates in the Lewy body (Wakabayashi et al., 2000). Thus, the decreased levels of PSMC4 found in the present invention, together with the progressive inhibition of the UPS by the slow but persistent aggregation of α-synuclein, may cause their accumulation in Lewy bodies (Jellinger, 2003; Liani et al., 2004). The observed increase in the expression of the glycosaminoclycan heparan-sulphate gene in PD samples, may possible contribute as well to fibrillation and aggregation of α-synuclein, as suggested (Cohlberg et al., 2002). Each of the six ATPases is essential and they account for the ATP-requirement for proteolysis and for association of the 20S and 19S complexes to form the 26S complex (Ghislain et al., 1993; Gordon et al., 1993). Therefore, the decline in the expression of PSMC4/TBP7 in brains of Parkinsonian patients may contribute to decreased levels of 26S proteasome complex, abnormal accumulation of ubiquitinated proteins and reduced rates of degradation of short-lived proteins such as cyclins, which in turn may induce cell defects (for review see (Coux et al., 1996; Voges et al., 1999)). Indeed, we observed a decline in cyclin G associated kinase gene. Accumulation of cyclins and cyclin-dependent kinases, has been reported in post-mitotic neurons undergoing apoptosis (Padmanabhan et al., 1999; Copani et al., 2001) and in the MPTP model of PD (Grunblatt et al., 2001). It has been suggested that this event may represent an attempt at cell cycle re-entry (Verdaguer et al., 2002). Our findings are in line with a previous report of decreased expression of the 20S proteasome α-subunits, but not n-subunits, decreased protein expression levels of some 19S subunits and functional deficits in the 26/20S proteasome activity in the SNpc of patients with sporadic PD (McNaught et al., 2003). In this study however, the exact nature of the different subunits affected was not established, as the antibodies employed in the western blot analyses recognize sequences that are common to the various proteasome subunits and because of possible protein dimerization, leading to overestimation of the molecular weight. Thus, the present invention provides a wider view, revealing gene changes in the expression of specific proteasome components.
Heat Shock Protein Chaperone HSC-70
Proteotoxic insults to cells or several stress conditions can induce up regulation of molecular chaperones aimed at protecting cells by assisting to correct folding of wild type and mutated proteins. One of such is the 70 kDa heat-shock cognate protein (Hsc70) (Zinsmaier and Bronk, 2001), a member of the heat-shock protein 70 (Hsp70) family. It is suggested that Hsc70 normally mediates cellular processes such as protein folding, refolding, assembly, disassembly and protein translocation via biological membranes. It has been shown that overexpression of Hsp70 reduces the amount of misfolded, aggregated α-synuclein species in vivo and in vitro (Klucken et al., 2004), prevents the loss of dopaminergic neurons in a model of PD in transgenic flies that express either the wild-type α-synuclein, as well as the mutant forms A30P and A53T (Auluck et al., 2002) and suppresses the degeneration associated with polyglutamine pathogenicity (Bonini, 2002). Recently, a functional polymorphism in the 5′ promoter region of HSP70-1 has been reported in 274 PD patients, which may increase susceptibility to PD (Wu et al., 2004). We have found significant decreased gene expression of HSP8, coding for HSC70/HSC54, both in the SNpc and SNr of five out of six parkinsonian patients, as confirmed by real-time quantitative PCR, whereas no significant alterations were observed between cerebellum from PD patients and controls, suggestive of tissue specificity. HSPA8 exists in two isoforms, the truncated one lacking 153 amino acids residues in the protein binding and variable domain. (Tsukahara et al., 2000) It is thought to function as an endogenous inhibitory regulator of Hsc70 by competing the co-chaperones (Tsukahara et al., 2000). At present, we can not differentiate which of these are more avidly affected.
In addition to the involvement of Hsc70 in vesicular and nuclear trafficking they may play a role in protein ubiquitination through recognition of unfolded or aberrant proteins and delivery to a cochaperone, E3 ligase enzyme CHIP (carboxyl-terminus of Hsc70 interacting protein) (Murata et al., 2003). CHIP can cooperate with Hsp90 and/or Hsp70/Hsc70 and ubiquitinate their attached misfolded substrates. Thus the Hsc70-CHIP pair represents an E3 ligase for specifically recognizing unfolded proteins presented by the chaperones. This molecular chaperone-UPS quality control system is of extreme importance in neurodegenerative diseases, since proteins that evade refolding or degradation by the UPS form aggregates that accumulates into inclusion bodies. Indeed, molecular chaperones as well as ubiquitin and proteasome are recruited to inclusion bodies and Lewy Body (Stenoien et al., 1999; Sherman and Goldberg, 2001) indicating the attempt of the quality control system to degrade damaged proteins or prevent their removal. In fact, a number of neurodegenerative diseases appear to result from failure of the protein quality control system (Sherman and Goldberg, 2001).
Cell Adhesion Molecules, Iron and Oxidative Stress
An association between genes playing essential roles in neuronal development, such as cell migration and axonal elongation as well as components of the cytoskeleton and AD, has recently been suggested (De Ferrari and Inestrosa, 2000). A similar hypothesis can be assigned to PD where abnormal phosphorylation of cytoskeleton components (e.g. neurofilaments, microtubule-associated proteins) and ubiquitination of synaptic and other proteins, are encountered within the Lewy bodies (Jellinger, 2003).
Our present findings showing abnormal upregulation of genes with structural and phosphorylative molecular functions integral to membrane and extracellular matrix (cell adhesion process), support this assumption and extend our knowledge revealing potential new gene targets for future therapeutics.
During brain inflammation such as those reported in PD, AD and multiple sclerosis, the breakdown in the organization of key cell adhesion molecules and reduction in their signal transduction pathways may generate brain pathology, probably related to the recruitment of glia cells and macrophages and to the elevation in cytokines and OS. In addition, in five out of the six patients we observed a striking induction by more than 1.5 fold of EGLN1 (eg1 nine homolog 1) gene, a recently described proline hydroxylase enzyme belonging to the iron- and 2-oxoglutarate-dependent dioxygenase superfamily (Epstein et al., 2001). These enzymes act as key iron and oxygen sensor controlling the expression of the transcription factor hypoxia-inducible factor-1 alpha (HIF), a master regulator orchestrating the coordinated induction of an array of hypoxia-sensitive genes. The target genes of HIF are especially related to angiogenesis, cell proliferation/survival and glucose/iron metabolism (Lee et al., 2004). Upon high oxygen levels or iron overload, the EGLN hydroxylases targets HIF to proteasomal degradation.
Interestingly, the free iron-induced proteasomal-mediated degradation of iron regulatory protein (IRP2) involves also activation of 2-oxoglutarate-dependent dioxygenases and is inhibited by iron chelators (Hanson et al., 2003; Wang et al., 2004). Thus, it is possible that IRP2 is a substrate of EGLN1, which causes post-translational modification, signaling it for protein degradation. Excessive production of EGLN1 in the SNpc may lead to a fall in IRP2 and subsequent decrease in transferrin receptor (TfR) mRNA and increase in ferritin levels, both subjected to positive and negative transcriptional regulation by IRP2, respectively (Meyron-Holtz et al., 2004; Ponka, 2004) Recent studies in knock-out mice for IRP2, have revealed accumulation of iron in the striatum with substantial bradykinesia and tremor (LaVaute et al., 2001).
Increased expression of the iron and OS sensor protein may be directly responsible for the observed reduction in phosphofructokinase and also the angiogenic factor VEGF, both regulated by the HIF proteins (Minchenko et al., 2003). Phosphofructokinase is the key regulatory enzyme that controls the glucose flux through the glycolytic pathway. Similarly, VEGF activates genes involved in glucose transport and metabolism via activation of the PI3K and ras pathways. The expression of two major players in these pathways, ras homolog gene family, member B and phosphoinositide-3-kinase, catalytic, alpha polypeptide, were also down-regulated in parkinsonian SNpc. This finding support previous reports in human PD patients using positron emission tomography (PET) analysis, demonstrating a decrease in glucose uptake into the SN (Berding et al., 2001) and give a wider view of major survival pathways affected by the disease. These observations, together with the additional decrease in a number of energy pathways/glycolysis-related genes, as well as increases in iron/OS and inflammatory genes, is consistent with the hypothesis that mitochondrial dysfunction and reactive oxygen and nitrogen species contribute to the pathogenesis of PD. In line with this, a recent study linked a recessive mutation in the putative free radical sensor DJ-1 gene, with early-onset parkinsonism (Bonifati et al., 2003).
To our surprise, apoptosis-related gene constituted the smallest functional class in PD samples, questioning the relevance of program cell death in the neurodegenerative cascade of events occurring in the disease. The role of apoptosis in PD is highly controversial since the evidence for it does not correlate with the pathological findings and the rate of neurodegeneration. We cannot exclude the possibility that apoptotic changes occurred at earlier stages and they decline with disease progression.
Dopamine Neurotransmission and Metabolism
The disruption in these neuronal processes may be exacerbated by the observed reduction in the expression of cell signaling genes, which may be coupled to cell adhesion protein complexes and in genes belonging to the vesicle secretory pathway and dopaminergic neurotransmission and metabolism in the parkinsonian SN. In a highly stringent analysis, sought to detect crucial alterations in at least five out of the six PD patients, with 1.5 fold change in gene expression relative to the control group, we observed major decreases in the expression of 8 genes, 3 of them related to DA transmission, including the vesicular monoamine transporter VMAT2 (SCL18A2). The real-time PCR confirmatory analysis revealed a specific reduction of VMAT2 mRNA in the SNpc while the changes in the SNr were less pronounced. This finding is in agreement with previous reports on post mortem SN of controls and PD showing a marked reduction of VMAT2 mRNA in PD which was associated with a marked reduction in both dopamine transporter and VMAT2 signal per cell in the remaining pigmented neurons (Harrington et al., 1996; Brooks, 2003). These alterations were specific for the SN as no significant difference was seen in the cerebellum between the PD and control samples. Two other striking changes were related to the reduction in DA transmission and metabolism-related genes ALDH1A1 and ARPP-21 coding for aldehyde dehydrogenase (ALDH) and cAMP regulated phosphoprotein, respectively. ARPP-21 is specifically enriched in DA-innervated brain regions of the basal ganglia (e.g. caudate-putamen), and in the substantia nigra (Ouimet et al., 1989; Tsou et al., 1993). Similar to DARPP-32 (dopamine and adenosine 3′:5′-monophosphate-regulated phosphoprotein-32K), ARPP-21 is activated by DA receptor D1, and thus may represent an index of functional activity of D-1 neurotransmission. As far as we know, this is the first report of decreased ARPP-21 mRNA levels in PD SN. ALDH1A1 was found to be expressed highly and specifically in DA cells of the SN and ventral tegmental area (VTA) and to be markedly reduced in SNpc dopaminergic neurons but not in those of the VTA of PD brains (Gaiter et al., 2003). This is in line with observation that striatal ADH activity in 6-hydroxydopamine or electrical induced lesion in rats (Agid et al., 1973) or cats (Duncan et al., 1972), was significantly reduced.
ALDHs are involved in the degradation of DA to aldehyde derivatives (3,4-dihydroxyphenylacetaldehyde and 4-hydroxy-3-methoxyphenylacetaldehyde) (Mardh and Vallee, 1986), which are then metabolized to homovanilic acid and dihydroxyphenylacetic acid and in detoxification of aldehydes which are highly reactive and neurotoxic (Hjelle and Petersen, 1983). Thus, alteration in DA transmission may alter ALDH activity and/or, conversely, changes in ALDH-mediated metabolism may affect DA levels in nerve cell bodies and terminal fields in basal ganglia and the limbic system. These proteins, in conjunction with VMAT2 may now be considered new markers for PD.
The overall decrease in the levels of these DA neurotransmission and metabolism-related genes observed in the present invention, may seriously compromise neurotransmitter storage and correlate with DA neuron dysfunction. This may further be exacerbated by the reduced levels of SEC22L1 mRNA, a member of the SEC22 family of early vesicle trafficking proteins from the endoplasmic reticulum to (Hay et al., 1996), and by the observed reduction in many membrane carrier transporters. Impairment of cell transport and vesicle trafficking, carrying a wide repertoire of intracellular and membrane proteins, may well account for their aggregation and cell deposition into cytoplasmic inclusions (aggresomes) and in the Lewy body. One simple hypothesis explaining our observations is that loss of DA-containing neuronal bodies and neuronal synapses may account for the reduce levels of the gene products. This assumption would predict similar gene expression changes in a vast repertoire of vesicle traffic or synaptic-related genes. However, as emerged from our microarray and real-time PCR confirmations, the expression of many genes related to these categories such as Rab3b, syntaxin 6 and COP zeta 2, was not altered in PD.
Protein Modification/Phosphorylation
The expression of two pre-mRNA splicing related genes was inversely and highly affected in each of the five or six PD samples compared to control brains: the SRPK2 gene coding for the SFRS protein kinase 2, was down-regulated and SRRM2 coding for serine/arginine repetitive matrix 2, was significantly increased. SRPK2 is a kinase highly expressed in brain with a predicted sequence for specific phosphorylation of arginine/serine-rich (SR) domain-containing splicing factors, which in turn regulate RNA splicing in brain regions (Wang et al., 1998). Interestingly, SRRM2 (or SRm300), is a SR domain-containing splicing factor functioning as a co-activator within a major complex SRm 160/300, responsible for the processing of a subset of constitutively spliced pre-mRNAs (Blencowe et al., 2000).
If SRMM2 is a substrate of SRPK2, then it seems plausible that the robust decline in the kinase SRPK2 and the high expression of SRMM2 in PD SNpc, leads to abnormal hypophosphorylated SRMM2 accumulation. This in turn may affect the nuclear distribution of SRm160 and/or other associated SR proteins converging eventually in a broad impairment in splicing activity.
In spite of the vast information regarding the pathology and etiology of PD, it is still premature to assert what is the primary event(s) that triggers the development of PD. The present invention shows a reduction in the levels of selective components of the UPS together with a progressive misregulation of extracellular matrix/cytoskeleton components, concurrent with a state of OS and inflammation. These series of events may act independently or cooperatively during the course of the disease, leading eventually to the demise of dopaminergic neurons. Thus, subtle alterations in the kinetics of the different affected proteins may have, during the decades, a cumulative effect underlying the slowly progressive neurodegeneration of the DA-containing neurons in PD.
nidulans)
Minchenko O, Opentanova I and Caro J (2003) Hypoxic regulation of the 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase gene family (PFKFB-1-4) expression in vivo. FEBS Lett 554(3): 264-270.
Polymeropoulos M H, Lavedan C, Leroy E, Ide S E, Dehejia A, Dutra A, Pike B, Root H, Rubenstein J, Boyer R, Stenroos E S, Chandrasekharappa S, Athanassiadou A, Papapetropoulos T, Johnson W G, Lazzarini A M, Duvoisin R C, Di Iorio G, Golbe L I, Nussbaum R L. (1997) Mutation in the alpha-synuclein gene identified in families with Parkinson's disease. Science. 276(5321):2045-7.
Number | Date | Country | Kind |
---|---|---|---|
04000968.0 | Jan 2004 | EP | regional |
0401877.8 | Aug 2004 | EP | regional |
This is a continuation of application Ser. No. 10/586,523, filed Jun. 1, 2007, which is a national phase of international application PCT/IL05/00064, filed Jan. 19, 2005. The entire contents of the prior applications are hereby incorporated by reference.
Number | Date | Country | |
---|---|---|---|
Parent | 10586523 | Jun 2007 | US |
Child | 12716819 | US |