Next Article in Journal
Constitutive Activation of Guanylate Cyclase by the G86R GCAP1 Variant Is Due to “Locking” Cation-π Interactions that Impair the Activator-to-Inhibitor Structural Transition
Next Article in Special Issue
Local and Systemic Humoral Response to Autologous Lineage-Negative Cells Intrathecal Administration in ALS Patients
Previous Article in Journal
Platelet Adhesion and Thrombus Formation in Microchannels: The Effect of Assay-Dependent Variables
Previous Article in Special Issue
An Update on the Progress of Isolation, Culture, Storage, and Clinical Application of Human Bone Marrow Mesenchymal Stem/Stromal Cells
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 21
Issue 3
10.3390/ijms21030751
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Bortezomib on Global Gene Expression in PC12-Derived Nerve Cells

by
Karolina Łuczkowska
,
Dorota Rogińska
,
Zofia Ulańczyk
and
Bogusław Machaliński
*
Department of General Pathology, Pomeranian Medical University, 70-111 Szczecin, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(3), 751; https://doi.org/10.3390/ijms21030751
Submission received: 30 December 2019 / Revised: 17 January 2020 / Accepted: 21 January 2020 / Published: 23 January 2020
(This article belongs to the Special Issue Paracrine Effects of Human Stem/Progenitor Cells—From Bench to Bedside)
Browse Figures
Review Reports Versions Notes

Abstract

:
Peripheral neuropathy is one of the main side-effects of novel therapeutics used in oncohematological diseases, but the molecular basis underlying its development and progression as well as neurotoxicity mechanisms induced by the use of these therapeutics are still not fully elucidated. The aim of this study was to demonstrate the effect of bortezomib on global gene and miRNA expression on PC12-derived nerve cells. Microarray analysis showed that expression of 1383 genes was downregulated at least two fold and 671 genes were upregulated at least two fold in PC12-derived nerve cells treated with bortezomib compared to untreated/control cells. Analysis of functional annotations mainly identified downregulated processes (e.g., regulation of cell cycle, DNA replication and repair, regulation of cell migration, neuron projection morphogenesis and neurotransmitter secretion). The result of miRNA expression analysis demonstrated only 11 significantly downregulated miRNAs (at least two fold) in bortezomib-treated PC12-derived nerve cells vs. control cells. MiRNAs regulate gene expression, therefore we decided to conduct an analysis comparing the outcomes of miRNA microarray expression data to the obtained mRNA data. The most interesting miRNA–target gene correlation is downregulated expression of miR-130a-3p and miR-152-3p and as a result of this downregulation the expression of the Gadd45 increased. This gene is a member of a group of genes, the transcript expression of which is enhanced after stressful growth arrest conditions and treatment with DNA-damaging agents like drugs or mutagens.

1. Introduction

Multiple myeloma (MM), which makes up around 10% of all haematological malignancies, is a malignant plasma cell disorder [1]. A major step forward in MM treatment has been the introduction of so-called novel agents, including the immunomodulatory drug (IMID) thalidomide, the subsequent development of its potent analogue lenalidomide, and the first-in-class proteasome inhibitor bortezomib (BZT). These agents serve a variety of anti-MM functions, including inhibition of MM cell proliferation and angiogenesis, apoptosis induction, and regulation of interactions between MM cells and bone marrow stromal cells [2]. Proteasome inhibitors are drugs that inhibit the ubiquitin-proteasome pathway and affect multiple signalling cascades in the cells. Bortezomib contains boron in its chemical structure and reversibly inhibits the chymotrypsin-like activity of the 26S proteasome degrading ubiquitinated protein, thus enabling the cell to maintain the balance of regulatory proteins involved in the cell cycle and division. Therefore, inhibition of proteasome activity by bortezomib leads to loss of cell cycle control and cell death [3]. BZT also decreases adhesion of myeloma plasma cells to stromal cells and disrupts secretion of specific cytokines in the bone marrow microenvironment such as vascular endothelial growth factor (VEGF), insulin-like growth factor 1 (IGF-1), interleukin-6 (IL-6), and tumor necrosis factor-α (TNF-α) [4]. The additional effect of bortezomib is based on the inhibition of DNA repair and angiogenesis, impairment of osteoclast activity [5] and disruption of intracellular calcium metabolism. This can lead to mitochondrial calcium inflow and apoptosis induced by caspase activation [6]. In addition, mitochondrial damage leads to oxidative stress through activation of ion channel TRPA1 in nociceptors.
From the early days of anti-MM agents’ clinical use, it has become apparent that peripheral neuropathy (PN) is one of their main non-haematological dose-limiting side-effects [7]. Bortezomib-induced peripheral neuropathy (BiPN) occurs in up to 40% of patients [8] and is related to the dose, schedule, and mode of administration [9]; it is usually reversible. Neurotoxicity can significantly affect the quality of a patient’s life, and may require dose-reduction, delay, or even premature termination of treatment. One of the hypothesised developments of neuropathy is the specific neurotoxic effect of bortezomib. BZT affects the transient release of intracellular calcium supply, leading to mitochondrial calcium influx and caspase-induced apoptosis. It has also been suggested that bortezomib exposure is closely associated with post-transcriptional gene silencing mediated regulation by non-coding RNAs, mainly microRNAs (miRNAs, miRs) [10,11]. MiRNAs are 20–22 nucleotide long molecules, which regulate gene expression and are involved in practically all biological processes, including apoptosis, cell proliferation, as well as disease-associated mechanisms [12]. According to the available literature, changes in miRNA expression are characteristic as a response to drug use in a number of diseases. [13]. However, the wide panel of mRNAs/miRNAs has not yet been tested in bortezomib-treated PC-12-derived nerve cells.
To better understand the underlying molecular mechanisms of bortezomib-induced neurotoxicity, in this study we analyzed the global expression of mRNA and miRNA in bortezomib-treated PC12-derived nerve cells and examined the changed biological pathways and miRNA–mRNA target interactions.

2. Results

2.1. Gene Expression Profile in PC12-derived Nerve Cells

Microarray analysis showed that expression of 1383 genes was downregulated at least two fold (fold -42.49 to -2) and 671 genes were upregulated at least two fold (fold 78.3 to 2) in PC12-derived nerve cells treated with bortezomib compared to untreated/control cells (Figure 1). The detailed fold values of most 15 up- or downregulated genes are presented in Table 1 and Table 2, respectively. Next, up/downregulated genes were classified according to the Gene Ontology (GO) classification of biological processes (Figure 2) and the KEGG Database of pathways (Figure 3). Analysis of functional annotations mainly identified downregulated processes (e.g., regulation of cell cycle, DNA replication and repair, regulation of cell migration, neuron projection morphogenesis, and neurotransmitter secretion). In addition, a detailed analysis of two pathways was performed (Figure 4: Axon guidance; Figure 5: Cell cycle).
Since treatment-induced neurotoxicity can be caused by a number of factors, in our analysis of microarray data we focused on several important processes that affect nervous system homeostasis. First, we focused on dysregulated processes essential for neural tissue function, such as neuron projection morphogenesis and neurotransmitter secretion, and genes involved in these processes that showed significant fold change in our study (Figure 6). Figure 7 presents more generalized processes that were also significantly dysregulated in our study (DNA repair, regulation of cell cycle, cell division), and genes involved in these processes that showed expression changes in our study.

Gene Set Enrichment Analysis

The Gene Set Enrichment Analysis (GSEA) method focuses on gene sets, that is, groups of genes that participate in a common biological function, chromosomal location, or regulation. Therefore, to more precisely investigate the biological influence of bortezomib on PC12-derived nerve cells, we carried out GSEA to identify altered pathways in bortezomib-treated PC12-derived nerve cells vs. control cells. According to the gene set enrichment analysis, 10 gene sets with a positive correlation were detected in bortezomib-treated cells and 10 gene sets with a negative correlation. Cell cycle checkpoints, DNA repair, and neuronal system-related genes were significantly downregulated in PC12-derived nerve cells treated with bortezomib for 24 h (false discovery rate (FDR) = 0.61%) (Figure 8).

2.2. Validation of Dysregulated mRNAs

We chose significantly dysregulated genes (Ascl1, Top2a, Slc7a11, Ddit3) to perform validation by qRT-PCR. In all selected genes, the results of qRT-PCR confirmed the result of expression of selected mRNAs obtained by global gene the expression analysis (Figure 9).

2.3. miRNAs Expression Profile in Neural Cells

The result of miRNA expression analysis demonstrated only 11 significantly downregulated miRNAs (at least two fold) in bortezomib-treated PC12-derived nerve vs. control cells (Figure 10).
MiRNAs regulate gene expression, therefore we decided to conduct an analysis comparing the outcomes of miRNA microarray expression data to the obtained mRNA data. Only this method shows the real impact of miRNAs on their target genes. Figure 11 shows changes in the expression of 9 miRNAs; target genes and changes in their expression in bortezomib-treated PC12-derived nerve compared to control cells. The identified miRNAs regulate various processes: Cell cycle, apoptotic process, regulation of cell growth.

2.4. Validation of Dysregulated miRNAs

We chose significantly dysregulated miRNAs (miRNA-21-5p; miRNA-322-5p; miRNA-532-5p) to perform validation by qRT-PCR. In all selected miRNAs, the results of qRT-PCR confirmed the result of expression of selected miRNAs obtained by global gene expression analysis (Figure 12).

3. Discussion

Since the incidence of chemotherapy-induced peripheral neuropathy in patients with MM treated with novel agents is still significant, the issue of the adverse effects of these drugs remains the important clinical problem. Suggested mechanisms underlying drug-induced PN comprise direct damage to dorsal root ganglion (DRG), mitochondrial and endoplasmic reticulum damage, deregulation of Ca2+ homeostasis, inhibition of neurotrophins’ transcription, microtubule stabilization, and autoimmune or inflammatory factors [14]. In this study, we demonstrated for the first time the global effect of bortezomib on expression of genes and miRNAs associated with neurogenesis, cell cycle, cell division, and DNA repair in PC12-derived nerve cells.
The introduction of bortezomib has been a significant breakthrough in the therapy of MM. BZT has a lot of different anti-myeloma effects including disruption of cell cycle, induction of apoptosis, alteration of the bone marrow microenvironment, and inhibition of nuclear factor kappa B (NFκB) [15]. The additional effect of bortezomib is based on the inhibition of DNA repair and angiogenesis, impairment of osteoclast activity [5], and disruption of intracellular calcium metabolism [16]. Recent studies have shown that bortezomib causes acute cytoskeletal damage, which is relevant in neural cell integrity. In addition, it has been demonstrated that bortezomib causes excessive protein carbonylation and destabilization of actin filaments [17]. Our current study provided interesting results that indicate a significant dysregulation of genes in neuronal cells, which may be characteristic for bortezomib activity (Figure 2, Figure 3 and Figure 8).
The DNA damage, for example altered bases and single- and double-strand breaks, is an important neurotoxic factor suggested in the aging and pathogenesis of many neurological disorders, including amyotrophic lateral sclerosis (ALS), Parkinson’s disease, and Alzheimer’s disease. DNA double-strand break (DSB) repair occurs in two ways, via: a) Homologous recombination, which uses a sister chromatid, and b) non-homologous end joining (NHEJ). DSB in the DNA cannot be completely repaired by the aforementioned repair ways in mature neurons because they are unable to actively divide. Therefore, it is considered that neurons accumulate DNA damage, which could potentially be the basis for the pathogenesis of neurodegenerative disorders [18]. DNA of peripheral nervous system (PNS) is more exposed to endo- and exogenous factors compared to the central nervous system (CNS) because the blood brain barrier (BBB) provides good protection [19]. One of the most downregulated genes in our study after bortezomib-treatment is topoisomerase (Top2a) (Table 1), which encodes an enzyme necessary for normal development. Topoisomerase relaxes coiled DNA to eliminate helical constraints, which may impede DNA replication and transcription and consequently inhibit cell growth [20]. Mutations in TDP1 play an important role in the development of neurodegenerative diseases such as ataxia telangiectasia (A-T) and spinocerebellar ataxia with axonal neuropathy (SCAN-1) [21]. We also observed reduced expression of another topoisomerase, topoisomerase (DNA) II binding protein 1 (Topbp1) (fold change = −7.37), which prevents DNA double-strand breaks. Lee et al. confirmed that TopBP1 is crucial for preventing replication-associated DNA strand breaks, but it is not necessary for replication [22]. In addition, they concluded that TopBP1 is essential for providing genome integrity in the early progenitors that direct neurogenesis. The basis for single-strand breaks (SSBs) repair is to generate DNA ends that enable ligation. DNA ligases play an important role in this process because they allow production of a phosphodiester bond to join broken DNA strands together. They play an important role in DNA replication, nuclear DNA metabolism, nucleotide excision repair (NER), single-strand break repair, homology-mediated DSB repair, and base excision repair (BER) [23]. In our study we observed decreased expression of Lig1 gene (fold = −3.605), which contribution to the development of neuropathy has been described [24]. The next gene with decreased expression (fold = −2.052) involved in the DNA repair process is Ercc6, responsible for the peripheral repair of transcribed sequences. It participates in the regulation of transcription and DNA repair, and ensures the stability of chromosomes. Troelstra et al. have demonstrated in a patient with Cockayne’s syndrome (CS) that ERCC6 is not crucial for cell viability but is specific for preferential repair of transcribed sequences [25].
In clinical practice (especially cancer treatment), drugs, such as cisplatin, bortezomib, or methotrexate, that affect apoptosis, DNA repair, or the cell cycle are used. Despite the effectiveness of these drugs in therapy, they also affect healthy tissues, causing side effects. Cisplatin is an anticancer drug that causes intra- and inter-strand crosslinks in the DNA. Cisplatin has a neurotoxic effect on the peripheral nervous system which is directly related to DNA damage. Dzagnidze et al. observed in mice with dysfunctional deficient for NER that chronic exposure to cisplatin caused accelerated accumulation of unrepaired intrastrand cross-links in neuronal cells, and was significantly correlated to the degree of sensory impairment as measured by electroneurography [26].
We also observed low expression in Kif5A gene (fold change = −3.32) in our study. KIF5A plays a significant role in the transport of organelles in the nerve cells and transfer of mitochondria from the neural soma to the axon terminal [27]. Nam et al. demonstrated that mutations reducing KIF5A expression are important in the development of hereditary spastic paraplegia (HSP) and axonal Charcot–Marie–Tooth peripheral neuropathy type 2 (CMT2). Reduced expression of this gene leads to axon degeneration in the peripheral nervous system [28]. The lipoprotein lipase (Lpl) gene had the most reduced expression in our study (fold change = −42.49). LPL supplies the lipids for myelin sheath by hydrolysing the triglyceride rich lipoproteins. Recent research by Rachana et al. has shown that LPL is an important factor in the regeneration of myelin sheath during development of diabetic neuropathy [29].
Our next interesting observation is the inhibition of mitogen-activated protein kinase/extracellular signal-regulated kinase (MAPK/ERK) signaling (Figure 4). Pathways such as extracellular signal regulated kinase1/2 (ERK1/2), phoshpatidyloinositol-3-kinase (PI3K), or extracellular signal regulated kinase5 (ERK5) are responsible for neuronal survival signaling [30] and mediating the transmission of signals from the synapse to cytoplasmic and nuclear effectors [31]. Davis et al. presented on a rat model that MAPK/ERK signaling is crucial for neuronal plasticity and long-term memory formation [32].
MicroRNAs (miRNA) are small, conserved non-coding RNA molecules involved in the regulation of gene expression. MiRNAs are also markers of disease diagnosis and prognosis [33]. In our study, we observed a small number of miRNA–target gene correlations. Thus, we hypothesize that the observed significant changes in gene expression after bortezomib exposure could probably be attributed to other mechanisms. Nonetheless, we have seen a very significant connection. Namely, expression of rno-miR-130a-3p and rno-miR-152-3p was downregulated and as a result the expression of the Gadd45 increased (fold change = 5.91). This gene is a member of a group of genes, the transcript expression of which is enhanced after stressful growth arrest conditions and treatment with DNA-damaging agents like drugs or mutagens [34].

Potential Study Limitations

Our study showed significant results; however, its weakness is the use of one dose of bortezomib at one time point. In clinical practice, repeated administration of bortezomib is used; therefore, in order to extensively assess the role of bortezomib in the development of drug-induced peripheral neuropathy in the treated patient, it would be advantageous to also use multiple doses in in vitro studies. However, our study was not designed nor extensive enough to unequivocally show molecular basis of BiPN. Our study rather presents miRNA and gene expression changes after single bortezomib exposure in PC12-derived nerve cells. While the nature of this study was rather preliminary, we believe it provides a solid basis for future studies that could fully explain the effect of bortezomib on the development of neuropathy.

4. Materials and Methods

4.1. Cell Culture and Incubation

PC12 pheochromocytoma cells (rat, ATCC, Manassas, VA, USA) were used in this study. The selected cell line is often used as a reliable model of induced peripheral neuropathy and neurotoxicity [35,36,37,38,39]. The cells respond reversibly to nerve growth factor (NGF) by induction of the neuronal phenotype when plated on collagen IV coated culture flasks. PC12 cells were incubated for growth in Corning® T-75 flasks (Corning, MO, USA) in RPMI-1640 (ATCC, Manassas, VO, USA) medium modified to contain 2 mM L-glutamine, 10 mM HEPES, 1 mM sodium pyruvate, 4500 mg/L glucose, and 1500 mg/L sodium bicarbonate, streptomycin (100 µg/mL), penicillin (100 U/mL) with 5% fetal bovine serum (FBS) (Sigma Aldrich, St. Louis, MO, USA) and 10% heat-inactivated horse serum (HS) (Sigma Aldrich, St. Louis, MO, USA) at 37 °C in saturated humidity atmosphere containing 5% CO2. The medium was changed every 2–3 days. PC12 cell differentiation into neurons was conducted according to the protocol described by Popova et al. [39]. PC12 cells were incubated at a density of 2x106 cells/well in 6-well culture plates on collagen IV coated (Corning, MO, USA) in RPMI-1640 medium supplemented with 2 mM L-glutamine, 10 mM HEPES, 1 mM sodium pyruvate, 4500 mg/L glucose, and 1500 mg/L sodium bicarbonate, streptomycin (100 µg/mL), penicillin (100 U/mL), 1% HS, and 100 ng/mL NGF (Sigma Aldrich, St. Louis, MO, USA). Every second day, half of the medium/well was changed. After 10 days, differentiated PC12 cells into neurons were received (Figure 13). Next, neurons were not treated (control) or treated with bortezomib 50 nM/L (Cell Signalling Technology, Danvers, MA, USA) and collected after 24 h of incubation. The cells were used to isolate mRNA and miRNA.

4.2. RNA

RNA was isolated from three separate cell incubations for each group. Total RNA was isolated from PC12-derived nerve cells (1.5 × 106) using the mirVana™ miRNA Isolation Kit (Thermo Fisher, Waltham, MA, USA) following the manufacturer’s instructions. The used kit allows RNA and miRNA isolation simultaneously. Quality and concentration of the received RNA were evaluated using Epoch spectrophotometer (Biotek, Winooski, VT, USA).

4.3. Affymetrix GeneChip Microarray and Data Analysis

RNA was isolated from three separate cell incubations. Next, the isolated total RNA samples were pooled to perform microarray analysis. mRNA was used to produce a complementary single-stranded DNA (cDNA) using an Ambion WT Expression Kit (Thermo Fisher Scientific, Waltham, MA, USA) and labelling with the GeneChip WT Terminal Labeling Kit (Affymetrix, Santa Clara, CA). After this step of the procedure, hybridization was performed on an Affymetrix Rat Gene 2.1 ST Array Strip. Hybridization, washing, and scanning were performed with an Affymetrix GeneAtlas™ System (Affymetrix, Santa Clara, CA, USA). The generated raw data were imported into the BioConductor software based on the statistical R programming language. The Robust Multiarray Averaging (RMA) algorithm implemented in the “affy” package of BioConductor was used for background normalization, correction, and summation of raw data. The BioConductor “oligo” package was used to obtain biological annotation.

4.4. Affymetrix GeneChip miRNA Microarray and Data Analysis

Total RNA-enrichment in miRNAs was isolated from three separate cell incubations. Next, the isolated total RNA samples were pooled to perform microarray analysis. The first step of the procedure was a poly (A) tailing reaction and ligation of the biotinylated signal molecule to the target RNA. Next, the sample was hybridized onto an Affymetrix miRNA 4.1 Array Strip (Affymetrix, Santa Clara, CA, USA). Finally, streptavidin-PE was added and array was prepared for scanning with the Affymetrix GeneAtlas system (Affymetrix, Santa Clara, CA, USA). Microarray data analysis was performed using BioConductor [40]. Normalized data were combined with the “pd.mirna.4.1” description file, including: Types, names, and sequences of miRNAs. The linear models for microarray data implemented in the “limma” library was used to determine differential expression [41].
A database miRTarBase was used to download the list of experimentally validated miRNA target genes. Only targets genes for changed expressed miRNA in the study were subtracted from the whole rat miRNA–target dataset. The target gene list was subjected to functional annotation and clusterization using the DAVID (Database for Annotation, Visualization, and Integrated Discovery) [42]. Target gene symbols of differentially expressed miRNA were uploaded to DAVID by the “RDAVIDWebService” BioConductor library [43], where targets were assigned to suitable Gene Ontology (GO) terms.
The miRNA–target gene interaction network was performed using “networkD3” library. Given that miRNAs mainly function as the repressor of the transcription of target genes, only target genes that negatively correlated with miRNAs were presented. If the expression of miRNAs was upregulated, the expression of its target genes should be downregulated and inversely. With such an assumption, interaction network was performed for all changed miRNA expressions, simultaneously considering only target genes with significantly altered expression.

4.5. DAVID

DAVID Bioinformatics Resources (Database for Annotation, Visualization, and Integrated Discovery) at http://david.abcc.ncifcrf.gov is a functional annotation and enrichment analysis. Functional annotation charts generated by DAVID with overrepresented gene annotations are shown as bubble plots from the BACA BioConductor package (https://cran.r-project.org/web/packages/BACA/BACA.pdf) [42,43].

4.6. Validation of Data Obtained from Microarrays

To validate the expression of genes and miRNA the qRT-PCR method was used. Only genes and miRNAs with significantly changed expression were selected for validation (fold >2). Reaction was performed on a Bio-Rad CFX96 Real-Time PCR Detection System (Bio-Rad Inc., Hercules, CA, USA). All primers for gene validation were designed by BLAST PRIMER and for miRNAs by miRPrimer, and purchased from the Laboratory of DNA Sequencing and Oligonucleotide Synthesis, Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Warsaw, Poland. The qRT-PCR program consisted of: 10-min initial denaturation at 95 °C, denaturation at 95 °C for 15 s, annealing at 56 °C for 5 s, and extension at 72 °C for 10 s. The relative gene expression was quantified using the comparative Ct method (2ΔCt, where ΔCt = [Ct of target genes] − [Ct of endogenous control gene]). All products were characterized by high specificity, which was checked by determining melting points (0.1 °C/s transition rate).

4.6.1. Validation of Gene Expression

The validation method using the qRT-PCR reaction consists of two main steps: (1) Reverse transcription of mRNA (0.1 µg) using the First Strand cDNA Synthesis Kit (Thermo Fisher Scientific, Waltham, MA, USA) and (2) main qRT-PCR reaction. The reaction mixture per sample consisted of: 7.5 μL of SYBR Green PCR Master Mix (Bio-Rad, Hercules, CA, USA); 1.5 μL cDNA template; 1.8 μL specific primers (0.9 μL reverse primer and 0.9 μL forward primer), and 4.2 μL Nuclease-Free Water. Gapdh gene was used as an endogenous control gene.

4.6.2. Validation of Selected miRNA Expression

miRNA reverse-transcription was performed using qScript microRNA cDNA Synthesis Kit (Quanta Biosciences, Beverly, MA, USA) following the manufacturer’s instructions. The mixture for the qRT-PCR reaction contained: 1 μL microRNA cDNA; 5 μL PerfeCTa SYBR Green SuperMix; 0.2 μL microRNA specific primer; 0.2 μL PerfeCTa Universal PCR Primer; 4.6 μL Nuclease-Free Water. miR93 was used as an endogenous control gene.

4.7. Gene Set Enrichment Analysis (GSEA)

GSEA is a computational method focused on gene sets, that is, groups of genes that participate in common biological functions, chromosomal location, or regulation. The method uses for calculations the Kolmogorov–Smirnov (K-S) statistical test to identify significantly enriched or depleted groups of genes [44]. GSEA analysis was conducted using the GSEA Java Desktop Application from the Broad Institute (http://software.broadinstitute.org/gsea/index.jsp). Normalized data from all genes were transformed into a suitable format and imported into the application (GSEA Java Desktop Application, http://software.broadinstitute.org/gsea/index.jsp). Then, the Molecular Signatures Database (MsigDB) was selected to create a predefined gene set database [45]. The signal-to-noise ratio with 1000 permutations allowed classified genes into the selected set according to the difference in their expression level. Calculation of the enrichment score (ES) for each selected gene set was performed using a sum statistic [45]. Enrichment scores were normalized by their gene set size, and false-positive findings were corrected by FDR. Significant gene sets were considered to be those with an adjusted nominal p value < 0.01 and FRD q value < 25%.

4.8. Statistical Methods

Arithmetic means and standard deviations were calculated using MS Excel. Parameter comparisons between the two groups were performed using the unpaired Student’s t-test, and a p-value of <0.05 was considered statistically significant.

5. Conclusions

The study shows the effect of single exposure to bortezomib on gene and miRNA expression in PC12-derived nerve cells. We have shown that bortezomib negatively affects the cell cycle, DNA repair, nerve processes and cell division. Disorders of these processes could have a negative influence on peripheral nervous system homeostasis and may be associated with neurotoxicity of bortezomib. Further studies on the development of drug-induced peripheral neuropathy are required to provide in-depth explanation of bortezomib action in the neural system.

Author Contributions

K.Ł. and D.R. performed the experiments and analysis, K.Ł. and Z.U. wrote the manuscript, B.M. analyzed the data and edited the manuscript, supervised the study. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Science Centre grant OPUS 2016/23/N/NZ7/03363

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

ALSamyotrophic lateral sclerosis
A-Tataxia telangiectasia
BiPNbortezomib-induced peripheral neuropathy
BZTbortezomib
CNScentral nervous system
CSCockayne’s syndrome
DAVIDDatabase for Annotation, Visualization, and Integrated Discovery
DSBDNA double-strand break
ESenrichment score
FBSfetal bovine serum
GOGene Ontology
GSEAGene Set Enrichment Analysis
HSheat-inactivated horse serum
HSPhereditary spastic paraplegia
IGF-1insulin-like growth factor 1
IMIDincluding immunomodulatory drug
IL-6interleukin-6
MAPK/ERKmitogen-activated protein kinase/extracellular signal-regulated kinase
miRNAsmicroRNAs
MMmultiple myeloma
MsigDBMolecular Signatures Database
NERnucleotide excision repair
NFκBnuclear factor kappa B
NGFnerve growth factor
NHEJnon-homologous end joining
PNperipheral neuropathy
PNSperipheral nervous system
qRT-PCRquantitative polymerase chain reaction
SCAN-1spinocerebellar ataxia with axonal neuropathy
TNF-αtumor necrosis factor-α
VEGFendothelial growth factor

References

  1. Rajkumar, S.V.; Kumar, S. Multiple Myeloma: Diagnosis and Treatment. Mayo Clin. Proc. 2016, 91, 101–119. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Quach, H.; Ritchie, D.; Stewart, A.K.; Neeson, P.; Harrison, S.; Smyth, M.J.; Prince, H.M. Mechanism of action of immunomodulatory drugs (IMiDS) in multiple myeloma. Leukemia 2010, 24, 22–32. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Adams, J.; Kauffman, M. Development of the proteasome inhibitor Velcade (Bortezomib). Cancer Investig. 2004, 22, 304–311. [Google Scholar] [CrossRef] [PubMed]
  4. Hideshima, T.; Mitsiades, C.; Akiyama, M.; Hayashi, T.; Chauhan, D.; Richardson, P.; Schlossman, R.; Podar, K.; Munshi, N.C.; Mitsiades, N.; et al. Molecular mechanisms mediating antimyeloma activity of proteasome inhibitor PS-341. Blood 2003, 101, 1530–1534. [Google Scholar] [CrossRef] [Green Version]
  5. Rajkumar, S.V.; Richardson, P.G.; Hideshima, T.; Anderson, K.C. Proteasome inhibition as a novel therapeutic target in human cancer. J. Clin. Oncol. 2005, 23, 630–639. [Google Scholar] [CrossRef]
  6. Landowski, T.H.; Megli, C.J.; Dorr, R.T.; Dalton, W.S. Calcium-mediated mitochondrial perturbation is a critical determinant of Velcade (PS-341/bortezomib) cytotoxicity in myeloma cell lines. Proc. Am. Assoc. Cancer Res. 2004, 67, 520. [Google Scholar]
  7. Richardson, P.G.; Delforge, M.; Beksac, M.; Wen, P.; Jongen, J.L.; Sezer, O.; Terpos, E.; Munshi, N.; Palumbo, A.; Rajkumar, S.V.; et al. Management of treatment-emergent peripheral neuropathy in multiple myeloma. Leukemia 2012, 26, 595–608. [Google Scholar] [CrossRef] [Green Version]
  8. Argyriou, A.A.; Iconomou, G.; Kalofonos, H.P. Bortezomib-induced peripheral neuropathy in multiple myeloma: A comprehensive review of the literature. Blood 2008, 112, 1593–1599. [Google Scholar] [CrossRef] [Green Version]
  9. Chen, D.; Frezza, M.; Schmitt, S.; Kanwar, J.; Dou, Q.P. Bortezomib as the first proteasome inhibitor anticancer drug: Current status and future perspectives. Curr. Cancer Drug Targets 2011, 11, 239–253. [Google Scholar] [CrossRef] [Green Version]
  10. Chu, Y.-Y.; Ko, C.-Y.; Wang, S.-M.; Lin, P.-I.; Wang, H.-Y.; Lin, W.-C.; Wu, D.-Y.; Wang, L.-H.; Wang, J.-M. Bortezomib-induced miRNAs direct epigenetic silencing of locus genes and trigger apoptosis in leukemia. Cell Death Dis. 2017, 8, e3167. [Google Scholar] [CrossRef] [Green Version]
  11. Łuczkowska, K.; Rogińska, D.; Ulańczyk, Z.; Paczkowska, E.; Schmidt, A.C.; Machaliński, B. Molecular Mechanisms of Bortezomib Action: Novel Evidence for the miRNA–mRNA Interaction Involvement. Int. J. Mol. Sci. 2020, 21, 350. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Wang, J.; Chen, J.; Sen, S. MicroRNA as Biomarkers and Diagnostics. J. Cell. Physiol. 2016, 231, 25–30. [Google Scholar] [CrossRef] [PubMed]
  13. Schraml, E.; Hackl, M.; Grillari, J. MicroRNAs and toxicology: A love marriage. Toxicol. Rep. 2017, 4, 634–636. [Google Scholar] [CrossRef] [PubMed]
  14. Łuczkowska, K.; Litwińska, Z.; Paczkowska, E.; Machaliński, B. Pathophysiology of drug-induced peripheral neuropathy in multiple myeloma patients. J. Physiol. Pharmacol. 2018, 69, 165–172. [Google Scholar] [CrossRef]
  15. Field-Smith, A.; Morgan, G.J.; Davies, F.E. Bortezomib (Velcade™) in the Treatment of Multiple Myeloma. Ther. Clin. Risk Manag. 2006, 2, 271–279. [Google Scholar] [CrossRef] [Green Version]
  16. Landowski, T.H.; Megli, C.J.; Nullmeyer, K.D.; Lynch, R.M.; Dorr, R.T. Mitochondrial-mediated disregulation of Ca2+ is a critical determinant of Velcade (PS-341/bortezomib) cytotoxicity in myeloma cell lines. Cancer Res. 2005, 65, 3828–3836. [Google Scholar] [CrossRef] [Green Version]
  17. Karademir, B.; Sari, G.; Jannuzzi, A.T.; Musunuri, S.; Wicher, G.; Grune, T.; Mi, J.; Hacioglu-Bay, H.; Forsberg-Nilsson, K.; Bergquist, J.; et al. Proteomic approach for understanding milder neurotoxicity of Carfilzomib against Bortezomib. Sci. Rep. 2018, 8, 16318. [Google Scholar] [CrossRef]
  18. Su, Y.; Ming, G.L.; Song, H. DNA damage and repair regulate neuronal gene expression. Cell Res. 2015, 25, 993–994. [Google Scholar] [CrossRef] [Green Version]
  19. Englander, E.W. DNA Damage Response in Peripheral Nervous System: Coping with Cancer Therapy-Induced DNA Lesions. DNA Repair (Amst.) 2013, 12, 685–690. [Google Scholar] [CrossRef] [Green Version]
  20. McKinnon, P.J. Topoisomerases and the regulation of neural function. Nat. Rev. Neurosci. 2016, 17, 673–679. [Google Scholar] [CrossRef] [Green Version]
  21. Takashima, H.; Boerkoel, C.F.; John, J.; Saifi, G.M.; Salih, M.A.; Armstrong, D.; Mao, Y.; Quiocho, F.A.; Roa, B.B.; Nakagawa, M.; et al. Mutation of TDP1, encoding a topoisomerase I-dependent DNA damage repair enzyme, in spinocerebellar ataxia with axonal neuropathy. Nat. Genet. 2002, 32, 267–272. [Google Scholar] [CrossRef] [PubMed]
  22. Lee, Y.; Katyal, S.; Downing, S.M.; Zhao, J.; Russell, H.R.; McKinnon, P.J. Neurogenesis requires TopBP1 to prevent catastrophic replicative DNA damage in early progenitors. Nat. Neurosci. 2012, 15, 819–826. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Lu, G.; Duan, J.; Shu, S.; Wang, X.; Gao, L.; Guo, J.; Zhang, Y. Ligase I and ligase III mediate the DNA double-strand break ligation in alternative end-joining. Proc. Natl. Acad. Sci. USA 2016, 113, 1256–1260. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ellenberger, T.; Tomkinson, A.E. Eukaryotic DNA Ligases: Structural and Functional Insights. Annu. Rev. Biochem. 2008, 77, 313–338. [Google Scholar] [CrossRef] [Green Version]
  25. Troelstra, C.; van Gool, A.; de Wit, J.; Vermeulen, W.; Bootsma, D.; Hoeijmakers, J.H. ERCC6, a member of a subfamily of putative helicases, is involved in Cockayne’s syndrome and preferential repair of active genes. Cell 1992, 71, 939–953. [Google Scholar] [CrossRef] [Green Version]
  26. Dzagnidze, A.; Katsarava, Z.; Makhalova, J.; Liedert, B.; Yoon, M.S.; Kaube, H.; Limmroth, V.; Thomale, J. Repair capacity for platinum-DNA adducts determines the severity of cisplatin-induced peripheral neuropathy. J. Neurosci. 2007, 27, 9451–9457. [Google Scholar] [CrossRef] [Green Version]
  27. Bermudez-Guzman, L.; Leal, A. DNA repair deficiency in neuropathogenesis: When all roads lead to mitochondria. Transl. Neurodegener. 2019, 8, 14. [Google Scholar] [CrossRef]
  28. Nam, D.E.; Yoo, D.H.; Choi, S.S.; Choi, B.O.; Chung, K.W. Wide phenotypic spectrum in axonal Charcot-Marie-Tooth neuropathy type 2 patients with KIF5A mutations. Genes Genom. 2018, 40, 77–84. [Google Scholar] [CrossRef]
  29. Rachana, K.S.; Manu, M.S.; Advirao, G.M. Insulin-induced upregulation of lipoprotein lipase in Schwann cells during diabetic peripheral neuropathy. Diabetes Metab. Syndr. 2018, 12, 525–530. [Google Scholar] [CrossRef]
  30. Hetman, M.; Vashishta, A.; Rempala, G. Neurotoxic mechanisms of DNA damage: Focus on transcriptional inhibition. J. Neurochem. 2010, 114, 1537–1549. [Google Scholar] [CrossRef] [Green Version]
  31. Blüthgen, N.; van Bentum, M.; Merz, B.; Kuhl, D.; Hermey, G. Profiling the MAPK/ERK dependent and independent activity regulated transcriptional programs in the murine hippocampus in vivo. Sci. Rep. 2017, 7, 45101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Davis, S.; Vanhoutte, P.; Pages, C.; Caboche, J.; Laroche, S. The MAPK/ERK cascade targets both Elk-1 and cAMP response element-binding protein to control long-term potentiation-dependent gene expression in the dentate gyrus in vivo. J. Neurosci. 2000, 20, 4563–4572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Macfarlane, L.A.; Murphy, P.R. MicroRNA: Biogenesis, Function and Role in Cancer. Curr. Genom. 2010, 11, 537–561. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Schafer, A. Gadd45 proteins: Key players of repair-mediated DNA demethylation. Adv. Exp. Med. Biol. 2013, 793, 35–50. [Google Scholar] [PubMed]
  35. Keilbaugh, S.A.; Prusoff, W.H.; Simpson, M.V. The PC12 cell as a model for studies of the mechanism of induction of peripheral neuropathy by anti-HIV-1 dideoxynucleoside analogs. Biochem. Pharmacol. 1991, 42, R5–R8. [Google Scholar] [CrossRef]
  36. Zhao, S.; Zhang, L.; Xu, Z.; Chen, W. Neurotoxic effects of iron overload under high glucose concentration. Neural Regen. Res. 2013, 8, 3423–3433. [Google Scholar]
  37. Slotkin, T.A.; MacKillop, E.A.; Ryde, I.T.; Tate, C.A.; Seidler, F.J. Screening for Developmental Neurotoxicity Using PC12 Cells: Comparisons of Organophosphates with a Carbamate, an Organochlorine, and Divalent Nickel. Environ. Health Perspect. 2007, 115, 93–101. [Google Scholar] [CrossRef]
  38. Christen, V.; Rusconi, M.; Crettaz, P.; Fent, K. Developmental neurotoxicity of different pesticides in PC-12 cells in vitro. Toxicol. Appl. Pharmacol. 2017, 325, 25–36. [Google Scholar] [CrossRef]
  39. Popova, D.; Karlsson, J.; Jacobsson, S.O.P. Comparison of neurons derived from mouse P19, rat PC12 and human SH-SY5Y cells in the assessment of chemical-and toxin-induced neurotoxicity. BMC Pharmacol. Toxicol. 2017, 18, 42. [Google Scholar] [CrossRef]
  40. Gautier, L.; Cope, L.; Bolstad, B.M.; Irizarry, R.A. Affy—Analysis of Affymetrix GeneChip data at the probe level. Bioinformatics 2004, 20, 307–315. [Google Scholar] [CrossRef]
  41. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  42. Dennis, G., Jr.; Sherman, B.T.; Hosack, D.A.; Yang, J.; Gao, W.; Lane, H.C.; Lempicki, R.A. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4, R60. [Google Scholar] [CrossRef] [Green Version]
  43. Fresno, C.; Fernandez, E.A. RDAVIDWebService: A versatile R interface to DAVID. Bioinformatics 2013, 29, 2810–2811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Liberzon, A.; Birger, C.; Thorvaldsdottir, H.; Ghandi, M.; Mesirov, J.P.; Tamayo, P. The Molecular Signatures Database (MSigDB) hallmark gene set collection. Cell Syst. 2015, 1, 417–425. [Google Scholar] [CrossRef] [Green Version]
  45. Subramanian, A.; Tamayo, P.; Mootha, V.K.; Mukherjee, S.; Ebert, B.L.; Gillette, M.A.; Paulovich, A.; Pomeroy, S.L.; Golub, T.R.; Lander, E.S.; et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl. Acad. Sci. USA 2005, 102, 15545–15550. [Google Scholar] [CrossRef] [Green Version]
Ijms 21 00751 g001 550
Figure 1. The scatter plot of global gene expression in bortezomib-treated PC12-derived nerve cells compared to the control cells. Genes are represented by dots (red colour: Downregulation; green colour: Upregulation). The graph shows the genes with at least twofold change and p < 0.05. The graph also contains the symbol of the genes with the largest change in expression.
Figure 1. The scatter plot of global gene expression in bortezomib-treated PC12-derived nerve cells compared to the control cells. Genes are represented by dots (red colour: Downregulation; green colour: Upregulation). The graph shows the genes with at least twofold change and p < 0.05. The graph also contains the symbol of the genes with the largest change in expression.
Ijms 21 00751 g001
Ijms 21 00751 g002 550
Figure 2. The bubble plot with changed biological processes assigned according to Gene Ontology (GO) classification in bortezomib-treated PC12-derived nerve cells compared to the control cells. Genes assigned to individual processes fulfilling the criteria: Adjusted p < 0.05, method = Benjamini, and minimum number of genes per group = 5, are presented. The bubble size indicates the number of genes represented in the corresponding annotation.
Figure 2. The bubble plot with changed biological processes assigned according to Gene Ontology (GO) classification in bortezomib-treated PC12-derived nerve cells compared to the control cells. Genes assigned to individual processes fulfilling the criteria: Adjusted p < 0.05, method = Benjamini, and minimum number of genes per group = 5, are presented. The bubble size indicates the number of genes represented in the corresponding annotation.
Ijms 21 00751 g002
Ijms 21 00751 g003 550
Figure 3. The bubble plot with changed pathways assigned according to the KEGG Pathway Database in bortezomib-treated PC12-derived nerve cells compared to the control cells. Genes assigned to individual processes fulfilling the criteria: Adjusted p < 0.05, method = Benjamini, and minimum number of genes per group = 5, are presented. The bubble size indicates the number of genes represented in the corresponding annotation.
Figure 3. The bubble plot with changed pathways assigned according to the KEGG Pathway Database in bortezomib-treated PC12-derived nerve cells compared to the control cells. Genes assigned to individual processes fulfilling the criteria: Adjusted p < 0.05, method = Benjamini, and minimum number of genes per group = 5, are presented. The bubble size indicates the number of genes represented in the corresponding annotation.
Ijms 21 00751 g003
Ijms 21 00751 g004 550
Figure 4. Graph of genes and processes involved in the pathway 04360.Axon guidance. The scheme shows genes and relationships between them. Additionally, the diagram indicates the genes up- or downregulated in our analysis. Solid arrows indicate a direct effect on specific gene expression while dotted arrows indicate gene involvement in a process or pathway. This allows an in-depth understanding of the effects of specific gene dysregulation on the selected process.
Figure 4. Graph of genes and processes involved in the pathway 04360.Axon guidance. The scheme shows genes and relationships between them. Additionally, the diagram indicates the genes up- or downregulated in our analysis. Solid arrows indicate a direct effect on specific gene expression while dotted arrows indicate gene involvement in a process or pathway. This allows an in-depth understanding of the effects of specific gene dysregulation on the selected process.
Ijms 21 00751 g004
Ijms 21 00751 g005 550
Figure 5. Graph of genes and processes involved in the pathway 04110: Cell cycle. The scheme shows genes and relationships between them. Additionally, the diagram indicates the genes up- or downregulated in our analysis. Solid arrows indicate a direct effect on specific gene expression while dotted arrows indicate gene involvement in a process or pathway. This allows an in-depth understanding of the effects of specific gene dysregulation on the selected process.
Figure 5. Graph of genes and processes involved in the pathway 04110: Cell cycle. The scheme shows genes and relationships between them. Additionally, the diagram indicates the genes up- or downregulated in our analysis. Solid arrows indicate a direct effect on specific gene expression while dotted arrows indicate gene involvement in a process or pathway. This allows an in-depth understanding of the effects of specific gene dysregulation on the selected process.
Ijms 21 00751 g005
Ijms 21 00751 g006 550
Figure 6. Circus plot shows the dysregulated processes (neuron projection morphogenesis and neurotransmitter secretion) and their associated genes after bortezomib-treatment in PC12-derived nerve cells. The level of expression for each gene is marked using red color (fold change −3 to −2).
Figure 6. Circus plot shows the dysregulated processes (neuron projection morphogenesis and neurotransmitter secretion) and their associated genes after bortezomib-treatment in PC12-derived nerve cells. The level of expression for each gene is marked using red color (fold change −3 to −2).
Ijms 21 00751 g006
Ijms 21 00751 g007 550
Figure 7. Circus plot shows the dysregulated processes (DNA repair, regulation of cell cycle, cell division) and their associated genes after bortezomib-treatment in PC12-derived nerve cells. The level of expression for each gene is marked using red color (fold change −3 to −2). Ribbons connecting areas of the Circus plots also indicate shared genes between groups.
Figure 7. Circus plot shows the dysregulated processes (DNA repair, regulation of cell cycle, cell division) and their associated genes after bortezomib-treatment in PC12-derived nerve cells. The level of expression for each gene is marked using red color (fold change −3 to −2). Ribbons connecting areas of the Circus plots also indicate shared genes between groups.
Ijms 21 00751 g007
Ijms 21 00751 g008 550
Figure 8. Gene sets enriched in bortezomib-treatment PC12-derived nerve cells. Gene sets are ranked according to the normalized enrichment score (NES).
Figure 8. Gene sets enriched in bortezomib-treatment PC12-derived nerve cells. Gene sets are ranked according to the normalized enrichment score (NES).
Ijms 21 00751 g008
Ijms 21 00751 g009 550
Figure 9. Real-time quantitation of selected genes Ascl1, Top2a, Slc7a11, Ddit3 in PC12-derived nerve cells after bortezomib-treatment, and in the control cells. Data are presented as the mean ± SD (n = 3). ** p < 0.002; *** p < 0.001; **** p < 0.0009 compared with the control cells.
Figure 9. Real-time quantitation of selected genes Ascl1, Top2a, Slc7a11, Ddit3 in PC12-derived nerve cells after bortezomib-treatment, and in the control cells. Data are presented as the mean ± SD (n = 3). ** p < 0.002; *** p < 0.001; **** p < 0.0009 compared with the control cells.
Ijms 21 00751 g009
Ijms 21 00751 g010 550
Figure 10. The scatter plot of global miRNAs expression in bortezomib-treated PC12-derived nerve cells compared to the control cells. Downregulated miRNAs are represented by red dots. The graph shows the miRNA with at least twofold change and p<0.05. The graph also contains the symbol of the miRNA with the largest change in expression.
Figure 10. The scatter plot of global miRNAs expression in bortezomib-treated PC12-derived nerve cells compared to the control cells. Downregulated miRNAs are represented by red dots. The graph shows the miRNA with at least twofold change and p<0.05. The graph also contains the symbol of the miRNA with the largest change in expression.
Ijms 21 00751 g010
Ijms 21 00751 g011 550
Figure 11. The diagram shows significantly downregulated miRNAs (at least fold −2). Target genes are assigned to each miRNA with a marked change in expression (at least fold −2) (red indicates upregulation and blue indicates downregulation).
Figure 11. The diagram shows significantly downregulated miRNAs (at least fold −2). Target genes are assigned to each miRNA with a marked change in expression (at least fold −2) (red indicates upregulation and blue indicates downregulation).
Ijms 21 00751 g011
Ijms 21 00751 g012 550
Figure 12. Real-time quantitation of selected miRNAs (miRNA-21-5p; miRNA322-5p; miRNA-532-5p) in PC12-derived nerve cells after bortezomib-treatment, and in control cells. Data are presented as the mean ± SD (n = 3). *p < 0.05; ** p < 0.01; *** p < 0.002 compared with the control cells.
Figure 12. Real-time quantitation of selected miRNAs (miRNA-21-5p; miRNA322-5p; miRNA-532-5p) in PC12-derived nerve cells after bortezomib-treatment, and in control cells. Data are presented as the mean ± SD (n = 3). *p < 0.05; ** p < 0.01; *** p < 0.002 compared with the control cells.
Ijms 21 00751 g012
Ijms 21 00751 g013 550
Figure 13. PC12 cells after 10 days of incubation with nerve growth factor (NGF).
Figure 13. PC12 cells after 10 days of incubation with nerve growth factor (NGF).
Ijms 21 00751 g013
Table
Table 1. The list of 15 the most downregulated genes in bortezomib-treated PC12-derived nerve cells compared to controls.
Table 1. The list of 15 the most downregulated genes in bortezomib-treated PC12-derived nerve cells compared to controls.
Gene SymbolGene NameGene FunctionFold Change
Lpllipoprotein lipaseRecruited to its site of action on the luminal surface of vascular endothelium by binding to GPIHBP1 and cell surface heparan sulfate proteoglycans−42.493
Ddcdopa decarboxylase (aromatic L-amino acid decarboxylase)important in the brain and nervous system. This enzyme takes part in the pathway that produces dopamine and serotonin, which are chemical messengers that transmit signals between nerve cells (neurotransmitters).−38.562
Lgr5leucine rich repeat containing G protein coupled receptor 5expressed across a diverse range of tissue such as in the muscle, placenta, spinal cord and brain and particularly as a biomarker of adult stem cells in certain tissues−32.995
Gpr22G protein-coupled receptor 22play an essential role in the regulation of cardiovascular function−30.754
Top2atopoisomerase (DNA) II alphacontrols and alters the topologic states of DNA during transcription; involved in processes such as chromosome condensation, chromatid separation, and the relief of torsional stress that occurs during DNA transcription and replication−29.831
Nusap1nucleolar and spindle associated protein 1binds and stabilizes microtubules (by similarity); can promote the organization of mitotic spindle microtubules around them−29.200
Ascl1achaete-scute family bHLH transcription factor 1 neuronal commitment and differentiation and in the generation of olfactory and autonomic neurons−28.366
HjurpHolliday junction recognition proteinplays a central role in the incorporation and maintenance of histone H3-like variant CENPA at centromeres−27.183
Mx2myxovirus (influenza virus) resistance 2exhibits antiviral activity; play a role in regulating nucleocytoplasmic transport and cell-cycle progression−24.350
RGD1564463similar to Mdes proteinable to upregulation of proliferation-related gene expression, such as that of CDX2 and PCNA−22.386
S100a5S100 calcium binding protein A5involved in the regulation of a number of cellular processes such as cell cycle progression and differentiation.−22.222
Bub1BUB1 mitotic checkpoint serine/threonine kinaseestablishment of the mitotic spindle checkpoint and chromosome congression−21.434
Ndufa4l2NADH dehydrogenase (ubiquinone) 1 alpha subcomplex, 4-like 2codes for a subunit of Complex I of the respiratory chain, which transfers electrons from NADH to ubiquinone−20.830
Nuf2NUF2, NDC80 kinetochore complex componentchromosome segregation−19.673
Faxdc2fatty acid hydroxylase domain containing 2catalyzes the synthesis of 2-hydroxysphingolipids−17.710
Table
Table 2. The list of 15 the most upregulated genes in bortezomib-treated PC12-derived nerve cells compared to controls.
Table 2. The list of 15 the most upregulated genes in bortezomib-treated PC12-derived nerve cells compared to controls.
Gene SymbolGene NameGene FunctionFold Change
Hamphepcidin antimicrobial peptidemaintaining iron balance in the body78.308
Ankrd1ankyrin repeat domain 1transcription factor involved in development and under conditions of stress36.103
Flrt3fibronectin leucine rich transmembrane protein 3cell adhesion and receptor signalling33.410
Fam171bfamily with sequence similarity 171, member Bassociated with neurodevelopmental disorders22.564
Slc7a11solute carrier family 7regulates synaptic activity by stimulating extrasynaptic receptors and performs nonvesicular glutamate release17.940
Fam129afamily with sequence similarity 129, member Amarker for certain cancers; regulates p53-mediated apoptosis16.321
Gata4GATA binding protein 4transcriptional regulator for many cardiac genes; regulates hypertrophic growth of the heart; promotes cardiac morphogenesis, cardiomyocytes survival, and maintains cardiac function in the adult heart.15.599
Sema6asema domain, transmembrane domain ™promotes reorganization of the actin cytoskeleton;plays an important role in axon guidance in the developing central nervous system15.091
Myo10myosin Xregulates cell shape, cell spreading and cell adhesion15.012
Lifrleukemia inhibitory factor receptor alphacontrols several cellular processes, including growth and division (proliferation), maturation (differentiation), and survival14.617
Prss46protease, serine, 46serine-type endopeptidase activity; involved in proteolysis13.349
Plaurplasminogen activator, urokinase receptorinfluences many normal and pathological processes related to cell-surface plasminogen activation and localized degradation of the extracellular matrix13.144
Dcxdoublecortininvolved in the movement of nerve cells to their proper locations in the developing brain, a process called neuronal migration13.099
Fosl1fos-like antigen 1regulates tumor cell proliferation and survival12.403
Ddit3DNA-damage inducible transcript 3adipogenesis, erythropoiesis, growth arrest and endoplasmic reticulum stress response7.306

Share and Cite

MDPI and ACS Style

Łuczkowska, K.; Rogińska, D.; Ulańczyk, Z.; Machaliński, B. Effect of Bortezomib on Global Gene Expression in PC12-Derived Nerve Cells. Int. J. Mol. Sci. 2020, 21, 751. https://doi.org/10.3390/ijms21030751

AMA Style

Łuczkowska K, Rogińska D, Ulańczyk Z, Machaliński B. Effect of Bortezomib on Global Gene Expression in PC12-Derived Nerve Cells. International Journal of Molecular Sciences. 2020; 21(3):751. https://doi.org/10.3390/ijms21030751

Chicago/Turabian Style

Łuczkowska, Karolina, Dorota Rogińska, Zofia Ulańczyk, and Bogusław Machaliński. 2020. "Effect of Bortezomib on Global Gene Expression in PC12-Derived Nerve Cells" International Journal of Molecular Sciences 21, no. 3: 751. https://doi.org/10.3390/ijms21030751

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