Pathway analysis for Gene Cluster s 1 and 2 ranked by p-value. The number of objects refers to the enrichment of genes/total gene related objects in the map.
1. Introduction
1.1. Primary Open Angle Glaucoma (POAG)
Elevated intraocular pressure (IOP) is the most important risk factor in POAG. Typically affecting older adults, in POAG the IOP exceeds the level that is tolerated by that individual’s optic nerve head (ONH). However, many individuals with clinical ocular hypertension do not exhibit glaucomatous changes in the optic disk; whereas, some individuals will develop glaucomatous changes at clinically normal IOP levels. These clinical findings indicate that individual variability in susceptibility of the ONH to IOP is an important factor in glaucomatous optic neuropathy. However, the molecular and cellular factors that may underlie variability in susceptibility of the ONH to elevated IOP have not been elucidated.
The target of the mechanical stress generated by elevated IOP is the lamina cribrosa in the ONH (Bellezza et al., 2003). In the glaucomatous ONH, compression, stretching, and remodeling of the cribriform plates of the lamina cribrosa occur. In many POAG patients, these elevated IOP related changes result in remodeling of the extracellular matrix (ECM), altering the quantity and composition of several ECM macromolecules that significantly affect the biomechanical properties of the tissue supporting the nerve fibers (Hernandez, 2000; Bellezza et al., 2003). Alterations of ECM components of the ONH in glaucoma perhaps sets the stage for further optic nerve damage from IOP during the progression of the disease.
1.2. POAG in African Americans
Epidemiological and genetic studies indicate that ethnic/genetic background plays an important role in susceptibility to POAG. POAG is more prevalent in Black Americans of African ancestry (AA) than in White Americans of European ancestry (CA), with reported frequencies of 3-4% in the AA population over the age of 40 years, as compared with approximately 1% in CA populations (Friedman et al., 2004). The disease is particularly frequent in Afro-Caribbean populations, with prevalence of 7% in Barbados and 8.8% in St. Lucia (Nemesure et al., 2001). On average AAs have increased duration (Quigley and Vitale, 1997) and progression of disease (Broman et al., 2008) compared to other populations. A positive family history of POAG is a major risk factor for the disease in AA (Leske et al., 2008). The Advanced Glaucoma Intervention Study (AGIS), which compared the glaucoma outcomes in AA and CA patients, concluded that after failure of medical therapy, surgical trabeculectomy delayed progression of glaucoma more effectively in CA than in AA patients (Beck, 2003; Ederer et al., 2004). AAs have significantly larger disc areas, larger cup areas, larger cup-to-disc ratios and smaller neural rim area-to-disc area ratios compared with CAs (Varma et al., 1994; Quigley et al., 1999). A morphometric study determined that AAs have a larger total area of the lamina cribrosa and a greater number of pores than CAs (Dandona et al., 1990). The cellular and molecular bases for these anatomic differences and disease prevalences have not been explored.
Our studies focus on how astrocytes in the lamina cribrosa of the human optic nerve head contribute to glaucomatous optic neuropathy that is associated with elevated IOP. There is mounting evidence that astrocytes are responsible for many of the pathological changes in the glaucomatous ONH and cells isolated from donors of different ethnic backgrounds vary in their gene expression (Miao et al., 2008; Hernandez et al., 2008). Cellular explanations for these population-related differences may be realized using primary ONH astrocytes as a source of differing cellular phenotypes and an
1.3. Astrocytes in the ONH
Astrocytes are the major glial cell type in the non-myelinated human ONH and provide cellular support functions to the axons while interfacing between connective tissue surfaces and surrounding blood vessels. In the normal ONH, astrocytes are quiescent, terminally differentiated cells. In the lamina cribrosa, quiescent astrocytes form lamellae oriented perpendicular to the axons surrounding a core of fibroelastic extracellular matrix (Hernandez, 2000). Astrocytes supply energy substrates to axons in the optic nerve and maintain extracellular pH and ion homeostasis in the periaxonal space (Fields and Stevens-Graham, 2002). ONH astrocytes express ECM proteins such as laminin, and proteoglycans, as well as bone morphogenetic proteins (Zode et al., 2007; Wordinger et al., 2002) neurotrophins and receptors (Lambert et al., 2004; Yang et al., 2007). Several of these may serve as neuroprotective factors for retinal ganglion cells (RGC). In addition to astrocytes, other cell types exist in the lamina cribrosa of humans and non human primates, including microglia, vascular endothelia and the lamina cribrosa cell. Lamina cribrosa cells can be distinguished from astrocytes because they do not express glial fibrillary acidic protein (GFAP) and they do not express vascular or microglial markers (Hernandez et al., 1988; Kirwan et al., 2005; Hernandez et al., 1989).
1.4. Reactive ONH astrocytes in glaucoma
Adult, quiescent astrocytes become "reactive" after injury or disease and participate in formation of a glial scar, which does not support axonal survival or growth (Sofroniew, 2005; Mc Graw et al., 2001; Hatten et al., 1991).The major hallmarks of a reactive astrocyte are an enlarged cell body and a thick network of processes with increased expression of GFAP and vimentin (Hatten et al., 1991). Similarly, reactive astrocytes in the glaucomatous ONH are large rounded cells with many thick processes and express increased amounts of GFAP, vimentin and HSP27 (Hernandez et al., 2008). The astrocytes are motile and migrate either to the edge of the laminar plates or to inside the nerve bundles (Hernandez et al., 2008). In glaucoma, reactive astrocytes exhibit putative neurodestructive and neuroprotective cellular cascades in the ONH (Hernandez and Pena, 1997). Previous studies demonstrated that the genomic responses of astrocytes in the glaucomatous ONH serve as the basis for these cascades (Hernandez et al., 2002). Damage to retinal ganglion cell axons and remodeling of the connective tissue plates in the lamina cribrosa appear to be mediated by astrocytes (Hernandez et al., 2008; Nickells, 2007). Reactive ONH astrocytes increase expression of various cell surface and extracellular matrix-related proteins such as laminin, tenascin C, and proteoglycans that play important roles in cell-cell recognition and in cell adhesion (Hernandez et al., 2002; Hernandez and Pena, 1997; Pena et al., 1999b). As might be expected there are also changes in signal transduction in reactive ONH astrocytes (Review (Hernandez et al., 2008)). For example, EGF receptors (Liu et al., 2006) and endothelin B receptors (Wang et al., 2006) are upregulated
1.5. Susceptibility to glaucoma: Differential expression in the ONH astrocyte transcriptomes from glaucomatous AA and CA donors
Gene expression in primary cultures of ONH astrocytes obtained from age-matched normal and glaucomatous donors of CA and AA populations was done using Affymetrix GeneChip microarrays. Gene Ontology analysis and networks of interacting proteins were constructed using the BioGRID database (Stark et al., 2006). The differential gene expression data were distributed among three networks that include regulation of myosin, actin, and protein trafficking (Lukas et al., 2008). Remarkably, cultured glaucomatous astrocytes retain differential expression of genes that promote cell motility and migration, regulate cell adhesion, and are associated with structural tissue changes during neurodegeneration. Similar changes in gene expression were observed in glaucomatous optic nerve head tissues as assessed by immunohistochemistry (Lukas et al., 2008). Thus, the
1.6. Model systems for glaucomatous astrocytes
One of the significant problems associated with work on glaucomatous astrocytes is that their availability is quite limited as the number of donors with glaucoma is much smaller than donors without eye disease. Therefore, we investigated methods to “transform” normal ONH astrocytes into a preglaucoma or glaucoma-like phenotype. In principle, mechanical stress applied to the cells by way of physical contact (as in the ONH) should mimic the
2. Results-gene expression studies in ONH astrocytes
2.1. Elevated hydrostatic pressure model system to simulate elevated IOP
We investigated the molecular bases of elevated hydrostatic pressure (HP) responses in primary cultures of ONH astrocytes derived from AA donors compared to CA donors. This model system simulates the mechanical stresses placed upon the optic nerve head by elevated IOP. Using high density microarrays we investigated global changes in gene expression in a cohort of three astrocyte cell lines from each group. We validated changes in expression induced by elevated HP in three to five additional AA and CA astrocyte lines using quantitative RT-PCR and/or Western Blotting. Using a global phosphoproteome approach, we also identified changes in protein phosphorylation associated with elevated HP. In both gene expression and proteome experiments, ONH astrocytes were subjected to elevated HP for periods of 3, 6, 24, and 48 hours. Control cells were cultured under ambient pressure (CP) for 6 or 48 hr. A custom made pressure chamber was used to subject cultured cells (6 well plates) to 60 mm (above atmospheric pressure) of HP for the desired periods of time. In order to compensate for potential changes in dissolved oxygen or pH, the gas mixture in the chamber was adjusted to 8% CO2 and buffering capacity in the media increased (Yang et al., 2004). Under these conditions changes in dissolved oxygen and pH are negligible during the course of the experiment (Ricard et al., 2000). Messenger RNA and protein were extracted from the cells immediately after the experiment using standard published techniques (Lukas et al., 2008; Miao et al., 2010; Chen et al., 2009). Gene expression was measured using Illumina Human-6 beadchips following the manufacturer’s protocols. These microarray chips have ~47,000 features representing ~31,000 human genes and expressed isoforms.
2.2. Microarray data analysis
To best utilize the unique features of Illumina BeadArray technology, we used the Bioconductor lumi package (Du et al., 2008; Lin et al., 2008; Du et al., 2007) to preprocess the raw data output by Illumina Beadstudio software. The data was preprocessed using variance stabilization transformation method (Lin et al., 2008) followed by quantile normalization. Probes with all samples “Absent” (lower or around background levels) were removed from further analysis to reduce false positives.
Two phases of analysis were done. First, to identify differentially expressed genes common to both the AA and CA groups, we applied routines implemented in limma package (Du et al., 2008) to fit linear models to the normalized expression values. The variance used in the t-score calculation was corrected by an empirical Bayesian method (Smyth, 2004) for better estimation with a small sample size. For each time point (3h, 6h, 24h and 48h), we compared all HP samples with the control samples, and defined the differentially expressed genes as fold-change higher than 1.3-fold (p<0.01). After this filtering we found 352 genes differentially expressed at least at one time point. In order to reduce variations across different clinical samples, we first normalized the pressure samples at each time point against the control sample from the same eye, then we standardized the normalized expression profile of each gene as one standard deviation and without centering (the fold-change direction is unchanged).
Hierarchical clustering was first applied to get an overview of the cluster distribution, and then we defined initial clusters with visually selected cutoff. Using these initial sets, we did K- Means clustering resulting in 4 clusters that include 67, 137, 56 and 93 genes, respectively. The average changes in expression for each cluster as a function of time of HP treatment are shown in Figure 1. The cluster profile is defined as the average of all gene time profiles in the cluster (solid lines in Figure 1). Clusters 1 and 4 exhibit a decreasing trend of expression, while clusters 2 and 3 have increased expression as a function of time of exposure to elevated HP. The profiles of clusters 3 and 4 exhibit a flattening at the 3 and 6 hr times followed by a continued increase or decrease at 24 and 48 hr. These data suggest that ONH astrocytes respond to elevated HP in a dynamic fashion that involves a fast (3-6 hr) and slower (or protracted) response that follows at 24-48 hr. These responses involve genes in different signaling pathways that define a transition to new cellular phenotypes. As we will demonstrate in subsequent sections, there is overlap between the pressure-induced phenotype and glaucomatous cellular profiles that have been described previously (Hernandez et al., 2002; Nickells, 2007). Similarly, we found that extended HP treatment of ONH astrocytes (72-120 hr) induces further changes in cellular phenotype with respect to the expression of genes that affect cell motility (Miao et al., 2010) and cell morphology (Yang et al., 1993).
2.3. HP-induced cellular processes in ONH astrocytes
For consistent comparison to the glaucomatous phenotype for ONH astrocytes (Nickells, 2007) we used the GeneGo suite to analyze the differential expression data for our elevated HP data. Using the gene clusters identified above, the genes were mapped onto existing canonical pathways that have been defined within GeneGo (Nikolsky et al., 2005). There are 650 such maps within the GeneGo library. The overrepresentation of genes within a pathway gives a p-value that estimates the significance of the mapped genes. The four top scoring maps for each cluster are summarized in Tables 1 and 2.
Map | Cell Process | p-value | Objects | |
Estradiol Metabolism | 0.00007 | 3/19 | ||
Regulation of Triiodothryronine and thyroxine signaling | Transcription | 0.0054 | 2/31 | |
Ligand dependent transcription of retinoid target genes | Transcription | 0.0047 | 2/32 | |
Androgen receptor nuclear signaling | Transcription | 0.0094 | 2/64 | |
Cell cycle-ESR1 regulation of G1/S transition | Response to hormone, extracellular stimulus | 0.00001 | 6/32 | |
BMP Signaling/ TGFβ inhibition | Cytokine mediated signaling | 0.00033 | 3/33 | |
Cell cycle regulation | Cell cycle | 0.0063 | 3/29 | |
Chemokines and cell adhesion | Cell Adhesion | 0.0079 | 3/50 |
In the first two clusters, the activation of transcription and cytokine mediated signaling are among the most significant pathways. The transcriptional processes decline with increased time of HP suggesting that this is an early response that stabilizes to near control levels at 48 hr. As will be shown later this response also correlates with changes in protein phosphorylation. The other significant pathways in Clusters 1 and 2 are cytokine-related or other extracellular stimuli. These increase with time of HP treatment suggesting a transition to a phenotype associated with increased sensitivity to extracellular agents. In this regard, this cluster includes serotonin receptor 5B (Table 5), a receptor known to be upregulated by mechanical stress (Liang et al., 2006; Sanden et al., 2000).
Cluster 3 genes are found in pathways leading to cytoskeleton remodeling, transport, and other signaling processes. ONH astrocytes undergo changes in morphology after being subjected to elevated HP (Yang et al., 1993) and up-regulation of cytoskeletal signaling processes is consistent with these properties. Moreover, altered cytoskeletal, transport, and cell adhesion gene expression is also found in populations of glaucomatous astrocytes (Lukas et al., 2008).
Map | Cell Process | p-value | Objects | |
Cell contraction mediated by G-protein coupled receptors | Cytoskeleton | 0.0019 | 3/15 | |
Transport_RAB1A regulation pathway | ER-Golgi transport | 0.0004 | 2/12 | |
Cytoskeleton remodeling by Rho GTPases | Cytoskeleton | 0.0053 | 4/23 | |
Regulation of insulin pathway signaling | Receptor-mediated signaling | 0.0050 | 2/42 | |
Transcription – Sin3 and NuRD in transcriptional regulation | Transcription | 0.00059 | 3/38 | |
G-protein receptor signaling via cyclic AMP pathways | Extracellular mediated signaling | 0.0022 | 2/52 | |
Regulation of lipid metabolism by transcription factor RXR via PPAR, RAR, and VDR pathways | Transcription | 0.0075 | 2/30 | |
Slit-Robo signaling in regulation of neuronal axon guidance | Development | 0.0075 | 2/30 |
Cluster 4 genes populate a distinct set of pathways in transcription, metabolism, and signal transduction. Inhibition of axonal guidance is altered under elevated HP due to secretion of Slit2/3 (Table 3) by astrocytes decreased with elevated HP over time. Further decreases in transcriptional pathways are evident in Cluster 4 with changes in the retinoid (RXR) based regulatory pathways. The latter involves changes in transcription through the post-translational modification of histones (Figure 2). Thus, besides the decrease in expression histones can undergo acetylation and methylation of lysine residues and phosphorylation of serine/threonine residues that alter chromatin structure and the transcription of specific DNA regions (Oki et al., 2007; Ito, 2007). As shown below we detected phosphorylation of lysine methyltransferase enzymes in ONH astrocytes at early time points of treatment with elevated HP. This is particularly interesting because there are few reports of epigenetic changes being induced by mechanical stress (Illi et al., 2005) and none by elevated HP. Thus, the transcriptome of ONH astrocytes may be undergoing reprogramming to accommodate this presence of this stress.
Other pathways represented in Cluster 4 genes that are altered by elevated HP include cyclic-AMP signaling. Our previous studies of short term elevated HP (30 min – 3 hr) treatment of ONH astrocytes revealed that cyclic-AMP signaling is differentially affected in astrocytes from African American donors (Chen et al., 2009). Upregulation of adenylate cyclases and parathyroid hormone-like hormone (PTHLH) that activates G-protein signaling to adenylate cyclases were detected (Chen et al., 2009). In the current work with longer term (24-48 hr) treatment both AA and CA astrocytes down regulate phosphodiesterases (PDE3A, PDE7B (Table 3.) which may result in higher levels of cyclic nucleotides and increased activity of transcription factors such as CREB that are activated by cyclic AMP-dependent phosphorylation (Johannessen et al., 2004; Delghandi et al., 2005).
Symbol | Description | HP3 | HP6 | HP24 | HP48 |
interleukin 6 (interferon, beta 2) | 2.22 | 1.41 | -1.15 | -1.67 | |
TAF5-like RNA polymerase II, p300/CBP-associated factor (PCAF) | 1.59 | 1.24 | 1.44 | 1.38 | |
phosphoprotein associated with glycosphingolipid microdomains 1 | 1.44 | 1.11 | -1.11 | 1.06 | |
wingless-type MMTV integration site family, member 5B | 1.38 | 1.26 | 1.22 | 1.08 | |
E74-like factor 4 (ets domain transcription factor) | 1.37 | 1.15 | 1.10 | -1.05 | |
itchy E3 ubiquitin protein ligase homolog | 1.35 | 1.35 | 1.18 | 1.10 | |
sulfotransferase family, cytosolic, 1A 3 | 1.33 | 1.18 | 1.24 | 1.03 | |
nuclear receptor corepressor 2 | 1.18 | 1.34 | -1.02 | 1.02 | |
FBJ osteosarcoma viral oncogene homolog | -1.21 | -1.74 | -1.58 | -1.44 | |
v-myc myelocytomatosis viral oncogene homolog | -1.39 | -1.37 | -1.08 | -1.03 | |
S-phase kinase-associated protein 2 (p45) | -1.30 | -1.09 | -1.00 | 1.03 | |
cAMP responsive element binding protein 1 | -1.32 | -1.42 | -1.21 | -1.16 | |
oncostatin M receptor | -1.16 | -1.23 | -1.53 | -1.14 | |
integrin, alpha 6 | -1.39 | -1.34 | -1.76 | -1.13 | |
paired-like homeodomain 1 | -1.44 | -1.25 | -1.27 | -1.37 | |
bone morphogenetic protein 4 | -1.79 | -1.35 | 1.04 | 1.06 | |
SMAD family member 5 | -1.34 | -1.04 | -1.00 | -1.22 | |
TGFB-induced factor homeobox 1 | -1.10 | -1.37 | 1.05 | 1.19 | |
collagen, type IV, alpha 4 | -1.56 | -1.63 | -1.32 | -1.88 | |
cofilin 2 | -1.12 | -1.67 | 1.00 | 1.24 | |
5-hydroxytryptamine receptor 2B | -1.61 | -1.69 | 1.05 | 1.16 | |
myosin, light chain 9, regulatory | 1.05 | 1.04 | 1.11 | 1.42 | |
phosphoinositide-3-kinase interacting protein 1 | -1.04 | 1.22 | 1.04 | 1.39 | |
insulin receptor substrate 2 | -1.04 | -1.35 | 1.36 | 1.43 | |
myosin ID | 1.06 | 1.12 | 1.34 | 1.22 | |
N-ethylmaleimide-sensitive factor attachment protein, alpha | -1.02 | -1.02 | 1.21 | 1.31 | |
RUN and SH3 domain containing 2 | 1.17 | 1.03 | 1.39 | 1.23 | |
ribosomal protein S6 kinase, 70kDa | 1.10 | 1.17 | 1.15 | 1.39 | |
protein kinase N1 | -1.28 | -1.08 | 1.03 | 1.32 | |
peroxisome proliferator-activated receptor gamma | -1.46 | -1.28 | -1.46 | -1.28 | |
phosphodiesterase 7B | -1.27 | -1.20 | -2.09 | -2.06 | |
histone cluster 1, H2bg | -1.26 | -1.19 | -1.50 | -1.31 | |
thyroid hormone receptor, alpha | -1.24 | -1.23 | -1.30 | -1.56 | |
histone cluster 1, H2bd | -1.23 | -1.34 | -1.63 | -1.77 | |
histone cluster 1, H4h | -1.18 | -1.31 | -1.38 | -1.50 | |
histone cluster 1, H4e | -1.16 | -1.17 | -1.33 | -1.54 | |
retinoic acid receptor, beta | -1.14 | 1.03 | -1.39 | -1.28 | |
histone cluster 2, H4a | -1.11 | -1.29 | -1.45 | -1.48 | |
slit homolog 2 (Drosophila) | -1.08 | -1.07 | -1.22 | -1.56 | |
transforming growth factor, beta receptor | -1.07 | -1.08 | -1.26 | -1.37 | |
histone cluster 1, H3d | -1.05 | -1.25 | -1.52 | -1.49 | |
histone cluster 2, H2aa3 | 1.04 | -1.01 | -1.31 | -1.45 | |
phosphodiesterase 3A, cGMP-inhibited | 1.00 | -1.62 | -1.24 | -1.54 |
2.4. Differential gene expression in ONH astrocytes from African American compared to Caucasian donors after treatment with elevated HP
Using the same gene selection and filtering criteria described above for the combined AA and CA data sets, we performed a cross comparison of gene expression at the same time points for genes differentially expressed by AA astrocytes. The 63 genes in Figure 3 were selected by satisfying these criteria: 1. a fold-change > 1.3, p-value < 0.01 of HP-AA vs. Ctrl at least at one time point, 2. a p-value < 0.01 of normalized HP-AA vs. normalized HP-CA at least at one time point (each time point was normalized against control).
To cluster the genes, we first standardized the normalized expression profile of each gene to one standard deviation and without centering (to keep the fold-change direction unchanged), then applied hierarchical cluster to get an overview of the cluster distribution. Based on the hierarchical clustering results, we defined initial clusters with visually selected cutoff value (5.5 for our case). Using these initial clusters, we did K-Means clustering. A total of 6 clusters were obtained. Two representative gene expression cluster patterns over the time course are shown in Figure 4. The early response genes (Figure 4A) increased significantly in AA (in red dotted line) after 3h exposure to HP compared to CA (black solid line). The late response genes (Figure 4B) increased significantly in AA (red dotted line) after 48h exposure to HP compared to CA (black solid line).
2.5. Validation of gene expression changes in ONH astrocytes subjected to elevated HP
Gene validation experiments were done primarily on the genes differentially expressed in the AA compared to CA ONH astrocytes as these were most likely to define some of the differences in sensitivity of AA cells to elevated HP. For validation, we tested GPNMB, CTSK, GCLM, HBEGF, and PLOD2. GPNMB is a late response gene. It has a requisite mutation (Anderson et al., 2008) in the hereditary glaucoma mouse model DBA-2J and is differentially expressed in ONH astrocytes in a primate ocular hypertension model (Kompass et al., 2008). Increased mRNA levels of GPNMB in AA correlate with the length of HP exposure (Figure 5A). Another late HP response gene is CTSK which is a lysosomal cysteine proteinase with strong degradative activity against the extracellular matrix and is suggested to play a role in tumor invasiveness. CTSK was significantly upregulated in AA
astrocytes after 48 hr exposure to HP (Figure 5A). Two down regulated genes in the AA astrocytes compared to CA are HBEGF and PLOD2. In earlier work we found that HBEGF and PLOD2 were expressed significantly lower in native AA astrocytes compared to CA (Miao et al., 2008). When exposed to HP, HBEGF expression was further downregulated in AA, while upregulated in CA astrocytes (Figure 5B). PLOD2 encodes an enzyme which catalyzes the hydroxylation of lysyl residues in collagen-like peptides. It is important for the stability of intermolecular crosslinks in basement membranes. These data support our idea that the two populations of ONH astrocytes respond differently to elevated hydrostatic pressure with respect to changes in gene expression. Some responses may be due to intrinsically lower expression in the AA astrocyte population that is enhanced by elevated HP, while others appear to be specifically up-regulated in these cells.
2.6. Phosphoproteome profiling of ONH astrocytes subjected to elevated HP
Global phosphoproteome profiling was done using standard affinity-based methods to trap phosphopeptides coupled to high performance liquid chromatography-mass spectrometry (LC-MS) for detection and identification similar to our previous studies of changes in protein phosphorylation in the retina following optic nerve crush (Lukas et al., 2009). The strategy here was to determine whether changes in gene expression can be correlated with changes in protein phosphorylation in ONH astrocytes subjected to elevated HP. These studies were done on two pairs of ONH astrocytes- one pair from AA donors and the other
from CA donors. As with the initial gene expression analysis, we combined all of the data at each time point to generate lists of the detected phosphoproteins. This was not a quantitative study, but rather a sampling of the phosphoproteome at each time point. This unbiased LC-MS method detects the most abundant phosphopeptides in the sample resulting in a limited survey of phosphoproteins. Phosphoproteome data were obtained at 3, 6, 12, and 24hr of elevated HP. These sets were individually submitted to Gene Ontology analysis using the GoMiner program (Zeeberg et al., 2003). During the time course of elevated HP, the number of detected nuclear and cell-cycle related phosphoproteins decreases and the cytoskeletal-associated phosphoproteins increase so as to make the differential enrichment significant (Table 4).
HP 3- 6hr | |
Accession # | |
IPI00448121 | |
IPI00855998 | |
IPI00103595 | |
IPI00297851 | |
IPI00009724 | |
IPI00015526 | |
IPI00171903 | |
IPI00456887 | |
IPI00299904 | |
IPI00332499 | |
IPI00470429 | |
IPI00031627 | |
IPI00783392 | |
IPI00024710 | |
IPI00004312 | |
IPI00044681 | |
IPI00783017 | |
IPI00003003 | |
IPI00022177 | |
IPI00166972 | |
IPI00747060 | |
IPI00640320 | |
IPI00009724 | |
IPI00410079 | |
IPI00171903 | |
IPI00307733 | |
IPI00025158 | |
IPI00749005 | |
IPI00170921 | |
IPI00217407 | |
IPI00796117 | |
IPI00298731 | |
Summarized in Tables 5 and 6 are selected phosphoproteins identified in the early (3-6 hr) and late (24-48 hr) phases of elevated HP. Within the nuclear/cell cycle groups of phosphoproteins we find transcription factors, nuclear riboproteins, and ubiquitin ligases. Also notable are histone methyl transferase proteins that may be functioning to reprogram the ONH astrocyte gene expression in response to the stress induced by elevated HP. Phosphorylation of these proteins may function to activate epigenetic reprogramming of the ONH astrocyte in response to the stress induced by elevated HP.
HP 3-6 hr | |
Accession# | |
IPI00472779 | |
IPI00011219 | |
IPI00008756 | |
IPI00306929 | |
IPI00479962 | |
IPI00307155 | |
IPI00289639 | |
IPI00010448 | |
IPI00457243 | |
IPI00792788 | |
IPI00216408 | |
IPI00642259 | |
IPI00260090 | |
IPI00452247 | |
IPI00848334 | |
IPI00306929 | |
IPI00479962 | |
IPI00303335 | |
IPI00071509 |
Within the cytoskeletal group were several proteins associated with focal adhesions (Table 6). These include for example, DOCK5/9 and ELMO1/2. ELMO was detected at the 3hr HP time point and the detected phosphopeptide, RIAFDAESEPNNSSGpSMEKR, corresponds to residues 329-348 of ELMO1. ELMO1 interacts with DOCK-180 and other proteins in a large complex found at focal adhesions (Beausoleil et al., 2006). Specific antiphosphopeptide antibodies to ELMO1 were not available, but there are phosphorylation site antibodies to p130CAS, another known phosphoprotein in the DOCK protein complex (Menniti et al., 2006). p130CAS is known to become highly phosphorylated in response to mechanical stress (Sawada et al., 2006; Geiger, 2006). We found that p130CAS phosphorylation increases with HP, particularly at the 3 and 6 hr time points in AA astrocytes and at 24-48 hr in CA astrocytes (Figure 6). Thus, the AA astrocytes respond earlier to elevated HP than the CA astrocytes with respect to p130CAS phosphorylation. Phosphorylation of p130CAS and ELMO may modulate the dynamics of focal adhesions. To support this idea, we found that
phosphorylation of focal adhesion kinase (FAK) was increased in CA astrocytes at 24-48 hr of elevated HP (Figure 7). Although FAK was phosphorylated in AA cells, its phosphorylation level was constant at each time point during elevated HP treatment (not shown).
2.7. Correlations of gene expression and protein phosphorylation in ONH astrocytes
The pattern of cytoskeletal protein phosphorylation is consistent with the differences in gene expression found in the late (24-48 hr) time of elevated HP treatment (Tables 1-3). Similarly, the pattern of protein phosphorylation at the early (3-6 hr) time of HP treatment is consistent with the transcriptional activity associated with this phase (Tables 1-3). Thus, the elevated HP system provides a potential “preglaucoma” state in ONH astrocytes because differential changes in expression of genes in the TGFβ pathways and cytoskeletal regulation are found in glaucomatous AA and CA astrocytes and in the glaucomatous ONH (Lukas et al., 2008). As shown earlier, elevated HP increases the expression of GPNMB protein in AA astrocytes compared to CA astrocytes (Figure 5C) consistent with the gene expression data. GPNMB is thought to function as a cell adhesion protein that connects the extracellular environment to downstream cellular signaling linked to gene transcription via TGFβ pathways that modulates the synthesis of extracellular matrix proteins in ONH astrocytes (Fuchshofer et al., 2005). Therefore, these pathways may be larger contributors to changes in cell adhesion in AA astrocytes compared to CA. On the other hand, the activation of focal adhesions complexes is linked to cytoskeletal remodeling and transcription through multiple pathways (Koyama et al., 2000). Thus, in signal transduction and gene expression, AA and CA astrocytes exhibit differential responses to elevated HP that can produce altered cell adhesion and morphology.
3. Conclusions
The types of experiments described in this report fall into the realm of systems biology. The guiding hypothesis is that subjecting a system such as cultured ONH astrocytes, to a stressor (elevated HP) will cause changes in gene expression and protein phosphorylation in a fashion that impacts multiple signaling networks and processes. We found that, indeed, elevated HP induces changes in gene expression patterns that are setting up new cellular phenotypes. Thus, the system is a model for what elevated IOP may induce in ONH astrocytes
In the same realm, changes in protein phosphorylation may define dynamic cellular phenotypes which have characteristics of different cellular populations. The few differences we uncovered in protein phosphorylation of cytoskeletal proteins between CA and AA astrocytes after treatment with elevated IOP are likely just a small sampling of the entire signaling system. However, one of the exciting prospects revealed in this study is that elevated IOP might be inducing epigenetic changes in ONH astrocytes that are responsible, in part, for altering gene transcription programs. This opens up unexplored areas for future glaucoma research.
Acknowledgments
This work is dedicated to the memory of our colleague Dr. M. Rosario Hernandez. It was supported, in part, by a grant from the National Institutes of Health (EY06416) and unrestricted funds from Research to Prevent Blindness.
References
- 1.
Anderson M. G. Nair K. S. Amonoo L. A. Mehalow A. Trantow C. M. Masli S. John S. W. 2008 GpnmbR150X allele must be present in bone marrow derived cells to mediate DBA/2J glaucoma. BMC Genet. 9,30 EOF - 2.
Beausoleil S. A. Villen J. Gerber S. A. Rush J. Gygi S. P. 2006 A probability-based approach for high-throughput protein phosphorylation analysis and site localization. Nat. Biotechnol.24 1285 1292 - 3.
Beck,A.D. 2003 Review of recent publications of the Advanced Glaucoma Intervention Study. Curr. Opin. Ophthalmol.14 83 85 - 4.
Bellezza A. J. Rintalan C. J. Thompson H. W. Downs J. C. Hart R. T. Burgoyne C. F. 2003 Deformation of the lamina cribrosa and anterior scleral canal wall in early experimental glaucoma. Invest Ophthalmol. Vis. Sci.44 623 637 - 5.
Broman A. T. Quigley H. A. West S. K. Katz J. Munoz B. Bandeen-Roche K. Tielsch J. M. Friedman D. S. Crowston J. Taylor H. R. Varma R. Leske M. C. Bengtsson B. Heijl A. He M. Foster P. J. 2008 Estimating the rate of progressive visual field damage in those with open-angle glaucoma, from cross-sectional data. Invest Ophthalmol. Vis. Sci.49 66 76 - 6.
Chen L. Lukas T. J. Hernandez M. R. 2009 Hydrostatic pressure-dependent changes in cyclic AMP signaling in optic nerve head astrocytes from Caucasian and African American donors. Mol. Vis.15 1664 1672 - 7.
Dandona L. Quigley H. A. Brown A. E. Enger C. 1990 Quantitative regional structure of the normal human lamina cribrosa. A racial comparison. Arch. Ophthalmol.108 393 398 - 8.
Delghandi M. P. Johannessen M. Moens U. 2005 The cAMP signalling pathway activates CREB through PKA,38 and MSK1 in NIH 3T3 cells. Cell Signal. 17, 1343-1351. - 9.
Du P. Kibbe W. A. Lin S. M. 2007 nuID: a universal naming scheme of oligonucleotides for illumina, affymetrix, and other microarrays. Biol. Direct. 2,16 EOF - 10.
Du P. Kibbe W. A. Lin S. M. 2008 lumi: a pipeline for processing Illumina microarray. Bioinformatics.24 1547 1548 - 11.
Ederer F. Gaasterland D. A. Dally L. G. Kim J. Van Veldhuisen P. C. Blackwell B. Prum B. Shafranov G. Allen R. C. Beck A. 2004 The Advanced Glaucoma Intervention Study (AGIS): 13. Comparison of treatment outcomes within race: 10-year results. Ophthalmology111 651 664 - 12.
Feng Y. Yang J. H. Huang H. Kennedy S. P. Turi T. G. Thompson J. F. Libby P. Lee R. T. 1999 Transcriptional profile of mechanically induced genes in human vascular smooth muscle cells. Circ. Res.85 1118 1123 - 13.
Fields R. D. Stevens-Graham B. 2002 New insights into neuron-glia communication. Science298 556 562 - 14.
Friedman D. S. Wolfs R. C. O’Colmain B. J. Klein B. E. Taylor H. R. West S. Leske M. C. Mitchell P. Congdon N. Kempen J. 2004 Prevalence of open-angle glaucoma among adults in the United States. Arch. Ophthalmol.122 532 538 - 15.
Fuchshofer R. Birke M. Welge-Lussen U. Kook D. Lutjen-Drecoll E. 2005 Transforming growth factor-beta 2 modulated extracellular matrix component expression in cultured human optic nerve head astrocytes. Invest Ophthalmol. Vis. Sci.46 568 578 - 16.
Geiger,B. 2006 A role for130Cas in mechanotransduction. Cell , 879-881. - 17.
Hashimoto K. Parker A. Malone P. Gabelt B. T. Rasmussen C. Kaufman P. S. Hernandez M. R. 2005 Long-term activation of c-Fos and c-Jun in optic nerve head astrocytes in experimental ocular hypertension in monkeys and after exposure to elevated pressure in vitro. Brain Res.1054 103 115 - 18.
Hatten M. E. Liem R. K. Shelanski M. L. Mason C. A. 1991 Astroglia in CNS injury. Glia4 233 243 - 19.
Hernandez,M.R. 2000 The optic nerve head in glaucoma: role of astrocytes in tissue remodeling. Prog. Retin. Eye Res.19 297 321 - 20.
Hernandez M. R. Agapova O. A. Yang P. Salvador-Silva M. Ricard C. S. Aoi S. 2002 Differential gene expression in astrocytes from human normal and glaucomatous optic nerve head analyzed by cDNA microarray. Glia38 45 64 - 21.
Hernandez M. R. Igoe F. Neufeld A. H. 1988 Cell culture of the human lamina cribrosa. Invest Ophthalmol. Vis. Sci.29 78 89 - 22.
Hernandez M. R. Luo X. X. Andrzejewska W. Neufeld A. H. 1989 Age-related changes in the extracellular matrix of the human optic nerve head. Am. J. Ophthalmol.107 476 484 - 23.
Hernandez M. R. Miao H. Lukas T. 2008 Astrocytes in glaucomatous optic neuropathy. Prog. Brain Res.173 353 373 - 24.
Hernandez M. R. Pena J. D. 1997 The optic nerve head in glaucomatous optic neuropathy. Arch. Ophthalmol.115 389 395 - 25.
Hernandez M. R. Pena J. D. Selvidge J. A. Salvador-Silva M. Yang P. 2000 Hydrostatic pressure stimulates synthesis of elastin in cultured optic nerve head astrocytes. Glia32 122 136 - 26.
Illi B. Scopece A. Nanni S. Farsetti A. Morgante L. Biglioli P. Capogrossi M. C. Gaetano C. 2005 Epigenetic histone modification and cardiovascular lineage programming in mouse embryonic stem cells exposed to laminar shear stress. Circ. Res.96 501 508 - 27.
Ito,T. 2007 Role of histone modification in chromatin dynamics. J. Biochem. (Tokyo)141 609 614 - 28.
Johannessen M. Delghandi M. P. Moens U. 2004 What turns CREB on? Cell Signal.16 1211 1227 - 29.
Johnson E. C. Jia L. Cepurna W. O. Doser T. A. Morrison J. C. 2007 Global changes in optic nerve head gene expression after exposure to elevated intraocular pressure in a rat glaucoma model. Invest Ophthalmol. Vis. Sci.48 3161 3177 - 30.
Kirwan R. P. Crean J. K. Fenerty C. H. Clark A. F. O’Brien C. J. 2004 Effect of cyclical mechanical stretch and exogenous transforming growth factor-beta1 on matrix metalloproteinase-2 activity in lamina cribrosa cells from the human optic nerve head. J. Glaucoma.13 327 334 - 31.
Kirwan R. P. Leonard M. O. Murphy M. Clark A. F. O’Brien C. J. 2005 Transforming growth factor-beta-regulated gene transcription and protein expression in human GFAP-negative lamina cribrosa cells. Glia52 309 324 - 32.
Kobayashi S. Vidal I. Pena J. D. Hernandez M. R. 1997 Expression of neural cell adhesion molecule (NCAM) characterizes a subpopulation of type 1 astrocytes in human optic nerve head. Glia20 262 273 - 33.
Kompass K. S. Agapova O. Li W. Kaufman P. L. Rasmussen C. Hernandez M. R. 2008 Bioinformatic and statistical analysis of the optic nerve head in a primate model of ocular hypertension. BMC Neuroscience in press. - 34.
Koyama Y. Yoshioka Y. Hashimoto H. Matsuda T. Baba A. 2000 Endothelins increase tyrosine phosphorylation of astrocytic focal adhesion kinase and paxillin accompanied by their association with cytoskeletal components. Neuroscience101 219 227 - 35.
Lambert W. S. Clark A. F. Wordinger R. J. 2004 Effect of exogenous neurotrophins on Trk receptor phosphorylation, cell proliferation, and neurotrophin secretion by cells isolated from the human lamina cribrosa. Mol. Vis.10 289 296 - 36.
Leske M. C. Wu S. Y. Hennis A. Honkanen R. Nemesure B. 2008 Risk Factors for Incident Open-angle Glaucoma The Barbados Eye Studies. Ophthalmology115 85 93 - 37.
Liang Y. J. Lai L. P. Wang B. W. Juang S. J. Chang C. M. Leu J. G. Shyu K. G. 2006 Mechanical stress enhances serotonin 2B receptor modulating brain natriuretic peptide through nuclear factor-kappaB in cardiomyocytes. Cardiovasc. Res.72 303 312 - 38.
Lin S. M. Du P. Huber W. Kibbe W. A. 2008 Model-based variance-stabilizing transformation for Illumina microarray data. Nucleic Acids Res. 36,e11 EOF e11 EOF - 39.
Liton P. B. Liu X. Challa P. Epstein D. L. Gonzalez P. 2005 Induction of TGF-beta1 in the trabecular meshwork under cyclic mechanical stress. J. Cell Physiol205 364 371 - 40.
Liu B. Chen H. Johns T. G. Neufeld A. H. 2006 Epidermal growth factor receptor activation: an upstream signal for transition of quiescent astrocytes into reactive astrocytes after neural injury. J. Neurosci.26 7532 7540 - 41.
Lukas T. J. Miao H. Chen L. Riordan S. Li W. Crabb A. M. Wise A. Du P. Lin S. M. Hernandez M. R. 2008 Susceptibility to glaucoma: differential comparison of the astrocyte transcriptome from glaucomatous African American and Caucasian American donors. Genome Biol. 9,R111 EOF - 42.
Lukas T. J. Wang A. L. Yuan M. Neufeld A. H. 2009 Early cellular signaling responses to axonal injury. Cell Commun. Signal. 7,5 EOF - 43.
Malone P. Miao H. Parker A. Juarez S. Hernandez M. R. 2007 Pressure induces loss of gap junction communication and redistribution of connexin 43 in astrocytes. Glia55 1085 1098 - 44.
Mc Graw J. Hiebert G. W. Steeves J. D. 2001 Modulating astrogliosis after neurotrauma. J. Neurosci. Res.63 109 115 - 45.
Menniti F. S. Faraci W. S. Schmidt C. J. 2006 Phosphodiesterases in the CNS: targets for drug development. Nat. Rev. Drug Discov.5 660 670 - 46.
Miao H. Chen L. Riordan S. M. Li W. Juarez S. Crabb A. M. Lukas T. J. Du P. Lin S. M. Wise A. Agapova O. A. Yang P. Gu C. C. Hernandez M. R. 2008 Gene expression and functional studies of the optic nerve head astrocyte transcriptome from normal african americans and caucasian americans donors. PLoS. ONE. 3,e2847 EOF - 47.
Miao H. Crabb A. W. Hernandez M. R. Lukas T. J. 2010 Modulation of factors affecting optic nerve head astrocyte migration. Invest Ophthalmol. Vis. Sci.51 4096 4103 - 48.
Neary J. T. Kang Y. Willoughby K. A. Ellis E. F. 2003 Activation of extracellular signal-regulated kinase by stretch-induced injury in astrocytes involves extracellular ATP and2 purinergic receptors. J. Neurosci. 23, 2348-2356. - 49.
Nemesure B. He Q. Mendell N. Wu S. Y. Hejtmancik J. F. Hennis A. Leske M. C. 2001 Inheritance of open-angle glaucoma in the Barbados family study. Am. J. Med. Genet.103 36 43 - 50.
Ni W. Rajkumar K. Nagy J. I. Murphy L. J. 1997 Impaired brain development and reduced astrocyte response to injury in transgenic mice expressing IGF binding protein-1. Brain Res.769 97 107 - 51.
Nickells,R.W. 2007 From ocular hypertension to ganglion cell death: a theoretical sequence of events leading to glaucoma. Can. J. Ophthalmol.42 278 287 - 52.
Nikolsky Y. Ekins S. Nikolskaya T. Bugrim A. 2005 A novel method for generation of signature networks as biomarkers from complex high throughput data. Toxicol. Lett.158 20 29 - 53.
Oki M. Aihara H. Ito T. 2007 Role of histone phosphorylation in chromatin dynamics and its implications in diseases. Subcell. Biochem.41 319 336 - 54.
Pena J. D. Netland P. A. Vidal I. Dorr D. A. Rasky A. Hernandez M. R. 1998 Elastosis of the lamina cribrosa in glaucomatous optic neuropathy. Exp. Eye Res.67 517 524 - 55.
Pena J. D. Taylor A. W. Ricard C. S. Vidal I. Hernandez M. R. 1999a Transforming growth factor beta isoforms in human optic nerve heads. Br. J. Ophthalmol.83 209 218 - 56.
Pena J. D. Varela H. J. Ricard C. S. Hernandez M. R. 1999b Enhanced tenascin expression associated with reactive astrocytes in human optic nerve heads with primary open angle glaucoma. Exp. Eye Res.68 29 40 - 57.
Quigley H. A. Varma R. Tielsch J. M. Katz J. Sommer A. Gilbert D. L. 1999 The relationship between optic disc area and open-angle glaucoma: the Baltimore Eye Survey. J. Glaucoma.8 347 352 - 58.
Quigley H. A. Vitale S. 1997 Models of open-angle glaucoma prevalence and incidence in the United States. Invest Ophthalmol. Vis. Sci.38 83 91 - 59.
Ricard C. S. Kobayashi S. Pena J. D. Salvador-Silva M. Agapova O. Hernandez M. R. 2000 Selective expression of neural cell adhesion molecule (NCAM)-180 in optic nerve head astrocytes exposed to elevated hydrostatic pressure in vitro. Brain Res. Mol. Brain Res.81 62 79 - 60.
Ricard C. S. Pena J. D. Hernandez M. R. 1999 Differential expression of neural cell adhesion molecule isoforms in normal and glaucomatous human optic nerve heads. Brain Res. Mol. Brain Res.74 69 82 - 61.
Salvador-Silva M. Aoi S. Parker A. Yang P. Pecen P. Hernandez M. R. 2004 Responses and signaling pathways in human optic nerve head astrocytes exposed to hydrostatic pressure in vitro. Glia45 364 377 - 62.
Salvador-Silva M. Ricard C. S. Agapova O. A. Yang P. Hernandez M. R. 2001 Expression of small heat shock proteins and intermediate filaments in the human optic nerve head astrocytes exposed to elevated hydrostatic pressure in vitro. J. Neurosci. Res.66 59 73 - 63.
Sanden N. Thorlin T. Blomstrand F. Persson P. A. Hansson E. 2000 Hydroxytryptamine2B receptors stimulate Ca2+ increases in cultured astrocytes from three different brain regions. Neurochem. Int.36 427 434 - 64.
Sawada Y. Tamada M. Dubin-Thaler B. J. Cherniavskaya O. Sakai R. Tanaka S. Sheetz M. P. 2006 Force sensing by mechanical extension of the Src family kinase substrate130Cas Cell , 1015-1026. - 65.
Smyth,G.K. 2004 Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat. Appl. Genet. Mol. Biol. 3, Article3. - 66.
Sofroniew,M.V. 2005 Reactive astrocytes in neural repair and protection. Neuroscientist.11 400 407 - 67.
Stark C. Breitkreutz B. J. Reguly T. Boucher L. Breitkreutz A. Tyers M. 2006 BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 34, D535 D539. - 68.
Varma R. Tielsch J. M. Quigley H. A. Hilton S. C. Katz J. Spaeth G. L. Sommer A. 1994 Race-, age-, gender-, and refractive error-related differences in the normal optic disc. Arch. Ophthalmol.112 1068 1076 - 69.
Wang L. Fortune B. Cull G. Dong J. Cioffi G. A. 2006 Endothelin B receptor in human glaucoma and experimentally induced optic nerve damage. Arch. Ophthalmol.124 717 724 - 70.
Wordinger R. J. Agarwal R. Talati M. Fuller J. Lambert W. Clark A. F. 2002 Expression of bone morphogenetic proteins (BMP), BMP receptors, and BMP associated proteins in human trabecular meshwork and optic nerve head cells and tissues. Mol. Vis.8 241 250 - 71.
Yang J. L. Neufeld A. H. Zorn M. B. Hernandez M. R. 1993 Collagen type I mRNA levels in cultured human lamina cribrosa cells: effects of elevated hydrostatic pressure. Exp. Eye Res.56 567 574 - 72.
Yang P. Agapova O. Parker A. Shannon W. Pecen P. Duncan J. Salvador-Silva M. Hernandez M. R. 2004 DNA microarray analysis of gene expression in human optic nerve head astrocytes in response to hydrostatic pressure. Physiol Genomics17 157 169 - 73.
Yang Z. Quigley H. A. Pease M. E. Yang Y. Qian J. Valenta D. Zack D. J. 2007 Changes in gene expression in experimental glaucoma and optic nerve transection: the equilibrium between protective and detrimental mechanisms. Invest Ophthalmol. Vis. Sci.48 5539 5548 - 74.
Zeeberg B. R. Feng W. Wang G. Wang M. D. Fojo A. T. Sunshine M. Narasimhan S. Kane D. W. Reinhold W. C. Lababidi S. Bussey K. J. Riss J. Barrett J. C. Weinstein J. N. 2003 GoMiner: a resource for biological interpretation of genomic and proteomic data. Genome Biol. 4,R28 EOF - 75.
Zode G. S. Clark A. F. Wordinger R. J. 2007 Activation of the BMP canonical signaling pathway in human optic nerve head tissue and isolated optic nerve head astrocytes and lamina cribrosa cells. Invest Ophthalmol. Vis. Sci.48 5058 5067