Figures
Abstract
Botulinum neurotoxin type A1 (BoNT/A1) is a potent protein toxin responsible for the potentially fatal human illness botulism. Notwithstanding, the long-lasting flaccid muscle paralysis caused by BoNT/A has led to its utility as a powerful and versatile bio-pharmaceutical. The flaccid paralysis is due to specific cleavage of neuronal SNAREs by BoNTs. However, actions of BoNTs on intoxicated neurons besides the cleavage of SNAREs have not been studied in detail. In this study we investigated by microarray analysis the effects of BoNT/A and a catalytically inactive derivative (BoNT/A ad) on the transcriptome of human induced pluripotent stem cell (hiPSC)-derived neurons at 2 days and 2 weeks after exposure. While there were only minor changes in expression levels at 2 days post exposure, at 2 weeks post exposure 492 genes were differentially expressed more than 2-fold in BoNT/A1-exposed cells when compared to non-exposed populations, and 682 genes were differentially expressed in BoNT/A ad-exposed cells. The vast majority of genes were similarly regulated in BoNT/A1 and BoNT/A ad-exposed neurons, and the few genes differentially regulated between BoNT/A1 and BoNT/A ad-exposed neurons were differentially expressed less than 3.5 fold. These data indicate a similar response of neurons to BoNT/A1 and BoNT/A ad exposure. The most highly regulated genes in cells exposed to either BoNT/A1 or BoNT/A ad are involved in neurite outgrowth and calcium channel sensitization.
Citation: Scherf JM, Hu XS, Tepp WH, Ichtchenko K, Johnson EA, Pellett S (2014) Analysis of Gene Expression in Induced Pluripotent Stem Cell-Derived Human Neurons Exposed to Botulinum Neurotoxin A Subtype 1 and a Type A Atoxic Derivative. PLoS ONE 9(10): e111238. https://doi.org/10.1371/journal.pone.0111238
Editor: Michel R. Popoff, Institute Pasteur, France
Received: July 8, 2014; Accepted: September 19, 2014; Published: October 22, 2014
Copyright: © 2014 Scherf et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. Data have been deposited into NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE58149.
Funding: This work was supported by the National Institute of Allergy and Infectious Diseases Award Number R01 AI093504 and the University of Wisconsin-Madison Food Research Institute (FRI). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: Dr Eric A. Johnson is serving as an editor for PLOS ONE. This does not alter the authors' adherence to PLOS ONE Editorial policies and criteria.
Introduction
Botulinum neurotoxins (BoNTs), the causal agents of botulism, are considered the most potent toxins to humankind [1], [2]. Of the seven currently known serotypes, BoNT/A has been shown to be particularly potent in humans; an ingested dose of 1–2 µg/kg body weight and an intravenous dose of 1–2 ng/kg is estimated to be fatal unless prompt medical treatment is received [3]–. BoNTs prevent neurotransmitter release by targeting the neuromuscular junction (NMJ), leading to flaccid muscular paralysis, with death usually caused by respiratory paralysis [2]. The severe effects of BoNTs are also due to the particularly long persistence of symptoms, lasting several months and up to 1 year in the case of BoNT/A, after which time the patients usually recover completely [2]. Despite the potentially fatal nature of BoNTs, the potency, duration, and reversible action of these neuroparalytic toxins has led to extensive applications in medicine. BoNT/A1, and to a lesser extend BoNT/B1, is now widely used as a unique and important bio-pharmaceutical to treat a variety of neuromuscular disorders, for cosmetic purposes, and to treat a variety of nervous-system-related disorders, such as migraine headaches, muscle spasticity, and various types of dystonia [6]–[8].
BoNTs specifically enter neuronal cells via receptor-mediated cell entry [9]. After binding to polysialogangliosides, BoNTs associate with synaptic vesicle proteins (in the case of BoNT/A SV2C) and are endocytosed. Once inside the endocytic vesicle, acidification of the vesicle lumen leads to a conformational change in the BoNT protein, which ultimately leads to translocation of the light chain (LC) into the cell cytosol [10], [11]. The disulfide bridge connecting the heavy chain (HC) and LC is reduced in the cell's cytosol [10]–[12], and the enzymatically active LC is released. BoNTs interrupt signal transduction by specifically cleaving a soluble N-ethylmaleimide-sensitive-factor attachment protein receptor (SNARE) protein; in the case of BoNT/A, the t-SNARE synaptosomal-associated protein 25 (SNAP-25) [13].
Upon entering the neuronal cell and cleaving SNAP-25, the toxin's LC remains active in the cytosol and cleaves newly synthesized SNAP-25, preventing recovery for an extended period of time [14]–[16]. Recovery is slow, and despite widespread use of BoNT/A as a therapeutic, cellular processes activated inside intoxicated neurons have largely remained unstudied in detail. Due to the specific neuronal cell entry, BoNT has been considered as a potential vehicle for delivering beneficial agents directly to neuronal cells [17]. Owing to the complicated and only partially understood mechanism of cell entry by BoNTs, a recently developed BoNT based specific neuronal delivery vehicle was designed based on a full length atoxic derivative of BoNT [18]. This derivative, BoNT/A ad, has been shown to be catalytically inactive at pM concentrations as used here, and has been proposed as a candidate for a neuron-targeting delivery mechanism [18], [19]. Although the intoxication pathway of BoNT/A ad has been investigated using in situ and neuronal cell models [19], the effect of the catalytically inactive toxin derivative on the human neuronal transcriptome remains unknown, especially compared to the effect of catalytically active BoNT/A. In addition, the effect of physiologically relevant doses of BoNT/A1, the BoNT serotype most commonly used for medical applications today, on the human neuronal transcriptome has not been investigated either.
The few functional studies in neurons exposed to BoNT/A1 have shown that there is post-intoxication neurite outgrowth in various human and animal models [20]–[25]. Other studies have investigated the effect of BoNTs on human epithelial cells and neuronal cell lines [26], as well as human skeletal muscle [27] and keloid fibroblasts [28]. These studies observed differential regulations of a wide range of genes involved in cellular functions such as inflammation, calcium (Ca2+) signaling, and response to oxidative stress after exposure times ranging from 6 hours to 1 year post-intoxication. However, there have been no studies to our knowledge examining the global gene expression response in non-cancerous and non-proliferating human neurons after BoNT/A exposure. In addition, it is not known whether previously observed physiological changes, such as neurite sprouting, are due to SNARE cleavage or to an alternative action of BoNT.
In this study, cDNA microarrays were used to analyze the effects of BoNT/A1 and a catalytically inactive BoNT/A1 derivative (BoNT/A ad) [18] on human neurons derived from human induced pluripotent stem cells (hiPSC). Our data indicate that the hiPSC-derived neurons used in this study and intoxicated with either BoNT/A or BoNT/A ad show few transcriptionally-linked changes after 2 days of intoxication. However, significant transcriptional changes were observed at 14 days post intoxication. Interestingly, gene expression changes were similar in both BoNT/A and BoNT/A ad-intoxicated cells, with strong differential changes in genes related to Ca2+ channel signaling and neurite sprouting, indicating that catalytic activity on SNAP-25 is not necessary for these changes to occur.
Materials and Methods
Botulinum neurotoxin
Isolated pure 150 kDa BoNT/A1 toxin heterodimer was obtained from the Clostridium botulinum strain Hall A-hyper as previously described [29]. Specific activity was determined by mouse bioassay (MBA) [30], [31] to be 1.25×108 mouse LD50 Units/mg. Full-length BoNT/A ad (atoxic derivative), which had an over 100,000-fold reduced toxicity in mice compared to BoNT/A1 (9×102 mouse LD50 Units/mg), was prepared as a recombinant protein expressed in insect cells via a baculovirus expression system as previously described [18].
Neuronal cell culture
The neurons used for this study were cryopreserved hiPSC-derived neurons (iCell Neurons) purchased from Cellular Dynamics International (Madison, WI). The cell population was a 97% pure population consisting of GABAergic and glutamatergic neurons, with a small percentage of dopaminergic neurons (information provided by CDI, Madison, WI). Thawing and plating of the cells was performed as recommended by the manufacturer. Approximately 650,000 live cells were plated per well in a 6-well plate (Techno Plastic Products [TPP]) treated with 0.01% poly-L-ornithine (PLO) (Sigma) and coated with 8.3 µg/ml Matrigel (Biosciences). In order to ensure same cell populations throughout the experiment, the thawed vials of cells were pooled and mixed, before cells were counted and seeded for the entire study. Cells were maintained in media supplied by Cellular Dynamics in a 37°C, 5% CO2 incubator.
Cell exposure to botulinum neurotoxin
Cells were exposed to 10 pM of BoNT/A or BoNT/A ad suspended in 3 ml of neuronal maintenance medium and incubated for 48 h. All cells were gently washed twice with 5 ml of fresh medium to remove extracellular toxin, and RNA was harvested as described below for the 2 d post exposure; the cells were maintained for another 12 d before harvesting RNA for the 14 d post exposure. All experiments were performed in triplicate. Non-exposed control samples were concurrently grown with exposed cells and harvested at the same time, and all mechanical manipulations including addition and removal of toxin to exposed cells was duplicated in the non-exposed samples with fresh medium not containing toxin. In parallel, two more triplicate sets of cells were exposed to the same toxin concentrations to confirm SNAP-25 cleavage. These samples were lysed in 1× lithium dodecyl sulfate (LDS) sample buffer (Invitrogen) and analyzed for SNAP-25 cleavage using Western immunoblotting as previously described [32], [33]. Briefly, analysis was performed via densitometry using a Foto/Analyst FX system and TotalLab Quant software (Fotodyne), correcting for background signal, with data plots being generated with Excel (Microsoft).
Preparation of samples and microarray analysis
Total RNA was isolated from exposed and non-exposed iCell neurons at 2 d and 14 d post-initial exposure using an RNeasy mini kit (Qiagen) with addition of β-mercaptoethanol (βME) to the lysis buffer as recommended by the manufacturer. DNA was removed during isolation using an RNeasy RNase-free DNase kit (Qiagen); for samples analyzed by qPCR an additional DNase treatment was performed (Turbo DNase kit, Ambion). Samples were stored at −20°C after isolation.
RNA labeling and cDNA hybridization was performed by the University of Wisconsin-Madison Biotechnology Center Gene Expression Center (GEC, Madison, WI). After quantification (NanoDrop 2000) and determination of sample integrity (Agilent 2100 Bioanalyzer, RNA 6000 Pico chip), 400 ng of each sample was labeled with an Ambion MessageAmp Premier IVT kit. After cDNA fragmentation, 10 µg of each sample was hybridized to an individual Affymetrix HG U133 Plus 2.0 gene chip according to the methods described in the Affymetrix “Genechip Expression Analysis Technical Manual.” Samples were hybridized for 16 h at 45°C in an Affymetrix Hybridization Oven 640. Post-processing was conducted using an Affymetrix GeneChip Fluidics Station 450 following HG U133 Plus 2.0-specific protocols provided by Affymetrix, and each GeneChip was scanned using an Affymetrix GC3000 G7 scanner. Data was obtained from the scanned images using the Affymetrix Expression Console v 1.2.0.20 software.
Gene enrichment was conducted on all samples using GeneSpring 12.6 (Agilent) software. Robust Multi-array Analysis with correction for oligo GC content (GCRMA) was used for summarization with quantile normalization, and background reduction was performed on all samples by filtering for probes flagged as present by the MAS 5.0 algorithm. Normalized (log 2) data were then subjected to a two-way analysis of variance (ANOVA) with Benjamini Hochberg FDR multiple testing correction and filtered for primary interactions (time, exposure) with P-values <0.05. Relevant interactions were used in determining genes meeting fold expression criteria among the various conditions, with an expression change greater than 2.0 being used to determine significance. Entities with an Affymetrix probe ID but no official gene symbol or designated title were removed from analysis and not counted in the reported data. Analysis of gene enrichment was performed using DAVID Bioinformatics Resources [34], [35], using the GOTERM_BP_FAT, GOTERM_CC_FAT, and GOTERM_MF_FAT datasets. The data discussed in this publication have been deposited in NCBI's Gene Expression Omnibus [36], [37] and are accessible through GEO Series accession number GSE58149.
Validation by real-time quantitative PCR analysis
In order to confirm the expression values from the microarray analysis, qPCR was conducted with all 18 samples used for microarray analysis. Primers were designed using the web-based Primer-BLAST tool (NCBI) [38]. Primers were tested for correct size and amplification of a single product by PCR using non-exposed iCell Neuron RNA as a template, followed by visual confirmation of product using gel electrophoresis (2% agarose gel). Reverse transcription of RNA into cDNA was done using Superscript II (Life Technologies) according to the manufacturer's instructions. Samples were diluted 1∶10 and efficiency tests were carried out on primers using a temperature gradient and DNA melt curve analysis; primers, product lengths, and efficiency values are given in Table 1. H2AFY and RTN4 were chosen as optimal reference genes based upon geNorm algorithm computation (qbase + software, Biogazelle).
QPCR was performed with a Roche LightCycler 480 in triplicate with SYBR Green fluorescent dye (SsoAdvanced SYBR Green, Bio-Rad) in triplicate according to the Roche LightCycler 480 SYBR Green I Master protocol. Fold change was determined using qbase+.
Results
BoNT/A1 but not BoNT/A ad exposure of hiPSC derived neurons resulted in SNAP-25 cleavage
The sensitivity of iCell Neurons to BoNT/A1 has previously been established [39]. In order to confirm catalytic activity of BoNT/A1 and absence of catalytic activity of BoNT/A ad in the cells at the concentrations used for the microarray analysis, iCell Neurons were exposed to 10 pM BoNT/A1 or BoNT/A ad for 48 h. Toxin was removed by washing, and at day 2 (immediately after 48 h exposure) and day 14 cells were lysed and SNAP-25 cleavage was determined by Western immunoblot. As expected, exposure of the neurons to BoNT/A1 resulted in detection of fully cleaved SNAP-25, while only uncleaved SNAP-25 was visible in both the BoNT/A ad exposed cells or the non-exposed control cells (Fig. 1).
A representative Western blot is shown for both 2 days and 14 days exposure. Percentage of SNAP-25 cleavage, conducted via densitometry and averaged for all three replicates after correcting for background values, is represented as a bar graph; error bars represent standard deviation of the three replicates.
There were minimal differences in expression profiles of neurons exposed to BoNT/A1 and BoNT/A ad for 2 days
Comparison of gene expression profiles of BoNT/A1 and BoNT/A ad-exposed cells at 2 d post-initial exposure against their accompanying non-exposed cells, which had undergone all of the same manipulation steps as the exposed cells but without toxin added, revealed only minimal differential gene expression profiles. There were no genes differentially regulated greater than 2-fold between the two exposed cell populations, and only 6 total genes were differentially regulated more than 2-fold in BoNT/A1 or BoNT/A ad-exposed cells compared to non-exposed cells (Table 2).
The general expression profile was significantly altered after 14 days post exposure in both BoNT/A1 and BoNT/A ad-exposed neurons
Comparison of the expression profile of non-exposed iCell Neurons at 2 d and 14 d revealed no differentially expressed genes (P<0.05). This suggests that these hiPSC-derived neurons are transcriptionally stable over the course of 14 d in culture, which is similar to data received from an expression study by the producer of cells (CDI, personal communication).
In contrast, at 14 d post-BoNT/A1 and BoNT/A ad exposures, significant changes were detected in the expression profiles of the exposed neurons as compared to non-exposed neurons. Of 3320 genes correlated with exposure (P<0.05), 354 and 503 genes were differentially regulated 2-fold or greater in BoNT/A1 and BoNT/A ad-exposed cells compared to non-exposed populations. However, the vast majority of these genes showed only minor differences in expression levels (below 5 fold), and only ∼12% were differentially regulated more than 5-fold. Surprisingly, the most highly regulated genes showed increased expression in BoNT/A1 and BoNT/A ad-exposed neurons (Table 3, Table S1). This was unexpected because the BoNT/A ad did not cleave SNAP-25 in the conditions used in this assay.
Most notably, the highly upregulated genes in both BoNT/A1 and BoNT/A ad-exposed neurons encoded proteins involved in neurite outgrowth (FGFR2, RDH10, ADAMTS3, SPARCL1, SALL3, SPON1) and calcium (Ca2+) channel activity (CACNA2D3, CHRM3, NRGN). In addition, several other genes were significantly upregulated, including genes associated with neuronal survival (PDE5A) and certain classes of migraine headaches (ATP1A2). Table 3 lists the most significantly expressed genes and their function. Of the genes that were significantly downregulated many were associated with cellular scaffolding and matrices, important processes involved in cellular proliferation, spatial organization, and adhesion. Examples include COL1A1 [40],COL8A1 [41], and COL11A1 [42], genes encoding subunits of types I, VIII, and XI collagen. Functional annotation clustering analysis of differentially expressed (2-fold or greater) genes in BoNT/A1 and BoNT/A ad-exposed populations at day 14 compared to non-exposed cells using DAVID revealed significant enrichment of genes involved in cell motion, adhesion, and signaling, as well as specific neuron differentiation, development, and projection (Table S2). The differential expression of several genes (CHRM3, DAPK3, MLL, NRGN, RDH10, and RTN4 and H2AFY as reference genes) was confirmed by qPCR (Tables 2 and 3). Overall, the qPCR data indicated an ∼2-fold lower differential gene regulation than indicated by microarray, which suggests that changes detected by microarray that are below 2-fold or even 3-fold may be less significant. There were 36 and 190 genes that are regulated greater than 2-fold exclusively when comparing BoNT/A1 BoNT/A ad to non-exposed cells, respectively. However, many of these genes were still regulated similarly in BoNT/A1 and BoNT/A ad exposed neurons, which are visualized generally (Fig 2) and in greater detail (Table S1). Direct comparison of genes differentially regulated in neurons exposed to BoNT/A1 or BoNT/A ad-exposed at 14 days post exposure showed only 60 genes differentially expressed greater than 2-fold (Table 4, Table S3). The most significant difference was observed for NEGR1, a neuronal growth regulator [43] that was upregulated 3.5-fold by BoNT/A1 compared to BoNT/A ad. Interestingly, when differential expression of NEGR1 in BoNT/A1 or BoNT/A ad exposed neurons was compared to non-exposed neurons, NEGR1 was downregulated significantly in BoNT/A ad-exposed cells, but not in BoNT/A1-exposed cells (-1.6-fold), thus indicating that this gene is not significantly regulated in BoNT/A1-exposed neurons. HSPA6, a heat shock/general stress protein [44], was the only gene that is downregulated (-2.0-fold) in BoNT/A1-exposed cells compared to BoNT/A ad-exposed cells.
All genes differentially expressed greater than 2-fold in either BoNT/A1 or BoNT/A ad-exposed cells are represented, sorted by normalized expression level as determined by the first replicate of BoNT/A1-exposed cells. All three replicates of each condition are shown as individual columns. Blue represents lower expression, while red represents higher expression on a scale of −3.0 to 3.0.
Discussion
BoNT/A1 is widely used as a pharmaceutical in medicine, with more than 3.5 million uses in the US in 2013 alone [45]. One of the hallmarks of these toxins is their long duration of action, lasting for 3–6 months post local injection in humans [2]. In a recent paper we described that this long duration of action was partly due to a persistence of active BoNT/A LC inside the neurons [14], which supported the results of Montecucco's group in seminal studies [12]. However, despite their widespread use and persistence in neurons, little is known of the transcriptional changes that take place in intoxicated neurons. In this report, we describe RNA expression profiles in neurons exposed to BoNT/A1 or to the recombinant BoNT/A1 derivative (BoNT/A ad) which did not show SNAP-25 cleavage in the assay used here. BoNT/A ad is being studied as a delivery platform to target drugs specifically to neuronal cells [18], and results have indicated that it can enter neurons in an analogous mechanism to BoNT/A1 [19].
In this study we showed that there were only minor changes in expression profiles of the neurons exposed to BoNT/A1 or BoNT/A ad compared to non-exposed neurons at 2 days post toxin addition, while there were significant changes in expression profiles after 2 weeks. Interestingly, the more highly regulated genes (FC>5.0) were similarly regulated in BoNT/A1 and BoNT/A ad in comparison to non-exposed cells (Table 3). This suggests abundant transcriptional changes in neurons, independent of SNAP-25 cleavage following BoNT exposure. BoNT/A ad has been shown to be at least 100,000-fold less toxic in mice than BoNT/A1 [19] and does not cleave SNAP-25 in neurons at the concentrations used in this study, and its effects on transcriptional changes are due to unknown mechanisms. Since BoNT/A ad binds to SNAP-25, it is possible that the observed effects are due to an interruption of intracellular signaling by SNAP-25 that is independent of catalytic cleavage, which may be reflected in changes to Ca2+ channel activity [46]. It has previously been observed that neuritogenesis in response to BoNT/A1 intoxication in cultured mouse spinal cord cells is dependent on the binding of the toxin to the neuronal membrane and independent of SNARE cleavage [20]. The majority of highly regulated genes in BoNT/A1 and BoNT/A ad-exposed cells compared to non-exposed cells were involved in neuritogenesis, supporting that the transcriptional changes are independent of SNAP-25 cleavage. The observation of distinct Ca2+ channel-related genes undergoing upregulation in human neurons upon BoNT exposure appears to be novel. However, the hypothesis that changes in SNAP-25 activity affects Ca2+ channel activity, particularly among astrocyte signaling, has been observed in mouse hippocampal neurons [46]. Furthermore, the inhibitory effect of SNAP-25 on Ca2+ channel activity in mouse hippocampal neurons has been reported [47]. The direct effect of BoNT intoxication on Ca2+ channel signaling indicates both an increase [48] and a decrease in Ca2+ utilization [49]. Future studies are required to further investigate these observations on a functional level and to determine whether association or cleavage of SNAP-25 or another property of BoNT/A1 and BoNT/A ad is responsible for the observed transcriptional changes. Finally, it is unlikely that the transcriptional changes were caused by an indirect effect of a contaminant common to both, BoNT/A1 and BoNT/A ad is unlikely, as the two proteins were produced and purified by completely independent procedures. The 150 kDa BoNT/A1 was purified from C. botulinum by published biochemical methods [29], whereas the BoNT/A ad was expressed in and purified from insect cells via a his- and strep-tag [18].
The neuronal cell population used in this study is a comprised primarily of GABAergic and glutamatergic neurons, resembling human forebrain neurons. This system differs from another microarray study that investigated gene expression in both a continuous cell line and epithelial cells [26]. Differences between the results of the two studies may be attributed to the different cell models, as well as differences in exposure parameters to BoNT/A1, including exposure time, toxin concentrations, and the use of BoNT/A1 complex compared to purified neurotoxin. Irrespective of the differences in methodologies, both studies demonstrated upregulation of genes involved in neurite outgrowth. While our cell model provides a human non-cancerous neuronal population, it does not encompass other cell types that are also components of the central nervous system (CNS). Microglial cells, for instance, form a reactive protection network that protects the CNS against harmful stimuli [50]. These cells not only react to environmental stimuli, but also communicate with neuronal cells and exert influence on neuron function [51]. Astrocytes together with neuronal pre- and post-synaptic terminals make up the tripartite synapse in the CNS [52], and these also respond to environmental stimuli and interact with neurons by Ca2+ channel variations [53]. Signals can be transferred along astrocyte chains, delivering messages to and from active and inactive synapses [53]–[55]. The absence of these interacting cells in the neuronal population used in this study may affect the overall transcriptional response of BoNT-exposed neurons. In addition, other neuronal subpopulations, such as motor neurons or nociceptors may have specific responses to BoNTs that could not be revealed in this study.
In summary, these data show that intoxication by BoNT/A1 and BoNT/A ad elicit similar responses in human neurons, particularly in genes that are related to neuritogenesis and Ca2+ channel sensitization. These findings suggest that cleavage of SNAP25 is not necessarily the key process that leads to the transcriptomic changes observed. It is unknown whether cellular entry or physical binding of the toxins to neurons are responsible for the observations described here. Further transcription-based research on the effects of BoNTs as well as their complexing proteins in different cell models is required to determine the mechanisms affecting gene expression. BoNT/A ad has previously been suggested as a general delivery vehicle to target compounds to neurons [18], and the similarities in genes differentially regulated in response to both BoNT/A1 and BoNT/A ad support such use.
Supporting Information
Table S1.
List of all named genes correlated with exposure at 14 d post-exposure with a differential expression greater than 2-fold in either BoNT/A1 or BoNT/A ad-exposed cells when compared to unexposed cells. The data are represented by two worksheets; the first is sorted by differential expression as determined by BoNT/A1 expression compared to control, and the second is sorted alphabetically by gene symbol. Duplicate entires (distinct probes that target identical genes) are grouped together by outline.
https://doi.org/10.1371/journal.pone.0111238.s001
(XLSX)
Table S2.
Functional Annotation Clustering output on genes that are differentially expressed in either 14 d BoNT/A1 or BoNT/A ad-exposed cells compared to non-exposed cells.
https://doi.org/10.1371/journal.pone.0111238.s002
(XLSX)
Table S3.
List of all named genes correlated with exposure at 14 d post-exposure with a differential expression greater than 2-fold in BoNT/A1-exposed cells when compared to BoNT/A ad-exposed cells.
https://doi.org/10.1371/journal.pone.0111238.s003
(XLSX)
Acknowledgments
The authors would like to thank the University of Wisconsin-Madison Biotechnology Center Gene Expression Center for Affymetrix GeneChip labeling and hybridization and qPCR services.
Author Contributions
Conceived and designed the experiments: EAJ SP. Performed the experiments: JMS WHT SP. Analyzed the data: JMS XSH. Contributed reagents/materials/analysis tools: KI. Wrote the paper: JMS XSH WHT EAJ SP. Produced and purified BoNT/A ad: KI.
References
- 1.
Johnson EA (2005) Clostridium botulinum and Clostridium tetani. In: S. P. Borriello PRMaGF, editor.Topley and Wilson's microbiology and microbial infections.London, United Kingdom: Hodder Arnold. pp. 1035–1088.
- 2.
Johnson EA, Montecucco C (2008) Chapter 11 Botulism. In: Andrew GE, editor. Handbook of Clinical Neurology: Elsevier. pp. 333–368.
- 3. Arnon SS, Schechter R, Inglesby TV, Henderson DA, Bartlett JG, et al. (2001) Botulinum toxin as a biological weapon: medical and public health management. JAMA: the journal of the American Medical Association 285: 1059–1070.
- 4. Gimenez DF, Gimenez JA (1995) The typing of botulinal neurotoxins. International journal of food microbiology 27: 1–9.
- 5. Smith TJ, Lou J, Geren IN, Forsyth CM, Tsai R, et al. (2005) Sequence variation within botulinum neurotoxin serotypes impacts antibody binding and neutralization. Infection and immunity 73: 5450–5457.
- 6. Brashear A (2010) Botulinum toxin type A: Exploring new indications. Drugs of today (Barcelona, Spain: 1998) 46: 671–682.
- 7.
Dressler D (2012) Clinical applications of botulinum toxin. Curr Opin Microbiol. England. pp. 325–336.
- 8. Chaddock JA, Marks PM (2006) Clostridial neurotoxins: structure-function led design of new therapeutics. Cellular and molecular life sciences: CMLS 63: 540–551.
- 9. Rummel A (2013) Double receptor anchorage of botulinum neurotoxins accounts for their exquisite neurospecificity. Curr Top Microbiol Immunol 364: 61–90.
- 10. Fischer A, Montal M (2007) Crucial role of the disulfide bridge between botulinum neurotoxin light and heavy chains in protease translocation across membranes. The Journal of biological chemistry 282: 29604–29611.
- 11. Fischer A, Montal M (2007) Single molecule detection of intermediates during botulinum neurotoxin translocation across membranes. Proceedings of the National Academy of Sciences of the United States of America 104: 10447–10452.
- 12. Pantano S, Montecucco C (2014) The blockade of the neurotransmitter release apparatus by botulinum neurotoxins. Cell Mol Life Sci 71: 793–811.
- 13. Blasi J, Chapman ER, Link E, Binz T, Yamasaki S, et al. (1993) Botulinum neurotoxin A selectively cleaves the synaptic protein SNAP-25. Nature 365: 160–163.
- 14. Whitemarsh RC, Tepp WH, Johnson EA, Pellett S (2014) Persistence of botulinum neurotoxin a subtypes 1–5 in primary rat spinal cord cells. PLoS One 9: e90252.
- 15. Foran PG, Mohammed N, Lisk GO, Nagwaney S, Lawrence GW, et al. (2003) Evaluation of the therapeutic usefulness of botulinum neurotoxin B, C1, E, and F compared with the long lasting type A. Basis for distinct durations of inhibition of exocytosis in central neurons. The Journal of biological chemistry 278: 1363–1371.
- 16. Keller JE, Cai F, Neale EA (2004) Uptake of botulinum neurotoxin into cultured neurons. Biochemistry 43: 526–532.
- 17. Bade S, Rummel A, Reisinger C, Karnath T, Ahnert-Hilger G, et al. (2004) Botulinum neurotoxin type D enables cytosolic delivery of enzymatically active cargo proteins to neurones via unfolded translocation intermediates. J Neurochem 91: 1461–1472.
- 18. Band PA, Blais S, Neubert TA, Cardozo TJ, Ichtchenko K (2010) Recombinant derivatives of botulinum neurotoxin A engineered for trafficking studies and neuronal delivery. Protein Expr Purif 71: 62–73.
- 19. Pellett S, Tepp WH, Stanker LH, Band PA, Johnson EA, et al. (2011) Neuronal targeting, internalization, and biological activity of a recombinant atoxic derivative of botulinum neurotoxin A. Biochemical and biophysical research communications 405: 673–677.
- 20. Coffield JA, Yan X (2009) Neuritogenic actions of botulinum neurotoxin A on cultured motor neurons. The Journal of pharmacology and experimental therapeutics 330: 352–358.
- 21. Angaut-Petit D, Molgó J, Comella JX, Faille L, Tabti N (1990) Terminal sprouting in mouse neuromuscular junctions poisoned with botulinum type A toxin: morphological and electrophysiological features. Neuroscience 37: 799–808.
- 22. Bonner PH, Friedli AF, Baker RS (1994) Botulinum A toxin stimulates neurite branching in nerve-muscle cocultures. Brain researchDevelopmental brain research 79: 39–46.
- 23. Juzans P, Comella JX, Molgo J, Faille L, Angaut-Petit D (1996) Nerve terminal sprouting in botulinum type-A treated mouse levator auris longus muscle. Neuromuscul Disord 6: 177–185.
- 24. Diaz J, Molgó J, Pécot-Dechavassine M (1989) Sprouting of frog motor nerve terminals after long-term paralysis by botulinum type A toxin. Neurosci Lett 96: 127–132.
- 25. Holland RL, Brown MC (1981) Nerve growth in botulinum toxin poisoned muscles. Neuroscience 6: 1167–1179.
- 26. Thirunavukkarasusx N, Ghosal KJ, Kukreja R, Zhou Y, Dombkowski A, et al. (2011) Microarray analysis of differentially regulated genes in human neuronal and epithelial cell lines upon exposure to type A botulinum neurotoxin. Biochem Biophys Res Commun 405: 684–690.
- 27.
Mukund K, Mathewson M, Minamoto V, Ward SR, Subramaniam S, et al.. (2014) Systems Analysis of Transcriptional Data Provides Insights Into Muscle's Biological Response to Botulinum Toxin. Muscle Nerve.
- 28. Xiaoxue W, Xi C, Zhibo X (2014) Effects of botulinum toxin type A on expression of genes in keloid fibroblasts. Aesthet Surg J 34: 154–159.
- 29. Malizio CJ, Goodnough MC, Johnson EA (2000) Purification of Clostridium botulinum type A neurotoxin. Methods in molecular biology (Clifton, NJ) 145: 27–39.
- 30. Schantz EJ, Kautter DA (1978) Standardized assay for Clostridium botulinum toxins. Journal - Association of Official Analytical Chemists 61: 96–99.
- 31.
Hatheway CL (1988) Botulism. In: Balows A, Hausler WH, Ohashi M, Turano MA, editors. Laboratory diagnosis of infectious diseases: principles and practice. New York: Springer-Verlag. pp. 111–133.
- 32. Pellett S, Tepp WH, Toth SI, Johnson EA (2010) Comparison of the primary rat spinal cord cell (RSC) assay and the mouse bioassay for botulinum neurotoxin type A potency determination. Journal of pharmacological and toxicological methods 61: 304–310.
- 33. Pellett S, Tepp WH, Clancy CM, Borodic GE, Johnson EA (2007) A neuronal cell-based botulinum neurotoxin assay for highly sensitive and specific detection of neutralizing serum antibodies. FEBS letters 581: 4803–4808.
- 34. Huang dW, Sherman BT, Lempicki RA (2009) Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res 37: 1–13.
- 35. Huang dW, Sherman BT, Lempicki RA (2009) Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 4: 44–57.
- 36. Edgar R, Domrachev M, Lash AE (2002) Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res 30: 207–210.
- 37. Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, et al. (2013) NCBI GEO: archive for functional genomics data sets–update. Nucleic Acids Res 41: D991–995.
- 38. Ye J, Coulouris G, Zaretskaya I, Cutcutache I, Rozen S, et al. (2012) Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics 13: 134.
- 39. Whitemarsh RC, Strathman MJ, Chase LG, Stankewicz C, Tepp WH, et al. (2012) Novel application of human neurons derived from induced pluripotent stem cells for highly sensitive botulinum neurotoxin detection. Toxicological sciences: an official journal of the Society of Toxicology 126: 426–435.
- 40. Rossert J, Terraz C, Dupont S (2000) Regulation of type I collagen genes expression. Nephrol Dial Transplant 15 Suppl 6: 66–68.
- 41. Muragaki Y, Mattei MG, Yamaguchi N, Olsen BR, Ninomiya Y (1991) The complete primary structure of the human alpha 1 (VIII) chain and assignment of its gene (COL8A1) to chromosome 3. Eur J Biochem 197: 615–622.
- 42. Bernard M, Yoshioka H, Rodriguez E, Van der Rest M, Kimura T, et al. (1988) Cloning and sequencing of pro-alpha 1 (XI) collagen cDNA demonstrates that type XI belongs to the fibrillar class of collagens and reveals that the expression of the gene is not restricted to cartilagenous tissue. J Biol Chem 263: 17159–17166.
- 43. Hashimoto T, Maekawa S, Miyata S (2009) IgLON cell adhesion molecules regulate synaptogenesis in hippocampal neurons. Cell Biochem Funct 27: 496–498.
- 44. Noonan EJ, Place RF, Giardina C, Hightower LE (2007) Hsp70B' regulation and function. Cell Stress Chaperones 12: 393–402.
- 45.
American Society for Aesthetic Plastic Surgery 2012 National Data Bank Statistics. Available: www.surgery.org.
- 46. Tafoya LC, Shuttleworth CW, Yanagawa Y, Obata K, Wilson MC (2008) The role of the t-SNARE SNAP-25 in action potential-dependent calcium signaling and expression in GABAergic and glutamatergic neurons. BMC Neurosci 9: 105.
- 47. Pozzi D, Condliffe S, Bozzi Y, Chikhladze M, Grumelli C, et al. (2008) Activity-dependent phosphorylation of Ser187 is required for SNAP-25-negative modulation of neuronal voltage-gated calcium channels. Proc Natl Acad Sci U S A 105: 323–328.
- 48. Fujikawa R, Ishii T, Komori S, Nishimura M (2008) Improved calcium utilization at motor nerve terminals exposed to botulinum neurotoxin in mice. J Physiol Sci 58: 419–424.
- 49. Molgo J, Comella JX, Angaut-Petit D, Pecot-Dechavassine M, Tabti N, et al. (1990) Presynaptic actions of botulinal neurotoxins at vertebrate neuromuscular junctions. J Physiol (Paris) 84: 152–166.
- 50. Hanisch UK, Kettenmann H (2007) Microglia: active sensor and versatile effector cells in the normal and pathologic brain. Nat Neurosci 10: 1387–1394.
- 51. Kettenmann H, Hanisch UK, Noda M, Verkhratsky A (2011) Physiology of microglia. Physiol Rev 91: 461–553.
- 52. Araque A, Parpura V, Sanzgiri RP, Haydon PG (1999) Tripartite synapses: glia, the unacknowledged partner. Trends Neurosci 22: 208–215.
- 53. Carmignoto G (2000) Reciprocal communication systems between astrocytes and neurones. Prog Neurobiol 62: 561–581.
- 54. Nedergaard M (1994) Direct signaling from astrocytes to neurons in cultures of mammalian brain cells. Science 263: 1768–1771.
- 55. Finkbeiner S (1992) Calcium waves in astrocytes-filling in the gaps. Neuron 8: 1101–1108.
- 56. Noureddini M, Verdi J, Mortazavi-Tabatabaei SA, Sharif S, Azimi A, et al. (2012) Human endometrial stem cell neurogenesis in response to NGF and bFGF. Cell Biol Int 36: 961–966.
- 57. Terwisscha van Scheltinga AF, Bakker SC, Kahn RS, Kas MJ (2013) Fibroblast growth factors in neurodevelopment and psychopathology. Neuroscientist 19: 479–494.
- 58. Siegenthaler JA, Ashique AM, Zarbalis K, Patterson KP, Hecht JH, et al. (2009) Retinoic acid from the meninges regulates cortical neuron generation. Cell 139: 597–609.
- 59. Lee SA, Belyaeva OV, Wu L, Kedishvili NY (2011) Retinol dehydrogenase 10 but not retinol/sterol dehydrogenase(s) regulates the expression of retinoic acid-responsive genes in human transgenic skin raft culture. J Biol Chem 286: 13550–13560.
- 60. Bolz H, Ramírez A, von Brederlow B, Kubisch C (2001) Characterization of ADAMTS14, a novel member of the ADAMTS metalloproteinase family. Biochim Biophys Acta 1522: 221–225.
- 61. Kähler AK, Djurovic S, Kulle B, Jönsson EG, Agartz I, et al. (2008) Association analysis of schizophrenia on 18 genes involved in neuronal migration: MDGA1 as a new susceptibility gene. Am J Med Genet B Neuropsychiatr Genet 147B: 1089–1100.
- 62. Harrison SJ, Parrish M, Monaghan AP (2008) Sall3 is required for the terminal maturation of olfactory glomerular interneurons. J Comp Neurol 507: 1780–1794.
- 63. Cheng YC, Chen TA, Chen CY, Liang CM, Liang SM (2012) 3′poly-G-tailed ODNs inhibit F-spondin to induce cell death and neurite retraction in rat embryonic neurons. Mol Neurobiol 45: 536–549.
- 64. Peterziel H, Sackmann T, Strelau J, Kuhn PH, Lichtenthaler SF, et al. (2011) F-spondin regulates neuronal survival through activation of disabled-1 in the chicken ciliary ganglion. Mol Cell Neurosci 46: 483–497.
- 65. Takita J, Chen Y, Okubo J, Sanada M, Adachi M, et al. (2011) Aberrations of NEGR1 on 1p31 and MYEOV on 11q13 in neuroblastoma. Cancer Sci 102: 1645–1650.
- 66. Neely GG, Hess A, Costigan M, Keene AC, Goulas S, et al. (2010) A genome-wide Drosophila screen for heat nociception identifies α2δ3 as an evolutionarily conserved pain gene. Cell 143: 628–638.
- 67. Pirone A, Kurt S, Zuccotti A, Rüttiger L, Pilz P, et al. (2014) α2δ3 is essential for normal structure and function of auditory nerve synapses and is a novel candidate for auditory processing disorders. J Neurosci 34: 434–445.
- 68. Huh JW, Kim YH, Lee SR, Kim H, Kim DS, et al. (2009) Gain of new exons and promoters by lineage-specific transposable elements-integration and conservation event on CHRM3 gene. Mol Cells 28: 111–117.
- 69. Lin HJ, Wan L, Chen WC, Lin JM, Lin CJ, et al. (2012) Muscarinic acetylcholine receptor 3 is dominant in myopia progression. Invest Ophthalmol Vis Sci 53: 6519–6525.
- 70. Krug A, Krach S, Jansen A, Nieratschker V, Witt SH, et al. (2013) The effect of neurogranin on neural correlates of episodic memory encoding and retrieval. Schizophr Bull 39: 141–150.