|Home | About | Journals | Submit | Contact Us | Français|
The molecular mechanisms underlying the changes in the nigrostriatal pathway in Parkinson’s disease (PD) are not completely understood. Here, we use mass spectrometry and microarrays to study the proteomic and transcriptomic changes in the striatum of two mouse models of PD, induced by the distinct neurotoxins 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) and methamphetamine (METH). Proteomic analyses resulted in the identification and relative quantification of 912 proteins with two or more unique peptides and 86 proteins with significant abundance changes following neurotoxin treatment. Similarly, microarray analyses revealed 181 genes with significant changes in mRNA, following neurotoxin treatment. The combined protein and gene list provides a clearer picture of the potential mechanisms underlying neurodegeneration observed in PD. Functional analysis of this combined list revealed a number of significant categories, including mitochondrial dysfunction, oxidative stress response, and apoptosis. These results constitute one of the largest descriptive data sets integrating protein and transcript changes for these neurotoxin models with many similar end point phenotypes but distinct mechanisms.
Parkinson’s disease (PD) is a progressive, neurodegenerative disorder primarily affecting the extrapyramidal motor system of the central nervous system (CNS). The disorder is marked by the loss of nigrostriatal dopamine (DA) neurons in the substantia nigra pars compacta (SN), resulting in a significant decrease of DA content in the striatum. The diminished DA content is thought to be responsible for the effects of the disease on movement, including akinesia, rigidity, and tremors.1,2 Increased DA production using 3,4-dihydroxy-l-phenylalanine (l-DOPA) is the standard treatment for PD patients but has adverse long-term effects. Despite intensive research, the exact cause of PD and the mechanisms that lead to nigrostriatal dopaminergic neuron death remain largely elusive; however, evidence suggests the potential role of mitochondrial dysfunction and oxidative stress in the pathogenesis of PD.3,4 Thus, a deeper understanding of disease pathology at the molecular level is urgently needed.
Technological advances have allowed new assessments of disease pathology at a genomic scale, particularly in mouse models of brain disease. These developments have provided new insights into the associated molecular pathways of neurodegeneration.5-9 Large-scale analysis is also possible for proteins in brain disorders using mass-spectrometry (MS)-based proteomics.10 Proteomic profiling using 2D gel electrophoresis coupled with MS has been applied to Parkin-deficient mice and revealed several proteins involved in mitochondrial dysfunction and oxidative stress.11 However, there have been few studies that compare in parallel the transcriptomic and proteomic changes in neurodegeneration.
Here, we investigate the striatal protein abundance and gene expression changes for two different toxicological mouse models of PD created using 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) and methamphetamine (METH). MPTP has been shown to induce Parkinsonism in mice with recapitulation of the hallmark cellular pathology, death of dopaminergic neurons in the SN, which decreases DA input to the target organ, the striatum.12,13 MPTP probably acts through the inhibition of oxidative phosphorylation in dopaminergic neurons. 14 METH-treated mice also exhibit cell death of dopaminergic neurons in the SN, and although the parallels to the human disease are less strong, the behavioral traits in the mice are akin to PD patients.15,16 METH stimulates catecholamine release, although the mechanisms by which large doses result in a PD model are unclear.17,18
We evaluated protein and transcript abundance changes in the striatum of the two mouse models using liquid chromatography (LC)–MS and microarrays, respectively. To date, this is the most comprehensive, integrated data set examining these two models of PD. The data provide insights into the physiological consequences of neurotoxicity and the potential molecular changes in the striatum of PD patients. In particular, we observed downregulation of a number of proteins and genes involved in mitochondrial dysfunction and oxidative response in the striatum. Furthermore, downregulation of several antiapoptotic and upregulation of pro-apoptotic proteins/genes was observed, also implicating programmed cell death in the disorder.
Adult C57BL/6J male mice (8 week, 21–27 g) received four i.p. injections of MPTP-HCl (15 mg/kg per injection), METH-HCl (10 mg/kg per injection), or 200 μL of 0.9% sterile saline (Con) at 2 h intervals.18-20 Striata were bilaterally removed 7 days following injections, and samples were frozen immediately in liquid nitrogen.
Neurochemical parameters were obtained from a small portion of the left striatum of each mouse used for microarray analyses. Striatum samples were weighed wet, sonicated with 250 μL of 0.1 M perchloric acid, and centrifuged at 14 000 rpm for 15 min (4 °C). The supernatant was filtered through a 0.2 μm polytetrafluoroethylene (PTFE) filter, and an aliquot was diluted with water (1:2) for high-performance liquid chromatography (HPLC) analysis. The solid pellet was suspended in 1.0 mL of 0.2 M NaOH for the protein assay. For HPLC, we used an ESA Biosciences HPLC model 580 solvent delivery module (dual-piston pump) and an ESA Biosciences Coulochem II electrochemical detector with an analytical cell operating at +350 mV and 500 nA. The mobile phase consisted of acetonitrile/sodium phosphate monobasic buffer [75 mM sodium phosphate, 1.8 mM 1-octanesulfonic acid sodium salt (OSA), and 12 μM ethylenediaminetetraacetic acid, disodium salt dehydrate (EDTA)] at 9.5:90.5 (v/v) and pH 3.1 (aqueous phase). The guard column was Adsorbosphere HS, C18, 7.5 × 4.6 mm, 3 μm. Both columns were from Alltech Associates, Inc. The flow rate was 0.8 mL/min, and the injector loop volume was 20 μL.
For each mouse, the isolated striatal samples from the left and right brain were pooled, which provides an average ~20 μg of wet tissue. All striatal samples were homogenized and digested into peptides as previously described.10 Peptides from control mice were individually labeled with 16O, and peptide samples from MPTP- or METH-treated mice were individually labeled with 18O using postdigestion immobilization trypsin-catalyzed labeling as previously described (n 5, each group).21 After labeling, an aliquot of the control sample (~150 μg/aliquot) was mixed with an equal amount of either an 18O-labeled MPTP- or METH-treated sample. A total of 10 pairs of 16O/18O-labeled peptide samples were thus generated with 5 replicates for each neurotoxin. The labeled peptide samples were further subjected to cysteinyl-peptide enrichment (CPE) as described previously22 to generate a cysteinyl (Cys) fraction and a noncysteinyl (non-Cys) fraction for each sample. Both the Cys and non-Cys peptide samples were analyzed using a custom-built capillary LC system coupled online using an in-house manufactured electrospray ionization (ESI) interface to an 11.5 T Fourier transform ion cyclotron resonance (FTICR) mass spectrometer.
The analyses of quantitative LC–FTICR data sets were as previously described.21 The detected LC–MS features were identified as peptides by matching the normalized elution time (NET) and accurate mass measurements to the pre-established accurate mass and time (AMT) tag database within a 5 ppm mass error and a 2% NET error and were quantified on the basis of 16O/18O intensity ratios. The AMT tag database was created by global 2D LC–MS/MS profiling of the whole mouse brain as previously described.21 All abundance ratios were log2-transformed, and a global normalization between data sets was applied by shifting the mean log ratio of each data set to zero. The peptide identification and quantification results from Cys and non-Cys data sets were combined to achieve increased proteome coverage as well as protein sequence coverage. All quantified peptides were rolled up to nonredundant protein groups using ProteinProphet,23 and the abundance ratio for each protein group was calculated by averaging the ratio of multiple unique peptides stemming from the same protein group.24
A one-sample t test was applied for both protein abundance ratios and peptide abundance ratios against a mean value of zero to determine whether the protein displayed significant abundance changes following each neurotoxin treatment. A two-sample t test was applied to detect abundance differences between the two models. A protein is considered to have a significant abundance change when the protein has a p value less than 0.05 with at least two different peptides, with both showing a relatively good peptide level p value (average p value < 0.2).
Gene expression data were obtained from the pooled left and right striata of the three groups of mice: MPTP, METH, and Con (n 3, each group). Gene expression levels were obtained using GeneChip Mouse Genome 430A 2.0 array (Affymetrix, Inc., Santa Clara, CA). Striatum cDNA samples were prepared using the One-Cycle cDNA synthesis kit (Affymetrix, Inc.) according to the protocol of the manufacturer, with RNA from one mouse being applied to one microarray. Our complete data set is deposited in the Gene Expression Omnibus database (accession number GSE8030; http://www.ncbi.nlm.nih.gov/geo/).
A one-sample t test was applied for transcript abundance ratios against a mean value of zero to determine whether the transcript was significantly different, following each neurotoxin treatment. A two-sample t test was applied to detect transcript abundance differences between the two models. A false discovery rate (FDR) was used to adjust for multiple hypothesis testing in R (http://faculty.washington.edu/jstorey/qvalue/).25,26
To detect groups of proteins having statistically significant concordant changes, we used gene set enrichment analysis (GSEA).27 The data set was first screened against Gene Ontology (GO), Kyoto encyclopedia of genes and genomes (KEGG), and other annotations using http://www.babelomics.org.28 For proper handling of the missing values and estimation of the confidence of the changes for the selected groups, we performed a sign test using the BSDA package of the R environment for statistical computing (www.r-project.org).29
C57/BL6J mice were treated with MPTP, METH, or saline (control) (n 5, each group, MS; n 3, each group, microarrays). A small portion of the left striatum from each mouse was sampled using HPLC. A significant loss of DA and its metabolites was observed in MPTP- (55% decrease, t test 1.81 × 10−4, df 1, p < 0.0002) and METH-treated animals (71% decrease, t test 2.21 × 10−4, df 1, p < 0.0003), while levels of a different neurotransmitter, serotonin, and its metabolites remained unchanged (Supplementary Figure 1 in the Supporting Information).
To analyze the potential neurotoxin-induced protein abundance changes, we applied a global quantitative proteomic approach (Figure 1), which integrates 16O/18O labeling and Cys-peptide fractionation with the AMT tag strategy to achieve relatively good proteome coverage with quantification.21 Striata from MPTP, METH, and control mice were individually processed and labeled, which led to the generation of 10 16O/18O-paired samples. Each sample was fractionated into Cys and non-Cys samples, which were individually analyzed by LC–MS to identify statistically significant changes. An extensive mouse brain peptide/protein database was recently developed from a global characterization of the mouse brain proteome by liquid chromatography coupled with tandem mass spectrometry (LC–MS/MS).10
The analyses resulted in the identification of ~4600 unique peptides corresponding to 1614 proteins, with all proteins quantified in at least 4 of the 10 paired biological samples. Relative protein abundance is expressed as a ratio of neurotoxin-treated to control sample levels. Figure 2A shows the overall reproducibility of the analyses by comparing the correlation of the raw peptide intensities (only 18O intensities were used) between any two biological samples. As shown, the Pearson correlation coefficient is 0.94 0.02 within the same METH or MPTP models (intratoxin), while the observed correlation coefficient 0.89 0.03 between the two models (intertoxin) is slightly lower, as expected. The results suggest overall good reproducibility of the quantitative data.
Statistical analyses confidently identified 912 proteins with at least two unique peptides among the total 1614 proteins. Using a Student’s t test, 199 and 149 proteins were identified with expression ratios significantly different from a log2 ratio of zero (p values < 0.05) in MPTP and METH models, respectively. To reduce false positives, this list of proteins was further filtered by requiring the log2 ratio to be at least 0.3 (25% change) and good agreement between multiple peptides detected from the same peptides (average peptide p value < 0.2). A final list of 86 proteins had significant abundance changes following neurotoxin treatment. The lists of the total 1614 proteins and the 86 proteins with significant abundance changes are supplied as Supplementary Table 1 and Table 1/Supplementary Table 2 (Supplementary Tables are in the Supporting Information), respectively.
Figure 2B shows a heatmap of the abundance patterns for the proteins with significant changes. A large percentage of proteins displayed consistent abundance regulation in both the MPTP and METH-treated mice. The consistency between samples extended to multiple peptides detected from the same protein (Figure 2C). Parts A and B of Figure 3 show examples of significantly up- and downregulated proteins in MPTP- and METH-treated striata, respectively. Consistent with this observation, a significant correlation (R 0.6563, p < 10−20) was found in protein abundance changes in response to both neurotoxins (Figure 3C). One example illustrating the sensitivity of detection of quantitative protein changes is Purkinje cell protein 4 (Pcp4), for which we detected an average log2 downregulation of −0.44 0.18 (p 0.004) in MPTP (a 28% decrease) and −0.25 0.13 (p 0.035) in METH. This agrees well with a recent report of a ~30% decrease following MPTP.30
Transcript profiling of the striatum was performed for each of the three conditions using Affymetrix 430A 2.0 microarrays, with the pooled left and right striatum of one animal analyzed using one array. Following robust multichip analysis (RMA) normalization of the data across all experiments, we identified 34 significant differentially expressed genes (FDR < 0.05) that were upregulated and 29 downregulated genes in the MPTP experiment compared to control samples. We also identified 51 and 40 genes up- and downregulated in the METH-treated striata, respectively (FDR < 0.05). The two gene lists with transcript abundance changes are in Supplementary Table 3 in the Supporting Information. A significant overlap was also found between the striatal responses to MPTP and METH, with 17 upregulated and 10 downregulated genes shared between the two treatment groups (x2 104, df 1, p < 0.001). These genes likely represent a common response of the striatum to the loss of dopaminergic afferents caused by the two different neurotoxins (parts A and B of Figure 4).
Consistent with the overlap among the significant genes between MPTP and METH, significant correlation (R 0.718, p < 10−16) was also found between all gene expression changes for both neurotoxin treatment conditions (Figure 4C). Such consistency between MPTP and METH agrees well with the proteomic observation, which is anticipated from the similar biological consequences of the two neurotoxin treatments. One of the genes upregulated in common between MPTP and METH was glial fibrillary acidic protein (Gfap), where similar abundance changes at the protein level were also observed. Gfap is an astrocyte marker known to be upregulated in acute traumatic brain injury and neurological diseases.31,32
We used qRT-PCR analyses of MPTP- and METH-treated striata to confirm the expression changes for two of the most significant differentially regulated genes common to the two toxins, Gfap and related RAS viral (r-ras) oncogene homologue 2 (Rras2) (Figure 4D). Similar to the microarray results, Gfap and Rras2 were significantly upregulated compared to controls in response to both MPTP and METH treatments as judged using qRT-PCR (Gfap, t 0.489, df 1, p 0.035, MPTP and t 0.487, df 1, p 0.048, METH; Rras2, t 0.491, df 1, p 0.029, MPTP and t 0.489, df 1, p 0.036, METH).
When the genes showing significant transcript and protein abundance changes were compared, only two, Gfap and glutathione peroxidase 4 (Gpx4), displayed significant abundance changes at both the transcript and protein level. Gpx4, an antioxidant defense enzyme previously reported to have a protection role in oxidative stress-induced apoptosis,33,34 is downregulated for METH-treated mice based on both transcript and protein levels. Consistent with the lack of overlap, there was no significant correlation between relative protein and transcript changes for either MPTP (R 0.034, p 0.333) or METH (R 0.043, p 0.216). One of the explanations for this low correlation is that neither technique offers high reproducibility to measure relatively small changes (<30%), while the majority of the observed genes do not exhibit large changes at either of the protein or peptide levels. Additionally, translational and post-translational controls may also play a role.
To evaluate the relation between absolute transcript and protein levels, we pooled all three experimental groups (Figure 4E). We estimated absolute protein abundance based on the total number of observations for a given protein by summing the observation peptide count across all five replicates from MPTP and all five from METH. This count measure is an indirect semiquantitative indicator of protein abundance.24 Transcript levels were determined by averaging the signal intensities across all nine microarray data sets, that is, MPTP, METH, and controls. Despite the lack of overlap between transcripts and proteins regulated by the two neurotoxins, there was a significant positive correlation between absolute protein and transcript abundance for the 1200 genes in common (R 0.2889, p < 10−21). These results suggest that transcript levels can serve as a partial indicator for protein abundance, despite post-transcriptional regulation.
Taking advantage of the protein and mRNA expression data, we examined whether codon usage influenced translational efficiency. Efficiently expressed genes (46 genes) were identified, whose residuals were greater than 2 SD above the LOWESS fitted curve (Figure 4E and Supplementary Table 4 in the Supporting Information). An approximately equal number of genes (43 genes) were identified as inefficiently expressed, with residuals less than 1.4 SD below the curve. A comparison of codon usage revealed that the efficiently translated genes tended to employ more commonly used codons than the inefficiently translated proteins (Supplementary Figure 2 in the Supporting Information).
MicroRNAs (miRNAs) are also thought to play a role in influencing translational efficiency. In both the MPTP- and METH-treated animals, we identified miRNAs (7 and 23 miRNAs, respectively), whose targeted genes showed significantly different expression levels between treated and untreated striata (Supplementary Table 5 in the Supporting Information). A large number of target genes for each miRNA (hundreds) were modestly repressed (1–2%) in the treated samples compared to the control. This subtle repression of many target genes is thought to be characteristic of miRNA regulation.35 One of the miRNAs (has_mir_570) was significant in both treatment conditions, suggesting a possible shared translational regulatory mechanism.
Despite the low overlap between the transcript and protein data, the combined data provide a total of 195 altered genes following the two neurotoxin treatments, giving a more complete picture of expression regulation in the striatum. These regulated genes/proteins are implicated in a number of functional categories, including mitochondrial dysfunction/oxidative phosphorylation, oxidative stress, apoptosis/cell death, neurotransmission, and signal transduction (Table 2).
Downregulation for a number of genes involved in mitochondrial dysfunction was found in both MPTP- and METH-treated mice, including five subunits of complex I and several ATP synthase subunits, suggesting mitochondrial dysfunction in the striatum of both MPTP- and METH-treated mice. Cytochrome C1, a known pro-apoptotic protein for neuronal death, was also upregulated in MPTP-treated mice. A previous study demonstrated that complex I deficiency increases the releasable soluble pool of cytochrome C1 in mitochondria.36 Additionally, both F- and V-type ATPase complexes were significantly downregulated using GSEA (Supplementary Figure 3A in the Supporting Information). Downregulation of F-type ATPase subunits is indeed in concordance with mitochondrial dysfunction and a decrease in ATP production. At the same time, although V-type ATPase is not involved in ATP biosynthesis, it is involved in the charging of synaptic vesicles with neurotransmitters.
Oxidative stress is widely considered a major consequence of mitochondrial dysfunction, thus potentially playing an important role in neurodegeneration. 4 Genes with antioxidant activity (Gpx4 and Gstm5), and possibly regulated in response to oxidative damage (Dusp1 and Fos) were downregulated (Table 2). In addition, upregulation for several other oxidative stress response proteins, PI3K, Usp14, and Rras2, was observed. Rras2 has been previously suggested to regulate oxide species production. 37 PI3K has also been shown to have a role in regulating the toxic levels of reactive oxygen species generated by oxidative stress.38
Another group of proteins showing a statistically significant increase in abundance were related to ubiquitin-mediated protein degradation (Supplementary Figure 3B in the Supporting Information). An interesting feature of this group is that it contains both proteins associated with ubiquitin–ligation activity and proteins from the proteasome subunits. While most of the ubiquitin ligases showed upregulation, all four proteasome subunit proteins showed concordant abundance decreases following drug treatment. Overall, the data suggested that generalized cell destruction or degradation of misfolded proteins plays a prominent role in these mouse models of PD.
Defining the molecular changes in the striatum of PD brains is necessary to understand the disorder and design of new therapeutic approaches. Microarrays have been used in other papers to analyze the striatum in mouse models of PD.6-8,39-42 However, differing toxins, doses, and time points in these studies complicate the identification of a reliable profile specific for PD in mice. We designed an approach using two neurotoxins, MPTP and METH, which each induce Parkinsonism in mice but have distinct mechanisms of action. Changes in common to both neurotoxins may represent a more confident molecular signature of PD, while changes specific to each neurotoxin may reflect their individual toxicology. While it is not entirely known whether the 7 day time point is best for measuring the differences between the two drugs, a significant loss of DA in the target organ, the striatum, is observed in both the MPTP and METH models of PD.
To gain deeper insights into the cellular response of the two neurotoxins, we assayed transcript and protein levels for both agents using microarrays and global quantitative proteomics. These techniques identified a number of significantly up- and downregulated transcripts and proteins in response to MPTP and METH. Some of the significantly regulated proteins may represent a protein group containing several isoforms with the common peptide detected. The current approach cannot pinpoint exactly which isoform(s) contributed to the observed changes, and orthogonal approaches will be necessary to validate the differences. The regulated genes/proteins may be attractive candidates for biomarker development in PD.
High correlations were found between proteins only or between transcripts only in response to MPTP and METH, suggesting that the general cellular response to both neurotoxins is similar. Somewhat surprisingly though, a low correlation was observed between the proteomics and microarray data in each of the neurotoxin models. Only two genes (Gpx4 and Gfap) overlapped between the regulated gene lists of proteomics and microarrays. The low concordance between microarray and proteomic data could be due to multiple factors, including translational and post-translational regulations. It is also possible that the differences in the quantification accuracy and the ability to detect low abundance genes for the two technologies play a major role in the observed low concordance. Most genes detected by global proteomics are relatively abundant because of the dynamic range limitation of current proteomic technologies;43 however, our results demonstrate the ability to detect as low as ~30% abundance changes exemplified by the Pcp4 protein. Low percentage changes in protein abundances could mean even smaller mRNA changes, which may not be detectable.
Another possibility for the low correlation may be the different mechanisms of the two neurotoxins. For example, we observed nearly a 6-fold decrease in protein abundances for aldehyde dehydrogenase gene Aldh1a1, a gene specifically expressed in dopamine neurons,44 in MPTP-treated mice but only a 3-fold decrease in METH-treated mice. MPTP causes cell death in neurons that project to the striatum, while METH suspends cell function and DA production without killing the cells. This difference may account for the more severe decrease in protein observed in the MPTP-treated mice.
Our investigation also evaluated proteins to better understand the relationship between transcript and protein levels in the striatum following the loss of dopaminergic afferents. A total of 1614 proteins were confidently identified, and very good consistency was found between replicates and neurotoxins. Differentially regulated proteins in response to drug treatment were also similar.
Interestingly, the Gfap protein was found to be upregulated in response to both MPTP and METH treatment, congruent with the microarray results. There were a number of other proteins that showed a similar response to both toxins. One protein with increased levels was spermidine synthase (Srm), which enhances Af3-induced neurotoxicity through increased free-radical levels.45 As mentioned before, the Aldh1a1 protein was downregulated in response to both neurotoxins and is of particular interest because decreasing Aldh1a1 in human lens epithelial cells using siRNA increases the susceptibility of the cells to oxidative damage and apoptosis.46
To explore the relationship between transcript and protein abundance, we examined the absolute intensity for transcript expression across all experiments and compared them to the protein abundance levels for the 1200 genes in common between the MPTP and METH experiments. Although there were individual cases in which as much as a 26-fold abundance difference could be found between protein and transcript levels, a significant correlation (p < 10−21) between proteins and transcripts was found, suggesting that transcript levels can be used as a general indicator of protein abundance.
Our results also identified a number of proteins functionally associated with apoptosis and cell death. In general, downregulation of antiapoptotic proteins and upregulation of pro-apoptotic proteins were observed, although apoptosis is a process in which cellular localization may play as much of a role as the expression level of these factors.47 Increased apoptosis and cell death factors are likely the consequence of mitochondrial dysfunction and oxidative damage. Indeed, many of the proteins involved in mitochondrial dysfunction and oxidative damage are also suggested to play a role in apoptosis. For example, proteins Gpx4 and Gstm5, which function primarily as antioxidants also have a protective effect against apoptosis.33,34 Both proteins were downregulated. In addition, antiapoptotic mitochondrial heat-shock protein 10 kDa (Hspe1) had significant reduced protein abundances in both mouse models, again suggesting reduced protection against cell death. Conversely, the pro-apoptotic protein cytochrome C1 (Cyc1) was upregulated in MPTP-treated mice, also as a result of mitochondrial dysfunction. Consistent with these observations, oxidative stress and mitochondrial dysfunction are implicated as major contributors to Parkinsonism in both MPTP-48 and METH-treated49,50 mice.
A number of proteins (Table 2) previously implicated in pro-apoptotic activity were observed with increased protein or transcript abundances in both mouse models. For instance, the classic glycolytic protein glyceraldehyde-3-phosphate dehydrogenase (Gapdh) has been identified as a general mediator of one or more apoptotic cascades and promotes Lewy body formation.51 Cyc1 and Gapdh were only observed with significant abundance increases in MPTP-treated mice, suggesting a potentially higher level of oxidative damage and cell death. The upregulation of the calpain-2 (Capn2) protein in both MPTP- and METH-treated mice is indicative of increased endoplasmic reticulum (ER) stress. Capn2 protein is known to be pro-apoptotic through activation of caspase-12.52 Increased Capn2 expression and neuronal death in MPTP-treated mice has been previously observed.53
Importantly, microarray and proteomic measurements seem to be complementary by providing two different sets of regulated genes. The combined regulated gene and protein list provides a much clearer picture of the biological changes occurring in the striatum following the two drug treatments. Overall, our data provide clear evidence on mitochondrial dysfunction, increased oxidative stress and damage, misregulated protein degradation, increased apoptosis and cell death, and the potential activation of the astrocytic response. While we have focused on aspects of neurotoxicity, other protein or gene expression changes may be related to the loss of dopamine signaling in the striatum. Many of the novel proteins and genes may represent interesting targets for therapeutic intervention or further mechanistic studies. The discovery of common changes helped illuminate pathways relevant to PD, but pathways specific for each toxin were also uncovered. The use of transcript and protein profiling to identify the molecular changes occurring in toxins with similar pathophysiological end effects may be a general approach for differentiating the molecular pathology of disease models from agent-specific effects.
Our data shows that understanding mRNA changes is not enough, and further investigation using multimodal analysis will provide a clearer picture of what is happening in the mouse brain. Because the current AMT method may fail to identify post-translational modifications, subsequent proteomic analyses will use newer technology to include potential modifications into its AMT tag database.54 In addition, a time-course study will illuminate the evolving changes to the neurotoxins in the mouse brain, while a chronic drug-delivery system might better replicate damage caused by drug abuse in humans.
This work was supported by the NIH (RO1 DA015802 and U24 NS052108 to D.J.S. and RO1 NS050148 to D.J.S. and R.D.S.), Staglin Music Festival and NARSAD Young Investigator Award, Tobacco-Related Disease Research Program (11RT-0172), and Alzheimer’s Association (IIRG-02-3609 to D.J.S.). PNNL efforts were supported by the NIH National Center for Research Resources (RR18522 to R.D.S.), and proteomics measurements were performed in the Environmental Molecular Sciences Laboratory, a U.S. Department of Energy national scientific user facility located at PNNL in Richland, Washington. PNNL is a multiprogram national laboratory operated by Battelle Memorial Institute for the DOE under Contract DE-AC05-76RL01830.
Supporting Information Available: Significantly decreased DA levels in the striata of treated versus untreated animals (Supplementary Figure 1), gene set enrichment analysis for three significant categories from the LC–FTICR–MS data: ATP metabolism, protein degradation, and cell movement (Supplementary Figure 2), codon usage analysis for efficiently and inefficiently translated genes (Supplementary Figure 3), all quantified proteins (Supplementary Table 1), annotated list of significantly changed proteins in response to neurotoxin (Supplementary Table 2), significantly up- and downregulated transcripts in both MPTP- and METH-treated mice (Supplementary Table 3), efficiently translated genes (Supplementary Table 4), and miRNAs whose targeted genes had significantly different expression levels between treated and untreated striata (Supplementary Table 5). This material is available free of charge via the Internet at http://pubs.acs.org.