Microarray analysis of expression of cell death-associated genes in rat spinal cord cells exposed to cyclic tensile stresses in vitro
© Uchida et al; licensee BioMed Central Ltd. 2010
Received: 7 December 2009
Accepted: 22 July 2010
Published: 22 July 2010
The application of mechanical insults to the spinal cord results in profound cellular and molecular changes, including the induction of neuronal cell death and altered gene expression profiles. Previous studies have described alterations in gene expression following spinal cord injury, but the specificity of this response to mechanical stimuli is difficult to investigate in vivo. Therefore, we have investigated the effect of cyclic tensile stresses on cultured spinal cord cells from E15 Sprague-Dawley rats, using the FX3000® Flexercell Strain Unit. We examined cell morphology and viability over a 72 hour time course. Microarray analysis of gene expression was performed using the Affymetrix GeneChip System®, where categorization of identified genes was performed using the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) systems. Changes in expression of 12 genes were validated with quantitative real-time reverse transcription polymerase chain reaction (RT-PCR).
The application of cyclic tensile stress reduced the viability of cultured spinal cord cells significantly in a dose- and time-dependent manner. Increasing either the strain or the strain rate independently was associated with significant decreases in spinal cord cell survival. There was no clear evidence of additive effects of strain level with strain rate. GO analysis identified 44 candidate genes which were significantly related to "apoptosis" and 17 genes related to "response to stimulus". KEGG analysis identified changes in the expression levels of 12 genes of the mitogen-activated protein kinase (MAPK) signaling pathway, which were confirmed to be upregulated by RT-PCR analysis.
We have demonstrated that spinal cord cells undergo cell death in response to cyclic tensile stresses, which were dose- and time-dependent. In addition, we have identified the up regulation of various genes, in particular of the MAPK pathway, which may be involved in this cellular response. These data may prove useful, as the accurate knowledge of neuronal gene expression in response to cyclic tensile stress will help in the development of molecular-based therapies for spinal cord injury.
Mechanical stresses applied to the spinal cord can potentially induce profound and irreversible paresis, secondary to induced pathological changes such as dysfunction and loss of neurons, impairment of neuronal cell survival mechanisms and protein synthesis, neuronal cell necrosis and apoptosis [1, 2]. Examples of such mechanically induced spinal cord damage include not only spinal cord compression but distraction insult [3, 4]; however, it is likely that tensile stresses form an important part of many injuries of the spinal cord. The primary mechanical event, which may occur in less than a second, can initiate a cascade of molecular and cellular events such as changes in gene expression, which may then influence cell function over minutes to hours or a much longer period. For example, transient disruption of Ca2+ homeostasis may be an early event in a series of aberrant signaling cascades that ultimately lead to cellular dysfunction or cell death. Extended consequences of this molecular cascade include changes in gene expression levels that are necessary for cell recovery or cell death [5–8].
Changes in the expression of several immediate early response genes have been documented in various in vivo models of spinal cord injury using microarray analysis [9–11]. These included the up regulation of transcription factors, suggesting that the expression of many other genes is potentially regulated after traumatic insult. Indeed, the differential and post-traumatic expression levels of several genes have been explored in vivo in an attempt to stabilize, both biologically and functionally, the spinal cord once injured [12, 13]. However, these in vivo experimental settings for studying the response of neuronal systems to mechanical injury suffer from several disadvantages over in vitro experimentation. For example, mechanical stress can cause additional or unexpected tissue or cell reactions such as activation of resident inflammatory cells or invasion of foreign cells from the periphery . The complexity of the in vivo situation may also result in a limited accessibility to specific areas of tissue or cell type of interest, preventing real-time and spatial measurement of biological or mechanical parameters . Thus, in vitro models of the spinal cord stimuli can be useful to gain a better understanding of the specific neuronal response to mechanical stress.
One approach to determine the pathophysiology of mechanical-stress-related spinal cord damage is to investigate the in vitro response of neuronal cells to loading. The use of neuronal cell culture models allows for better control of the extracellular environment, is relative easy to manipulate, and permits for repeated access to neural cells for specific analysis. The spinal cord and neurons are always subjected to mechanical stress including tensile stresses, during spine movement. Longitudinal vertebral distraction and the physiological tension zone  of the spinal cord are closely correlated each other when the spine is subjected to flexural positioning [17, 18] and excessive kyphosis in the thoracic vertebrae . The Flexercell Strain Unit (FX3000®, Flexercell International, Hillsborough, NC) is a cell-stretching apparatus that allows application of cyclic tensile force to cultured cells. The system has been used to elucidate the mechanism of mechanical signaling in various types of cells [20–22]. In our previous study using this equipment , we investigated the in vitro effects of cyclic tensile stress on cultured spinal cord cells, with a special focus on the expressions of neurotrophins and their receptor genes. The results of that study showed that the application of tensile stress increased the expression levels of nerve growth factor, brain-derived neurotrophic factor, trkB, p75 neurotrophin receptor (p75NTR), glial cell line-derived neurotrophic factor, and caspase-9 mRNAs in the acute phase, followed by increased lactate dehydogenase release and induction of necrotic cell death. These findings led us to investigate further the expression of several genes related to cell death in cultured spinal cord cells under cyclic tensile stress.
The present study was thus designed to examine further the molecular changes and gene expression profiles in cultured spinal cord cells using the above cell-stretching apparatus and DNA microarray technology in order to provide a more complete picture of the changes in the expression of specific genes involved in neuronal response to cyclic tensile stress.
Cyclic tensile stress induced spinal cord cell death
Cluster analysis of gene expression profiles
Identification of upregulated genes through Gene Ontology (GO) analysis
List of "apoptosis" genes selectively upregulated in cluster 3
Entrez Gene No.
nuclear receptor subfamily 4, group A, member 2
fos-like antigen 1
integrin alpha V
twist homolog 2 (Drosophila)
retinoic acid receptor, beta
dual specificity phosphatase 1
heme oxygenase (decycling) 1
apoptotic chromatin condensation inducer 1
Platelet derived growth factor receptor, beta polypeptide
glutamate-cysteine ligase, catalytic subunit
inositol hexaphosphate kinase 2
death associated protein kinase 1
adhesion molecule with Ig like domain 2
matrix metallopeptidase 2
Homeodomain interacting protein kinase 2
baculoviral IAP repeat-containing 3
myeloid differentiation primary response gene 116
interleukin 1 receptor antagonist
homeodomain interacting protein kinase 1
similar to death effector domain-containing DNA binding protein 2
receptor (TNFRSF)-interacting serine-threonine kinase 2
death associated protein kinase 1
Tissue inhibitor of metalloproteinase 3
ceroid lipofuscinosis, neuronal 3, juvenile (Batten, Spielmeyer-Vogt disease)
excision repair cross-complementing rodent repair deficiency, complementation group 5
phosphatidylinositol 3-kinase, catalytic, alpha polypeptide
v-raf-leukemia viral oncogene 1
endothelin receptor type B
homeodomain interacting protein kinase 2
Tissue inhibitor of metalloproteinase 3
vascular endothelial growth factor A
furin (paired basic amino acid cleaving enzyme)
E1A binding protein p300
SH3-domain kinase binding protein 1
PTK2 protein tyrosine kinase 2 beta
similar to Tnf receptor-associated factor 1
neutrophil cytosolic factor 1
amyloid beta (A4) precursor protein-binding, family B, member 2
ceroid-lipofuscinosis, neuronal 8
SH3-domain kinase binding protein 1
List of "response to stimulus" genes selectively upregulated in cluster 1
Entrez Gene No.
serine (or cysteine) peptidase inhibitor, clade A, member 3N
prostaglandin-endoperoxide synthase 2
chemokine (C-C motif) ligand 20
chemokine (C-X-C motif) ligand 3
heat shock protein 4 like
nuclear protein 1
crystallin, alpha B
cytochrome P450, family 4, subfamily b, polypeptide 1
chemokine (C-X-C motif) ligand 10
solute carrier family 12 (sodium/potassium/chloride transporters), member 2
Identification of upregulated genes using Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis
High frequency pathways* identified by KEGG/pathway analysis in clusters 1, 2, and 3
pathway (total gene count)
Glutamate metabolism (26)
D-Glutamine and D-glutamate metabolism (3)
Taurine and hypotaurine metabolism (10)
Pantothenate and CoA biosynthesis (11)
T cell receptor signaling pathway (95)
Glycosphingolipid biosynthesis - ganglio series (15)
Cell adhesion molecules (CAMs) (151)
Focal adhesion (187)
Bladder cancer (36)
TGF-beta signaling pathway (81)
Pathways in cancer (313)
Prostate cancer (92)
Small cell lung cancer (92)
ECM-receptor interaction (74)
Renal cell carcinoma (69)
PPAR signaling pathway (70)
mTOR signaling pathway (54)
Gap junction (90)
Glycerophospholipid metabolism (46)
MAPK signaling pathway (255)
Non-small cell lung cancer (52)
Pancreatic cancer (71)
GnRH signaling pathway (91)
Adherens junction (73)
Pyruvate metabolism (38)
Regulation of actin cytoskeleton (204)
Leukocyte transendothelial migration (116)
Chronic myeloid leukemia (79)
Glutathione metabolism (42)
Gene-specific real-time reverse transcription polymerase chain reaction (RT-PCR) in MAPK signaling pathway
In vitro systems offer several advantages over in vivo systems in the analysis of cellular responses to their mechanical environment, including the precise specification of loading parameters (strain levels and strain rate), control of the extracellular environment (temperature, ion concentration, partial pressure of gases), the relative ease and repeated access of the cells, and simplified administration of pharmacological compounds. By precisely controlling the loading conditions, the quantitative relationship between the severity of mechanical injury and response to the injury can be evaluated . Previous mechanical stimulation studies of neuronal cells using immortalized cell lines reported that the physiological strain conditions of neuronal cells are not as high as osteoblasts or skeletal muscle cells . Rat motoneurons stretched at 11% strain showed ischemic changes but no mechanical damage , whereas at a 6% strain level there were no ischemic changes or mechanical damage . PC12 cells were subjected to cyclic tensile strain levels ranging from 4 to 16% at strain frequencies of 1-2 Hz as physiological mechanical conditions . Based on these reports and considering primary cultured spinal cord cells , we selected the range of cyclic tensile stresses most appropriate to our culture system in the present study, and the cells were observed morphologically following application of various tensile strains levels (5%, 10%, and 15%) and strain frequencies (0.5 Hz and 1 Hz) using the Flexercell system. Interestingly, we found a significant decrease in the cell survival rate (%) when cultures were subjected to either a cyclic tensile stress of 5% at 1 Hz or 15% at 0.5 Hz, compared with a cyclic tensile stress at 10% at 0.5 Hz, which was taken as a control level. Our results suggest that both the level of strain applied and frequency of its application influence cell viability. Furthermore the results demonstrate that a higher strain level at a lower strain rate can have a similar effect as a lower strain level at a higher strain rate in neuron-rich spinal cord cells.
To understand the molecular mechanism of neuronal responses to cyclic tensile stress, DNA microarray was employed to identify the specific gene expression patterns in cultured spinal cord cells, and the hierachical clustering algorithm was used in the analysis. Consequently, we found that 44 genes related to "apoptosis" among cluster 3 included the majority of progressively upregulated genes by enrichment analysis using biological process in the GO system. Furthermore, 17 genes related to "response to stimulus" were also identified among those of cluster 1, including most upregulated genes in the late phase after application of cyclic tensile stress. On the other hand, the KEGG/pathway analysis identified different pathways in the upregulated genes in a time-dependent manner in the clusters 1, 2, and 3 which had been identified via the GO system.
A number of investigators have attempted to characterize or determine changes in specific genes and signaling pathways of direct mechanical stress or load using cells from diverse tissues and in vitro models of mechanical trauma. These systems are multifold and can be separated into detection, response, and modulation pathways . A load stimulus is detected by a mechano-electrochemical sensory system, including mechanically sensitive ion channels (stretch-activated channels) linked to the cytoskeleton . The principle detection system is the matrix-integrin-mechanosensory protein complex-cytoskeleton machinery [30, 31]. This system is linked to the kinase cascade (tyrosine or nontyrosine kinase cascade or the JAK/STAT kinase cascade) system, which provides amplification, diversity, selectivity, and modulation capabilities and contains focal adhesion kinase (FAK) . The Src homology protein complex (SHC) and growth factor receptor binding adaptor protein linking receptor (GRB2) link to the Ras-raf signaling pathway, and these are regulated through FAK. Ras and GTPase activation regulates the activation of Raf, MAPK and extracellular regulated kinase (ERK) [33, 34]. Activated ERKs enter the nucleus and upregulate the expression of various transcription factors (jun, fos, myc, erg-1) and activate nuclear binding proteins such as NF-κB . Specific regulation may occur at load response elements in promoters of certain genes such as platelet-derived growth factor (PDGF) . The mechanically sensitive ion channel could also be a gap junction channel, which propagates a Ca2+ wave from cell to cell after a mechanical signal is detected through a cell adhesion molecule (CAM).
It is clear that systems involved in evoking the response of neuronal cells to applied mechanical load are highly complex, which requires an expansive approach to their elucidation. Our results in primary neuronal-rich culture have analyzed gene expression profiles using the Affymatrix system and thereby have identified 3,412 genes with altered expression as a result of applied cyclic tensile stress, many of which are consistent with previous studies. These data will be useful for future examination of how the effects of mechanical insult on the spinal cord, in particular, are regulated. Signalling pathways that may be involved in this process can be modulated by molecules such as G proteins, kinase phosphorylation and activation, and kinase inhibitors or phosphatases . Indeed, the MAPK kinase pathway was identified by KEGG analysis of the alterations in gene expression seen in the mechanically loaded spinal cord cells. MAPKs are a family of related serine/threonine protein kinases that transduce several signals responsible for cell proliferation or cellular stress, and are also intracellular signaling systems that induce optimal stress responses. In general, MAPK is activated by phosphorylation of tyrosine and threonine residues by MAPKK, which is activated through phosphorylation by MAPKKK. In KEGG analysis, the MAPK signaling pathway can be classified into three main groups : the classical MAP pathway, the c-Jun N-terminal kinase or stress-activated protein kinase (JNK), the p38 MAPK signaling pathway, and the ERK pathway. In the present study, the mRNA expression levels of PDGFR and G12 significantly increased during the mid period of cyclic tensile stress application, while CACN, Raf1, MKP, JIP3, HSP72, GAD D45, and GAD D153 mRNA levels increased during the late phase of the cyclic tensile stress application. These genes, which were further identified by real time RT-PCR analysis, may play an important role in the response of spinal cord cells to neuronal injury. Previous in vivo studies of spinal cord injuries using microarray have also demonstrated the participation of similar genes in acute phage after injury [9, 13]. Gene expression of trkA/B, FGFR, PDGFR, and Raf1 could facilitate neuronal survival, while gene expression of G12, HSP72, GAD D45, GAD D153 could be involved in DNA damage. Further studies that target specific gene expression pathways will help determine their precise role in neuronal responses to mechanical load.
Our study is perhaps the first to attempt investigating the expression of genes related to apoptotic cell death during the application of cyclic tensile stress to neuronal rich spinal cord cells; however, it has several limitations. These include (i) imperfect validity of microarray results, with a CV of 5-15% for quantitative signals ; (ii) the neuronal culture was not 100% pure (the rate of NeuN-positive cells was 71%), and embryonic cells were used, not adults cells; and (iii) the lack of immunohistological or in situ hybridization data to identify the exact types of cells in which alterations in gene expression occurred. Nevertheless, we believe that the present study might be the first to comprehensively profile changes in gene expression involved in neuronal responses to cyclic tensile stress in cultured spinal cord cells using DNA microarray. This is a considerable improvement in examining the specific response of neuronal cells to mechanical load. Considered together with our previous findings , we can conclude that certain apoptosis-specific genes are activated in neuronal cell rich cultures during the application of cyclic tensile stress. The clinical relevance of tensile stress may specifically include the tethering effect with the developmental ascensus medullaris , cervical myelopathy in association with kyphotic deformity , and complicated spinal cord distraction injury. However, it is perhaps intuitive to consider that abnormal tensile stresses are involved in many mechanical insults of the spinal cord. Thus, we believe that our study provides new insights into the pathophysiology of spinal cord damage in various disease entities. Furthermore, the current study may help understand the response of neuronal cells to cyclic tensile stress and therapeutic issues related to the mechanically damaged spinal cord.
We have investigated the effects of cyclic tensile stresses on cultured spinal cord cells and demonstrate that cell death was induced depending on the level and duration of strain applied. Furthermore, we have performed a comprehensive analysis of alterations in gene expression profiles that occur following this mechanical stress, and identified in particular an upregulation of members of the MAPK pathway. Knowledge of the specific response of neuronal cells to mechanical insult could be a potentially useful tool for molecular-based therapy for spinal cord injury.
Cell isolation and culture
Primary cultures were established using the method described previously by our group . In brief, the spinal cords of Sprague-Dawley rat embryos were dissected out at post-coital day 15 and stripped of the dorsal root ganglia and meninges. Dissected tissues were rinsed with cold Ca2+- and Mg2+-free Hanks balanced salt solution (HBSS) supplemented with 4 g/L glucose, and incubated at 37°C for 20 minutes with 0.03% (w/v) trypsin solution in HBSS under mild shaking. They were transferred into HBSS containing 0.1% (w/v) soyabean trypsin inhibitor (Sigma, St. Louis, MO) and 0.2% (w/v) bovine serum albumin (BSA), and triturated very mildly. The cell suspension was filtered through nylon mesh (70 μm, Cell Strainer; Becton Dickinson, Bedford, MA). The culture medium consisted of 75 mL Leibovitz's L-15 medium supplemented with glucose (4 g/L), 1.0 mL N2 supplement, 15 mL 0.15 M sodium bicarbonate, 10 mL heat-inactivated horse serum, 1 mL of 100 mM L-cysteine and 1 mL penicillin G 104 U/mL and neutralized with CO2. After centrifugation at 400 × g for 15 minutes at 4°C, the precipitated cells were gently re-suspended in a fresh culture medium and plated at a density of 4.0 × 105 cells/well onto a 6-well culture plate with a flexible-polystyrene bottom coated with type IV collagen (BioFlex® Baseplate, Flexercell International).
The experiment was carried out in the Orthopaedic Spinal Cord Laboratory of our University Medical Faculty. The experimental protocol strictly followed the Ethics Review Committee Guidelines for Animal Experimentation of our University.
Application of cyclic tensile stress to the cultured spinal cord cells
The cell stretching device used in this study was the Flexercell FX3000® (Flexercell International). The device consists of a computer-controlled vacuum unit, a culture plate with a flexible-polystyrene well-bottom coated with type IV collagen (BioFlex® Baseplate), and another culture plate with a non-deformable culture well bottom constructed of the same material as control. The culture plates consisted of 6-well (35 mm diameter) flexible-bottomed culture plate with a hydrophilic surface. The application of a vacuum provides a hemispherically downward deforming force onto the flexible bottom, resulting in a non-homogenous strain profile with a maximum at the periphery and a geometric decline toward zero at the center of the culture well bottom. For these experiments, the cultured spinal cord cells were subjected to various conditions. The flexible-bottomed culture plates including the control plates were then placed on the vacuum baseplate in the incubator (37°C in 5% CO2). Three days after cell seeding, the cells were subjected to cyclic tensile stress for up to 72 hours. Previously we have documented  that ~71% of cells isolated from the spinal cord using this methodology are positive for the neuronal marker NeuN (Chemicon International, Temecula, CA). Repeated examinations by phase contrast microscopy showed that the cells remained attached to the substratum during elongation of the flexible-polystyrene well plates as tensile stresses were applied. Panjabi and White  demonstrated that the elastic properties of the cervical spinal cord are dramatically altered when it is subjected to approximately 10% elongation. In our previous in vivo study , the amplitude of epidurally-recorded spinal cord evoked potentials began to diminish, especially the second component (N2 spike), when the longitudinal extension of the cord shortened by 10-17%. Thus, in the present study, we set the tensile stress levels at a maximum 15% elongation and set the strain rate to estimate an appropriate frequency of spine movement in everyday life. Therefore, the cells were observed morphologically following application of various tensile stresses resulting in strains of 5%, 10%, and 15% applied at frequencies of 0.5 Hz and 1 Hz. Analyses of cell survival, DNA microarrays and real-time RT-PCR were conducted at 0, 2, 6, 12, 24, 48 and 72 hours after the application of the cyclic tensile load. DNA microarray analysis and real time RT-PCR were performed to compare the levels of gene expression at time 0 with levels of gene expression thereafter after the cyclic tensile loading at o.5 Hz, resulting in 10% strain.
Quantification of cell survival under cyclic tensile stress
Cell survival was investigated by scoring the number of living cells after tensile stress application using the Live/Dead Assay (Live/Dead Assay, Molecular Probes, Eugene, OR), according to the manufacturers' instruction. This assay if based on the differential staining of cells with calcein-acetoxymethyl ester (calcein AM: 3.00 μM) to identify living cells and ethidium homodimer-1 to identify dead cells [39, 40]. Calcein-AM is a membrane-permeable dye that is cleaved by intracellular esterase to produce an impermeant green-wavelength fluorophore in living cells. Ethidium homodimer-1 cannot penetrate live cells, but it can enter dead cells which have a porous membrane and hence bind to DNA to produce red fluorescence. The culture medium was removed and the cells were then washed twice with PBS, and stained for 75 minutes at 32°C. The numbers of attached living cells (green) in at least 6 high-power fields (each containing at least 100 cells) were counted using fluoromicroscopy (IX70, Olympus, Tokyo) and a color image analyzer (MacSCOPE, Mitani, Fukui, Japan) in more than three wells for each time point. There was no evidence of spinal cord cell proliferation during the 3 day period prior to treatment with cyclic tensile stress, i.e. cell counts were almost uniform and at a density between 3.3 × 105 and 4.8 × 105 cells/well after dissemination on Bioflex® Baseplate in the absence of mechanical stimuli. The cell survival rate (%) at each time point for cultures which were subjected to 10% strain at 0.5 Hz frequency was calculated relative to the cell number at 0-hour. These cultures (subjected to 10% tensile strain at 0.5 Hz frequency) were then set as a standard, to which the cell viability of other levels of strain and frequency were compared.
All values were expressed as mean ± standard error of the mean (SEM). Differences between values of the loaded and control cultures were tested at each point by one-way ANOVA and Tukey posthoc test using the SPSS software version 11.0 (SPSS, Chicago, IL). P values of less than 0.05 denoted the presence of a statistically significant difference.
The presence of DNA fragmentation was examined via TEM examination. After application of tensile stresses, cultured spinal cord cells were washed twice with PBS and fixed with 2.5% glutaraldehyde and 2.5% paraformaldehyde, followed by late fixation in 1% osmium tetroxide for 2 hours. Fixed specimens were dehydrated in a graded series of alcohol, embedded in epoxy resin and polymerized at 60°C for 2 days. Ultrathin sections were cut by ultramicrotome, stained with uranyl acetate and lead citrate, and examined with a Hitachi H-7650 TEM (Hitachi, Tokyo).
RNA preparation and DNA microarray hybridization
The cultured cells on each well at 0, 2, 6, 12, 24, 48 and 72 hours were disrupted in a lysis buffer containing β-mercaptoethanol and the total RNA from five animals was pooled at each time point, and further purified using RNeasy® Mini Kit (Qiagen, Valencia, CA) and treated with DNase I (TaKaRa Bio, Ohtsu, Japan). The quality of RNA was initially assessed by electrophoresis on a 1.5% agarose gel, and further determined by using the RNA 6000 Nano LabChip Kit® and Agilent Bioanalyzer 2100 (Agilent, Palo Alto, CA). cDNAs were synthesized by GeneChip T7-Oligo (dT) Promoter Primer Kit® (Affymetrix, Santa Clara, CA) and TaKaRa cDNA Synthesis Kit® (TaKarRa Bio) from 3 μg total RNA. Biotinylated cRNA were synthesized using the IVT Labeling Kit (Affymetrix). Following fragmentation, 20 μg of cRNA were hybridized for 16 hours at 45°C on the GeneChip® Rat Genome 230 2.0 Array. GeneChips were washed and stained in the Affymetrix Fluidics Station 450, and scanned using GeneChip Scanner 3000 7G.
Data analysis of GeneChip® expression array
Microarray data was initially processed using GeneChip Operating Software (GCOS, Affymetrix). Signal intensity was calculated by Microarray Suite version 5.0 (MAS5.0) with Affymetrix default setting and global scaling as the normalization method. The trimmed mean target intensity of each array was arbitrarily set to 500. The ratio and difference of intensity between two corresponding genes on each array was calculated. Significance of the difference in levels of gene expression was set at a ratio and difference of intensity of ≧2 and ≧200, respectively. Spots that could not be interpreted were excluded, resulting in a list of 3,412 genes available for subsequent analysis. After filtering, the hierarchical clustering was applied to the axis using the complete linkage method, as implemented in the program Spotfire 8.1.1. The distance matrices used were uncentered, with correlation calculated on the basis of expression signals for clustered genes.
GO and KEGG analysis
where N is the total number of genes, M is the number of genes that annotated for a given term or pathway, n is the size of the list of genes of interest and k is the number of genes within the list that annotated for a given term or pathway.
Microarray including GeneChip® loads some different prove sequences in order to identify one gene. Strictly speaking, microarray data based on signal intensity represents the loading common probe set domain designed by Affymetrix. Therefore, it is possible that the number of gene symbols is different from that of the identified gene by the probe set.
Real-time RT-PCR analysis
Sequences of primers used for real-time PCR
PCR Product Size (bp)
Sequence Accession No.
We are indebted to Yoshimasa Tsujimoto (Dragon Genomics Center, TAKARA BIO INC.) for support of analysis the data. This work was supported in part by Grants-in-Aid to HB, KU and HN for General Scientific Research of the Ministry of Education, Science and Culture of Japan (grants numbers 16390435, 18390411, 19791023, 21591895, 21791389, and 22390287). This work was also supported in part by grants from the Investigation Committees on Ossification of the Spinal ligaments (2005-2009) and Research Program of Spinal Cord Intractable Pain (2009), The Public Health Bureau of the Japanese Ministry of Health and Welfare.
- Baptiste DC, Fehling MG: Pathophysiology of cervical myelopathy. Spine J. 2006, 6 (6 Suppl): 190S-197S. 10.1016/j.spinee.2006.04.024.View ArticlePubMedGoogle Scholar
- Uchida K, Baba H, Maezawa Y, Furukawa S, Omiya M, Kokubo Y, Kubota C, Nakajima H: Increased expression of neurotrophins and their receptors in the mechanically compressed spinal cord of the spinal hyperostotic mouse (twy/twy). Acta Neuropathol. 2003, 106: 29-36.PubMedGoogle Scholar
- Ghafoor AU, Martin TW, Gopalakrishnan S, Viswamitra S: Caring for the patients with cervical spine injuries: what have we learned?. J Clin Anesth. 2005, 17: 640-649. 10.1016/j.jclinane.2005.04.003.View ArticlePubMedGoogle Scholar
- Uchida K, Nakajima H, Sato R, Yayama T, Erisa S, Mwaka ES, Kobayashi S, Baba Hisatoshi: Cervical spondylotic myelopathy associated with kyphosis or sagittal sigmoid alignment: outcome after anterior or posterior decompression. J Neurosurg Spine. 11: 521-528.
- Hayes RL, Yang K, Raghupathi R, McIntosh TK: Changes in gene expression following traumatic brain injury in the rat. J Neurotrauma. 1995, 12: 779-790. 10.1089/neu.1995.12.779.View ArticlePubMedGoogle Scholar
- Raghupathi R, McIntosh TK: Regionally and temporally distinct patterns of induction of c-fos, c-jun and junB mRNAs following experimental brain injury in the rat. Brain Res Mol Brain Res. 1996, 37: 134-144. 10.1016/0169-328X(95)00289-5.View ArticlePubMedGoogle Scholar
- Sall JM, Morehead M, Murphy S, Goldman H, Walker PD: Alterations in CNS gene expression in a rodent model of moderate traumatic brain injury complicated by acute alcohol intoxication. Exp Neurol. 1996, 139: 257-268. 10.1006/exnr.1996.0099.View ArticlePubMedGoogle Scholar
- Yakovlev AG, Faden AI: Molecular biology of CNS injury. J Neurotrauma. 1995, 12: 767-777. 10.1089/neu.1995.12.767.View ArticlePubMedGoogle Scholar
- Carmel JB, Galante A, Soteropoulos P, Tolias P, Recce M, Young W, Hart RP: Gene expression profiling of acute spinal cord injury reveals spreading inflammatory signals and neuron loss. Physiol Genomic. 2001, 7: 201-213.View ArticleGoogle Scholar
- Song G, Cechvala C, Resnick DK, Dempsey RJ, Rao VL: GeneChip analysis after acute spinal cord injury in rat. J Neurochem. 2001, 79: 804-815. 10.1046/j.1471-4159.2001.00626.x.View ArticlePubMedGoogle Scholar
- Tachibana T, Noguchi K, Ruda MA: Analysis of gene expression following spinal cord injury in rat using complementary DNA microarray. Neurosci Lett. 2002, 327: 133-137. 10.1016/S0304-3940(02)00375-0.View ArticlePubMedGoogle Scholar
- Aimone JB, Leasure JL, Perreau VM, Thallmair M: Spatial and temporal gene expression profiling of the contused rat spinal cord. Exp Neurol. 2004, 189: 204-221. 10.1016/j.expneurol.2004.05.042.View ArticlePubMedGoogle Scholar
- Di Giovanni S, Knoblach SM, Brandoli C, Aden SA, Hoffman EP, Faden AI: Gene profiling in spinal cord injury shows role of cell cycle in neuronal death. Ann Neurol. 2003, 53: 454-468. 10.1002/ana.10472.View ArticlePubMedGoogle Scholar
- Popovich PG, Strokes BT, Whitacre CC: Concept of autoimmunity following spinal cord injury: possible roles for T lymphocytes in the traumatized central nerve system. J Neurosci Res. 1996, 45: 349-363. 10.1002/(SICI)1097-4547(19960815)45:4<349::AID-JNR4>3.0.CO;2-9.View ArticlePubMedGoogle Scholar
- Morrison B, Saatman KE, Meaney DF, McIntosh TK: In vitro central nervous system models of mechanically induced trauma. A review. J Neurotrauma. 1998, 15: 911-928. 10.1089/neu.1998.15.911.View ArticlePubMedGoogle Scholar
- Panjabi MM, White AA: Physical properties and functional biomechanics of the spine: the spinal cord. Clinical Biomechanics of the Spine. Edited by: White AA III, Panjabi MM. 1996, Philadelphia: J.B. Lippincott, 67-71. 2Google Scholar
- Kawahara N, Baba H, Nagata S, Kikuchi Y, Tomita K, Nomura S, Yugami H: Experimental studies on the spinal cord evoked potentials in cervical spine distraction injuries. Spinal Cord Monitoring and Electrodiagnosis. Edited by: Shimoji K, Kurokawa T, Tamaki T, Willis WD Jr. 1991, Tokyo: Springer-Verlag, 107-115.View ArticleGoogle Scholar
- Bilston LE, Thibault LE: The mechanical properties of the human cervical spinal cord in vitro. Ann Biomed Eng. 1996, 24: 67-74. 10.1007/BF02770996.View ArticlePubMedGoogle Scholar
- Shimizu K, Nakamura M, Nishikawa Y, Nishikawa Y, Hijikata S, Chiba K, Toyama Y: Spinal kyphosis causes demyelination and neuronal loss in the spinal cord: a new model of kyphotic deformity using juvenile Japanese small game fowls. Spine. 2005, 30: 2388-2392. 10.1097/01.brs.0000184378.67465.5c.View ArticlePubMedGoogle Scholar
- Gilbert JA, Weinhold PS, Banes AJ, Link GW, Jones GL: Strain profiles for circular cell culture plates containing flexible surfaces employed to mechanically deform cells in vitro. J Biomech. 1994, 27: 1169-1177. 10.1016/0021-9290(94)90057-4.View ArticlePubMedGoogle Scholar
- Matsumoto T, Kawakami M, Kuribayashi K, Takenaka T, Tamaki T: Cyclic mechanical stretch stress increases the growth rate and collagen synthesis of nucleus pulposus cells in vitro. Spine. 1999, 24: 315-319. 10.1097/00007632-199902150-00002.View ArticlePubMedGoogle Scholar
- Nakatani T, Marui T, Hirota T, Doita M, Nishida K, Kurosaka M: Mechanical stretching force promotes collagen synthesis by cultured cells from human ligamentum flavum via transforming growth factor-β1. J Orthop Res. 2002, 20: 1380-1386. 10.1016/S0736-0266(02)00046-3.View ArticlePubMedGoogle Scholar
- Uchida K, Nakajima H, Takamura T, Furukawa S, Kobayashi S, Yayama T, Baba H: Gene expression profiles of neurotrophic factors in rat cultured spinal cord cells under cyclic tensile stress. Spine. 2008, 33: 2596-2604. 10.1097/BRS.0b013e31818917af.View ArticlePubMedGoogle Scholar
- Morrison B, Meaney DF, Margulies SS, McIntosh TK: Dynamic mechanical stretch of organotypic brain slice cultures induces differential genomic expression: relationship to mechanical parameters. J Biomech Eng. 2000, 122: 224-230. 10.1115/1.429650.View ArticlePubMedGoogle Scholar
- Haq F, Keith C, Zhang G: Neurite development in PC12 cells on flexible micro-textured substrates under cyclic stretch. Biotechnol Prog. 2006, 22: 133-140. 10.1021/bp0501625.View ArticlePubMedGoogle Scholar
- Morrison B, Cater HL, Wang CC, Thomas FC, Hung CT, Ateshian GA, Sundstrom LE: A tissue level tolerance criterion for living brain developed with an in vitro model of traumatic mechanical loading. Stapp Car Crash J. 2003, 47: 93-105.PubMedGoogle Scholar
- Tanoue M, Yamaga M, Ide J, Takagi K: Acute stretching of peripheral nerves inhibits retrograde axonal transport. J Hand Surg [Br]. 1996, 21: 358-363.View ArticleGoogle Scholar
- Banes AJ, Tsuzaki M, Yamamoto J, Fischer T, Brigman B, Brown T, Miller L: Mechanoreception at the cellular level; The detection, interpretation, and the diversity of responses to mechanical signals. Biochem Cell Biology. 1995, 73: 349-365. 10.1139/o95-043.View ArticleGoogle Scholar
- Guharay F, Sachs F: Stretch-activated single ion channel currents in tissue-cultured embryonic chick skeletal muscle. J Physiol. 1984, 352: 685-701.PubMed CentralView ArticlePubMedGoogle Scholar
- Burridge K, Fath K, Kelly T, Nuckolls G, Turner C: Focal adhesions: transmembrane junctions between the extracellular matrix and the cytoskeleton. Ann Rev Cell Biol. 1988, 4: 487-525.View ArticlePubMedGoogle Scholar
- Clark EA, Brugge JS: Integrins and signal transduction pathways: the road taken. Science. 1995, 268: 233-239. 10.1126/science.7716514.View ArticlePubMedGoogle Scholar
- Shyy JY, Chien S: Role of integrins in cellular responses to mechanical stress and adhesion. Curr Opin Cell Biol. 1997, 9: 707-713. 10.1016/S0955-0674(97)80125-1.View ArticlePubMedGoogle Scholar
- Li C, Hu Y, Mayr M, Xu Q: Cyclic strain stress-induced mitogen-activated protein kinase (MAPK) phosphatase 1 expression in vascular smooth muscle cells is regulated by Ras/Rac-MAPK pathways. J Biol Chem. 1999, 274: 25273-25280. 10.1074/jbc.274.36.25273.View ArticlePubMedGoogle Scholar
- Li C, Xu Q: Mechanical stress-initiated signal transductions in vascular smooth muscle cells. Cell Signal. 2000, 12: 435-445. 10.1016/S0898-6568(00)00096-6.View ArticlePubMedGoogle Scholar
- Copland IB, Post M: Stretch-activated signaling pathways responsible for early response gene expression in fetal lung epithelial cells. J Cell Physiol. 2007, 210: 133-143. 10.1002/jcp.20840.View ArticlePubMedGoogle Scholar
- Resnick N, Collins T, Atkinson W, Bonthron DT, Dewey CF, Gimbrone MA: Platelet-derived growth factor B chain promoter contains a cis-acting fluid shear-stress-responsive element. Proc Natl Acad Sci USA. 1993, 90: 4591-4595. 10.1073/pnas.90.10.4591.PubMed CentralView ArticlePubMedGoogle Scholar
- Davis RJ: The mitogen-activated protein kinase signal transduction pathway. J Biol Chem. 1993, 268: 14553-14556.PubMedGoogle Scholar
- Shi L, Reid LH, Jones WD, Shippy R, Warrington JA, Baker SC, Collins PJ, de Longueville F, Kawasaki ES, Lee KY, Luo Y, Sun YA, Willey JC, Setterquist RA, Fischer GM, Tong W, Dragan YP, Dix DJ, Frueh FW, Goodsaid FM, Herman D, Jensen RV, Johnson CD, Lobenhofer EK, Puri RK, Schrf U, Thierry-Mieg J, Wang C, Wilson M, Wolber PK, Zhang L, Amur S, Bao W, Barbacioru CC, Lucas AB, Bertholet V, Boysen C, Bromley B, Brown D, Brunner A, Canales R, Cao XM, Cebula TA, Chen JJ, Cheng J, Chu TM, Chudin E, Corson J, Corton JC, Croner LJ, Davies C, Davison TS, Delenstarr G, Deng X, Dorris D, Eklund AC, Fan XH, Fang H, Fulmer-Smentek S, Fuscoe JC, Gallagher K, Ge W, Guo L, Guo X, Hager J, Haje PK, Han J, Han T, Harbottle HC, Harris SC, Hatchwell E, Hauser CA, Hester S, Hong H, Hurban P, Jackson SA, Ji H, Knight CR, Kuo WP, LeClerc JE, Levy S, Li QZ, Liu C, Liu Y, Lombardi MJ, Ma Y, Magnuson SR, Maqsodi B, McDaniel T, Mei N, Myklebost O, Ning B, Novoradovskaya N, Orr MS, Osborn TW, Papallo A, Patterson TA, Perkins RG, Peters EH, Peterson R, Philips KL, Pine PS, Pusztai L, Qian F, Ren H, Rosen M, Rosenzweig BA, Samaha RR, Schena M, Schroth GP, Shchegrova S, Smith DD, Staedtler F, Su Z, Sun H, Szallasi Z, Tezak Z, Thierry-Mieg D, Thompson KL, Tikhonova I, Turpaz Y, Vallanat B, Van C, Walker SJ, Wang SJ, Wang Y, Wolfinger R, Wong A, Wu J, Xiao C, Xie Q, Xu J, Yang W, Zhang L, Zhong S, Zong Y, Slikker W: The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006, 24: 1151-1161. 10.1038/nbt1239.View ArticlePubMedGoogle Scholar
- Takeno K, Kobayashi S, Negoro K, Uchida K, Miyazaki T, Yayama T, Shimada S, Baba H: Physical limitations to tissue engineering of intervertebral disc cells: effect of extracellular osmotic change on glycosaminoglycan production and cell metabolism. Laboratory investigation. J Neurosurg Spine. 2007, 7: 637-644. 10.3171/SPI-07/12/637.View ArticlePubMedGoogle Scholar
- Vaughan PJ, Pike CJ, Cotman CW, Cunningham DD: Thrombin receptor activation protects neurons and astrocytes from cell death produced by environmental insults. J Neurosci. 1995, 15: 5389-5401.PubMedGoogle Scholar