Susceptibility to glaucoma: differential comparison of the astrocyte transcriptome from glaucomatous African American and Caucasian American donors
- Thomas J Lukas†1Email author,
- Haixi Miao†2,
- Lin Chen2,
- Sean M Riordan2,
- Wenjun Li2,
- Andrea M Crabb2,
- Alexandria Wise3,
- Pan Du4,
- Simon M Lin4 and
- M Rosario Hernandez2
© Lukas et al.; licensee BioMed Central Ltd. 2008
Received: 9 May 2008
Accepted: 9 July 2008
Published: 09 July 2008
Epidemiological and genetic studies indicate that ethnic/genetic background plays an important role in susceptibility to primary open angle glaucoma (POAG). POAG is more prevalent among the African-descent population compared to the Caucasian population. Damage in POAG occurs at the level of the optic nerve head (ONH) and is mediated by astrocytes. Here we investigated differences in gene expression in primary cultures of ONH astrocytes obtained from age-matched normal and glaucomatous donors of Caucasian American (CA) and African American (AA) populations using oligonucleotide microarrays.
Gene expression data were obtained from cultured astrocytes representing 12 normal CA and 12 normal AA eyes, 6 AA eyes with POAG and 8 CA eyes with POAG. Data were normalized and significant differential gene expression levels detected by using empirical Bayesian shrinkage moderated t-statistics. Gene Ontology analysis and networks of interacting proteins were constructed using the BioGRID database. Network maps included regulation of myosin, actin, and protein trafficking. Real-time RT-PCR, western blots, ELISA, and functional assays validated genes in the networks.
Cultured AA and CA glaucomatous astrocytes retain differential expression of genes that promote cell motility and migration, regulate cell adhesion, and are associated with structural tissue changes that collectively contribute to neural degeneration. Key upregulated genes include those encoding myosin light chain kinase (MYLK), transforming growth factor-β receptor 2 (TGFBR2), rho-family GTPase-2 (RAC2), and versican (VCAN). These genes along with other differentially expressed components of integrated networks may reflect functional susceptibility to chronic elevated intraocular pressure that is enhanced in the optic nerve head of African Americans.
Glaucoma comprises a group of diseases that are characterized by optic neuropathy associated with optic disc cupping and loss of visual field and, in many patients, with elevated intraocular pressure (IOP) . There are several types of glaucoma, including juvenile and adult-onset types, primary open angle glaucoma (POAG), narrow-angle glaucoma, and secondary glaucoma, with different pathogenic mechanisms. POAG is more prevalent in Black Americans of African American (AA) ancestry than in Caucasian American (CA) populations 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 . The disease is particularly frequent in Afro-Caribbean persons, with a prevalence of 7% in Barbados and 8.8% in St Lucia . On average, African Americans have the longest duration  and higher progression of disease  compared to other populations. In addition to racial differences, a positive family history of POAG is a major risk factor for the disease in African Americans . 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 [7, 8].
Abnormally elevated IOP elicits a complex sequence of putative neurodestructive and neuroprotective cellular responses in the optic nerve head (ONH) . Previous studies demonstrated that gene expression in astrocytes of the glaucomatous ONH serve as the basis for these responses . Here we present evidence that primary cultures of AA and CA astrocytes derived from POAG donors exhibit differential gene expression of genes that relate to reactive astrocytes and to pathological changes that occur in the glaucomatous ONH. Validations of changes in expression of selected genes were done by quantitative real-time RT-PCR, western blots, enzyme-linked immunosorbent assay (ELISA) and various functional assays. Network analysis of gene product interactions focused our findings on specific functional pathways. Our data indicate that both normal and glaucomatous astrocytes from AA donors exhibit differential expression in genes that regulate signal transduction, cell migration, intracellular trafficking and secretory pathways.
Results and discussion
Primary cultures of ONH astrocytes from normal and glaucomatous donors
Demographics and clinical history
Demographic characteristics of the normal AA and CA donors used in this study are detailed in Additional data file 2. Demographic and clinical data for AA donors with glaucoma (AAGs) and CA donors with glaucoma (CAGs) included in the microarray analyses and other assays are detailed in Additional data file 1. Twelve eyes from ten CAG donors and six eyes from AAG donors were used in this study. Glaucoma drug treatment history was available for some POAG donors. None of the drug treatments are known to affect astrocytes in the ONH. The degree of glaucomatous damage in donors with POAG was assessed using histories when available and by evaluating axon degeneration in cross-sections of the myelinated optic nerve (Additional data file 1). A limitation of this study is that only six eyes from three AAG donors were available due to the extreme rarity of these samples. Consequently, we used all six eyes to generate primary cultures for all experiments in our study. Primary cultures of samples from AAG and CAG donors were fully characterized as ONH astrocytes as described in detail earlier .
Identification of differentially expressed genes in ONH astrocytes from AA and CA donors with POAG
For the comparisons amongst the four groups, our primary focus was to establish the differentially expressed genes between AAG and CAG donors (Additional data file 7); our secondary focus was the comparison between normal and glaucomatous astrocytes and our tertiary focus was to identify differentially expressed genes within each population: AAG versus AA and CAG versus CA.
Common genes significantly decreased in glaucomatous ONH astrocytes compared to their normal counterparts
Adhesion molecule with Ig-like domain 2
Bone morphogenetic protein 1
Cysteine-rich protein 2
Diacylglycerol kinase, alpha 80 kDa
Dystrophia myotonica-protein kinase
EF-hand domain family, member D1
Microtubule associated monoxygenase, calponin and LIM domain containing 2
NIPA-like domain containing 3
Platelet-derived growth factor alpha polypeptide
Solute carrier family 12, member 2
Solute carrier family 12, member 4
WW domain containing E3 ubiquitin protein ligase 2
Differentially expressed genes in glaucomatous astrocytes*
Genes associated with myosin regulation
Myosin, heavy chain 2
Myosin, light polypeptide kinase
Phosphoinositide-3-kinase, subunit (p85-alpha)
Protein phosphatase 1, regulator subunit 12A (PPP1R12A)
Ras-related 2 (Rho family, Rac2)
Ribosomal protein S6 kinase, 90 kDa, polypeptide 3
Genes associated with actin regulation
Rho guanine nucleotide exchange factor (GEF) 7
NCK adaptor protein 1
PDZ and LIM domain 1 (elfin, CLP36)
Phosphoinositide-3-kinase, regulatory subunit 1
Plectin 1, intermediate filament binding protein
Protein tyrosine phosphatase, non-receptor type 11
Ras-related 2 (Rho family, Rac2)
SMAD, mothers against DPP homolog 3
Transforming growth factor, beta receptor I
Transforming growth factor, beta receptor II
Genes associated with protein trafficking
Amyloid beta precursor protein binding protein 1
Chemokine (C-C motif) ligand 5
Cadherin 2, type 1, N-cadherin (neuronal)
Collagen, type IV, alpha 4
Catenin (cadherin-associated protein), beta 1, 88 kDa
Catenin (cadherin-associated protein), delta 1
Golgi autoantigen, golgin subfamily a, 1
Golgi autoantigen, golgin subfamily a, 2
Golgi autoantigen, golgin subfamily a, 3
Hyaluronan and proteoglycan link protein 1
Protease, serine, 3 (mesotrypsin)
RAB1A, member RAS oncogene family
RAB4A, member RAS oncogene family
RAB5B, member RAS oncogene family
RAB9A, member RAS oncogene family
RAB9 effector protein with kelch motifs
Rab geranylgeranyltransferase, beta subunit
Versican (chondroitin sulfate proteoglycan 2, CSPG2)
Eight eyes from six CAG donors were used to generate astrocytes for eight Hu95v2 chips. Six eyes from three AAG were used to generate astrocytes for six Hu95Av2 chips and six Hu133A 2.0 chips. Eighteen Hu133 2.0 chips from nine normal AA and nine normal CA donors, and seven Hu95v2 chips from six normal CA donors were used for comparisons within the appropriate platform. All microarray data have been deposited in the NCBI GEO database under the series accession number GSE9963.
The data measured by the two types of chips were normalized separately by RMA normalization as described in Materials and methods. Differentially expressed genes required an up or down fold-change of more than 1.5-fold (p < 0.01, false discovery rate < 0.05). A total of 618 genes were differentially expressed in AAG-CAG comparisons, 484 upregulated and 134 downregulated (Additional data file 7); 509 genes were differentially expressed in AAG compared to normal AA astrocytes, 167 upregulated and 342 downregulated (Additional data file 5); and 195 genes were differentially expressed in the CAG-CA comparison, 132 upregulated and 63 downregulated (Additional data file 6). We used empirical Bayesian methods to identify differentially expressed genes; both our results (not shown) and previous studies [12, 13] have suggested that the empirical Bayesian method has performance similar to statistical analysis of microarrays (SAM). To reduce batch effects, we added fold-change criteria because genes with larger fold-change are less likely to be affected by such effects.
Gene Ontology (GO) analysis of differential expression in glaucomatous astrocytes was done with GoMiner . There were 33 significant categories for CAG-CA, 80 for AAG-AA, and 67 for AAG-CAG comparisons (p < 0.01). The significant genes in selected categories were mined using GOstats in Bioconductor (Additional data file 9). The phosphorylation category (GoID: 16310) was significant in the three datasets. The percent distribution of the genes common to all of the datasets in this category was determined (Additional data file 10). For example, the genes encoding myosin light chain kinase (MYLK) and calcium/calmodulin-dependent serine protein kinase (CASK1) were found in all three glaucoma comparisons. Those encoding the regulatory subunit of phosphatidylinositol-3-kinase (PIK3R1), transforming growth factor (TGF)β-receptor 2 (TGFBR2), ERBB2, and Ephrin receptor A5 were some of the genes found in two datasets (AAG-CAG and AAG-AA). Similarly, another category with overlaps between the datasets was cell-cell signaling (Additional data file 10). Some of the genes in this category include those encoding latent transforming growth factor beta binding protein 4 (LTBP4), the glutamate receptor subunit (GRIK2), and parathyroid hormone-like protein (PTHLH). As we show below, expansion of these and other GO categories using network-protein interaction software yielded three networks that include differentially expressed GTPases, protein kinases, transmembrane receptors, and proteins involved in trafficking at cellular membranes. Altogether, the GO analysis suggests that alterations in the signaling networks that regulate cell motility, polarity, adhesion, and trafficking are present in glaucomatous astrocytes. Moreover, the overlap among the datasets in multiple categories suggests that there is a spectrum of changes in gene expression in glaucoma.
Common genes significantly increased in glaucomatous ONH astrocytes compared to their normal counterparts
ATP-binding cassette, sub-family A, member 8
Chromosome 5 open reading frame 30
Calcium/calmodulin-dependent serine protein kinase
Caspase 4, apoptosis-related cysteine peptidase
Glutathione S-transferase A4
GULP, engulfment adaptor PTB domain containing 1
Potassium channel, subfamily K, member 2
LIM domain only 4
Myosin, heavy polypeptide 10, non-muscle
Phosphorylase, glycogen; liver
Retinol binding protein 1, cellular
Serpin peptidase inhibitor, clade G, member 1
SH3-domain binding protein 5
Slit homolog 2
TGF beta-inducible nuclear protein 1
Cellular motility and migration in AAG astrocytes
Migration of reactive astrocytes is an important component in the remodeling of the ONH in glaucoma [15, 16]. In glaucoma, reactive astrocytes migrate from the cribriform plates into the nerve bundles [9, 17] and synthesize neurotoxic mediators such as nitric oxide and tumor necrosis factor (TNF)α, which may be released near the axons, causing neuronal damage [18, 19]. Previous work in our laboratory demonstrated that human ONH astrocytes in vitro respond to elevated pressure predominantly with an increase in cell migration that may be relevant to axonal degeneration and tissue remodeling in glaucomatous optic neuropathy .
Here we provide in vitro data of differential astrocyte migration in astrocytes from AAG donors using a standardized migration assay. As shown in Figure 1a, migration of AAG astrocytes is significantly increased compared to CAG astrocytes and migration is faster in AA compared to CA astrocytes. Because multiple cellular processes impact cell motility and migration, we divided our analysis between two interacting networks that regulate myosin and actin.
Myosin-dependent astrocyte migration
From the microarray and quantitative RT-PCR (qRT-PCR) data, the following genes related to myosin regulation were differentially expressed in AAG: MYLK, MYPT1, RAC2, CALM1, RPS6KA3, MYH10, and PIK3R1. Shown in Figure 1b is the network of proteins associated with the phosphorylation of the regulatory light chain of myosin II and activation of myosin-ATPase (MYH10). Two network nodes are critical for the regulation of myosin. These include MYLK, a calmodulin-activated protein kinase that phosphorylates Ser19 on the myosin regulatory light chain and MYPT1, the regulatory subunit of myosin-light chain phosphatase, which dephosphorylates the myosin light chain. We found that both genes were expressed in AAG astrocytes at significantly higher levels than in CAG astrocytes (Table 3). Similarly, calmodulin (CALM1), the activator of MYLK is also upregulated in AAG astrocytes (Table 3)
The upregulation of MYLK suggests that the myosin regulatory system may exhibit increased responsiveness towards modulation by various cellular second messenger signaling systems such as Ca2+, diacylglycerol, and cyclic nucleotides . Similarly, changes in expression of RAC2 indicate that other members of the Rho-family signaling network are altered in AAG astrocytes (Figure 1c). These changes allow us to predict that the myosin-regulated motility may be sensitized to signals from Ca2+, Rho GTPase, and growth/trophic factors coupled to the activation of phosphoinositides. Within the phosphoinositide pathway, PIK3R1 is upregulated in AAG astrocytes (Figure 1c). The PIK3R1 pathway is important for the motility of ONH astrocytes  and their responses to increased hydrostatic pressure . PIK3R1 is the regulatory subunit of the lipid kinase that transforms phosphoinositide (4,5) biphosphate (PIP2) into the triphosphate (PIP3). PIP3 in turn mediates activation of several of the Rho GTPases as well as selected protein kinases. Thus, in AAG astrocytes, lipid-activated pathways that modulate astrocyte motility are altered.
ERK1 potentiates MYLK activity through phosphorylation  and interacts with PEA15 (Phosphoprotein enriched in astrocytes) . The increased expression of the S6-family kinase (RPS6KA3) may compete with ERK1 for binding to the phosphoprotein PEA15 , potentially increasing the pool of active ERK1. Consistent with this finding, we have shown that ERK1 is activated in normal CA ONH astrocytes, under increased hydrostatic pressure and in experimental glaucoma in primates . Thus, myosin-based motility may be influenced by changes in MYLK expression and potentiation through ERK1 activation under hydrostatic pressure.
The MYLK gene has multiple genes within its locus . In some tissues up to three transcripts are expressed, including for long and short forms of the kinase and a protein identical to the carboxyl-terminal sequence . ONH astrocytes express both the 130 kDa (MYLK-130) and 210 kDa (MYLK-210) kinase isoforms and we quantified changes in both using standard densitometry measurements. Western blots (Figure 4b) show that the fraction of MYLK-210 in ONH astrocytes is higher in AAG and CAG compared to normal astrocytes, while the fraction of the MYLK-130 isoform decreases (Figure 4b). These differences were quantified using densitometry (Figure 4c, d). Thus, in glaucoma there appears to be MYLK isoform switching towards the larger protein. The difference between the two proteins is the presence of an amino-terminal extension in the 210 kDa species that contains additional actin binding domains. Other studies have shown that MYLK-210 displays enhanced interaction with the actin cytoskeleton compared to the 130 kDa isoform [28, 29]. These results are consistent with the enhanced migration of ONH astrocytes mediated in part by increased expression of MYLK-210.
MYLK variants have been found to confer risk of lung injury , asthma or sepsis , particularly in African Americans . Some of the common polymorphisms in MYLK affect its expression . Therefore, in some populations, it is possible that the effects of increased expression of MYLK may be further modified by genetic polymorphisms.
Actin-dependent astrocyte migration
From the microarray and qRT-PCR data the following genes were differentially expressed in AAG: TGFBR2, TGFBR1, SMAD3, NCK1, PTPN11, ARHGEF7, PDLIM1, LM04, and PLEC1. Figure 2a shows several signal transduction networks that participate in the regulation of actin. Remodeling or redistribution of actin at cellular edges is an essential part of establishing cell polarity  and the formation of processes in astrocytes . Actomyosin interactions and actin polymerization are regulated by intracellular proteins such as α-actinin (ACTN4) and the ARP protein complex (ACTR2, WASP: Figure 2a). These networks involve the Rho GTPase signaling pathway. Therefore, we used a pull-down Rho activation assay to measure activated Rho in cell lysates. ONH astrocytes from CAG and AAG donors exhibited significantly higher Rho activity compared to those from normal AA and CA donors (Figure 2b, c), consistent with the differential expression of Rho regulatory components. Rho activity was also increased in astrocytes exposed to elevated hydrostatic pressure . Thus, increased Rho activity is another contributor towards increased migration of AAG astrocytes. We suspect that Rho activity may be altered by changes in the signaling proteins included in these networks. For example, RAC2 and ARGEF7 are upregulated in AAG. The Rho-family GTPase, RAC2, is downstream of TGFβ signaling  and ARHGEF7 stimulates guanine nucleotide exchange on Rho family GTP-binding proteins. We further elaborated changes in TGFβ signaling as a driver to changes in Rho activity.
TGFβ signaling in AAG astrocytes
TGFβ regulates cellular motility through two components. One is through the expression of extracellular matrix (ECM) proteins, which will be discussed in detail below. Contractile forces are transmitted to the ECM through actin-based stress fibers via focal adhesions, which are assemblies of ECM proteins, transmembrane receptors, and cytoplasmic structural and signaling proteins, such as integrins. TGFβ modulates integrin-mediated cellular migration, where FYN is one of the primary signal transducing proteins. A second component of TGFβ signaling is the regulation of cell polarity. For example, PARD3 and PARD6 are part of a multi-component polarity complex that controls polarized cell migration . These complexes involve the Rho, CDC42, and RAC signaling pathways, which provide the means to remodel actin during migration [33, 41]
As shown in Figure 2a, NCK1 was upregulated in AAG (Table 3). The Nck1 SH2/SH3 adaptor couples phosphotyrosine signals to the actin cytoskeleton and receptor signaling to the regulatory machinery of the cytoskeleton . The enigma family member PDLIM1 was upregulated in AAG astrocytes (Table 3) and functions by allowing interactions among cytoskeletal proteins through PDZ and amino LIM domains [43, 44]. Downregulation of other actin binding proteins such as PLEC1 (Table 3) may alter actin dynamics with respect to cytoskeletal changes induced by Rho-GTPase, phospholipids, and tyrosine kinase (Src) mediated signaling .
TGFBR2 receptors in optic nerve head astrocytes
Figure 5b illustrates immunohistochemistry of the TGFBR2 on astrocytes in normal and glaucomatous ONH tissue. GFAP positive astrocytes in the lamina cribrosa of AAG exhibit higher expression of TGFBR2 compared with astrocytes in normal ONH tissue. Consistent with these findings, western blots of lysates of ONH astrocytes from AAG indicate higher levels of TGFBR2 protein compared to the normal tissue and CAG (Figure 5c).
To further investigate alterations in TGFβ signaling in ONH astrocytes, we examined the production of TGFβ1 and TGFβ2. As seen in Figure 5d, TGFβ2 is the primary form of TGFβ produced by ONH astrocytes . There are significantly increased levels of secreted TGFβ1 in AA compared to CA astrocyte supernatants but the increases in AAG and CAG astrocytes were not significant compared to normal astrocytes. These data suggest that most of the changes in TGFβ signaling are due to alterations at the level of TGFβ receptors in astrocytes from AAG.
Mutations in TGFBR2 are associated with Marfan syndrome type 2 [47–49]. Ocular abnormalities, including glaucoma, are associated with Marfan syndrome type 1 in which there are mutations in the gene for fibrillin (FBN1) . However, it has not been established that mutations of TGFBR2 are associated with ocular problems in Marfan syndrome type 2 [48, 49].
Intracellular trafficking and the endoplasmic reticulum/Golgi compartments
From the microarray and quantitative RT-PCR data the following genes were differentially expressed in AAG. Endosome group, RAB4A, RAB5B, RAB9P40, RAB9A; plasma membrane group, PRSS3, APPB1, CTNND1, CTNNB1, CDH2, VCAN, HAPLN1, CCL5, COL4A4, TGM2, SLIT2, GPC1; Golgi group, GOLGA1, GOLGA3, GOLGA2, RAB1A, RABGGTB (Figure 3a). Six Rab family signaling genes involved in intracellular transport of organelles were differentially regulated (Table 3). Three small GTPases, RAB4A, RAB5B, and RAB9A, were upregulated (Table 3, Figure 3b), suggesting increased endosomal transport and processing. RAB4A and RAB5B selectively regulate intracellular trafficking and signaling of G protein-coupled receptors, such as the angiotensin receptor and adrenergic receptors (β2-AR and α2B-AR) from the cell surface [51, 52]. RAB9A participates in late endosomal events leading to fusion with the lysosomal compartment .
In AAG astrocytes there was a predominant increase in transcription of Golgi-resident protein transcripts (Additional data file 7). These include RAB1A, and three members of the golgin family, GOLGA1, GOLGA2 and GOLGA3 (Table 3), which may function in the stacking of Golgi cisternae and in vesicular transport . GOLGA3 promotes cell surface expression of the beta adrenergic receptors . Thus, the increased expression of Golgi proteins may further enhance adrenergic receptor signaling. Note that the RAB proteins upregulated in the endosomal pathway (above) also affect trafficking of these receptors.
Included in the protein trafficking network are plasma membrane associated proteins involved in cell-cell communication from the junctional matrix (Figure 3a). Catenins (CTNNB1, CTNND1) form membrane trafficking complexes that integrate other cadherins (CDH2), and members of the amyloid precursor protein complex (presenilin, APPBP1, PRSS3). In particular, CTNND1 functions to regulate membrane trafficking either through blocking cadherin interactions, or through Rho-GTPases such as Rho A, Rac and CDC42 . As with the myosin and actin motility networks, the change in expression of GTPase regulatory proteins will likely impact plasma membrane trafficking. The upregulation of chondroitin sulfate proteoglycan 2 (versican; VCAN), transglutaminase 2 (TGM2), and hyaluronan and proteoglycan link protein 1 (HAPLN1) are significant modifiers of the ECM . Both HAPLN1 and VCAN mRNA levels were upregulated in AAG compared to CAG astrocytes by qRT-PCR (Figure 3b). VCAN immunoreactivity was observed in the ECM of the cribriform plates, the perivascular matrix and a few astrocytes in the lamina cribrosa of normal AA and CA donors (Figure 3c). In glaucomatous tissues there was a marked increase in VCAN staining in astrocytes in the cribriform plates and hypertrophied reactive astrocytes in the nerve bundles in both populations (Figure 3c). TGFβ2 signaling upregulates VCAN  in astrocyte cell types and expression of collagen type 4 and transglutaminase 2 in ONH astrocytes . Our data on changes in TGFβ receptor expression and ECM proteins are similar to those found in microarray profiling of ONH tissue from a rat model of glaucoma . Expression of ECM proteins is also modulated by TGFβ in GFAP-negative lamina cribrosa cells in culture .
There is substantial evidence that ONH astrocytes are responsible for the normal maintenance of the ECM in normal tissue and that reactive astrocytes remodel the ECM in response to elevated IOP in human and experimental glaucoma [10, 60, 61]. Reactive astrocytes in the ONH express abnormal ECM in glaucoma, leading to loss of resiliency and deformability in response to elevated IOP. Alterations in TGFβ2 levels and TGFβ receptors and abnormal synthesis of ECM in AAG may convey connective tissue components of susceptibility to elevated IOP to this population.
cAMP signaling in glaucomatous ONH astrocytes
We searched the expression data for differentially expressed genes that might explain the difference in basal cAMP levels between glaucomatous and normal astrocytes. One potential candidate for increasing basal cAMP is PTHLH, a parathyroid hormone-like protein that is upregulated in glaucomatous astrocytes (Figure 6b). This protein binds to ubiquitous PTH receptors that are coupled to stimulation of adenylate cyclase and elevated cyclic AMP . Thus, upregulation of PTHLH provides an autocrine pathway leading to increased basal cyclic AMP levels in glaucomatous astrocytes. Another gene that might also contribute to the activity of adenylate cyclases is CAP2 . However, we found that CAP2 was not differentially expressed in glaucomatous ONH astrocytes by qRT-PCR (Figure 6b).
Other disease-associated genes differentially regulated in glaucomatous OHN astrocytes
In this list we found that specific cell-surface-associated proteins are downregulated in glaucoma. These include BMP1, which activates cleavage of LTBP1 proteins that release nascent TGFβ1 , and AMIGO2, a type I transmembrane protein that regulates axon extension . Down-regulation of BMP1 may reduce the levels of free TGFβ1 and thus unbalance signaling between TGFβ isoforms. A decrease in AMIGO2 might negatively impact axon survival.
Two differentially expressed genes that are involved in reactive astrocyte responses to neuronal injury are SLIT2 and GPC1 (glypican-1). SLIT2 serves as a chemorepellant for multiple types of axons , while GPC1 is a proteoglycan that binds SLIT2 . Upregulation of expression of SLIT2 and a reduction of GPC1 by glaucomatous astrocytes suggest an inhibitory microenvironment for RGC axons in the ONH. These data are consistent with the idea that the enhanced migratory properties of glaucomatous astrocytes coupled with the release of factors that negatively impact upon axon survival are part of the pathophysiology of the disease.
Finally, lectin-like oxidized-LDL receptor (OLR1; also known as LOX-1) is highly upregulated (Additional data files 5 and 11) in AAG astrocytes. OLR1 expression is induced by TGFβ1 signaling and is known to be a component of the fluid shear stress response of endothelial cells in early atherosclerotic lesions . These data are further confirmation of enhanced TGFβ signaling in AAG astrocytes as suggested by the differential receptor expression described earlier.
Intracellular calcium signaling/transport systems in ONH astrocytes
Two genes directly involved in Ca2+ homeostasis are differentially regulated in ONH astrocytes of AAG (Additional data files 5 and 7). CACNB4 encodes a beta subunit of the voltage-dependent calcium channel complex. CACNB4 plays an important role in calcium channel function by modulating G protein inhibition, increasing peak calcium current, controlling the alpha-1 subunit targeting to the membrane and shifting the voltage dependence of activation and inactivation. The second gene, ATP2C1 (Additional data file 7), encodes a protein that belongs to the family of P-type primary ion transport ATPases, which pump Ca2+ into the endoplasmic reticulum.
Transcripts encoding the calcium/calmodulin-related signaling proteins calmodulin 1 (CALM1) and Ca2+/calmodulin-dependent membrane-associated kinase (CASK1) are differentially expressed in one or more glaucoma groups. CALM1 was increased in AAG compared to AA donors (Additional data file 5), while CASK1 was increased in glaucomatous astrocytes from both AA and CA donors (Table 1, Figure 7). Calmodulin is the Ca2+ sensor of key signaling molecules, such as adenylyl cyclase, CAMKII, CAMKIV, and MYLK discussed above. CASK1 is a member of the membrane-associated guanylate kinase proteins (MAGUKs), a prominent family of scaffolding molecules associated with intercellular junctions. CASK1 targets Ca2+ and K+ channels  and/or the Ca2+ pump 4b/CI  to the plasma membrane, interacts with liprins  and regulates transcription by interacting with transcription factors in the nucleus . Interestingly, CASK is a candidate gene for X-linked optic atrophy . The differential expression of genes in Ca2+ signaling pathways could be a common theme in glaucomatous astrocytes that may have a higher impact in optic nerves from AA donors due to increased sensitivity to elevated IOP in these donors.
Glaucomatous ONH astrocytes share many characteristics of reactive astrocytes in the central nervous system; however, certain properties may be specific to the pathophysiology of glaucoma. The current work and previous studies demonstrate that cultured glaucomatous ONH astrocytes exhibit differential expression of genes that promote cell motility and migration, downregulate cell adhesion, are associated with structural tissue changes, and contribute to neural degeneration. Our data further strengthen the idea that reprogramming of transcription in glaucomatous astrocytes shifts signaling towards TGFβ, Rho GTPase and Ca2+ systems, which impact the multiple networks described earlier.
Our demonstration of this wide variety of genes that remain differentially expressed after weeks in culture suggests that glaucomatous ONH astrocytes have an altered phenotype. In the current study, using microarray analysis, we identified a number of genes (for example, MYLK, TGFBR2, VCAN, and RAC2) whose expression may underlie higher susceptibility of astrocytes of AA individuals to elevated IOP and that may be relevant to reactive astrocyte responses in glaucoma. Some limitations of our approach should be noted. First, ONH astrocyte derived from human glaucomatous eyes during the disease process does not allow assessment of changes or the identification of early mechanisms of disease that might be available from animal models. In addition, the difficulty to obtain and include more AA glaucomatous eyes limited our ability to identify differentially expressed genes in this group. However, stringent filters allowed the selection of a group of genes with functional significance. For each comparison, selected genes were validated by qRT-PCR and relevant gene products were confirmed by western blots in the four groups.
We propose that part of the increased susceptibility to elevated IOP in AAG relates directly to astrocyte functions in the ONH. Astrocytes in AAG, which are reactive astrocytes, may have increased responsiveness to TGFβ signaling and enhanced migratory abilities, which may impact the remodeling of the ECM, inhibit axon survival, and alter vascular permeability in the glaucomatous ONH. Any one of these changes may represent a susceptibility risk factor in the AA population to withstand abnormally elevated IOP.
This study provides an initial survey of the molecular differences of ONH astrocytes from AA and CA donors with glaucoma. Genes encoding many potential therapeutic targets, such as motility genes, ion channels, adhesion molecules, and signaling pathways, are selectively expressed in glaucomatous astrocytes, making them interesting as potential targets for astrocyte-specific therapeutics. Additional applications of these data include identification and characterization of signaling pathways involved in astrocyte function and further exploration of the role of selected identified genes in experimental animal and in vitro models of glaucoma.
Materials and methods
Thirteen eyes from eleven CA donors (age 73 ± 9 years) with POAG (referred to as CAG) and six eyes from three AA donors (age 62 ± 13 years) with POAG (referred to as AAG) were used to generate ONH astrocyte cultures as described (Additional data file 1). Myelinated optic nerves were fixed in 4% paraformaldehyde, post-fixed in osmium, embedded in epoxy resin and stained with paraphenylendiamine to detect axon degeneration as described earlier to confirm glaucoma and to assess optic nerve damage (Additional data file 1). Normal eyes were from 12 CA donors (age 60 + 11 years) and 12 AA donors (age 58 + 12 years) with no history of eye disease, diabetes, or chronic central nervous system disease (Additional data file 2).
Primary cultures of human ONH astrocytes were generated as described previously . Briefly, four explants from each lamina cribrosa were dissected and placed into 25 cm2 Primaria tissue culture flasks (Falcon, Lincoln Park, NJ, USA). Explants were maintained in DMEM-F12 supplemented with 10% fetal bovine serum (Biowhittaker, Walkerswille, MD, USA) and 10 μl/ml of PSFM (10,000 U/ml penicillin, 10,000 μg/ml streptomycin and 25 μg/ml amphotericin B; Gibco/BRL, Gaithersburg, MD, USA). Cells were kept in a 37°C, 5% CO2 incubator. Primary cultures were purified by using an immunopanning procedure . Purified cells were expanded after characterization by immunostaining for astrocyte markers GFAP and NCAM (Neural cell adhesion molecule) as described . Second passage cell cultures were stored in RPMI 1640 with 10% DMSO in liquid nitrogen until use. For each set of experiments, cells were thawed and cultured so that sufficient cells from the same batch were available for multiple experiments.
An affinity purified rabbit polyclonal antibody to MYLK was a gift from Dr Linda van Eldik (Northwestern University). It was used in western blotting (1:10,000) and immunohistochemistry (1:50). Another MYLK antibody (M7905) is a mouse monoclonal antibody (Sigma-Aldrich, St Louis, MO, USA) and it was used in western blotting (1:10,000) and immunohistochemistry(1:50). TGFBRII (L-21) is a rabbit polyclonal antibody (Santa Cruz Biotechnology Inc., Santa Cruz, CA, USA). It was used in western blotting (1: 200) and immunohistochemistry (1:50). VCAN is a goat polyclonal antibody (R&D Systems, Minneapolis, MN, USA). It was used in immunohistochemistry (1:20).
Oligonucleotide microarray analysis
Total RNA was extracted using Qiagen RNeasy mini kits (Qiagen, Valencia, CA, USA). RNA was then purified and quantified by measuring absorbance at 260 nm. Quality and intactness of the RNA was assessed by capillary electrophoresis analysis using an Agilent 2100 Bioanalyzer (Agilent, Palo Alto, CA, USA). cDNA was synthesized from 2-5 μg purified RNA by using Superscript Choice system (Gibco BRL Life Technologies, Gaithersburg, MD, USA) and T7-(dT)24 primer (GENSET, La Jolla, CA, USA). Using Bioassay High Yield RNA Transcript Labeling Kit (Enzo Diagnostics, Farmingdale, NY, USA), in vitro transcription was carried out with the cleaned double-stranded cDNA as a template in the presence of biotinylated UTP and CTP. Purified biotin-labeled cRNA was fragmented before the hybridization. Hybridization of the labeled cRNA to Human Genome U95Av2, U133A, U133A 2.0 chips (Affymetrix, Santa Clara, CA, USA) was carried out by using Genechip Instrument System (Affymetrix) at the Genechip Core Facility of Washington University School of Medicine. The arrays were washed and stained with streptavidin-phycoerythrin (Molecular Probes, Eugene, OR, USA) followed by scanning on an Agilent GeneArray Scanner G2500A (Agilent Technologies, Palo Alto, CA, USA).
Pretreatment of data
The first step in the analysis of the microarray data was to determine which genes to consider 'present' or 'absent'. We estimated the probe-set present/absent calls by using the Wilcoxon signed rank-based algorithm. In order to reduce false positives, we removed the probe-sets with all samples as 'absent' (Additional data file 3).
Comparison between glaucomatous ONH astrocytes from AA and CA normal donors
As the experiments were done at different times, two types of Affymetrix microarrays (Human Genome U95Av2 array and Human Genome U133A 2.0 array) were used. Samples from eight CAG donors, seven CA normal donors and six AAG samples were measured using a Human Genome U95Av2 array. Eighteen CA samples, eighteen AA samples and six AAG samples were measured using a Human Genome U133A 2.0 array. The data measured by two types of arrays were normalized separately by RMA normalization [74, 75]. We defined common glaucoma-related genes as genes differentially expressed in both CAG versus CA and AAG versus AA, and did comparisons of CAG versus CA and AAG versus AA separately. The differentially expressed genes were identified by the empirical Bayesian shrinkage moderated t-statistics in the limma Bioconductor package . A mixed effects model was used to account for the effect of technical replicates. Genes exhibiting a fold-change >1.5 and p-value < 0.01 were considered significant. To reduce false positives because the AAG has only three biological replicates, we applied the Benjamini and Hochberg false discovery rate multiple testing correction with a false discovery rate of 0.05 (AAG versus CAG and AAG versus AA).
To compare the significant gene list based on two types of microarray platforms, the Affymetrix probeset IDs were transferred as Entrez IDs based on the Bioconductor library. Genes whose Entrez IDs appear in both the differentially expressed gene list from CAG versus CA (using the Human Genome U95Av2 array) and AAG versus AA comparisons (using the Human Genome U133A 2.0 array) and change in the same direction were considered as common glaucoma-related genes. Genes that are differentially expressed for AAG versus CAG (using the Human Genome U95Av2 array), but without significant changes for AA versus CA (using Human Genome U133A 2.0 array), were considered as the glaucoma race-related genes. Here we considered a p-value > 0.05 as indicative of changes that were not significant.
GO analysis of differential expression in glaucomatous astrocytes was done with GoMiner . Briefly, gene lists of up- and downregulated genes (p < 0.01 as described above) were normalized to 1 and -1, respectively, for genes that exhibited at least a 1.5-fold change in either direction. These lists were then loaded into GoMiner using local GO databases accessed using the 'Derby' module. GoMiner output was analyzed with a significance cutoff of p < 0.01 and at least four genes per category.
Initially, we scanned the differentially expressed gene lists for AAG-CAG, AAG-AA, and CAG-CA comparisons for groups of genes that were either in common GO categories, or were highly over- or underexpressed (>1.5-fold, p < 0.01). These short lists were then used as a source of nodes for each network group. Networks of interacting proteins were constructed using the BioGRID database . BioGRID is a freely accessible database of physical and genetic interactions. BioGRID release version 2.0 includes more than 116,000 interactions from Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster and Homo sapiens. Graphs with embedded protein, gene and interaction attributes were constructed with a visualization program, Osprey , that is dynamically linked to the BioGRID. Each network was begun using a single gene or node. Then more interactions were added using The BioGrid Database lookup function. These were curated to simplify the graphs, and non-expanded nodes were minimized. In general, nodes that were not differentially expressed required at least two connections or edges to remain in the network. Expression of genes depicted in the networks were checked for a 'present call' in the microarray data or otherwise validated by quantitative real time RT-PCR.
Real time qRT-PCR was done as previously described . To compare expression of specific genes amongst the four groups included in this study (AA, CA, AAG and CAG), we used 12 ONH astrocyte cultures from normal CA and 12 cultures from normal AA donors. cDNA of eight eyes from eight CAG donors and of six eyes from three AAG donors were used. Individual samples were processed simultaneously under the same conditions and the data were analyzed for significance using a two-tailed t-test on sample pairs (Prism 3.0 GraphPad software, San Diego, CA, USA). Primers used in this study are listed in Additional data file 4.
Protein lysates from four samples of each group of ONH astrocytes were processed together in the appropriate combinations: four AAG and four AA; four CA and four CAG. Western blots were run in triplicate to accommodate all available samples. Protein lysates containing 3-10 μg were used depending on the specific antibody. β-Actin was used as a loading control. Films of blots were scanned using a flatbed scanner in 8-bit gray scale mode. ImageJ (National Institutes of Health) was used to quantify band intensities on the blots.
Detection of TGFβ1 and TGFβ2 by ELISA
TGFβ1 and TGFβ2 were measured in cell culture supernatants using ELISA kits (R & D Systems) specific for each protein. Briefly, astrocytes (70-80% confluent) were incubated for 24 h in 6 ml of cell culture medium without serum. Media was harvested and divided into 1 ml aliquots and frozen at -80°C until analysis. For each sample, cell counts were made and recorded. Media samples were thawed on ice and 200 μl aliquots activated by incubating with 40 μl of 0.1 N HCl at room temperature for 40 minutes. The reactions were quenched by adding 40 μl of 0.1 N NaOH in 0.5 M HEPES and mixed. Samples were diluted with the appropriate ELISA assay buffer to 400 μl. Aliquots of these solutions (50 μl TGFβ1:100 μl TGFβ2) were then assayed according to the manufacturers' protocol. Experiments were performed in duplicate and each astrocyte cell culture (n = 5-7 samples per each group) was assayed at least twice. Expressed protein values in picograms of TGFβ1/2 per ml were normalized to 106 cells using the cell counts obtained at harvest. The means of the content were considered significantly different if p < 0.05 (two-tailed t-test; Prism 3.0 GraphPad software.).
Cyclic AMP assay
Primary ONH astrocyte cultures obtained from six normal AA, six normal CA, eight CAG and three AAG were grown in 60 mm dishes until 80% confluence. Growth media was replaced with serum free media and the cells incubated for an additional 24 h. After washing with ice-cold phosphate-buffered saline (PBS), cells were lysed in 95% chilled ethanol for 1 h and then centrifuged at 2000 × g for 15 minutes at 4°C. The supernatant was evaporated using a Speed Vac concentrator and resuspended in 100 μl of the assay buffer and analyzed as described in the cAMP Biotrak Enzyme Immunoassay Kit (Amersham Bioscience RPN225, Piscataway, NJ, USA). cAMP concentration per well was expressed as pmol/mg of protein. Each value represents the mean cAMP level (± standard deviation) of independent experiments using primary astrocyte cultures from each donor and performed in triplicate. Sample pairs were analyzed by two-tailed t-test (Prism 3.0 GraphPad software) for significance (p < 0.05).
CytoSelect™ 24-well cell migration assay (Cell BioLabs, San Diego, CA USA) was used to measure the migratory properties of cells. The assay was performed according to the manufacture's protocol. Briefly, media with 10% fetal bovine serum was placed in the lower wells followed by 50,000 cells in 300 μl of serum free media in each of the well inserts. After incubation at 37°C in a 5% CO2 atmosphere for 24 h, the media was removed from the inserts. Cells that did not migrate were removed from the inserts using a cotton swab. The inserts were stained with 400 μl of cell staining solution and washed three times with water. Cells were treated with 200 μl of extraction solution and the solution transferred to individual wells of a new plate. The absorbance of the extracted samples was measured at 560 nm by a Thermo Multiskan Spectrum plate reader. Six astrocyte cultures from each group (AA, CA, AAG and CAG) were used in the assay and data were analyzed for significance with ANOVA (Prism 3.0 GraphPad software).
Rho activation assay
Rho activation assay kit (Upstate Biotechnology Billerica, MA, USA) was used to detect activated Rho in cell lysates. Unstimulated cells were cultured in 60 mm dishes until 85-90% confluence and then harvested in ice cold 1 × Mg2+ Lysis/Wash (MLB) buffer (according to the manufacturer's protocol). Protein concentration was determined by the Bradford method. Protein lysate (200 μg) were mixed with 32 μl of Rho assay reagent slurry containing GST-Rhotekin-RBD fusion protein, and incubated for 45 minutes at 4°C with gentle agitation. After pelleting and washing three times with 1 × MLB, the beads were resuspended in 2 × NuPage LDS sample buffer (Invitrogen Carlsbad, CA, USA) supplemented with 0.075 M DTT and boiled at 95°C for 5 minutes. Samples were subjected to western blot analysis. An anti-Rho antibody that recognizes Rho-A, Rho-B and Rho-C was used for detection. Four cultures from each group (AA, CA, AAG and CAG) were used in the assay. Western blots were performed in duplicate. Representative blots are shown in the results and the mean optical density was used in density analysis. Statistical significance was based upon two-tailed t-test (Prism 3.0 GraphPad software) and p-value < 0.05
Six eyes from normal CA donors, six eyes from normal AA donors, six eyes from CAG donors and four eyes from AAG donors were used. All donors were age matched. Tissues were fixed with 4% paraformaldehyde in 0.1 M phosphate-buffered saline pH 7.4 and processed for paraffin embedding. Two slides were stained per donor containing at least two 6 μm optic nerve sections each. In double labeling experiments we used monoclonal or polyclonal antibodies against human glial acidic fibrillar protein (GFAP) as an astrocyte marker. Secondary antibodies labeled with Alexa 488 and Alexa 568 (1:800) were from Molecular Probes. For negative controls, the primary antibody was replaced with non-immune serum. Serial sections used in comparisons (AAG versus CAG) were stained simultaneously to control for variations in immunostaining. Slides were examined in a Nikon Eclipse 80 i microscope (Tokyo, Japan) equipped with epifluorescent illumination and digital cameras (CoolSnap ES and CF, Photometrics). The images were processed using MetaMorph software (Molecular Devices Sunnyvale, CA, USA).
Additional data files
The following additional data files are available in the online version of the paper. Additional data file 1 is a table listing clinical information about CAG and AAG eyes used to generate primary cultures of ONH astrocytes. Additional data file 2 is a table listing demographic information of CA and AA normal donor eyes used to generate primary cultures of ONH astrocytes. Additional data file 3 is a table that summarizes the number of probe-sets on the chip and used in analysis. Additional data file 4 is spreadsheet listing the primers used for qRT-PCR. Additional data file 5 is a spreadsheet listing genes differentially expressed in glaucomatous ONH astrocytes and including the comparison between AAG versus normal AA. Additional data file 6 is a spreadsheet listing differentially expressed genes between CAG and normal CA. Additional data file 7 is a spreadsheet listing differentially expressed genes between AAG and CAG. Additional data file 8 is a spreadsheet listing genes differentially expressed in ONH astrocytes from AAG compared to both normal AA and CAG. Additional data file 7 is a spreadsheet summarizing Gene Ontology for the comparisons between AAG and AA data. Additional data file 8 is a spreadsheet with Gene ontology comparisons for CAG and CA. Additional data file 9 is a spreadsheet with GO comparisons for AAG versus CAG expression sets. Additional data file 10 is a figure showing the distribution of genes in two GO categories. Additional data file 11 is a figure showing qRT-PCR data that confirm additional differentially expressed genes from the CAG-CA comparison.
AA donor with glaucoma
CA donor with glaucoma
enzyme-linked immunosorbent assay
glial acidic fibrillar protein
myosin light chain kinase
optic nerve head
primary open angle glaucoma
transforming growth factor
tumor necrosis factor.
We thank Jose Bongolan for immunohistochemistry staining and Ping Yang and Marina Vracar-Grabar for generating cell cultures. This work was supported in part by NIH grant EY 06416 and an unrestricted grant from Research to Prevent Blindness.
- Weinreb RN, Khaw PT: Primary open-angle glaucoma. Lancet. 2004, 363: 1711-1720. 10.1016/S0140-6736(04)16257-0.PubMedView ArticleGoogle Scholar
- Friedman DS, Wolfs RC, O'Colmain BJ, Klein BE, Taylor HR, West S, Leske MC, Mitchell P, Congdon N, Kempen J: Prevalence of open-angle glaucoma among adults in the United States. Arch Ophthalmol. 2004, 122: 532-538. 10.1001/archopht.122.7.1019.PubMedView ArticleGoogle Scholar
- Nemesure B, He Q, Mendell N, Wu SY, Hejtmancik JF, Hennis A, Leske MC: Inheritance of open-angle glaucoma in the Barbados family study. Am J Med Genet. 2001, 103: 36-43. 10.1002/ajmg.1498.PubMedView ArticleGoogle Scholar
- Quigley HA, Vitale S: Models of open-angle glaucoma prevalence and incidence in the United States. Invest Ophthalmol Vis Sci. 1997, 38: 83-91.PubMedGoogle Scholar
- Broman AT, Quigley HA, West SK, Katz J, Munoz B, Bandeen-Roche K, Tielsch JM, Friedman DS, Crowston J, Taylor HR, Varma R, Leske MC, Bengtsson B, Heijl A, He M, Foster PJ: Estimating the rate of progressive visual field damage in those with open-angle glaucoma, from cross-sectional data. Invest Ophthalmol Vis Sci. 2008, 49: 66-76. 10.1167/iovs.07-0866.PubMedPubMed CentralView ArticleGoogle Scholar
- Leske MC, Wu SY, Hennis A, Honkanen R, Nemesure B: Risk factors for incident open-angle glaucoma: the Barbados eye studies. Ophthalmology. 2007, 115: 85-93. 10.1016/j.ophtha.2007.03.017.PubMedView ArticleGoogle Scholar
- Beck AD: Review of recent publications of the Advanced Glaucoma Intervention Study. Curr Opin Ophthalmol. 2003, 14: 83-85. 10.1097/00055735-200304000-00005.PubMedView ArticleGoogle Scholar
- Ederer F, Gaasterland DA, Dally LG, Kim J, VanVeldhuisen PC, Blackwell B, Prum B, Shafranov G, Allen RC, Beck A: The Advanced Glaucoma Intervention Study (AGIS): 13. Comparison of treatment outcomes within race: 10-year results. Ophthalmology. 2004, 111: 651-664. 10.1016/j.ophtha.2004.07.004.PubMedView ArticleGoogle Scholar
- Hernandez MR, Pena JD: The optic nerve head in glaucomatous optic neuropathy. Arch Ophthalmol. 1997, 115: 389-395.PubMedView ArticleGoogle Scholar
- Hernandez MR, Agapova OA, Yang P, Salvador-Silva M, Ricard CS, Aoi S: Differential gene expression in astrocytes from human normal and glaucomatous optic nerve head analyzed by cDNA microarray. Glia. 2002, 38: 45-64. 10.1002/glia.10051.PubMedView ArticleGoogle Scholar
- Yang P, Hernandez MR: Purification of astrocytes from adult human optic nerve heads by immunopanning. Brain Res Brain Res Protoc. 2003, 12: 67-76. 10.1016/S1385-299X(03)00073-4.PubMedView ArticleGoogle Scholar
- Kooperberg C, Aragaki A, Strand AD, Olson JM: Significance testing for small microarray experiments. Stat Med. 2005, 24: 2281-2298. 10.1002/sim.2109.PubMedView ArticleGoogle Scholar
- Ritchie ME, Silver J, Oshlack A, Holmes M, Diyagama D, Holloway A, Smyth GK: A comparison of background correction methods for two-colour microarrays. Bioinformatics. 2007, 23: 2700-2707. 10.1093/bioinformatics/btm412.PubMedView ArticleGoogle Scholar
- Zeeberg BR, Feng W, Wang G, Wang MD, Fojo AT, Sunshine M, Narasimhan S, Kane DW, Reinhold WC, Lababidi S, Bussey KJ, Riss J, Barrett JC, Weinstein JN: GoMiner a resource for biological interpretation of genomic and proteomic data. Genome Biol. 2003, 4: R28-10.1186/gb-2003-4-4-r28.PubMedPubMed CentralView ArticleGoogle Scholar
- Varela HJ, Hernandez MR: Astrocyte responses in human optic nerve head with primary open-angle glaucoma. J Glaucoma. 1997, 6: 303-313. 10.1097/00061198-199710000-00007.PubMedView ArticleGoogle Scholar
- Hernandez MR, Pena JD, Selvidge JA, Salvador-Silva M, Yang P: Hydrostatic pressure stimulates synthesis of elastin in cultured optic nerve head astrocytes. Glia. 2000, 32: 122-136. 10.1002/1098-1136(200011)32:2<122::AID-GLIA20>3.0.CO;2-J.PubMedView ArticleGoogle Scholar
- Pena JD, Varela HJ, Ricard CS, Hernandez MR: Enhanced tenascin expression associated with reactive astrocytes in human optic nerve heads with primary open angle glaucoma. Exp Eye Res. 1999, 68: 29-40. 10.1006/exer.1998.0577.PubMedView ArticleGoogle Scholar
- Liu B, Neufeld AH: Expression of nitric oxide synthase-2 (NOS-2) in reactive astrocytes of the human glaucomatous optic nerve head. Glia. 2000, 30: 178-186. 10.1002/(SICI)1098-1136(200004)30:2<178::AID-GLIA7>3.0.CO;2-C.PubMedView ArticleGoogle Scholar
- Yuan L, Neufeld AH: Tumor necrosis factor-alpha: a potentially neurodestructive cytokine produced by glia in the human glaucomatous optic nerve head. Glia. 2000, 32: 42-50. 10.1002/1098-1136(200010)32:1<42::AID-GLIA40>3.0.CO;2-3.PubMedView ArticleGoogle Scholar
- Salvador-Silva M, Aoi S, Parker A, Yang P, Pecen P, Hernandez MR: Responses and signaling pathways in human optic nerve head astrocytes exposed to hydrostatic pressure in vitro. Glia. 2004, 45: 364-377. 10.1002/glia.10342.PubMedView ArticleGoogle Scholar
- Lukas TJ: A signal transduction pathway model prototype II: Application to Ca2+-calmodulin signaling and myosin light chain phosphorylation. Biophys J. 2004, 87: 1417-1425. 10.1529/biophysj.104.042721.PubMedPubMed CentralView ArticleGoogle Scholar
- Tezel G, Hernandez MR, Wax MB: In vitro evaluation of reactive astrocyte migration, a component of tissue remodeling in glaucomatous optic nerve head. Glia. 2001, 34: 178-189. 10.1002/glia.1052.PubMedView ArticleGoogle Scholar
- Nguyen DH, Catling AD, Webb DJ, Sankovic M, Walker LA, Somlyo AV, Weber MJ, Gonias SL: Myosin light chain kinase functions downstream of Ras/ERK to promote migration of urokinase-type plasminogen activator-stimulated cells in an integrin-selective manner. J Cell Biol. 1999, 146: 149-164.PubMedPubMed CentralView ArticleGoogle Scholar
- Callaway K, Abramczyk O, Martin L, Dalby KN: The anti-apoptotic protein PEA-15 is a tight binding inhibitor of ERK1 and ERK2, which blocks docking interactions at the D-recruitment site. Biochemistry. 2007, 46: 9187-9198. 10.1021/bi700206u.PubMedView ArticleGoogle Scholar
- Vaidyanathan H, Ramos JW: RSK2 activity is regulated by its interaction with PEA-15. J Biol Chem. 2003, 278: 32367-32372. 10.1074/jbc.M303988200.PubMedView ArticleGoogle Scholar
- Hashimoto K, Parker A, Malone P, Gabelt BT, Rasmussen C, Kaufman PS, Hernandez MR: 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. 2005, 1054: 103-115. 10.1016/j.brainres.2005.06.050.PubMedView ArticleGoogle Scholar
- Birukov KG, Schavocky JP, Shirinsky VP, Chibalina MV, Van Eldik LJ, Watterson DM: Organization of the genetic locus for chicken myosin light chain kinase is complex: multiple proteins are encoded and exhibit differential expression and localization. J Cell Biochem. 1998, 70: 402-413. 10.1002/(SICI)1097-4644(19980901)70:3<402::AID-JCB13>3.0.CO;2-M.PubMedView ArticleGoogle Scholar
- Smith L, Parizi-Robinson M, Zhu MS, Zhi G, Fukui R, Kamm KE, Stull JT: Properties of long myosin light chain kinase binding to F-actin in vitro and in vivo. J Biol Chem. 2002, 277: 35597-35604. 10.1074/jbc.M206483200.PubMedView ArticleGoogle Scholar
- Kudryashov DS, Stepanova OV, Vilitkevich EL, Nikonenko TA, Nadezhdina ES, Shanina NA, Lukas TJ, Van Eldik LJ, Watterson DM, Shirinsky VP: Myosin light chain kinase (210 kDa) is a potential cytoskeleton integrator through its unique N-terminal domain. Exp Cell Res. 2004, 298: 407-417. 10.1016/j.yexcr.2004.04.025.PubMedView ArticleGoogle Scholar
- Gao L, Grant A, Halder I, Brower R, Sevransky J, Maloney JP, Moss M, Shanholtz C, Yates CR, Meduri GU, Shriver MD, Ingersoll R, Scott AF, Beaty TH, Moitra J, Ma SF, Ye SQ, Barnes KC, Garcia JG: Novel polymorphisms in the myosin light chain kinase gene confer risk for acute lung injury. Am J Respir Cell Mol Biol. 2006, 34: 487-495. 10.1165/rcmb.2005-0404OC.PubMedPubMed CentralView ArticleGoogle Scholar
- Gao L, Grant AV, Rafaels N, Stockton-Porter M, Watkins T, Gao P, Chi P, Munoz M, Watson H, Dunston G, Togias A, Hansel N, Sevransky J, Maloney JP, Moss M, Shanholtz C, Brower R, Garcia JG, Grigoryev DN, Cheadle C, Beaty TH, Mathias RA, Barnes KC: Polymorphisms in the myosin light chain kinase gene that confer risk of severe sepsis are associated with a lower risk of asthma. J Allergy Clin Immunol. 2007, 119: 1111-1118. 10.1016/j.jaci.2007.03.019.PubMedView ArticleGoogle Scholar
- Flores C, Ma SF, Maresso K, Ober C, Garcia JG: A variant of the myosin light chain kinase gene is associated with severe asthma in African Americans. Genet Epidemiol. 2007, 31: 296-305. 10.1002/gepi.20210.PubMedView ArticleGoogle Scholar
- Gutjahr MC, Rossy J, Niggli V: Role of Rho, Rac, and Rho-kinase in phosphorylation of myosin light chain, development of polarity, and spontaneous migration of Walker 256 carcinosarcoma cells. Exp Cell Res. 2005, 308: 422-438. 10.1016/j.yexcr.2005.05.001.PubMedView ArticleGoogle Scholar
- Kalman D, Gomperts SN, Hardy S, Kitamura M, Bishop JM: Ras family GTPases control growth of astrocyte processes. Mol Biol Cell. 1999, 10: 1665-1683.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang P, Agapova O, Parker A, Shannon W, Pecen P, Duncan J, Salvador-Silva M, Hernandez MR: DNA microarray analysis of gene expression in human optic nerve head astrocytes in response to hydrostatic pressure. Physiol Genomics. 2004, 17: 157-169. 10.1152/physiolgenomics.00182.2003.PubMedView ArticleGoogle Scholar
- Kim JS, Kim JG, Moon MY, Jeon CY, Won HY, Kim HJ, Jeon YJ, Seo JY, Kim JI, Kim J, Lee JY, Kim PH, Park JB: Transforming growth factor-beta1 regulates macrophage migration via RhoA. Blood. 2006, 108: 1821-1829. 10.1182/blood-2005-10-009191.PubMedView ArticleGoogle Scholar
- Fuchshofer R, Birke M, Welge-Lussen U, Kook D, Lutjen-Drecoll E: Transforming growth factor-beta 2 modulated extracellular matrix component expression in cultured human optic nerve head astrocytes. Invest Ophthalmol Vis Sci. 2005, 46: 568-578. 10.1167/iovs.04-0649.PubMedView ArticleGoogle Scholar
- Flanders KC, Ren RF, Lippa CF: Transforming growth factor-betas in neurodegenerative disease. Prog Neurobiol. 1998, 54: 71-85. 10.1016/S0301-0082(97)00066-X.PubMedView ArticleGoogle Scholar
- Lu Z, Lam KS, Wang N, Xu X, Cortes M, Andersen B: LMO4 can interact with Smad proteins and modulate transforming growth factor-beta signaling in epithelial cells. Oncogene. 2006, 25: 2920-2930. 10.1038/sj.onc.1209318.PubMedView ArticleGoogle Scholar
- Warner DR, Pisano MM, Roberts EA, Greene RM: Identification of three novel Smad binding proteins involved in cell polarity. FEBS Lett. 2003, 539: 167-173. 10.1016/S0014-5793(03)00155-8.PubMedView ArticleGoogle Scholar
- Etienne-Manneville S, Hall A: Integrin-mediated activation of Cdc42 controls cell polarity in migrating astrocytes through PKCzeta. Cell. 2001, 106: 489-498. 10.1016/S0092-8674(01)00471-8.PubMedView ArticleGoogle Scholar
- Frese S, Schubert WD, Findeis AC, Marquardt T, Roske YS, Stradal TE, Heinz DW: The phosphotyrosine peptide binding specificity of Nck1 and Nck2 Src homology 2 domains. J Biol Chem. 2006, 281: 18236-18245. 10.1074/jbc.M512917200.PubMedView ArticleGoogle Scholar
- Bauer K, Kratzer M, Otte M, de Quintana KL, Hagmann J, Arnold GJ, Eckerskorn C, Lottspeich F, Siess W: Human CLP36, a PDZ-domain and LIM-domain protein, binds to alpha-actinin-1 and associates with actin filaments and stress fibers in activated platelets and endothelial cells. Blood. 2000, 96: 4236-4245.PubMedGoogle Scholar
- Kotaka M, Kostin S, Ngai S, Chan K, Lau Y, Lee SM, Li H, Ng EK, Schaper J, Tsui SK, Fung K, Lee C, Waye MM: Interaction of hCLIM1, an enigma family protein, with alpha-actinin 2. J Cell Biochem. 2000, 78: 558-565. 10.1002/1097-4644(20000915)78:4<558::AID-JCB5>3.0.CO;2-I.PubMedView ArticleGoogle Scholar
- Andra K, Nikolic B, Stocher M, Drenckhahn D, Wiche G: Not just scaffolding: plectin regulates actin dynamics in cultured cells. Genes Dev. 1998, 12: 3442-3451. 10.1101/gad.12.21.3442.PubMedPubMed CentralView ArticleGoogle Scholar
- Pena JD, Taylor AW, Ricard CS, Vidal I, Hernandez MR: Transforming growth factor beta isoforms in human optic nerve heads. Br J Ophthalmol. 1999, 83: 209-218.PubMedPubMed CentralView ArticleGoogle Scholar
- Akutsu K, Morisaki H, Takeshita S, Sakamoto S, Tamori Y, Yoshimuta T, Yokoyama N, Nonogi H, Ogino H, Morisaki T: Phenotypic heterogeneity of Marfan-like connective tissue disorders associated with mutations in the transforming growth factor-beta receptor genes. Circ J. 2007, 71: 1305-1309. 10.1253/circj.71.1305.PubMedView ArticleGoogle Scholar
- Singh KK, Rommel K, Mishra A, Karck M, Haverich A, Schmidtke J, Arslan-Kirchner M: TGFBR1 and TGFBR2 mutations in patients with features of Marfan syndrome and Loeys-Dietz syndrome. Hum Mutat. 2006, 27: 770-777. 10.1002/humu.20354.PubMedView ArticleGoogle Scholar
- Mizuguchi T, Collod-Beroud G, Akiyama T, Abifadel M, Harada N, Morisaki T, Allard D, Varret M, Claustres M, Morisaki H, Ihara M, Kinoshita A, Yoshiura K, Junien C, Kajii T, Jondeau G, Ohta T, Kishino T, Furukawa Y, Nakamura Y, Niikawa N, Boileau C, Matsumoto N: Heterozygous TGFBR2 mutations in Marfan syndrome. Nat Genet. 2004, 36: 855-860. 10.1038/ng1392.PubMedPubMed CentralView ArticleGoogle Scholar
- Izquierdo NJ, Traboulsi EI, Enger C, Maumenee IH: Glaucoma in the Marfan syndrome. Trans Am Ophthalmol Soc. 1992, 90: 111-117.PubMedPubMed CentralGoogle Scholar
- Dale LB, Seachrist JL, Babwah AV, Ferguson SS: Regulation of angiotensin II type 1A receptor intracellular retention, degradation, and recycling by Rab5, Rab7, and Rab11 GTPases. J Biol Chem. 2004, 279: 13110-13118. 10.1074/jbc.M313333200.PubMedView ArticleGoogle Scholar
- Seachrist JL, Ferguson SS: Regulation of G protein-coupled receptor endocytosis and trafficking by Rab GTPases. Life Sci. 2003, 74: 225-235. 10.1016/j.lfs.2003.09.009.PubMedView ArticleGoogle Scholar
- Hicks SW, Machamer CE: Golgi structure in stress sensing and apoptosis. Biochim Biophys Acta. 2005, 1744: 406-414. 10.1016/j.bbamcr.2005.03.002.PubMedView ArticleGoogle Scholar
- Hicks SW, Horn TA, McCaffery JM, Zuckerman DM, Machamer CE: Golgin-160 promotes cell surface expression of the beta-1 adrenergic receptor. Traffic. 2006, 7: 1666-1677. 10.1111/j.1600-0854.2006.00504.x.PubMedView ArticleGoogle Scholar
- Xiao K, Oas RG, Chiasson CM, Kowalczyk AP: Role of p120-catenin in cadherin trafficking. Biochim Biophys Acta. 2007, 1773: 8-16. 10.1016/j.bbamcr.2006.07.005.PubMedView ArticleGoogle Scholar
- Shi S, Grothe S, Zhang Y, O'Connor-McCourt MD, Poole AR, Roughley PJ, Mort JS: Link protein has greater affinity for versican than aggrecan. J Biol Chem. 2004, 279: 12060-12066. 10.1074/jbc.M310091200.PubMedView ArticleGoogle Scholar
- Arslan F, Bosserhoff AK, Nickl-Jockschat T, Doerfelt A, Bogdahn U, Hau P: The role of versican isoforms V0/V1 in glioma migration mediated by transforming growth factor-beta2. Br J Cancer. 2007, 96: 1560-1568. 10.1038/sj.bjc.6603766.PubMedPubMed CentralView ArticleGoogle Scholar
- Johnson EC, Jia L, Cepurna WO, Doser TA, Morrison JC: Global changes in optic nerve head gene expression after exposure to elevated intraocular pressure in a rat glaucoma model. Invest Ophthalmol Vis Sci. 2007, 48: 3161-3177. 10.1167/iovs.06-1282.PubMedPubMed CentralView ArticleGoogle Scholar
- Kirwan RP, Leonard MO, Murphy M, Clark AF, O'Brien CJ: Transforming growth factor-beta-regulated gene transcription and protein expression in human GFAP-negative lamina cribrosa cells. Glia. 2005, 52: 309-324. 10.1002/glia.20247.PubMedView ArticleGoogle Scholar
- Agapova OA, Kaufman PL, Lucarelli MJ, Gabelt BT, Hernandez MR: Differential expression of matrix metalloproteinases in monkey eyes with experimental glaucoma or optic nerve transection. Brain Res. 2003, 967: 132-143. 10.1016/S0006-8993(02)04234-8.PubMedView ArticleGoogle Scholar
- Hernandez MR: The optic nerve head in glaucoma: role of astrocytes in tissue remodeling. Prog Retin Eye Res. 2000, 19: 297-321. 10.1016/S1350-9462(99)00017-8.PubMedView ArticleGoogle Scholar
- Schubert P, Morino T, Miyazaki H, Ogata T, Nakamura Y, Marchini C, Ferroni S: Cascading glia reactions: a common pathomechanism and its differentiated control by cyclic nucleotide signaling. Ann NY Acad Sci. 2000, 903: 24-33. 10.1111/j.1749-6632.2000.tb06346.x.PubMedView ArticleGoogle Scholar
- Hastings RH, Araiza F, Burton DW, Bedley M, Deftos LJ: Parathyroid hormone-related protein regulates apoptosis in lung cancer cells through protein kinase A. Am J Physiol Cell Physiol. 2004, 287: C1616-C1622. 10.1152/ajpcell.00300.2004.PubMedView ArticleGoogle Scholar
- Yu G, Swiston J, Young D: Comparison of human CAP and CAP2, homologs of the yeast adenylyl cyclase-associated proteins. J Cell Sci. 1994, 107: 1671-1678.PubMedGoogle Scholar
- Ge G, Greenspan DS: BMP1 controls TGFbeta1 activation via cleavage of latent TGFbeta-binding protein. J Cell Biol. 2006, 175: 111-120. 10.1083/jcb.200606058.PubMedPubMed CentralView ArticleGoogle Scholar
- Kuja-Panula J, Kiiltomaki M, Yamashiro T, Rouhiainen A, Rauvala H: AMIGO, a transmembrane protein implicated in axon tract development, defines a novel protein family with leucine-rich repeats. J Cell Biol. 2003, 160: 963-973. 10.1083/jcb.200209074.PubMedPubMed CentralView ArticleGoogle Scholar
- Hagino S, Iseki K, Mori T, Zhang Y, Hikake T, Yokoya S, Takeuchi M, Hasimoto H, Kikuchi S, Wanaka A: Slit and glypican-1 mRNAs are coexpressed in the reactive astrocytes of the injured adult brain. Glia. 2003, 42: 130-138. 10.1002/glia.10207.PubMedView ArticleGoogle Scholar
- Minami M, Kume N, Kataoka H, Morimoto M, Hayashida K, Sawamura T, Masaki T, Kita T: Transforming growth factor-beta(1) increases the expression of lectin-like oxidized low-density lipoprotein receptor-1. Biochem Biophys Res Commun. 2000, 272: 357-361. 10.1006/bbrc.2000.2778.PubMedView ArticleGoogle Scholar
- Leonoudakis D, Conti LR, Radeke CM, McGuire LM, Vandenberg CA: A multiprotein trafficking complex composed of SAP97, CASK, Veli, and Mint1 is associated with inward rectifier Kir2 potassium channels. J Biol Chem. 2004, 279: 19051-19063. 10.1074/jbc.M400284200.PubMedView ArticleGoogle Scholar
- Schuh K, Uldrijan S, Gambaryan S, Roethlein N, Neyses L: Interaction of the plasma membrane Ca2+ pump 4b/CI with the Ca2+/calmodulin-dependent membrane-associated kinase CASK. J Biol Chem. 2003, 278: 9778-9783. 10.1074/jbc.M212507200.PubMedView ArticleGoogle Scholar
- Olsen O, Moore KA, Fukata M, Kazuta T, Trinidad JC, Kauer FW, Streuli M, Misawa H, Burlingame AL, Nicoll RA, Bredt DS: Neurotransmitter release regulated by a MALS-liprin-alpha presynaptic complex. J Cell Biol. 2005, 170: 1127-1134. 10.1083/jcb.200503011.PubMedPubMed CentralView ArticleGoogle Scholar
- Hsueh YP, Wang TF, Yang FC, Sheng M: Nuclear translocation and transcription regulation by the membrane-associated guanylate kinase CASK/LIN-2. Nature. 2000, 404: 298-302. 10.1038/35005118.PubMedView ArticleGoogle Scholar
- Dimitratos SD, Stathakis DG, Nelson CA, Woods DF, Bryant PJ: The location of human CASK at Xp11.4 identifies this gene as a candidate for X-linked optic atrophy. Genomics. 1998, 51: 308-309. 10.1006/geno.1998.5404.PubMedView ArticleGoogle Scholar
- Irizarry RA, Bolstad BM, Collin F, Cope LM, Hobbs B, Speed TP: Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003, 31: e15-10.1093/nar/gng015.PubMedPubMed CentralView ArticleGoogle Scholar
- Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19: 185-193. 10.1093/bioinformatics/19.2.185.PubMedView ArticleGoogle Scholar
- 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-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView ArticleGoogle Scholar
- Stark C, Breitkreutz BJ, Reguly T, Boucher L, Breitkreutz A, Tyers M: BioGRID: a general repository for interaction datasets. Nucleic Acids Res. 2006, 34: D535-D539. 10.1093/nar/gkj109.PubMedPubMed CentralView ArticleGoogle Scholar
- Breitkreutz BJ, Stark C, Tyers M: Osprey: a network visualization system. Genome Biol. 2003, 4: R22-10.1186/gb-2003-4-3-r22.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.