Molecular Vision 2008; 14:1050-1063 <http://www.molvis.org/molvis/v14/a126>
Received 03 January 2008 | Accepted 23 April 2008 | Published 03 June 2008

Time course profiling of the retinal transcriptome after optic nerve transection and optic nerve crush

Marta Agudo,1,2 Maria Cruz Pérez-Marín,1 Ulrika Lönngren,3 Paloma Sobrado,1 Ana Conesa,4 Isabel Cánovas,1 Manuel Salinas-Navarro,1 Jaime Miralles-Imperial,2 Finn Hallböök,3 Manuel Vidal-Sanz1

The second and third authors contributed equally to this work.

1Departamento Oftalmología, Facultad de Medicina, Universidad de Murcia, Murcia Spain; 2Hospital Universitario Virgen de la Arrixaca, Murcia, Spain; 3Department of Neuroscience, Unit for Developmental Neuroscience, Uppsala University, Uppsala, Sweden; 4Bioinformatics Department Centro de Investigación Príncipe Felipe, Valencia, Spain

Correspondence to: Marta Agudo, Departamento de Oftalmología, Facultad Medicina,Campus Espinardo, Universidad de Murcia. 30100 Murcia, Spain; Phone: +34968363996; FAX: +34968363962; email: martabar@um.es

Abstract

Purpose: A time-course analysis of gene regulation in the adult rat retina after intraorbital nerve crush (IONC) and intraorbital nerve transection (IONT).

Methods: RNA was extracted from adult rat retinas undergoing either IONT or IONC at increasing times post-lesion. Affymetrix RAE230.2 arrays were hybridized and analyzed. Statistically regulated genes were annotated and functionally clustered. Arrays were validated by means of quantative reverse transcription polymerase chain reaction (qRT–PCR) on ten regulated genes at two times post-lesion. Western blotting and immunohistofluorescence for four pro-apoptotic proteins were performed on naïve and injured retinas. Finally, custom signaling maps for IONT- and IONC-induced death response were generated (MetaCore, Genego Inc.).

Results: Here we show that over time, 3,219 sequences were regulated after IONT and 1,996 after IONC. Out of the total of regulated sequences, 1,078 were commonly regulated by both injuries. Interestingly, while IONT mainly triggers a gene upregulation-sustained over time, IONC causes a transitory downregulation. Functional clustering identified the regulation of high interest biologic processes, most importantly cell death wherein apoptosis was the most significant cluster. Ten death-related genes upregulated by both injuries were used for array validation by means of qRT–PCR. In addition, western blotting and immunohistofluorescence of total and active Caspase 3 (Casp3), tumor necrosis factor receptor type 1 associated death domain (TRADD), tumor necrosis factor receptor superfamily member 1a (TNFR1a), and c-fos were performed to confirm their protein regulation and expression pattern in naïve and injured retinas. These analyses demonstrated that for these genes, protein regulation followed transcriptional regulation and that these pro-apoptotic proteins were expressed by retinal ganglion cells (RGCs). MetaCore-based death-signaling maps show that several apoptotic cascades were regulated in the retina following optic nerve injury and highlight the similarities and differences between IONT and IONC in cell death profiling.

Conclusions: This comprehensive time course retinal transcriptome study comparing IONT and IONC lesions provides a unique valuable tool to understand the molecular mechanisms underlying optic nerve injury and to design neuroprotective protocols.

Introduction

Optic nerve injury triggers retinal ganglion cell (RGC) death [1]. The progression of this death depends on the type of lesion, crush or transection, and on its distance from the eye [2-5].

The effects of intraorbital nerve transection (IONT) and intraorbital nerve crush (IONC) on retinal degeneration, and more specifically on RGC death, have been thoroughly studied by our group [5-8]. Both injuries trigger a massive RGC loss, which is slower and less acute after crush than after transection [3,4]. Thus, anatomically, it is observed that 38% and 20% of RGCs are lost seven days post-optic nerve transection or crush, respectively. Our IONT injury cleanly severs all the RGC axons, sparing the blood supply. The IONC injury compresses the blood vessels and axons in 10 s. Axons are thus severed by the acutely exerted pressure and by the glial reaction that occurs later in the retina and around the lesion site [1,9-13]. During the crush, the retinal blood supply is transiently hampered but in the present study we did check that our surgical manipulation avoided an ischemic insult.

Numerous efforts have been made to slow down RGC death triggered by optic nerve injury. So far none of them have successfully delayed RGC death beyond 15 days post-optic nerve injury [6-8,14-19]. Therefore, we purposed to investigate the molecular signals triggered by optic nerve injury in the retina. We have performed an exhaustive time-course analysis of the retinal transcriptome profile comparing mRNA expression from IONC and IONT retinas to naïve retinas using Affymetrix RAE230.2 arrays.

Sequences significantly regulated after IONT and after IONC at each time point were extracted and analyzed. Data were compared to learn the differences and similarities between both injuries. An extensive functional clustering was done, revealing that optic nerve injury alters highly important biologic functions in the retina as cell death, visual perception, and cytoskeleton. Due to the large number of regulated sequences and because RGC death is the most demolishing effect of optic nerve injury, we focused this work on the study of regulated genes involved in this process.

Methods

Animal handling and surgery

One hundred and forty adult female Sprague-Dawley rats (180–220 g bodyweight) were used for the array study and 48 for western and immunohistofluorescence studies. Rats were obtained from the university breeding colony. For anesthesia, we used xylazine (10 mg/kg bodyweight; Rompun; Bayer, Kiel, Germany) and ketamine (60 mg/kg bodyweight; Ketolar; Pfizer, Alcobendas, Madrid, Spain). All experimental procedures were performed in accordance with the Association for Research in Vision and Ophthalmology guidelines for the use of animals in research. Sterile precautions were maintained for all surgical procedures.

The array experiment-- Animals were divided into a control group (naïve), which did not undergo any experimental manipulation (n=20) and two experimental groups, one receiving an intraorbital nerve transection (IONT, n=60) and the other an intraorbital nerve crush (IONC, n=60). The left optic nerve was intraorbitally axotomized according to procedures that are standard in our laboratory [3-5,11]. Briefly, an incision was made in the superior orbital rim, the superoexternal orbital contents were dissected, and the superior and external rectus muscles were removed. When performing IONT, the dura mater was opened longitudinally to spare the blood supply, and the optic nerve was transected 0.5 mm from the optic disc [8,14]. The purpose of the IONC injury was to crush the entire population of optic nerve fibers. This lesion was performed by crushing the optic nerve 3 mm from the optic disc for 10 s with a pair of watchmaker’s forceps. Before and after the procedure, the eye fundus was observed through the operating microscope to assess the integrity of the retinal blood flow. In the very few instances in which blood supply did not restore within the following minutes, the animal was discarded.

Western blotting analyses-- Animals were subjected to surgery as above (IONT and IONC n=4 per time point and lesion).

Immunohistofluorescence analyses-- To identify RGCs in the sectioned retinas, they were retrogradely traced by applying Fluorogold (3% diluted in 10% DMSO-saline, Fluorochrome, LLC, Denver, CO) to the superior colliculus (SC) one week before the optic nerve injury was performed, as previously described [11-13]. Briefly the rats were anesthetized, each midbrain was exposed and, after removing the piamater overlying the SC, a small piece of gelatine sponge (Spongostan Film, Ferronsan, Denmark) soaked in 3% FG, diluted as aforementioned, was laid over the SC to label RGCs by retrograde axonal transport. Three animal groups were made: naïve, IONT, and IONC (n=5 per group). Optic nerve injury was performed as aforementioned, and animals were processed 48 h post lesion (hpl).

Array experimental procedures

Tissue processing and RNA extraction-- Animals were kept the appropriate time post-injury (12 h, 24 h, 48 h, 3 days, 7 days, and 15 days, n=8–12 per time point) and then sacrificed by an intraperitoneal overdose of sodium pentobarbital. Left retinas were freshly dissected and immediately frozen in liquid nitrogen. Four retinas from each time point, animal group, and biologic replica were pooled, and RNA was extracted using Trizol (Invitrogen, Barcelona, Spain). RNA was further cleaned through RNeasy mini kit columns (Qiagen, Izasa, Madrid, Spain). The RNA integrity was checked using a Bioanalyzer (Agilent Technologies, Santa Clara, CA), and concentrations were determined using a NanoDrop (NanoDrop Technologies, Wilmington, DE).

Array experiment-- RNA was labeled, arrays were hybridized, and a signal was acquired at the Genomic Service (Universidad Complutense de Madrid, Madrid, Spain) as described in [20]. Briefly, cDNA, which was created from 5 μg of each RNA sample, was synthesized with One-Cycle cDNA synthesis kit of Affymetrix (Santa Clara, CA) according to the Gene Chip Expression Analysis Technical Manual (2003; Affymetrix). From this cDNA, cRNA was synthesized according to the protocol of IVT Labeling kit of Affymetrix. Finally, this cRNA was purified with the Gene Chip Sample Cleanup Module of Affymetrix, recovering it in 22 μl of water. cRNA (15 μg) was fragmented to prepare the hybridization cocktails. Samples were hybridized to Rat genome 230 arrays (rat RAE230.2, Affymetrix), which are two gene chips that contain 31,000 probe-sets (genes or expressed sequence tags [ESTs]). Following hybridization, the gene chips were washed, stained, scanned, and analyzed with Gene Chip Operating Software (GCOS 1. 2; Affymetrix). The quality of the experiments was controlled by the visual inspection of the gene chips and by the presence of the spike controls and housekeeping control genes. A total of 35 gene chips were analyzed.

Each array replica was hybridized with an RNA pool of four retinas from independent RNA extractions and different animals (i.e., Array 1 was hybridized with naïve RNA from animals 1–4, array 2 with naïve RNA from animals 5–8 and so on). See Table 1 for a summary.

Array normalization and analysis-- Expression data in the CEL file format were imported directly into Bioconductor software for normalization and analysis [21]. Normalization was performed using the package justGCRMA (Robust Multi-Array expression measure using sequence information). To check the linearity and reproducibility that represent an adequate normalization of the arrays, scatter plots were performed between the log2 intensities of the injured retina arrays and the naïve retina arrays. Raw and normalized data for all hybridizations have been deposited as data series in NCBIs Gene Expression Omnibus (GEO, and are accessible through GEO Series accession number GSE9918).

The results from the IONT- and IONC-array analyses were compared to the results from the naïve-array ones. This allowed the identification of IONT- and IONC-regulated sequences. Statistically significant regulated sequences at each time point were extracted using the Limma package (Data analysis, linear models, and differential expression for micro array data; bioinf) with a cut-off for those sequences with a B value >0 and an adjusted p value false discovery rate (FDR) <0.01 and with a minimum regulation of 20% when compared to control. However, sequences modestly regulated along time would be missed with the Limma analysis since they would not be considered significant. To avoid this bias we also performed an analysis of the temporal regulation using the maSigPro package [21]. This package is specifically designed to analyze and compare time course experiments (p value FDR less than 0.01 and R2 threshold equal to 0.6). Log2 output values of regulated sequences were transformed to their net values (fold change). It is important to highlight here that when a sequence is not expressed, its mRNA level is 0. For this reason, the control level of expression was set to 1, any value less than 1 and more than 0 will then indicate downregulation and everything above 1 will indicate upregulation. IONT- and IONC-regulated sequences (p value <0,01, B value >0) were annotated and functionally clustered using the Database for Annotation, Visualization and Integrated Discovery, (DAVID [22]), the Kyoto Encyclopedia of Genes and Genomes database (KEGG), and PubMed reports. Clustering (molecular function and biologic process) was done using both Affymetrix RAE230.2 and the Rattus norvegicus genome as the background. Clusters that were significantly regulated were detected using the EASE tool [23] (p value FDR <1.00E-04). Finally, all regulated genes were imported to MetaCoreTM software (GeneGoLtd .) to generate “signaling maps” that include molecular interactions between the regulated sequences that can be associated to certain biologic processes such as cell death. Significantly regulated MetacoreTM cell-death signaling maps (pEASE <1.00E-04) together with their corresponding curated interactions were merged, and custom maps for IONT and IONC gene regulation were generated using the MapEditorTM software (GeneGo Ltd).

Quantitative real time polymerase chain reaction

One microgram of total RNA was used for cDNA synthesis using random hexamer primers (TaqMan Reverse Transcriptase, Applied Biosystems, Foster City, CA). Analysis of mRNA levels was performed using iQ SYBR Green Supermix (Bio-Rad Laboratories, Hercules, CA) in combination with a MyiQ Single Color Real-Time PCR Detection System (Bio-Rad Laboratories) and sequence specific primers (Table 2). Primers were designed using the software Primer Express 2.0 (Applied Biosystems) to target the same regions as the Affymetrix probes (predominantly in the 3′−UTR region). Melting temperature for the primers were approximately 59 °C and amplicon length approximately 80 bp. Primer concentrations were optimized for each primer pair and ranged between 150 and 500 nM final concentration. Polymerase chain reactions (PCRs) were performed in triplicate in a 96 well plate according to the manufacturer’s protocol. For each series, a threshold cycle value (Ct) corresponding to the PCR cycle at which the fluorescent emission reaches a threshold above the baseline emission was determined. Based on this, initial mRNA levels were calculated. To confirm specificity of amplified PCR products, dissociation curve analyses and agarose gel electrophoresis were performed. As control genes RPL30, RPL31, RPL37, β-actin, and hypoxanthine guanine phosphoribosyl transferase (HPRT) were measured. These control genes were chosen because no change among samples was observed in the array experiment. Initial mRNA levels were normalized to the average mRNA levels of RPL30, RPL31, and RPL37 to correct for variations in the cDNA synthesis. Relative mRNA levels in naïve retinas were used as a reference.

Immunohistofluorescence

Naïve, IONT-, and IONC-injured animals whose retinal ganglion cells had been previously traced with Fluorogold (n=5 per group, 48 h post-axotomy) were deeply anesthetized with pentobarbital and perfused transcardially with 4% paraformaldehyde (PFA in 0.1 M phosphate buffer) after a saline rinse. Retinas were dissected and post-fixed for 2 h at 4 °C before switching to 30% sucrose (Sigma, Alcobendas, Madrid, Spain) for at least 12 h. Tissue was embedded in optimal cutting temperature (OCT) compound (Sakura Finetek, Torrance, CA), frozen on liquid nitrogen, and kept at −70 °C. Retinas were sectioned at a thickness of 15 μm on a cryostat. Sections were blocked in 2% donkey serum in phosphate buffered saline (PBS) with 0.1% Triton (PBST) and incubated overnight a 4 °C with the appropriate antibody. Immunoreactivity was detected using fluorescence-conjugated secondary antibodies (4 h at room temperature in PBST). Images were taken with a CCD camera using the ImageProPlus software (Image-Pro® Plus 5.1 for Windows®; Media Cybernetics, Inc., Silver Spring, MD) and further processed with Adobe Photoshop 7.0 (Adobe Systems Inc., San Jose, CA).

Western blot analysis

Freshly dissected retinas from naïve, IONT-, and IONC-injured animals (12 h, 48 h, and 3 days or 7 days post-axotomy, n=4 per group and time point) were immediately frozen in liquid nitrogen and stored at −70 °C until protein was extracted. Tissue was homogenized in 300 μl of lysis buffer (1% Nonidet-p40, 20mM Hepes pH 7.4, 100 mM NaCl, 100 mM NaF, 1 mM Na3VO4, and 5 mM EDTA with 1X protease inhibitor cocktail, Roche Diagnosis, Barcelona, Spain) and then incubated on ice for 2 h. After centrifugation at 18928x g for 15 min at 4 °C, the supernatant was collected to establish total protein values. The amount of protein was determined using a bicinchoninic acid (BCA) protein assay kit (Pierce, ThermoFisher Scientific, Cultek SL, Madrid, Spain) and β-actin detection. β-actin western blots were scanned, and intensities were measured (Gene Tool, Syngene, Synoptics Ltd, Cambridge, UK). If necessary, the loading volume was corrected to equal the amount of loaded protein among all samples. Protein (20–60 μg) was loaded on 16% or 12% SDS-polyacrylamide gels for standard electrophoresis. Proteins were transferred to a polyvinylidene fluoride (PVDF; Amersham, GE Healthcare Europe GmBH, Barcelona, Spain) membrane (Amersham Pharmacia Biotech). Membranes were blocked with 5% fat free milk in PBS with 0.1% Tween 20 and incubated with the appropriate antibody overnight at 4 °C. Membranes were then incubated with horseradish peroxidase (HRP)-conjugated secondary antibodies. HRP activity was visualized by applying chemiluminescent substrate (Enhanced Chemo Luminiscence (ECL), Amersham, GE Healthcare Europe GmBH, Barcelona, Spain) followed by exposure of the membrane to an X-ray film. The exposed films were analyzed with the Gen Tools program (Syngene). The signal intensity of treated retinas was referred to the signal in naïve ones, which were arbitrarily considered as 100%. To avoid biologic variability, extracts from four animals were loaded in parallel, and to avoid technical variability, each western analysis was performed three times. Data shown are the averaged values of these replicas with their standard error of the mean (SEM) value.

Antibodies and working dilutions

All antibodies were diluted in phosphate buffered saline with 0.1% Triton (PBSTx) or phosphate buffered saline with 0.1% Tween 20 (PBSTw) for immunohistofluorescence and western blotting experiments, respectively.

Primary antibodies-- Rabbit anti-β-actin (Santa Cruz Biotechnologies, Heidelberg, Germany) was diluted at 1:200 for western blotting analysis. Rabbit anti-TRADD (Abcam, Cambridge, UK) was diluted 1:70 for western blotting and 1:50 for immnunofluorescence experiments. Mouse anti-TNF- R1 (H-5, Santa Cruz Biotechnologies, Heidelberg, Germany) was diluted: 1:200 for western blotting and 1:50 for immunohistofluorescence studies. Rabbit anti-c-Fos (4, Santa Cruz Biotechnologies, Heidelberg, Germany) was diluted 1:200 for western blotting and 1:50 for immunohistofluorescence experiments. Rabbit anti-active Caspase 3 (H-277, Santa Cruz Biotechnologies, Heidelberg, Germany) was diluted 1:200 for western blotting analysis. Rabbit anti-total Caspase 3: (BD Biosciences PharMingen, Franklin Lakes, NJ) was diluted 1:200 for western blotting and immunohistofluorescence experiments.

Secondary antibodies-- Peroxidase-conjugated AffiniPure Goat anti-Mouse IgG (H+L) and peroxidase- conjugated AffiniPure Goat anti-Rabbit IgG (H+L; both from Jackson ImmunoResearch, Suffolk, UK) were diluted 1:5000 in PBSTw for secondary detection in western blot analyses. In immunohistofluorescence experiments secondary detection was performed using Alexa Fluor 568-goat anti-mouse IgG (H+L; Molecular Probes, Invitrogen, Barcelona, Spain) and Cy3-goat anti-Rabbit IgG (H+L; Jackson ImmunoResearch) diluted 1:500 in PBSTx. These antibodies were previously checked for unspecific labeling (4 h incubation onto retinal sections without primary antibody incubation).

Results

General overview: number of sequences that are regulated in the retina by optic nerve injury

We performed a time-course analysis of the retinal transcriptome in two different models of pathological RGC death, optic nerve transection and optic nerve crush, and compared them to naïve retinas. Based on anatomic data from our group [3,5], we choose five time points post-lesion (pl), 12 h, 24 h, 48 h, 3 day, 7 day, and 15 day (further details on number and array replicas see Methods). At early time points (from 12 hpl to 3 days post lesion [dpl]) no RGC death is observed. At 7 dpl, 38% (IONT) and 20% (IONC) of the RGCs have died. Finally, 90% and 68% of the RGCs have disappeared two weeks after IONT and IONC, respectively.

The global number of sequences (genes and ESTs) that were regulated in this study (p value FDR <0,01, B value >0) is shown in Figure 1A. The first observation is that IONT triggers a higher retinal transcriptome regulation (by 36.2%) than IONC does (3,219 sequences regulated after IONT and 1,996 after IONC). This is not due to the extra time point analyzed after IONT (day 15) since if these sequences are removed, IONT regulation is still higher by 29.8% (2,846 IONT sequences regulated from 12 hpl until 7 dpl). Out of all the regulated sequences, 2,071 sequences correspond to genes (Figure 1B), which can be divided into: a) genes regulated only by IONT (1,117, IONT-specific genes); b) genes regulated only by IONC (293, IONC-specific genes); and c) genes regulated after both injuries (661, IONT and IONC common genes).

Number of sequences regulated after intraorbital nerve transection and intraorbital nerve crush in the adult rat retina and their time course expression

Genes-- The number of upregulated and downregulated genes and their temporal course of expression are summarized in Figure 1C. Gene regulation is mainly transient after IONC (Figure 1C, right) where 59.6% out of the total IONC-triggered genes are regulated only at one time point, namely at 12 hpl. This differs from IONT where most genes (68.7%) are regulated at more than one time point (Figure 1C, left). Thus, IONT induces mainly a sustained response and IONC a transitory response. In Figure 1C, it is also noticeable that after IONT, 60% of the genes are upregulated whereas after IONC 67% are down-regulated. Genes that are regulated after both lesions follow the same pattern of upregulation or downregulation with just a few exceptions (Appendix 1). However, though similarly regulated after both injuries, they do not always follow the same time course. This is observed in Figure 1C where the gray bars (common genes) in the IONT graph do not correlate to the same time points shown in the IONC graph and vice versa.

Expressed sequence tags-- The pattern of temporal distribution and regulation of ESTs resembles that of the annotated sequences as shown in Figure 1D. Since the optic nerve injury triggers a high number of gene regulation responses in the retina, we focused this study on the regulation of annotated sequences, leaving aside the ESTs.

Functional clustering of optic nerve injury-regulated genes

Genes regulated after IONC and IONT were clustered using the web-based tools, DAVID and EASE [24,25]. The main clusters (pEASE <1,00E-04) and cell pathways regulated after IONT and IONC and the number of genes in each one are summarized in Figure 2. A detailed and comprehensive list can be found in Appendix 2 and Appendix 3. Each main cluster was hierarchically divided into sub-clusters based on more specific functions (i.e., main cluster included cell death and sub-clusters included apoptosis, autophagy, and inflammatory response linked to cell death). The highest regulated biologic processes were signaling (cluster 1, pathways 1–9), cytoskeleton and associated processes (cluster 13, pathways 12–16), primary metabolism (clusters 2 and 3 and pathways 22–25 and 27–33 where generation of cellular energy, pathways 30–33, is altered mainly after IONT, although glycolysis and gluconeogenesis are affected by IONC as well), protein metabolism (cluster 3, pathways 21 and 23), immune response and inflammation (cluster 7, pathways 10 and 11, which are mainly upregulated after both injuries), ribosomal protein genes (pathway 34, regulated after IONT but not after IONC), and RNA metabolism, processing, and translation (cluster 8, affected after both lesions but mainly upregulated after IONT and downregulated after IONC).

As expected, cell death was highly regulated (cluster 5 pEASE value 5.60E-08 for IONT and 2.20E-7 for IONC). Our data indicate that the retinal cell death caused by optic nerve injury is apoptotic (pathway 4, pEASE value 2,40E-19 for IONT and 4,80E-06 for IONC) since no other cell death mechanism was found statistically regulated. In total, 186 genes were clustered under cell death out of which 89 and 29 were specific for IONT or IONC, respectively, and 68 were regulated by both lesions (Appendix 3). Some of these genes have been further analyzed.

Array validation

Once genes were annotated and clustered, 10 IONT and IONC upregulated genes related with death and inflammation were chosen for array validation using quantitative-reverse transcriptase-polymerase chain reaction (qRT–PCR). The results are shown in Table 3. The analyzed genes follow the same trend after both experiments. However, differences in net regulation between the results that were obtained after both techniques were observed. These differences can be explained by the variations in experimental set-up and the sensitivity of each technique. In the array experiments, all transcripts hybridize simultaneously and the hybridization is not optimized for each single sequence. In the qRT–PCR, each reaction is optimized for measuring the single cDNA of a transcript that the primers were designed to bind to. For the majority of genes studied using qRT–PCR, their regulation was very similar to that observed in the array experiment. The only exception is Eef2K that is regulated after IONT when analyzed by qRT–PCR but not in the array analysis. This was probably due to the low level of upregulation observed in the array experiment 12 h after IONT (1.14 fold change compared to naïve levels). The low level of upregulation together with the variability among the array replicas (SEM: 0.15) resulted in a not significant change (p value FDR=0.07). In conclusion, there is a correlation between both experimental techniques, which validates the array experiment.

Temporal expression level and retinal pattern distribution of a few death-related proteins

The mRNA of tumor necrosis factor receptor 1a (TNFR1a), Caspase 3 (Casp3) and c-fos is upregulated after both optic nerve injuries (Appendix 3). We used western blotting to study whether their protein levels followed their mRNA regulation. In addition, we measured active Casp3 to know whether its upregulation was followed by its activation. Casp3 implication in IONT-induced RGC death has been described previously [24,26,27] but to our knowledge not in the IONC model. In addition, the temporal expression analysis of this executioner protein will provide insight as when the RGCs are committed to death in our models. TNFR1 has been implicated in RGC death in the mouse retina following optic nerve crush [25]. TNFR1a, depending on its downstream adaptor, could be linked to apoptosis or survival [28]. For this reason, we also measured the regulation of tumor necrosis factor associated death domain (TRADD), which links TNFR1a and TNFR12A (Fn14) signaling to apoptosis [29,30]. Clarify here that the arrays did not contain probes for TRADD mRNA and therefore it was not possible to asses for its mRNA regulation. The results are shown in Figure 3, TNFR1a, Casp3, and c-fos proteins are upregulated along time after both injuries. This is in agreement with the array results. c-fos upregulation after both injuries contrasts with other reports [31,32] where this protein does not show regulation. This may be due to the different time points analyzed together with the different sensitivity of the techniques used. The levels of TRADD protein also increased after both injuries in our study. As for the active form of Casp3, it peaked at 48 h after IONT and IONC, which was earlier than previously described [24]. This temporal difference may be a reflection of the distance from the eye where the injury is performed, which is closer in our model. This supports previous reports, which demonstrate that the closer to the eye the injury occurs, the more rapid the RGC degeneration is [5].

To find out whether these upregulated pro-apoptotic proteins were expressed by RGCs, we performed immunohistofluorescence on radial sections from naïve, IONT-, and IONC-injured retinas. To identify RGCs, these were traced with Fluorogold (FG) one week before the injury was performed. Animals were kept 48 h post-optic nerve injury since it was at this time point when the protein levels peaked as shown by the western blot analyses. Results are shown in Figure 4. Naïve and injured RGCs expressed TRADD and c-fos as seen by colocalization of the fluorogold tracer and red fluorescence immunoreactivity. Casp3 and TNFR1a proteins were detected in injured RGCs but not in naïve RGCs.

Discussion

We have compared the temporal profiling of the retinal transcriptome after two types of optic nerve injury that trigger RGC death, intraorbital nerve transection and intraorbital nerve crush [1-3,5]. The results show that optic nerve injury triggers a plethora of gene regulation responses in the retina as soon as 12 hpl. Injury-regulated sequences fall in three categories, those uniquely regulated after IONT, those uniquely regulated after IONC, and those regulated in a similar way after either of both injuries. There are two different temporal responses, genes that are transiently regulated and genes that have a sustained regulation over time. Three main differences are observed between both injuries, 1) the number of the regulated sequences, 2) the temporal profiles of the regulated sequences, and 3) the majority of upregulation or downregulation. The number of regulated sequences is higher after IONT compared to the number of sequences after IONC. While IONC mainly triggers transient responses that are particularly evident at 12 hpl, IONT mainly triggers a sustained response. Finally, after IONC, the overall response is downregulation whereas after IONT, the main response is upregulation. These differences are likely to reflect the progression of retinal cell death and may account for the different time-courses of RGC death observed after IONT and IONC [2-5].

Fifteen days post-IONT, hundreds of genes are regulated in agreement with previous reports [33]. At this time post-lesion, just a small proportion of RGCs is still alive. Therefore, these changes might be attributable to other retinal cells. However, this does not mean that these cells are dying but rather that they are responding to the lesion inflicted to the RGCs.

Functional clustering of regulated genes revealed that many biologic processes and signaling pathways are involved in the response to optic nerve injury. In agreement with previous studies [34,35] and according to our results, the main cell death mechanism is apoptosis. IONT regulates a higher number of death-related genes than IONC. This concurs with the more rapid course of RGC death observed after IONT than after IONC. We have focused this study on genes directly involved in cell death processes since RGC death is the most pronounced effect after optic nerve injury. There is a possibility that RGC death is not a primary effect of the injuries but rather a secondary one. Thus, it could not be ruled out that the optic nerve injury may primarily affect other processes in the retina such as cell-cell communication or even basic metabolic alterations that will eventually trigger cell death.

Since both injuries lead to RGC death, it is likely that the commonly regulated cell death-related genes are responsible for retinal degeneration associated to axonal injury. Therefore, we chose ten of those genes to validate our array data by qRT–PCR. This is the first time that the retinal regulation of these genes by both models of optic nerve injury is reported with the exception of clusterin, which has been shown to be upregulated at 15 dpl after IONT [33]. Interestingly, clusterin, lipocalin 2, litaf, and TNFR1a are also upregulated by other injuries as elevation of intraocular pressure (IOP) or ischemia [36-40]. An increase of Clusterin and Lipocalin 2 expression together with that of ceruloplasmin and early growth response 1 (Egr1) is associated with increased apoptosis and immune-inflammatory responses [41,42]. All of those genes, which were upregulated after IONT and IONC, are also upregulated after ischemia, retinal scraping, experimental glaucoma, and phototoxic lesions (in the case of Egr1). These injuries trigger degeneration of retinal neurons, and so the similar upregulation of these genes after different retinal insults indicates that there might be a common response to injury. Whether such common responses are involved in or may modulate the progression of injury-associated cell death is subject for further studies.

To bring the validation of the arrays one step further and to gain some additional information about the translatability of the transcriptional changes to the protein level, we performed western blot and immunohistofluorescence analyses in naïve, IONT-, and IONC-injured retinas. These results demonstrate that for these genes, protein regulation follows transcriptional regulation. In addition, we showed that these pro-apoptotic proteins are expressed by RGCs, and this result gives support that the identified networks are active in injured RGCs. In conclusion, the expression pattern and regulation of TNFR1a, TRADD, and the executioner protease, Caspase 3, indicate that the extrinsic pathway to apoptosis [43] is active in RGCs undergoing IONT or IONC.

Signaling maps

IONT- and IONC-regulated genes were loaded into MetaCoreTM to search for regulated gene networks and signaling maps. Death signaling maps that were significantly regulated by optic nerve injury (pEASE <1.00E-04) were selected and merged for further analysis and visualization purposes using MapEditorTM. The merged map reflects the main gene networks regulated by optic nerve injury that lead to cell death. (Figure 5 and Appendix 4). In addition to the array results, regulation after IONT and IONC for many of the genes shown in the map has been validated by both qRT–PCR and western blot. Besides, some of these proteins have been shown to be expressed by injured RGCs. These maps visualize the retinal responses to IONC and IONT and are based on a significant over-representation of regulated genes. However, we are aware that these are partial maps of the retinal response to optic nerve injury since just 47 out of the 189 death-related genes which are regulated have been included. Furthermore, the data that support the maps belong to “known gene interactions.” Therefore, these maps are not irrefutable and should be considered as a step to unravel the molecular mechanisms underlying optic nerve injury- induced retinal degeneration. As shown in Figure 5 and in detail in the Appendix 4, eight main cellular processes, which may lead to cell death, were identified as being involved in the retinal response to optic nerve injury: stress response, induction of apoptosis by extracellular signals (extrinsic pathway), by intracellular signals (pro-apoptotic genes, intrinsic pathway), regulation of survival signals, activation of immune-inflammatory pathways, DNA damage, cell cycle deregulation, and the lysosomal death pathway. Interestingly, both IONT and IONC trigger TNFR1-TRADD signaling, which leads to the activation of the extrinsic apoptosis pathway [29,30,43,44]. This is strengthened by the upregulation after IONT and IONC of caspases 3 and 11 (caspase 4 in humans) and of another member of the tumor necrosis factor family, TNRSF12A (Fn14), which is also linked to apoptosis through TRADD signaling [29,30,45]. This is further supported by the upregulation of caspase 8 in axotomized RGCs reported by Weishaupt et al., 2003. [46]. Pro-apoptotic genes, associated to the intrinsic apoptotic pathway, are regulated after IONT and IONC (14–3–3, RAD9, RASSF4, Rb protein, Smad4, PPM1B), although IONT (BMF, Bcl.10, Bax, Bcl-B, DLC1, MCL1) seems to trigger a higher response than IONC does. Some of these genes are linked to cell cycle control as 14–3–3, RAD9, and PPM1B. Cell cycle deregulation has been reported as a threshold for cell death [47], and genes involved in the cell cycle control are regulated after both injuries (Cyclin G1, GADD45α, Rb protein, and aforementioned). Again, there is a higher response after IONT (Cyclins, CDK1, Aurora, p27KIP1, p21) compared to IONC. Interestingly, an upregulation of cell cycle-related genes is observed as well after spinal cord lesion, another model of central nervous system axonal injury and neuronal degeneration [48]. DNA damage is another trigger of apoptosis [49] and GADD45α and nibrin, involved in this process, are regulated after optic nerve injury (the former after both lesions and the latter only after IONT). Increased levels of caspase 2 are important in cell death activated by DNA damage [50], and we found this gene to be downregulated by both injuries, which suggests that DNA damage is not the main pathway for the optic nerve injury-induced retinal degeneration. Survival and anti-apoptotic signals [51] such as XIAP, ERK2, SOS, c-Raf-1, and YY1 are downregulated. It is worth highlighting here that among the genes linked to apoptosis triggered by mitochondrial proteins (Bad and Bid among others) only Bax is regulated. Interestingly, it is regulated only by transection and quite late after the injury at 3 dpl. This would agree with previous reports where administration of Bcl-2 simultaneously to IONT did not delay RGC-associated death [6]. Stress response to injury is evidenced by the regulation of calcineurin and RPS6K after IONT and C/EBP, JNK, and heat shock protein [52,53] genes after IONT and IONC. There is a strong immune-inflammatory response [54] as seen by the upregulation of JAK2, Litaf, Stat 1, and caspase 11 among others. It is important to highlight here that some of these genes such as caspase 11, stat 1, and litaf are also related to the apoptotic pathway [55-57]. Finally, optic nerve injury regulates several cathepsins, (H, L, and S after IONT and C and Z after IONT and IONC). Cathepsins are proteases located in the lysosomes, which have a role in apoptosis through the “lysosomal pathway of apoptosis” [57-59]. Lysosomal permeabilization and release of cathepsins are facilitated by sphingosine, a lipid whose overexpression is linked to apoptosis [60]. Our array study reveals that the RNA of some enzymes linked to sphingosine metabolism are soon upregulated after IONT (Sphk1 and Sgpl1) and after IONC (Sphk1 and 2) [58-60]. Cathepsins’ role in apoptosis has two sides. First, they act as proteases cleaving and thus activating caspases. Second, they act on the mitochondria to induce mitochondrial dysfunction [58,61].

Together, all this indicates that IONT triggers a more acute and more pronounced apoptotic response than IONC [62]. First, more genes are regulated after IONT, and those, which are commonly regulated, show in general a longer lasting regulation after transection than after crush. Second, the regulation of the extrinsic apoptotic pathway and immune-inflammatory response is similar after IONT and IONC. However, lysosomal cell death, DNA damage, cell cycle control, and pro-apoptotic signals that may activate the intrinsic apoptotic pathway are mostly deregulated after IONT. All these may account for the slower degeneration of RGC after optic nerve crush than after transection and, most importantly, it highlights the common responses that may cause the RGC degeneration triggered by optic nerve injury.

Appendix 1. Genes commonly regulated by intraorbital nerve crush and intraorbital nerve transection that follow opposite regulation

Appendix 2. Functional clusters regulated in the retina after intraorbital nerve transection and intraorbital nerve crush.

Appendix 3. Time course gene regulation in adult rat retinas after intraorbital nerve transection and intraorbital nerve crush and functional clustering

Appendix 4. Cell-death signaling maps based on retinal gene regulation after intraorbital nerve crush and intraorbital nerve transection.

Acknowledgments

We thank A. Harrison for his tips on array analysis. This work is funded by ISCII: CP003/00119 and BIOCARM 2005/01–64692

Reference

  1. Kermer P, Klocker N, Weishaupt JH, Bahr M. Transection of the optic nerve in rats: studying neuronal death and survival in vivo. Brain Res Brain Res Protoc. 2001; 7:255-60. [PMID: 11431127]
  2. Berkelaar M, Clarke DB, Wang YC, Bray GM, Aguayo AJ. Axotomy results in delayed death and apoptosis of retinal ganglion cells in adult rats. J Neurosci. 1994; 14:4368-74. [PMID: 8027784]
  3. Parrilla-Reverter G. Sobrado, P., Mayor, S., Canovas, I., Aguilera, M.E., Lambrou, G.N., Latour. E., Villegas-Pérez, M.P., Vidal-Sanz, M (2004) Intraorbital optic nerve crush induces progressive retinal ganglion cell loss ARVO Annual Meeting; 2004 April 25–29; Fort Lauderdale (FL).
  4. Parrilla-Reverter G. Sobrado, P., Mayor, S., Canovas, I., Bernal, J.M. Soro, M.I., Aguilera, M.E. Villegas-Pérez, M.P., Vidal-Sanz, M (2006) Retinal ganglion cell axotzomy induced by intraorbital nerve crush or optic nerve transection results in different time course degeneration and expression of neurofilaments ARVO Annual Meeting; 2006 April 30-May 4; Fort Lauderdale (FL).
  5. Villegas-Perez MP, Vidal-Sanz M, Rasminsky M, Bray GM, Aguayo AJ. Rapid and protracted phases of retinal ganglion cell loss follow axotomy in the optic nerve of adult rats. J Neurobiol. 1993; 24:23-36. [PMID: 8419522]
  6. Vidal-Sanz M, Lafuente M, Sobrado-Calvo P, Selles-Navarro I, Rodriguez E, Mayor-Torroglosa S, Villegas-Perez MP. Death and neuroprotection of retinal ganglion cells after different types of injury. Neurotox Res. 2000; 2:215-27. [PMID: 16787842]
  7. Peinado-Ramon P, Salvador M, Villegas-Perez MP, Vidal-Sanz M. Effects of axotomy and intraocular administration of NT-4, NT-3, and brain-derived neurotrophic factor on the survival of adult rat retinal ganglion cells. A quantitative in vivo study. Invest Ophthalmol Vis Sci. 1996; 37:489-500. [PMID: 8595949]
  8. Villegas-Perez MP, Vidal-Sanz M, Bray GM, Aguayo AJ. Influences of peripheral nerve grafts on the survival and regrowth of axotomized retinal ganglion cells in adult rats. J Neurosci. 1988; 8:265-80. [PMID: 2448429]
  9. Panagis L, Thanos S, Fischer D, Dermon CR. Unilateral optic nerve crush induces bilateral retinal glial cell proliferation. Eur J Neurosci. 2005; 21:2305-9. [PMID: 15869529]
  10. Salvador-Silva M, Vidal-Sanz M, Villegas-Perez MP. Microglial cells in the retina of Carassius auratus: effects of optic nerve crush. J Comp Neurol. 2000; 417:431-47. [PMID: 10701865]
  11. Sobrado-Calvo P, Vidal-Sanz M, Villegas-Perez MP. Rat retinal microglial cells under normal conditions, after optic nerve section, and after optic nerve section and intravitreal injection of trophic factors or macrophage inhibitory factor. J Comp Neurol. 2007; 501:866-78. [PMID: 17311318]
  12. Lafuente MP, Villegas-Pérez MP, Sellés-Navarro I, Mayor-Torroglosa S, Miralles de Imperial J, Vidal-Sanz M. Retinal ganglion cell death after acute retinal ischemia is an ongoing process whose severity and duration depends on the duration of the insult. Neuroscience. 2002; 109:157-68. [PMID: 11784707]
  13. Selles-Navarro I, Villegas-Perez MP, Salvador-Silva M, Ruiz-Gomez JM, Vidal-Sanz M. Retinal ganglion cell death after different transient periods of pressure-induced ischemia and survival intervals. A quantitative in vivo study. Invest Ophthalmol Vis Sci. 1996; 37:2002-14. [PMID: 8814140]
  14. Vidal-Sanz M, Bray GM, Villegas-Perez MP, Thanos S, Aguayo AJ. Axonal regeneration and synapse formation in the superior colliculus by retinal ganglion cells in the adult rat. J Neurosci. 1987; 7:2894-909. [PMID: 3625278]
  15. Cheng L, Sapieha P, Kittlerova P, Hauswirth WW, Di Polo A. TrkB gene transfer protects retinal ganglion cells from axotomy-induced death in vivo. J Neurosci. 2002; 22:3977-86. [PMID: 12019317]
  16. Cui Q, So KF, Yip HK. Major biological effects of neurotrophic factors on retinal ganglion cells in mammals. Biol Signals Recept. 1998; 7:220-6. [PMID: 9730581]
  17. Cui Q, Cho KS, So KF, Yip HK. Synergistic effect of Nogo-neutralizing antibody IN-1 and ciliary neurotrophic factor on axonal regeneration in adult rodent visual systems. J Neurotrauma. 2004; 21:617-25. [PMID: 15165369]
  18. Mansour-Robaey S, Clarke DB, Wang YC, Bray GM, Aguayo AJ. Effects of ocular injury and administration of brain-derived neurotrophic factor on survival and regrowth of axotomized retinal ganglion cells. Proc Natl Acad Sci USA. 1994; 91:1632-6. [PMID: 8127857]
  19. Zhang CW, Lu Q, You SW, Zhi Y, Yip HK, Wu W, So KF, Cui Q. CNTF and BDNF have similar effects on retinal ganglion cell survival but differential effects on nitric oxide synthase expression soon after optic nerve injury. Invest Ophthalmol Vis Sci. 2005; 46:1497-503. [PMID: 15790921]
  20. Arbillaga L, Azqueta A, van Delft JH, López de Cerain A. In vitro gene expression data supporting a DNA non-reactive genotoxic mechanism for ochratoxin A. Toxicol Appl Pharmacol. 2007; 220:216-24. [PMID: 17316727]
  21. Conesa A, Nueda MJ, Ferrer A, Talon M. maSigPro: a method to identify significantly differential expression profiles in time-course microarray experiments. Bioinformatics. 2006; 22:1096-102. [PMID: 16481333]
  22. Dennis G, , Jr Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003; 4(5)P:3. 12734009
  23. Hosack DA, Dennis G, , Jr Sherman BT, Lane HC, Lempicki RA. Identifying biological themes within lists of genes with EASE. Genome Biol. 2003; 4:R70 [PMID: 14519205]
  24. Cheung ZH, Chan YM, Siu FK, Yip HK, Wu W, Leung MC, So KF. Regulation of caspase activation in axotomized retinal ganglion cells. Mol Cell Neurosci. 2004; 25:383-93. [PMID: 15033167]
  25. Tezel G, Yang X, Yang J, Wax MB. Role of tumor necrosis factor receptor-1 in the death of retinal ganglion cells following optic nerve crush injury in mice. Brain Res. 2004; 996:202-12. [PMID: 14697498]
  26. Kermer P, Klocker N, Labes M, Thomsen S, Srinivasan A, Bahr M. Activation of caspase-3 in axotomized rat retinal ganglion cells in vivo. FEBS Lett. 1999; 453:361-4. [PMID: 10405176]
  27. Kermer P, Klocker N, Labes M, Bahr M. Insulin-like growth factor-I protects axotomized rat retinal ganglion cells from secondary death via PI3-K-dependent Akt phosphorylation and inhibition of caspase-3 In vivo. J Neurosci. 2000; 20:2-8. [PMID: 10632601]
  28. Dempsey PW, Doyle SE, He JQ, Cheng G. The signaling adaptors and pathways activated by TNF superfamily. Cytokine Growth Factor Rev. 2003; 14:193-209. [PMID: 12787559]
  29. Ashkenazi A, Dixit VM. Death receptors: signaling and modulation. Science. 1998; 281:1305-8. [PMID: 9721089]
  30. Ashkenazi A, Dixit VM. Apoptosis control by death and decoy receptors. Curr Opin Cell Biol. 1999; 11:255-60. [PMID: 10209153]
  31. Robinson GA. Immediate early gene expression in axotomized and regenerating retinal ganglion cells of the adult rat. Brain Res Mol Brain Res. 1994; 24:43-54. [PMID: 7968376]
  32. Yang Z, Quigley HA, Pease ME, Yang Y, Qian J, Valenta D, Zack DJ. Changes in gene expression in experimental glaucoma and optic nerve transection: the equilibrium between protective and detrimental mechanisms. Invest Ophthalmol Vis Sci. 2007; 48:5539-48. [PMID: 18055803]
  33. Piri N, Kwong JM, Song M, Elashoff D, Caprioli J. Gene expression changes in the retina following optic nerve transection. Mol Vis. 2006; 12:1660-73. [PMID: 17200666]
  34. Reme CE, Grimm C, Hafezi F, Marti A, Wenzel A. Apoptotic cell death in retinal degenerations. Prog Retin Eye Res. 1998; 17:443-64. [PMID: 9777646]
  35. Weishaupt JH, Bahr M. Degeneration of axotomized retinal ganglion cells as a model for neuronal apoptosis in the central nervous system - molecular death and survival pathways. Restor Neurol Neurosci. 2001; 19:19-27. [PMID: 12082226]
  36. Ahmed F, Brown KM, Stephan DA, Morrison JC, Johnson EC, Tomarev SI. Microarray analysis of changes in mRNA levels in the rat retina after experimental elevation of intraocular pressure. Invest Ophthalmol Vis Sci. 2004; 45:1247-58. [PMID: 15037594]
  37. Fontaine V, Mohand-Said S, Hanoteau N, Fuchs C, Pfizenmaier K, Eisel U. Neurodegenerative and neuroprotective effects of tumor Necrosis factor (TNF) in retinal ischemia: opposite roles of TNF receptor 1 and TNF receptor 2. J Neurosci. 2002; 22:RC216 [PMID: 11917000]
  38. Gwon JS, Kim IB, Lee MY, Oh SJ, Chun MH. Expression of clusterin in Muller cells of the rat retina after pressure-induced ischemia. Glia. 2004; 47:35-45. [PMID: 15139011]
  39. Tezel G, Li LY, Patil RV, Wax MB. TNF-alpha and TNF-alpha receptor-1 in the retina of normal and glaucomatous eyes. Invest Ophthalmol Vis Sci. 2001; 42:1787-94. [PMID: 11431443]
  40. Yoshimura N, Kikuchi T, Kuroiwa S, Gaun S. Differential temporal and spatial expression of immediate early genes in retinal neurons after ischemia-reperfusion injury. Invest Ophthalmol Vis Sci. 2003; 44:2211-20. [PMID: 12714663]
  41. Krones-Herzig A, Mittal S, Yule K, Liang H, English C, Urcis R, Soni T, Adamson ED, Mercola D. Early growth response 1 acts as a tumor suppressor in vivo and in vitro via regulation of p53. Cancer Res. 2005; 65:5133-43. [PMID: 15958557]
  42. Lee S, Lee J, Kim S, Park JY, Lee WH, Mori K, Kim SH, Kim IK, Suk K. A dual role of lipocalin 2 in the apoptosis and deramification of activated microglia. J Immunol. 2007; 179:3231-41. [PMID: 17709539]
  43. Jin Z, El Deiry WS. Overview of cell death signaling pathways. Cancer Biol Ther. 2005; 4:139-63. [PMID: 15725726]
  44. Khosravi-Far R, Esposti MD. Death receptor signals to mitochondria. Cancer Biol Ther. 2004; 3:1051-7. [PMID: 15640619]
  45. Yepes M. Tweak and FN14 in central nervous system health and disease. Front Biosci. 2007; 12:2772-81. [PMID: 17485258]
  46. Weishaupt JH, Diem R, Kermer P, Krajewski S, Reed JC, Bahr M. Contribution of caspase-8 to apoptosis of axotomized rat retinal ganglion cells in vivo. Neurobiol Dis. 2003; 13:124-35. [PMID: 12828936]
  47. Greene LA, Liu DX, Troy CM, Biswas SC. Cell cycle molecules define a pathway required for neuron death in development and disease. Biochim Biophys Acta. 2007; 1772:392-401. [PMID: 17229557]
  48. 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-68. [PMID: 12666113]
  49. Norbury CJ, Zhivotovsky B. DNA damage-induced apoptosis. Oncogene. 2004; 23:2797-808. [PMID: 15077143]
  50. Zhivotovsky B, Orrenius S. Caspase-2 function in response to DNA damage. Biochem Biophys Res Commun. 2005; 331:859-67. [PMID: 15865942]
  51. Frebel K, Wiese S. Signalling molecules essential for neuronal survival and differentiation. Biochem Soc Trans. 2006; 34:1287-90. [PMID: 17073803]
  52. Beere HM. “The stress of dying”: the role of heat shock proteins in the regulation of apoptosis. J Cell Sci. 2004; 117:2641-51. [PMID: 15169835]
  53. Beere HM. Death versus survival: functional interaction between the apoptotic and stress-inducible heat shock protein pathways. J Clin Invest. 2005; 115:2633-9. [PMID: 16200196]
  54. Martinon F, Tschopp J. Inflammatory caspases and inflammasomes: master switches of inflammation. Cell Death Differ. 2007; 14:10-22. [PMID: 16977329]
  55. Bromberg J, Darnell JE, Jr. The role of STATs in transcriptional control and their impact on cellular function. Oncogene. 2000; 19:2468-73. [PMID: 10851045]
  56. Kim HS, Lee MS. STAT1 as a key modulator of cell death. Cell Signal. 2007; 19:454-65. [PMID: 17085014]
  57. Tang X, Molina M, Amar S. p53 short peptide (p53pep164) regulates lipopolysaccharide-induced tumor necrosis factor-alpha factor/cytokine expression. Cancer Res. 2007; 67:1308-16. [PMID: 17283168]
  58. Guicciardi ME, Leist M, Gores GJ. Lysosomes in cell death. Oncogene. 2004; 23:2881-90. [PMID: 15077151]
  59. Ferri KF, Kroemer G. Organelle-specific initiation of cell death pathways. Nat Cell Biol. 2001; 3:E255-63. [PMID: 11715037]
  60. Kagedal K, Zhao M, Svensson I, Brunk UT. Sphingosine-induced apoptosis is dependent on lysosomal proteases. Biochem J. 2001; 359:335-43. [PMID: 11583579]
  61. Turk B, Turk D, Turk V. Lysosomal cysteine proteases: more than scavengers. Biochim Biophys Acta. 2000; 1477:98-111. [PMID: 10708852]
  62. Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J. Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004; 5:R80 [PMID: 15461798]