|
|
||||||||
From the Casey Eye Institute, Oregon Health and Science University, Portland, Oregon.
| Abstract |
|---|
|
|
|---|
METHODS. Porcine TM cells were subjected to sustained mechanical stretching, and RNA was isolated after 12, 24, or 48 hours. Changes in gene expression were evaluated with microarrays containing approximately 8000 cDNAs. Select mRNA changes were then compared by quantitative reverse transcriptionpolymerase chain reaction (qRT-PCR). Western immunoblots were used to determine whether some of these changes were associated with changes in protein levels.
RESULTS. On the microarrays, 126 genes were significantly upregulated, and 29 genes were significantly downregulated at one or more time points, according to very conservative statistical and biological criteria. Of the genes that changed, several ECM regulatory genes, cytoskeletal-regulatory genes, signal-transduction genes, and stress-response genes were notable. These included several proteoglycans and matricellular ECM proteins composed of common repetitive binding domains. The results of analysis of mRNA changes in more than 20 selected genes by qRT-PCR supported the findings in the microarray analysis. Western immunoblots of several proteins demonstrated protein level changes associated with changes in the level of mRNA.
CONCLUSIONS. The expression of a variety of TM genes is significantly affected by mechanical stretching. These include several ECM proteins that contain multiple binding sites and may serve organizational roles in the TM. Several proteins that could contribute to the homeostatic modification of aqueous humor outflow resistance are also upregulated or downregulated.
Much of the aqueous humor outflow resistance appears to reside within the deepest portion of the TM.11 Consequently, this resistance is functionally stretched like a diaphragm across Schlemms canal. Alterations in the outflow resistance that are sufficient to affect IOP change the degree of stretching and distortion of the juxtacanalicular TM. A likely sensing mechanism in TM cells, indicating that an increase or decrease in outflow resistance is needed, is that of alterations in mechanical tension or stretching. This sensing could be mediated and transduced by integrinECM or similar types of interactions.12 13 We and others have found that trabecular cells can sense changes in IOP and mechanical stretching and produce several distinctive responses.14 15 16 17 18 19 20 21 22 23 24 In support of our overall hypothesis of IOP homeostasis, trabecular cells respond to pressure elevations in perfused organ culture or to mechanical stretching in cell culture by increasing MMP-2 (gelatinase A) and MMP-14 (membrane-type-1-MMP) activity or levels, while dramatically reducing levels of their primary tissue inhibitor, TIMP-2.14 Increased trabecular MMP activity results in reduced outflow resistance7 and thus should restore IOP to normal levels. Therefore, the components of a self-contained trabecular IOP homeostasis mechanism are present and functional. In perfused anterior segment organ culture, increased flow rates produce initial elevations of IOP, which return to normal over several days, even with sustained increases in perfusion rate.14 23
An important component of this putative IOP homeostasis mechanism is the nature of the changes in the ECM that adjust the outflow resistance. The ECM turnover process is initiated by secretion/activation of the MMPs, which partially degrades select ECM proteins. This phase is followed by the removal of proteolytic ECM fragments, presumably facilitated by TM cell endocytosis. New materials must then be synthesized by TM cells to replace degraded ECM components. To adjust the outflow resistance, changes in the amount, composition, or organization of this ECM are likely to be instituted. This process must occur within the active pathway of aqueous humor outflow without allowing excessive structural disorganization. To begin unraveling the details of this complex process, microarray analysis of TM gene expression changes occurring after sustained mechanical stretching were conducted.
| Materials and Methods |
|---|
|
|
|---|
TM Cell Culture, Mechanical Stretching, Treatments, and Extractions
Porcine TM cells were cultured as previously described.5 25 26 27 By passages 3 to 5, cells were plated at a density of approximately 90% confluence onto cell culture insert membranes in six-well culture plates.14 28 When cells were densely confluent, 3 to 5 days later, they were placed in serum-free medium for 24 hours before and during stretching experiments. To apply mechanical stretch/distortion, a smooth, 5.25-mm glass bead was placed in the dish beneath the center of the insert membrane. A weight was then applied to the lid of the plate, which forced the insert down onto the bead.14 28 This process created a defined upward distortion of the membrane, which increased the surface area and produced mechanical stretching of the cells and their ECM. More than 20 different porcine cell lines, each comprising cells from 10 to 20 eyes, were studied. DNA analysis (PicoGreen; Molecular Probes) was conducted in parallel wells in some studies, to estimate cell density. Because the differences between wells were always less than 10%, the analysis was not conducted in all studies. At the indicated times after initiation of stretching, medium was collected and stored in aliquots frozen at 20°C until use. For cellular or ECM protein analysis, membranes were immediately rinsed with ice-cold phosphate-buffered saline and the cells lysed with a modified RIPA buffer29 30 31 (2 mM EDTA, 2 mM EGTA, 1% NP-40, 0.5% sodium deoxycholate, 0.1% sodium dodecyl sulfate [SDS], 2 mM DTT, 1 mM PMSF, proteinase inhibitor cocktail, and 50 mM Tris [pH 7.5]) on ice. For RNA analysis, the membranes were cut from the inserts after stretching and placed in denaturation buffer (RNeasy kit; Ambion). They were then vortexed and the lysate removed. The remaining steps were as described in the kit protocol, and total RNA was processed according the manufacturers instructions. The extracted RNA was subjected to DNase treatment with 20 units of RNase-free DNase (Promega, Madison, WI) for 20 minutes. This was followed by two phenol-chloroform extraction steps to obtain high-quality RNA with an A260/A280 optical density ratio of 1.8 to 2.0.
Gene Chip Microarray
From each control or stretch sample, 9 µg of total purified RNA was provided to the Oregon Health and Science University Spotted Microarray Core (OHSU SMC) for processing (http://www.ohsu.edu/gmsr/smc/index.shtml). The stretch or the nonstretch control RNA was labeled by reverse transcription using either fluorescein-modified or biotinylated nucleotides and an oligo-dT primer. These labeled probes were mixed and hybridized to the spotted chips at 65°C for 16 hours (M-series LifterSlips; Erie Scientific, Portsmouth, NH) in a deep well hybridization chamber (TeleChem, Sunnyvale, CA). Hybridized arrays were then probed with Tyramide signal amplification (TSA) kits (PerkinElmer, Boston, MA), with horseradish peroxidase-conjugated anti-fluorescein antibody or streptavidin. Slides were developed with Cy3- and Cy5-tyramine and scanned (ScanArray 4000 XL with ScanArray Express software; Perkin-Elmer).
The two types of OHSU SMC microarray chips that were used had duplicate separated spottings of approximately 5700 or 8400 human cDNA clones per slide. The data presented herein are from the 8400 chips, but very similar results were obtained earlier with the 5700 chips. The libraries used were prepared by the IMAGE consortium32 and distributed in sequence-verified form by Invitrogen (Carlsbad, CA). Amplified PCR products were printed onto microarray slides (UltraGAPS; Corning Costar, Corning, NY) using microarray printing pins (TeleChem 16 CMP-3) and an array printer (Cartesian PixSys 5500 XL; Genomic Solutions, Irvine, CA). Twenty-three plant genes (Arabidopsis thaliana) spotted twice per chip, were used for standard reference providing both a positive and a negative control. In addition, 176 spots were left blank on each chip to provide another type of negative control. Scanned images were analyzed by computer (ImaGene; BioDiscovery, Marina del Rey, CA) to determine relative fluorescence intensity for each label; to identify spot locations, shapes, sizes, and anomalies; and to determine spot-specific background intensity for each label. At each of the three time points (12, 24, and 48 hours), three totally independent experiments were conducted. Each sample from each experiment was analyzed on two separate identical microchips. On each chip, each of the 8400 genes was spotted twice in separate regions. Consequently, at each of the three time points, there were 12 separate stretch and control data pairs for each gene.
Microarray spot fluorescence intensities for control and stretch samples with associated background values were then analyzed separately at each of the three time points. To determine statistical significance, the data from each chip with the duplicate spot values were separated into half chips, with one stretch and one control value and their respective specific backgrounds for each of the approximately 8400 cDNAs. The stretch and control data from each half-chip were then corrected for background and normalized via Lowess33 (GeneSpring software; Silicon Genetics, Redwood City, CA). These normalized fluorescence intensities were then combined into a single parallel data set for each time point, to identify those genes that achieved statistical significance using significance analysis of microarrays (SAM).34 Twelve normalized stretch and 12 normalized control fluorescence intensities for each gene at each time point were thus used to determine which changes were statistically significant. The number of permutations was set at 1000; the analysis was in a paired format; and the nearest number imputer was set to 10 in the SAM software (http://www.stat.stanford.edu/
tibs/SAM/ provided in the public domain by Stanford University, Stanford, CA). The false-discovery rate
parameter was then adjusted to produce a median value of less than 0.5 falsely significant genes from the complete set. A separate evaluation was conducted by using two criteria to determine which gene changes achieved biological significance. To do this, all 12 of the stretch and control fluorescence spot intensity pairs with their associated backgrounds were combined into a single long data set. They were corrected for background and then normalized via Lowess (GeneSpring; Silicon Genetics), as though they were replicates on one large microchip. These data were then filtered to select genes that exhibited at least a 1.5-fold average increase or at least a 50% average decrease. These were then further filtered to eliminate genes expressed at low levels both before and after mechanical stretching. Spots with average normalized fluorescence intensities below 300 for both control and stretch, after the background was subtracted, were eliminated. This lower expression threshold was selected because it was more than 2 standard deviations above the average specific background on these chips. The subset of genes that passed all three of these selection criteria (i.e., exhibited both statistical and biological significance for at least one time point) was retained.
Quantitative RT-PCR
Quantitative RT-PCR was performed (LightCycler and one-step Lightcycler-RNA amplification Kit, SYBR Green I; Roche Diagnostics, Indianapolis, IN).35 Primer pairs were designed in adjacent exons so that amplification of potential traces of contaminating genomic DNA would be easily identified. After the RT cycle (55°C, 10 minutes), 40 PCR cycles (94°C, 30 seconds; 50°C [primer dependent], 10 seconds; and 72°C, 14 seconds) were used. Fluorescence was acquired at the end of each extension step, and a melting curve was obtained at the end of each run. In a few cases in which multiple peaks were observed in the melting curve, the fluorescence was acquired in a separate step after each cycle at a temperature above the melting temperature of the additional nonspecific peak. After the PCR, products were analyzed on gels to verify band sizes and purity. In some cases, the product was also sequenced. Data were analyzed by the system software (LightCycler; Roche Diagnostics), by using a 10-point dilution standard curve and the crossing-point method to determine the relative template concentration of samples. Because fibronectin exhibits alternative splicing, several additional primers were designed to determine which splice forms were expressed. Twenty genes were evaluated by qRT-PCR, including several that increased, several that decreased, and several that did not show any change on the microarrays at any time after stretching.
Western Immunoblots
Western immunoblots, transferred electrophoretically from standard SDS-PAGE gels to nitrocellulose or PVDF membranes, were probed with the indicated primary antibodies. Detection used the appropriate secondary antibodies with conjugated horseradish peroxidase and chemiluminescence according to the manufacturers instructions (SuperSignal; Pierce)5 14 31 36 Autoradiographs were scanned and relative band density analyzed37 with a densitometry program (UVP, Upland, CA). For both the qRT-PCR and the Western immunoblot data, Students t-test or Mann-Whitney ranked sum analysis was used to determine significance, when comparing treatment results. All experiments presented were repeated at least three times, and typical gels were selected for presentation.
| Results |
|---|
|
|
|---|
parameter decreased incrementally, a lower plateau was reached at approximately a median number of 0.5 falsely significant genes per data set. Additional increments reduced the number of genes without reducing the false-discovery rate. Therefore, the highest number of genes that achieved this false-discovery rate was accepted. Further filtering of this group for biological significance, in terms of expression level in the TM and magnitude of changes with mechanical stretch, gave 126 genes that achieved a 1.5-fold or larger increase and 29 genes that achieved a 50% or greater decrease at one or more of the three time points. Analysis of the results from the SMC5700 chips showed good agreement for the genes common to both arrays (data not shown). The OHSU SMC microarray facility also conducted a detailed data analysis in which somewhat different specific normalization and significance testing approaches were used. Their lists of changed genes were very similar to those we obtained, as detailed earlier.
|
Effects of Mechanical Stretch on Specific TM Gene Expression
Twenty-seven transcripts that code for ECM molecules or for molecules integrally involved in ECM regulation or in cell adhesion were upregulated, and two were downregulated, according to the stringent statistical and biological criteria described earlier (Table 1) . Several were elevated at all three time points, but most showed significant increases at only one or two time points. Several of those that did not reach the significance criteria were nonetheless clearly affected. As an example, fibronectin 1 was elevated at 24 hours based on our statistical significance criteria, but only 1.33-fold, and thus is labeled NS (not significant) in Table 1 at this time point.
|
|
|
|
The chromosomal location of each gene that reached both statistical and biological significance in our studies is included in the last columns in Tables 1 2 3 and 4 . Analysis of the actual chromosomal location of each possibly relevant gene relative to the markers used to map these glaucoma loci as localized on the most recent human chromosomal maps, showed that most of the changed genes in Tables 1 2 3 and 4 are not actually within any of the mapped glaucoma loci. A summary of changed genes that are located within several mapped loci is shown in Table 5 . None were within the GLC1B, GLC1C, GLC1D, GLC1F, or GLC1G loci.40 41 42 43 44 Several changed genes were identified in some of the other loci that have been mapped but not named, as indicated in Table 5 .45 46 47
|
|
|
|
|
| Discussion |
|---|
|
|
|---|
The 8000 genes on this microarray represent approximately 25% of the human genome. Although it seems possible that a few of the porcine probe sequences do not hybridize with high affinity to the human cDNAs spotted on the microarrays, this number is probably small. Because human and porcine sequences are generally similar, microarray hybridization conditions are of only moderate stringency, and both the labeled porcine probes and the spotted human cDNAs are several hundred base pairs in length, lack of cross-species hybridization should be minimal in these studies. This selected subset of 126 increased and 29 decreased mRNAs should provide useful insights into this homeostatic process.
The genes that were changed by stretching displayed varied temporal patterns, with some changing at only one time point and others at two or all three points. One gene, IGF-binding protein 4 proteinase, was elevated 1.5-fold at 12 hours, not significantly changed at 24 hours, and significantly decreased (to 0.54-fold) at 48 hours. This temporal variation in gene expression is not surprising and suggests that the TM cell response to stretching is both complex and intricately coordinated. The distribution of the changes suggests a bias favoring only a few gene functional groups. Genes involved in ECM, cytoskeleton, or stress responses or in their regulation comprised a disproportionately large fraction of those we identified. The genes in this microarray are biased toward genes with known functions and include a wide representation of genes involved in most major biological processes. A fairly predictable and significant number of genes involved in transcription, transcriptional regulation, translation, or overall cellular regulation were significantly affected by stretching. Perhaps not too surprisingly, few genes involved in apoptosis, cellular proliferation, most aspects of cellular metabolism, or most other major cellular functions were affected. Stretching is presumably a normal homeostatic adjustment signal requiring only a moderate and focused response from these cells.
In the current literature, microarray and cDNA library studies of human TM cells subjected to elevated pressures in perfused organ culture are somewhat comparable to our results.22 24 48 The TM cells in these studies, which mostly involved experimental conditions different from ours, were probably responding to mechanical stretching that is functionally somewhat similar to the process we used. Matrix Gla protein, several metallothioneins, alpha B crystallin, alpha tubulin, transgelin, chaperonin-containing TCP1, and periostin were among the genes showing similar changes. Because of the very complex structural organization of the cells and the ECM within the TM in vivo, the relative degree to which different TM cells are exposed to stretching and distortion with increasing IOP is difficult to predict. It seems probable that each TM cell will experience different degrees of stretching and distortion with elevated IOP. In our cell culture stretch model, more uniformity is likely to be attained. However, different regions of the membranes in our model may also experience different degrees of stretch/distortion. We estimate that many of the cells in our model are experiencing degrees and types of stretching and distortion that are similar to that of many of the cells in the TM in vivo. However, it is very clear that additional detailed biophysical analysis of the TM is needed to clarify this contention.
In terms of our working hypothesis of IOP homeostasis, the changes in ECM or ECM-related genes listed in Table 1 are of particular interest. In addition to the initiation of ECM turnover by the concerted action of MMP-2, TIMP2, and MMP-14,14 28 subsequent biosynthetic replacement of the degraded ECM components is likely to be an integral aspect of the homeostatic IOP adjustment. This probably entails modifying the composition or organization of the ECM to change the outflow resistance and thus modulate IOP. One group of ECM genes that are affected by stretching includes: NELL2, tenascin C, SPARC, fibronectin, laminin
1 chain, and collagen XIV. A common feature of this group is that they are all modular matricellular ECM proteins composed of numerous repeat domains containing motifs, such as EGF-like, calcium-binding EGF-like, thrombospondin 1, fibronectin type III, von Willebrand factor C, heparin-binding or integrin-binding RGD motif.49 50 51 52 53 54 55 56 57 58 59 60 Domains of these types commonly serve as specific binding cassettes for proteinprotein, proteinglycosaminoglycan, or proteincell interactions and function in ECM and tissue structure and organization. Their elevation during stretching is likely to be involved in maintaining and facilitating the tissue and ECM reorganization associated with the ECM turnover process initiated by the MMPs. Because they bind glycosaminoglycans (GAGs) and proteoglycans, which are thought to provide at least a portion of the outflow resistance,3 4 61 they could also directly affect outflow by changing resistance component orientation.
The increase in fibronectin levels and the splicing changes, which increase the fraction of the protein that contains V1 and V3 domains, are intriguing. Perfusion of anterior segment organ cultures with the second heparin binding (Hep II) fragment of fibronectinthe 14th to 16th type III domains as indicated in Figure 4 reversibly increases outflow facility.57 In another tissue, it has been shown that alternative splicing of the V1 and V3 domains affect several specific biological activities of the adjacent Hep II fragment without affecting other activities.56 Fibronectin has an RGDS cell-binding site in the 11th type III domain that binds to
5ß1 integrins (Fig. 4) . An
4ß1 integrin-binding site is found in the distal portion of the Hep II region. In the alternatively spliced variable regions of IIICS, V1 contains an LDV and V3 contains an REDV, both of which are also
4ß1 integrin-binding motifs.62 63 Thus, TM cells switch from a mix of splice variants, with and without these two sites, to exclusively expressing the form that contains these two additional cell-binding motifs. In addition, the V2 domain of fibronectin contains a GAG-binding site, and these splice variants may impact binding to CD44 and syndecan 2 (both discussed later). Although we did not see changes in fibronectin EIIIA or EIIIB exon splicing with stretching, we detected very low levels of EIIIB in both stretched and control cells. A previous report did not detect any of either of these exons in serum-free control TM cells, although both were detected after TGF-ß2 or serum treatment.64
Of interest, laminin
1 was increased 1.5-fold at 12 hours, but
2,
4,
5, ß1, and the laminin receptor 1 were totally unaffected at all three time points. The laminins are a central component of cellular basement membranes, acting as scaffolds for the recruitment of other ECM components.65 Because laminins 1 to 15 are different heterotrimeric combinations of the
1 to 5, ß1 to 4, and
1 to 3 subunits, an increase in
1 may reflect a partial shift from one laminin form to another,58 at a time when the ECM is undergoing remodeling and change. The
1-subunit of laminin is also involved in binding to nidogen and in polymer formation.66 Although the primary cell-binding sites on laminin are in the laminin
-subunit, cell-binding sites have been identified on the
1 subunit of laminin.67 68 69 Based on relative fluorescence intensities in our microarrays, the
2 subunit was expressed at low to negligible levels, ß1 only at moderately low levels; but
4 and
5,
1, and laminin receptor 1 were all expressed at relatively high levels by the TM.
A second group of ECM genes that changed expression with stretching includes several proteoglycans, some of which may serve a direct role in outflow resistance.3 4 6 61 Chondroitin sulfate proteoglycan 4, fibromodulin, biglycan, syndecan 2, and the part-time proteoglycan and hyaluronan GAG receptor, CD44, increased with stretching, whereas mimecan decreased. Presumably, one likely way to adjust the outflow resistance in response to mechanical stretching would be to adjust the proteoglycan/GAG levels or composition in the TMs ECM. Although it is tempting to speculate on the molecular mode by which increasing one or more of these proteoglycans or of decreasing mimecan, could decrease the outflow resistance, more information is needed to understand the process. Changes in the GAG side chains are also probable modifiers of the outflow resistance.
The increase in CD44, the hyaluronan receptor, is very likely to be important in this process. A relationship between hyaluronan levels and the trabecular outflow resistance is generally accepted.70 Several recent studies have demonstrated the involvement of CD44 in TM cellular function and have provided hypotheses explaining relationships with the outflow resistance.71 72 CD44 is also integrally involved and interactive with MMP-2 and -14,73 74 75 both of which are upregulated at the translational but not the transcriptional level in TM cells by mechanical stretching.14 28 CD44 exists in forms, with and without GAG side chains, and the form with chondroitindermatan sulfates binds to collagen XIV, which is also increased with stretching.
Syndecan 2, a transmembrane heparan sulfate proteoglycan,76 77 was upregulated at 24 hours (Table 1) . It interacts with
5ß1 integrin in cellular signaling and with the Hep II, V1, and V2 regions of fibronectin. This upregulation also affects cytoskeletal organization in an ezrin-dependent step via Rho A GTPase activation.78 79 Syndecan 2 is phosphorylated on two C-terminal tyrosines by the kinase domain of EphB2 receptors as a step in signaling to the cytoskeleton.77 Although EphB2 levels were not upregulated, its ligand, ephrin-B2, was upregulated by stretching. Thus, several proteins that have important interactions with each other are changed, further increasing the relevance of this process. CD44 also contains a cytoplasmic ezrin-binding domain, linking it to this same signaling system and to the Rho A-Rho kinasemediated cytoskeletal regulation.80 81 82 83 Therefore, the syndecan 2, ephrin-B2, and CD44 upregulation by stretching may be closely linked in maintaining the TMs ECM and may be involved in, among other events, the Rho Acytoskeletal modulatory changes, which also occurred with mechanical stretching of TM cells (Table 2) .
The cytoskeletal changes observed (Table 2) involve components of all three filament systemsmicrofilaments, microtubules, and intermediate filamentswhich suggests significant cross-talk between systems.84 Early cytoskeletal changes after mechanical stretching of TM cells have been reported,15 and sustained responses and long-term differences can be expected. Agents that manipulate the TM cell cytoskeleton have been shown to affect outflow facility.85 86 87 88 89 Whether the cytoskeletal changes are a direct cause or an effect of the ECM changes and outflow resistance adjustments remains unclear.
Several other intriguing gene expression changes are notable, such as increases in the chemokines, metallothionines, and metastasis-associated protein 1. The magnitude of some of these changes suggests important, if only partially understood, regulatory processes. Several metallothionein mRNAs increase from low or negligible levels to quite high levels, increasing between 20- to more than 100-fold at one time point. Functions of these genes in other tissues include heavy metal, particularly zinc, detoxification90 and homeostasis.91 This suggests that TM cells are responding to elevated or released zinc or to oxidative stress in association with the stretching.92 A zinc transporter and several other stress genes are also induced. The several chemokine ligand increases and the large metastasis-associated 1 increases are similarly intriguing and similarly enigmatic.
Although the number of genes with expression levels that are significantly changed by mechanical stretch is not large, the potential involvement of many of them in maintaining and adjusting the TMs ECM and putatively the outflow resistance is likely. Understanding this complex process that appears to be responsible for maintaining IOP homeostasis clearly requires additional studies focused on the function of these proteins in outflow pathway biology.
| Footnotes |
|---|
Submitted for publication January 21, 2005; revised April 11, 2005; accepted April 22, 2005.
Disclosure: V. Vittal, None; A. Rose, None; K.E. Gregory, None; M.J. Kelley, None; T.S. Acott, Alcon Labs (F)
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Ted S. Acott, Casey Eye Institute (CERES), Oregon Health and Science University, 3375 SW Terwilliger, Portland, OR 97239-4197; acott{at}ohsu.edu.
| References |
|---|
|
|
|---|
. Invest Ophthalmol Vis Sci. 2000;41:422430.
B-crystallin: an early cellular response to mechanical stretch. Biochem Biophys Res Commun. 1997;235:6973.[CrossRef][ISI][Medline][Order article via Infotrieve]
regulation of trabecular matrix metalloproteinases and TIMPs. Invest Ophthalmol Vis Sci. 2001;42:28312838.