The Longitudinal Transcriptomic Response of the Substantia Nigra to Intrastriatal 6-Hydroxydopamine Reveals Significant Upregulation of Regeneration-Associated Genes

Nicholas M. Kanaan; Timothy J. Collier; Allyson Cole-Strauss; Tessa Grabinski; Zachary R. Mattingly; Mary E. Winn; Kathy Steece-Collier; Caryl E. Sortwell; Fredric P. Manfredsson; Jack W. Lipton

doi:10.1371/journal.pone.0127768

Abstract

We hypothesized that the study of gene expression at 1, 2, 4, 6 and 16 weeks in the substantia nigra (SN) after intrastriatal 6-OHDA in the Sprague-Dawley rat (rattus norvegicus) would identify cellular responses during the degenerative process that could be axoprotective. Specifically, we hypothesized that genes expressed within the SN that followed a profile of being highly upregulated early after the lesion (during active axonal degeneration) and then progressively declined to baseline over 16 weeks as DA neurons died are indicative of potential protective responses to the striatal 6-OHDA insult. Utilizing a κ-means cluster analysis strategy, we demonstrated that one such cluster followed this hypothesized expression pattern over time, and that this cluster contained several interrelated transcripts that are classified as regeneration-associated genes (RAGs) including Atf3, Sprr1a, Ecel1, Gadd45a, Gpnmb, Sox11, Mmp19, Srgap1, Rab15,Lifr, Trib3, Tgfb1, and Sema3c. All exemplar transcripts tested from this cluster (Sprr1a, Ecel1, Gadd45a, Atf3 and Sox11) were validated by qPCR and a smaller subset (Sprr1a, Gadd45a and Sox11) were shown to be exclusively localized to SN DA neurons using a dual label approach with RNAScope in situ hybridization and immunohistochemistry. Upregulation of RAGs is typically associated with the response to axonal injury in the peripheral nerves and was not previously reported as part of the axodegenerative process for DA neurons of the SN. Interestingly, as part of this cluster, other transcripts were identified based on their expression pattern but without a RAG provenance in the literature. These "RAG-like" transcripts need further characterization to determine if they possess similar functions to or interact with known RAG transcripts. Ultimately, it is hoped that some of the newly identified axodegeneration-reactive transcripts could be exploited as axoprotective therapies in PD and other neurodegenerative diseases.

Citation: Kanaan NM, Collier TJ, Cole-Strauss A, Grabinski T, Mattingly ZR, Winn ME, et al. (2015) The Longitudinal Transcriptomic Response of the Substantia Nigra to Intrastriatal 6-Hydroxydopamine Reveals Significant Upregulation of Regeneration-Associated Genes. PLoS ONE 10(5): e0127768. https://doi.org/10.1371/journal.pone.0127768

Academic Editor: Richard Jay Smeyne, St. Jude Children's Research Hospital, UNITED STATES

Received: February 26, 2015; Accepted: April 20, 2015; Published: May 20, 2015

Copyright: © 2015 Kanaan et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited

Data Availability: All CEL files are available from the NCBI Gene Expression Omnibus Repository (http://www.ncbi.nlm.nih.gov/geo [#GSE58710]).

Funding: This study was supported by the National Institute of Neurological Disorders and Stroke (US), P50 NS058830 (to TJC, JWL, KSC, CES), the Office of the Vice President of Research Graduate Studies, Michigan State University (to JWL), and Saint Mary's Foundation (Doran Foundation), Mercy Health Saint Mary's, Grand Rapids, MI (to NMK). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

Introduction

The 6-hydroxydopamine (6-OHDA) model of Parkinson’s disease (PD) has been extensively characterized and employed to study the consequences of dopamine neuron loss and to screen therapeutics [1–9]. While there is continued debate as to the etiology of PD, it is clear that 6-OHDA rodent models remain relevant preclinical tools for examining the behavioral and biologic consequences of nigrostriatal dopamine (DA) neuron degeneration during disease progression. In this regard, our lesioning paradigm using intrastriatal administration of 6-OHDA [6–8,10] causes retrograde axonal degeneration that results in a rapid and near complete denervation of the striatum (STR) with a slower, more progressive, 76–80% loss of substantia nigra (SN) DA neurons by approximately six weeks [6,8]. Notably, a recent postmortem evaluation of nigrostriatal pathology in people who died in varying stages of PD found that near complete striatal DA axon denervation occurs within the first four years after clinical diagnosis, followed later by progressive DA neuron loss [11]. The key finding of this study was that striatal DA denervation is an early and fairly rapid process in PD. Hence, the intrastriatal 6-OHDA lesion model replicates the basic ordinal progression of the PD degenerative process, albeit in an accelerated timeframe.

Using this model, we hypothesized that the study of progressive gene expression changes over time in the SN after 6-OHDA insult would identify protective (albeit insufficient) cellular responses that are engaged during the neurodegenerative process. Subsequent studies could then assess whether artificial induction or silencing of the observed gene expression changes could be exploited to protect against DA neuronal or axonal damage. In this way, we are not attempting to identify etiologic factors of PD, but rather endeavor to exploit the natural response of DA neurons to degenerative cues in order to identify new manipulatable targets to repair, stabilize, or interfere with the neurodegenerative process of SN DA neurons.

The current study examined the progressive gene expression changes over 16 weeks in the SN after an intrastriatal 6-OHDA lesion. We hypothesized that patterns of gene expression within the SN characterized by early, highly upregulated expression would be indicative of potential protective responses to a striatal 6-OHDA insult. We additionally hypothesized that the use of k-means cluster analysis to group transcripts based upon expression patterns would aid in the identification of previously unknown neuroprotective candidates. Our results demonstrated that intrastriatal 6-OHDA treatment upregulated several interrelated transcripts classified as members of the regeneration-associated gene (RAGs) family that segregated into one cluster. Upregulation of RAGs is historically associated with the response to axonal injury in the peripheral nerves [12–15] and was not previously associated with axonal degeneration of SN DA neurons. In addition, other transcripts displaying a similar regulation pattern were identified and should be examined for similar function or interaction with the RAG family. These clustered RAG and RAG-like transcripts should be examined for their role in neurodegenerative and neuroprotective processes in SN DA neurons.

Materials and Methods

Subjects

Male, Sprague-Dawley rats (rattus norvegicus, Harlan, Indianapolis, IN) weighing between 200–225 g (approximately 60–75 days old) were given access to food and water ad libitum, and housed in an AAALAC accredited animal facility on a standard light-dark cycle. Housing conditions and animal procedures were performed according to the Association for Assessment and Accreditation of Laboratory Animal Care International guidelines and with the approval of the Michigan State University Institutional Animal Care and Use Committee. This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. All surgery was performed under sodium pentobarbital anesthesia, and all efforts were made to minimize suffering. Animals were divided into the following experimental conditions: No surgery control (WK 0), unilateral intrastriatal 6-hydroxydopamine lesion (6-OHDA) or sham surgical controls (SHAM). SHAM and 6-OHDA subjects were sacrificed at 5 different times post-lesion: 1 (WK1), 2 (WK2), 4 (WK4), 6 (WK6) or 16 (WK16) weeks post-surgery. Three subjects from each of the 11 groups (5 time points x 2 treatments + the 1 NAÏVE group) who met all inclusion criteria (see below) were subjected to microarray analyses. The quantitative real-time PCR (qPCR) validation was conducted in the same cohort that was subjected to microarray analyses. Subjects utilized for immunohistochemistry and in situ hybridization (ISH) were sacrificed at 1 week post-lesion.

Unilateral intrastriatal 6-OHDA lesion

Rats were anesthetized with equithesin (0.3ml/100 g of weight, i.p; chloral hydrate 42.5 mg/ml + sodium pentobarbital 9.72 mg/ml) and placed into a stereotaxic frame. Each rat received two injections of 6-OHDA HBr [Sigma-Aldrich, St. Louis, MO; 5 μg/μl, calculated as free base in 0.2% ascorbic acid and a 0.9% saline solution] at 2 μl per site. The two injection coordinates were (i) AP +1.6 mm, ML +2.4 mm, DV −4.2mm, (ii) MP +0.2 mm, ML +2.6 mm, DV −7.0mm. For all rats, the needle was zeroed at the skull directly above the injection site in order to target the DV coordinate. For each injection, the needle was lowered slowly to the injection site and 1 minute elapsed before injection commenced, 6-OHDA was injected at 0.5 μl/min and at the end of the injection the needle remained in place for an additional 4 minutes before retraction.

Tissue collection

At the appropriate post-surgical time point, rats were anesthetized with pentobarbital (50mg/kg intraperitoneally), and decapitated. Brains were removed rapidly and submerged for 30 s in a 250ml beaker of isopentane chilled in powdered dry ice. The brains were then wrapped in foil and stored at -80°C until dissected. Frozen brains were slabbed on an inverted petri dish over a bed of crushed ice. Slabs containing the striatum were dissected using scalpels. A small portion (~ 2mm³) of the striatum from each side of the brain was separately reserved for confirmation of lesion status using HPLC (see below). Tissue from the substantia nigra (SN) was placed in 1 ml of trizol, (Life Technologies, Carlsbad, CA), homogenized by hand with a disposable plastic pestle, frozen on dry ice and stored at -80°C in preparation for RNA isolation

High performance liquid chromatography (HPLC)

Striatal DA levels (i.e. 6-OHDA lesion status) were quantified by HPLC as described previously[16–18]. Briefly, samples were sonicated into 150 μl (SN) or 250 μl (STR) of a 0.4 N perchlorate, 1.34 mM EDTA and 0.53 mM sodium metabisulfite solution. A 20 μl aliquot of the homogenate was reserved for protein determination and the remaining homogenate was centrifuged at 10,500 rpm for 10 minutes at 4°C. The supernatant was stored in a separate tube at-80°C. Sample separation was performed on a 250x4.6mm Microsorb MV C18 100–5 column (Agilent, Santa Clara, CA). DA levels were detected and quantitated using a 12-channel CoulArray 5200 coulometric array detector (ESA, Chelmsford, MA). The mobile phase consisted of 100mM Citric Acid, 75mM Na₂HPO₄Na, 80°M 1-heptanesulfonate monohydrate, sodium salt, 5% MeOH, pH 4.25. Samples values were interpolated against a 6 point standard curve. The final values were standardized based on protein content (BCA Protein Assay Kit, Pierce Inc., Rockford, IL). Striatal DA depletion of >95% in the lesioned hemisphere as compared to the unlesioned hemisphere was used as a criterion for inclusion in the study.

RNA isolation and quality evaluation

RNA extraction was performed using the RNA Clean and Concentrator kit (Zymo Research, Irvine, CA) and eluted into 15μl H₂O. RNA quality was assessed using the RNA Nano 6000 Assay on an Agilent Bioanalyzer, (Santa Clara, CA). RNA quality was measured using the 10-point scale associated with the RNA Integrity Number (RIN). Only samples with RIN values ≥ 7 qualified for inclusion in microarray analyses. The mean and standard deviation of RIN values for all of the samples was 8.7 ± 0.71 (n = 33).

Microarray sample processing and hybridization

Isolated RNA from tissue samples (n = 33) were processed for microarray hybridization on the Rat Gene 1.0 ST Array at the Gene Expression Microarray Core of Cincinnati Children’s Hospital Medical Center, Cincinnati, OH. 50-120ng of total RNA was converted to biotin–labeled sense-strand cDNA for hybridization using the Ambion WT Expression Kit (Life Technologies, Carlsbad, CA) combined with the GeneChip WT Terminal Labeling Kit (Affymetrix, Santa Clara, CA). Chips were incubated at 45°C for 17 hours in the GeneChip Hybridization Oven 640, washed and stained in the Fluidics Station 450 (Affymetrix, Santa Clara, CA), and scanned using an Affymetrix Gene Chip Scanner 3000 7G (Affymetrix, Santa Clara, CA).

Microarray image analysis and quality control

Only array images meeting all of the quality control measures defined by the Affymetrix Expression Control Program were included in this study. Specific quality control metrics included signal histogram, relative log expression signal, Pearson’s correlation, PM mean (average signal intensity of probes), and positive and negative area under the curve (AUC). Also measured were the expression values of spiked-in poly-A RNA controls, and values of spiked-in hybridization controls. Raw data are available on the NCBI Gene Expression Omnibus Repository (http://www.ncbi.nlm.nih.gov/geo [#GSE58710])

Analysis of microarray data

Data processing and analysis were performed using R v3.0.2 (http://www.R-project.org) and Bioconductor [19]. Affymetrix-formatted CEL files were imported into R, assessed for quality, and Robust-Multi Array (RMA) normalized at the transcript level using the oligo package (v1.26.0). Array quality was evaluated via visual assessment of the array images, intensity distributions, multi-array plots, and principal component analysis. Differential expression between 6-OHDA and sham rats was assessed at each time point using limma (v3.18.6). The resulting moderated t-statistics with an associated Benjamini-Hochberg false discovery rate (FDR) adjusted p-value less than 0.05 were used to identify expression patterns of interest using ConsensusClusterPlus (v1.16.0) [20]. Consensus clustering is a quantitative, unsupervised method used to assign objects into groups (i.e. clusters). Clustering is performed using the intrinsic relationships within the data and without external information. Repeated clustering using subsets of the data provides information regarding confidence in the number of groups and the genes within each group. The moderated t-statistic is the ratio of the log2-fold change to its standard error. T-statistics were utilized due to improved stability of consensus clustering as compared to log2-fold change. Cluster number (k) was chosen for cluster stability and a maximum mean cluster consensus score after an initial drop.

Validation of microarray gene expression changes by qPCR

To validate expression differences observed from microarray analysis, qPCR was performed to quantify levels of Sprr1a, Ecel1, Gadd45, Atf3, and Sox11. The mRNA levels in frozen SN micropunches were derived from the same subjects at the WK2 post-surgical time point. Total RNA was converted to cDNA using Superscript VILO Mastermix (Invitrogen/Life Technologies, Carlsbad, CA). PCR reactions were run in 30 μl using target specific, FAM labeled Taqman hydrolysis probes (Applied Biosystems/Life Technologies, Carlsbad, CA) multiplexed with a VIC labeled primer-limited probe set to the Gapdh reference gene, which was confirmed to not vary with respect to the different treatment conditions (data not shown). Normalized gene expression was determined by differences in the cycle thresholds (Ct) between genes of interest and Gapdh (∆Ct) on a Viia7 qPCR System (Life Technologies). Differences between these normalized values across treatment groups (∆∆Ct) were then linearized to determine fold change differences (2^-∆∆Ct) and compared to fold changes observed from samples subjected to microarray hybridization and quantitation.

Pathway analysis of Cluster-1

The 94 probe sets representing 88 genes from Cluster-1 and their corresponding fold changes from WK1 (see S1 Table) were imported into MetaCore version 6.19 build 65960 (Thomson Reuters, New York, NY) for pathway and network analysis. Pathway analysis was performed using the Pathway Maps One-Click Analysis. MetaCore Pathway Maps are rigorously curated pathways supported by experimental evidence reported in the literature. An algorithm based on the hypergeometric distribution is used to calculate enrichment p-values. Pathway Maps with a FDR adjusted p-value less than 0.05 were considered significant. Network analysis was performed using the “shortest paths” algorithm with a maximum number of steps in the path set to 1 without the use of canonical pathways. Nodes and edges were exported from MetaCore and imported into Cytoscape v3.1.1 for final visualization.

Tyrosine hydroxylase (TH) immunostaining and in situ hybridization (ISH) of Cluster-1 transcripts to determine colocalization

Animals were sacrificed 1 week after 6-OHDA (the peak of transcript expression in Cluster -1 and prior to significant loss of SN DA neurons). We utilized a combination of immunohistochemical (IHC) staining for TH and ISH for determining whether the highest fold-change Cluster-1 transcript, Sprr1a, and two transcripts known to influence Sprr1a expression Sox11 and Gadd45a were localized to DA neurons in the SN. Briefly, the brains were cut into 20μm thick free-floating sections and sections were first processed for RNAscope ISH for Sprr1a, Sox11 or Gadd45a. RNAscope is a commercially available branched or “tree” hybridization technology provided as a kit by Advanced Cell Diagnostics. Custom RNAscope target probes were generated against Sprr1a (407411, NM_021864.1, probe region 3–801), Sox11 (414971, NM_053349.1, probe region 304–1955) and Gadd45a (414961, NM_024127.2, probe region 5–995). The ISH signal was detected using 3, 3’-diaminobenzidine to produce brown punctate staining. Following RNAscope ISH labeling, the sections were processed for TH (a DA neurons marker, Millipore, MAB318, 1:50,000) IHC using Vector slate-grey (SG) to produce blue-grey labeling of DA neurons (for a detailed protocol see [21]). DA neurons were identified both by their expression of TH and their morphology within the SN pars compacta (SNpc) as TH neurons can lose their phenotype when stressed or undergoing degeneration.

Results

Lesion Status Confirmation

Lesion status in the STR was quantified by the proportion of DA (ng/mg protein) from the ipsilateral (operated) to the contralateral (non-operated) side and converted to a percentage of the contralateral side. The 6-OHDA subjects that met exclusion criteria had mean DA levels of 0.8% ± 2.1% in the operated striatum (0.9 ng/mg protein) as compared to the intact side (119.1 ng/mg protein). Sham subjects had mean DA levels of 101.8% ±3.2% in the operated striatum (126.3 ng/mg protein) as compared to the non-operated striatum (124.1 ng/mg protein).

Consequences of striatal 6-OHDA injections on transcriptomic changes in the SN over 16 weeks

Microarray analysis.

Using k-means cluster analysis we identified eight clusters of genes in the SN that covaried over time as a result of intrastriatal 6-OHDA lesion (Fig 1). For data from all clusters see S2 Table. Transcripts in Cluster-1 exhibited a robust increase one week following the 6-OHDA lesion, progressively declined at 2, 4 and 6 weeks post-lesion, and returned to basal or near basal levels by 16 weeks. The top 20 transcript changes for Cluster-1 are depicted in Fig 2 and enumerated in Table 1; significantly enriched pathways for the Cluster-1 are listed in Table 2. Interestingly, 13 transcripts within Cluster-1 fit the definition of “regeneration-associated genes” or RAGs, which identify a functional class of genes that are up-regulated in neuronal perikarya following axotomy or axonal damage [12–15,22–35]. These genes include Atf3, Sprr1a, Ecel1, Gadd45a, Gpnmb, Sox11, Mmp19, Srgap1, Rab15,Lifr, Trib3, Tgfb1, and Sema3c. The changes in expression of these RAGs in comparison to other members of the cluster are shown in Fig 2. Since the remainder of the transcripts in Cluster-1 demonstrated a similar differential expression profile, henceforth we will refer to them as “RAG-like” transcripts.

Download:

Fig 1. K-means cluster analysis identifies 88 genes acutely upregulated in response to lesioning.

Eight clusters were identified through consensus clustering analysis. We hypothesized that a cluster of genes would be highly upregulated due to 6-OHDA administration that would fall-off over time due to ongoing DA neuron degeneration. Cluster-1 followed this hypothesized pattern.

https://doi.org/10.1371/journal.pone.0127768.g001

Download:

Fig 2. Thirteen known regeneration-associated genes are part of the 88 genes of Cluster-1.

The expression pattern of the non-RAG genes follows a “RAG-like” expression pattern suggesting that they may too be part of the RAG family.

https://doi.org/10.1371/journal.pone.0127768.g002

Download:

Table 1. Top twenty differentially expressed genes from Cluster-1.

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

Download:

Table 2. Significantly enriched pathways in Cluster-1.

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

While the remaining seven clusters did not fit with our a priori hypothesis, these data are available in S2 Table. Transcripts within these other seven clusters showed varying patterns in the degree of change and time course of gene expression. For example, Cluster-2 genes exhibited progressive downregulation from WK1 to WK2 that persisted through WK16. This cluster contained a number of transcripts associated with the DAergic phenotype. However, clusters 2–8 should be treated with caution due to a lack of a priori hypotheses regarding their patterns and the reality that consensus clustering can produce false clusters and inaccurate identification of the true k [36]. Due to the large number of genes being assessed and repeated sampling of the data there is a high probability that many of these same patterns could be found among a completely unrelated set of samples. Regardless, the differences in expression between lesioned and sham rats at each time point remain valid and should be further explored.

qPCR validation of genes from Cluster-1.

The changes in expression levels for Sprr1a, Ecel1, Gadd45, Atf3, and Sox11 were validated by qPCR in samples from the same 6-OHDA-WK1 and SHAM-WK1 subjects evaluated by microarray. These five RAG transcripts were chosen for validation due to the highly enriched presence of RAGs in Cluster-1 and were among the top six FDR adjusted p-values in Cluster-1. The qPCR confirmed a similar fold change in expression levels for, Ecel1, Gadd45a, Atf3, and Sox11. However Sprr1a, which demonstrated a 19 fold increase by microarray, exhibited an 80-fold increase by qPCR, suggesting that the microarray may have reached signal saturation (see Fig 3).

Download:

Fig 3. Exemplar Transcripts from Cluster-1 are validated by qPCR.

The top 5 RAG transcripts from Cluster-1 were subjected to validation by qPCR. There was a high degree of agreement between the microarray and qPCR differential expression with the exception of Sprr1a. Differential expression of Sprr1a in the SN of 6-OHDA lesioned subjects was much greater when quantified by qPCR (80.3 fold) than originally observed by microarray (19.2 fold). This is likely due to saturation of the Sprr1a signal on the microarray chip, resulting in an underestimate of the differential upregulation.

https://doi.org/10.1371/journal.pone.0127768.g003

Demonstration that exemplar Cluster-1 transcripts were localized to DA neurons of the SNpc.

Next, we sought to confirm our hypothesis that the upregulation of genes attempting to confer axoprotection (i.e. RAGs) in Cluster-1 occurred specifically in DA neurons of the SN. We therefore chose three transcripts from Cluster-1 to pursue further using a combination of ISH to identify transcript changes in the RAGs and IHC to identify SN DA neurons. We observed a robust upregulation of Sprr1a, Sox11 and Gadd45a transcripts (as indicated by RNAscope ISH) in TH+ SNpc neurons (Fig 4). These data suggest that, at least for the exemplar transcripts examined here, the progressive DA neuron loss from this intrastriatal 6-OHDA lesion is likely responsible for the progressive loss of signal over time in the transcripts of Cluster-1. These data also directly confirm that these RAGs are upregulated specifically in SN DA neurons in response to axonal degeneration by 6-OHDA.

Download:

Fig 4. Upregulation of Sprr1a, Gadd45a and Sox11 early after intrastriatal 6-OHDA lesion occurs specifically in substantia nigra (SN) dopamine (DA) neurons. A-F).

The mRNA for Sprr1a (A and B), Gadd45a (C and D) and Sox11 (E and F) was identified using RNAscope in situ hybridization (ISH; brown puncta) and DA neurons in the SN were identified with immunohistochemistry for tyrosine hydroxylase (blue staining, a DA neuron marker) in animals one week after an intrastriatal 6-OHDA lesion. The DA neurons on the intact side (i.e. contralateral to the 6-OHDA lesion) had little to no ISH signal due to the normally low basal expression of these RAGs (A, C, and E). In contrast, the DA neurons on the lesioned side had robust ISH signal (B, D and F). In the lesioned side, DA neurons appeared unhealthy and three cellular phenotypes were present, 1) little to no RAG ISH but strong TH indicating an apparently healthy neuron (arrow), 2) mild-moderate RAG ISH with moderate-weak TH staining indicating a DA neuron that is likely enroute to degenerating (arrowhead), and 3) robust ISH signal and little-no TH staining likely indicating DA neurons closer to degenerating (asterisk; DA neuron determination is based on neuroanatomical location within the SN and cell morphology). These data indicate that the upregulation of Sprr1a, Gadd45a and Sox11 identified by gene arrays and confirmed by qPCR occurred specifically in SN DA neurons. Scale bar = 50μm.

https://doi.org/10.1371/journal.pone.0127768.g004

Network analysis of cCluster-1

MetaCore network analysis demonstrated that transcripts from Cluster-1 were highly connected, with 62 of the transcripts showing evidence of interaction (Fig 5). This network included four significantly overconnected nodes (Ddit3; p = 4.14e-06, Atf5; p = 5.12e-05, Atf3; p = 0.001, Cebpg; p = 0.003, Mthdf3; p = 0.002) that are indicative of a response to stress [37,38].

Download:

Fig 5. Pathway analysis demonstrates a high degree of interaction among the transcripts of Cluster-1.

MetaCore Pathway analysis was conducted on the differential gene expression of Cluster-1 at 1 week post-surgery. The data were imported into Cytoscape for visualization. Diamonds represent transcription factors and ovals represent all other proteins. Arrows represent the direction of the interaction. The intensity of red coloring is indicative of the magnitude of upregulation. Values ranged from 1.23 (lightest pink) to 19.2 (darkest red). Of the 88 genes in Cluster-1, 62 had known interactions. The significant level of interaction among the transcripts along with the a priori hypothesis for the existence of Cluster-1 suggests this is a stable, biologically relevant cluster.

https://doi.org/10.1371/journal.pone.0127768.g005

Discussion

In this study, we examined transcriptomic changes in the SN over 16 weeks after a unilateral, intrastriatal 6-OHDA lesion that results in progressive nigrostriatal DA axon degeneration followed by SN DA neuron loss [10]. Thus, this model replicates the early, profound nigrostriatal axon degeneration recently observed in post-mortem PD brains [11]. By using a novel k-means cluster analysis strategy, we generated clusters of transcripts defined by their covariation over time. This clustering suggests that groups of transcripts are responding to disconnection from the striatum in the same way or are interacting with each other in a coordinated fashion to respond to the 6-OHDA insult. Our a priori hypothesis was that genes expressed after induction of axonal degeneration would include potential axoprotective responses that are highly upregulated in the immediate/acute phases after injury but would then fall as overt loss of DA neurons occurs. In support of this hypothesis, we observed that transcripts with axon regenerative properties in the peripheral nervous system [12–15,24,30,39], including several members of the RAG family, are prominently upregulated early after insult in the SN followed by a rapid waning of expression (Cluster-1). The robust upregulation of a number of the RAGs identified by microarray were confirmed using qPCR, and in the case of Sprr1a, an increase of around 80 fold was observed.

The rapid induction and subsequent reduction in RAG and RAG-like expression over time indicates that these transcripts likely are (1) expressed in DA neurons of the SN and (2) their reduction over time is a result of the progressive degeneration of those DA neurons. Indeed, combing TH IHC with ISH for either Sprr1a, Sox11 or Gadd45a confirmed that expression of at least these three RAG transcripts were substantially increased specifically in DA neurons of the SNpc (see Fig 4). While we cannot definitively conclude that the origin of all differentially expressed Cluster-1 transcripts in the SN come from DA neurons, our evidence supports the hypothesis that their rapid up-regulation is most likely the result of the DA axon damage that occurs in the initial stages of STR terminal damage. Future studies are needed to determine the cellular origin of the remaining member transcripts from Cluster-1. Additionally, we need to determine, prior to the one week time point, the temporal order of induction of Cluster-1 members in order to better understand the ordinal response of this gene cluster to axonal insult.

The specific RAG transcripts that we found upregulated in Cluster-1 include Atf3, Sprr1a, Ecel1, Gadd45a, Gpnmb, Sox11, Mmp19, Srgap1, Rab15,Lifr, Trib3, Tgfb1, and Sema3c. By definition, all of these transcripts are associated with a response to axonal damage and exhibit potential axoprotective or regenerative properties. Moreover, the majority of transcripts from Cluster-1 have significant interactions as illustrated in the network pathway analysis depicted in Fig 5, suggesting that other transcript members of the cluster may have RAG-like properties.

Atf3 is among the more well-characterized RAGs in response to axonal damage. Atf3 is upregulated in damaged dorsal root ganglion (DRG) neurons concurrent with neurite elongation in vitro [12]. Similarly, Sprr1a is relatively well-studied in peripheral neurons, and following DRG neuron axon damage it is significantly upregulated. The Sprr1a protein is associated with the axonal growth cone and is postulated to be integral to the growth and repair process of optic nerve and DRG neurons. Overexpression of Sprr1a enhances axonal recovery after injury, while reducing expression of Sprr1a impairs regenerative capabilities [14,15,24]. Sprr1a also known covaries with the transcription factor Atf3 after axonal injury in cultured neurons [14]. Furthermore, neurite outgrowth is enhanced and Sprr1a expression is induced in cultured DRG neurons transfected with Sox11 [15], suggesting that Sox11 is a transcription factor for Sprr1a. In this regard, the present TH/Sox11 and TH/Sprr1a immunostaining analysis demonstrates that both of these RAGs are highly upregulated in damaged DA neurons (Fig 4). This novel observation indicates that the response of SN DA neurons to axonal injury is under the control of many of the same factors observed in DRG neurons. Functionally, Sox11 may regulate expression of cytoskeletal proteins (e.g. βIII-tubulin, MAP2 and doublecortin), while Sprr1a may directly interact with structural cytoskeleton proteins critical (e.g. actin; for a review see [40]).

Another member of the RAG family induced by axonal injury in nigrostriatal DA neurons is Gadd45a. Spinal cord injury induces Gadd45a expression [41,42] in both sensory and motor neurons of the spinal cord. Neuronal damage activates the JNK-p53-GADD45α cascade and Gadd45a inhibition may be a potential neuroprotective strategy [25], whereas another study suggests that Gadd45a may be neuroprotective as its overexpression can prevent neurodegeneration [13]. The underlying cause of these discrepant findings is unclear and whether modulation of Gadd45a expression will affect DA neuron degeneration awaits future studies.

Intriguingly, parkin overexpression induces Gadd45a expression in vitro suggesting a link between some forms of parkinsonism and this transcript [43]. With respect to other RAGs in Cluster-1, Ecel1 (aka, DINE) is upregulated in response to DRG injury [22] and injury to retinal ganglion cells induces Ecel1, Atf3, Gadd45a and Sprr1a [30]. Once again, these data demonstrate the highly interrelated response of RAGs to injury among a wide variety of neuronal cell types.

Other transcripts that behaved similar to established RAGs are less well characterized yet may play a significant role in the response to axonal damage in DA neurons. For instance, Mmp19 expression was only recently found in neurons and is speculated to promote axonal growth and remodeling [31]. Lifr, the alpha receptor for leukemia inhibitory factor, also has axon regenerative properties in retinal ganglion cells [39]. Srgap1 interacts with axonal guidance molecules during development [39] and plays a potentially important role in neurite outgrowth [32]. Likewise, the small GTPase Rab15 is implicated in neuronal differentiation [39]. Tgfb1, has anti-apoptotic, anti-inflammatory and pro-myelin maintenance and production properties [35,44,45]. Gpnmb may show similar functions as it was recently identified as a potential ALS-related factor in the spinal cords of SOD1^G93A mice [26]. Another RAG-like interaction of note is between Trib3 and Atf4. Trib3 is a pseudokinase that is upregulated in response to neuronal death and increases the activity of Atf4 [33,34]. Atf4, in turn, is implicated in the axodegenerative process [46], suggesting that in our model Trib3 upregulation might signal antagonistic axoprotective responses such as RAG induction. Finally, Sema3c is an axonal guidance molecule that is part of the semaphorin family of genes which are involved in axonal navigation and pathfinding in development and neuronal injury [23,27–29,47]. The remainder of Cluster-1 warrants thorough examination for potential roles in the coordinated response of the SNc to DA axon damage, due to their similar RAG-like regulation profile.

Exactly what signals cause the upregulation of RAG and RAG-like transcripts requires further investigation. Interestingly, several transcription factors that are implicated in the ER stress response show significant pathway interaction as indicated by MetaCore analysis (Ddit3, Atf5, Atf3, Cebpg)[37,38]. An important gene connected to these transcription factors is Bcl2l11, a known important mediator of apoptosis [48]. Understanding which of these stress-induced transcripts are signals of cellular trauma, and which are transducers of axoprotective responses will be critical in clarifying their complex, coordinated interrelationship. This understanding could lead to new targets for the development of axoprotective therapeutics.

Notably, the SN RAG transcripts up-regulated in response to striatal 6-OHDA administration were not previously described in DA neurons. For many of these transcripts, this is the first time they were observed to be upregulated in response to CNS injury [12,13,24,30,39,49–53]. This striking observation suggests a conservation of genetic axoprotective programs within the nervous system. Furthermore, these data underscore the importance of these gene pathways as critical mediators of axonal maintenance and regrowth after injury in nigrostriatal DAergic neurons.

Finally, it should be restated here that while there is a wealth of data demonstrating differential gene expression across genes from all eight clusters, one must be cautious in drawing conclusions about the association of genes within any particular cluster without an a priori hypotheses due to a high probability of false positive associations. As a result, new hypotheses could be proposed regarding the longitudinal response to intrastriatal lesioning as long as the appropriate caveats are acknowledged.

Recent discoveries regarding the early and robust nature of axonal degeneration in the DAergic nigrostriatal pathway that precedes overt DA neuron loss in the SN highlights the importance of targeting axonal degeneration in treating PD. In light of this recent observation [11], the data presented here provide potential new avenues to pursue toward target development to prevent or reverse nigrostriatal axon degeneration. The fact that many of the RAG changes are attributable to the DA neurons themselves very early after 6-OHDA insult suggest that transcript specific responses of DA neurons are attempting to resist and repair axonal degeneration. Moving forward, the critical next step is to perform studies aimed at determining whether modulation of the identified DA neuron-derived RAGs and RAG-like transcripts can provide axoprotection and resultant neuroprotection in SN DA neurons. These novel transcript responses identified in DA neurons may result in new strategies for preventing the early axonal loss observed in PD patients [11].

Supporting Information

S1 Table. Fold changes at each time point, F-statistics and p values for the transcripts within Cluster-1.

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

(TXT)

S2 Table. Cluster assignment (1–8) and t-statistic at each time point for differentially expressed transcripts from the microarray analysis.

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

(TXT)

Acknowledgments

We would like to thank Dr. Scott E. Counts for his editorial advice on the manuscript and the Mattingly Family for the distinct privilege of allowing us the opportunity to work with Mr. Zachary R. Mattingly.

Author Contributions

Conceived and designed the experiments: NMK TJC KSC CES FPM JWL. Performed the experiments: NMK TJC ACS TG ZRM KSC CES FPM JWL. Analyzed the data: NMK MEW ACS JWL. Contributed reagents/materials/analysis tools: NMK MEW JWL. Wrote the paper: NMK MEW JWL.

References

1. Mendez JS, Finn BW. Use of 6-hydroxydopamine to create lesions in catecholamine neurons in rats. J Neurosurg. 1975 Feb;42(2):166–73. pmid:234524
- View Article
- PubMed/NCBI
- Google Scholar
2. Henry B, Crossman AR, Brotchie JM. Characterization of enhanced behavioral responses to L-DOPA following repeated administration in the 6-hydroxydopamine-lesioned rat model of Parkinson's disease. Exp Neurol. 1998 Jun;151(2):334–42. pmid:9628768
- View Article
- PubMed/NCBI
- Google Scholar
3. Vu TQ, Ling ZD, Ma SY, Robie HC, Tong CW, Chen EY, et al. Pramipexole attenuates the dopaminergic cell loss induced by intraventricular 6-hydroxydopamine. J Neural Transm. 2000;107(2):159–76. pmid:10847557
- View Article
- PubMed/NCBI
- Google Scholar
4. Ojeda V, Fuentealba JA, Galleguillos D, Andres ME. Rapid increase of Nurr1 expression in the substantia nigra after 6-hydroxydopamine lesion in the striatum of the rat. J Neurosci Res. 2003 Sep 1;73(5):686–97. pmid:12929136
- View Article
- PubMed/NCBI
- Google Scholar
5. Ling ZD, Chang Q, Lipton JW, Tong CW, Landers TM, Carvey PM. Combined toxicity of prenatal bacterial endotoxin exposure and postnatal 6-hydroxydopamine in the adult rat midbrain. Neuroscience. 2004;124(3):619–28. pmid:14980732
- View Article
- PubMed/NCBI
- Google Scholar
6. Spieles-Engemann AL, Steece-Collier K, Behbehani MM, Collier TJ, Wohlgenant SL, Kemp CJ, et al. Subthalamic nucleus stimulation increases brain derived neurotrophic factor in the nigrostriatal system and primary motor cortex. J Parkinsons Dis. 2011;1(1):123–36. pmid:22328911
- View Article
- PubMed/NCBI
- Google Scholar
7. Gombash SE, Lipton JW, Collier TJ, Madhavan L, Steece-Collier K, Cole-Strauss A, et al. Striatal pleiotrophin overexpression provides functional and morphological neuroprotection in the 6-hydroxydopamine model. Mol Ther. 2012 Mar;20(3):544–54. pmid:22008908
- View Article
- PubMed/NCBI
- Google Scholar
8. Gombash SE, Manfredsson FP, Mandel RJ, Collier TJ, Fischer DL, Kemp CJ, et al. Neuroprotective potential of pleiotrophin overexpression in the striatonigral pathway compared with overexpression in both the striatonigral and nigrostriatal pathways. Gene Ther. 2014 May 8.
9. Zhang Y, Meredith GE, Mendoza-Elias N, Rademacher DJ, Tseng KY, Steece-Collier K. Aberrant restoration of spines and their synapses in L-DOPA-induced dyskinesia: involvement of corticostriatal but not thalamostriatal synapses. J Neurosci. 2013 Jul 10;33(28):11655–67. pmid:23843533
- View Article
- PubMed/NCBI
- Google Scholar
10. Spieles-Engemann AL, Behbehani MM, Collier TJ, Wohlgenant SL, Steece-Collier K, Paumier K, et al. Stimulation of the rat subthalamic nucleus is neuroprotective following significant nigral dopamine neuron loss. Neurobiol Dis. 2010 Jul;39(1):105–15. pmid:20307668
- View Article
- PubMed/NCBI
- Google Scholar
11. Kordower JH, Olanow CW, Dodiya HB, Chu Y, Beach TG, Adler CH, et al. Disease duration and the integrity of the nigrostriatal system in Parkinson's disease. Brain. 2013 Aug;136(Pt 8):2419–31. pmid:23884810
- View Article
- PubMed/NCBI
- Google Scholar
12. Seijffers R, Mills CD, Woolf CJ. ATF3 increases the intrinsic growth state of DRG neurons to enhance peripheral nerve regeneration. J Neurosci. 2007 Jul 25;27(30):7911–20. pmid:17652582
- View Article
- PubMed/NCBI
- Google Scholar
13. Lin CR, Yang CH, Huang CE, Wu CH, Chen YS, Sheen-Chen SM, et al. GADD45A protects against cell death in dorsal root ganglion neurons following peripheral nerve injury. J Neurosci Res. 2011 May;89(5):689–99. pmid:21337369
- View Article
- PubMed/NCBI
- Google Scholar
14. Starkey ML, Davies M, Yip PK, Carter LM, Wong DJ, McMahon SB, et al. Expression of the regeneration-associated protein SPRR1A in primary sensory neurons and spinal cord of the adult mouse following peripheral and central injury. J Comp Neurol. 2009 Mar 1;513(1):51–68. pmid:19107756
- View Article
- PubMed/NCBI
- Google Scholar
15. Jing X, Wang T, Huang S, Glorioso JC, Albers KM. The transcription factor Sox11 promotes nerve regeneration through activation of the regeneration-associated gene Sprr1a. Exp Neurol. 2012 Jan;233(1):221–32. pmid:22024412
- View Article
- PubMed/NCBI
- Google Scholar
16. Perez SE, Lazarov O, Koprich JB, Chen EY, Rodriguez-Menendez V, Lipton JW, et al. Nigrostriatal dysfunction in familial Alzheimer's disease-linked APPswe/PS1DeltaE9 transgenic mice. J Neurosci. 2005 Nov 2;25(44):10220–9. pmid:16267229
- View Article
- PubMed/NCBI
- Google Scholar
17. Koprich JB, Chen EY, Kanaan NM, Campbell NG, Kordower JH, Lipton JW. Prenatal 3,4-methylenedioxymethamphetamine (ecstasy) alters exploratory behavior, reduces monoamine metabolism, and increases forebrain tyrosine hydroxylase fiber density of juvenile rats. Neurotoxicol Teratol. 2003 Sep;25(5):509–17. pmid:12972064
- View Article
- PubMed/NCBI
- Google Scholar
18. Davis JF, Tracy AL, Schurdak JD, Tschop MH, Lipton JW, Clegg DJ, et al. Exposure to elevated levels of dietary fat attenuates psychostimulant reward and mesolimbic dopamine turnover in the rat. Behav Neurosci. 2008 Dec;122(6):1257–63. pmid:19045945
- View Article
- PubMed/NCBI
- Google Scholar
19. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80. pmid:15461798
- View Article
- PubMed/NCBI
- Google Scholar
20. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010 Jun 15;26(12):1572–3. pmid:20427518
- View Article
- PubMed/NCBI
- Google Scholar
21. Grabinski T, Kneynsberg A, Manfredsson FP, Kanaan NM. A method for combining RNAscope in situ hybridization with immunohistochemistry in thick free-floating brain sections and primary neuronal cultures. PLoS One. 2015;in press.
22. Kiryu-Seo S. Identification and functional analysis of damage-induced neuronal endopeptidase (DINE), a nerve injury associated molecule. Anat Sci Int. 2006 Mar;81(1):1–6. pmid:16526590
- View Article
- PubMed/NCBI
- Google Scholar
23. Gonthier B, Nasarre C, Roth L, Perraut M, Thomasset N, Roussel G, et al. Functional interaction between matrix metalloproteinase-3 and semaphorin-3C during cortical axonal growth and guidance. Cereb Cortex. 2007 Jul;17(7):1712–21. pmid:17021275
- View Article
- PubMed/NCBI
- Google Scholar
24. Marklund N, Fulp CT, Shimizu S, Puri R, McMillan A, Strittmatter SM, et al. Selective temporal and regional alterations of Nogo-A and small proline-rich repeat protein 1A (SPRR1A) but not Nogo-66 receptor (NgR) occur following traumatic brain injury in the rat. Exp Neurol. 2006 Jan;197(1):70–83. pmid:16321384
- View Article
- PubMed/NCBI
- Google Scholar
25. Choi HJ, Kang KS, Fukui M, Zhu BT. Critical role of the JNK-p53-GADD45alpha apoptotic cascade in mediating oxidative cytotoxicity in hippocampal neurons. Br J Pharmacol. 2011 Jan;162(1):175–92. pmid:20955365
- View Article
- PubMed/NCBI
- Google Scholar
26. Tanaka H, Shimazawa M, Kimura M, Takata M, Tsuruma K, Yamada M, et al. The potential of GPNMB as novel neuroprotective factor in amyotrophic lateral sclerosis. Sci Rep. 2012;2:573. pmid:22891158
- View Article
- PubMed/NCBI
- Google Scholar
27. Steup A, Lohrum M, Hamscho N, Savaskan NE, Ninnemann O, Nitsch R, et al. Sema3C and netrin-1 differentially affect axon growth in the hippocampal formation. Mol Cell Neurosci. 2000 Feb;15(2):141–55. pmid:10673323
- View Article
- PubMed/NCBI
- Google Scholar
28. Moreno-Flores MT, Martin-Aparicio E, Martin-Bermejo MJ, Agudo M, McMahon S, Avila J, et al. Semaphorin 3C preserves survival and induces neuritogenesis of cerebellar granule neurons in culture. J Neurochem. 2003 Nov;87(4):879–90. pmid:14622119
- View Article
- PubMed/NCBI
- Google Scholar
29. Ara J, Bannerman P, Hahn A, Ramirez S, Pleasure D. Modulation of sciatic nerve expression of class 3 semaphorins by nerve injury. Neurochem Res. 2004 Jun;29(6):1153–9. pmid:15176472
- View Article
- PubMed/NCBI
- Google Scholar
30. Fernandes KA, Harder JM, Kim J, Libby RT. JUN regulates early transcriptional responses to axonal injury in retinal ganglion cells. Exp Eye Res. 2013 Jul;112:106–17. pmid:23648575
- View Article
- PubMed/NCBI
- Google Scholar
31. Fudge NJ, Mearow KM. Extracellular matrix-associated gene expression in adult sensory neuron populations cultured on a laminin substrate. BMC Neurosci. 2013;14:15. pmid:23360524
- View Article
- PubMed/NCBI
- Google Scholar
32. Ma Y, Mi YJ, Dai YK, Fu HL, Cui DX, Jin WL. The inverse F-BAR domain protein srGAP2 acts through srGAP3 to modulate neuronal differentiation and neurite outgrowth of mouse neuroblastoma cells. PLoS One. 2013;8(3):e57865. pmid:23505444
- View Article
- PubMed/NCBI
- Google Scholar
33. Ohoka N, Yoshii S, Hattori T, Onozaki K, Hayashi H. TRB3, a novel ER stress-inducible gene, is induced via ATF4-CHOP pathway and is involved in cell death. EMBO J. 2005 Mar 23;24(6):1243–55. pmid:15775988
- View Article
- PubMed/NCBI
- Google Scholar
34. Ord D, Ord T. Characterization of human NIPK (TRB3, SKIP3) gene activation in stressful conditions. Biochem Biophys Res Commun. 2005 Apr 29;330(1):210–8. pmid:15781252
- View Article
- PubMed/NCBI
- Google Scholar
35. Nie X, Tian W, Zhang Y, Chen X, Dong R, Jiang M, et al. Induction of transforming growth factor-beta 1 on dentine pulp cells in different culture patterns. Cell Biol Int. 2006 Apr;30(4):295–300. pmid:16458025
- View Article
- PubMed/NCBI
- Google Scholar
36. Senbabaoglu Y, Michailidis G, Li JZ. Critical limitations of consensus clustering in class discovery. Sci Rep. 2014;4:6207. pmid:25158761
- View Article
- PubMed/NCBI
- Google Scholar
37. Ghosh AP, Klocke BJ, Ballestas ME, Roth KA. CHOP potentially co-operates with FOXO3a in neuronal cells to regulate PUMA and BIM expression in response to ER stress. PLoS One. 2012;7(6):e39586. pmid:22761832
- View Article
- PubMed/NCBI
- Google Scholar
38. Jousse C, Deval C, Maurin AC, Parry L, Cherasse Y, Chaveroux C, et al. TRB3 inhibits the transcriptional activation of stress-regulated genes by a negative feedback on the ATF4 pathway. J Biol Chem. 2007 May 25;282(21):15851–61. pmid:17369260
- View Article
- PubMed/NCBI
- Google Scholar
39. Leibinger M, Muller A, Andreadaki A, Hauk TG, Kirsch M, Fischer D. Neuroprotective and axon growth-promoting effects following inflammatory stimulation on mature retinal ganglion cells in mice depend on ciliary neurotrophic factor and leukemia inhibitory factor. J Neurosci. 2009 Nov 11;29(45):14334–41. pmid:19906980
- View Article
- PubMed/NCBI
- Google Scholar
40. Reiprich S, Wegner M. From CNS stem cells to neurons and glia: Sox for everyone. Cell Tissue Res. 2015 Jan;359(1):111–24. pmid:24894327
- View Article
- PubMed/NCBI
- Google Scholar
41. Di GS, Knoblach SM, Brandoli C, Aden SA, Hoffman EP, Faden AI. Gene profiling in spinal cord injury shows role of cell cycle in neuronal death. Ann Neurol. 2003 Apr;53(4):454–68. pmid:12666113
- View Article
- PubMed/NCBI
- Google Scholar
42. Befort K, Karchewski L, Lanoue C, Woolf CJ. Selective up-regulation of the growth arrest DNA damage-inducible gene Gadd45 alpha in sensory and motor neurons after peripheral nerve injury. Eur J Neurosci. 2003 Aug;18(4):911–22. pmid:12925017
- View Article
- PubMed/NCBI
- Google Scholar
43. Unschuld PG, Dachsel J, Darios F, Kohlmann A, Casademunt E, Lehmann-Horn K, et al. Parkin modulates gene expression in control and ceramide-treated PC12 cells. Mol Biol Rep. 2006 Mar;33(1):13–32. pmid:16636914
- View Article
- PubMed/NCBI
- Google Scholar
44. Nie X, Deng M, Yang M, Liu L, Zhang Y, Wen X. Axonal regeneration and remyelination evaluation of chitosan/gelatin-based nerve guide combined with transforming growth factor-beta1 and Schwann cells. Cell Biochem Biophys. 2014 Jan;68(1):163–72. pmid:23740553
- View Article
- PubMed/NCBI
- Google Scholar
45. Scherer SS, Kamholz J, Jakowlew SB. Axons modulate the expression of transforming growth factor-betas in Schwann cells. Glia. 1993 Aug;8(4):265–76. pmid:8406683
- View Article
- PubMed/NCBI
- Google Scholar
46. Baleriola J, Walker CA, Jean YY, Crary JF, Troy CM, Nagy PL, et al. Axonally synthesized ATF4 transmits a neurodegenerative signal across brain regions. Cell. 2014 Aug 28;158(5):1159–72. pmid:25171414
- View Article
- PubMed/NCBI
- Google Scholar
47. Puschel AW, Adams RH, Betz H. Murine semaphorin D/collapsin is a member of a diverse gene family and creates domains inhibitory for axonal extension. Neuron. 1995 May;14(5):941–8. pmid:7748561
- View Article
- PubMed/NCBI
- Google Scholar
48. Gong Y, Somwar R, Politi K, Balak M, Chmielecki J, Jiang X, et al. Induction of BIM is essential for apoptosis triggered by EGFR kinase inhibitors in mutant EGFR-dependent lung adenocarcinomas. PLoS Med. 2007 Oct 9;4(10):e294. pmid:17927446
- View Article
- PubMed/NCBI
- Google Scholar
49. Bonilla IE, Tanabe K, Strittmatter SM. Small proline-rich repeat protein 1A is expressed by axotomized neurons and promotes axonal outgrowth. J Neurosci. 2002 Feb 15;22(4):1303–15. pmid:11850458
- View Article
- PubMed/NCBI
- Google Scholar
50. Li S, Strittmatter SM. Delayed systemic Nogo-66 receptor antagonist promotes recovery from spinal cord injury. J Neurosci. 2003 May 15;23(10):4219–27. pmid:12764110
- View Article
- PubMed/NCBI
- Google Scholar
51. Tanabe K, Bonilla I, Winkles JA, Strittmatter SM. Fibroblast growth factor-inducible-14 is induced in axotomized neurons and promotes neurite outgrowth. J Neurosci. 2003 Oct 22;23(29):9675–86. pmid:14573547
- View Article
- PubMed/NCBI
- Google Scholar
52. Kim JE, Liu BP, Park JH, Strittmatter SM. Nogo-66 receptor prevents raphespinal and rubrospinal axon regeneration and limits functional recovery from spinal cord injury. Neuron. 2004 Oct 28;44(3):439–51. pmid:15504325
- View Article
- PubMed/NCBI
- Google Scholar
53. Huebner EA, Strittmatter SM. Axon regeneration in the peripheral and central nervous systems. Results Probl Cell Differ. 2009;48:339–51. pmid:19582408
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Mendez JS, Finn BW. Use of 6-hydroxydopamine to create lesions in catecholamine neurons in rats. J Neurosurg. 1975 Feb;42(2):166–73. pmid:234524
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Henry B, Crossman AR, Brotchie JM. Characterization of enhanced behavioral responses to L-DOPA following repeated administration in the 6-hydroxydopamine-lesioned rat model of Parkinson's disease. Exp Neurol. 1998 Jun;151(2):334–42. pmid:9628768
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Vu TQ, Ling ZD, Ma SY, Robie HC, Tong CW, Chen EY, et al. Pramipexole attenuates the dopaminergic cell loss induced by intraventricular 6-hydroxydopamine. J Neural Transm. 2000;107(2):159–76. pmid:10847557
View Article
PubMed/NCBI
Google Scholar

[10] View Article

[11] PubMed/NCBI

[12] Google Scholar

[ref4] 4. Ojeda V, Fuentealba JA, Galleguillos D, Andres ME. Rapid increase of Nurr1 expression in the substantia nigra after 6-hydroxydopamine lesion in the striatum of the rat. J Neurosci Res. 2003 Sep 1;73(5):686–97. pmid:12929136
View Article
PubMed/NCBI
Google Scholar

[14] View Article

[15] PubMed/NCBI

[16] Google Scholar

[ref5] 5. Ling ZD, Chang Q, Lipton JW, Tong CW, Landers TM, Carvey PM. Combined toxicity of prenatal bacterial endotoxin exposure and postnatal 6-hydroxydopamine in the adult rat midbrain. Neuroscience. 2004;124(3):619–28. pmid:14980732
View Article
PubMed/NCBI
Google Scholar

[18] View Article

[19] PubMed/NCBI

[20] Google Scholar

[ref6] 6. Spieles-Engemann AL, Steece-Collier K, Behbehani MM, Collier TJ, Wohlgenant SL, Kemp CJ, et al. Subthalamic nucleus stimulation increases brain derived neurotrophic factor in the nigrostriatal system and primary motor cortex. J Parkinsons Dis. 2011;1(1):123–36. pmid:22328911
View Article
PubMed/NCBI
Google Scholar

[22] View Article

[23] PubMed/NCBI

[24] Google Scholar

[ref7] 7. Gombash SE, Lipton JW, Collier TJ, Madhavan L, Steece-Collier K, Cole-Strauss A, et al. Striatal pleiotrophin overexpression provides functional and morphological neuroprotection in the 6-hydroxydopamine model. Mol Ther. 2012 Mar;20(3):544–54. pmid:22008908
View Article
PubMed/NCBI
Google Scholar

[26] View Article

[27] PubMed/NCBI

[28] Google Scholar

[ref8] 8. Gombash SE, Manfredsson FP, Mandel RJ, Collier TJ, Fischer DL, Kemp CJ, et al. Neuroprotective potential of pleiotrophin overexpression in the striatonigral pathway compared with overexpression in both the striatonigral and nigrostriatal pathways. Gene Ther. 2014 May 8.

[ref9] 9. Zhang Y, Meredith GE, Mendoza-Elias N, Rademacher DJ, Tseng KY, Steece-Collier K. Aberrant restoration of spines and their synapses in L-DOPA-induced dyskinesia: involvement of corticostriatal but not thalamostriatal synapses. J Neurosci. 2013 Jul 10;33(28):11655–67. pmid:23843533
View Article
PubMed/NCBI
Google Scholar

[31] View Article

[32] PubMed/NCBI

[33] Google Scholar

[ref10] 10. Spieles-Engemann AL, Behbehani MM, Collier TJ, Wohlgenant SL, Steece-Collier K, Paumier K, et al. Stimulation of the rat subthalamic nucleus is neuroprotective following significant nigral dopamine neuron loss. Neurobiol Dis. 2010 Jul;39(1):105–15. pmid:20307668
View Article
PubMed/NCBI
Google Scholar

[35] View Article

[36] PubMed/NCBI

[37] Google Scholar

[ref11] 11. Kordower JH, Olanow CW, Dodiya HB, Chu Y, Beach TG, Adler CH, et al. Disease duration and the integrity of the nigrostriatal system in Parkinson's disease. Brain. 2013 Aug;136(Pt 8):2419–31. pmid:23884810
View Article
PubMed/NCBI
Google Scholar

[39] View Article

[40] PubMed/NCBI

[41] Google Scholar

[ref12] 12. Seijffers R, Mills CD, Woolf CJ. ATF3 increases the intrinsic growth state of DRG neurons to enhance peripheral nerve regeneration. J Neurosci. 2007 Jul 25;27(30):7911–20. pmid:17652582
View Article
PubMed/NCBI
Google Scholar

[43] View Article

[44] PubMed/NCBI

[45] Google Scholar

[ref13] 13. Lin CR, Yang CH, Huang CE, Wu CH, Chen YS, Sheen-Chen SM, et al. GADD45A protects against cell death in dorsal root ganglion neurons following peripheral nerve injury. J Neurosci Res. 2011 May;89(5):689–99. pmid:21337369
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref14] 14. Starkey ML, Davies M, Yip PK, Carter LM, Wong DJ, McMahon SB, et al. Expression of the regeneration-associated protein SPRR1A in primary sensory neurons and spinal cord of the adult mouse following peripheral and central injury. J Comp Neurol. 2009 Mar 1;513(1):51–68. pmid:19107756
View Article
PubMed/NCBI
Google Scholar

[51] View Article

[52] PubMed/NCBI

[53] Google Scholar

[ref15] 15. Jing X, Wang T, Huang S, Glorioso JC, Albers KM. The transcription factor Sox11 promotes nerve regeneration through activation of the regeneration-associated gene Sprr1a. Exp Neurol. 2012 Jan;233(1):221–32. pmid:22024412
View Article
PubMed/NCBI
Google Scholar

[55] View Article

[56] PubMed/NCBI

[57] Google Scholar

[ref16] 16. Perez SE, Lazarov O, Koprich JB, Chen EY, Rodriguez-Menendez V, Lipton JW, et al. Nigrostriatal dysfunction in familial Alzheimer's disease-linked APPswe/PS1DeltaE9 transgenic mice. J Neurosci. 2005 Nov 2;25(44):10220–9. pmid:16267229
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref17] 17. Koprich JB, Chen EY, Kanaan NM, Campbell NG, Kordower JH, Lipton JW. Prenatal 3,4-methylenedioxymethamphetamine (ecstasy) alters exploratory behavior, reduces monoamine metabolism, and increases forebrain tyrosine hydroxylase fiber density of juvenile rats. Neurotoxicol Teratol. 2003 Sep;25(5):509–17. pmid:12972064
View Article
PubMed/NCBI
Google Scholar

[63] View Article

[64] PubMed/NCBI

[65] Google Scholar

[ref18] 18. Davis JF, Tracy AL, Schurdak JD, Tschop MH, Lipton JW, Clegg DJ, et al. Exposure to elevated levels of dietary fat attenuates psychostimulant reward and mesolimbic dopamine turnover in the rat. Behav Neurosci. 2008 Dec;122(6):1257–63. pmid:19045945
View Article
PubMed/NCBI
Google Scholar

[67] View Article

[68] PubMed/NCBI

[69] Google Scholar

[ref19] 19. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, et al. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004;5(10):R80. pmid:15461798
View Article
PubMed/NCBI
Google Scholar

[71] View Article

[72] PubMed/NCBI

[73] Google Scholar

[ref20] 20. Wilkerson MD, Hayes DN. ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010 Jun 15;26(12):1572–3. pmid:20427518
View Article
PubMed/NCBI
Google Scholar

[75] View Article

[76] PubMed/NCBI

[77] Google Scholar

[ref21] 21. Grabinski T, Kneynsberg A, Manfredsson FP, Kanaan NM. A method for combining RNAscope in situ hybridization with immunohistochemistry in thick free-floating brain sections and primary neuronal cultures. PLoS One. 2015;in press.

[ref22] 22. Kiryu-Seo S. Identification and functional analysis of damage-induced neuronal endopeptidase (DINE), a nerve injury associated molecule. Anat Sci Int. 2006 Mar;81(1):1–6. pmid:16526590
View Article
PubMed/NCBI
Google Scholar

[80] View Article

[81] PubMed/NCBI

[82] Google Scholar

[ref23] 23. Gonthier B, Nasarre C, Roth L, Perraut M, Thomasset N, Roussel G, et al. Functional interaction between matrix metalloproteinase-3 and semaphorin-3C during cortical axonal growth and guidance. Cereb Cortex. 2007 Jul;17(7):1712–21. pmid:17021275
View Article
PubMed/NCBI
Google Scholar

[84] View Article

[85] PubMed/NCBI

[86] Google Scholar

[ref24] 24. Marklund N, Fulp CT, Shimizu S, Puri R, McMillan A, Strittmatter SM, et al. Selective temporal and regional alterations of Nogo-A and small proline-rich repeat protein 1A (SPRR1A) but not Nogo-66 receptor (NgR) occur following traumatic brain injury in the rat. Exp Neurol. 2006 Jan;197(1):70–83. pmid:16321384
View Article
PubMed/NCBI
Google Scholar

[88] View Article

[89] PubMed/NCBI

[90] Google Scholar

[ref25] 25. Choi HJ, Kang KS, Fukui M, Zhu BT. Critical role of the JNK-p53-GADD45alpha apoptotic cascade in mediating oxidative cytotoxicity in hippocampal neurons. Br J Pharmacol. 2011 Jan;162(1):175–92. pmid:20955365
View Article
PubMed/NCBI
Google Scholar

[92] View Article

[93] PubMed/NCBI

[94] Google Scholar

[ref26] 26. Tanaka H, Shimazawa M, Kimura M, Takata M, Tsuruma K, Yamada M, et al. The potential of GPNMB as novel neuroprotective factor in amyotrophic lateral sclerosis. Sci Rep. 2012;2:573. pmid:22891158
View Article
PubMed/NCBI
Google Scholar

[96] View Article

[97] PubMed/NCBI

[98] Google Scholar

[ref27] 27. Steup A, Lohrum M, Hamscho N, Savaskan NE, Ninnemann O, Nitsch R, et al. Sema3C and netrin-1 differentially affect axon growth in the hippocampal formation. Mol Cell Neurosci. 2000 Feb;15(2):141–55. pmid:10673323
View Article
PubMed/NCBI
Google Scholar

[100] View Article

[101] PubMed/NCBI

[102] Google Scholar

[ref28] 28. Moreno-Flores MT, Martin-Aparicio E, Martin-Bermejo MJ, Agudo M, McMahon S, Avila J, et al. Semaphorin 3C preserves survival and induces neuritogenesis of cerebellar granule neurons in culture. J Neurochem. 2003 Nov;87(4):879–90. pmid:14622119
View Article
PubMed/NCBI
Google Scholar

[104] View Article

[105] PubMed/NCBI

[106] Google Scholar

[ref29] 29. Ara J, Bannerman P, Hahn A, Ramirez S, Pleasure D. Modulation of sciatic nerve expression of class 3 semaphorins by nerve injury. Neurochem Res. 2004 Jun;29(6):1153–9. pmid:15176472
View Article
PubMed/NCBI
Google Scholar

[108] View Article

[109] PubMed/NCBI

[110] Google Scholar

[ref30] 30. Fernandes KA, Harder JM, Kim J, Libby RT. JUN regulates early transcriptional responses to axonal injury in retinal ganglion cells. Exp Eye Res. 2013 Jul;112:106–17. pmid:23648575
View Article
PubMed/NCBI
Google Scholar

[112] View Article

[113] PubMed/NCBI

[114] Google Scholar

[ref31] 31. Fudge NJ, Mearow KM. Extracellular matrix-associated gene expression in adult sensory neuron populations cultured on a laminin substrate. BMC Neurosci. 2013;14:15. pmid:23360524
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref32] 32. Ma Y, Mi YJ, Dai YK, Fu HL, Cui DX, Jin WL. The inverse F-BAR domain protein srGAP2 acts through srGAP3 to modulate neuronal differentiation and neurite outgrowth of mouse neuroblastoma cells. PLoS One. 2013;8(3):e57865. pmid:23505444
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref33] 33. Ohoka N, Yoshii S, Hattori T, Onozaki K, Hayashi H. TRB3, a novel ER stress-inducible gene, is induced via ATF4-CHOP pathway and is involved in cell death. EMBO J. 2005 Mar 23;24(6):1243–55. pmid:15775988
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref34] 34. Ord D, Ord T. Characterization of human NIPK (TRB3, SKIP3) gene activation in stressful conditions. Biochem Biophys Res Commun. 2005 Apr 29;330(1):210–8. pmid:15781252
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref35] 35. Nie X, Tian W, Zhang Y, Chen X, Dong R, Jiang M, et al. Induction of transforming growth factor-beta 1 on dentine pulp cells in different culture patterns. Cell Biol Int. 2006 Apr;30(4):295–300. pmid:16458025
View Article
PubMed/NCBI
Google Scholar

[132] View Article

[133] PubMed/NCBI

[134] Google Scholar

[ref36] 36. Senbabaoglu Y, Michailidis G, Li JZ. Critical limitations of consensus clustering in class discovery. Sci Rep. 2014;4:6207. pmid:25158761
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref37] 37. Ghosh AP, Klocke BJ, Ballestas ME, Roth KA. CHOP potentially co-operates with FOXO3a in neuronal cells to regulate PUMA and BIM expression in response to ER stress. PLoS One. 2012;7(6):e39586. pmid:22761832
View Article
PubMed/NCBI
Google Scholar

[140] View Article

[141] PubMed/NCBI

[142] Google Scholar

[ref38] 38. Jousse C, Deval C, Maurin AC, Parry L, Cherasse Y, Chaveroux C, et al. TRB3 inhibits the transcriptional activation of stress-regulated genes by a negative feedback on the ATF4 pathway. J Biol Chem. 2007 May 25;282(21):15851–61. pmid:17369260
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref39] 39. Leibinger M, Muller A, Andreadaki A, Hauk TG, Kirsch M, Fischer D. Neuroprotective and axon growth-promoting effects following inflammatory stimulation on mature retinal ganglion cells in mice depend on ciliary neurotrophic factor and leukemia inhibitory factor. J Neurosci. 2009 Nov 11;29(45):14334–41. pmid:19906980
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref40] 40. Reiprich S, Wegner M. From CNS stem cells to neurons and glia: Sox for everyone. Cell Tissue Res. 2015 Jan;359(1):111–24. pmid:24894327
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref41] 41. Di GS, Knoblach SM, Brandoli C, Aden SA, Hoffman EP, Faden AI. Gene profiling in spinal cord injury shows role of cell cycle in neuronal death. Ann Neurol. 2003 Apr;53(4):454–68. pmid:12666113
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref42] 42. Befort K, Karchewski L, Lanoue C, Woolf CJ. Selective up-regulation of the growth arrest DNA damage-inducible gene Gadd45 alpha in sensory and motor neurons after peripheral nerve injury. Eur J Neurosci. 2003 Aug;18(4):911–22. pmid:12925017
View Article
PubMed/NCBI
Google Scholar

[160] View Article

[161] PubMed/NCBI

[162] Google Scholar

[ref43] 43. Unschuld PG, Dachsel J, Darios F, Kohlmann A, Casademunt E, Lehmann-Horn K, et al. Parkin modulates gene expression in control and ceramide-treated PC12 cells. Mol Biol Rep. 2006 Mar;33(1):13–32. pmid:16636914
View Article
PubMed/NCBI
Google Scholar

[164] View Article

[165] PubMed/NCBI

[166] Google Scholar

[ref44] 44. Nie X, Deng M, Yang M, Liu L, Zhang Y, Wen X. Axonal regeneration and remyelination evaluation of chitosan/gelatin-based nerve guide combined with transforming growth factor-beta1 and Schwann cells. Cell Biochem Biophys. 2014 Jan;68(1):163–72. pmid:23740553
View Article
PubMed/NCBI
Google Scholar

[168] View Article

[169] PubMed/NCBI

[170] Google Scholar

[ref45] 45. Scherer SS, Kamholz J, Jakowlew SB. Axons modulate the expression of transforming growth factor-betas in Schwann cells. Glia. 1993 Aug;8(4):265–76. pmid:8406683
View Article
PubMed/NCBI
Google Scholar

[172] View Article

[173] PubMed/NCBI

[174] Google Scholar

[ref46] 46. Baleriola J, Walker CA, Jean YY, Crary JF, Troy CM, Nagy PL, et al. Axonally synthesized ATF4 transmits a neurodegenerative signal across brain regions. Cell. 2014 Aug 28;158(5):1159–72. pmid:25171414
View Article
PubMed/NCBI
Google Scholar

[176] View Article

[177] PubMed/NCBI

[178] Google Scholar

[ref47] 47. Puschel AW, Adams RH, Betz H. Murine semaphorin D/collapsin is a member of a diverse gene family and creates domains inhibitory for axonal extension. Neuron. 1995 May;14(5):941–8. pmid:7748561
View Article
PubMed/NCBI
Google Scholar

[180] View Article

[181] PubMed/NCBI

[182] Google Scholar

[ref48] 48. Gong Y, Somwar R, Politi K, Balak M, Chmielecki J, Jiang X, et al. Induction of BIM is essential for apoptosis triggered by EGFR kinase inhibitors in mutant EGFR-dependent lung adenocarcinomas. PLoS Med. 2007 Oct 9;4(10):e294. pmid:17927446
View Article
PubMed/NCBI
Google Scholar

[184] View Article

[185] PubMed/NCBI

[186] Google Scholar

[ref49] 49. Bonilla IE, Tanabe K, Strittmatter SM. Small proline-rich repeat protein 1A is expressed by axotomized neurons and promotes axonal outgrowth. J Neurosci. 2002 Feb 15;22(4):1303–15. pmid:11850458
View Article
PubMed/NCBI
Google Scholar

[188] View Article

[189] PubMed/NCBI

[190] Google Scholar

[ref50] 50. Li S, Strittmatter SM. Delayed systemic Nogo-66 receptor antagonist promotes recovery from spinal cord injury. J Neurosci. 2003 May 15;23(10):4219–27. pmid:12764110
View Article
PubMed/NCBI
Google Scholar

[192] View Article

[193] PubMed/NCBI

[194] Google Scholar

[ref51] 51. Tanabe K, Bonilla I, Winkles JA, Strittmatter SM. Fibroblast growth factor-inducible-14 is induced in axotomized neurons and promotes neurite outgrowth. J Neurosci. 2003 Oct 22;23(29):9675–86. pmid:14573547
View Article
PubMed/NCBI
Google Scholar

[196] View Article

[197] PubMed/NCBI

[198] Google Scholar

[ref52] 52. Kim JE, Liu BP, Park JH, Strittmatter SM. Nogo-66 receptor prevents raphespinal and rubrospinal axon regeneration and limits functional recovery from spinal cord injury. Neuron. 2004 Oct 28;44(3):439–51. pmid:15504325
View Article
PubMed/NCBI
Google Scholar

[200] View Article

[201] PubMed/NCBI

[202] Google Scholar

[ref53] 53. Huebner EA, Strittmatter SM. Axon regeneration in the peripheral and central nervous systems. Results Probl Cell Differ. 2009;48:339–51. pmid:19582408
View Article
PubMed/NCBI
Google Scholar

[204] View Article

[205] PubMed/NCBI

[206] Google Scholar

Figures

Abstract

Introduction

Materials and Methods

Subjects

Unilateral intrastriatal 6-OHDA lesion

Tissue collection

High performance liquid chromatography (HPLC)

RNA isolation and quality evaluation

Microarray sample processing and hybridization

Microarray image analysis and quality control

Analysis of microarray data

Validation of microarray gene expression changes by qPCR

Pathway analysis of Cluster-1

Tyrosine hydroxylase (TH) immunostaining and in situ hybridization (ISH) of Cluster-1 transcripts to determine colocalization

Results

Lesion Status Confirmation

Consequences of striatal 6-OHDA injections on transcriptomic changes in the SN over 16 weeks

Microarray analysis.

qPCR validation of genes from Cluster-1.

Demonstration that exemplar Cluster-1 transcripts were localized to DA neurons of the SNpc.

Network analysis of cCluster-1

Discussion

Supporting Information

S1 Table. Fold changes at each time point, F-statistics and p values for the transcripts within Cluster-1.

S2 Table. Cluster assignment (1–8) and t-statistic at each time point for differentially expressed transcripts from the microarray analysis.

Acknowledgments

Author Contributions

References