Next Article in Journal
Transfer of Anti-Rotavirus Antibodies during Pregnancy and in Milk Following Maternal Vaccination with a Herpes Simplex Virus Type-1 Amplicon Vector
Next Article in Special Issue
Distinct Antigen Delivery Systems Induce Dendritic Cells’ Divergent Transcriptional Response: New Insights from a Comparative and Reproducible Computational Analysis
Previous Article in Journal
Synthesis of 8-Aryl-O-methylcyanidins and Their Usage for Dye-Sensitized Solar Cell Devices
Previous Article in Special Issue
SpidermiR: An R/Bioconductor Package for Integrative Analysis with miRNA Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Epidermal Growth Factor Receptor Related Gene Expression Changes in a Cellular and Animal Model of Parkinson’s Disease

Department of Biotechnology, Konkuk University, Chungju 380-701, Korea
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2017, 18(2), 430; https://doi.org/10.3390/ijms18020430
Submission received: 24 October 2016 / Revised: 27 January 2017 / Accepted: 8 February 2017 / Published: 16 February 2017
(This article belongs to the Special Issue Transcriptome Profiling in Human Diseases)

Abstract

:
We employed transcriptome analysis of epidermal growth factor receptor related gene expression changes in cellular and animal models of Parkinson’s disease (PD). We used a well-known Parkinsonian toxin 1-methyl-4-phenylpyridine (MPP+) to induce neuronal apoptosis in the human neuroblastoma SH-SY5Y cell line. The MPP+-treatment of SH-SY5Y cells was capable of inducing neuro-apoptosis, but it remains unclear what kinds of transcriptional genes are affected by MPP+ toxicity. Therefore the pathways that were significantly perturbed in MPP+ treated human neuroblastoma SH-SY5Y cells were identified based on genome-wide gene expression data at two time points (24 and 48 h). We found that the Epidermal Growth Factor Receptor (EGFR) pathway-related genes showed significantly differential expression at all time points. The EGFR pathway has been linked to diverse cellular events such as proliferation, differentiation, and apoptosis. Further, to evaluate the functional significance of the altered EGFR related gene expression observed in MPP+-treated SH-SY5Y cells, the EGFR related GJB2 (Cx26) gene expression was analyzed in an MPP+-intoxicated animal PD model. Our findings identify that the EGFR signaling pathway and its related genes, such as Cx26, might play a significant role in dopaminergic (DAergic) neuronal cell death during the process of neuro-apoptosis and therefore can be focused on as potential targets for therapeutic intervention.

Graphical Abstract

1. Introduction

Parkinson’s disease (PD) is the second most long-term degenerative disorder of the central nervous system and is associated with a progressive loss of dopaminergic (DAergic) neurons in the substantia nigra. Although, the pathogenesis of PD is not yet clear, evidence suggests that oxidative stress-induced mitochondrial dysfunction, associated with inhibition of respiratory chain complex I, is one of the main reasons for the onset of PD [1,2]. Experimentally, the symptoms of PD have been achieved using neurotoxin 1-methyl-4-phenylpyridinium (MPP+), which selectively and potently inhibits complex I of the mitochondrial electron transport chain [3,4,5,6]. The activation of neuronal cell death pathways involving the mitochondrial dysfunction, as well as increased oxidative stress, most likely represents the process of PD neurodegeneration [7]. Mitochondrial dysfunction results in reduced intracellular adenosine triphosphate (ATP), lower mitochondrial membrane potential, and less endogenous cellular respiration with an augmentation of oxidative stress [8].
Several toxin-induced in vitro models such as 1-methyl-4-phenylpyridinium (MPP+), 6-hydroxydopamine (6-OHDA), or rotenone have been widely used in studies to understand the pathogenesis of PD, and, in many aspects, the DAergic neuronal death observed in these conditions is similar to that observed in PD. A reliable cellular model enables large-scale functional screening of genes and interpretation of the molecular signatures through transcriptional profiling. Conn et al. [9] have employed cDNA microarray analysis to characterize the transcriptional response of SH-SY5Y cells to MPP+. However, their findings were limited because they examined only 1185 genes for differential expression. Thus, a genome-wide gene expression study is required to identify biological pathways affected by MPP+ at the cellular level. Identifying the molecular pathways involved in MPP+ toxicity could provide a better understanding of the molecular mechanisms involved in the cell death of PD. Therefore, we compared the expression profiles of MPP+ induced genes of an untreated group with a MPP+ treatment group, using whole genome expression array (Human Transcriptome Pico Assay 2.0, Affymetrix, Inc., San Diego, CA, USA). Bioinformatics analysis revealed many differentially expressed genes involved in important biological processes. Expression analysis was used to confirm the differential expression of various genes and proteins.

2. Results

2.1. MPP+ Induces Dopaminergic Neuronal Cell Death in SH-SY5Y

Cell viability was measured with MTT (5-diphenyl-tetrazolium bromide) assay based on the conversion of MTT to dark blue formazan crystals by mitochondrial dehydrogenase [10]. Exposure to MPP+ showed that concentration and time decreased cell viability when compared to untreated control cells. MPP+ at 0.5, 1, and 2 mM significantly (p < 0.05 and p < 0.01) reduced cell viability when compared to control cells. The higher the concentration of MPP+ used, the greater the decrease in cell viability (Figure 1A). The final concentration of MPP+ at 1 mM was selected for these trials. As shown in Figure 1B, MPP+ was significantly (p < 0.05 at 24 h and p < 0.01 at 48 h) cytotoxic to SH-SY5Y cells in a time-dependent manner. The percentage of cell viability was approximately 50% after 48 h treatment with 1 mM MPP+. Exposure level of cell death was also assayed by LDH (Lactate dehydrogenase) assay, which detected the release of LDH in culture medium. The exposure to 1 mM MPP+ for 48 h resulted in a 2.01 ± 0.11-fold increase in the release of LDH when compared to the control cultures (p < 0.01 vs. control group) (Figure 1D).
The Bcl-2 family plays a key role in the mitochondrial apoptotic pathway [11]. PARP (Poly ADP ribose polymerase) is an abundant nuclear enzyme with normal function in DNA repair. However, extensive PARP activation can promote cell death [12]. PARP is a biochemical marker for the execution phase of apoptosis. Therefore, we investigated whether MPP+ had an apoptotic effect on the expressions of the Bcl-2 family and cleaved PARP in SH-SY5Y cells using RT-PCR and immunoblot analysis. As shown in Figure 1C, the anti-apoptotic molecule Bcl-2 mRNA and the cleaved PARP protein in MPP+-treated SH-SY5Y cells were markedly reduced in a time-dependent manner, whereas pro-apoptotic molecule Bax mRNA was increased in a time-dependent manner. Further, there is a time-dependent increase in the Bax/Bcl-2 ratio in cells exposed to 1 mM MPP+ (Figure 1E). The ratio of anti-apoptotic versus pro-apoptotic Bcl-2 family proteins regulates apoptosis sensitivity and could be partly responsible in the cell fate decisions [12].

2.2. Differential Gene Expression and Pathway from MPP+-Treated SH-SY5Y Cells at Two Time Points (24 and 48 h)

To investigate the responses of SH-SY5Y cells exposed to MPP+ and understand the underlying molecular mechanisms of MPP+ action in SH-SY5Y cells, we used microarrays to determine the gene expression in the two experimental groups; MPP+ treatment for 24 and 48 h. Statistically significant (p < 0.05) changes in gene expression (fold change ≥2 relative to stimulated cells) were identified by microarray analysis. Many genes were up- or down- regulated by more than two fold by MPP+. A total of 6.4% of genes were differentially expressed by MPP+ induced cytotoxicity. Among these genes, MPP+ treatment t for 24 and 48 h up-regulated 0.5% and 1.3% of the genes in SH-SY5Y cells, respectively, evoking the common expression of 40 genes. MPP+ treatment for 24 and 48 h also down-regulated 2.8% and 1.8% of genes, respectively, showing common expression of 0.8% of genes, including the following mitochondrial controlled apoptosis-related genes; BAG1 (BCL-2 binding athanogene-1), BAG5 (B-cell CLL/lymphoma 7C), and BCL-2. Anti-apoptotic molecules were down-regulated by MPP+ induced cell death (Table 1).
Differentially expressed genes (DEGs) were sorted based on the direction of regulation and corresponding Gene Ontology (GO) categories (Figure 2). The use of the DAVID bioinformatics resources (NIAID/NIH, 2008, available online: http://david. abcc.ncifcrf.gov) enabled us to detect gene ontology biological process categories for these differentially expressed genes. The results revealed proteins encoded by these genes associated with several categories (Figure 2). The biggest GO biological process classification of down-regulated genes after MPP+ treatment for 24 and 48 h was for genes encoding proteins associated with signal transduction. Between down-regulated genes after MPP+ treatment for 24 and 48 h, the common GO biological processes were signal transduction, cell communication, other intracellular signaling cascades, co-enzyme and prosthetic group metabolism, vitamin/co-factor transport, and steroid hormone-mediated signaling.
GO biological process classification of up-regulated genes was distinct between the two MPP+ treatment times. The biggest GO biological process classification of up-regulated genes after MPP+ treatment for 24 h was for genes encoding proteins associated with signal transduction. However, the GO biological process classification of up-regulated genes after MPP+ treatment for 48 h was for genes encoding proteins associated with nucleoside, nucleotide, and nucleic acid metabolism (Figure 2C). These results suggest that signal transduction plays a leading role in MPP+ induced DAergic neuronal cell death.

2.3. Validation of Microarray Results

To increase the fidelity of the microarray data, all gene expression values were measured using independent samples at each time point. We selected twelve genes from Table 1 and subjected them to internal validation of expression. The expression patterns of these genes agreed very well with the microarray results (Figure 3). Some of the identified genes were of interest from the context of apoptosis. The TNFRSF12A (fibroblast growth factor-inducible 14 receptor) binds to the soluble TNF-like weak inducer of apoptosis (TWEAK) and induces multiple biological responses, such as proliferation, migration, differentiation, angiogenesis, and apoptosis [13]. In particular, Amyloid β A4 protein precursor (APP) has been recently identified as an integral membrane protein implicated as a regulator of synapse formation, neural plasticity, and iron export [14]. Abnormal metabolism of APP results in β-amyloid accumulation, causing neuronal degeneration in the brains of AD (Alzheimer's disease) patients [15]. Lie et al. [16] reported that over expression of APP will disrupt mitochondrial membrane potentials and redox states.

2.4. Analysis of EGFR Pathway-Related Gens in MPP+-Induced DAergic Neuronal Cell Death

The use of the KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway database (available online: http://www.genome.jp/kegg/pathway) and NetPath (available online: http://www.netpath.org) enabled us to detect pathways for the differentially expressed genes identified in this study. Various pathways including the EGFR pathway, the Notch pathway, a two-component system, and neuroactive ligand-receptor interaction were found for those differentially expressed genes. In particular, the EGFR signaling pathway-related genes showed significant expression at all time points. We selected twelve genes associated with MPP+ induced dopaminergic neuronal cell death (Table 2). Among these genes, three genes (HSPB1, IGFBP3, and GPRC5A) were down regulated (Figure 4). HSPB1, IGFBP3, and ZFP36L2 were associated with the p38 MAPK pathway in the EGFR pathway involved in various stresses including reactive oxygen species (ROS), ultra-violet (UV), and toxin [17,18,19]. Other genes (MELK, IFI44L, SHMT2, SAA4, ZFP36L2, CAV1, GJB2, LAD1, and ACOT7) associated with the EGFR signaling pathway were up-regulated in MPP+ induced apoptosis (Figure 4). Also two genes (GPRC5A and CJB2) were associated with the G-protein coupled receptor (GPCR) protein signaling pathway. The EGF signaling pathway regulated the GPCR-related genes [20].

2.5. Expression of GJB2 (Cx26) in In Vitro and In Vivo Pakinson’s Model

To investigate the expression of Cx26 in toxin induced apoptosis and the MPTP (1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine) animal model, we analyzed expression of Cx26 for selected treated times. As shown in Figure 5, there is a significant (p < 0.05) and time-dependent increase in the protein expression of Cx26 with the increase in neuronal cell death. After treatment with other toxins such as hydrogen peroxide (H2O2) and 6-OHDA (100 µM H2O2 and 25 µM 6-OHDA) for up to 48 h, almost half of the neuronal cell loss was observed, which was similar to MPP+-induced cell death. This can result in the generation of reactive oxygen species (ROS). It was well documented that ROS is involved in the pathophysiology of several neurodegenerative diseases, including PD [21,22]. There is a time-dependent increase in the expression of Cx26 protein on toxin treated SH-SY5Y cells, which was similar to that in MPP+ treated cells.
The injection of the MPTP neurotoxin has become a widely-used method to investigate the neurodegenerative process. To investigate whether Cx26 levels may be altered in Parkinsonism, we examined the Cx26 protein level in our MPTP-induced mouse PD model. In this model, chronic exposure to MPTP remarkably reduced tyrosine hydroxylase (TH) expression in the substantia nigra pars compacta (SNpc), the same area where a loss of DA neurons occurs in human PD [23]. Figure 6 shows that Cx26 was found in the midbrain and striatum and that the Cx26 protein level was up-regulated in a time-dependent manner. Results indicated that expression of Cx26 was induced by toxins; its role, especially linked on neuronal cell death and neurodegeneration, was evident in both in vitro and in vivo models.

3. Discussion

It is well documented that that DAergic neuronal cell death underlies the symptoms of Parkinson’s diseases [24]. The development of a stable and reliable DAergic neuronal injury model is particularly necessary for studying the pathogenesis of PD. Human neuroblastoma SH-SY5Ycells mimic many aspects of the DAergic neuronal cell death observed in PD when treated with neurotoxins such as MPP+, 6-hydroxydopamine, and rotenone. We investigated the global gene expression pattern in MPP+ treated SH-SY5Y cells. Microarray analysis revealed that 6.4% genes of cellular transcripts were differentially regulated in MPP+ induced DAergic neuronal cell death. These differentially expressed genes were sorted based on the direction of regulation. Their corresponding GO categories were identified. The GO terms of down-regulated genes in MPP+ induced apoptosis were mostly signal transduction. For up-regulated genes, the GO terms are different. After MPP+ treatment for 24 h, most of those up-regulated genes were involved in signal transduction. After 48 h of treatment with MPP+, most of the up-regulated genes were involved in nucleoside, nucleotide, and nucleic acid metabolism. The expression results of sixteen differentially-expressed genes identified by microarray analysis were confirmed by RT-PCR, confirming the reliability of the microarray approach.
Genes in the related EGFR signal pathway of DEGs were found to be differentially expressed with different exposure time to MPP+. Multiple alterations of gene expressions and signal transduction pathways lie downstream of activated EGFRs, and some of these genes are responsible for cellular defense or apoptosis through expressing either anti- or pro-apoptotic proteins [24,25,26,27]. Some of these genes are responsible for cellular defense or apoptosis through expressing either anti- or pro-apoptotic proteins. Closer inspection of the genes associated with this pathway after MPP+ treatment showed that two genes encoding receptor tyrosine kinase (RTK), such as MELK and SHMT2, were differentially expressed (Figure 4). These genes are located very upstream in the EGFR pathway. They affect the entry points controlling this pathway. Therefore, their alteration can be widely propagated throughout the pathway. Also, a previous report indicated increased HSPB1 in PD under conditions of ischemic or thermal stress [28]. HSPB1, which acts through caspase inactivation signaling mechanisms, is known to be neuroprotective and anti-apoptotic [29]. Thus, the increased in expression observed in MPP+ induced apoptosis may represent another compensatory response.
Among DEGs of the EGFR related pathway, the connexin proteins are encoded by a multi-gene family, and so far 23 different human Cx genes have been identified [30,31]. The modulation of connexins affected the cell viability or growth, implying that connexins may have an important role in maintaining homeostasis in various organs [32,33]. The connexin is very crucial for the homeostatic regulation in the brain cells and helps with the communication of various electrical and chemical signals for maintaining the neuronal activities [34,35]. Alteration in Cx protein expression may lead to various congenital diseases including hearing loss, which was linked to the Cx protein family [36]. The connexins were regulated by the EGF signaling pathway and apoptosis [37]. Among Cx genes, we have found a Cx26 gene that markedly increased during the MPP+ treatment time (Table 2). The alteration in the expression level of Cx26 during MPP+-induced apoptosis and whether a similar response was shown in H2O2- and 6-OHDA-treated SH-SY5Y cells was investigated (Figure 5). Further, we examined the expression pattern of Cx26 in MPTP-intoxicated mice. The Cx26 protein level was significantly increased after MPTP injection, compared with that of saline controls. It is interesting to note that Cx26 that is expressed to colorectal cancer may correlate with Bcl-xL and Bax [38]. Considering the data obtained, our findings suggest that Cx26 might play significant role in DAergic neuronal cell death during the process of neuro-apoptosis. EGFR signaling pathways might also be responsible for DAergic neuronal cell death and therefore can be focused on as potential targets for therapeutic intervention.

4. Materials and Methods

4.1. Reagents

MPP+, MPTP, and 5-diphenyl-tetrazolium bromide (MTT) were obtained from Sigma Aldrich (St. Louis, MO, USA). Six-well and 100 mm culture dishes were purchased from Nunc Inc. (North Aurora Road, IL, USA). Dulbecco’s modified Eagle’s medium (DMEM), fetal bovine serum (FBS), and other cell culture reagents was purchased from Gibco-BRL Technologies (Rockville, CA, USA). Antibodies against cleaved PARP and β-actin were obtained from Cell Signaling Co., (Boston, MA, USA). All other chemicals used in this study were of analytical grade and obtained from Sigma Chemical Co. (St. Louis, MO, USA) unless otherwise noted.

4.2. Cell Culture and Treatments

SH-SY5Y neuroblastoma cells were obtained as described previously [39]. Briefly, cells were grown and maintained in DMEM containing 10% heat-inactivated fetal bovine serum (FBS) and 100 μg/mL penicillin-streptomycin and maintained in a humidified incubator at 37 °C with 5% CO2. Following the induction of differentiation with 10 µM all-trans retinoic acid (Sigma-Aldrich, St. Louis, MO, USA) for 7 days prior to treatment, 1 mM MPP+, 100 µM H2O2, and 25 µM 6-OHDA at final concentration were used to treat the SH-SY5Y cells. The total RNAs and proteins were harvested from these cells treated with MPP+ for 0, 6, 12, 24, and 48 h, respectively.

4.3. Assessment of Cell Viability

Cell viability was measured by using a previously described quantitative colorimetric MTT assay and LDH assay [39]. Briefly, cells seeded at 2.5 × 105 cells/mL in a 96-well plate were treated with 1 mM MPP+ for 48 h. The supernatants were removed from each well and the formazan crystals in viable cells were dissolved in DMSO. Optical density was measured at 550 nm using a microplate reader (Molecular Devices, Sunnyvale, CA, USA). Each cell suspension was centrifuged (4000× g, 5 min, 4 °C) (Mega 17R, Hanil Scientific Inc., Kimpo-si, Korea) and then the supernatant was collected. LDH activity in the supernatant was performed by using a cytotoxicity assay kit according to the manufacturer’s instructions (Takara Bio, Shiga, Japan). Absorbance was read at 440 nm and cytotoxicity (%) was calculated as:
{(supernatant value − blank value)/[(supernatant value − blank value) + (upper control value − blank value)]} × 100.
Each experiment was performed in triplicate.

4.4. Animals and Treatments

The animals were cared for as described previously [40]. Male C57BL/6 mice (age, 8 weeks; weight, 25–30 g) were obtained from Samtako Bio Korea (Gyeonggi-do, Korea) and were used in the present study. All experiments were performed in accordance with the Principles of Laboratory Animal Care (NIH publication No. 85-23, revised 1985) and the experimental procedures were approved by the Institutional Animal Care and Use Committee of Konkuk University (No. KU16171, 6 May 2016). The animals were housed in a controlled environment (23 ± 1 °C and 50% ± 5% humidity) and allowed food and water ad libitum. The room lights were on between 8:00 a.m. and 20:00 p.m. 25 animals were taken and divided into five groups containing five animals each. MPTP sub-chronic (30 mg/kg/day for five days) was administerd to the animals. MPTP was dissolved in saline and prepared just prior to dosing. In the sub-chronic group, the animals were sacrificed at 1, 2, 4, and 7 days to characterize the expression pattern of Cx26.

4.5. Microarray Analysis

The experiment was performed using a Human Transcriptome Pico Assay 2.0 (Agilent Technologies, Santa Clara, CA, USA). Microarray analyses were performed as described previously [41]. Total RNA was isolated from cells using TRIzol (Invitrogen, Carlsbad, CA, USA). Extracted total RNA was inspected with a Bioanalyzer 2100 (Agilent Technologies). Microarray hybridization was done according to the manufacturer’s instructions. Fluorescence-labeled cDNA probes were prepared from 5 μg of total RNA by RT primer (Genisphere, Hatfield, PA, USA)-primed polymerization using Super-Script II reverse transcriptase (Invitrogen) in a total reaction volume of 10.5 mL. After reverse transcription, the sample RNA was degraded by adding 1 mL of stop solution (0.5 M NaOH/50 mM EDTA) and incubating 65 μL for 10 min. After two labeled cDNAs were mixed, the mixture was denaturized at 95 °C for 2 min and incubated in a 45 °C water chamber for 20 min. The cDNA mixture was then placed on the cDNA chip and covered by a hybridization chamber. The slides were hybridized in a 62 °C hybridization oven for 12 h. Hybridized slides were scanned with a GenePix 4200 A scanner (Molecular Devices, Sunnyvale, CA, USA). Scanned images were analyzed using GeneSpring 7.3.1 software (Agilent Technologies). To allow algorithmic elimination of bad spots, no data points were eliminated by visual inspection from the initial GenePix image. For signal normalization, positive control genes were spotted onto each slide. The signals of these spots were used for normalization. To determine the background signal intensity, spotting solution was spotted onto each slide. To filter out unreliable data, spots with a signal-to-noise ratio (signal background-background, standard deviation) below 100 were not included in the data. Data were normalized by global, lowess, print-tip, and scaled normalization for data reliability. Genes of interest were those that exhibited a 2-fold change between test and control samples. Data were clustered in groups of genes that behaved similarly across the experiment using GeneSpring 7.3.1 (Agilent Technologies). An algorithm based on Euclidean distance was utilized to separate the genes with similar patterns. The distance cutoff was considered as statistically significant when there was a 2-fold change between time course experiments. The correlation cutoff was 0.95.

4.6. Expression Analysis

SH-SY5Y cells (1 × 106 cells/mL) were cultured in six well plates. Total RNA was extracted with TRIzol (Invitrogen). For the reverse transcription-polymerase chain reaction (RT-PCR), 2.5 µg of total RNA from the SH-SY5Y cells was used using First Strand cDNA synthesis kit (Invitrogen). RT-PCR was performed using the cDNA as a template. Nucleotide sequences of the primers used for RT-PCR are shown in Table 3. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as an internal control to evaluate the relative expressions of specific genes.
The protein expression analyses were performed as described previously [39]. Aliquot (0.1 mL) of RIPA buffer (1× PBS, 1% NP-40, 0.5% sodium deoxycholate, 0.1% SDS containing freshly added protease inhibitor cocktail (Calbiochem, San Diego, CA, USA) was added to the cells cultured in 100 mm plates to obtain the total cell lysate. Cells were scraped, incubated for 10 min on ice, and centrifuged at 14,000 × rpm for 10 min at 4 °C. For the animal experiment, tissues were washed two times with PBS, placed at 4 °C, and homogenized using a 1 mL syringe in lysis buffer (1x RIPA lysis buffer, protease inhibitor cocktail, and phosphatase inhibitor cocktail) and then finally passed through a 311/2 gauge syringe needle and centrifuged at 14,000 rpm at 4 °C for 15 min. The supernatants were collected for further analysis. Protein concentration was determined by the Bio-Rad DC Protein Assay (Hercules, CA, USA). 15 µg of whole cell lysates (cell experiments), or 40 µg of protein (animal experiments) were separated electrophoretically by 10% sodium dodecyl sulfate-polyacrylamide electrophoresis (SDS-PAGE), and the resolved proteins were transferred to polyvinylidene difluoride membranes (Millipore, Bedford, MA, USA). The membranes were incubated for 1 h with 5% skim milk in PBS buffer to block non-specific binding and then incubated with primary antibodies to anti-cleaved PARP, anti-Cx26 (1:1000; Cell Signaling Technology, Boston, MA, USA), and anti-β-actin (1:2000; Cell Signaling Technology). The blots were visualized using the SuperSignal West Pico Chemiluminescent Substrate Detection System (Thermo Scientific, Rockford, IL, USA), according to the manufacturer’s procedure. The optical densities of the antibody-specific bands were analyzed using a luminescent image analyzer, LAS-3000 (Fuji, Tokyo, Japan).

4.7. Statistical Analysis

All measurements were performed in triplicates and repeated three times using different samples. Results are presented as the mean ± S.E.M., (n = 3). The level of statistical significance was determined by analysis of variance (ANOVA) followed by Dunnett’s t-test for multiple comparisons. p-values of less than 0.05 were considered statistically significant.

5. Conclusions

SH-SY5Yneuroblastoma cells show great potential as a useful research model for research on various CNS diseases such as neurodegenerative diseases, in which DAergic neuronal injury is prominent in the pathophysiology. The DAergic neuronal injury model is particularly necessary for studying the pathogenesis of Parkinson’s disease (PD). We investigated a global pattern of gene expression analysis in MPP+ treated SH-SY5Y cells. DEGs were sorted into lists based on the direction of regulation, and corresponding GO categories were identified. The GO term of MPP+ induced apoptosis was common in most Signal transductions. Within signal transduction, we investigated the expression of EGFR pathway-related genes. Moreover, we have shown that expression of Cx26 was increased with DAergic neuronal cell death. Therefore EGFR pathway-related genes and Cx26 associated genes might be considered to play an important role, and their regulation may mitigate the neuronal cell death seen in neurodegenerative disorders including PD.

Acknowledgments

This work was funded by Konkuk University 2016.

Author Contributions

In-Su Kim and Dong-Kug Choi conceived and designed the experiment. In-Su Kim, Sushruta Koppula, and Shin-Young Park performed experiments and analyzed the data. In-Su Kim, Sushruta Koppula, and Dong-Kug Choi wrote the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pieczenik, S.R.; Neustadt, J. Mitochondrial dysfunction and molecular pathways of disease. Exp. Mol. Pathol. 2007, 83, 84–92. [Google Scholar] [CrossRef] [PubMed]
  2. Rodriguez-Arribas, M.; Yakhine-Diop, S.M.; Pedro, J.M.; Gomez-Suaga, P.; Gomez-Sanchez, R.; Martinez-Chacon, G.; Fuentes, J.M.; Gonzalez-Polo, R.A.; Niso-Santano, M. Mitochondria-associated membranes (MAMs): Overview and its role in Parkinson’s disease. Mol. Neurobiol. 2016. [Google Scholar] [CrossRef] [PubMed]
  3. Singer, T.P.; Ramsay, R.R. Mechanism of the neurotoxicity of MPTP. An update. FEBS Lett. 1990, 274, 1–8. [Google Scholar] [PubMed]
  4. Eberhardt, O.; Schulz, J.B. Apoptotic mechanisms and antiapoptotic therapy in the MPTP model of Parkinson’s disease. Toxicol. Lett. 2003, 139, 135–151. [Google Scholar] [CrossRef]
  5. Przedborski, S.; Tieu, K.; Perier, C.; Vila, M. MPTP as a mitochondrial neurotoxic model of Parkinson’s disease. J. Bioenerg. Biomembr. 2004, 36, 375–379. [Google Scholar] [CrossRef] [PubMed]
  6. Blesa, J.; Przedborski, S. Parkinson’s disease: Animal models and dopaminergic cell vulnerability. Front. Neuroanat. 2014, 8, 155. [Google Scholar] [CrossRef] [PubMed]
  7. Fiskum, G.; Starkov, A.; Polster, B.M.; Chinopoulos, C. Mitochondrial mechanisms of neural cell death and neuroprotective interventions in Parkinson’s disease. Ann. N. Y. Acad. Sci. 2003, 991, 111–119. [Google Scholar] [CrossRef] [PubMed]
  8. Chaturvedi, R.K.; Flint Beal, M. Mitochondrial diseases of the brain. Free Radic. Biol. Med. 2013, 63, 1–29. [Google Scholar] [CrossRef] [PubMed]
  9. Conn, K.J.; Ullman, M.D.; Larned, M.J.; Eisenhauer, P.B.; Fine, R.E.; Wells, J.M. cDNA microarray analysis of changes in gene expression associated with MPP+ toxicity in SH-SY5Y cells. Neurochem. Res. 2003, 28, 1873–1881. [Google Scholar] [CrossRef] [PubMed]
  10. Datki, Z.; Juhasz, A.; Galfi, M.; Soos, K.; Papp, R.; Zadori, D.; Penke, B. Method for measuring neurotoxicity of aggregating polypeptides with the MTT assay on differentiated neuroblastoma cells. Brain Res. Bull. 2003, 62, 223–229. [Google Scholar] [CrossRef] [PubMed]
  11. Cory, S.; Adams, J.M. The BCL2 family: Regulators of the cellular life-or-death switch. Nat. Rev. Cancer 2002, 2, 647–656. [Google Scholar] [CrossRef] [PubMed]
  12. Martire, S.; Mosca, L.; d’Erme, M. PARP-1 involvement in neurodegeneration: A focus on Alzheimer’s and Parkinson’s diseases. Mech. Ageing Dev. 2015, 146–148, 53–64. [Google Scholar] [CrossRef] [PubMed]
  13. Perez, J.G.; Tran, N.L.; Rosenblum, M.G.; Schneider, C.S.; Connolly, N.P.; Kim, A.J.; Woodworth, G.F.; Winkles, J.A. The tweak receptor FN14 is a potential cell surface portal for targeted delivery of glioblastoma therapeutics. Oncogene 2016, 35, 2145–2155. [Google Scholar] [CrossRef] [PubMed]
  14. Crouch, P.J.; Harding, S.M.; White, A.R.; Camakaris, J.; Bush, A.I.; Masters, C.L. Mechanisms of a β mediated neurodegeneration in Alzheimer’s disease. Int. J. Biochem. Cell. Biol. 2008, 40, 181–198. [Google Scholar] [CrossRef] [PubMed]
  15. Hoe, H.S.; Lee, H.K.; Pak, D.T. The upside of APP at synapses. CNS Neurosci. Ther. 2012, 18, 47–56. [Google Scholar] [CrossRef] [PubMed]
  16. Lu, J.; Esposito, G.; Scuderi, C.; Steardo, L.; Delli-Bovi, L.C.; Hecht, J.L.; Dickinson, B.C.; Chang, C.J.; Mori, T.; Sheen, V. S100b and APP promote a gliocentric shift and impaired neurogenesis in down syndrome neural progenitors. PLoS ONE 2011, 6, e22126. [Google Scholar] [CrossRef] [PubMed]
  17. Lin, N.Y.; Lin, T.Y.; Yang, W.H.; Wang, S.C.; Wang, K.T.; Su, Y.L.; Jiang, Y.W.; Chang, G.D.; Chang, C.J. Differential expression and functional analysis of the tristetraprolin family during early differentiation of 3T3-L1 preadipocytes. Int. J. Biol. Sci. 2012, 8, 761–777. [Google Scholar] [CrossRef] [PubMed]
  18. Chen, S.F.; Nieh, S.; Jao, S.W.; Liu, C.L.; Wu, C.H.; Chang, Y.C.; Yang, C.Y.; Lin, Y.S. Quercetin suppresses drug-resistant spheres via the P38 MAPK-HSP27 apoptotic pathway in oral cancer cells. PLoS ONE 2012, 7, e49275. [Google Scholar] [CrossRef] [PubMed]
  19. Leibowitz, B.J.; Cohick, W.S. Endogenous IGFBP-3 is required for both growth factor-stimulated cell proliferation and cytokine-induced apoptosis in mammary epithelial cells. J. Cell. Physiol. 2009, 220, 182–188. [Google Scholar] [CrossRef] [PubMed]
  20. Rash, J.E.; Yasumura, T.; Davidson, K.G.; Furman, C.S.; Dudek, F.E.; Nagy, J.I. Identification of cells expressing Cx43, Cx30, Cx26, Cx32 and Cx36 in gap junctions of rat brain and spinal cord. Cell. Commun. Adhes. 2001, 8, 315–320. [Google Scholar] [CrossRef] [PubMed]
  21. Halliwell, B. Reactive oxygen species and the central nervous system. J. Neurochem. 1992, 59, 1609–1623. [Google Scholar] [CrossRef] [PubMed]
  22. Hirsch, E.C. Does oxidative stress participate in nerve cell death in Parkinson’s disease? Eur. Neurol. 1993, 33, 52–59. [Google Scholar] [CrossRef] [PubMed]
  23. Jackson-Lewis, V.; Przedborski, S. Protocol for the MPTP mouse model of Parkinson’s disease. Nat. Protoc. 2007, 2, 141–151. [Google Scholar] [CrossRef] [PubMed]
  24. Riederer, P.; Jellinger, K. Morphological and biochemical changes in the aging brain: Pathophysiological and possible therapeutic consequences. Exp. Brain Res. 1982, 5, 158–166. [Google Scholar] [PubMed]
  25. Jorissen, R.N.; Walker, F.; Pouliot, N.; Garrett, T.P.; Ward, C.W.; Burgess, A.W. Epidermal growth factor receptor: Mechanisms of activation and signalling. Exp. Cell. Res. 2003, 284, 31–53. [Google Scholar] [CrossRef]
  26. Van der Geer, P.; Hunter, T.; Lindberg, R.A. Receptor protein-tyrosine kinases and their signal transduction pathways. Annu. Rev. Cell. Biol. 1994, 10, 251–337. [Google Scholar] [CrossRef] [PubMed]
  27. Daub, H.; Weiss, F.U.; Wallasch, C.; Ullrich, A. Role of transactivation of the EGF receptor in signalling by G-protein-coupled receptors. Nature 1996, 379, 557–560. [Google Scholar] [CrossRef] [PubMed]
  28. Kato, H.; Kogure, K.; Liu, X.H.; Araki, T.; Kato, K.; Itoyama, Y. Immunohistochemical localization of the low molecular weight stress protein HSP27 following focal cerebral ischemia in the rat. Brain Res. 1995, 679, 1–7. [Google Scholar] [CrossRef]
  29. Bruey, J.M.; Ducasse, C.; Bonniaud, P.; Ravagnan, L.; Susin, S.A.; Diaz-Latoud, C.; Gurbuxani, S.; Arrigo, A.P.; Kroemer, G.; Solary, E.; et al. HSP27 negatively regulates cell death by interacting with cytochrome C. Nat. Cell. Biol. 2000, 2, 645–652. [Google Scholar] [PubMed]
  30. Belluardo, N.; White, T.W.; Srinivas, M.; Trovato-Salinaro, A.; Ripps, H.; Mudo, G.; Bruzzone, R.; Condorelli, D.F. Identification and functional expression of HCX31.9, a novel gap junction gene. Cell. Commun. Adhes. 2001, 8, 173–178. [Google Scholar] [CrossRef] [PubMed]
  31. Willecke, K.; Eiberger, J.; Degen, J.; Eckardt, D.; Romualdi, A.; Guldenagel, M.; Deutsch, U.; Sohl, G. Structural and functional diversity of connexin genes in the mouse and human genome. Biol. Chem. 2002, 383, 725–737. [Google Scholar] [CrossRef] [PubMed]
  32. Hayashi, T.; Matesic, D.F.; Nomata, K.; Kang, K.S.; Chang, C.C.; Trosko, J.E. Stimulation of cell proliferation and inhibition of gap junctional intercellular communication by linoleic acid. Cancer Lett. 1997, 112, 103–111. [Google Scholar] [CrossRef]
  33. Ogawa, T.; Hayashi, T.; Tokunou, M.; Nakachi, K.; Trosko, J.E.; Chang, C.C.; Yorioka, N. Suberoylanilide hydroxamic acid enhances gap junctional intercellular communication via acetylation of histone containing connexin 43 gene locus. Cancer Res. 2005, 65, 9771–9778. [Google Scholar] [CrossRef] [PubMed]
  34. Anderson, C.M.; Swanson, R.A. Astrocyte glutamate transport: Review of properties, regulation, and physiological functions. Glia 2000, 32, 1–14. [Google Scholar] [CrossRef]
  35. Ransom, B.R.; Brown, A.M. Intracellular Ca2+ release and ischemic axon injury: The Trojan horse is back. Neuron 2003, 40, 2–4. [Google Scholar] [CrossRef]
  36. Kelly, J.J.; Simek, J.; Laird, D.W. Mechanisms linking connexin mutations to human diseases. Cell. Tissue Res. 2015, 360, 701–721. [Google Scholar] [CrossRef] [PubMed]
  37. Dube, E.; Dufresne, J.; Chan, P.T.; Cyr, D.G. Epidermal growth factor regulates connexin 43 in the human epididymis: Role of gap junctions in azoospermia. Hum. Reprod. 2012, 27, 2285–2296. [Google Scholar] [CrossRef] [PubMed]
  38. Kanczuga-Koda, L.; Sulkowski, S.; Koda, M.; Skrzydlewska, E.; Sulkowska, M. Connexin 26 correlates with Bcl-xl and Bax proteins expression in colorectal cancer. World J. Gastroenterol. 2005, 11, 1544–1548. [Google Scholar] [CrossRef] [PubMed]
  39. Kim, I.S.; Ko, H.M.; Koppula, S.; Kim, B.W.; Choi, D.K. Protective effect of Chrysanthemum indicum Linne against 1-methyl-4-phenylpridinium ion and lipopolysaccharide-induced cytotoxicity in cellular model of Parkinson’s disease. Food Chem. Toxicol. 2009, 49, 963–973. [Google Scholar] [CrossRef] [PubMed]
  40. Kim, B.W.; Koppula, S.; Kumar, H.; Park, J.Y.; Kim, I.W.; More, S.V.; Kim, I.S.; Han, S.D.; Kim, S.K.; Yoon, S.H.; Choi, D.K. α-Asarone attenuates microglia-mediated neuroinflammation by inhibiting NF κB activation and mitigates MPTP-induced behavioral deficits in a mouse model of Parkinson’s disease. Neuropharmacology 2015, 97, 46–57. [Google Scholar] [CrossRef] [PubMed]
  41. Kim, I.S.; Park, J.H.; Do, J.H.; Ko, H.M.; Park, E.J.; Suk, K.; Choi, D.K. Expression profiling of neuroinflammatory genes in LPS and LPS/IFNγ-stimulated BV-2 microglia cells. Biochip J. 2009, 3, 299–305. [Google Scholar]
Figure 1. Cell viability (A,B,D) and biochemical markers of apoptosis (C,E,F) in SH-SY5Y cells after treatment for up to 48 h. Cell viability was assessed by the MTT method (A,B) and LDH assay (D) presented as a percentage of untreated controls. The biochemical markers were analyzed by 1mM 1-methyl-4-phenylpyridine (MPP+) induced apoptosis after treatment for up to 48 h. (A) SH-SY5Y cells were treated with different concentrations (0.5, 1, 2 mM) of MPP+ for 48 h; (B,D) SH-SY5Y cells were incubated for 48 h with 1 mM MPP+; (C) The biochemical markers of apoptosis were shown by expression analysis; (E) The Bax/Bcl-2 ratio was determined; (F) The levels of cleaved PARP were quantified by densitometric analysis. Data are expressed as a percentage of values compared to untreated control cultures. Values are mean ± standard error (* p < 0.05, ** p < 0.01 and *** p < 0.005 vs. control group).
Figure 1. Cell viability (A,B,D) and biochemical markers of apoptosis (C,E,F) in SH-SY5Y cells after treatment for up to 48 h. Cell viability was assessed by the MTT method (A,B) and LDH assay (D) presented as a percentage of untreated controls. The biochemical markers were analyzed by 1mM 1-methyl-4-phenylpyridine (MPP+) induced apoptosis after treatment for up to 48 h. (A) SH-SY5Y cells were treated with different concentrations (0.5, 1, 2 mM) of MPP+ for 48 h; (B,D) SH-SY5Y cells were incubated for 48 h with 1 mM MPP+; (C) The biochemical markers of apoptosis were shown by expression analysis; (E) The Bax/Bcl-2 ratio was determined; (F) The levels of cleaved PARP were quantified by densitometric analysis. Data are expressed as a percentage of values compared to untreated control cultures. Values are mean ± standard error (* p < 0.05, ** p < 0.01 and *** p < 0.005 vs. control group).
Ijms 18 00430 g001aIjms 18 00430 g001b
Figure 2. Gene Ontology (GO) biological process classification of differentially expressed genes (DEGs) induced by MPP+ treatment. (A) Down-regulated genes in SH-SY5Y cells treated by MPP+ for 24 h; (B) Down-regulated genes in SH-SY5Y cells treated by MPP+ for 48 h; (C) Up-regulated genes in SH-SY5Y cells treated by MPP+ for 24 h; (D) Up-regulated genes in SH-SY5Y cells treated by MPP+ for 48 h.
Figure 2. Gene Ontology (GO) biological process classification of differentially expressed genes (DEGs) induced by MPP+ treatment. (A) Down-regulated genes in SH-SY5Y cells treated by MPP+ for 24 h; (B) Down-regulated genes in SH-SY5Y cells treated by MPP+ for 48 h; (C) Up-regulated genes in SH-SY5Y cells treated by MPP+ for 24 h; (D) Up-regulated genes in SH-SY5Y cells treated by MPP+ for 48 h.
Ijms 18 00430 g002
Figure 3. Validation of microarray results. Cells were treated with 1 mM of MPP+ for 6 to 48 h. (A) Total RNA was extracted and converted to cDNA followed by RT-PCR; (B) Bar graphs show microarray results for selected DEGs; (C) The levels of mRNA expression were quantified by densitometric analysis. Values are mean ± standard error.
Figure 3. Validation of microarray results. Cells were treated with 1 mM of MPP+ for 6 to 48 h. (A) Total RNA was extracted and converted to cDNA followed by RT-PCR; (B) Bar graphs show microarray results for selected DEGs; (C) The levels of mRNA expression were quantified by densitometric analysis. Values are mean ± standard error.
Ijms 18 00430 g003
Figure 4. Expression analysis of EGFR related genes in MPP+ induced apoptosis. Cells were treated with 1 mM MPP+ for 6 to 48 h; (A) Total RNA was extracted and converted to cDNA followed by RT-PCR; (B) Bar graphs show microarray results for selected DEG; (C) The levels of mRNA expression were quantified by densitometric analysis. Values are mean ± standard error.
Figure 4. Expression analysis of EGFR related genes in MPP+ induced apoptosis. Cells were treated with 1 mM MPP+ for 6 to 48 h; (A) Total RNA was extracted and converted to cDNA followed by RT-PCR; (B) Bar graphs show microarray results for selected DEG; (C) The levels of mRNA expression were quantified by densitometric analysis. Values are mean ± standard error.
Ijms 18 00430 g004
Figure 5. Expression of Cx26 on toxin treated SH-SY5Y cells. SH-SY5Y cells were treated with (A) 1 mM MPP+, (B) 100 µM H2O2, and (C) 25 µM 6-OHDA for 6 to 48 h. Cell lysate samples (20 µg protein/lane) were measured by immunoblot analysis using Cx26 antibody. (D) The levels of CX26 protein were quantified by densitometric analysis. Values are mean ± standard error (*** p < 0.01 vs. control group).
Figure 5. Expression of Cx26 on toxin treated SH-SY5Y cells. SH-SY5Y cells were treated with (A) 1 mM MPP+, (B) 100 µM H2O2, and (C) 25 µM 6-OHDA for 6 to 48 h. Cell lysate samples (20 µg protein/lane) were measured by immunoblot analysis using Cx26 antibody. (D) The levels of CX26 protein were quantified by densitometric analysis. Values are mean ± standard error (*** p < 0.01 vs. control group).
Ijms 18 00430 g005
Figure 6. Immunoblot analysis of Cx26 in the midbrain and striatum of a C57BL/6 mouse after sub-chronic MPTP injection. Expression pattern of Cx26 was performed after sub-chronic (A) injection of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) in C57BL/6 mice. Bar graphs show quantitative data for Cx26 signals that are normalized to the β-actin signal (n = 3–4 per group); (B) midbrain (C) striatum. Values are mean ± standard error (* p < 0.05 vs. control group).
Figure 6. Immunoblot analysis of Cx26 in the midbrain and striatum of a C57BL/6 mouse after sub-chronic MPTP injection. Expression pattern of Cx26 was performed after sub-chronic (A) injection of 1-methyl-4-phenyl-1,2,3,6-tetrahydropyridine (MPTP) in C57BL/6 mice. Bar graphs show quantitative data for Cx26 signals that are normalized to the β-actin signal (n = 3–4 per group); (B) midbrain (C) striatum. Values are mean ± standard error (* p < 0.05 vs. control group).
Ijms 18 00430 g006
Table 1. List of genes with significant (p < 0.05) mRNA level changes in MPP+ induced apoptosis.
Table 1. List of genes with significant (p < 0.05) mRNA level changes in MPP+ induced apoptosis.
GenBank Accession No.Gene SymbolDescriptionFold ChangeKEGG Pathwaysp-Value
24 h48 h
NM_001775CD38ADP-ribosyl cyclase 10.490.47Calcium signaling pathway6.48 × 10−4
NM_000957PTGER3Prostaglandin E2 receptor, EP3 subtype0.480.441.58- × 10−4
NM_000833GRIN2AGlutamate receptor subunit epsilon 1 precursor0.490.092.62 × 105
NM_004160PYYPeptide YY precursor0.490.50Neuroactive ligand-receptor interaction2.36 × 10−4
BC045651P2RY5P2Y purinoceptor 50.430.471.54 × 10−4
BC036030GABRG2Gamma-aminobutyric-acid receptor γ-2 subunit precursor0.460.435.56 × 105
AY358893CCL26Small inducible cytokine A26 precursor0.450.49Cytokine-cytokine receptor interaction4.66 × 105
NM_003811TNFSF9Tumor necrosis factor ligand superfamily member 90.440.492.05 × 10−4
BM555452TNFRSF12ATumor necrosis factor receptor superfamily member Fn14 precursor4.330.988.36 × 105
X06374PDGFASplice Isoform Short of Platelet-derived growth factor, A chain precursor0.690.222.43 × 10−4
NM_000565IL6RInterleukin-6 receptor α chain precursor0.440.422.23 × 10−4
NM_006917RXRGRetinoic acid receptor RXR-γ0.460.44Adipocytokine signaling pathway1.63 × 10−4
NM_004457ACSL3Long-chain-fatty-acid--CoA ligase 30.460.063.37 × 10−4
BC030529PSMA4Proteasome subunit α type 42.256.75Proteasome8.35 × 10−4
NM_002790PSMA5Proteasome subunit α type 50.450.474.39 × 10−4
AF097935DSG1Desmoglein-1 precursor0.460.44Cell Communication2.58 × 108
BC034761ADMRAdrenomedullin receptor0.840.11G-protein coupled receptor protein signaling pathway5.33 × 10−4
NM_012202GNG3Guanine nucleotide-binding protein G(I)/G(S)/G(O) γ-3 subunit.0.730.473.68 × 10−4
BC055286A4GALTLactosylceramide 4-α-galactosyltransferase0.430.45Globoside metabolism4.71 × 10−4
AY643499HEXBβ-Hexosaminidase β chain precursor0.430.498.27 × 10−4
CR627398MAML2Mastermind-like 2 (Drosophila)5.902.96Notch signaling pathway6.79 × 10−4
BC018937APPAmyloid β A4 protein precursor2.211.272.11 × 10−4
AK057791METTL6Hypothetical protein MGC241324.602.07Histidine metabolism4.55 × 106
NM_020997LEFTY1Left-right determination factor B precursor0.250.30TGF-β signaling pathway8.39 × 10−4
L04751CYP4A11Cytochrome P450 4A11 precursor0.480.45Fatty acid metabolism6.27 × 10−4
BC073828PPM1JProtein phosphatase 2a, catalytic subunit, isoform0.460.45Tight junction2.61 × 10−4
NM_001015049BAG5B-cell CLL/lymphoma 7C0.490.43Anti-apoptosis7.38 × 10−4
NM_004323BAG1BCL-2 binding athanogene-10.480.326.83 × 10−4
NM_000633BCL2B-cell CLL/lymphoma 20.480.212.1 × 106
AK001497DEDDDeath effector domain-containing protein0.540.47Induction of apoptosis7.85 × 10−4
NM_031882PCDHAC1Protocadherin α subfamily C, 10.850.23Cell adhesion3.67 × 105
CR622836NALP10NACHT-, LRR- and PYD-containing protein 100.812.13Unknown8.90 × 106
AK075382ATP6AP2Vacuolar ATP synthase membrane sector associated protein M8-90.880.263.47 × 10−4
NM_016630SPG21Spastic paraplegia 211.813.881.24 × 104
BC039110SMNDC1Survival of motor neuron-related splicing factor 3022.71.54 × 105
Table 2. List of genes identified in EGFR signaling pathway related genes.
Table 2. List of genes identified in EGFR signaling pathway related genes.
GenBank Accession No.Gene SymbolDescriptionFold Changep-Value
24 h48 h
BM907768HSPB1Heat-shock protein β-10.48 ± 0.090.58 ± 0.163.47 × 104
NM_001013398IGFBP3Insulin-like growth factor binding protein 3 precursor0.49 ± 0.080.64 ± 0.141.54 × 105
NM_014791MELKMaternal embryonic leucine zipper kinase1.96 ± 0.242.15 ± 0.091.18 × 104
AL832618IFI44Linterferon induced protein 44-like1.76 ± 0.092.31 ± 0.331.20 × 104
AK055053SHMT2Serine hydroxymethyltransferase1.59 ± 0.122.12 ± 0.293.47 × 104
NM_004004GJB2Connexin 262.1 ± 0.309.5 ± 0.464.12 × 105
NM_006887ZFP36L2Butyrate response factor 20.95 ± 0.242.77 ± 0.424.51 × 104
NM_001753CAV1Caveolin-11.33 ± 0.062.38 ± 0.181.20 × 104
NM_003979GPRC5AG-protein coupled receptor family C group 5 member A0.42 ± 0.040.41 ± 0.023.66 × 104
NM_005558LAD1Ladinin 11.86 ± 0.143.7 ± 0.295.54 × 104
NM_181865ACOT7Cytosolicacyl coenzyme A thioester hydrolase1.36 ± 0.192.14 ± 0.161.20 × 104
AY358893SAA4Serum amyloid A-4 protein precursor2.09 ± 0.261.02 ± 0.164.62 × 104
Table 3. List of primers used for RT-PCR.
Table 3. List of primers used for RT-PCR.
Gene SymbolPrimer
ForwardReverse
Bcl-2ACGACTTCTCCCGCCGCTACCCCAGCCTCCGTTATCCTGG
BaxCACCAAGGTGCCGGAACTGAAATGCCCATGTCCCCCAATC
TNFSF9ACTCGGCCTTCGGTTTCCAGAAGGGGAGGTCAGCAGCAGG
ADMRGACCGCTATGTCACCCTCACGGTGCTGTACGTTTCAAAAGGT
APPAGCTTCTGGCCACTGGGAGGAGTGGGACATGCGGGGAAGT
TNFRSF12ACTCCTCCAACCACACAAGGGGTCCCCTCCAAACTCTCCCCA
PDGFACCAGCGACTCCTGGAGATAGACGTCCTGGTCTTGCAGACAG
PCDHAC1TGTGGGGTGGCAGTTTTATGTGAAAGGAAGCGAAAGTTCCG
NALP10CTACTTACGGGATATGACCCTGTGCTGGTCCACAAGTTCCAACA
GNG3CCGGTGAACAGCACTATGAGTGGCATCACAGTAAGTCATCAGG
ATP6AP2AGCTCCGTAATCGCCTGTTTCGATGCTTATGACGAGACAGCAA
SPG21TCTGACTGGATGGGGTTACCGGAAGCGCCAAAAAGATGAACTT
DEDDGCACCGCATGTTTGACATCGTCACGTCCATTTCGGATGAGT
BAG5AATGCAAACCACCCACACCGATCAGCCCGAGAGCACACA
HSPB1GTCCCTGGATGTCAACCACTCTTTACTTGGCGGCAGTCTC
IGFBP3CCAGCTCCAGGTGAGCCGGTCATGTCCTTGGCAGTCT
MELKCCATGTGCTAGAGACAGCCAATGCTACTGGGAGAGAGCCA
IFI44LTCCTAGCCATGTGTCCTTCCGGCCTACCTCCTCAATTTCC
SHMT2AGACTCAGACTGGGGAAGCATCAGTGCCAGCTGGTTGTAG
GJB2CTACTTCCCCATCTCCCACAGGGAAATGCTAGCGACTGAG
ZFP36L2CACCTTTCATACCATCGGCTGCAAAGTTGTGGGTCTGGAT
CAV1TGAGATTCAGTGCATCAGCCAACTTGAAATTGGCACCAGG
GPRC5ATCGTCGCAATGAAGACTTTGTGGTTCTGCAGCTGAAAATG
LAD1GACCATCTCCTTTCGGATGACCCTTGAGACGAAGACTTGC
ACOT7CGTCCTCAGAAAGGAAGTCGGTTCCTCCACTTGGTCTCCA
SAA4CAGCACAATGAGGCTTTTCAAGTGACCCTGTGTCCCTGTC
GAPDHGTCAGTGGTGGACCTGACCTTGTGAGGAGGGGAGATTCAG

Share and Cite

MDPI and ACS Style

Kim, I.-S.; Koppula, S.; Park, S.-Y.; Choi, D.-K. Analysis of Epidermal Growth Factor Receptor Related Gene Expression Changes in a Cellular and Animal Model of Parkinson’s Disease. Int. J. Mol. Sci. 2017, 18, 430. https://doi.org/10.3390/ijms18020430

AMA Style

Kim I-S, Koppula S, Park S-Y, Choi D-K. Analysis of Epidermal Growth Factor Receptor Related Gene Expression Changes in a Cellular and Animal Model of Parkinson’s Disease. International Journal of Molecular Sciences. 2017; 18(2):430. https://doi.org/10.3390/ijms18020430

Chicago/Turabian Style

Kim, In-Su, Sushruta Koppula, Shin-Young Park, and Dong-Kug Choi. 2017. "Analysis of Epidermal Growth Factor Receptor Related Gene Expression Changes in a Cellular and Animal Model of Parkinson’s Disease" International Journal of Molecular Sciences 18, no. 2: 430. https://doi.org/10.3390/ijms18020430

APA Style

Kim, I. -S., Koppula, S., Park, S. -Y., & Choi, D. -K. (2017). Analysis of Epidermal Growth Factor Receptor Related Gene Expression Changes in a Cellular and Animal Model of Parkinson’s Disease. International Journal of Molecular Sciences, 18(2), 430. https://doi.org/10.3390/ijms18020430

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop