|
|
||||||||
From the Michael Panitch Macular Degeneration Laboratory, Wilmer Eye Institute, Johns Hopkins Medical Institutes, Baltimore, Maryland.
| Abstract |
|---|
|
|
|---|
METHODS. Morphologically normal RPE cells from 15 human globes from donors aged 52 to 82 years old were laser capture microdissected. Total RNA from 5000 cells was SMART amplified, [33]P-labeled, and hybridized to a cDNA array containing 4325 known genes. Expression profiles were analyzed by hierarchical cluster analysis, Prediction Analysis of Microarrays (PAM), and Significance Analysis for Microarrays (SAM). Differentially expressed genes were evaluated further by real time RT-PCR.
RESULTS. The overall expression profiles of RPE cells from the macula and periphery were similar. Unsupervised and supervised hierarchical cluster analysis showed that patient genotype was a stronger separating factor than topographical location. SAM analysis identified 11 genes that were underexpressed by macular RPE cells. The expression patterns of these 11 genes were confirmed by real time RT-PCR, with 5 genes reaching statistical significance.
CONCLUSIONS. Whereas the overall expression profiles were similar between cells from the macula and periphery, subtle differential expression of five genes could contribute to RPE phenotypic differences based on topographic location.
A hypothesis that the mRNA phenotype of the RPE varies by topographical location was made. The onset of microarray technology has resulted in a rapid and more comprehensive molecular characterization than previously for a number of diseases.10 11 12 Technical issues such as difficulty in dissecting macular from peripheral RPE cells, the inherent heterogeneous phenotype of RPE cells within a macula, particularly in elderly eyes,2 and the resultant small amount of tissue available for molecular studies, have been impediments that prevent determining accurate mRNA phenotypes. The development of laser capture microdissection allows for the removal of pure cell populations from small regions of tissue, such as the macula, for microarray studies.13 Since little is known about the expression profiles of RPE cells in vivo, the expression profile of morphologically normal RPE cells microdissected from the macula and periphery of human donor globes was characterized.
| Materials and Methods |
|---|
|
|
|---|
|
Laser Capture Microdissection
Sample sections from each sample were examined before laser capture microdissection, and three observers confirmed that RPE morphologic abnormalities, drusen, or basal deposits were not observed. Cells of interest were dissected with an Arcturus PixCell II laser capture microdissector (Arcturus Engineering, Inc., Mountain View, CA) using transfer film (Cap-Sure TF-100; Arcturus Engineering) according to our previously published protocol.16 Morphologically normal RPE cells attached to nonthickened Bruchs membrane were defined using the criteria established by Curcio et al. and Sarks.17 18 Normal macular RPE were defined as having regular cuboidalcolumnar cell shape, homogeneous melanin pigmentation, and a height estimated at 10 to 15 µm19 using the 7.5 and 15 µm spot size of the laser aiming beam. Peripheral RPE typically have a cuboidal cell shape and a lower cell height. Normal peripheral RPE were defined as having a cuboidal cell shape, homogenous melanin pigmentation, and a height of
7.5 µm. In addition, cells were included as normal only if they were overlying a nonthickened Bruchs membrane. For this age group, Okubo et al.18 20 have reported a thickness of 3 to 5 µm, so we defined a nonthickened Bruchs membrane as <
RPE cell height and not associated with drusen or other Bruchs membrane abnormality. After dissection, the transfer cap was inspected with the microscope for contaminating tissue, which verified a cleavage plane at the RPEBruchs membrane junction, before being placed in 200 µL denaturing buffer that contained 4 M guanidine isothiocyanate, 0.02 M sodium citrate, 0.5% sarcosyl, and 2 µL ß-mercaptoethanol (14.5 M; Qiagen Inc, Valencia, CA).
RNA Extraction
Total RNA was extracted from laser captured RPE cells using the RNeasy Mini-kit (Qiagen Inc.) according to the manufacturers recommendations. RNA was treated with DNase I (Qiagen, Inc.) during RNA purification. A sample of RPE cells was obtained from a peripheral calotte that was not used for microarray analysis from each donor which showed preserved 28S and 18S rRNA bands. Before synthesizing probe, RNA quality was assessed by the expression of GAPDH from 100 cells using real time RTPCR with primers designed at the 5' end of the gene.
Probe Synthesis
Probe was prepared from total RNA of 5000 laser captured RPE cells using the Super SMART PCR cDNA synthesis kit (BD Biosciences Clontech, Palo Alto, CA) according to the manufacturers recommendations. Total RNA was reverse transcribed, and column-purified first strand cDNA was PCR amplified. The PCR cycle number was optimized using a small aliquot (5 µL) after 15 cycles and every three cycles thereafter, with 1.2% agarose/ethidium bromide gel electrophoresis. Typically, 23 cycles produced ds cDNA that remained in the exponential phase of amplification that produced a smear from 0.5 to 5 kb. The cDNA was column purified with a NucleoSpin Extraction kit (BD Biosciences Clontech) and labeled using the BD Atlas SMART probe amplification kit (BD Biosciences Clontech) in the presence of 50 µCi [
-33P]dATP, 1 µg random hexamers, and 2 units Klenow fragment at 50°C for 30 minutes. The probe was purified by passage through a Bio-Spin 6 chromatography column (BioRad Laboratories, Hercules, CA).
Microarray Analysis
The labeled cDNA was denatured and hybridized to the cDNA GeneFilter "Named genes" human array (4325 genes; Invitrogen/Research Genetics, Inc., Huntsville, AL) using the manufacturers protocol. This array contains genes with known function that are an insert DNA from a sequence-verified IMAGE/LLNL clone from the 3' end of the gene. Arrays were exposed for 3 days to a high density phosphorimager screen (BioRad Laboratories) and scanned at 50 µm resolution in a phosphorimager instrument (FX Pro-Plus; BioRad Laboratories).
Image and Statistical Analysis
The TIFF images acquired from the phosphorimager were imported into the image analysis software (Pathways 3; Invitrogen/Research Genetics, Inc.). This software aligns the images, quantifies a signal intensity for each gene, and normalizes the different hybridization signals on the basis of the 75% average signal intensity of the entire array. To allow statistical comparison of the arrays, the gene expression signals were scaled according to the method of Tusher et al., and our previously published protocol.21 22
Hierarchical cluster analysis of macular and peripheral RPE cells was determined with cluster analysis and visualized with TreeView.23 A class prediction from gene expression profiling based on the "nearest prototype (centroid) classifier" approach, was performed with the Prediction Analysis of Microarrays (PAM) software, to identify gene subsets that best characterize each class, i.e., macula and periphery.24 This program uses soft thresholding rather than screening, and focuses on misclassification error.24 The Significance Analysis of Microarrays method (SAM, version 1.12) was used to determine individual gene expression differences by topographical distribution.21 SAM calculates a significance score for each gene based on the gene expression change relative to the SD of repeated values.
Real Time Reverse TranscriptionPolymerase Chain Reaction (RT-PCR)
An aliquot (1 ng) from the same ds cDNA used for microarray analysis was assayed using the LightCycler apparatus (Roche Diagnostics, Nutley, NJ). The primer sequences used in this study were designed using Primer 3 (Whitehead Institute/MIT, Cambridge, MA) or the LightCycler Probe Design software, and sequences were verified using NCBI Unigene (www.ncbi.nlm.hih.gov/) (Table 2) . The standard curve consisted of PCR products for the gene of interest using serial dilutions of 1 pg106 pg. Thermo cycling of each reaction was performed in a final volume of 20 µL containing SYBR Green PCR Master Mix (10 µL; Qiagen, Inc.), Primer A and B (10 µM each), and 2 µL template DNA in a concentration of 2.5 mM MgCl2. The cDNA was denatured at 95°C for 15 minutes followed by PCR settings of 94°C for 15 seconds, Tm-5°C for 20 seconds, and 72°C for x seconds where x = PCR product length (bp/20). PCR products were quantified using the second derivate maximum values calculated by the Light-Cycler analysis software. Negative controls without template were produced for each run. Expression levels of all genes were normalized to GAPDH mRNA levels. All PCR products were checked by melting point analysis. For each sample, the experiment was repeated once, and the average expression of each sample was used to calculate the expression ratio. The Wilcoxon signed rank test was used to compare the differential gene expression between macular and peripheral RPE. P < 0.05 was considered significant.
|
| Results |
|---|
|
|
|---|
|
|
|
|
|
Fifty-six genes that are involved in the main oxidative stress defense systems are represented on this array. The normalized average expression values are listed as seen in Table 4 . The most highly expressed antioxidant genes include superoxide dismutase (SOD), and genes from the glutathione antioxidant system including glutathione-S-transferase (GST), glutathione peroxidase (GPX),
-glutamyltransferase-like activity 1, peroxiredoxin, and glutathione synthesis. A complete set of isoforms for several of these genes were represented on the array which allows a characterization of the relevant isoenzymes expressed by the RPE. For example, SOD1 and SOD2, which detoxifies superoxide, had higher expression than SOD3. GPX 3, GST A4, M1, M4, O1 and peroxiredoxin-1 were the most highly expressed isoforms. Supervised cluster analysis using this gene set also showed no difference between RPE cells derived from the macula or periphery.
|
Gene Expression Differences between Macular and Peripheral RPE
Individual gene expression differences between macula and periphery were assessed by SAM analysis. With a false discovery rate of 9%, two-class unpaired SAM analysis identified 11 genes that had significant differential expression (Table 5) . The fold-differences in expression between macular and peripheral cells were small with a maximum of 4.2-fold. All genes were downregulated by macular RPE cells, and have functions related to apoptosis, differentiation, oxidative stress, metabolism, and matrix regulation.
|
| Discussion |
|---|
|
|
|---|
The ability to assess RNA quality and a sample quantity that misrepresents the intended cell population are potential limitations of laser capture microdissection. Donors were chosen who had limited premortem illness since the rapidity of death is a greater influence on RNA quality than postmortem factors.14 The arrays were selected because they contain 700 to 1000 BP cDNA inserts from the 3' end of the gene, the region which is most resistant to RNA degradation. We found similar expression of GAPDH on the arrays and by real-time RT-PCR using primers designed from the 5' end, a wide size distribution of amplified cDNA products, and no signal degradation of total arrays or individually selected genes when plotted against postmortem factors. While these analyses do not quantify the degree of RNA degradation, they do indicate that this set of globes had similar quality RNA. Amplification can potentially induce bias. Seth et al. however, showed that the relative expression of low, medium, and high abundance genes was retained after SMART amplification independent of transcript abundance, coding region, and PCR product size.27 Amplification can allow a reproducible microarray signal from as few as 10 cells. Since several studies have demonstrated significant variability in gene expression across RPE cells in vivo,6 8 9 we were reluctant to reduce the cell number to this level over concern of obtaining an expression profile that misrepresents the intended cell population. By using 5000 cells for our analysis, we acknowledge that the expression profiles identified in this study may not be generalizable to all RPE cells or alternatively, the cell population studied may not represent cells that are vulnerable to preferential aging changes within the macula or periphery.
These data provide insights into the genes that maintain RPE homeostasis. The most highly abundant genes were common to RPE cells regardless of topographical location, and illustrate the wide functional diversity of the RPE cell. Of the most highly expressed genes, 26% are involved in protein synthesis and degradation. Sharon et al., 28 using SAGE analysis of native RPE from a single donor, found that 30% of the most abundant tags were related to protein synthesis and degradation, possibly due to the RPEs high phagocytic activity. Buraczynska et al.29 characterized a native RPE library of over 1100 genes and expressed sequence tags (ESTs) from two donors. One hundred and sixty-seven genes from this data set were identified on the array, and all were expressed by RPE cells, regardless of location. While a complete comparison with these libraries is difficult due to differences in technique, the present results along with these previous studies, help to define the mRNA phenotype of native RPE cells.
The RPE is exposed to significant oxidative stress due to its high metabolic activity, high oxygen fluxes from the choriocapillaris, and marked light exposure, and as a result, has a significant antioxidant defense system. A complete isoform set for several important antioxidant enzymes were evaluated in this study. The high expression suggests that SOD1 and SOD2, GPX3, GST A4, M1, M4, O1, and peroxiredoxin 1 are relevant isoenzymes in the RPE cells oxidative defense system. While high catalase expression has been previously reported,30 our low catalase expression is more consistent with reports of an age-related decline by the RPE.30 31 The RPE is susceptible to oxidative damage by lipid peroxidation products from phagocytosed outer segments.32 Lipid peroxidation products are detoxified not by SOD and catalase, but instead by the GST alpha class and GPXs.33 34 35 36 The high expression of GSTA4 and GPX3 suggests that these isoforms could be relevant lipid peroxidation detoxifiers for the RPE.
Our hierarchical cluster analysis found that the overall expression profiles between macular and peripheral RPE were similar, and suggested that the genotype was a stronger factor than topography. Alternatively, race, gender, or exposure to similar environmental conditions such as UV exposure, diet, associated co-morbidities, or medications, could have influenced the gene expression profiles. Oleksiak et al. 37 showed that, despite considerable expression profile variability among individuals within the same population, many important expression differences were small. The authors reasoned that genes relevant to a phenotype are tightly regulated and have constant expression, so that small expression changes will produce biologically important differences.37 This supposition has particular relevance when determining topographically-related changes since we would expect to find small expression changes by regional location. In the present study, the SAM program identified eleven genes with < 4.2-fold expression differences. Due to this small expression difference, statistically validated real time RT-PCR results were used to determine genes that may distinguish macular from peripheral cells. While the real time PCR experiments validated the expression pattern of all 11 genes tested, only five genes were statistically validated. Besides the relatively small expression differences between topographical location, the variable expression of GAPDH which was used for normalization, could be a confounding factor. Recent studies have suggested that there is no housekeeping gene that is optimally normalizes gene expression, and that using total RNA may be a more accurate method to normalize gene expression.38 39 The small starting material from laser captured material prevents accurate measurement of total RNA from samples. While no statistical difference between macular and peripheral RPE GAPDH expression by the arrays and RT-PCR (data not shown) was found, the variability of GAPDH expression across donors could have influenced the differential expression results for real time RT-PCR validation experiments. Regardless, these results demonstrate the importance of statistical assessment for RT-PCR validation studies, the need for finding the optimal method for normalizing expression, and suggest that the genes with statistically significant differential expression are worth exploration to determine a role in distinguishing macular from peripheral RPE. The translation of an mRNA into a protein can be highly variable, and post-translational protein modifications are important changes that influence aging. Ultimately, these gene expression changes need to be correlated with their respective protein levels, and more importantly, how these changes influence their biological function.
The macula experiences more apoptotic RPE cell loss than the periphery with aging.3 4 5 The SAM and RT-PCR analysis revealed that the cell cycle gene c-KIT was underexpressed by macular cells. c-KIT has been linked to bcl-2 upregulation, and when downregulated, makes cells susceptible to apoptosis.40 Likewise, since cysteine-rich, angiogenic inducer 61 is a cytokine that promotes cell proliferation, its under-expression by macular cells would also support preferential cell loss in the macula.41
Decreased macular GSTM1 expression was revealed by SAM and RT-PCR analyses. The GST Mu class is highly polymorphic. For example, 50% of the white population lacks GSTM1 activity due to two GSTM1 null alleles.42 Patients with the GSTM1, but not the P1, T1, or Z1 null genotype, do not neutralize photo-oxidative stress, and are highly susceptible to solar keratosis,43 or if they smoke, are at increased risk for atherosclerosis from an impaired ability to detoxify tobacco smoke, and have increased oxidatively-induced DNA damage within atherosclerotic lesions.44 45 It is possible that reduced GSTM1 activity, whether from genetic susceptibility or age-related decline, could make macular RPE susceptible over time, to either tobacco related or photo-oxidative stress, which coincidently, are risk factors for age-related macular degeneration.46 47
Compared to peripheral cells, Watzke et al. 2 observed a preferential degeneration in macular RPE cell morphology with aging. Our analysis identified reduced expression of aldehyde dehydrogenase 6 (ALD6). ALD6 is an essential enzyme in the synthesis retinoic acid, which has an established role in epithelial cell differentiation.48 The reduced macular expression of ALD6 could support preferential aging related morphologic changes to macular RPE cells via alterations in the retinoic acid pathway.
Genes were examined that, in general, contribute to the overall homeostasis of a number of cell types. This approach may be justified since to date, little is known about the overall RPE mRNA phenotype. A shortcoming of this strategy, of course, is that RPE specific genes were not examined. Currently, one available human RPE specific library contains 1100 nonredundant genes.29 We are currently evaluating in detail, the expression signature of macular RPE cells using an expanded RPE specific array. Whereas the functional effects of these differentially expressed genes remain unexplored, the results of this study provide a foundation for studying differences in macular RPE cells as a function of topography.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication February 18, 2004; revised April 29, 2004; accepted June 1, 2004.
Disclosure: K. Ishibashi, None; J. Tian, None; J.T. Handa, None
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: James T. Handa, 3109 Jefferson St. Building, 600 N. Wolfe St, Johns Hopkins Medical Institutes, Baltimore, MD 21287; jthanda{at}jhmi.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
N. Nagai, A. Klimava, W.-H. Lee, K. Izumi-Nagai, and J. T. Handa CTGF Is Increased in Basal Deposits and Regulates Matrix Production through the ERK (p42/p44mapk) MAPK and the p38 MAPK Signaling Pathways Invest. Ophthalmol. Vis. Sci., April 1, 2009; 50(4): 1903 - 1910. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. K. Chen, P. J. Patel, G. S. Uppal, G. S. Rubin, P. J. Coffey, G. W. Aylward, and L. Da Cruz A Comparison of Macular Translocation with Patch Graft in Neovascular Age-Related Macular Degeneration Invest. Ophthalmol. Vis. Sci., April 1, 2009; 50(4): 1848 - 1855. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Bowes Rickman, J. N. Ebright, Z. J. Zavodni, L. Yu, T. Wang, S. P. Daiger, G. Wistow, K. Boon, and M. A. Hauser Defining the Human Macula Transcriptome and Candidate Retinal Disease Genes Using EyeSAGE Invest. Ophthalmol. Vis. Sci., June 1, 2006; 47(6): 2305 - 2316. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Yang, J. L. Wiser, J. J. Peairs, J. N. Ebright, Z. J. Zavodni, C. Bowes Rickman, and G. J. Jaffe Human RPE Expression of Cell Survival Factors Invest. Ophthalmol. Vis. Sci., May 1, 2005; 46(5): 1755 - 1764. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |