|
|
||||||||
From the Departments of Ophthalmology and Visual Sciences and Human Genetics, W. K. Kellogg Eye Center, University of Michigan, Ann Arbor, Michigan.
| Abstract |
|---|
|
|
|---|
METHODS. Gene microarray slides containing 2400 human genes (primarily neuronal) were hybridized to biotin or dinitrophenyl (DNP)-labeled target cDNAs that were synthesized using total RNAs from young (1314 years) and elderly (6274 years) human retinas. Hybridization signals were visualized with cyanine (Cy)-5 or Cy-3 fluorescent reporter molecules, and the fluorescence intensities of the images were analyzed by computer. Northern blot analysis and real-time quantitative reverse transcription PCR (qRT-PCR) were performed to validate the microarray results.
RESULTS. Of the 2400 genes represented on the microarray slides, more than 50% hybridized to the retinal cDNA targets. Expression of a majority of these genes was not altered during aging; nonetheless, changes in the expression of 24 genes were detected between young and elderly retinas. These genes could be clustered into four categories: energy metabolism, stress response, cell growth, and neuronal transmission/signaling. Northern blot analysis and qRT-PCR results confirmed the changes in expression of 8 of 10 genes examined.
CONCLUSIONS. Using commercially available slide microarrays, the authors show that aging of the human retina is associated with changes in patterns of gene expression. This analysis suggests that pathways involved in stress response and energy metabolism play key roles in retinal aging. These studies demonstrate the utility of gene microarrays in identifying global patterns of retinal gene expression and lay the foundation for future studies defining the genetic basis of aging-associated retinal diseases, such as age-related macular degeneration.
| Introduction |
|---|
|
|
|---|
Aging in humans is associated with progressive and perhaps irreversible impairment of physiological functions, including vision. Age-related macular degeneration (AMD) is a major cause of untreatable vision loss in the elderly. The origin of this disease is dependent on complex interactions between genetic and environmental risk factorsthe two strongest risk factors being age and family history.3 Vision loss in AMD is attributed to photoreceptor dysfunction, which is caused by abnormalities of the retinal pigment epithelium (RPE), Bruchs membrane, and/or choriocapillaris.4 The physiological changes in the aging retina are remarkably similar, albeit less severe, to the pathologic changes in AMD tissue. In both aging and AMD, multiple factors, including DNA damage and oxidative stress, are believed to contribute to the pathology. In light of these studies, we hypothesized that a broader understanding of the molecular events that modulate aging of the retina would provide insights into the pathogenesis of AMD.
Although a change in the activity of a single gene can be directly correlated to a cellular phenotype (which forms the theoretical basis for linkage analysis of Mendelian traits), the molecular definition of complex phenotypes (e.g., aging and AMD) depends on the interaction of multiple genes and cellular pathways. Microarray technology offers an opportunity to delineate the complex gene expression profiles of specialized cell types.5 6 7 8 9 10 A unique strength of this methodology is that gene discovery is hypothesis independent. We have used commercially produced gene microarrays to examine changes in retinal gene expression, with the goal of identifying genes and cellular events that contribute to retinal aging. In this report, we demonstrate the utility of this technology, report that aging of the human retina is accompanied by discrete changes in gene expression patterns, and present evidence for aging-associated pathways in the retina.
| Materials and Methods |
|---|
|
|
|---|
Microarray Experimentation
Microarray slides (Micromax) were purchased from NEN Life Science Products, Inc. (Boston, MA) and contained 2400 human genes (called probes). Of these, 80% are expressed in the brain. (A complete description of the genes on the microarray slides is available at http://www.umich.edu/
retina/yoshida-iovs.html)
Target Labeling.
Four micrograms of total RNA samples and 350 pg of plant control RNA were converted into either biotin or dinitrophenyl (DNP) cDNA targets using reverse transcriptase (Micromax cDNA Microarray System; NEN Life Science Products; see http://www.nen.com/pdf/MMMANUALPDF.pdf/ for details). To evaluate target labeling, aliquots of labeled cDNAs and a control cDNA of known concentration were serially diluted, individually spotted onto a membrane (GeneScreen; NEN Life Science Products) and fixed by UV cross-linking at 1200 µJ. The membrane filters were developed by chemiluminescence, as described (the enhanced chemiluminescence protocol, NEN Life Science Products). The concentration of test cDNA was determined by identifying the dilution spot with an intensity that most closely matched the least dense intensity spot of the control cDNA, and the cDNA concentrations were calculated by multiplying the target sample dilution factor by the concentration of the control cDNA producing the least dense equivalent spot intensity (data not shown). Labeled cDNAs were used as targets in the hybridization reactions if a 1:800 dilution spot of the target cDNA was visible and the difference between the paired biotin cDNA and the DNP cDNA mass yields was less than fourfold.
Hybridization.
A mixture of equal specific activity of biotin- and DNP-labeled cDNA targets was simultaneously hybridized to the microarray slide at 65°C for 16 hours. After hybridization, slides were washed at room temperature in 0.5x SSC-0.01% SDS for 5 minutes and 0.06x SSC-0.01% SDS for 5 minutes, followed by a 2-minute wash in 0.06x SSC. Biotin and DNP hybridization signals were then detected with Cy-5 and -3 fluorescent reporter molecules, respectively (Tyramide Signal Amplification [TSA protocol]; NEN Life Science Products).
Image Acquisition and Analysis.
Fluorescence intensities were determined by NEN from images taken with The ScanArray 3000 (GSI Lumonics, Inc., Watertown, MA), which was equipped with laser excitation sources and interference filters appropriate for the Cy-3 and -5 fluors. Data were then filtered to remove values from poorly hybridized spots, so that the intensity levels of more than 60% of the pixels in a single spot were greater than the local background.11
The young-to-elderly ratio measurements for each remaining spot on each array were calculated by computer (ScanAlyze software; http://genome-www.stanford.edu/software). A median ratio for each hybridized slide was then calculated based on the ratio of fluors in the young and elderly target populations (see the Results section).8
12
These values were converted to 1, and this conversion factor was used to normalize the results of the hybridization.
Northern Blot Analysis
Total RNA (5 µg) was separated on a 1% formaldehyde-agarose gel and transferred to nylon membrane (Hybond N+; Amersham Pharmacia Biotech, Piscataway, NJ). The cDNA clones for the human KIAA 0120, TGF-ß receptorinteracting protein 1 (TRIP1) and IFN-responsive transcription factor subunit (ISGF3G) were obtained from Research Genetics (Huntsville, AL). Human cDNAs for creatine kinase B (CKB), transferrin (TF) and glyceralaldehyde-3-phosphate dehydrogenase (GAPDH) had been isolated in our laboratory from a subtracted retinal library.13
Purified cDNA inserts were labeled with [
-32P] dCTP using a DNA-labeling system (Multiprime; Amersham Pharmacia Biotech) and hybridized to Northern blot analysis. The blots were normalized against ribosomal RNA, and image intensities were quantified by densitometry of computer (Image Beta 3b software; Scion Corp., Frederick, MD).
Real-Time qRT-PCR
qRT-PCR was performed with gene-specific primers using sequences derived from GenBank (Table 1)
. Purified retinal RNA was prepared from three human samples from young individuals (ages 16, 16, and 18 years) and three from elderly persons (ages 70, 74, and 78 years). These retinal samples were different from those used in the microarray analyses. First-strand cDNA was synthesized with reverse transcriptase (Invitrogen, Carlsbad, CA) and oligo-dT priming. Real-time PCR was performed in triplicate for each individual, with a commercial system (iCycler; Bio-Rad, Hercules, CA) and fluorescence detection (SYBR Green; PE-Applied Biosystems, Foster City, CA), as described herein. The 50-µL reactions were performed in 96-well plates with optical sealing tape (Bio-Rad) and contained the following components: 50 mM Tris (pH 8.3), 3 mM MgCl2, 0.5 mg/mL BSA, 200 µM each dNTP, 0.33x green fluorescence, 0.5 U Taq polymerase (AmpliTaq; PE-Applied Biosystems), two primers at an optimal concentration, and the cDNA template or water. To confirm the specificity of the PCR reactions, products were separated by electrophoresis in 1.0% agarose gels and visualized by ethidium bromide staining. Changes in gene expression were quantified by calculating the average value of the triplicate reactions from each of the three young individuals and comparing that to the average derived from the triplicate analysis of each of the three elderly individuals. This generated a total of nine pair-wise comparisons, from which the mean ± SEM were calculated. Each reaction was normalized against human ß-actin (ACTB).
|
| Results |
|---|
|
|
|---|
After filtration, 52% of the 2400 genes were found to hybridize strongly with the retinal cDNA targets, suggesting that approximately 1250 genes on the microarray slide are expressed in the human retina. Data were further normalized based on the expectation that the expression levels of the majority of retinal genes are not altered during aging (see the Materials and Methods section). (The raw data for all experiments can be accessed at the web site www.umich.edu/
retina/yoshida-iovs.html.) We then defined a level of significance for variations in the expression of individual genes; if a two-fold change in expression (i.e., probes whose Cy3-to-Cy5 ratio is greater than +2 or <0.5) was observed in three of four experiments, a gene was considered to be differentially expressed. Based on this criterion, the level of expression of most of the genes on the chip was constant between the young and elderly retinas. However, a small subset of the genes demonstrated significant and reproducible alterations in mRNA levels.
Our comparative hybridization studies demonstrated that retinal aging is associated with alterations in mRNA levels, which reflect a change in gene expression and/or mRNA stability. The data revealed that 17 (1.4%) genes were expressed at higher levels in retinas from young individuals (referred as young-dominant), whereas 7 (0.6%) genes were designated as elderly-dominant (Table 2) . To establish the validity of the microarray data, we analyzed the patterns of expression of six genes by Northern blot analysis, using the same sources of young and elderly retinal RNAs (Fig. 1) . A comparison of the hybridization intensities indicated that mRNA levels of KIAA0120, TRIP1, and ISGF3G were elevated in retinas from young individuals by 3.5-, 5.0-, and 3.8-fold, respectively (Figs. 1A 1B 1C) , whereas CKB expression in the elderly retinas was nearly twice that in the young retina (Fig. 1D) . Two control genes, TF and GAPDH, showed nearly equivalent levels of gene expression in the young and elderly retinas (Figs. 1E 1F) . In each case, the levels of gene expression on the Northern analyses were consistent with those observed in the microarray experiments.
|
|
To identify potential biological relationships among similarly upregulated genes, we used the recently developed PubGene database (http://www.pubgene.org)14
to identify citation-based gene network associations (Figs. 2A
2B)
. This database includes 10 million MEDLINE records (provided in the public domain from the National Center for Biotechnology Information, Bethesda, MD). We submitted gene symbols of all extracted genes (except ISGF3G and COL7A1) to PubGene and obtained cocitation neighborhoods. The network of eight of the young-dominant genes defined a cocitation neighborhood with somatostatin (SST) at the center, whereas another network of four elderly-dominant genes centered at interleukin-1
(IL1A).
|
| Discussion |
|---|
|
|
|---|
Although microarray technology has the potential to define global patterns of gene expression, the completeness of this catalog is dependent on the arrays that are analyzed. Because we used a commercially produced array containing less than 10% of the genes in the human genome, the analysis was expected to identify only a restricted number of differentially expressed genes. Despite the limitation, this analysis has identified several interesting genes, including a candidate retinal aging gene, CKB. Creatine kinases, including the cytosolic and mitochondrial isoforms, are thought to play a central role in cellular energy metabolism by transporting energy from the site of production in the mitochondrion to that of utilization. Previous studies demonstrated that CKB is expressed in the retina,13 with the highest concentrations in inner segments of the photoreceptors and in the plexiform layers.16 Studies of mouse brain aging have also found evidence for age-related increase in the expression of CKB.7 Higher expression of CKB has been linked to cellular energy stress and may reflect the extent of brain damage.17 These associations are particularly intriguing because expression of several specific stress response genes is also enhanced in the aging retina. Finally, the ability of our unselected array to identify candidate aging gene(s) underscores the power of this technology to generate a hypothesis-independent expression profile of the aging retina.
The analysis of data from microarray studies presents a major challenge,18
and the interpretation of the biological characteristics of genes in each cluster has remained primarily a manual and subjective task. In an attempt to perform a more objective analysis, we used the recently launched PubGene database.14
Although the database is biased toward well-studied genes that are extensively reported in the literature relative to newly discovered genes, it offers a method for rapidly establishing potential associations between genes and functional pathways. Our analysis positions interleukin-1
(IL1A), a potent mediator of inflammation and immunity, in the center of the literature network of elderly-dominant genes (Fig. 2B)
. This is consistent with the accumulating evidence that aging is associated with inflammatory response.7
19
The observation that Somatostatin positions at the center of the young-dominant network (Fig. 2A)
suggests that this neuropeptide may play a more significant role in regulating retinal function than previously envisaged.
Although cDNA microarray technology can provide considerable new insights into gene expression, many aspects require additional development. These include better methods for target labeling, image acquisition, and processing; reduction of intra- and interslide variations, and clustering of gene expression data.10 18 20 We used the TSA system in this analysis, because this was the only available labeling method for small amounts of human RNA when the studies were initiated. Because this is an amplification-based method, the TSA system may introduce some level of bias during label incorporation, reverse transcription, and signal amplification. Our recent studies have demonstrated that another microarray detection system (3DNA Submicro Expression Array Detection System; Genisphere, Hatfield, PA) provides more consistent hybridization data in slide arrays.21 Another possible approach is to complement the slide microarrays with oligonucleotide-based microarrays (e.g., GeneChips from Affymetrix, Santa Clara, CA), in which a two-probe pair strategy is used to minimize cross-hybridization and background signal.22
In addition to the variables associated with the cDNA microarrays, a high sample-to-sample variation is inherent in human tissue samples.23 An optimized method of donor eye preservation24 may reduce such variations. In identical (isogenic) biological systems, the expression of a single gene can fluctuate and exhibit gene-specific patterns of noise.25 To control for both biological and methodological noise, we used multiple sources of young and elderly retinal tissue and exchanged the label (Cy-3 and Cy-5) between the tissue sources. In addition, we increased the significance of our analysis by performing multiple hybridizations.26 This allowed us to generate results with a certain level of significance and to focus only on genes that show larger changes in RNA levels. However, the ability to collect more profiles in parallel may directly influence the extraction of useful data, especially in investigations that use human tissue.27 The variations in gene expression patterns that were observed by using complementary methods with different sets of retinal samples emphasize the importance of analyzing a sufficiently large pool of samples to minimize sample-to-sample variations. Additional investigations with custom gene microarrays of expressed sequence tags (ESTs) generated from retina-RPE libraries13 28 29 30 31 and with an increased number of retinas at various ages are necessary to obtain a comprehensive profile of aging-associated changes in gene expression. Our studies have laid the foundation for future global profiling of retinal and RPE gene expression during development, aging, and disease.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication August 23, 2001; revised March 7, 2002; accepted April, 3, 2002.
Commercial relationships policy: N.
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: Anand Swaroop, Department of Ophthalmology and Visual Sciences, W. K. Kellogg Eye Center, University of Michigan, 1000 Wall Street, Ann Arbor, MI 48105-0714; swaroop{at}umich.edu.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
E. Carrasco, C. Hernandez, A. Miralles, P. Huguet, J. Farres, and R. Simo Lower Somatostatin Expression Is an Early Event in Diabetic Retinopathy and Is Associated With Retinal Neurodegeneration Diabetes Care, November 1, 2007; 30(11): 2902 - 2908. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. P. Gruber, C. D. Coldren, M. D. Woolum, G. P. Cosgrove, C. Zeng, A. E. Baron, M. D. Moore, C. D. Cool, G. S. Worthen, K. K. Brown, et al. Human Lung Project: Evaluating Variance of Gene Expression in the Human Lung Am. J. Respir. Cell Mol. Biol., July 1, 2006; 35(1): 65 - 71. [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] |
||||
![]() |
J. L. Haines, N. Schnetz-Boutaud, S. Schmidt, W. K. Scott, A. Agarwal, E. A. Postel, L. Olson, S. J. Kenealy, M. Hauser, J. R. Gilbert, et al. Functional Candidate Genes in Age-Related Macular Degeneration: Significant Association with VEGF, VLDLR, and LRP6 Invest. Ophthalmol. Vis. Sci., January 1, 2006; 47(1): 329 - 335. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Yu, S. He, J. S. Friedman, M. Akimoto, D. Ghosh, A. J. Mears, D. Hicks, and A. Swaroop Altered Expression of Genes of the Bmp/Smad and Wnt/Calcium Signaling Pathways in the Cone-only Nrl-/- Mouse Retina, Revealed by Gene Profiling Using Custom cDNA Microarrays J. Biol. Chem., October 1, 2004; 279(40): 42211 - 42220. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Zareparsi, A. Hero, D. J. Zack, R. W. Williams, and A. Swaroop Seeing the Unseen: Microarray-Based Gene Expression Profiling in Vision Invest. Ophthalmol. Vis. Sci., August 1, 2004; 45(8): 2457 - 2462. [Full Text] [PDF] |
||||
![]() |
F. Ahmed, K. M. Brown, D. A. Stephan, J. C. Morrison, E. C. Johnson, and S. I. Tomarev Microarray Analysis of Changes in mRNA Levels in the Rat Retina after Experimental Elevation of Intraocular Pressure Invest. Ophthalmol. Vis. Sci., April 1, 2004; 45(4): 1247 - 1258. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. I. Dorrell, E. Aguilar, C. Weber, and M. Friedlander Global Gene Expression Analysis of the Developing Postnatal Mouse Retina Invest. Ophthalmol. Vis. Sci., March 1, 2004; 45(3): 1009 - 1019. [Abstract] [Full Text] [PDF] |
||||
![]() |
I. Chowers, D. Liu, R. H. Farkas, T. L. Gunatilaka, A. S. Hackam, S. L. Bernstein, P. A. Campochiaro, G. Parmigiani, and D. J. Zack Gene expression variation in the adult human retina Hum. Mol. Genet., November 15, 2003; 12(22): 2881 - 2893. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. Nielsen, K. Birkenkamp-Demtroder, N. Ehlers, and T. F. Orntoft Identification of Differentially Expressed Genes in Keratoconus Epithelium Analyzed on Microarrays Invest. Ophthalmol. Vis. Sci., June 1, 2003; 44(6): 2466 - 2476. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |