- Open Access
A microarray study of MPP+-treated PC12 Cells: Mechanisms of toxicity (MOT) analysis using bioinformatics tools
BMC Bioinformatics volume 6, Article number: S8 (2005)
This paper describes a microarray study including data quality control, data analysis and the analysis of the mechanism of toxicity (MOT) induced by 1-methyl-4-phenylpyridinium (MPP+) in a rat adrenal pheochromocytoma cell line (PC12 cells) using bioinformatics tools. MPP+ depletes dopamine content and elicits cell death in PC12 cells. However, the mechanism of MPP+-induced neurotoxicity is still unclear.
In this study, Agilent rat oligo 22K microarrays were used to examine alterations in gene expression of PC12 cells after 500 μM MPP+ treatment. Relative gene expression of control and treated cells represented by spot intensities on the array chips was analyzed using bioinformatics tools. Raw data from each array were input into the NCTR ArrayTrack database, and normalized using a Lowess normalization method. Data quality was monitored in ArrayTrack. The means of the averaged log ratio of the paired samples were used to identify the fold changes of gene expression in PC12 cells after MPP+ treatment. Our data showed that 106 genes and ESTs (Expressed Sequence Tags) were changed 2-fold and above with MPP+ treatment; among these, 75 genes had gene symbols and 59 genes had known functions according to the Agilent gene Refguide and ArrayTrack-linked gene library. The mechanism of MPP+-induced toxicity in PC12 cells was analyzed based on their genes functions, biological process, pathways and previous published literatures.
Multiple pathways were suggested to be involved in the mechanism of MPP+-induced toxicity, including oxidative stress, DNA and protein damage, cell cycling arrest, and apoptosis.
DNA microarrays have been increasingly applied as a tool for the simultaneous monitoring of relative expression levels of thousands of genes for samples under various conditions, e.g., normal versus disease, and control versus drug or toxicant treatment [1–3], and offer a promising means to better understand how cells react to environmental perturbations. Their popularity, in part, is reflected by the number of microarray-related publications indexed in PubMed, which have been increasing exponentially. With the improvement of microarray chip quality and tools for data quality assurance, the focus in this field has gradually switched from determination of the gene numbers with altered expression levels to the analysis of the biological mechanism by categorizing significantly changed genes into functional groups and pathways [5–7]. Toxicogenomics, a emerging field combining genomics and bioinformatics to identify and characterize mechanisms of toxicity (MOT) for drugs and various compounds, has been developed quickly during the last several years [8–13]. Microarrays, with their power to examine genome-wide transcriptional responses, have become a key technology in toxicogenomics.
During the last two decades, MPTP (1-methyl-4-phenyl-1,2,5,6-tetrahydropyridine) – induced neurotoxicity has attracted a great deal of attention because of the similarity of its toxic effects to the biochemical changes in the brains of patients with Parkinson's disease (PD) [14, 15]. PD is a progressive neurodegenerative disorder that results in degeneration of dopaminergic neurons in the substantia nigra (SN) and dopamine depletion in the striatum . Although MPTP does not exactly reproduce PD, it has been an extremely valuable tool to model many features of PD in animals, and has led to a better understanding of crucial aspects of the sub-cellular events participating in the evolution of the PD clinical syndrome and neurotoxicity [17–20]. In order to be active, MPTP requires monoamine oxidase B to be converted into MPP+ (1-methyl-4-phenylpyridinium) . MPP+ is selectively taken up by dopaminergic neurons via the dopamine transporter of the plasma membrane [22, 23] and produces neuronal loss in substantia nigra (SN), striatal dopamine (DA) depletion and behavioral impairments in humans , primates , and mice [15, 25–28].
PC12 cells, a rat clonal pheochromocytoma cell line , possess dopamine synthesis, metabolism and transporter systems [30–32], and therefore have been used as a model for studies of MPP+ neurotoxicity and PD. A plethora of evidence has demonstrated that MPP+ depletes dopamine and elicits cell death in PC12 cells [33–36]. Previously, we demonstrated that MPP+-induced DA depletion and cell loss in PC12 cells is dose and time-course responsive . However, the mechanism of MPP+ neurotoxicity in PC12 cells is still unclear. Generally, it is believed that MPP+ directly and/or indirectly inhibits mitochondrial complex I, causing abnormal energy metabolism and increased production of reactive oxygen species (ROS), resulting in cell death [34, 36, 38]. Our study indicated that MPP+ may compromise heat shock protein (HSP) cell defense systems and cause apoptosis in NGF-differentiated PC12 cells and C56, but not CD1, mice [37, 39]. However, some studies have suggested that the effect of MPP+ on PC12 cells might be independent of ROS [35, 40]. Recent studies have also demonstrated that endoplasmic reticulum , PI-3  and cAMP pathways  may be possible targets for MPP+-induced neurotoxicity. Thus, MPP+ may cause PC12 cell injury and death via multiple complex mechanisms. Therefore, DNA microarray analysis, by measuring the expression of large numbers of genes, is a promising tool to help elucidate the mechanism of MPP+-induced toxicity in PC12 cells.
With the development of bioinformatics, various toxicogenomic databases and bioinformatics tools are available for data manipulation and mining [9, 44, 45]. The Center for Toxicoinformatics in the FDA's National Center for Toxicological Research (NCTR/FDA) has developed a public microarray data management and analysis software, called ArrayTrack, for FDA-wide microarray data storage and preliminary data analysis [10, 46]. In this study, alterations of gene expression were examined using the Agilent microarray platform and the data were analyzed using ArrayTrack and a software package called IRIDESCENT.
Methods and materials
Cell culture and MPP+ treatment
PC12 cells at passage 15–20 (ATCC, Manassas, VA) were grown in 75 cm2tissue culture flasks at 37°C under an atmosphere of 5% CO2/95% air in RPMI 1640 medium (Sigma, St. Louis, MO) containing 10% horse serum and 5% fetal bovine serum (complete media). At about 80% confluence, cells were placed into 24-well plates and cultured for another 2 days, and then were treated with 500 μM MPP+ iodide (Sigma, 100% purity by HPLC), a dose that induces significant decrease in dopamine content and cell viability .
Twenty-four hours after MPP+ treatment, 12 wells of PC12 cells from same treatment group in a 24-well plate were pooled into one single sample. Two control and two MPP+-treated samples were prepared for microarray hybridization. Total RNA was extracted using Qiagen™ RNeasy Mini Kits (Qiagen, Valencia, CA) following the RNeasy mini protocols for isolation from animal cells (I. spin protocol). RNA quantity was measured using the NanoDrop® ND-1000 UV-Vis Spectrophotometer (Wilmington, DE), while RNA quality was monitored using Agilent 2100 bioanalyzer and expressed as RNA integrity number (RIN) value using Agilent 2100 RIN Beta Version Software (Palo Alto, CA). For all 4 RNA samples, the ratios of 260/280 (indication of protein contamination) and 260/230 (indication of reagent contamination) determined by the NanoDrop spectrophotometer were above 2, and the RIN values determined by an Agilent bioanalyzer were 10 (Best RNA quality according to the grading of Agilent 2100 RIN Software).
Agilent 22K 60-mer oligonucleotide microarray slides (Palo Alto, CA) were used and a dye-swap experimental design was applied. High-quality RNA samples (200 ng each) acquired from PC12 cells were amplified and labeled with Cy5-and Cy3-CTP (Amersham Biosciences) to produce labeled cRNA using Agilent low RNA input fluorescent linear amplification kits following the manufacturers protocol. After the amplification and labeling, the dye-incorporation ratio was determined using a Nanodrop spectrophotometer and the ratios were within 10 to 20 pmol per μg cRNA, the range the manufacturer suggests prior to hybridization. For hybridization, the Agilent 60-mer oligo microarray processing protocol (Rev. 7, SSC Wash/6-screw hybridization chamber) was strictly followed. Briefly, 750 ng Cy3-labeled control and 750 ng Cy5-labeled MPP+-treated sample were mixed and incubated with an Agilent microarray slide for 17 hours using an Agilent in situ hybridization kit following SSC buffer washing. Sample pairs were dye-swapped and processed at the same time. The washed slides were immediately dried using an ultra pure filtered N2 stream in an ozone-free Biobubble. After drying, the slides were scanned using an Axon GenePix 4000B scanner with the PMT settings at 770 for Cy5 and 670 for Cy3, and the raw data were generated using GenePix Pro 6.0 software (Axon Instruments, Union City, CA).
Data analysis and MOT interpretation
Four datasets were acquired in this experiment, Slide 230: control 1 (Cy3)/treated 1(Cy5); Slide 231: control 2(Cy3)/treated 2(Cy5); Slide 232 (dye-swap of slide 230): control 1 (Cy5)/treated 1 (Cy3); Slide 233 (dye-swap of slide 231): control 2(Cy5)/treated 2(Cy3). Raw data generated from GenePix Pro 6.0 software was input into the ArrayTrack database. ArrayTrack is logically constructed of three linked components: a) a database (MicroarrayDB) that stores microarray experiment information; b) tools (TOOL) for data visualization and analysis; and c) libraries (LIB) that provide curated functional data from public databases for data interpretation. Specifically, ArrayTrack is MIAME (Minimum Information About A Microarray Experiment) supportive for storing both microarray data and experiment parameters associated with a toxicogenomics study. A quality control mechanism is implemented to assess the quality of entire expression data as well as quality of each spot on the chip. In addition, many data analysis and visualization tools are available with ArrayTrack, including four normalization methods, several statistical approaches for identification of differentially expressed genes, clustering/classification methods. ArrayTrack also provides a rich collection of functional information about genes, proteins and pathways drawn from various public biological databases for facilitating data interpretation [10, 46]. After data input, an automated data quality control was made by serial criteria of data quality control parameters in ArrayTrack to assure the data quality (Fig. 1A). Subsequently, data were normalized using the Lowess normalization method  to correct intensity-dependent ratio bias in ArrayTrack (Fig. 1B). Raw intensity data were logarithm (base 2) transformed and the log ratio of treated/control was calculated in ArrayTrack. The data set was then exported and the subsequent data analysis was performed employing JMP statistical software 5.0.1 and Spotfire DecisionSite bioinformatics software. The genes were categorized based on their functions derived from OMIM http://www.ncbi.nlm.nih.gov and PubMed literature http://www.pubmed.com as well as GO biological process and molecular functions.
A software package called IRIDESCENT [48–50] was also used to tie responding genes back to the published literature and identify commonalities among and between genes, diseases, chemical compounds, ontological categories and FDA-approved drugs. Commonalities are scored based upon how many genes within the entire set of microarray responders are observed to be related to each ''object'' via their co-occurrence within MEDLINE abstracts versus how many relationships would be expected statistically, by chance. This enables validation of experimental data by comparing the microarray response to the published literature to verify that at least some previously documented relationships are being reproduced by the experiment. It also enables the identification of common themes and potential new leads. Graphs were drawn using GraphViz, a software package written by AT&T Labs.
Microarray data consistency
Data consistency was examined employing both ArrayTrack tools and JMP software. Figure 2 shows the correlation ratios between sample pairs (r = 0.898 for slide 230 versus 231; r = 0.922 for slide 232 versus 233), indicating data consistency between biological samples. However, poor correlation ratios of dye-swapped slides (Slide 230/232 and slides 231/233) were observed in this study (Fig. 3) due to anti-correlation, a common phenomenon in dye-swapped slides. Further analysis using ArrayTrack tools indicated that these spots are consistent in the two swapped-slide sets (shown red) (Fig. 3A) and most of them are low-intensity spots and internal Cy3 controls (shown red) (Fig. 3B).
Identification of genes altered by MPP+ treatment
Genes with the lowest average intensity in this study (an arbitrary threshold at log 2 transformed intensity = 8) were excluded to reduce variance. Log ratios of the dye-swapped slides were averaged and the mean of two sample pairs were subsequently used to determine the fold change of gene expression in MPP+-treated groups in contrast to control. Figure 4 shows the correlation of the treated/control log (base 2) of the two sample pairs after the dye-swapped slides were averaged (r = 0.90), indicating low variance between the two samples sets. In this study, a 2-fold cutoff was applied as significant. Genes with log 2 transformed ratio (MPP+ treated/control) = 1 (representing 2-fold and above increase) showed red and ≤ 1 (representing 2-fold and above decrease) showed green. In total, there were 106 genes (44 genes induced and 62 genes repressed) with equal to and greater than a 2-fold change with MPP+ treatment.
Based on GO term information, OMIM and PubMed literature, ArrayTrack was used to categorize genes into 6 functional groups, which were: Oxidative stress, DNA and protein damage, metabolic process, neurotransmission and neuronal growth, cell arrest and apoptosis, and cell cycle (Tables 1, 2, 3, 4, 5, 6). Additionally, genes that cannot be categorized or with unclear functions were also listed (Table 7).
Genes associated with oxidative stress (Table 1)
Nine genes were functionally categorized in this group. Vmp1 is a stress-induced gene which may promote formation of intracellular vacuoles followed by cell death . Although it is still not clear the exact function of isg12(b), it is believed to play a role in resisting cellular or environmental stress. Herpud1 encodes an endoplasmic reticulum stress-inducible protein . Heat shock protein is generally induced by oxidative stress. Ephrin A1 acts as an early response protein in response to oxidative stress. This gene along with its receptor increased in response to DNA damage and is regulated by p53 in apoptosis . Ftl1 is a gene that is often induced due to iron-related oxidative stress, which is important for Parkinson's disease . The induction of these genes indicates that MPP+ may induce oxidative stress in PC12 cells. However, leukotriene B4 12-hydroxydehydrogenase which also plays a role in antioxidative function , CD24 involved in defense response  and Gstm1, a general stress-response gene, were down-regulated in response to MPP+ treatment, an outcome in opposition to the enhanced oxidative stress hypothesis. However, it should be noted that the above 3 genes are also involved in other functions. Alterations of these 3 genes may provide other mechanistic information for MPP+-induced toxicity.
Genes associated with DNA and protein damage (Table 2)
Both Myd116  and Ddit3 (OMIM: +126337) are DNA-damage inducible genes, and their induction provides direct evidence for DNA damage with MPP+ treatment. Dut is involved in purine metabolism, thus, the decrease of this enzyme may compromise the DNA-repair system. In addition, as mentioned above, increased expression of ephrin A1 is also an indication of DNA damage. Limited evidence was found to indicate protein damage. Tryptophanyl-tRNA synthetase catalyzes the first step of protein synthesis and is an essential function in the cell's protein synthesis machinery (OMIM: *191050). Cathepsin L is a protein associated with proteolysis and peptidolysis and induced by oxidative stress [58, 59]. These genes, along with heat shock protein 70, indicate cellular response to protein damage. Tcp1 functions as a cytosolic chaperone in the biosynthesis of tubulin [60, 61]. Although there is no evidence to indicate that this protein is involved in protein damage and repair, it has been reported that tubulin is involved in MPP+ elicited toxicity in PC12 cells . The decreased expression of this gene may affect the synthesis of tubulin, whose assembly into microtubules is critical to many cellular processes. It should be noted that LOC296368 has a sequence similar to ubiquitin-conjugation enzyme E2C, and is down-regulated with MPP+ treatment. However, because the exact function of this sequence is still unknown, the effect of these data on our hypothesis is unclear.
Genes associated with metabolic processes (Table 3)
If cellular functions are severely damaged, it is expected that most metabolic activities, especially fuel-utilization and energy metabolism, will be compromised. As Table 3 shows, MPP+ induces suppression of chromogranin A, a gene that modulates glucose, lipid and protein metabolism, and other genes associated with energy production, except acsl4. Therefore, we can primarily conclude that metabolic processes were down-regulated by MPP+ in PC12 cells. In addition, chromogranin A is also co-stored and co-released with catecholamines from storage granules in the adrenal medulla (OMIM: *118910). Although the function is still not known, the decrease of chromogranin A may indicate decreased dopamine content with MPP+ treatment as previously reported [33, 37].
Genes associated with neurotransmission and neuronal growth (Table 4)
MPP+-induces dopamine depletion, however, our preliminary study using real-time PCR has demonstrated that there are only minor gene expression changes with MPP+ treatment in PC12 cells (data not shown), unlike that in vivo . In this study, instead of tyrosine hydroxylase, dopa decarboxylase (Ddc) was down-regulated following MPP+ treatment. In addition, cytochrome b561, a major transmembrane protein that is specific to catecholamine and neuropeptide secretory vesicles of the adrenal medulla (OMIM: *600019), was down-regulated. As mentioned above, chromogranin A may be involved in dopamine metabolism in PC12 cells.
In PC12 cells, genes involved in neurogenesis and neuronal development were down-regulated following MPP+ treatment. These include Vgf (OMIM: *602186), Nnat , Plk2 , Hmgb2 (GO terms, biological process) and Bhlhb3 (GO terms, biological process). The decrease of these growth-related genes may suggest cell damage, cell growth arrest and a tendency toward apoptosis.
Genes associated with general cell growth arrest and apoptosis (Table 5)
With MPP+ treatment,, the expression of Gadd45a and Gas5 , genes associated with cell growth and responsive to DNA-damage stimulus , were induced. As shown above, induction of ephrin A1 also indicates p53-dependent activation in response to DNA damage . Furthermore, Trib3 (Kegg pathway) and Vmp1, which are involved in apoptosis, are both induced. Adenosine A2a receptor is also included in this group because it has been reported that Activation of A2A adenosine receptors (Adora2a) was found to prevent reactive oxygen species (ROS) formation and apoptosis in pheochromocytoma PC12 cells . The repression of this gene may facilitate apoptosis and ROS damage. In addition, this receptor is cAMP mediated, thus the repression of this receptor may be involved in subsequent cAMP regulated functions, such as CREB. ADM is one of the proteins that are cAMP pathway regulated. Although the function is still not clear, it has been reported to be highly expressed in pheochromocytoma and adrenal medulla (OMIM: *103275).
Genes associated with cell-cycle (Table 6)
It is evident that if cell growth is arrested, cell cycle and proliferation should be ceased. Thus, it is not surprising to see that genes for anti-proliferation (Btg1, Copeb) were induced and genes for cell cycle and replications were reduced. As shown in Table 6, Mcmd6 (mini chromosome maintenance deficient 6), which belongs to a family of early S-phase proteins required for DNA replication, may play a role in cell-cycle progression and DNA replication . Calcyclin (S100A6) is a calcium-binding protein that belongs to the family of S100 proteins. Its gene was discovered on the basis of its cell cycle-dependent expression. This gene is expressed at its maximal level during the transition between G0 to S phase of the cell cycle [69, 70]. Snf1lk encodes a protein kinase that may function in cell cycle regulation . Tcf19 is a transcription factor that may be involved in cell growth . H2Afz is a protein synthesized throughout the cell cycle (OMIM: *142763). These genes, along with some typical genes of the cell cycle (cdc2a, S-phase kinase-associated protein, Rfc2, Ret, Pcna), were all down-regulated by MPP+ treatment. Although the oncogene c-Myc has been considered a gene for proliferation, it is usually up-regulated during apoptosis [73, 74]. Thus, the induction of this gene is consistent with apoptotic processes induced by MPP+.
Genes with unknown functions (Table 7)
There are also some genes that cannot be categorized into gene groups. Generally this is due to their unknown functions. However, GPNMB, expressed in melanoma cell lines, was preferentially expressed in low-metastatic cell lines (OMIM: *604368). There was an inverse relationship between the expression of GPNMB compared to calcyclin (OMIM: *604368). In PC12 cells, this pattern was also observed. Acp5 is a gene related to bone growth (OMIM: *171640), which made it difficult to correlate with its function in PC12 cells.
Thus far, genes have been discussed individually and in terms of their known ontological processes. Gene categories shown above indicated that MPP+-induced toxicity might be related to the production of oxidative stress, DNA and protein damage, cell growth arrest, cell cycle/proliferation repression and apoptosis. IRIDESCENT was used to gain a broader perspective on how the 59 responding genes with known functions were related and whether these gene expression changes resonated with MPTP-induced neurotoxicity and Parkinson's disease. Only 47 of these 59 gene names, however, appeared in the literature at least once and could be used in the analysis based on the IRIDESCENT search. Genes with sequence similarity to other genes were not used (e.g. the LOC* entries). When examining how each of these genes was related to one another within the scientific literature, 34 genes had no detected relationships, suggesting that at least a few of these genes may be the result of noise or a non-specific response. It is also possible, of course, that their relationship to MPP+ and PD has not yet been documented and are thus valuable leads.
Discarding the genes without known connections, the most apparent pathway identified by the literature analysis was the DNA damage pathway (Figure 5), with GADD45 and GADD153 (Ddit3) in a relatively central role. Interestingly, these genes have been documented in other microarray experiments to play a role in dopamine-induced toxicity , yet have been only weakly associated with PD thus far. Both MPTP and PD are included in this graph so that literature relationships to the responding genes are apparent. Note the strong relationship between MPTP and PD and that the genes related to MPTP are also strongly related to PD. Roughly, the graph also segregates between a responding (up-regulated) and repressed (down-regulated) group of genes. Because we expect genes related by function or process in the literature to behave in a similar manner to a greater degree than unrelated genes, this helps confirm the specificity of the microarray experiments conducted. This graph also helps visualize to a greater degree the central players in the transcriptional response – in this case, PCNA, which has been called the "ringmaster of the genome" and can lead to apoptosis when downregulated . GADD45 is also upregulated, which is known to occur after DNA damage . Collectively, these factors within this analysis suggest that the MOT is likely due to DNA damage. As part of the shared relationship analysis, IRIDESCENT also identified Methyl Methanesulfonate (MMS) as being highly related to the responding genes. Because MMS is known to induce DNA damage through DNA alkylation, it would be interesting to compare the MPTP and MMS responses in future microarray experiments to identify both similarities and differences in their gene expression patterns.
During the last decade, DNA microarrays have been one of the hottest topics in biological research and biotechnology industry. The first paper about a DNA microarray was published from Pat Brown's laboratory in Stanford University in 1995 . By the end of 2004, the number of microarray-related publications was 11,625 (http://www.pubmed.com, keyword searching: DNA microarray). However, data assurance and data quality are still great challenges in microarray studies [4, 79]. In this study, the experiment was performed under stringent quality control, with minimum standards sets for RNA purity, RNA degradation, and dye-incorporation ratios. The correlation ratio between the two sample-sets was up to 0.9, thus, we assume that variance introduced by the experimental procedure was small. The anti-correlation effect exists more or less in all dye-swapped microarray experiments, which is due to natural differences between the Cy5 and Cy3 dyes. As our data indicated, the anti-correlation was intensity-dependent, and occurred primarily with low-intensity spots. The cut-off threshold applied in this study (log intensity ≥ 8) can, therefore, largely exclude these low intensity spots. This cutoff may cause false-negative discovery because the low-expressed genes, which may be altered by treatment, are also eliminated. However, it was interesting to note that none of the genes with expression levels lower than the cutoff level were changed up to 2-fold with MPP+ treatment. This is mainly due to the dye anti-correlation effect as shown in Figure 2.
In this study, the good correlation of sample set 1 and sample set 2 (Fig. 4) indicated low variance, which made it justifiable to use fold-change criterion to determine genes that were altered following MPP+ treatment. Actually, the correlation ratio of the two sample sets was up to 0.98 if only those genes with 2-fold change were counted (data not shown). Two-fold is a somewhat arbitrary fold-change criterion that has been applied in many microarray studies [80–83], but nonetheless necessary to delineate a tentative set of responders. We have to admit that it is an arbitrary criterion, and it does not mean that the genes altered lower than 2-fold change are not meaningful. Certainly, more genes will be determined to be significant if the threshold is decreased to 1.5 or lower, but this will produce more false discovery. Currently, this criterion is still the most acceptable criterion regarding the magnitude of fold change and biological significance.
Data interpretation (mechanism of toxicity) is the ultimate and most complicated aspect of a microarray study, in part because there is no standard approach. Although there are more and more tools, websites and biological libraries to assist MOT analysis, biological knowledge and understanding the pertinent literature are generally more important for MOT analysis. Our study indicated that MPP+ may produce its toxicity via the production of oxidative stress as indicated by induction of associated genes (Vmp1, Hspa5, Herpud1, Ftl1, Efna1, Ctsl, etc.), and at the same time suppression of general defense systems (Gstm1, Cd24, Adora2a, etc.). Previous studies have demonstrated that MPP+ elicits its toxicity in PC12 cells via the production of oxidative stress [84, 85], and this study provides evidence of oxidative stress production at the gene level. Oxidative stress may be responsible for the DNA and protein damage as represented by the induction of DNA and protein damage-inducible genes. It has been reported that MPP+ elicits DNA damage and subsequent apoptosis [86, 87] and induced heat shock proteins such as HSP70, which is an indicator of protein damage and repair [88, 89]. DNA and protein damage are involved in cell growth arrest and cell cycle repression and eventually apoptosis [90, 91]. A plethora of evidence has demonstrated that MPP+ can elicit cell death in PC12 cells [34, 36, 88]. Therefore, this study confirmed the previous findings at the gene expression level, and more importantly the genes identified in this study suggested specific mechanistic pathways of MPP+-induced toxicity.
In summary, this study provides new insight on the mechanism of toxicity induced by MPP+. Our data indicate that MPP+ induces oxidative stress, elicits damage in DNA and proteins, and causes cell growth arrest, cell cycle suppression and apoptosis. Literature analysis by IRIDESCENT has revealed some genes such as PCNA and GADD45 that are likely central players in this toxic response. In addition, IRIDESCENT analysis demonstrated that some of the genes altered by MPP+ in PC12 cells are involved in MPTP induced neurotoxicity in mice and Parkinson's disease. We have also identified MMS as a compound related to the genes in our experimental transcriptional response, and since no studies on the effects of MMS on PC12 cells have been published thus far, this is a potentially promising avenue for future experimentation the mechanism of toxicity for MPTP.
Gershon D: Microarray technology: an array of opportunities. Nature 2002, 416(6883):885–891. 10.1038/416885a
Lockhart DJ, Barlow C: Expressing what's on your mind: DNA arrays and the brain. Nat Rev Neurosci 2001, 2(1):63–68. 10.1038/35049070
Barlow C, Lockhart DJ: DNA arrays and neurobiology – what's new and what's next? Curr Opin Neurobiol 2002, 12(5):554–561. 10.1016/S0959-4388(02)00353-7
Shi L, Tong W, Goodsaid F, Frueh FW, Fang H, Han T, Fuscoe JC, Casciano DA: QA/QC: challenges and pitfalls facing the microarray community and regulatory agencies. Expert Rev Mol Diagn 2004, 4(6):761–777. 10.1586/14737184.108.40.2061
Moggs JG, Tinwell H, Spurway T, Chang HS, Pate I, Lim FL, Moore DJ, Soames A, Stuckey R, Currie R, et al.: Phenotypic anchoring of gene expression changes during estrogen-induced uterine growth. Environ Health Perspect 2004, 112(16):1589–1606.
Jung JW, Park JS, Hwang JW, Kang KS, Lee YS, Song BS, Lee GJ, Yeo CD, Kang JS, Lee WS, et al.: Gene expression analysis of peroxisome proliferators-and phenytoin-induced hepatotoxicity using cDNA microarray. J Vet Med Sci 2004, 66(11):1329–1333. 10.1292/jvms.66.1329
Lee GJ, Lee WS, Jeon KS, Um CH, Kim YS, Kim SJ, Lee CH, Yoon HK, Hwang SY, Park JS, et al.: cDNA microarray gene expression analysis and toxicological phenotype for anticancer drug. J Vet Med Sci 2004, 66(11):1339–1345. 10.1292/jvms.66.1339
Freeman K: Toxicogenomics data: the road to acceptance. Environ Health Perspect 2004, 112(12):A678–685.
Mattes WB, Pettit SD, Sansone SA, Bushel PR, Waters MD: Database development in toxicogenomics: issues and efforts. Environ Health Perspect 2004, 112(4):495–505.
Tong W, Cao X, Harris S, Sun H, Fang H, Fuscoe J, Harris A, Hong H, Xie Q, Perkins R, et al.: ArrayTrack – supporting toxicogenomic research at the U.S. Food and Drug Administration National Center for Toxicological Research. Environ Health Perspect 2003, 111(15):1819–1826.
Luyendyk JP, Mattes WB, Burgoon LD, Zacharewski TR, Maddox JF, Cosma GN, Ganey PE, Roth RA: Gene expression analysis points to hemostasis in livers of rats cotreated with lipopolysaccharide and ranitidine. Toxicol Sci 2004, 80(1):203–213. 10.1093/toxsci/kfh146
Lutz W, Kur B: [Toxicogenomics. New perspectives for the molecular toxicology]. Med Pr 2004, 55(2):193–202.
Waters MD, Fostel JM: Toxicogenomics and systems toxicology: aims and prospects. Nat Rev Genet 2004, 5(12):936–948. 10.1038/nrg1493
Burns RS, Chiueh CC, Markey SP, Ebert MH, Jacobowitz DM, Kopin IJ: A primate model of parkinsonism: selective destruction of dopaminergic neurons in the pars compacta of the substantia nigra by N-methyl-4-phenyl-1,2,3,6-tetrahydropyridine. Proc Natl Acad Sci U S A 1983, 80(14):4546–4550. 10.1073/pnas.80.14.4546
Heikkila RE, Hess A, Duvoisin RC: Dopaminergic neurotoxicity of 1-methyl-4-phenyl-1,2,5,6-tetrahydropyridine in mice. Science 1984, 224(4656):1451–1453.
Agid Y: Parkinson's disease: pathophysiology. Lancet 1991, 337: 1321–1324. 10.1016/0140-6736(91)92989-F
Dauer W, Kholodilov N, Vila M, Trillat AC, Goodchild R, Larsen KE, Staal R, Tieu K, Schmitz Y, Yuan CA, et al.: Resistance of alpha-synuclein null mice to the parkinsonian neurotoxin MPTP. Proc Natl Acad Sci U S A 2002, 99(22):14524–14529. 10.1073/pnas.172514599
Speciale S: MPTP: insights into parkinsonian neurodegeneration. Neurotoxicol Teratol 2002, 24(5):607–620. 10.1016/S0892-0362(02)00222-2
Krieglstein K, Suter-Crazzolara C, Unsicker K: Development of mesencephalic dopaminergic neurons and the transforming growth factor-β superfamily. J Neural Transm Suppl 1995, 46: 209–216.
Kuhn K, Wellen J, Link N, Maskri L, Lubbert H, Stichel CC: The mouse MPTP model: gene expression changes in dopaminergic neurons. Eur J Neurosci 2003, 17(1):1–12. 10.1046/j.1460-9568.2003.02408.x
Przedborski S, Jackson-Lewis V: Mechanisms of MPTP toxicity. Mov Disord 1998, 13(1):35–38.
Markey SP, Johannessen JN, Chiueh CC, Burns RS, Herkenham MA: Intraneuronal generation of a pyridinium metabolite may cause drug-induced parkinsonism. Nature 1984, 311(5985):464–467. 10.1038/311464a0
Chiba K, Trevor A, Castagnoli NJ: Metabolism of the neurotoxic tertiary amine, MPTP, by brain monoamine oxidase. Biochem Biophys Res Commun 1984, 120(2):574–578. 10.1016/0006-291X(84)91293-2
Davis GC, Williams AC, Markey SP, Ebert MH, Caine ED, Reichert CM, Kopin IJ: Chronic Parkinsonism secondary to intravenous injection of meperidine analogues. Psychiatry Res 1979, 1(3):249–254. 10.1016/0165-1781(79)90006-4
Tipton KF, Singer TP: Advances in our understanding of the mechanisms of the neurotoxicity of MPTP and related compounds [Review]. J Neurochem 1993, 61(4):1191–1206.
Ali SF, David SN, Newport GD, L CJ, Slikker WJ: MPTP-induced oxidative stress and neurotoxicity are age-dependent: evidence from measures of reactive oxygen species and striatal dopamine levels. Synapse 1994, 18(1):27–34. 10.1002/syn.890180105
Ali SF, David SN, Newport G: Age-related susceptibility to MPTP-induced neurotoxicity in mice. Neurotoxicology 1993, 14(1):29–34.
Xu Z, Cawthon D, McCastlain KA, Slikker W Jr, Ali SF: Selective alterations of gene expression in mice induced by MPTP. Synapse 2005, 55(1):45–51. 10.1002/syn.20089
Greene LA, Tischler AS: Establishment of a noradrenergic clonal line of rat adrenal pheochromocytoma cells which respond to nerve growth factor. Proc Nat Acad Sci 1976, 73: 2424–2428. 10.1073/pnas.73.7.2424
Tuler SM, Hazen AA, Bowen JM: Release and metabolism of dopamine in a clonal line of pheochromocytoma (PC12) cells exposed to fenthion. Fundam Appl Toxicol 1989, 13: 484–492. 10.1016/0272-0590(89)90284-4
Hatanaka H: Nerve growth factor-mediated stimulation of tyrosine hydroxylase activity in a clonal rat pheochromocytoma cell line. Brain Res 1981, 222(2):225–233. 10.1016/0006-8993(81)91029-5
Rebois RV, Reynolds EE, Toll L, Howard BD: Storage of dopamine and acetylcholine in granules of PC12, a clonal pheochromocytoma cell line. Biochemistry 1980, 19(6):1240–1248. 10.1021/bi00547a031
Itano Y, Kitamura Y, Nomura Y: Biphasic effects of MPP+, a possible parkinsonism inducer, on dopamine content and tyrosine hydroxylase mRNA expression in PC12 cells. Neurochem Int 1995, 26(2):165–171. 10.1016/0197-0186(94)00122-B
Sheng GQ, Zhang JR, Pu XP, Ma J, Li CL: Protective effect of verbascoside on 1-methyl-4-phenylpyridinium ion-induced neurotoxicity in PC12 cells. Eur J Pharmacol 2002, 451(2):119–124. 10.1016/S0014-2999(02)02240-9
Fonck C, Baudry M: Toxic effects of MPP(+) and MPTP in PC12 cells independent of reactive oxygen species formation. Brain Res 2001, 905(1–2):199–206. 10.1016/S0006-8993(01)02551-3
Hartley A, Stone JM, Heron C, Cooper JM, Schapira AH: Complex I inhibitors induce dose-dependent apoptosis in PC12 cells: relevance to Parkinson's disease. J Neurochem 1994, 63(5):1987–1990.
Xu Z, Cawthon D, McCastlain KA, Duhart HM, D NG, Fang H, Patterson TA, Slikker WJ, Ali S: Selective alterations of transcription factors in MPP+-induced neurotoxicity in PC12 cells. Nuerotoxicology 2005, in press.
Chalmers-Redman RM, Fraser AD, Carlile GW, Pong A, Tatton WG: Glucose protection from MPP+-induced apoptosis depends on mitochondrial membrane potential and ATP synthase. Biochem Biophys Res Commun 1999, 257(2):440–447. 10.1006/bbrc.1999.0487
Freyaldenhoven TE, Ali SF, Hart RW: MPTP-and MPP(+)-induced effects on body temperature exhibit age-and strain-dependence in mice. Brain Res 1995, 688(1–2):161–170. 10.1016/0006-8993(95)00529-Y
Lee CH, Han JH, Jang YY, Song JH, Han ES: Differential effect of catecholamines and MPP+on membrane permeability in brain mitochondria and cell viability in PC12 cells. Neurochem Int 2002, 40: 361–369. 10.1016/S0197-0186(01)00069-9
Ryu EJ, Harding HP, Angelastro JM, Vitolo OV, Ron D, Greene LA: Endoplasmic reticulum stress and the unfolded protein response in cellular models of Parkinson's disease. J Neurosci 2002, 22(24):10690–10698.
Salinas M, Martin D, Alvarez A, Cuadrado A: Akt1/PKBalpha protects PC12 cells against the parkinsonism-inducing neurotoxin 1-methyl-4-phenylpyridinium and reduces the levels of oxygen-free radicals. Mol Cell Neurosci 2001, 17(1):67–77. 10.1006/mcne.2000.0921
Ribeiro L, Azevedo I, Martel F: Comparison of the effect of cyclic AMP on the content and release of dopamine and 1-methyl-4-phenylpyridinium (Mpp+) in PC12 cells. Auton Autacoid Pharmacol 2002, 22(5–6):277–289. 10.1046/j.1474-8673.2002.00273.x
Waters MD, Olden K, Tennant RW: Toxicogenomic approach for assessing toxicant-related disease. Mutat Res 2003, 544(2–3):415–424.
Wetmore BA, Merrick BA: Toxicoproteomics: proteomics applied to toxicology and pathology. Toxicol Pathol 2004, 32(6):619–642. 10.1080/01926230490518244
Tong W, Harris S, Cao X, Fang H, Shi L, Sun H, Fuscoe J, Harris A, Hong H, Xie Q, et al.: Development of public toxicogenomics software for microarray data management and analysis. Mutat Res 2004, 549(1–2):241–253.
Yang YH, Dudoit S, Luu P, Lin DM, Peng V, Ngai J, Speed TP: Normalization for cDNA microarray data: a robust composite method addressing single and multiple slide systematic variation. Nucleic Acids Res 2002, 30(4):e15. 10.1093/nar/30.4.e15
Wren JD: Extending the mutual information measure to rank inferred literature relationships. BMC Bioinformatics 2004, 5(1):145. 10.1186/1471-2105-5-145
Wren JD, Garner HR: Shared relationship analysis: ranking set cohesion and commonalities within a literature-derived relationship network. Bioinformatics 2004, 20(2):191–198. 10.1093/bioinformatics/btg390
Wren JD, Bekeredjian R, Stewart JA, Shohet RV, Garner HR: Knowledge discovery by automated identification and ranking of implicit relationships. Bioinformatics 2004, 20(3):389–398. 10.1093/bioinformatics/btg421
Jiang PH, Motoo Y, Vaccaro MI, Iovanna JL, Okada G, Sawabu N: Expression of vacuole membrane protein 1 (VMP1) in spontaneous chronic pancreatitis in the WBN/Kob rat. Pancreas 2004, 29(3):225–230. 10.1097/00006676-200410000-00008
Kokame K, Agarwala KL, Kato H, Miyata T: Herp, a new ubiquitin-like membrane protein induced by endoplasmic reticulum stress. J Biol Chem 2000, 275(42):32846–32853. 10.1074/jbc.M002063200
Dohn M, Jiang J, Chen X: Receptor tyrosine kinase EphA2 is regulated by p53-family proteins and induces apoptosis. Oncogene 2001, 20(45):6503–6515. 10.1038/sj.onc.1204816
Gotz ME, Double K, Gerlach M, Youdim MB, Riederer P: The relevance of iron in the pathogenesis of Parkinson's disease. Ann N Y Acad Sci 2004, 1012: 193–208. 10.1196/annals.1306.017
Dick RA, Kwak MK, Sutter TR, Kensler TW: Antioxidative function and substrate specificity of NAD(P)H-dependent alkenal/one oxidoreductase. A new role for leukotriene B4 12-hydroxydehydrogenase/15-oxoprostaglandin 13-reductase. J Biol Chem 2001, 276(44):40803–40810. 10.1074/jbc.M105487200
Challen GA, Martinez G, Davis MJ, Taylor DF, Crowe M, Teasdale RD, Grimmond SM, Little MH: Identifying the molecular phenotype of renal progenitor cells. J Am Soc Nephrol 2004, 15(9):2344–2357. 10.1097/01.ASN.0000136779.17837.8F
Doutheil J, Althausen S, Gissel C, Paschen W: Activation of MYD116 (gadd34) expression following transient forebrain ischemia of rat: implications for a role of disturbances of endoplasmic reticulum calcium homeostasis. Brain Res Mol Brain Res 1999, 63(2):225–232. 10.1016/S0169-328X(98)00276-9
Kagedal K, Johansson U, Ollinger K: The lysosomal protease cathepsin D mediates apoptosis induced by oxidative stress. Faseb J 2001, 15(9):1592–1594.
Garcia-Arumi E, Andreu AL, Lopez-Hellin J, Schwartz S: Effect of oxidative stress on lymphocytes from elderly subjects. Clin Sci (Lond) 1998, 94(4):447–452.
Bourke GJ, El Alami W, Wilson SJ, Yuan A, Roobol A, Carden MJ: Slow axonal transport of the cytosolic chaperonin CCT with Hsc73 and actin in motor neurons. J Neurosci Res 2002, 68(1):29–35. 10.1002/jnr.10186
Roobol A, Holmes FE, Hayes NV, Baines AJ, Carden MJ: Cytoplasmic chaperonin complexes enter neurites developing in vitro and differ in subunit composition within single cells. J Cell Sci 1995, 108(Pt 4):1477–1488.
Cappelletti G, Incani C, Maci R: Involvement of tubulin in MPP+neurotoxicity on NGF-differentiated PC12 cells. Cell Biol Int 1995, 19(8):687–93. 10.1006/cbir.1995.1118
Zheng S, Chou AH, Jimenez AL, Khodadadi O, Son S, Melega WP, Howard BD: The fetal and neonatal brain protein neuronatin protects PC12 cells against certain types of toxic insult. Brain Res Dev Brain Res 2002, 136(2):101–110. 10.1016/S0165-3806(02)00348-6
Myer DL, Bahassi el M, Stambrook PJ: The Plk3-Cdc25 circuit. Oncogene 2005, 24(2):299–305. 10.1038/sj.onc.1208278
Shugart EC, Levenson AS, Constance CM, Umek RM: Differential expression of gas and gadd genes at distinct growth arrest points during adipocyte development. Cell Growth Differ 1995, 6(12):1541–1547.
Rees WD, Hay SM, Fontanier-Razzaq NC, Antipatis C, Harries DN: Expression of the growth arrest genes (GAS and GADD) changes during organogenesis in the rat fetus. J Nutr 1999, 129(8):1532–1536.
Huang NK: Adenosine A2A receptors regulate oxidative stress formation in rat pheochromocytoma PC12 cells during serum deprivation. Neurosci Lett 2003, 350(2):127–131. 10.1016/S0304-3940(03)00860-7
Sykes DE, Weiser MM: Rat intestinal crypt-cell replication factor with homology to early S-phase proteins required for cell division. Gene 1995, 163(2):243–247. 10.1016/0378-1119(95)00297-J
Nowotny M, Bhattacharya S, Filipek A, Krezel AM, Chazin W, Kuznicki J: Characterization of the interaction of calcyclin (S100A6) and calcyclin-binding protein. J Biol Chem 2000, 275(40):31178–31182. 10.1074/jbc.M001622200
Lesniak W, Swart GW, Bloemers HP, Kuznicki J: Regulation of cell specific expression of calcyclin (S100A6) in nerve cells and other tissues. Acta Neurobiol Exp (Wars) 2000, 60(4):569–575.
Stephenson A, Huang GY, Nguyen NT, Reuter S, McBride JL, Ruiz JC: snf1lk encodes a protein kinase that may function in cell cycle regulation. Genomics 2004, 83(6):1105–1115. 10.1016/j.ygeno.2003.12.007
Teraoka Y, Naruse TK, Oka A, Matsuzawa Y, Shiina T, Iizuka M, Iwashita K, Ozawa A, Inoko H: Genetic polymorphisms in the cell growth regulated gene, SC1 telomeric of the HLA-C gene and lack of association of psoriasis vulgaris. Tissue Antigens 2000, 55(3):206–211. 10.1034/j.1399-0039.2000.550303.x
Jiang B, Liu JH, Bao YM, An LJ: Hydrogen peroxide-induced apoptosis in pc12 cells and the protective effect of puerarin. Cell Biol Int 2003, 27(12):1025–1031. 10.1016/j.cellbi.2003.09.007
Tarabin V, Schwaninger M: The role of NF-kappaB in 6-hydroxydopamine-and TNFalpha-induced apoptosis of PC12 cells. Naunyn Schmiedebergs Arch Pharmacol 2004, 369(6):563–569. 10.1007/s00210-004-0938-1
Stokes AH, Freeman WM, Mitchell SG, Burnette TA, Hellmann GM, Vrana KE: Induction of GADD45 and GADD153 in neuroblastoma cells by dopamine-induced toxicity. Neurotoxicology 2002, 23(6):675–684. 10.1016/S0161-813X(02)00093-1
Paunesku T, Mittal S, Protic M, Oryhon J, Korolev SV, Joachimiak A, Woloschak GE: Proliferating cell nuclear antigen (PCNA): ringmaster of the genome. Int J Radiat Biol 2001, 77(10):1007–1021. 10.1080/09553000110069335
Xiao G, Chicas A, Olivier M, Taya Y, Tyagi S, Kramer FR, Bargonetti J: A DNA damage signal is required for p53 to activate gadd45. Cancer Res 2000, 60(6):1711–1719.
Schena M, Shalon D, Davis RW, Brown PO: Quantitative monitoring of gene expression patterns with a complementary DNA microarray. Science 1995, 270(5235):467–470.
Johnson K, Lin S: QA/QC as a pressing need for microarray analysis: meeting report from CAMDA'02. Biotechniques 2003, (Suppl):62–63.
Ivanov D, Dvoriantchikova G, Pestova A, Nathanson L, Shestopalov VI: Microarray analysis of fiber cell maturation in the lens. FEBS Lett 2005, 579(5):1213–1219. 10.1016/j.febslet.2005.01.016
Ludwikow A, Gallois P, Sadowski J: Ozone-induced oxidative stress response in Arabidopsis: transcription profiling by microarray approach. Cell Mol Biol Lett 2004, 9(4B):829–842.
Noh HS, Lee HP, Kim DW, Kang SS, Cho GJ, Rho JM, Choi WS: A cDNA microarray analysis of gene expression profiles in rat hippocampus following a ketogenic diet. Brain Res Mol Brain Res 2004, 129(1–2):80–87.
King HC, Sinha AA: Gene expression profile analysis by DNA microarrays: promise and pitfalls. Jama 2001, 286(18):2280–2288. 10.1001/jama.286.18.2280
Lee CS, Han ES, Lee WB: Antioxidant effect of phenelzine on MPP+-induced cell viability loss in differentiated PC12 cells. Neurochem Res 2003, 28(12):1833–1841. 10.1023/A:1026119708124
Virmani A, Gaetani F, Binienda Z, Xu A, Duhart H, Ali SF: Role of Mitochondrial Dysfunction in Neurotoxicity of MPP+: Partial Protection of PC12 Cells by Acetyl-l-Carnitine. Ann N Y Acad Sci 2004, 1025: 267–273. 10.1196/annals.1316.033
Qing H, Xu H, Wei Z, Gibson K, Li XM: The ability of atypical antipsychotic drugs vs. haloperidol to protect PC12 cells against MPP+-induced apoptosis. Eur J Neurosci 2003, 17(8):1563–1570. 10.1046/j.1460-9568.2003.02590.x
Sheng G, Pu X, Lei L, Tu P, Li C: Tubuloside B from Cistanche salsa rescues the PC12 neuronal cells from 1-methyl-4-phenylpyridinium ion-induced apoptosis and oxidative stress. Planta Med 2002, 68(11):966–970. 10.1055/s-2002-35667
Quigney DJ, Gorman AM, Samali A: Heat shock protects PC12 cells against MPP+toxicity. Brain Res 2003, 993(1–2):133–139. 10.1016/j.brainres.2003.09.004
Iwaki K, Chi SH, Dillmann WH, Mestril R: Induction of HSP70 in cultured rat neonatal cardiomyocytes by hypoxia and metabolic stress. Circulation 1993, 87: 2023–2032.
Mestril R, Dillmann WH: Heat shock and adaptive response to ischemia. Trend Cardiovasc Med 1991, 1: 240–244. 10.1016/1050-1738(91)90028-D
Bagchi D, Bagchi M, Hassoun EA, Stohs SJ: In vitro and in vivo generation of reactive oxygen species, DNA damage and lactate dehydrogenase leakage by selected pesticides. Toxicology 1995, 104: 129–140. 10.1016/0300-483X(95)03156-A
About this article
Cite this article
Xu, Z., Patterson, T.A., Wren, J.D. et al. A microarray study of MPP+-treated PC12 Cells: Mechanisms of toxicity (MOT) analysis using bioinformatics tools. BMC Bioinformatics 6 (Suppl 2), S8 (2005). https://doi.org/10.1186/1471-2105-6-S2-S8
- PC12 Cell
- Protein Damage
- Cell Growth Arrest
- Elicit Cell Death