A Cross-Study Transcriptional Analysis of Parkinson's Disease

Abstract

The study of Parkinson's disease (PD), like other complex neurodegenerative disorders, is limited by access to brain tissue from patients with a confirmed diagnosis. Alternatively the study of peripheral tissues may offer some insight into the molecular basis of disease susceptibility and progression, but this approach still relies on brain tissue to benchmark relevant molecular changes against. Several studies have reported whole-genome expression profiling in post-mortem brain but reported concordance between these analyses is lacking. Here we apply a standardised pathway analysis to seven independent case-control studies, and demonstrate increased concordance between data sets. Moreover data convergence increased when the analysis was limited to the five substantia nigra (SN) data sets; this highlighted the down regulation of dopamine receptor signaling and insulin-like growth factor 1 (IGF1) signaling pathways. We also show that case-control comparisons of affected post mortem brain tissue are more likely to reflect terminal cytoarchitectural differences rather than primary pathogenic mechanisms. The implementation of a correction factor for dopaminergic neuronal loss predictably resulted in the loss of significance of the dopamine signaling pathway while axon guidance pathways increased in significance. Interestingly the IGF1 signaling pathway was also over-represented when data from non-SN areas, unaffected or only terminally affected in PD, were considered. Our findings suggest that there is greater concordance in PD whole-genome expression profiling when standardised pathway membership rather than ranked gene list is used for comparison.

About the Authors:

Greg T. Sutherland

Contributed equally to this work with: Greg T. Sutherland, Nicholas A. Matigian

Affiliation: National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia

Nicholas A. Matigian

Contributed equally to this work with: Greg T. Sutherland, Nicholas A. Matigian

Affiliation: National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia

Alistair M. Chalk

Affiliation: National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia

Matthew J. Anderson

Affiliation: National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia

Peter A. Silburn

Affiliations National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia, The School of Medicine, University of Queensland, Brisbane, Queensland, Australia

Alan Mackay-Sim

Affiliation: National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia

Christine A. Wells

Affiliation: National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia

George D. Mellick

* E-mail: [email protected]

Affiliations National Centre for Adult Stem Cell Research, Eskitis Institute for Cell and Molecular Therapies, Griffith University, Brisbane, Queensland, Australia, Department of Neurology, Princess Alexandra Hospital, Brisbane, Queensland, Australia, Department of Neurology, Royal Brisbane and Women's Hospital, Brisbane, Queensland, Australia

Introduction

Parkinson's disease (PD, OMIM: #168600) is a uniquely human disease that is clinically characterised by cardinal motor symptoms such as postural instability, bradykinesia and resting tremor [1]. In the PD brain there is pathognomonic loss of more than 50% of the dopaminergic neurons in the substantia nigra pars compacta (SNc) located in the midbrain [2], [3]. However, the disease is also characterised by non-motor symptoms such as sleep disorders, depression, hyposmia and autonomic dysfunction [4]–[7]. Accordingly, pathology in PD is not just restricted to the SNc but also affects the olfactory pathway, spinal cord and dorsal cranial nuclei of the medulla [8]–[11].

High throughput discovery platforms, such as microarrays, that assume no a priori aetiological hypotheses, promise much in elucidating the pathogenesis of complex diseases such as PD. Moreover, one would hope that these microarray data would reveal clues previously inaccessible via other means. Several microarray-based studies have used human tissue to look for differentially expressed genes in Parkinson's disease [12]–[21]. The majority of these used post-mortem whole brain tissue from the substantia nigra (SN) [14]–[17] and although most authors emphasised differential expression in the ubiquitin-proteasome system or cellular energy pathways, their published gene lists appeared quite discordant. Others extended their studies to include pathologically normal brain regions [12], [18], [20], [21] and these highlighted other biological mechanisms such as G-Protein-coupled receptor signaling and transcriptional regulation.

One study assayed SNc dopaminergic neurons only, following isolation by laser capture microscopy (LCM) [13]. Here gender differences were more pronounced than PD versus control differences. Uniquely, one study compared the transcriptomes of whole blood samples [19] and reported expression differences in a number of unrelated genes.

Given an apparent lack of concordance in published data sets one might ask what relevance these transcriptional approaches can have to PD pathogenesis. Certainly the utilisation of post mortem brain tissue appears to represent the best opportunity for finding PD-specific changes in gene expression. Furthermore such ‘benchmarks’ facilitate the evaluation of clinical samples and model systems for their utility in PD research. However the approaches to generating and analysing microarray data are not standardised and therefore could account for much of the apparent discrepancy between reported gene lists.

Here we apply a common analytical approach to the available human transcriptomic data in an attempt to find greater data convergence and generate new insight into the pathways systematically altered during PD pathogenesis. We have also generated an online search tool and extend an invitation to other researchers to explore the data themselves http://ncascr.griffith.edu.au/pdreview/2008/.

Results

Demographics of microarray studies using PD case-control tissue

Ten published transcriptomic studies met the initial criterion of comparing primary tissues derived from PD patients and controls (Table 1). Of these, seven were included because they used the same common gene expression array. These seven studies provided 13 case-control data sets comparing tissue from 119 PD cases and 74 controls (Table 2). Three studies did not meet the inclusion criteria: (1) Grunblatt et al. [15] who used a Focus array (Affymetrix) with limited gene coverage; (2) Miller et al [22] who used the Codelink™ bioarray platform and (3) Bossers et al. [12] who used Agilent technology. We were also unable to source raw data from the study of Papapetropoulos and colleagues [18]. Although the raw data from these latter experiments was not included in our combined re-analysis, their published gene lists were used when initially evaluating data convergence between all studies.

[Figure omitted. See PDF.]

Table 1. Summary of tissue used, size and array feature of published Parkinson's disease microarray publications.

https://doi.org/10.1371/journal.pone.0004955.t001

[Figure omitted. See PDF.]

Table 2. Summary of each data set used in the re-analysis and the number of differentially expressed probes before and after neuronal correction.

https://doi.org/10.1371/journal.pone.0004955.t002

Lack of convergence of pathways between studies from published ranked gene lists

We postulated, as others had done previously, that genes and pathways that appear consistently as differentially expressed in multiple studies and different source tissues are likely to be important in PD [20]. To look for such convergence we compared the published lists of differentially expressed genes using the Ingenuity Pathways Analysis (IPA) package. There was little overlap between published ranked gene lists and there were no genes consistently identified in more than three datasets (Table S1 ). We then asked whether we would see greater convergence using pathway analysis of the originally reported ranked gene lists. The ERK/MAPK was the most highly represented pathway, although it was over-represented in only four of the 11 data sets (Table 3).

[Figure omitted. See PDF.]

Table 3. Summary of over-represented IPA pathways from the published array data.

https://doi.org/10.1371/journal.pone.0004955.t003

Application of common data analysis methodology

Our common analysis method was applied to the 13 datasets (from seven studies) that met our platform inclusion criteria, and new ranked probe lists were produced. These are listed in Table S2, while the number of differentially-expressed probes is shown in Table 2. Additionally these data can be viewed online http://ncascr.griffith.edu.au/pdreview/2008/. This improved the convergence of the genes between the datasets with 20 genes now consistently differently regulated across six of 13 datasets (data not shown). Ranked IPA pathways from these gene lists can be found in Table S3.

Pathway analysis of SN data sets reveals a common dysregulation of dopamine signaling

Dopaminergic neuron loss in the SNc is the prominent neuropathological entity in PD so we initially focused on the SN data sets for their convergence and reproducibility. Over-representation of the dopamine receptor signaling pathway was consistently and significantly observed in all SN data sets (p-values >0.003–0.026) suggesting that not only can PD-related pathways be dissected out of complex transcriptomic data but that these changes are robustly reproducible between comparable studies (Table 4 and Table S2 A–E for full ranked pathway lists).

[Figure omitted. See PDF.]

Table 4. Summary of over-represented IPA pathways in the SN data sets before and after neuronal correction.

https://doi.org/10.1371/journal.pone.0004955.t004

Neuronal loss in PD is more severe in the lateral SN compared to the medial SN [23]–[25] while areas such as the superior frontal gyrus (SFG) are largely unaffected. A comparison of these three anatomical areas from the data in Moran et al. (LSN v MSN v SFG) showed a direct correlation between fold changes in genes of dopamine receptor signaling pathway and the severity of PD neuronal loss [14]. Therefore, the prominence of this pathway appeared to represent the disparate numbers of residual dopaminergic neurons in PD compared to control brains.

Correcting for dopaminergic neuronal loss changes ranked pathways

In order to bias our analysis towards underlying pathogenic mechanisms rather than terminal pathology, we devised a correction paradigm based on Moran's observations on neuronal loss. Our rationale is described in detail in Table S4. 217 genes with fold changes LSN>MSN>SFG from the Moran et al. data sets were defined as potentially “neuronal-loss-associated”, as these were the genes whose differential expression most likely resulted from relative loss of dopaminergic neurons from the brain regions sampled, rather than transcriptomic differences in residual cells. However eight of these genes could be shown to be differentially expressed in residual dopaminergic neurons (Cantuti-Castelvetri et al. study)[13] and so were retained in the analysis. The removal of 209 genes (Table S5) from the SN ranked lists resulted in two alternative pathways gaining prominence: ephrin receptor signaling (p-values <0.003–0.04) and the axonal guidance pathway (p-values <0.004–0.049) (Table 4 and Table S6 A–E for full ranked pathway lists).

Assessment of glial contribution to the neuronal-loss corrected datasets

Given the loss of dopaminergic cells in the PD SN the major contribution to the expression profile in the SN PD samples would presumably now come from the non-dopaminergic cells. Furthermore neuropathology in the PD SN is characterised by a reactive gliosis or “glial inflammation” (reviewed by Orr et al., 2002) [26]. Therefore glial markers and in particular reactive microglia markers such as CD68 and ICAM-1 might be expected to be upregulated in the SN data sets [26], [27]. However these genes and those of glial markers in general were largely indifferent between PD and controls. This potential anomaly is illustrated further using a selection of PD-related glial markers (Table S7).

Analysis of non-SN brain tissue highlights growth factor signaling pathways

We also analysed non-SN tissues as they are not subject to cytoarchitectural changes seen in the SN or are only affected late in the disease. Five IPA pathways were overrepresented in three out of five non-SN data sets (Table 5). Of these, IGF-1 and VEGF signaling were also found to be dysregulated in two of the neuronal loss-corrected SN datasets. Whereas, the four pathways (PTEN signaling, JAK/STAT signaling, ephrin receptor signaling, axonal guidance signaling) found to be significance in three or more of the corrected SN datasets were only significance in two of the non-SN datasets.

[Figure omitted. See PDF.]

Table 5. Comparison of over-represented IPA pathways in non-SN data sets and neuronal-loss corrected SN data sets.

https://doi.org/10.1371/journal.pone.0004955.t005

Differential expression in whole blood

Finally, given its clinical accessibility, we also re-analysed a whole blood dataset [19]. Our analysis revealed that the inositol phosphate metabolism and VEGF signaling pathways were the most differentially expressed in this dataset. The prominence of these pathways was quite distinct from SN tissues but commonality with non-SN tissues was observed with the VEGF signaling pathway (Table S2 L).

Discussion

Microarrays promise much in elucidating the pathogenesis of complex diseases such as PD but the lack of concordance in published data sets to date certainly questions their relevance. Here we have shown that a standardised approach to analysing PD-related microarray data can account for a considerable proportion of the discordance. We used a common analytical approach which improved data convergence and uncovered new leads for PD pathogenesis. We also recognised a potential anatomical bias in the datasets derived from brain regions with high neuronal loss. Our approach therefore provided an improved comparative analysis between existing datasets and further considered ‘tissue-of-origin’ effects.

Pathway analysis may uncover concordance between datasets not found in gene lists

Complex phenotypes, by their very nature, are aetiologically heterogeneous. This implies that single gene signatures may not be shared by all affected individuals. However, the identification of particularly relevant genetic pathways, have a higher probability of being revealed as convergent across multiple individuals and multiple studies than individual genes per se. Moreover, the differences reflected in a pathway or network of genes may be robust enough to overcome the effects of experimental noise and inter-study variability prone to bias single gene expression values. Therefore, this approach results in an increase in sensitivity to detect interesting and novel patterns in gene expression between multiple samples of cases and controls.

Common data analysis highlights dopamine signaling pathway in SN

Context is very important in gene expression studies and as expected the analysis of SN tissue-derived data sets further improved our concordance and highlighted the ‘dopamine receptor signaling’ pathway. However rather than representing a primary pathogenic effect the extensive down regulation of genes such as DOPA decarboxylase (DDC), dopamine receptor 2 (DRD2), dopamine transporter (SLC6A3 or DAT) and tyrosine hydroxylase (TH) was probably entirely due to a disproportionate number of SN dopaminergic neurons between cases and controls. On a positive note, the microarray data was providing accurate molecular fingerprints of the comparative tissue being examined; a modern pixelated analogy to the histological section. However the relative neuronal loss in the SN also reduced our signal to noise ratio for pathogenic relevance.

Accordingly we pursued two alternative approaches to maximise potentially useful information on the underlying biological processes. First, the implementation of a correction factor for dopaminergic neuronal loss in the SN data sets and second, the analysis of non-SN or unaffected tissues for data convergence. The purpose of the correction factor was not to magically recreate the early disease landscape but to remove ‘red herrings’ that solely reflected the relative numbers of neurons between cases and controls. The retention of genes differentially regulated in the residual dopaminergic neuron data set should have improved the overall specificity of this approach.

Following our ‘neuronal loss’ correction, two alternative pathways gained prominence: ephrin receptor signaling and the axonal guidance pathway. The latter is consistent with the findings of recent studies [12], [28] including one that combined gene expression data with genotype data from two genome-wide association studies [28]. There is actually considerable overlap between the axon guidance and ephrin signaling pathways with ephrins along with netrins, slits and semaphorins being the main families of guidance molecules in the developing nervous system [29]. It remains to be clarified whether the differential expression of axon guidance genes in PD represent neurodevelopmental manifestations, compensatory attempts at rewiring, or dysregulated expression patterns induced by a devastated environment.

As discussed above microarray data is very powerful in illustrating cytoarchitectural differences between cases and controls such as dopaminergic neuron loss. Given the considerable literature supporting the involvement of glia in PD pathogenesis [26], [27], [30], [31] we might have expected glial markers to inversely differentiated in the SN data sets. The absence of significant fold changes in activated microglial markers in particular argues against a substantial glial component to neurodegeneration in the terminal PD brain.

Growth factor signaling pathways are prominent in non-SN tissues

The distinct gene expression pattern of brain areas that are not overtly affected by PD pathology may be less confounded than the SN with respect to the cell death associated with PD. It could be argued that the transcriptomes of unaffected tissues might be too divergent from those of predilection sites such as the SN, such that they provide very little informative data. Our analysis, which has highlighted consistent differences in growth-factor signaling in non-SN datasets, argues that areas affected late in the disease, such as the prefrontal cortex [8] could at the time of post mortem exhibit similar mechanisms of degeneration as initially occurred in the SN. Our approach highlights the differential expression of IGF-1 and VEGF signaling pathways. Importantly the IGF1 signaling pathway was also over-represented in two ‘corrected’ SN datasets (Table 5).

This pathway has been largely unexamined for associations with PD although IGF1 signaling is reported to have neuroprotective effects on dopaminergic neurons [32], [33] and it has recently been suggested that excessive IGF1 signaling accelerates ageing through deleteriously effects on protective mechanisms against proteotoxicity such as Lewy body formation in PD [34]. Furthermore we have recently showed that a polymorphism in the 3′ untranslated region of the IGF2 gene, a homologue of IGF1 was protective against PD [35].

Similarly VEGF is known to promote the growth and survival of dopaminergic neurons [36]–[38] Interestingly both IGF1 and insulin enhance VEGF expression in vitro [39], [40] providing a plausible mechanism that might underlie the co-prominence of these signaling pathways in our re-analysis.

Challenges and Future Directions

Case-control expression analysis in a degenerative disease like PD poses difficult issues when attempting to uncover pathways contributing to disease initiation. It would be advantageous to target tissues that express the proteins that are fundamental to the disease process and are different in individuals who are at risk of the disease. At the same time we need to account for any influences of the pathological process on these profiles. Microarray data of predilection sites such as the SN illustrates cytoarchitectural differences between cases and controls but to understand some of the early pathogenic processes, we would ideally want to assay a brain region very similar to SN but that is only belatedly affected.

An additional consideration is the ability of the pathway approach, used in our analyses, to provide adequate specificity for PD over other neurodegenerative conditions. This issue remains to be clarified, and requires further investigation. It is important to recognise that there may be genetic expression patterns common to neurodegenerative diseases, generally. These may reflect common pathological changes (such as cell death, markers of oxidative stress or neuro-inflammation etc) or shared risk factors influencing neurodegeneration.

There are still inherent difficulties in obtaining reproducible gene expression data from post mortem brain, even if an optimal region of the brain could be assayed [41]–[44]. Furthermore this information can only be used retrospectively for the potential benefit of future PD patients. Therefore there is considerable interest in developing strategies to obtain human RNA from more accessible sources such as blood or neuronal-like cell lines. The finding of down regulation of α-synuclein in microarray analysis of whole blood samples from PD patients versus controls is exciting because it implies that PD-specific changes can be found ante mortem in readily accessible tissues [19]. It is also of note that these tissues also revealed differential expression of the VEGF signaling pathway (observed in other non-SN tissue samples). However, there is a real risk that peripheral tissues, such as whole blood examined here, may express few proteins fundamental to the disease process and therefore be of limited ability to demonstrate case - control differences relevant to nervous system disorders [45]. The lack of available gene expression data from multiple tissues in PD patients at various stages of the disease prevents such an analysis but highlights the need for ongoing research efforts in this area.

Interestingly one peripherally accessible neural tissue, the olfactory mucosa, has been used to demonstrate significant differences in functional assays and gene expression between schizophrenics, bipolar affective disorder and controls [46]. Such cells from PD patients and controls may yet provide an opportunity to interrogate neuronal mechanisms without relying on post mortem tissue [47].

In this paper we have presented a summary of the available microarray data from PD case-control studies and have suggested some potential strategies for uncovering primary pathogenic mechanisms. For others who wish to use and explore these data we have constructed an online database which enables rapid evaluation on a single gene or pathway basis.

Materials and Methods

We conducted literature searches in National Center for Biotechnology Information (NCBI) PubMed and dataset searches in NCBI gene expression omnibus and ArrayExpress (EBI) [48], [49] to identify all reported microarray studies that explored differential gene expression in Parkinson disease. Studies satisfied the inclusion criteria if they: 1) compared tissues from PD patients and controls; 2) assessed transcripts on a genome-wide basis; and 3) used Affymetrix gene expression arrays. For the selected studies we obtained raw microarray data (CEL files) from public microarrays repositories [13], [14], [19] or from the study authors [16, Vogt, 2006 #1529,21,]. In one case the data was publicly available from the follow up study [28] rather than the primary study [18].

Meta-analysis data summarisation, normalisation and analysis of variance

All studies used Affymetrix arrays, the probes on the arrays and the experimentally chosen fluorescence thresholds varied. Consequently, the data could not be simply combined without avoiding study bias and the effects of probe-level sequence information [50]. To overcome this problem, the raw data (as CEL files) for each dataset were imported individually into GeneSpring® ™ 7.3.1(Agilent) and the probes sets were summarised and normalised by the Robust Multichip Average (RMA) algorithm[51]. Some studies included non-PD disease controls but our analysis was performed on PD and control patients only (Table 1). For studies which used multiple brain regions, each area was treated as a separate data set. Differentially expressed genes between PD and controls were determined by an analysis of variances (ANOVA) using a Welch t-test with a p-value cut-off of ≤0.01.

Pathway over-representation analysis

The ranked genes lists for each study were assessed by integrating the data at a pathway level. Each ranked list was imported into the Ingenuity Pathways Analysis 6.3 (IPA, from Ingenuity® Systems, www.ingenuity.com) which incorporates an extensive literature-derived knowledge base from which to assign pathway affiliation. The significance value for pathway over-representation was calculated using a right-tailed Fisher's exact test. Each pathway was ranked by assessing the number of studies that were statistically over-represented (p-value ≤0.05). This pathway over-representation ranking was performed individually on data sets utilising substantia nigra (SN) (Hauser, Moran LSN+MSN, Zhang, Lesnick)[14], [16], [21], [28] and non-SN tissues (Zhang-Putamen, Zhang-BA9, Moran-SFG, Vogt-OCT, Vogt-Putamen, Vogt-CB)[14], [20], [21]. Additionally this analysis was performed independently on data sets derived from whole blood (Scherzer)[19] and dopaminergic neurons (Cantuti-Castelvetri)[13].

Correction for dopaminergic neuronal loss

The substantia nigra pars compacta of PD patients is characterised by the loss of neuromelanin-containing dopaminergic neurons [2], [3]. Furthermore neurons are lost in a particular pattern; severity decreasing from ventrolateral to dorsomedial [23]–[25]. In an attempt to correct for the effects on expression arising directly from the neuronal loss associated with PD, we devised the following correction paradigm. We first used data from the Moran et al. study[14] which compared gene expression profiles in three brain regions, lateral SN (LSN), medial SN (MSN) and superior frontal gyrus (SFG), from the same patients. The actual neuronal loss in PD is known to be greater in LSN compared to MSN, with the SFG relatively spared. Therefore significant probes (p≤0.01) with fold changes in LSN>MSN>SFG were defined as potentially “neuronal-loss-associated”. For example the tyrosine hydroxylase (TH) gene showed a fold change pattern of −14.5 (LSN), −4.9 (MSN) and no change (SFG). Any probes which were defined as “neuronal-loss-associated” and were not differentially expressed in residual laser-captured dopaminergic neurons (Cantuti-Castelvetri et al. study)[13] were removed from the subsequent analyses of all SN data sets.

Online Database

The differentially expressed gene list generated for each study by this re-analysis can be found with their respective p-values and fold changes can be found at http://ncascr.griffith.edu.au/pdreview/2008/. A search can be performed individually using Entrez gene ID, gene symbol, or collectively by publicly available lists/pathways.

Accession Numbers

The National Center for Biotechnology Information Entrez Gene website (http://www.ncbi.nlm.nih.gov/sites/entrez?dbgene) accession numbers (GeneIDs) for the genes named in the paper include: DDC(1644), DRD2(1813), SLC6A3(6531), TH(7054), IGF2(3481).

Supporting Information

[Figure omitted. See PDF.]

Table S1.

Comparison of overlap in genes between PD-related transcriptomic studies. The enclosed table illustrates the increase in data convergence between PD-related transcriptomic studies following the implementation of our common analysis methodology.

https://doi.org/10.1371/journal.pone.0004955.s001

(0.02 MB PDF)

Table S2.

The differentially expressed probes generated by common analysis for each study. Ranked probe lists for each study generated by common analysis method with fold change and p-value.

https://doi.org/10.1371/journal.pone.0004955.s002

(1.08 MB XLS)

Table S3.

Over-represented pathway categories from IPA analysis. Over-represented pathway categories from IPA analysis of the differentally expressed probes of PD patients compared to controls: A- SN - Hauser study; B- SN - Zhang study; C- lateral SN - Moran study; D- medial SN - Moran study; E- SN - Lesnick study; F- Brodmann Area 9 - Zhang study; G- Putamen - Zhang study; H- Superior Frontal Gyrus - Moran study; I- Occipital Cortex - Vogt study; J- Putamen - Vogt study; K- Cerebellum - Vogt study; L- Whole Blood - Scherzer study; M- Laser Captured SN dopaminergic neurons - Castelvetri study.

https://doi.org/10.1371/journal.pone.0004955.s003

(0.07 MB XLS)

Table S4.

Selection of neuronal loss-associated genes. A work flow diagram and hypothetical examples illustrate the selection of neuronal loss-associated genes that were removed from the SN datasets prior to pathway analysis.

https://doi.org/10.1371/journal.pone.0004955.s004

(0.12 MB PDF)

Table S5.

‘Neuronal-loss’ associated genes removed by correction paradigm. Identified probes with fold change's for the three brain regions in the Moran study that followed the pattern, LSN>MSN>SFG.

https://doi.org/10.1371/journal.pone.0004955.s005

(0.05 MB XLS)

Table S6.

Identified probes with fold change's for the three brain regions in the Moran study that followed the pattern, LSN>MSN>SFG. Over-represented pathway categories from IPA analysis after neuronal correction of the differentally expressed probes from SN tissue: A- Hauser study; B- Zhang study; C- lateral SN; Moran study; D- medial SN Moran study; E- Lesnick study.

https://doi.org/10.1371/journal.pone.0004955.s006

(0.04 MB XLS)

Table S7.

Fold changes in PD-related glial markers. The lack of differential expression of PD-related glial markers is illustrated in the enclosed table.

https://doi.org/10.1371/journal.pone.0004955.s007

(0.07 MB PDF)

Acknowledgments

The authors would like to acknowledge our PD colleagues who have kindly made their gene expression data available through public databases or forwarding it directly to us. We would like to thank Amanda Miotto and Othmar Korn for informatics support.

Author Contributions

Conceived and designed the experiments: GS NAM AC CAW GM. Performed the experiments: GS NAM. Analyzed the data: GS NAM AC MA AMS CAW GM. Contributed reagents/materials/analysis tools: GS NAM AC MA PS AMS CAW GM. Wrote the paper: GS NAM AC PS AMS CAW GM.

Citation: Sutherland GT, Matigian NA, Chalk AM, Anderson MJ, Silburn PA, Mackay-Sim A, et al. (2009) A Cross-Study Transcriptional Analysis of Parkinson's Disease. PLoS ONE4(3): e4955. https://doi.org/10.1371/journal.pone.0004955

References

1. Gelb DJ, Oliver E, Gilman S (1999) Diagnostic criteria for Parkinson disease. Arch Neurol 56: 33–39.DJ GelbE. OliverS. Gilman1999Diagnostic criteria for Parkinson disease.Arch Neurol563339

2. Hughes AJ (1997) Clinicopathological aspects of Parkinson's disease. Eur Neurol 38: Suppl 213–20.AJ Hughes1997Clinicopathological aspects of Parkinson's disease.Eur Neurol38Suppl 21320

3. Hughes AJ, Daniel SE, Kilford L, Lees AJ (1992) Accuracy of clinical diagnosis of idiopathic Parkinson's disease: a clinico-pathological study of 100 cases. J Neurol Neurosurg Psychiatry 55: 181–184.AJ HughesSE DanielL. KilfordAJ Lees1992Accuracy of clinical diagnosis of idiopathic Parkinson's disease: a clinico-pathological study of 100 cases.J Neurol Neurosurg Psychiatry55181184

4. Comella CL (2007) Sleep disorders in Parkinson's disease: an overview. Mov Disord 22: Suppl 17S367–373.CL Comella2007Sleep disorders in Parkinson's disease: an overview.Mov Disord22Suppl 17S367373

5. Dubow JS (2007) Autonomic dysfunction in Parkinson's disease. Dis Mon 53: 265–274.JS Dubow2007Autonomic dysfunction in Parkinson's disease.Dis Mon53265274

6. Hawkes C (2006) Olfaction in neurodegenerative disorder. Adv Otorhinolaryngol 63: 133–151.C. Hawkes2006Olfaction in neurodegenerative disorder.Adv Otorhinolaryngol63133151

7. Ziemssen T, Reichmann H (2007) Non-motor dysfunction in Parkinson's disease. Parkinsonism Relat Disord 13: 323–332.T. ZiemssenH. Reichmann2007Non-motor dysfunction in Parkinson's disease.Parkinsonism Relat Disord13323332

8. Braak H, Del Tredici K, Rub U, de Vos RA, Jansen Steur EN, et al. (2003) Staging of brain pathology related to sporadic Parkinson's disease. Neurobiol Aging 24: 197–211.H. BraakK. Del TrediciU. RubRA de VosEN Jansen Steur2003Staging of brain pathology related to sporadic Parkinson's disease.Neurobiol Aging24197211

9. Braak H, Muller CM, Rub U, Ackermann H, Bratzke H, et al. (2006) Pathology associated with sporadic Parkinson's disease–where does it end? J Neural Transm Suppl 89–97.H. BraakCM MullerU. RubH. AckermannH. Bratzke2006Pathology associated with sporadic Parkinson's disease–where does it end?J Neural Transm Suppl8997

10. Hubbard PS, Esiri MM, Reading M, McShane R, Nagy Z (2007) Alpha-synuclein pathology in the olfactory pathways of dementia patients. J Anat 211: 117–124.PS HubbardMM EsiriM. ReadingR. McShaneZ. Nagy2007Alpha-synuclein pathology in the olfactory pathways of dementia patients.J Anat211117124

11. Kalaitzakis ME, Graeber MB, Gentleman SM, Pearce RK (2008) Controversies over the staging of alpha-synuclein pathology in Parkinson's disease. Acta Neuropathol 116: 125–128.ME KalaitzakisMB GraeberSM GentlemanRK Pearce2008Controversies over the staging of alpha-synuclein pathology in Parkinson's disease.Acta Neuropathol116125128

12. Bossers K, Meerhoff G, Balesar R, van Dongen JW, Kruse CG, et al. (2008) Analysis of Gene Expression in Parkinson's Disease: Possible Involvement of Neurotrophic Support and Axon Guidance in Dopaminergic Cell Death. Brain Pathol. K. BossersG. MeerhoffR. BalesarJW van DongenCG Kruse2008Analysis of Gene Expression in Parkinson's Disease: Possible Involvement of Neurotrophic Support and Axon Guidance in Dopaminergic Cell Death.Brain Pathol

13. Cantuti-Castelvetri I, Keller-McGandy C, Bouzou B, Asteris G, Clark TW, et al. (2007) Effects of gender on nigral gene expression and parkinson disease. Neurobiol Dis 26: 606–614.I. Cantuti-CastelvetriC. Keller-McGandyB. BouzouG. AsterisTW Clark2007Effects of gender on nigral gene expression and parkinson disease.Neurobiol Dis26606614

14. Moran LB, Duke DC, Deprez M, Dexter DT, Pearce RK, et al. (2006) Whole genome expression profiling of the medial and lateral substantia nigra in Parkinson's disease. Neurogenetics 7: 1–11.LB MoranDC DukeM. DeprezDT DexterRK Pearce2006Whole genome expression profiling of the medial and lateral substantia nigra in Parkinson's disease.Neurogenetics7111

15. Grunblatt E, Mandel S, Jacob-Hirsch J, Zeligson S, Amariglo N, et al. (2004) Gene expression profiling of parkinsonian substantia nigra pars compacta; alterations in ubiquitin-proteasome, heat shock protein, iron and oxidative stress regulated proteins, cell adhesion/cellular matrix and vesicle trafficking genes. J Neural Transm 111: 1543–1573.E. GrunblattS. MandelJ. Jacob-HirschS. ZeligsonN. Amariglo2004Gene expression profiling of parkinsonian substantia nigra pars compacta; alterations in ubiquitin-proteasome, heat shock protein, iron and oxidative stress regulated proteins, cell adhesion/cellular matrix and vesicle trafficking genes.J Neural Transm11115431573

16. Hauser MA, Li YJ, Xu H, Noureddine MA, Shao YS, et al. (2005) Expression profiling of substantia nigra in Parkinson disease, progressive supranuclear palsy, and frontotemporal dementia with parkinsonism. Arch Neurol 62: 917–921.MA HauserYJ LiH. XuMA NoureddineYS Shao2005Expression profiling of substantia nigra in Parkinson disease, progressive supranuclear palsy, and frontotemporal dementia with parkinsonism.Arch Neurol62917921

17. Miller RM, Kiser GL, Kaysser-Kranich TM, Lockner RJ, Palaniappan C, et al. (2006) Robust dysregulation of gene expression in substantia nigra and striatum in Parkinson's disease. Neurobiol Dis 21: 305–313.RM MillerGL KiserTM Kaysser-KranichRJ LocknerC. Palaniappan2006Robust dysregulation of gene expression in substantia nigra and striatum in Parkinson's disease.Neurobiol Dis21305313

18. Papapetropoulos S, Ffrench-Mullen J, McCorquodale D, Qin Y, Pablo J, et al. (2006) Multiregional gene expression profiling identifies MRPS6 as a possible candidate gene for Parkinson's disease. Gene Expr 13: 205–215.S. PapapetropoulosJ. Ffrench-MullenD. McCorquodaleY. QinJ. Pablo2006Multiregional gene expression profiling identifies MRPS6 as a possible candidate gene for Parkinson's disease.Gene Expr13205215

19. Scherzer CR, Eklund AC, Morse LJ, Liao Z, Locascio JJ, et al. (2007) Molecular markers of early Parkinson's disease based on gene expression in blood. Proc Natl Acad Sci U S A 104: 955–960.CR ScherzerAC EklundLJ MorseZ. LiaoJJ Locascio2007Molecular markers of early Parkinson's disease based on gene expression in blood.Proc Natl Acad Sci U S A104955960

20. Vogt IR, Lees AJ, Evert BO, Klockgether T, Bonin M, et al. (2006) Transcriptional changes in multiple system atrophy and Parkinson's disease putamen. Exp Neurol 199: 465–478.IR VogtAJ LeesBO EvertT. KlockgetherM. Bonin2006Transcriptional changes in multiple system atrophy and Parkinson's disease putamen.Exp Neurol199465478

21. Zhang Y, James M, Middleton FA, Davis RL (2005) Transcriptional analysis of multiple brain regions in Parkinson's disease supports the involvement of specific protein processing, energy metabolism, and signaling pathways, and suggests novel disease mechanisms. Am J Med Genet B Neuropsychiatr Genet 137: 5–16.Y. ZhangM. JamesFA MiddletonRL Davis2005Transcriptional analysis of multiple brain regions in Parkinson's disease supports the involvement of specific protein processing, energy metabolism, and signaling pathways, and suggests novel disease mechanisms.Am J Med Genet B Neuropsychiatr Genet137516

22. Miller RM, Federoff HJ (2005) Altered gene expression profiles reveal similarities and differences between Parkinson disease and model systems. Neuroscientist 11: 539–549.RM MillerHJ Federoff2005Altered gene expression profiles reveal similarities and differences between Parkinson disease and model systems.Neuroscientist11539549

23. Damier P, Hirsch EC, Agid Y, Graybiel AM (1999) The substantia nigra of the human brain. II. Patterns of loss of dopamine-containing neurons in Parkinson's disease. Brain 122(Pt 8): 1437–1448.P. DamierEC HirschY. AgidAM Graybiel1999The substantia nigra of the human brain. II. Patterns of loss of dopamine-containing neurons in Parkinson's disease.Brain122(Pt 8)14371448

24. Fearnley JM, Lees AJ (1991) Ageing and Parkinson's disease: substantia nigra regional selectivity. Brain 114(Pt 5): 2283–2301.JM FearnleyAJ Lees1991Ageing and Parkinson's disease: substantia nigra regional selectivity.Brain114(Pt 5)22832301

25. Gibb WR (1991) Neuropathology of the substantia nigra. Eur Neurol 31: Suppl 148–59.WR Gibb1991Neuropathology of the substantia nigra.Eur Neurol31Suppl 14859

26. Orr CF, Rowe DB, Halliday GM (2002) An inflammatory review of Parkinson's disease. Prog Neurobiol 68: 325–340.CF OrrDB RoweGM Halliday2002An inflammatory review of Parkinson's disease.Prog Neurobiol68325340

27. McGeer PL, McGeer EG (2008) Glial reactions in Parkinson's disease. Mov Disord 23: 474–483.PL McGeerEG McGeer2008Glial reactions in Parkinson's disease.Mov Disord23474483

28. Lesnick TG, Papapetropoulos S, Mash DC, Ffrench-Mullen J, Shehadeh L, et al. (2007) A genomic pathway approach to a complex disease: axon guidance and Parkinson disease. PLoS Genet 3: e98.TG LesnickS. PapapetropoulosDC MashJ. Ffrench-MullenL. Shehadeh2007A genomic pathway approach to a complex disease: axon guidance and Parkinson disease.PLoS Genet3e98

29. Dickson BJ (2002) Molecular mechanisms of axon guidance. Science 298: 1959–1964.BJ Dickson2002Molecular mechanisms of axon guidance.Science29819591964

30. Croisier E, Graeber MB (2006) Glial degeneration and reactive gliosis in alpha-synucleinopathies: the emerging concept of primary gliodegeneration. Acta Neuropathol 112: 517–530.E. CroisierMB Graeber2006Glial degeneration and reactive gliosis in alpha-synucleinopathies: the emerging concept of primary gliodegeneration.Acta Neuropathol112517530

31. Mirza B, Hadberg H, Thomsen P, Moos T (2000) The absence of reactive astrocytosis is indicative of a unique inflammatory process in Parkinson's disease. Neuroscience 95: 425–432.B. MirzaH. HadbergP. ThomsenT. Moos2000The absence of reactive astrocytosis is indicative of a unique inflammatory process in Parkinson's disease.Neuroscience95425432

32. Chung CY, Seo H, Sonntag KC, Brooks A, Lin L, et al. (2005) Cell type-specific gene expression of midbrain dopaminergic neurons reveals molecules involved in their vulnerability and protection. Hum Mol Genet 14: 1709–1725.CY ChungH. SeoKC SonntagA. BrooksL. Lin2005Cell type-specific gene expression of midbrain dopaminergic neurons reveals molecules involved in their vulnerability and protection.Hum Mol Genet1417091725

33. Offen D, Shtaif B, Hadad D, Weizman A, Melamed E, et al. (2001) Protective effect of insulin-like-growth-factor-1 against dopamine-induced neurotoxicity in human and rodent neuronal cultures: possible implications for Parkinson's disease. Neurosci Lett 316: 129–132.D. OffenB. ShtaifD. HadadA. WeizmanE. Melamed2001Protective effect of insulin-like-growth-factor-1 against dopamine-induced neurotoxicity in human and rodent neuronal cultures: possible implications for Parkinson's disease.Neurosci Lett316129132

34. Cohen E, Dillin A (2008) The insulin paradox: aging, proteotoxicity and neurodegeneration. Nat Rev Neurosci 9: 759–767.E. CohenA. Dillin2008The insulin paradox: aging, proteotoxicity and neurodegeneration.Nat Rev Neurosci9759767

35. Sutherland G, Mellick G, Newman J, Double KL, Stevens J, et al. (2008) Haplotype analysis of the IGF2-INS-TH gene cluster in Parkinson's disease. Am J Med Genet B Neuropsychiatr Genet 147B: 495–499.G. SutherlandG. MellickJ. NewmanKL DoubleJ. Stevens2008Haplotype analysis of the IGF2-INS-TH gene cluster in Parkinson's disease.Am J Med Genet B Neuropsychiatr Genet147B495499

36. Barcia C, Bautista V, Sanchez-Bahillo A, Fernandez-Villalba E, Faucheux B, et al. (2005) Changes in vascularization in substantia nigra pars compacta of monkeys rendered parkinsonian. J Neural Transm 112: 1237–1248.C. BarciaV. BautistaA. Sanchez-BahilloE. Fernandez-VillalbaB. Faucheux2005Changes in vascularization in substantia nigra pars compacta of monkeys rendered parkinsonian.J Neural Transm11212371248

37. Silverman WF, Krum JM, Mani N, Rosenstein JM (1999) Vascular, glial and neuronal effects of vascular endothelial growth factor in mesencephalic explant cultures. Neuroscience 90: 1529–1541.WF SilvermanJM KrumN. ManiJM Rosenstein1999Vascular, glial and neuronal effects of vascular endothelial growth factor in mesencephalic explant cultures.Neuroscience9015291541

38. Yasuhara T, Shingo T, Muraoka K, Kameda M, Agari T, et al. (2005) Neurorescue effects of VEGF on a rat model of Parkinson's disease. Brain Res 1053: 10–18.T. YasuharaT. ShingoK. MuraokaM. KamedaT. Agari2005Neurorescue effects of VEGF on a rat model of Parkinson's disease.Brain Res10531018

39. Miele C, Rochford JJ, Filippa N, Giorgetti-Peraldi S, Van Obberghen E (2000) Insulin and insulin-like growth factor-I induce vascular endothelial growth factor mRNA expression via different signaling pathways. J Biol Chem 275: 21695–21702.C. MieleJJ RochfordN. FilippaS. Giorgetti-PeraldiE. Van Obberghen2000Insulin and insulin-like growth factor-I induce vascular endothelial growth factor mRNA expression via different signaling pathways.J Biol Chem2752169521702

40. Beckert S, Farrahi F, Perveen Ghani Q, Aslam R, Scheuenstuhl H, et al. (2006) IGF-I-induced VEGF expression in HUVEC involves phosphorylation and inhibition of poly (ADP-ribose) polymerase. Biochem Biophys Res Commun 341: 67–72.S. BeckertF. FarrahiQ. Perveen GhaniR. AslamH. Scheuenstuhl2006IGF-I-induced VEGF expression in HUVEC involves phosphorylation and inhibition of poly (ADP-ribose) polymerase.Biochem Biophys Res Commun3416772

41. Preece P, Cairns NJ (2003) Quantifying mRNA in postmortem human brain: influence of gender, age at death, postmortem interval, brain pH, agonal state and inter-lobe mRNA variance. Brain Res Mol Brain Res 118: 60–71.P. PreeceNJ Cairns2003Quantifying mRNA in postmortem human brain: influence of gender, age at death, postmortem interval, brain pH, agonal state and inter-lobe mRNA variance.Brain Res Mol Brain Res1186071

42. Preece P, Virley DJ, Costandi M, Coombes R, Moss SJ, et al. (2003) An optimistic view for quantifying mRNA in post-mortem human brain. Brain Res Mol Brain Res 116: 7–16.P. PreeceDJ VirleyM. CostandiR. CoombesSJ Moss2003An optimistic view for quantifying mRNA in post-mortem human brain.Brain Res Mol Brain Res116716

43. Marcotte ER, Srivastava LK, Quirion R (2003) cDNA microarray and proteomic approaches in the study of brain diseases: focus on schizophrenia and Alzheimer's disease. Pharmacol Ther 100: 63–74.ER MarcotteLK SrivastavaR. Quirion2003cDNA microarray and proteomic approaches in the study of brain diseases: focus on schizophrenia and Alzheimer's disease.Pharmacol Ther1006374

44. Atz M, Walsh D, Cartagena P, Li J, Evans S, et al. (2007) Methodological considerations for gene expression profiling of human brain. J Neurosci Methods 163: 295–309.M. AtzD. WalshP. CartagenaJ. LiS. Evans2007Methodological considerations for gene expression profiling of human brain.J Neurosci Methods163295309

45. Matigian NA, McCurdy RD, Feron F, Perry C, Smith H, et al. (2008) Fibroblast and lymphoblast gene expression profiles in schizophrenia: are non-neural cells informative? PLoS ONE 3: e2412.NA MatigianRD McCurdyF. FeronC. PerryH. Smith2008Fibroblast and lymphoblast gene expression profiles in schizophrenia: are non-neural cells informative?PLoS ONE3e2412

46. McCurdy RD, Feron F, Perry C, Chant DC, McLean D, et al. (2006) Cell cycle alterations in biopsied olfactory neuroepithelium in schizophrenia and bipolar I disorder using cell culture and gene expression analyses. Schizophr Res 82: 163–173.RD McCurdyF. FeronC. PerryDC ChantD. McLean2006Cell cycle alterations in biopsied olfactory neuroepithelium in schizophrenia and bipolar I disorder using cell culture and gene expression analyses.Schizophr Res82163173

47. Murrell W, Wetzig A, Donnellan M, Feron F, Burne T, et al. (2008) Olfactory Mucosa is a Potential Source for Autologous Stem Cell Therapy for Parkinson's Disease. Stem Cells. W. MurrellA. WetzigM. DonnellanF. FeronT. Burne2008Olfactory Mucosa is a Potential Source for Autologous Stem Cell Therapy for Parkinson's Disease.Stem Cells

48. Edgar R, Domrachev M, Lash AE (2002) Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 30: 207–210.R. EdgarM. DomrachevAE Lash2002Gene Expression Omnibus: NCBI gene expression and hybridization array data repository.Nucleic Acids Res30207210

49. Parkinson H, Kapushesky M, Shojatalab M, Abeygunawardena N, Coulson R, et al. (2007) ArrayExpress–a public database of microarray experiments and gene expression profiles. Nucleic Acids Res 35: D747–750.H. ParkinsonM. KapusheskyM. ShojatalabN. AbeygunawardenaR. Coulson2007ArrayExpress–a public database of microarray experiments and gene expression profiles.Nucleic Acids Res35D747750

50. Hwang KB, Kong SW, Greenberg SA, Park PJ (2004) Combining gene expression data from different generations of oligonucleotide arrays. BMC Bioinformatics 5: 159.KB HwangSW KongSA GreenbergPJ Park2004Combining gene expression data from different generations of oligonucleotide arrays.BMC Bioinformatics5159

51. Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, et al. (2003) Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res 31: e15.RA IrizarryBM BolstadF. CollinLM CopeB. Hobbs2003Summaries of Affymetrix GeneChip probe level data.Nucleic Acids Res31e15