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

SPATIAL: A System-level PAThway Impact AnaLysis approach

The goal of pathway analysis is to identify the pathways that are significantly impacted when a biological system is perturbed, e.g. by a disease or drug. Current methods treat pathways as independent entities. However, many signals are constantly sent from one pathway to another, essentially linking all pathways into a global, system-wide complex. In this work, we propose a set of three pathway analysis methods based on the impact analysis, that performs a system-level analysis by considering all signals between pathways, as well as their overlaps. Briefly, the global system is modeled in two ways: (i) considering the inter-pathway interaction exchange for each individual pathways, and (ii) combining all individual pathways to form a global, system-wide graph. The third analysis method is a hybrid of these two models. The new methods were compared with DAVID, GSEA, GSA, PathNet, Crosstalk and SPIA on 23 GEO data sets involving 19 tissues investigated in 12 conditions. The results show that both the ranking and the P-values of the target pathways are substantially improved when the analysis considers the system-wide dependencies and interactions between pathways.

Efficient and biologically relevant consensus strategy for Parkinson's disease gene prioritization

BACKGROUND

The systemic information enclosed in microarray data encodes relevant clues to overcome the poorly understood combination of genetic and environmental factors in Parkinson's disease (PD), which represents the major obstacle to understand its pathogenesis and to develop disease-modifying therapeutics. While several gene prioritization approaches have been proposed, none dominate over the rest. Instead, hybrid approaches seem to outperform individual approaches.

METHODS

A consensus strategy is proposed for PD related gene prioritization from mRNA microarray data based on the combination of three independent prioritization approaches: Limma, machine learning, and weighted gene co-expression networks.

RESULTS

The consensus strategy outperformed the individual approaches in terms of statistical significance, overall enrichment and early recognition ability. In addition to a significant biological relevance, the set of 50 genes prioritized exhibited an excellent early recognition ability (6 of the top 10 genes are directly associated with PD). 40 % of the prioritized genes were previously associated with PD including well-known PD related genes such as SLC18A2, TH or DRD2. Eight genes (CCNH, DLK1, PCDH8, SLIT1, DLD, PBX1, INSM1, and BMI1) were found to be significantly associated to biological process affected in PD, representing potentially novel PD biomarkers or therapeutic targets. Additionally, several metrics of standard use in chemoinformatics are proposed to evaluate the early recognition ability of gene prioritization tools.

CONCLUSIONS

The proposed consensus strategy represents an efficient and biologically relevant approach for gene prioritization tasks providing a valuable decision-making tool for the study of PD pathogenesis and the development of disease-modifying PD therapeutics.

An Effective Method to Identify Shared Pathways and Common Factors among Neurodegenerative Diseases

Groups of distinct but related diseases often share common symptoms, which suggest likely overlaps in underlying pathogenic mechanisms. Identifying the shared pathways and common factors among those disorders can be expected to deepen our understanding for them and help designing new treatment strategies effected on those diseases. Neurodegeneration diseases, including Alzheimer's disease (AD), Parkinson's disease (PD) and Huntington's disease (HD), were taken as a case study in this research. Reported susceptibility genes for AD, PD and HD were collected and human protein-protein interaction network (hPPIN) was used to identify biological pathways related to neurodegeneration. 81 KEGG pathways were found to be correlated with neurodegenerative disorders. 36 out of the 81 are human disease pathways, and the remaining ones are involved in miscellaneous human functional pathways. Cancers and infectious diseases are two major subclasses within the disease group. Apoptosis is one of the most significant functional pathways. Most of those pathways found here are actually consistent with prior knowledge of neurodegenerative diseases except two cell communication pathways: adherens and tight junctions. Gene expression analysis showed a high probability that the two pathways were related to neurodegenerative diseases. A combination of common susceptibility genes and hPPIN is an effective method to study shared pathways involved in a group of closely related disorders. Common modules, which might play a bridging role in linking neurodegenerative disorders and the enriched pathways, were identified by clustering analysis. The identified shared pathways and common modules can be expected to yield clues for effective target discovery efforts on neurodegeneration.

Comparative study on gene set and pathway topology-based enrichment methods

Background

Enrichment analysis is a popular approach to identify pathways or sets of genes which are significantly enriched in the context of differentially expressed genes. The traditional gene set enrichment approach considers a pathway as a simple gene list disregarding any knowledge of gene or protein interactions. In contrast, the new group of so called pathway topology-based methods integrates the topological structure of a pathway into the analysis.

Methods

We comparatively investigated gene set and pathway topology-based enrichment approaches, considering three gene set and four topological methods. These methods were compared in two extensive simulation studies and on a benchmark of 36 real datasets, providing the same pathway input data for all methods.

Results

In the benchmark data analysis both types of methods showed a comparable ability to detect enriched pathways. The first simulation study was conducted with KEGG pathways, which showed considerable gene overlaps between each other. In this study with original KEGG pathways, none of the topology-based methods outperformed the gene set approach. Therefore, a second simulation study was performed on non-overlapping pathways created by unique gene IDs. Here, methods accounting for pathway topology reached higher accuracy than the gene set methods, however their sensitivity was lower.

Conclusions

We conducted one of the first comprehensive comparative works on evaluating gene set against pathway topology-based enrichment methods. The topological methods showed better performance in the simulation scenarios with non-overlapping pathways, however, they were not conclusively better in the other scenarios. This suggests that simple gene set approach might be sufficient to detect an enriched pathway under realistic circumstances. Nevertheless, more extensive studies and further benchmark data are needed to systematically evaluate these methods and to assess what gain and cost pathway topology information introduces into enrichment analysis. Both types of methods for enrichment analysis require further improvements in order to deal with the problem of pathway overlaps.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0751-5) contains supplementary material, which is available to authorized users.

14-3-3 inhibition promotes dopaminergic neuron loss and 14-3-3θ overexpression promotes recovery in the MPTP mouse model of Parkinson's disease

14-3-3s are a highly conserved protein family that plays important roles in cell survival and interact with several proteins implicated in Parkinson's disease (PD). Disruption of 14-3-3 expression and function has been implicated in the pathogenesis of PD. We have previously shown that increasing the expression level of 14-3-3θ is protective against rotenone and 1-methyl-4-phenylpyridinium (MPP(+)) in cultured cells. Here, we extend our studies to examine the effects of 14-3-3s in the 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) mouse model of PD. We first investigated whether targeted nigral 14-3-3θ overexpression mediated by adeno-associated virus offers neuroprotection against MPTP-induced toxicity. 14-3-3θ overexpression using this approach did not reduce MPTP-induced dopaminergic cell loss in the substantia nigra nor the depletion of dopamine (DA) and its metabolites in the striatum at three weeks after MPTP administration. However, 14-3-3θ-overexpressing mice showed a later partial recovery in striatal DA metabolites at eight weeks after MPTP administration compared to controls, suggesting that 14-3-3θ overexpression may help in the functional recovery of those dopaminergic neurons that survive. Conversely, we investigated whether disrupting 14-3-3 function in transgenic mice expressing the pan 14-3-3 inhibitor difopein exacerbates MPTP-induced toxicity. We found that difopein expression promoted dopaminergic cell loss in response to MPTP treatment. Together, these findings suggest that 14-3-3θ overexpression promotes recovery of DA metabolites whereas 14-3-3 inhibition exacerbates neuron loss in the MPTP mouse model of PD.

Proteotoxicity and cardiac dysfunction

Baseline physiological function of the mammalian heart is under the constant threat of environmental or intrinsic pathological insults. Cardiomyocyte proteins are thus subject to unremitting pressure to function optimally, and this depends on them assuming and maintaining proper conformation. This review explores the multiple defenses a cell may use for its proteins to assume and maintain correct protein folding and conformation. There are multiple quality control mechanisms to ensure that nascent polypeptides are properly folded and mature proteins maintain their functional conformation. When proteins do misfold, either in the face of normal or pathological stimuli or because of intrinsic mutations or post-translational modifications, they must either be refolded correctly or recycled. In the absence of these corrective processes, they may become toxic to the cell. Herein, we explore some of the underlying mechanisms that lead to proteotoxicity. The continued presence and chronic accumulation of misfolded or unfolded proteins can be disastrous in cardiomyocytes because these misfolded proteins can lead to aggregation or the formation of soluble peptides that are proteotoxic. This in turn leads to compromised protein quality control and precipitating a downward spiral of the cell's ability to maintain protein homeostasis. Some underlying mechanisms are discussed and the therapeutic potential of interfering with proteotoxicity in the heart is explored.

Evidence for Immune Response, Axonal Dysfunction and Reduced Endocytosis in the Substantia Nigra in Early Stage Parkinson's Disease

Subjects with incidental Lewy body disease (iLBD) may represent the premotor stage of Parkinson’s disease (PD). To elucidate molecular mechanisms underlying neuronal dysfunction and alpha-synuclein pathology in the premotor phase of PD, we investigated the transcriptome of the substantia nigra (SN) of well-characterized iLBD, PD donors and age-matched controls with Braak alpha-synuclein stage ranging from 0–6. In Braak alpha-synuclein stages 1 and 2, we observed deregulation of pathways linked to axonal degeneration, immune response and endocytosis, including axonal guidance signaling, mTOR signaling, EIF2 signaling and clathrin-mediated endocytosis in the SN. In Braak stages 3 and 4, we observed deregulation of pathways involved in protein translation and cell survival, including mTOR and EIF2 signaling. In Braak stages 5 and 6, we observed deregulation of dopaminergic signaling, axonal guidance signaling and thrombin signaling. Throughout the progression of PD pathology, we observed a deregulation of mTOR, EIF2 and regulation of eIF4 and p70S6K signaling in the SN. Our results indicate that molecular mechanisms related to axonal dysfunction, endocytosis and immune response are an early event in PD pathology, whereas mTOR and EIF2 signaling are impaired throughout disease progression. These pathways may hold the key to altering the disease progression in PD.

Moment based gene set tests

BACKGROUND

Permutation-based gene set tests are standard approaches for testing relationships between collections of related genes and an outcome of interest in high throughput expression analyses. Using M random permutations, one can attain p-values as small as 1/(M+1). When many gene sets are tested, we need smaller p-values, hence larger M, to achieve significance while accounting for the number of simultaneous tests being made. As a result, the number of permutations to be done rises along with the cost per permutation. To reduce this cost, we seek parametric approximations to the permutation distributions for gene set tests.

RESULTS

We study two gene set methods based on sums and sums of squared correlations. The statistics we study are among the best performers in the extensive simulation of 261 gene set methods by Ackermann and Strimmer in 2009. Our approach calculates exact relevant moments of these statistics and uses them to fit parametric distributions. The computational cost of our algorithm for the linear case is on the order of doing |G| permutations, where |G| is the number of genes in set G. For the quadratic statistics, the cost is on the order of |G|(2) permutations which can still be orders of magnitude faster than plain permutation sampling. We applied the permutation approximation method to three public Parkinson's Disease expression datasets and discovered enriched gene sets not previously discussed. We found that the moment-based gene set enrichment p-values closely approximate the permutation method p-values at a tiny fraction of their cost. They also gave nearly identical rankings to the gene sets being compared.

CONCLUSIONS

We have developed a moment based approximation to linear and quadratic gene set test statistics' permutation distribution. This allows approximate testing to be done orders of magnitude faster than one could do by sampling permutations. We have implemented our method as a publicly available Bioconductor package, npGSEA (www.bioconductor.org) .

Heat shock protein responses to aging and proteotoxicity in the olfactory bulb

The olfactory bulb is one of the most vulnerable brain regions in age-related proteinopathies. Proteinopathic stress is mitigated by the heat shock protein (Hsp) family of chaperones. Here, we describe age-related decreases in Hsc70 in the olfactory bulb of the female rat and higher levels of Hsp70 and Hsp25 in middle and old age than at 2-4 months. To model proteotoxic and oxidative stress in the olfactory bulb, primary olfactory bulb cultures were treated with the proteasome inhibitors lactacystin and MG132 or the pro-oxidant paraquat. Toxin-induced increases were observed in Hsp70, Hsp25, and Hsp32. To determine the functional consequences of the increase in Hsp70, we attenuated Hsp70 activity with two mechanistically distinct inhibitors. The Hsp70 inhibitors greatly potentiated the toxicity of sublethal lactacystin or MG132 but not of paraquat. Although ubiquitinated protein levels were unchanged with aging in vivo or with sublethal MG132 in vitro, there was a large, synergistic increase in ubiquitinated proteins when proteasome and Hsp70 functions were simultaneously inhibited. Our study suggests that olfactory bulb cells rely heavily on Hsp70 chaperones to maintain homeostasis during mild proteotoxic, but not oxidative insults, and that Hsp70 prevents the accrual of ubiquitinated proteins in these cells. The olfactory bulb is affected in the early phases of many age-related neurodegenerative disorders. Here, we described the impact of aging on multiple heat shock proteins (Hsps), such as Hsp70, in the female rat olfactory bulb in vivo. Using multiple proteasome and Hsp70 inhibitors (see schematic), we found that proteotoxicity elicited a compensatory increase in Hsp70 in primary olfactory bulb cells in vitro. Hsp70 then reduced the proteotoxic buildup of ubiquitinated proteins and robustly protected against cell death according to three independent viability assays. Thus, olfactory bulb neurons can mount impressive natural adaptations to proteotoxic injury, perhaps explaining why neurodegenerative disorders are so delayed in onset and so slow to progress.

Regulators of G-protein-signaling proteins: negative modulators of G-protein-coupled receptor signaling

Regulators of G-protein-signaling (RGS) proteins are a category of intracellular proteins that have an inhibitory effect on the intracellular signaling produced by G-protein-coupled receptors (GPCRs). RGS along with RGS-like proteins switch on through direct contact G-alpha subunits providing a variety of intracellular functions through intracellular signaling. RGS proteins have a common RGS domain that binds to G alpha. RGS proteins accelerate GTPase and thus enhance guanosine triphosphate hydrolysis through the alpha subunit of heterotrimeric G proteins. As a result, they inactivate the G protein and quickly turn off GPCR signaling thus terminating the resulting downstream signals. Activity and subcellular localization of RGS proteins can be changed through covalent molecular changes to the enzyme, differential gene splicing, and processing of the protein. Other roles of RGS proteins have shown them to not be solely committed to being inhibitors but behave more as modulators and integrators of signaling. RGS proteins modulate the duration and kinetics of slow calcium oscillations and rapid phototransduction and ion signaling events. In other cases, RGS proteins integrate G proteins with signaling pathways linked to such diverse cellular responses as cell growth and differentiation, cell motility, and intracellular trafficking. Human and animal studies have revealed that RGS proteins play a vital role in physiology and can be ideal targets for diseases such as those related to addiction where receptor signaling seems continuously switched on.

Heat shock protein defenses in the neocortex and allocortex of the telencephalon

The telencephalic allocortex develops protein inclusions before the neocortex in many age-related proteinopathies. One major defense mechanism against proteinopathic stress is the heat shock protein (Hsp) network. We therefore contrasted Hsp defenses in stressed primary neocortical and allocortical cells. Neocortical neurons were more resistant to the proteasome inhibitor MG132 than neurons from 3 allocortical subregions: entorhinal cortex, piriform cortex, and hippocampus. However, allocortical neurons exhibited higher MG132-induced increases in Hsp70 and heat shock cognate 70 (Hsc70). MG132-treated allocortical neurons also exhibited greater levels of protein ubiquitination. Inhibition of Hsp70/Hsc70 activity synergistically exacerbated MG132 toxicity in allocortical neurons more than neocortical neurons, suggesting that the allocortex is more reliant on these Hsp defenses. In contrast, astrocytes harvested from the neocortex or allocortex did not differ in their response to Hsp70/Hsc70 inhibition. Consistent with the idea that chaperones are maximally engaged in allocortical neurons, an increase in Hsp70/Hsc70 activity was protective only in neocortical neurons. Finally, the levels of select Hsps were altered in the neocortex and allocortex in vivo with aging.

Validity of the MPTP-Treated Mouse as a Model for Parkinson's Disease

Parkinson's disease (PD) is characterized by dopaminergic (DA) neuron death in the substantia nigra (SN) and subsequent striatal adaptations. Mice treated with the neurotoxin 1-methyl-4-phenyl-1,2,3,6-tetrahydropyrimidine (MPTP) are widely used as a model for PD. To assess the validity of the MPTP mouse model for PD pathogenesis, we here identify the biological processes that are dysregulated in both human PD and MPTP-treated mice. Gene enrichment analysis of published differentially expressed messenger RNAs (mRNAs) in the SN of PD patients and MPTP-treated mice revealed an enrichment of gene categories related to motor dysfunction and neurodegeneration. In the PD striatum, a similar enrichment was found, whereas in the striatum of MPTP mice, acute processes linked to epilepsy were selectively enriched shortly following MPTP treatment. More importantly, we integrated the proteins encoded by the differentially expressed mRNAs into molecular landscapes showing PD pathogenesis-implicated processes only in the SN, including vesicular trafficking, exocytosis, mitochondrial apoptosis, and DA neuron-specific transcription, but not in the striatum. We conclude that the current use of the MPTP mouse as a model for studying the molecular processes in PD pathogenesis is more valid for SN than striatal mechanisms in PD. This novel insight has important practical implications for future studies using this model to investigate PD pathogenesis and evaluate the efficacy of new treatments.

Computational deconvolution of genome wide expression data from Parkinson's and Huntington's disease brain tissues using population-specific expression analysis

The characterization of molecular changes in diseased tissues gives insight into pathophysiological mechanisms and is important for therapeutic development. Genome-wide gene expression analysis has proven valuable for identifying biological processes in neurodegenerative diseases using post mortem human brain tissue and numerous datasets are publically available. However, many studies utilize heterogeneous tissue samples consisting of multiple cell types, all of which contribute to global gene expression values, confounding biological interpretation of the data. In particular, changes in numbers of neuronal and glial cells occurring in neurodegeneration confound transcriptomic analyses, particularly in human brain tissues where sample availability and controls are limited. To identify cell specific gene expression changes in neurodegenerative disease, we have applied our recently published computational deconvolution method, population specific expression analysis (PSEA). PSEA estimates cell-type-specific expression values using reference expression measures, which in the case of brain tissue comprises mRNAs with cell-type-specific expression in neurons, astrocytes, oligodendrocytes and microglia. As an exercise in PSEA implementation and hypothesis development regarding neurodegenerative diseases, we applied PSEA to Parkinson's and Huntington's disease (PD, HD) datasets. Genes identified as differentially expressed in substantia nigra pars compacta neurons by PSEA were validated using external laser capture microdissection data. Network analysis and Annotation Clustering (DAVID) identified molecular processes implicated by differential gene expression in specific cell types. The results of these analyses provided new insights into the implementation of PSEA in brain tissues and additional refinement of molecular signatures in human HD and PD.

Dysfunction of endocytic kinase AAK1 in ALS

Mechanisms of human mutant superoxide dismutase 1 (SOD1)-induced toxicity in causing the familial form of amyotrophic lateral sclerosis (ALS) remain elusive. Identification of new proteins that can selectively interact with mutant SOD1s and investigation of their potential roles in ALS are important to discover new pathways that are involved in disease pathology. Using the yeast two-hybrid system, we identified the adaptor-associated kinase 1 (AAK1), a regulatory protein in clathrin-coated vesicle endocytic pathway that selectively interacted with the mutant but not the wild-type SOD1. Using both transgenic mouse and rat SOD1-linked familial ALS (FALS) models, we found that AAK1 was partially colocalized with the endosomal and presynaptic protein markers under the normal physiological condition, but was mislocated into aggregates that contained mutant SOD1s and the neurofilament proteins in rodent models of ALS in disease. AAK1 protein levels were also decreased in ALS patients. These results suggest that dysfunction of a component in the endosomal and synaptic vesicle recycling pathway is involved in ALS pathology.

Differential Co-Expression between α-Synuclein and IFN-γ Signaling Genes across Development and in Parkinson's Disease

Expression patterns of the alpha-synuclein gene (SNCA) were studied across anatomy, development, and disease to better characterize its role in the brain. In this postmortem study, negative spatial co-expression between SNCA and 73 interferon-γ (IFN-γ) signaling genes was observed across many brain regions. Recent animal studies have demonstrated that IFN-γ induces loss of dopamine neurons and nigrostriatal degeneration. This opposing pattern between SNCA and IFN-γ signaling genes increases with age (rho = −0.78). In contrast, a meta-analysis of four microarray experiments representing 126 substantia nigra samples reveals a switch to positive co-expression in Parkinson’s disease (p<0.005). Use of genome-wide testing demonstrates this relationship is specific to SNCA (p<0.002). This change in co-expression suggests an immunomodulatory role of SNCA that may provide insight into neurodegeneration. Genes showing similar co-expression patterns have been previously linked to Alzheimer’s (ANK1) and Parkinson’s disease (UBE2E2, PCMT1, HPRT1 and RIT2).

Comparative pathway and network analysis of brain transcriptome changes during adult aging and in Parkinson's disease

Aging is considered as one of the main factors promoting the risk for Parkinson's disease (PD), and common mechanisms of dopamine neuron degeneration in aging and PD have been proposed in recent years. Here, we use a statistical meta-analysis of human brain transcriptomics data to investigate potential mechanistic relationships between adult brain aging and PD pathogenesis at the pathway and network level. The analyses identify statistically significant shared pathway and network alterations in aging and PD and an enrichment in PD-associated sequence variants from genome-wide association studies among the jointly deregulated genes. We find robust discriminative patterns for groups of functionally related genes with potential applications as combined risk biomarkers to detect aging- and PD-linked oxidative stress, e.g., a consistent over-expression of metallothioneins matching with findings in previous independent studies. Interestingly, analyzing the regulatory network and mouse knockout expression data for NR4A2, a transcription factor previously associated with rare mutations in PD and here found as the most significantly under-expressed gene in PD among the jointly altered genes, suggests that aging-related NR4A2 expression changes may increase PD risk via downstream effects similar to disease-linked mutations and to expression changes in sporadic PD. Overall, the analyses suggest mechanistic explanations for the age-dependence of PD risk and reveal significant and robust shared process alterations with potential applications in biomarker development for pre-symptomatic risk assessment or early stage diagnosis.

Analysis of Parkinson's disease pathophysiology using an integrated genomics-bioinformatics approach

The pathogenesis and pathophysiology of a disease determine how it should be diagnosed and treated. Yet, understanding the cause and mechanisms of progression often requires intensive human efforts, especially for diseases with complex etiology. The latest genomic technology coupled with advanced, large-scale data analysis in the field known as bioinformatics has promised a high-throughput approach that can quickly identify disease-affected genes and pathways by examining tissue samples collected from patients and control subjects. Furthermore, significant biological themes indicative of genomic events can be recognized on the basis of affected genes. However, given identified biological themes, it is not clear how to organize genomic events to arrive at a coherent pathophysiological explanation about the disease. To address this important issue, we have developed an innovative method named "Expression Data Up-Stream Analysis" (EDUSA) that can perform a bioinformatics analysis to identify and rank upstream processes effectively. We applied it to Parkinson's disease (PD) using a genomic data set available at a public data repository known as Gene Expression Omnibus (GEO). In this study, disease-affected genes were identified using GEO2R software, and disease-pertinent processes were identified using EASE software. Then the EDUSA program was used to determine the upstream versus downstream hierarchy of the processes. The results confirmed the current misfolded protein theory about the pathogenesis of PD, and provided new insights as well. Particularly, our program discovered that RNA (ribonucleic acid) metabolism pathology was a potential cause of PD, which in fact, is an emerging theory of neurodegenerative disorders. In addition, it was found that the dysfunction of the transport system seemed to occur in the early phase of neurodegeneration, whereas mitochondrial dysfunction appeared at a later stage. Using this methodology, we have demonstrated how to determine the stages of disease development with single-point data collection.

Heat shock proteins in neurodegenerative disorders and aging

Many members of the heat shock protein family act in unison to refold or degrade misfolded proteins. Some heat shock proteins also directly interfere with apoptosis. These homeostatic functions are especially important in proteinopathic neurodegenerative diseases, in which specific proteins misfold, aggregate, and kill cells through proteotoxic stress. Heat shock protein levels may be increased or decreased in these disorders, with the direction of the response depending on the individual heat shock protein, the disease, cell type, and brain region. Aging is also associated with an accrual of proteotoxic stress and modulates expression of several heat shock proteins. We speculate that the increase in some heat shock proteins in neurodegenerative conditions may be partly responsible for the slow progression of these disorders, whereas the increase in some heat shock proteins with aging may help delay senescence. The protective nature of many heat shock proteins in experimental models of neurodegeneration supports these hypotheses. Furthermore, some heat shock proteins appear to be expressed at higher levels in longer-lived species. However, increases in heat shock proteins may be insufficient to override overwhelming proteotoxic stress or reverse the course of these conditions, because the expression of several other heat shock proteins and endogenous defense systems is lowered. In this review we describe a number of stress-induced changes in heat shock proteins as a function of age and neurodegenerative pathology, with an emphasis on the heat shock protein 70 (Hsp70) family and the two most common proteinopathic disorders of the brain, Alzheimer's and Parkinson's disease.

Transcription factors Foxa1 and Foxa2 are required for adult dopamine neurons maintenance

The proteins Foxa1 and Foxa2 belong to the forkhead family of transcription factors and are involved in the development of several tissues, including liver, pancreas, lung, prostate, and the neural system. Both Foxa1 and Foxa2 are also crucial for the specification and differentiation of dopamine (DA) neurons during embryonic development, while about 30% of mice with an embryonic deletion of a single allele of the Foxa2 gene exhibit an age-related asymmetric loss of DA neurons and develop locomotor symptoms resembling Parkinson's disease (PD). Notably, both Foxa1 and Foxa2 factors continue to be expressed in the adult dopamine system. To directly assess their functions selectively in adult DA neurons, we induced genetic deletions of Foxa1/2 transcription factors in mice using a tamoxifen inducible tissue-specific CreERT2 recombinase expressed under control of the dopamine transporter (DAT) promoter (DATCreERT2). The conditional DA neurons-specific ablation of both genes, but not of Foxa2 alone, in early adulthood, caused a decline of striatal dopamine and its metabolites, along with locomotor deficits. At early pre-symptomatic stages, we observed a decline in aldehyde dehydrogenase family 1, subfamily A1 (Aldh1a1) protein expression in DA neurons. Further analyses revealed a decline of aromatic amino acid decarboxylase (AADC) and a complete loss of DAT expression in these neurons. These molecular changes ultimately led to a reduction of DA neuron numbers in the substantia nigra pars compacta (SNpc) of aged cFoxa1/2^−/− mice, resembling the progressive course of PD in humans. Altogether, in this study, we address the molecular, cellular, and functional role of both Foxa1 and Foxa2 factors in the maintenance of the adult dopamine system which may help to find better approaches for PD treatment.

Integrated multi-cohort transcriptional meta-analysis of neurodegenerative diseases

Introduction

Neurodegenerative diseases share common pathologic features including neuroinflammation, mitochondrial dysfunction and protein aggregation, suggesting common underlying mechanisms of neurodegeneration. We undertook a meta-analysis of public gene expression data for neurodegenerative diseases to identify a common transcriptional signature of neurodegeneration.

Results

Using 1,270 post-mortem central nervous system tissue samples from 13 patient cohorts covering four neurodegenerative diseases, we identified 243 differentially expressed genes, which were similarly dysregulated in 15 additional patient cohorts of 205 samples including seven neurodegenerative diseases. This gene signature correlated with histologic disease severity. Metallothioneins featured prominently among differentially expressed genes, and functional pathway analysis identified specific convergent themes of dysregulation. MetaCore network analyses revealed various novel candidate hub genes (e.g. STAU2). Genes associated with M1-polarized macrophages and reactive astrocytes were strongly enriched in the meta-analysis data. Evaluation of genes enriched in neurons revealed 70 down-regulated genes, over half not previously associated with neurodegeneration. Comparison with aging brain data (3 patient cohorts, 221 samples) revealed 53 of these to be unique to neurodegenerative disease, many of which are strong candidates to be important in neuropathogenesis (e.g. NDN, NAP1L2). ENCODE ChIP-seq analysis predicted common upstream transcriptional regulators not associated with normal aging (REST, RBBP5, SIN3A, SP2, YY1, ZNF143, IKZF1). Finally, we removed genes common to neurodegeneration from disease-specific gene signatures, revealing uniquely robust immune response and JAK-STAT signaling in amyotrophic lateral sclerosis.

Conclusions

Our results implicate pervasive bioenergetic deficits, M1-type microglial activation and gliosis as unifying themes of neurodegeneration, and identify numerous novel genes associated with neurodegenerative processes.

Electronic supplementary material

The online version of this article (doi:10.1186/s40478-014-0093-y) contains supplementary material, which is available to authorized users.

Systems-based analyses of brain regions functionally impacted in Parkinson's disease reveals underlying causal mechanisms

Detailed analysis of disease-affected tissue provides insight into molecular mechanisms contributing to pathogenesis. Substantia nigra, striatum, and cortex are functionally connected with increasing degrees of alpha-synuclein pathology in Parkinson's disease. We undertook functional and causal pathway analysis of gene expression and proteomic alterations in these three regions, and the data revealed pathways that correlated with disease progression. In addition, microarray and RNAseq experiments revealed previously unidentified causal changes related to oligodendrocyte function and synaptic vesicle release, and these and other changes were reflected across all brain regions. Importantly, subsets of these changes were replicated in Parkinson's disease blood; suggesting peripheral tissue may provide important avenues for understanding and measuring disease status and progression. Proteomic assessment revealed alterations in mitochondria and vesicular transport proteins that preceded gene expression changes indicating defects in translation and/or protein turnover. Our combined approach of proteomics, RNAseq and microarray analyses provides a comprehensive view of the molecular changes that accompany functional loss and alpha-synuclein pathology in Parkinson's disease, and may be instrumental to understand, diagnose and follow Parkinson's disease progression.

Preventing α-synuclein aggregation: the role of the small heat-shock molecular chaperone proteins

Protein homeostasis, or proteostasis, is the process of maintaining the conformational and functional integrity of the proteome. The failure of proteostasis can result in the accumulation of non-native proteins leading to their aggregation and deposition in cells and in tissues. The amyloid fibrillar aggregation of the protein α-synuclein into Lewy bodies and Lewy neuritis is associated with neurodegenerative diseases classified as α-synucleinopathies, which include Parkinson's disease and dementia with Lewy bodies. The small heat-shock proteins (sHsps) are molecular chaperones that are one of the cell's first lines of defence against protein aggregation. They act to stabilise partially folded protein intermediates, in an ATP-independent manner, to maintain cellular proteostasis under stress conditions. Thus, the sHsps appear ideally suited to protect against α-synuclein aggregation, yet these fail to do so in the context of the α-synucleinopathies. This review discusses how sHsps interact with α-synuclein to prevent its aggregation and, in doing so, highlights the multi-faceted nature of the mechanisms used by sHsps to prevent the fibrillar aggregation of proteins. It also examines what factors may contribute to α-synuclein escaping the sHsp chaperones in the context of the α-synucleinopathies.

'Neuroinflammation' differs categorically from inflammation: transcriptomes of Alzheimer's disease, Parkinson's disease, schizophrenia and inflammatory diseases compared

'Neuroinflammation' has become a widely applied term in the basic and clinical neurosciences but there is no generally accepted neuropathological tissue correlate. Inflammation, which is characterized by the presence of perivascular infiltrates of cells of the adaptive immune system, is indeed seen in the central nervous system (CNS) under certain conditions. Authors who refer to microglial activation as neuroinflammation confuse this issue because autoimmune neuroinflammation serves as a synonym for multiple sclerosis, the prototypical inflammatory disease of the CNS. We have asked the question whether a data-driven, unbiased in silico approach may help to clarify the nomenclatorial confusion. Specifically, we have examined whether unsupervised analysis of microarray data obtained from human cerebral cortex of Alzheimer's, Parkinson's and schizophrenia patients would reveal a degree of relatedness between these diseases and recognized inflammatory conditions including multiple sclerosis. Our results using two different data analysis methods provide strong evidence against this hypothesis demonstrating that very different sets of genes are involved. Consequently, the designations inflammation and neuroinflammation are not interchangeable. They represent different categories not only at the histophenotypic but also at the transcriptomic level. Therefore, non-autoimmune neuroinflammation remains a term in need of definition.

Impact of aging on heat shock protein expression in the substantia nigra and striatum of the female rat

Many heat shock proteins are chaperones that help refold or degrade misfolded proteins and battle apoptosis. Because of their capacity to protect against protein misfolding, they may help keep diseases of aging at bay. A few reports have examined heat shock proteins (eg. Hsp25, Hsp60, Hsp70, and heat shock cognate 70 or Hsc70) as a function of age in the striatum and nigra. In the present study, we examined the impact of aging on Hsp25, heme oxygenase 1 (HO1 or Hsp32), Hsp40, Hsp60, Hsc70, Hsc/Hsp70 interacting protein (Hip), 78 kDa glucose-regulated protein (GRP78), Hsp90, and ubiquitinated proteins in the nigra and striatum of the female rat by infrared immunoblotting. Female animals are not typically examined in aging studies, adding further to the novelty of our study. Striatal HO1 and Hsp40 were both higher in middle-aged females than in the oldest group. Hsp60 levels were also highest in middle age in the nigra, but were highest in the oldest animals in the striatum. Striatal levels of Hsc70 and the co-chaperone Hip were lower in the oldest group relative to the youngest animals. In contrast, Hsp25 rose with advancing age in both regions. Hsp25 was also colocalized with tyrosine hydroxylase in nigral neurons. Ubiquitinated proteins exhibited a trend to rise in the oldest animals in both regions, and K48 linkage-specific ubiquitin rose significantly from 4-6 to 16-19 months in the striatum. Our study reveals a complex array of age-related changes in heat shock proteins. Furthermore, the age-related rises in some proteins, such as Hsp25, may reflect endogenous adaptations to cellular stress.

Tumor Necrosis Factor-Associated Protein 1 (TRAP1) is Released from the Mitochondria Following 6-hydroxydopamine Treatment

Parkinson's disease (PD) is a common neurodegenerative disorder characterized by progressive degeneration of dopaminergic neurons in the substantia nigra pars compacta. Most cases are sporadic and its etiology is incompletely understood. However, increasing evidence suggests that oxidative stress and mitochondrial dysfunction may be involved in the pathogenesis of Parkinson's disease. The aim of this study was to investigate changes in mitochondrial protein profiles during dopaminergic neuronal cell death using two-dimensional gel electrophoresis in conjunction with mass spectrometry. Several protein spots were found to be significantly altered following treatment of MN9D dopaminergic neuronal cells with 6-hydroxydopamine (6-OHDA). Among several identified candidates, TNF receptor-associated protein 1 (TRAP1), a mitochondrial molecular chaperone, was released from the mitochondria into the cytosol in MN9D cells as well as primary cultures of dopaminergic neurons following 6-OHDA treatment. This event was drug-specific in that such apoptotic inducers as staurosporine and etoposide did not cause translocation of TRAP1 into the cytosol. To our knowledge, the present study is the first to demonstrate the drug-induced subcellular translocation of TRAP1 during neurodegeneration. Further studies delineating cellular mechanism associated with this phenomenon and its functional consequence may provide better understanding of dopaminergic neurodegeneration that underlies PD pathogenesis.

Orchestrated increase of dopamine and PARK mRNAs but not miR-133b in dopamine neurons in Parkinson's disease

Progressive loss of substantia nigra dopamine neurons (SN DA) is a hallmark of aging and of Parkinson's disease (PD). Mutations in PARK genes cause familial PD forms. Increased expression of alpha-synuclein (PARK4) is a disease-triggering event in familial PD and also observed in SN DA neurons in sporadic PD but related transcriptional changes are unknown. With optimized single-cell quantitative real-time polymerase chain reaction analysis, we compared messenger RNA and microRNA levels in SN DA neurons from sporadic PD patients and controls. Non-optimally matched donor ages and RNA integrities are common problems when analyzing human samples. We dissected the influence of distinct ages and RNA integrities of our samples by applying a specifically-optimized, linear-mixed-effects model to quantitative real-time polymerase chain reaction-data. We identified that elevated alpha-synuclein messenger RNA levels in SN DA neurons of human PD brains were positively correlated with corresponding elevated levels of mRNAs for functional compensation of progressive SN DA loss and for enhanced proteasomal (PARK5/UCHL1) and lysosomal (PARK9/ATPase13A2) function, possibly counteracting alpha-synuclein toxicity. In contrast, microRNA miR-133b levels, previously implicated in transcriptional dysregulation in PD, were not altered in SN DA neurons in PD.

Genomic convergence and network analysis approach to identify candidate genes in Alzheimer's disease

Background

Alzheimer’s disease (AD) is one of the leading genetically complex and heterogeneous disorder that is influenced by both genetic and environmental factors. The underlying risk factors remain largely unclear for this heterogeneous disorder. In recent years, high throughput methodologies, such as genome-wide linkage analysis (GWL), genome-wide association (GWA) studies, and genome-wide expression profiling (GWE), have led to the identification of several candidate genes associated with AD. However, due to lack of consistency within their findings, an integrative approach is warranted. Here, we have designed a rank based gene prioritization approach involving convergent analysis of multi-dimensional data and protein-protein interaction (PPI) network modelling.

Results

Our approach employs integration of three different AD datasets- GWL,GWA and GWE to identify overlapping candidate genes ranked using a novel cumulative rank score (S_R) based method followed by prioritization using clusters derived from PPI network. S_R for each gene is calculated by addition of rank assigned to individual gene based on either p value or score in three datasets. This analysis yielded 108 plausible AD genes. Network modelling by creating PPI using proteins encoded by these genes and their direct interactors resulted in a layered network of 640 proteins. Clustering of these proteins further helped us in identifying 6 significant clusters with 7 proteins (EGFR, ACTB, CDC2, IRAK1, APOE, ABCA1 and AMPH) forming the central hub nodes. Functional annotation of 108 genes revealed their role in several biological activities such as neurogenesis, regulation of MAP kinase activity, response to calcium ion, endocytosis paralleling the AD specific attributes. Finally, 3 potential biochemical biomarkers were found from the overlap of 108 AD proteins with proteins from CSF and plasma proteome. EGFR and ACTB were found to be the two most significant AD risk genes.

Conclusions

With the assumption that common genetic signals obtained from different methodological platforms might serve as robust AD risk markers than candidates identified using single dimension approach, here we demonstrated an integrated genomic convergence approach for disease candidate gene prioritization from heterogeneous data sources linked to AD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-199) contains supplementary material, which is available to authorized users.

Integrating pathways of Parkinson's disease in a molecular interaction map

Parkinson's disease (PD) is a major neurodegenerative chronic disease, most likely caused by a complex interplay of genetic and environmental factors. Information on various aspects of PD pathogenesis is rapidly increasing and needs to be efficiently organized, so that the resulting data is available for exploration and analysis. Here we introduce a computationally tractable, comprehensive molecular interaction map of PD. This map integrates pathways implicated in PD pathogenesis such as synaptic and mitochondrial dysfunction, impaired protein degradation, alpha-synuclein pathobiology and neuroinflammation. We also present bioinformatics tools for the analysis, enrichment and annotation of the map, allowing the research community to open new avenues in PD research. The PD map is accessible at http://minerva.uni.lu/pd_map .

The mitochondrial rhomboid protease: its rise from obscurity to the pinnacle of disease-relevant genes

The Rhomboid proteases belong to a highly conserved family of proteins that are present in all branches of life. In Drosophila, the secretory pathway-localized rhomboid proteases are crucial for epidermal growth factor (EGF) signaling. The identification of a mitochondrial-localized rhomboid protease shed light on other functions of rhomboid proteases including the maintenance of mitochondrial morphology and the regulation of apoptosis. More recent work has revealed other functions of the mitochondrial rhomboid protease in mitochondrial and cellular biology, failure of which have been implicated in human diseases. In this review, we will summarize the current knowledge and disease relevance of the mitochondrial-localized rhomboid protease. This article is part of a Special Issue entitled: Intramembrane Proteases.

Gene expression profiling analysis of the putamen for the investigation of compensatory mechanisms in Parkinson's disease

BACKGROUND

Parkinson's disease (PD) is affecting 5 million people worldwide, but the response mechanisms of the striatum are still unclear. Therefore, identification of gene expression alterations in the striatum will greatly assist the development of novel therapy strategies.

METHODS

We performed a comprehensive gene expression analysis in 15 PD patients and 15 normal controls to identify differentially expressed genes (DEGs) using the expression profile GSE20291 from Gene Expression Omnibus (GEO). Gene Ontology (GO) analysis and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway enrichment analysis were used to define functions and pathways altered in PD. Protein-protein interaction network was constructed to find out the modules with close interactions.

RESULTS

Total715 DEGs including 268 up-regulated and 447 down-regulated genes were obtained. GO functional enrichment analysis indicated that the genes related with neurons function and cell morphogenesis might be changed upon PD. KEGG pathway enrichment analysis showed that most of the genes were enriched in the nodes of Gap junction, calcium signaling pathway, phosphatidylinositol signaling system, long-term potentiation, Alzheimer's disease and GnRH signaling pathway. Protein-protein interaction network and module analysis suggested that some apoptosis related genes, such as PRKCA, CDC42 and BCL2 may play critical roles in striatal neurons growth.

CONCLUSION

Intrinsic striatal tyrosine hydroxylase interneurons growth may be promoted by changes in several genes expression and thus reduce the functional excitatory synapses.

SPIEDw: a searchable platform-independent expression database web tool

Background

SPIEDw is a web tool designed to facilitate fast and simple quantitative querying of publically available gene expression data. The resource is motivated by the observation that transcriptional profiles can serve as effective means of comparing biological states across a wide set of experiments.

Results

Gene expression data for over 200,000 experiments across multiple species and platforms have been collected into a searchable database. The new resource is a development of the previously published SPIED, which was designed to be downloaded and queried locally with SPIED software. SPIEDw features three significant improvements over the original version. Firstly, the number of experiments covered has been doubled and now includes Agilent and Illumina technologies. Secondly, SPIEDw has an enhanced search algorithm for speedy web-based querying and lastly an abridged dataset comprising the most regulated genes has been included for a speedier search and searching for enrichment of gene sets.

Conclusions

SPIEDw is simple to use, not requiring any expertise in microarray analysis, and the output straightforward to interpret. It is hoped that this will open up gene expression data mining to the wider research community.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-14-765) contains supplementary material, which is available to authorized users.

Continuously tunable Ca(2+) regulation of RNA-edited CaV1.3 channels

CaV1.3 ion channels are dominant Ca(2+) portals into pacemaking neurons, residing at the epicenter of brain rhythmicity and neurodegeneration. Negative Ca(2+) feedback regulation of CaV1.3 channels (CDI) is therefore critical for Ca(2+) homeostasis. Intriguingly, nearly half the CaV1.3 transcripts in the brain are RNA edited to reduce CDI and influence oscillatory activity. It is then mechanistically remarkable that this editing occurs precisely within an IQ domain, whose interaction with Ca(2+)-bound calmodulin (Ca(2+)/CaM) is believed to induce CDI. Here, we sought the mechanism underlying the altered CDI of edited channels. Unexpectedly, editing failed to attenuate Ca(2+)/CaM binding. Instead, editing weakened the prebinding of Ca(2+)-free CaM (apoCaM) to channels, which proves essential for CDI. Thus, editing might render CDI continuously tunable by fluctuations in ambient CaM, a prominent effect we substantiate in substantia nigral neurons. This adjustability of Ca(2+) regulation by CaM now looms as a key element of CNS Ca(2+) homeostasis.

A common structural scaffold in CTD phosphatases that supports distinct catalytic mechanisms

The phosphorylation and dephosphorylation of the carboxyl-terminal domain (CTD) of the largest RNA polymerase II (RNAPII) subunit is a critical regulatory checkpoint for transcription and mRNA processing. This CTD is unique to eukaryotic organisms and it contains multiple tandem-repeats with the consensus sequence Tyr(1) -Ser(2) -Pro(3) -Thr(4) -Ser(5) -Pro(6) -Ser(7) . Traditionally, CTD phosphatases that use metal-ion-independent (cysteine-based) and metal-ion-assisted (aspartate-based) catalytic mechanisms have been considered to belong to two independent groups. However, using structural comparisons we have identified a common structural scaffold in these two groups of CTD phosphatases. This common scaffold accommodates different catalytic processes with the same substrate specificity, in this case phospho-serine/threonine residues flanked by prolines. Furthermore, this scaffold provides a structural connection between two groups of protein tyrosine phosphatases (PTPs): Cys-based (classes I, II, and III) and Asp-based (class IV) PTPs. Redundancy in catalytic mechanisms is not infrequent and may arise in specific biological settings. To better understand the activity of the CTD phosphatases, we combined our structural analyses with data on CTD phosphatase expression in different human and mouse tissues. The results suggest that aspartate- and cysteine-based CTD-dephosphorylation acts in concert during cellular stress, when high levels of reactive oxygen species can inhibit the nucleophilic function of the catalytic cysteine, as occurs in mental and neurodegenerative disorders like schizophrenia, Alzheimer's and Parkinson's diseases. Moreover, these findings have significant implications for the study of the RNAPII-CTD dephosphorylation in eukaryotes.

The impact of pathogenic mitochondrial DNA mutations on substantia nigra neurons

Mitochondrial defects within substantia nigra (SN) neurons are implicated in the pathogenesis of Parkinson's disease. SN neurons show increased mitochondrial defects, mitochondrial DNA deletion levels, and susceptibility to such dysfunction, although the role of mitochondria in neuronal degeneration remains uncertain. In this study, we addressed this important question by exploring changes within the mitochondria of SN neurons from patients with primary mitochondrial diseases to determine whether mitochondrial dysfunction leads directly to neuronal cell loss. We counted the pigmented neurons and quantified mitochondrial respiratory activity, deficiencies in mitochondrial proteins, and the percentage of pathogenic mutations in single neurons. We found evidence of defects of both complex I and complex IV of the respiratory chain in all patients. We found that marked neuronal cell loss was only observed in a few patients with mitochondrial disease and that all these patients had mutations in polymerase gamma (POLG), which leads to the formation of multiple mitochondrial DNA deletions over time, similar to aging and Parkinson's disease. Interestingly, we detected α-synuclein pathology in two mitochondrial patients with POLG mutations. Our observations highlight the complex relationship between mitochondrial dysfunction and the susceptibility of SN neurons to degeneration and α-synuclein pathology. Our finding that the loss of SN neurons was only severe in patients with POLG mutations suggests that acquired mitochondrial defects may be less well tolerated by SN neurons than by inherited ones.

Adaptation and sensitization to proteotoxic stress

Although severe stress can elicit toxicity, mild stress often elicits adaptations. Here we review the literature on stress-induced adaptations versus stress sensitization in models of neurodegenerative diseases. We also describe our recent findings that chronic proteotoxic stress can elicit adaptations if the dose is low but that high-dose proteotoxic stress sensitizes cells to subsequent challenges. In these experiments, long-term, low-dose proteasome inhibition elicited protection in a superoxide dismutase-dependent manner. In contrast, acute, high-dose proteotoxic stress sensitized cells to subsequent proteotoxic challenges by eliciting catastrophic loss of glutathione. However, even in the latter model of synergistic toxicity, several defensive proteins were upregulated by severe proteotoxicity. This led us to wonder whether high-dose proteotoxic stress can elicit protection against subsequent challenges in astrocytes, a cell type well known for their resilience. In support of this new hypothesis, we found that the astrocytes that survived severe proteotoxicity became harder to kill. The adaptive mechanism was glutathione dependent. If these findings can be generalized to the human brain, similar endogenous adaptations may help explain why neurodegenerative diseases are so delayed in appearance and so slow to progress. In contrast, sensitization to severe stress may explain why defenses eventually collapse in vulnerable neurons.

α-Synuclein overexpression represses 14-3-3θ transcription

Previous gene microarray studies have shown that expression of 14-3-3θ is significantly decreased in an α-synuclein transgenic mouse model. In this study, we tested whether α-synuclein can regulate 14-3-3θ transcription. We demonstrate that the 14-3-3θ mRNA level is decreased in SH-SY5Y cells overexpressing α-synuclein. Luciferase activity under the control of the 14-3-3θ promoter is reduced both in stable SH-SY5Y cells constitutively overexpressing α-synuclein and in doxycycline-inducible SH-SY5Y cells upon α-synuclein induction, suggesting that the regulation of 14-3-3θ by α-synuclein occurs at the transcriptional level. Knockdown of α-synuclein by RNA interference does not increase the 14-3-3θ mRNA level. These findings suggest that α-synuclein represses 14-3-3θ transcription under pathologic conditions, but that regulation of 14-3-3θ expression is not a function of endogenous α-synuclein at baseline.

Distinctive patterns of DNA methylation associated with Parkinson disease: identification of concordant epigenetic changes in brain and peripheral blood leukocytes

Parkinson disease (PD) is a multifactorial neurodegenerative disorder with high incidence in the elderly, where environmental and genetic factors are involved in etiology. In addition, epigenetic mechanisms, including deregulation of DNA methylation have been recently associated to PD. As accurate diagnosis cannot be achieved pre-mortem, identification of early pathological changes is crucial to enable therapeutic interventions before major neuropathological damage occurs. Here we investigated genome-wide DNA methylation in brain and blood samples from PD patients and observed a distinctive pattern of methylation involving many genes previously associated to PD, therefore supporting the role of epigenetic alterations as a molecular mechanism in neurodegeneration. Importantly, we identified concordant methylation alterations in brain and blood, suggesting that blood might hold promise as a surrogate for brain tissue to detect DNA methylation in PD and as a source for biomarker discovery.

Decreased parvalbumin mRNA expression in dorsolateral prefrontal cortex in Parkinson's disease

It has recently been shown that expression of the rate-limiting GABA-synthesizing enzyme glutamic acid decarboxylase (GAD) is decreased in Brodmann area 9 (BA9) of the dorsolateral prefrontal cortex (DLPFC) in Parkinson's disease (PD) compared to control brains (Lanoue, A.C., Dumitriu, A., Myers, R.H., Soghomonian, JJ., 2010. Exp. Neurol. 206 (1), 207-217). A subpopulation of cortical GABAergic interneurons expresses the calcium-binding protein parvalbumin and plays a critical role in the control of pyramidal neuron excitability and the generation of cortical gamma frequency oscillations. In view of its key role in the physiology of the cerebral cortex, we sought to determine whether the expression of parvalbumin and the number of parvalbumin-expressing neurons are altered in BA9 of PD brains. First, isotopic in situ hybridization histochemistry was used to examine mRNA expression of parvalbumin on post-mortem brain sections. Second, immunohistochemistry and design-based stereology were used to determine the density of parvalbumin-positive interneurons in BA9. Quantification of mRNA labeling at the single cell level showed a significant decrease in parvalbumin expression in PD cases. In contrast, neuronal density of parvalbumin-positive neurons was not significantly different between PD and controls. Results confirm that the GABAergic system is altered in the DLPFC in PD and identify the contribution of parvalbumin-expressing neurons in these alterations. We speculate that these effects could contribute to altered cortical excitability and oscillatory activity previously documented in PD.

A network view on Parkinson's disease

Network-based systems biology tools including Pathway Studio 9.0 were used to identify Parkinson's disease (PD) critical molecular players, drug targets, and underlying biological processes. Utilizing several microarray gene expression datasets, biomolecular networks such as direct interaction, shortest path, and microRNA regulatory networks were constructed and analyzed for the disease conditions. Network topology analysis of node connectivity and centrality revealed in combination with the guilt-by-association rule 17 novel genes of PD-potential interest. Seven new microRNAs (miR-132, miR-133a1, miR-181-1, miR-182, miR-218-1, miR-29a, and miR-330) related to Parkinson's disease were identified, along with more microRNA targeted genes of interest like RIMS3, SEMA6D and SYNJ1. David and IPA enrichment analysis of KEGG and canonical pathways provided valuable mechanistic information emphasizing among others the role of chemokine signaling, adherence junction, and regulation of actin cytoskeleton pathways. Several routes for possible disease initiation and neuro protection mechanisms triggered via the extra-cellular ligands such as CX3CL1, SEMA6D and IL12B were thus uncovered, and a dual regulatory system of integrated transcription factors and microRNAs mechanisms was detected.

Applying transcriptomic and proteomic knowledge to Parkinson's disease drug discovery

It is recognised that in both genetic and sporadic cases of Parkinson's disease (PD), the basis of its etiopathology resides in the particular vulnerability of the dopaminergic neurons of the substantia nigra pars compacta (SNpc) to oxidative stress and in the failure to adequately remove abnormal proteins. These observations have been confirmed recently by microarray transcriptomic studies in human SN from PD brains and have extended understanding of the molecular pathways underlying the PD pathology. This article reviews recent gene expression profiling studies in sporadic PD postmortem SN and highlights gene candidates as putative molecular signatures for early disease diagnosis. In addition, the application of transcriptomics and proteomics in the quest for multifunctional neuroprotective-neurorescue drugs that might possess disease-modifying action is discussed.

Phenotypic characterization of retinoic acid differentiated SH-SY5Y cells by transcriptional profiling

Multiple genetic and environmental factors play a role in the development and progression of Parkinson's disease (PD). The main neuropathological hallmark of PD is the degeneration of dopaminergic (DAergic) neurons in the substantia nigra pars compacta. To study genetic and molecular contributors to the disease process, there is a great need for readily accessible cells with prominent DAergic features that can be used for reproducible in vitro cellular screening. Here, we investigated the molecular phenotype of retinoic acid (RA) differentiated SH-SY5Y cells using genome wide transcriptional profiling combined with gene ontology, transcription factor and molecular pathway analysis. We demonstrated that RA induces a general neuronal differentiation program in SH-SY5Y cells and that these cells develop a predominantly mature DAergic-like neurotransmitter phenotype. This phenotype is characterized by increased dopamine levels together with a substantial suppression of other neurotransmitter phenotypes, such as those for noradrenaline, acetylcholine, glutamate, serotonin and histamine. In addition, we show that RA differentiated SH-SY5Y cells express the dopamine and noradrenalin neurotransmitter transporters that are responsible for uptake of MPP(+), a well known DAergic cell toxicant. MPP(+) treatment alters mitochondrial activity according to its proposed cytotoxic effect in DAergic neurons. Taken together, RA differentiated SH-SY5Y cells have a DAergic-like phenotype, and provide a good cellular screening tool to find novel genes or compounds that affect cytotoxic processes that are associated with PD.

Upregulation of a small vault RNA (svtRNA2-1a) is an early event in Parkinson disease and induces neuronal dysfunction

MicroRNAs (miRNAs) and other small non-coding RNAs (sncRNAs) are post-transcriptional regulators of gene expression, playing key roles in neuronal development, plasticity, and disease. Transcriptome deregulation caused by miRNA dysfunction has been associated to neurodegenerative diseases. Parkinson disease (PD) is the second most common neurodegenerative disease showing deregulation of the coding and small non-coding transcriptome. On profiling sncRNA in PD brain areas differently affected, we found that upregulation of a small vault RNA (svtRNA2-1a) is widespread in PD brains, occurring early in the course of the disease (at pre-motor stages). SvtRNA2-1a biogenesis was dependent on Dicer activity on its precursor (vtRNA2-1) but independent of Drosha endonuclease, unlike the canonical miRNAs. Although endogenous svtRNA2-1a was enriched in Ago-2 immunoprecipitates in differentiated SH-SY5Y neuronal cells, overexpression of svtRNA2-1a induced subtle transcriptomic changes, suggesting that gene expression regulation may involve other mechanisms than mRNA decay only. Function enrichment analysis of the genes deregulated by svtRNA2-1a overexpression or svtRNA2-1a predicted targets identified pathways related to nervous system development and cell type specification. The expression pattern of svtRNA2-1a during development and aging of the human brain and the detrimental consequences of a svtRNA2-1a mimic overexpression in neuronal cells further indicate that low svtRNA2-1a levels may be important for the maintenance of neurons. Our results suggest that early svtRNA2-1a upregulation in PD may contribute to perturbations of gene expression networks, underlying metabolic impairment and cell dysfunction. A better understanding of the pathways regulated by svtRNA2-a, and also the mechanisms regulating its expression should facilitate the identification of new targets for therapeutic approaches in PD.

Metal dyshomeostasis and inflammation in Alzheimer's and Parkinson's diseases: possible impact of environmental exposures

A dysregulated metal homeostasis is associated with both Alzheimer's (AD) and Parkinson's (PD) diseases; AD patients have decreased cortex and elevated serum copper levels along with extracellular amyloid-beta plaques containing copper, iron, and zinc. For AD, a putative hepcidin-mediated lowering of cortex copper mechanism is suggested. An age-related mild chronic inflammation and/or elevated intracellular iron can trigger hepcidin production followed by its binding to ferroportin which is the only neuronal iron exporter, thereby subjecting it to lysosomal degradation. Subsequently raised neuronal iron levels can induce translation of the ferroportin assisting and copper binding amyloid precursor protein (APP); constitutive APP transmembrane passage lowers the copper pool which is important for many enzymes. Using in silico gene expression analyses, we here show significantly decreased expression of copper-dependent enzymes in AD brain and metallothioneins were upregulated in both diseases. Although few AD exposure risk factors are known, AD-related tauopathies can result from cyanobacterial microcystin and β-methylamino-L-alanine (BMAA) intake. Several environmental exposures may represent risk factors for PD; for this disease neurodegeneration is likely to involve mitochondrial dysfunction, microglial activation, and neuroinflammation. Administration of metal chelators and anti-inflammatory agents could affect disease outcomes.

Alternative α-synuclein transcript usage as a convergent mechanism in Parkinson's disease pathology

α-Synuclein is implicated both in physiological functions at neuronal synaptic terminals as well as in pathological processes in the context of Parkinson's disease. However, the molecular mechanisms for these apparently diverse roles are unclear. Here we show that specific RNA transcript isoforms of α-synuclein with an extended 3' untranslated region, termed aSynL, appear selectively linked to pathological processes, relative to shorter α-synuclein transcripts. Common variants in the aSynL 3' untranslated region associated with Parkinson's disease risk promote the accumulation and translation of aSynL transcripts. The presence of intracellular dopamine can further enhance the relative abundance of aSynL transcripts through alternative polyadenylation site selection. We demonstrate that the presence of the extended aSynL transcript 3' untranslated region impacts accumulation of α-synuclein protein, which appears redirected away from synaptic terminals and towards mitochondria, reminiscent of Parkinson's disease pathology. Taken together, these findings identify a novel mechanism for aSyn regulation in the context of Parkinson's disease-associated genetic and environmental variations.

Knowledge Driven Variable Selection (KDVS) - a new approach to enrichment analysis of gene signatures obtained from high-throughput data

Background

High–throughput (HT) technologies provide huge amount of gene expression data that can be used to identify biomarkers useful in the clinical practice. The most frequently used approaches first select a set of genes (i.e. gene signature) able to characterize differences between two or more phenotypical conditions, and then provide a functional assessment of the selected genes with an a posteriori enrichment analysis, based on biological knowledge. However, this approach comes with some drawbacks. First, gene selection procedure often requires tunable parameters that affect the outcome, typically producing many false hits. Second, a posteriori enrichment analysis is based on mapping between biological concepts and gene expression measurements, which is hard to compute because of constant changes in biological knowledge and genome analysis. Third, such mapping is typically used in the assessment of the coverage of gene signature by biological concepts, that is either score–based or requires tunable parameters as well, limiting its power.

Results

We present Knowledge Driven Variable Selection (KDVS), a framework that uses a priori biological knowledge in HT data analysis. The expression data matrix is transformed, according to prior knowledge, into smaller matrices, easier to analyze and to interpret from both computational and biological viewpoints. Therefore KDVS, unlike most approaches, does not exclude a priori any function or process potentially relevant for the biological question under investigation. Differently from the standard approach where gene selection and functional assessment are applied independently, KDVS embeds these two steps into a unified statistical framework, decreasing the variability derived from the threshold–dependent selection, the mapping to the biological concepts, and the signature coverage. We present three case studies to assess the usefulness of the method.

Conclusions

We showed that KDVS not only enables the selection of known biological functionalities with accuracy, but also identification of new ones. An efficient implementation of KDVS was devised to obtain results in a fast and robust way. Computing time is drastically reduced by the effective use of distributed resources. Finally, integrated visualization techniques immediately increase the interpretability of results. Overall, KDVS approach can be considered as a viable alternative to enrichment–based approaches.

Neurochemistry and the non-motor aspects of PD

Parkinson disease (PD) is a systemic disease with variegated non-motor deficits and neurological symptoms, including impaired olfaction, autonomic failure, cognitive impairment and psychiatric symptoms, in addition to the classical motor symptoms. Many non-motor symptoms appear before or in parallel with motor deficits and then worsen with disease progression. Although there is a relationship, albeit not causal, between motor symptoms and the presence of Lewy bodies (LBs) and neurites filled with abnormal α-synuclein, other neurological alterations are independent of the amount of α-synuclein inclusions in neurons and neurites, thereby indicating that different mechanisms probably converge in the degenerative process. This may apply to complex alterations interfering with olfactory and autonomic nervous systemfunctions, emotions, sleep regulation, and behavioral, cognitive and mental performance. Involvement of the cerebral cortex leading to impaired behavior and cognition is related to several convergent altered factors including: a. dopaminergic, noradrenergic, serotoninergic and cholinergic cortical innervation; b. synapses; c. cortical metabolism; d. mitochondrial function and energy production; e. oxidative damage; f. transcription; g. protein expression; h. lipid composition; and i. ubiquitin–proteasome system and autophagy, among others. This complex situation indicates that multiple subcellular failure in selected cell populations is difficult to reconcilewith a reductionistic scenario of a single causative cascade of events leading to non-motor symptoms in PD. Furthermore, these alterationsmay appear at early stages of the disease and may precede the appearance of substantial irreversible cell loss by years. These observations have important implications in the design of therapeutic approaches geared to prevention and treatment of PD.

Biosignatures for Parkinson's disease and atypical parkinsonian disorders patients

Diagnosis of Parkinson' disease (PD) carries a high misdiagnosis rate due to failure to recognize atypical parkinsonian disorders (APD). Usually by the time of diagnosis greater than 60% of the neurons in the substantia nigra are dead. Therefore, early detection would be beneficial so that therapeutic intervention may be initiated early in the disease process. We used splice variant-specific microarrays to identify mRNAs whose expression is altered in peripheral blood of early-stage PD patients compared to healthy and neurodegenerative disease controls. Quantitative polymerase chain reaction assays were used to validate splice variant transcripts in independent sample sets. Here we report a PD signature used to classify blinded samples with 90% sensitivity and 94% specificity and an APD signature that resulted in a diagnosis with 95% sensitivity and 94% specificity. This study provides the first discriminant functions with coherent diagnostic signatures for PD and APD. Analysis of the PD biomarkers identified a regulatory network with nodes centered on the transcription factors HNF4A and TNF, which have been implicated in insulin regulation.

Gene expression profiling in human neurodegenerative disease

Transcriptome study in neurodegenerative disease has advanced considerably in the past 5 years. Increasing scientific rigour and improved analytical tools have led to more-reproducible data. Many transcriptome analysis platforms assay the expression of the entire genome, enabling a complete biological context to be captured. Gene expression profiling (GEP) is, therefore, uniquely placed to discover pathways of disease pathogenesis, potential therapeutic targets, and biomarkers. This Review summarizes microarray human GEP studies in the common neurodegenerative diseases amyotrophic lateral sclerosis (ALS), Parkinson disease (PD) and Alzheimer disease (AD). Several interesting reports have compared pathological gene expression in different patient groups, disease stages and anatomical areas. In all three diseases, GEP has revealed dysregulation of genes related to neuroinflammation. In ALS and PD, gene expression related to RNA splicing and protein turnover is disrupted, and several studies in ALS support involvement of the cytoskeleton. GEP studies have implicated the ubiquitin-proteasome system in PD pathogenesis, and have provided evidence of mitochondrial dysfunction in PD and AD. Lastly, in AD, a possible role for dysregulation of intracellular signalling pathways, including calcium signalling, has been highlighted. This Review also provides a discussion of methodological considerations in microarray sample preparation and data analysis.