(Investigative Ophthalmology and Visual Science. 2004;45:2503-2513.)
© 2004 by The Association for Research in Vision and Ophthalmology, Inc.
DOI: 10.1167/iovs.03-1391
Gene Expression Profiling of Purified Rat Retinal Ganglion Cells
Ronald H. Farkas,1,2
Jiang Qian,1,2
Jeffrey L. Goldberg,3,4
Harry A. Quigley,1,5 and
Donald J. Zack1,2,6,7,8
1From the Department of Ophthalmology, Wilmer Eye Institute, the
2Guerrieri Center for Genetic Engineering and Molecular Ophthalmology, and the
5Dana Center for Preventive Ophthalmology at the Wilmer Eye Institute, Baltimore, Maryland; the Departments of
6Molecular Biology and Genetics,
7Neuroscience, and the
8McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland; and the
3Department of Neurobiology, Stanford University, Stanford, California.
4Present affiliation: Bascom Palmer Eye Institute, McKnight Vision Research Center, Miami, Florida.
 |
Abstract
|
|---|
PURPOSE. The phenotype of specialized cells arises, in part, from their characteristic gene expression patterns. Retinal ganglion cells (RGCs) are of wide interest in neuroscience and die in glaucoma and other optic neuropathies. In this study the genes expressed by RGCs were profiled by expressed sequence tag (EST) analysis.
METHODS. ESTs were generated from a cDNA library constructed from RGCs isolated by immunopanning. The RGC genes were compared with published microarray expression profiles from 13 different neural regions. Immunohistochemistry was performed by standard methods.
RESULTS. Clustering of 4791 RGC ESTs identified 2360 unique gene clusters. Of these, 60% represented known genes, 27% uncharacterized genes/ESTs, and 13% novel sequence. Unexpectedly, one of the largest RGC clusters, RESP18, corresponded to a neuroendocrine-specific gene preferentially expressed in the hypothalamus. RESP18 immunoreactivity within the retina was found mainly in the RGC layer. DDAH1, a gene involved in nitric oxide metabolism, was localized to RGC and amacrine layers. Comparison of gene expression patterns across neuronal regions revealed a prominent subset of RGC genes that were overexpressed in dorsal root and trigeminal ganglia. To narrow the search for candidate disease-related genes, RGC genes were mapped to known disease loci for optic neuropathies.
CONCLUSIONS. This work is one of the first efforts to profile gene expression in a purified population of retinal neurons, the RGCs. The profiling, in addition to revealing both known and novel genes underlying the RGC phenotype, also uncovered common patterns of gene expression between RGCs and other sensory ganglia.
Retinal ganglion cells (RGCs) perform the final steps of retinal vision processing before carrying the visual signal to higher centers in the brain. The death of RGCs is a common cause of visual disability, most often occurring from glaucoma, a leading cause of irreversible blindness in developed countries and worldwide.1 2 RGC death also occurs in many other optic neuropathies, which can be of environmental and/or genetic origin, with inheritance patterns ranging from mitochondrial to autosomal dominant. RGCs are widely used as a model system for neurobiological studies of development, trauma, and neuronal degeneration and regeneration. Thus, knowledge of RGC functions is important to the understanding of many neuronal processes.
The behavior of RGCs, and of all cells on a fundamental level, arises in large part from the specific distribution of genes that the cells express. Relatively little is known about the genes expressed by RGCs. Gene expression studies have been performed on the retina as a whole,3 using a variety of methods, including expressed sequence tag (EST) sequencing,4 5 6 serial analysis of gene expression (SAGE),7 8 and microarrays.9 10 11 12 13 However, the retina is a complex tissue composed of neuronal, glial, and vascular cell types. The dominant cell type of the retina is the photoreceptor, with RGCs comprising 5% or less of all retinal cells.14 15 Thus, gene expression profiles from whole retina are unrepresentative of RGC gene expression. Previous efforts to identify the complete picture of retinal gene expression may have missed important RGC genes, due to the low relative frequency of RGCs and, hence, their mRNA transcripts.
Efforts to document gene expression of specific cell types in the eye have used EST sequencing from libraries derived separately from lens, iris, and retinal pigment epithelium-choroid.4 Similarly, optic nerve astrocytes isolated from human retina have been studied with microarrays.16 Efforts are currently underway to identify the expression pattern of single retinal cells.17 Gene expression profiling has been undertaken for different subregions of the brain18 19 20 21 22 23 24 25 and auditory system,26 with the use of message amplification for cells microdissected individually20 22 27 28 or purified by fluorescence-activated cell sorting.29
We describe herein an initial gene expression profile of rat RGCs that were separated from other cell types by immunopanning. RGCs purified by this method have provided information on the interactions that govern RGC physiology and development.30 Through EST sequencing of a cDNA library prepared from purified RGCs, we found expression of both known and novel genes and have begun further investigation of how these genes relate to RGC function.
 |
Materials and Methods
|
|---|
RGC cDNA Library
All experiments were conducted with adherence to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. RGCs from postnatal day (P)21 Sprague-Dawley rats (Simonsen Laboratories, Gilroy, CA) were purified by sequential immunopanning using a monoclonal antibody to Thy1.1, by a method that has been reported to yield approximately 99% purity.30 31 33 Briefly, whole retina tissue was removed from rats from four litters, incubated in a papain solution, dissociated into single-cell suspension, and subjected to negative selection with anti-macrophage antiserum (1:100; Axell Accurate Chemical and Scientific Corp, Westbury, NY) to remove macrophage and endothelial cells. Remaining cells were then placed on panning plates containing Thy1.1 antibody, and unbound cells were removed with washing. RGCs were released after incubation with a trypsin solution. RGCs were cultured overnight on poly-D-lysine (Sigma-Aldrich) and laminin (Invitrogen-Gibco, Grand Island, NY), in serum-free defined medium, as described.31 32 The yield of RGC cells was between 50% and 70%, with a survival of more than 90%.31 Total RNA was extracted and confirmed to be intact by denaturing gel electrophoresis. The total RNA was reverse transcribed, and a nonnormalized, nonsubtracted cDNA library was constructed from 1 µg total RNA using SMART technology and 16 rounds of PCR amplification (BD Biosciences-Clontech, Palo Alto, CA). This method is based on template-switching at the 5' end of the mRNA transcript, with binding of an anchor sequence that can serve as a priming site for PCR. The cDNA was then size fractionated and ligated into pDNR-LIB plasmid (BD Biosciences-Clontech). Escherichia coli (Electromax DH10B; Invitrogen, San Diego, CA) was transformed by electroporation and grown for 1 hour in SOC medium without further amplification. Clones were plated and randomly chosen for 5' sequencing (Michigan State University Genomics Facility, East Lansing, MI, or Integrated Genomics GmbH, Jena, Germany). RGC EST sequences were submitted to the National Center for Biotechnology Information (NCBI; accession numbers CF974889 to CF9749532).
EST Sequence Clustering
Our clustering procedure consisted of two steps. First, the CAP3 sequence assembly program34 was used with stringent criteria: overlap length of 200 bp and sequence identity percentage of 95% (command line options: -o 200 -p 95). Then, the sequence clusters resulting from the CAP3 program were aligned against the rat, mouse, and human Reference Sequence Project (RefSeq) databases35 36 with BLASTn.37 The e-value cutoff for the alignment was defined as 1e-50. The clusters aligning to the same RefSeq sequence were further merged into one cluster. The remaining clusters were then aligned against the rat, mouse, and human UniGene38 and genomic databases with BLASTn at an e-value defined as 1e-20 and then further aligned with BLASTx alignment at an e-value cutoff of 1e-5.
Cluster Annotation
All clusters were annotated to RefSeq preferentially35 36 and otherwise to UniGene. Rat, mouse, and human RefSeq were used as complementary annotation resources. If a cluster aligned to a RefSeq from a different species than rat, the RefSeq with lowest e-value alignment was selected as the annotation for the cluster. In cases in which the alignments between a cluster and a RefSeq from different species yielded the same e-value (e.g., e-value = 0.0), the priority order of annotation was rat>mouse>human.
dbEST Search
The sequences without a RefSeq or UniGene match were further aligned against rat, mouse, and human database for ESTs (dbEST)39 with the same e-value cutoff (1e-20) applied.
Comparative Gene Expression Clustering
The RGC genes used for hierarchical clustering composed the subset of RGC clusters represented in the microarray data compiled by Su et al.19 (http://expression.gnf.org/cgi-bin/index.cgi, Gene Expression Atlas, Genomics Institute of the Novartis Research Foundation, San Diego, CA) from 13 different neural regions. Replicate array data available for each region were first averaged. The Cluster and TreeView programs developed by Michael Eisen (http://rana.lbl.gov/EisenSoftware.htm; provided in the public domain by Eisen Lab, Lawrence Berkeley, national Lab, University of California at Berkeley) were then used for clustering and visualization, respectively. Both genes and neural regions were normalized by median centering and filtered for genes with expression in at least one neural region. Negative expression levels were set to zero. The parameters used for clustering were "absolute correlation" and "average linkage clustering."
Disease Loci
All sequence clusters were mapped to the human genome with an mRNA alignment program, SPIDEY.40 The program was run with the interspecies alignment option on. The clusters falling into the disease loci were identified as disease candidate genes.
Immunohistochemistry
Eight-micrometer paraffin sections of rat retina were collected onto slides (Superfrost Plus; Fisher Scientific, Pittsburgh, PA) before immunolabeling by the streptavidin-biotin peroxidase method (Vector-SG detection kit; Vector Laboratories Inc., Burlingame, CA). Primary antibody was to RESP18 (a gift from Betty Eipper, University of Connecticut Health Center, Farmington, CT) or DDAH1 (a gift of Milan Va
ák, Department of Biochemistry, Universität Zürich, Zürich, Switzerland). Negative controls included nonimmune serum of the same species as the primary antibody, or preimmune serum at the same protein concentration, or incubation buffer alone. Labeled sections were mounted in glycerol jelly and viewed by Nomarski optics (Axioskop; Carl Zeiss Meditec, Thornwood, NY).
 |
Results
|
|---|
RGC Library Clones: Known and Novel
A nonnormalized, nonsubtracted cDNA library was constructed from purified postnatal day (P)21 rat RGCs. Generation of the library incorporated 16 rounds of PCR amplification to compensate for the limiting quantities of RNA available. The library contained 1.2 x 105 primary transcripts, with an average insert size of 1.8 kb. No further amplification of the library was performed. Five thousand clones were randomly selected for 5' sequencing, generating 4791 usable EST sequences, with an average usable read length of approximately 500 bp.
Clustering of these ESTs to group together sequences representing the same transcript revealed 2360 unique gene clusters (complete listing available in Supplemental Table 1 at www.iovs.org/cgi/content/full/45/8/2503/DC1). Approximately 40% of the RGC library clones could be characterized as poorly characterized sequence, commonly termed ESTs. In a more detailed analysis, approximately 60% of the clusters (1408) represented known genes in the sense that they matched genes categorized in the NCBI Reference Sequence Project (RefSeq) sets for rat, mouse, and human genes. For the remaining 952 clusters a less stringent homology search was performed using UniGene for rat, mouse, and human. Approximately 60% of these clusters (n = 565) matched to UniGene clusters. Of the other 301 clusters, 92 matched to previously submitted EST sequences in dbEST. For each of the RGC clusters that matched only to otherwise unclassified dbEST ESTs, the tissue of origin of the matching dbEST ESTs was identified (Supplemental Table 1), to serve as a rough indication of the tissue types that express the gene.
The next step in the partitioning of RGC genes involved the 209 RGC clusters that had no match in the EST database by BLASTn. To characterize these remaining clusters further, and to select those most likely to represent rarely transcribed sequences, we aligned the clusters against the rat genome using BLAST. This alignment demonstrated that 82 of the clusters possessed two or more discontinuous areas of alignment on the genome, as would occur for exons of transcribed genes. In agreement with this, most of these sequences, 73%, could be readily amplified by PCR. Of these 82 clusters with intronexon structure, 74 showed no homology in the protein database or matched to predictedtheoretical protein records.
Library Representation of RGC Gene Expression
The distribution of RGC gene cluster size is shown in Figure 1 . Most (77%) of the ESTs had no significant sequence overlap with other RGC ESTs,and so were in a cluster of size 1, as reflects the fact that most cellular mRNAs are expressed at low levels. As cluster size increased, the number of clusters at each increasing size was fewer, in an exponential-like decline, reflecting the fact that a smaller percentage of cellular genes are expressed at moderate levels, with fewer still expressed very highly. Though genes that are very highly expressed are few in number, they still account for a relatively large percentage of total mRNA produced by a cell, with the 10 largest RGC gene clusters containing approximately 500 of the 4791 total ESTs sequenced, or nearly 10% of the total (Table 1) . We observed a strong inverse correlation between cluster size and the percentage of clusters that represented an unknown EST. For singleton clusters, as mentioned earlier, approximately 40% were ESTs. As cluster size increased, the percentage of ESTs decreased, from 32% for cluster size 2, to 16%, 8%, 6%, and 0% for cluster sizes between 3 and 6 or more, respectively.

View larger version (12K):
[in this window]
[in a new window]
|
FIGURE 1. Cluster analysis for 4791 RGC ESTs. The histogram groups the EST clusters by size, from clusters made of only 1 EST sequence, to clusters containing 34 or more overlapping ESTs that collectively represent a single expressed sequence. Most (77%) of the ESTs had no significant sequence overlap with other RGC ESTs and so were in a cluster of size 1. As cluster size increased, the number of clusters at a given size decreased exponentially, reflecting the fact that fewer cellular genes are expressed at moderate and high levels.
|
|
Highly Expressed RGC Genes
Gene function assignments were made for the 100 most frequently sequenced RGC genes (Fig. 2) . For this purpose, each gene was assigned to a single main category. Despite being highly represented, some genes have poorly understood functions. Nevertheless, functional assignments were possible for most genes in this group.

View larger version (28K):
[in this window]
[in a new window]
|
FIGURE 2. Gene ontogeny of the 100 most frequently sequenced RGC clones. Each gene was assigned to a single category corresponding to the main functional role of the gene, as determined by literature review.
|
|
The protein modification enzyme ubiquitin was the single most frequently sequenced RGC gene, second only to the combination of mitochondria-encoded proteins. Of the 13 proteins encoded by the mitochondrial genome, 12 were sequenced from the RGC library (Table 2) . The number of ESTs encountered for each of these genes varied greatly, from 18 for the largest cluster (ATP synthase F0 subunit 6) to NADH dehydrogenase subunit 3, which was not encountered at all. Mitochondrial mutations lead to selective RGC death in Lebers optic neuropathy, suggesting that mitochondrial gene expression levels may be particularly important in ganglion cell function. We therefore compared the cluster size found in the RGC library for each mitochondria-encoded gene to the cluster size in a non-normalized, nonsubtracted cDNA library (dbEST library ID.8847) from rat dorsal root ganglia, a cell type with similarities to RGCs that does not significantly degenerate in Lebers optic neuropathy. Although our sample size was small, there was a good overall correlation between cluster sizes for the mitochondria-encoded genes in the two cell types (R2 = 0.4; Fig. 3 ), and thus no evidence of large differences in expression levels of mitochondria-encoded genes between RGC and DRG cells.

View larger version (11K):
[in this window]
[in a new window]
|
FIGURE 3. Scatterplot showing similar expression levels (as represented by cluster size) of mitochondria-encoded genes in the RGC library compared with a library derived from dorsal root ganglion (DRG) cells (dbEST Library ID.8847,41 ). The 13 mitochondria-encoded genes were assigned rank-order values (113) based on the number of ESTs/cluster in each library (R2 = 0.4).
|
|
Also highly expressed in RGC cells were glycolytic and citric acid cycle enzymes, including malate dehydrogenase, lactate dehydrogenase, aldolase C, and aldolase A. These proteins are also highly represented in SAGE results from brain and retina,7 8 and represent large clusters in UniGene, reflecting broad expression in many different tissues. Calcium-binding proteins were frequently encountered in the library, including calmodulin-1, -2, and -3; calretinin; calpactin 1; and calcyclin. Antioxidant genes were also prominent, including peroxiredoxin, glutathione S-transferase, and the poorly characterized selenoprotein W. After the five most represented functional categories (energy, protein synthesis/modification, lysosomal/vesicular, synaptic function, cytostructure), the next largest group was ESTs without known function.
Photoreceptor-specific genes are among the most frequently encountered ESTs in cDNA libraries made from whole retina, whereas in contrast in the RGC ESTs they were almost entirely excluded. For example, no ESTs were encountered for rhodopsin, transducin, or rod photoreceptor cGMP phosphodiesterase subunits. Similarly, no clones were encountered in the RGC library for Müller cell genes, such as cellular retinaldehyde-binding protein (CRALBP), vimentin, or glutamine synthase.
Neural- and Disease-Associated Gene Expression in RGCs
Surprisingly, a neuroendocrine-specific gene with poorly understood function, RESP18, was one of the most frequently encountered genes in the RGC library. We used immunohistochemistry to correlate RESP18 message found in the library with protein expression in the retina, and to localize its cellular expression pattern. RESP18 immunoreactivity was observed mainly in the RGC layer, with faint labeling of the inner plexiform layer and inner nuclear layer (Fig. 4) . Staining was present mainly in cell bodies, in a perinuclear distribution consistent with its previously identified localization to the lumen of the endoplasmic reticulum.42

View larger version (156K):
[in this window]
[in a new window]
|
FIGURE 4. Expression of RESP18 in rat retina. RESP18 immunoreactivity was present mainly in large cell bodies in the ganglion cell layer, with faint immunoreactivity in the inner plexiform and inner nuclear layers and negligible labeling elsewhere. Immunoperoxidase staining with brown 3-amino-9-ethylcarbazole (AEC) reaction product. INL, inner nuclear layer; IPL, inner plexiform layer; RGC, retinal ganglion cell layer. Bar, 15 µm.
|
|
The functions of the remaining RGC clusters were then examined for genes that might be involved in pathways particularly important to RGC disease. One of the singleton genes, DDAH1, has been investigated for its involvement in nitric oxide metabolism43 44 45 46 and neuronal injury. We found immunoreactivity for DDAH1 to be present in RGC and amacrine cell layers (Fig. 5) . The sublaminae of the inner plexiform layer were particularly intensely immunoreactive, as was the retinal pigment epithelium.

View larger version (183K):
[in this window]
[in a new window]
|
FIGURE 5. Expression of DDAH1 in rat retina. DDAH1 immunoreactivity was present in the ganglion cell layer, multiple sublaminae of the inner plexiform layer, and amacrine cell regions of the inner nuclear layer. Strong immunoreactivity was also present in the retinal pigment epithelium. Immunoperoxidase labeling with brown 3-amino-9-ethylcarbazole (AEC) reaction product. INL, inner nuclear layer; IPL, inner plexiform layer; RGC, retinal ganglion cell layer. Bar, 15 µm.
|
|
Comparative Gene Expression Profiling of RGC Genes
Just as knowing which genes are expressed by RGCs should be informative of cell function, knowing the different neuronal types that express a gene should be informative of gene function. We therefore sought to determine where else in the nervous system genes expressed in RGCs were also expressed. Using hierarchical clustering to visualize expression patterns (Fig. 6) , we compared our RGC data set to gene expression profiles generated using oligonucleotide microarrays from 13 separate neuronal regions19 (http://expression.gnf.org/cgi-bin/index.cgi). Of the 2360 RGC gene clusters, 959 were represented on these arrays, with 89% (851/959) showing positive signal in at least one neural area. A major resultant feature of the clustering was a subgroup of 24 genes relatively overexpressed in both dorsal root and trigeminal ganglia (Fig. 6 , large oval). These sensory ganglia overexpressed genes, along with numerical expression level data are shown in Table 3 . Hierarchical clustering was also performed across neural regions, as indicated by the tree structure. For the subset of genes expressed in retinal ganglion cells, the dorsal root and trigeminal ganglia were more closely related to each other than to the other neural regions.

View larger version (18K):
[in this window]
[in a new window]
|
FIGURE 6. Comparative gene expression profiles of RGC genes. Comparison was made between the gene expression profile derived from RGCs through EST sequencing versus the gene expression profiles of 13 separate neuronal regions derived by Su et al.18 using oligonucleotide microarrays. Hierarchical clustering is shown of the microarray output from the resultant 851 genes present in both ganglion cell and microarray data sets. Mean-centered normalization was applied across rows, so that red represents expression higher than average across the neural regions and green represents below average expression level. Overall brightness represents microarray signal intensity. Genes with similar expression patterns across neuronal areas are placed in proximity to each other. Genes relatively overexpressed in both dorsal root and trigeminal ganglia are indicated by the large oval and genes relatively overexpressed in hypothalamus are indicated by the small oval.
|
|
To further explore the prominent RGC expression of RESP18, we focused on the smaller subset of RGC genes shown through clustering to be overexpressed in hypothalamus (Fig. 6 small oval, Table 4 ). One of these genes, necdin, has been reported to be overexpressed in hypothalamus,47 48 whereas another, secretogranin 2, has known neuroendocrine functions. Necdin is thought to play a role in cell cycle arrest in terminally differentiated neurons, and secretogranin 2 is one of the chromogranin family of proteins, involved in secretory granulogenesis and the sorting and processing of secretory proteins.
Candidate Optic Neuropathy Genes Expressed by RGCs
To aid in candidate gene approaches to finding optic neuropathy-associated genes, we mapped the RGC genes to syntenic regions of the human genome associated with the optic neuropathyassociated loci GLC1B,49 GLC1C,50 GLC1D,51 GLC1F,52 OPA2,53 OPA4,54 WFS2,55 SCABD,56 and ARTS57 (Table 5) .
 |
Discussion
|
|---|
We have established an initial gene expression profile of purified rat RGCs through the approach of EST sequencing. This has allowed us to ascribe to RGCs the expression of a large number of known and novel genes. Several hundred novel sequences were identified from our sampling of only a minority of total genes expressed by RGCs, suggesting that we still have much to discover about RGC gene expression, especially about those genes expressed at low levels.
Although novel sequences from cDNA libraries are often described as ESTs, the term has certain limitations. For example, it does not distinguish rare novel sequences from more common sequences that are as yet unnamed. We therefore undertook a more detailed analysis of RGC gene clusters based on the full information available in NCBI databases. The RefSeq collections represent stable, curated databases of genes that have been documented at a higher level of certainty than those found in other experimental gene collections such as UniGene. The RGC genes that matched RefSeq entries were therefore considered known genes for this analysis. It is important to note, however, that most genes in RefSeq remain largely uncharacterized and that some characterized genes have not, at present, been included in RefSeq. The UniGene sets are based on automated partitioning of EST data into an experimental set of unique genes, based on sequence alignment combined with evidence of polyadenylation. Thus, in general, the RGC genes that matched to UniGene, but not to RefSeq, represented genes that are poorly characterized, possibly because of low expression levels. Those RGC sequences so rare as not to be found even in UniGene may represent transcripts of particular interest for future study.
The RGC clones were randomly selected and sequenced from a non-normalized, nonsubtracted cDNA library. This methodology can provide approximate information about the relative expression level of genes, but must be considered with caution in this instance because PCR was used in library construction to compensate for the limited quantity of RGC RNA. Several observations suggest that the rank ordering of genes we obtained may still be a useful representation of RGC gene expression levels: (1) In general, within any population of total mRNA, a small number of genes are highly abundant, with each of their individual mRNAs accounting for up to 1% of total cellular mRNA transcripts. For example, rhodopsin transcripts account for 1% or more of total retinal mRNA, and it is often the case that the 10 most highly expressed genes in a tissue account for approximately 10% of total mRNA. Such a distribution of highly expressed genes was encountered in this library, and included polyubiquitin and tubulin, among other genes expected to be highly expressed. (2) Along these same lines, nearly all genes in a tissue are of medium (0.1%1%) or low abundance (0.1% or less of total mRNA).58 59 Single EST sequences represented 77% of clusters in the RGC library, and the overall cluster size distribution was similar to that of a non-PCR amplified retinal cDNA library11 as would be expected if the RGC library also contained a sampling of low-, medium-, and high-abundance mRNA species. (3) Mitochondria-encoded genes were expressed in the RGC library at ratios similar to levels in a non-PCR amplified library from dorsal root ganglion cells (Fig. 4) . Lastly, (4) a strong correlation was present between RGC cluster size and the likelihood of the clusters representing a known gene, as would be expected if highly expressed genes are more likely to be known, whereas low expression level genes are more likely to be unstudied ESTs.
Within the 35 genes most frequently sequenced from the RGC library (Table 1) , a surprising number have poorly understood functions, and/or have not been studied in the retina. One of these for which there is some interesting information available is RESP18, which is expressed in endocrine and neural cells that secrete peptide hormones60 61 and is suggested to be involved in intracellular signaling pathways.62 In the brain, expression of RESP18 is limited to very specific regions, including the paraventricular and supraoptic nuclei of the hypothalamus, and the neurointermediate and anterior lobes of the pituitary. RESP18 is homologous to a region of the luminal domain of neuroendocrine-specific receptor-type protein tyrosine phosphatases and appears to be involved in intracellular signaling from the Golgi to the nucleus.62 We speculate that further study of RESP18 may uncover a group of functionally interrelated genes that act within similar pathways in RGCs and the small subset of specialized neurons that preferentially express this gene.
Most mammalian genes remain poorly characterized, and similarly most of the genes we identified as singletons in the RGC library remain poorly characterized, both in general and in the retina. We selected DDAH1 for further study based on its involvement in nitric oxide metabolism, which has been particularly implicated in glaucomatous RGC death.63 64 65 66 Also of interest, DDAH1 is increased in neurons with neurofibrillary tangles in Alzheimer disease67 and in neuronal trauma.68 DDAH1 hydrolyzes arginine residues that are asymmetrically methylated and endogenously generated. Such residues are competitive inhibitors of all three isoforms of nitric oxide synthase (NOS). Thus, inhibiting DDAH1 leads to decreased NO synthesis. This action may provide a therapeutic avenue for modulating NO levels in disease. The observation that DDAH1 is present both in the RGC and amacrine cell layers, as well as in multiple sublaminae of the inner plexiform layer, suggests a role for DDAH1 in ganglion and amacrine cell signaling.
Several of the genes most frequently sequenced from the RGC library have been found to change expression in RGC injury. Synuclein gamma is increased in retinal glial cells in human glaucoma.69 In an optic nerve crush model of ganglion cell death, stathmin showed reduced expression, whereas malate dehydrogenase and heat shock protein cognate 70 were both increased.70
Human genetic studies have identified several loci for inherited optic neuropathies. Candidate gene approaches to finding the specific gene causing these RGC diseases will be aided by knowledge of which genes in mapped loci are in fact expressed by RGCs (recognizing that a gene causing RGC death need not be expressed by RGCs). Of the RGC-expressed genes in these regions (Table 5) , those known to be involved in neuronal degenerations or that are specifically expressed in neurons may be stronger candidates for RGC disease. Thus, the OPA4 region contained an "EST similar to ELAV (Embryonic lethal abnormal vision)-type RNA-binding protein 3" as annotated by NCBI. The ELAV family of genes is essential for Drosophila visual system development, and autoimmunity to these proteins can cause encephalomyelitis with sensory neuropathy.71 Also within the OPA4 region was Ras-like without CAAX/Rin/Rit2, a small guanosine triphosphate (GTP)binding protein that is expressed exclusively in neurons and in the retina is found only in ganglion cells and in the inner nuclear layer.72 73 The GLC1D region similarly contained a neuron-restricted gene, potassium channel alpha subunit (Kv8.1). This gene has been analyzed in a positional candidate approach for adult familial myoclonic epilepsy, but without findings of pathogenic changes.74
One way to understand how the genome influences cell function is to compare the gene expression profiles among different cell types. This approach is amenable to high-throughput methods, such as EST sequencing, microarrays, and SAGE, and is particularly useful for gaining an initial insight into the roles of the large number of mammalian genes for which no detailed studies have been undertaken. In general, genes that are expressed at high levels in one neuronal type, but not in others, are likely contributors to phenotypic differences. In contrast, genes expressed at similar levels across neuronal subtypes might be expected to perform general cellular or neuronal functions and not to have a large role in establishing phenotypic differences. We examined the expression profile of the RGC genes in 13 different neural regions, as determined by the oligonucleotide microarray analysis of Su et al.19 It is known that RGCs and other sensory ganglia such as the dorsal root and trigeminal ganglia share transcription factors such as the Brn-3 family.75 As a consequence, it has been proposed that these neurons share similar genetic regulatory hierarchies.75 Consequently, it might be expected that the cells would show downstream gene expression patterns with some significant similarities. This, in fact, appeared to be the case, with hierarchical clustering of gene expression revealing a pattern of RGC genes that were preferentially expressed in both the dorsal root and trigeminal ganglia. Several of these genes were among the largest RGC clusters, including tubulin, synuclein gamma, stathmin-like 2, and HSC70. Ion channels were also prominent in this group, including chloride intracellular channel 1, annexin 5, calpactin 1, and aquaporin 1.
There are many RGC subtypes, some of which have properties such as photoreception that set them apart from others.76 77 78 This suggests that diversity will be found in gene expression patterns, even among RGCs, and highlights the need for increasingly stringent cell-type separations for future gene expression studies. The method of RGC purification used in the present study did not specifically bias against particular RGC subtypes, and the yield and survival of RGCs suggests that a fairly representative sample of RGCs was obtained. As a consequence, those subtypes present in small numbers, such as photoreceptive RGCs, would have made little contribution to the overall gene expression profile obtained. It is important to consider that gene expression patterns can vary greatly between seemingly similar samples, even controlling for such factors as genetic background, age, sex, and environment.79 The RGC population examined in the current study was collected at a single time point in the midmorning, and the expression of some genes may be altered at different circadian time points, for example. Similarly, the specific conditions during RGC purification and overnight culture probably affected the expression levels of some genes. Thus, instead of a single, static gene expression profile defining a cell type, it appears that a range of profiles, when combined, better illustrates the more dynamic nature of transcriptional control of cellular phenotype.
 |
Acknowledgements
|
|---|
The authors thank Danielle Valenta and Mary Ellen Pease for technical assistance.
 |
Footnotes
|
|---|
Supported by grants from the National Eye Institute, Glaucoma Foundation, Glaucoma Research Foundation, the American Health Assistance Foundation, and Santen Pharmaceutical, and by generous gifts from Mr. and Mrs. Marshall and Stevie Wishnack and from Mr. and Mrs. Robert and Clarice Smith. HAQ is the A. Edward Maumenee Professor of Ophthalmology; DJZ is the Guerrieri Professor of Genetic Engineering and Molecular Ophthalmology and is a recipient of a Research to Prevent Blindness Senior Investigator Award.
Submitted for publication December 23, 2003; revised March 2, 2004; accepted March 12, 2004.
Disclosure: R.H. Farkas, None; J. Qian, None; J.L. Goldberg, None; H.A. Quigley, None; D.J. Zack, 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: Ronald H. Farkas, Department of Ophthalmology, Wilmer Eye Institute, Maumenee 813, Johns Hopkins University School of Medicine, Baltimore, MD 21287; rfarkas{at}jhmi.edu.
 |
References
|
|---|
- Sommer A, Tielsch JM, Katz J, et al. Racial differences in the cause-specific prevalence of blindness in east Baltimore. N Engl J Med. 1991;325:14121417.[Abstract]
- Quigley HA. Number of people with glaucoma worldwide. Br J Ophthalmol. 1996;80:389393.[Abstract/Free Full Text]
- Swaroop A, Zack DJ. Transcriptome analysis of the retina. Genome Biol. 2002;3:REVIEWS1022.[Medline][Order article via Infotrieve]
- Wistow G. A project for ocular bioinformatics: NEIBank. Mol Vis. 2002;8:161163.[ISI][Medline][Order article via Infotrieve]
- Sinha S, Sharma A, Agarwal N, Swaroop A, Yang-Feng TL. Expression profile and chromosomal location of cDNA clones, identified from an enriched adult retina library. Invest Ophthalmol Vis Sci. 2000;41:2428.[Abstract/Free Full Text]
- Katsanis N, Worley KC, Gonzalez G, Ansley SJ, Lupski JR. A computational/functional genomics approach for the enrichment of the retinal transcriptome and the identification of positional candidate retinopathy genes. Proc Natl Acad Sci USA. 2002;99:1432614331.[Abstract/Free Full Text]
- Sharon D, Blackshaw S, Cepko CL, Dryja TP. Profile of the genes expressed in the human peripheral retina, macula, and retinal pigment epithelium determined through serial analysis of gene expression (SAGE). Proc Natl Acad Sci USA. 2002;99:315320.[Abstract/Free Full Text]
- Blackshaw S, Fraioli RE, Furukawa T, Cepko CL. Comprehensive analysis of photoreceptor gene expression and the identification of candidate retinal disease genes. Cell. 2001;107:579589.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Chowers I, Liu D, Farkas RH, et al. Gene expression variation in the adult human retina. Hum Mol Genet. 2003;12:28812893.[Abstract/Free Full Text]
- Chowers I, Gunatilaka TL, Farkas RH, et al. Identification of novel genes preferentially expressed in the retina using a custom human retina cDNA microarray. Invest Ophthalmol Vis Sci. 2003;44:37323741.[Abstract/Free Full Text]
- Mu X, Zhao S, Pershad R, et al. Gene expression in the developing mouse retina by EST sequencing and microarray analysis. Nucleic Acids Res. 2001;29:49834993.[Abstract/Free Full Text]
- Livesey FJ, Furukawa T, Steffen MA, Church GM, Cepko CL. Microarray analysis of the transcriptional network controlled by the photoreceptor homeobox gene Crx. Curr Biol. 2000;10:301310.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Farjo R, Yu J, Othman MI, et al. Mouse eye gene microarrays for investigating ocular development and disease. Vision Res. 2002;42:463470.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Jeon CJ, Strettoi E, Masland RH. The major cell populations of the mouse retina. J Neurosci. 1998;18:89368946.[Abstract/Free Full Text]
- Ahmad KM, Klug K, Herr S, Sterling P, Schein S. Cell density ratios in a foveal patch in macaque retina. Vis Neurosci. 2003;20:189209.[CrossRef][ISI][Medline][Order article via Infotrieve]
- 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:4564.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Wahlin KJ, Grice EA, Hackam AS, et al. A method for analysis of gene expression in isolated mouse photoreceptor and Müller cells. Mol Vis. 2004;10:366375.[ISI][Medline][Order article via Infotrieve]
- Miki R, Kadota K, Bono H, et al. Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays. Proc Natl Acad Sci USA. 2001;98:21992204.[Abstract/Free Full Text]
- Su AI, Cooke MP, Ching KA, et al. Large-scale analysis of the human and mouse transcriptomes. Proc Natl Acad Sci USA. 2002;99:44654470.[Abstract/Free Full Text]
- Bonaventure P, Guo H, Tian B, et al. Nuclei and subnuclei gene expression profiling in mammalian brain. Brain Res. 2002;943:3847.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Zhao X, Lein ES, He A, Smith SC, Aston C, Gage FH. Transcriptional profiling reveals strict boundaries between hippocampal subregions. J Comp Neurol. 2001;441:187196.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Che S, Ginsberg SD. Amplification of RNA transcripts using terminal continuation. Lab Invest. 2004;84:131137.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Pavlidis P, Noble WS. Analysis of strain and regional variation in gene expression in mouse brain. Genome Biol. 2001;2:RESEARCH0042.[Medline][Order article via Infotrieve]
- Zirlinger M, Kreiman G, Anderson DJ. Amygdala-enriched genes identified by microarray technology are restricted to specific amygdaloid subnuclei. Proc Natl Acad Sci USA. 2001;98:52705275.[Abstract/Free Full Text]
- Lee CK, Weindruch R, Prolla TA. Gene-expression profile of the ageing brain in mice. Nat Genet. 2000;25:294297.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Cho Y, Gong TW, Stover T, Lomax MI, Altschuler RA. Gene expression profiles of the rat cochlea, cochlear nucleus, and inferior colliculus. J Assoc Res Otolaryngol. 2002;3:5467.[CrossRef][Medline][Order article via Infotrieve]
- Zirlinger M. Selection and validation of microarray candidate genes from subregions and subnuclei of the brain. Methods. 2003;31:290300.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Dent GW, ODell DM, Eberwine JH. Gene expression profiling in the amygdala: an approach to examine the molecular substrates of mammalian behavior. Physiol Behav. 2001;73:841847.[CrossRef][Medline][Order article via Infotrieve]
- Barrett T, Xie T, Piao Y, et al. A murine dopamine neuron-specific cDNA library and microarray: increased COX1 expression during methamphetamine neurotoxicity. Neurobiol Dis. 2001;8:822833.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Goldberg JL, Espinosa JS, Xu Y, Davidson N, Kovacs GT, Barres BA. Retinal ganglion cells do not extend axons by default: promotion by neurotrophic signaling and electrical activity. Neuron. 2002;33:689702.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Meyer-Franke A, Kaplan MR, Pfrieger FW, Barres BA. Characterization of the signaling interactions that promote the survival and growth of developing retinal ganglion cells in culture. Neuron. 1995;15:805819.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Barres BA, Silverstein BE, Corey DP, Chun LL. Immunological, morphological, and electrophysiological variation among retinal ganglion cells purified by panning. Neuron. 1988;1:791803.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Goldberg JL, Klassen MP, Hua Y, Barres BA. Amacrine-signaled loss of intrinsic axon growth ability by retinal ganglion cells. Science. 2002;296:18601864.[Abstract/Free Full Text]
- Huang X, Madan A. CAP3: A DNA sequence assembly program. Genome Res. 1999;9:868877.[Abstract/Free Full Text]
- Pruitt KD, Katz KS, Sicotte H, Maglott DR. Introducing RefSeq and LocusLink: curated human genome resources at the NCBI. Trends Genet. 2000;16:4447.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Pruitt KD, Maglott DR. RefSeq and LocusLink: NCBI gene-centered resources. Nucleic Acids Res. 2001;29:137140.[Abstract/Free Full Text]
- Altschul SF, Madden TL, Schaffer AA, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:33893402.[Abstract/Free Full Text]
- Schuler GD, Boguski MS, Stewart EA, et al. A gene map of the human genome. Science. 1996;274:540546.[Abstract/Free Full Text]
- Boguski MS, Lowe TM, Tolstoshev CM. dbEST: database for "expressed sequence tags.". Nat Genet. 1993;4:332333.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Wheelan SJ, Church DM, Ostell JM. Spidey: a tool for mRNA-to-genomic alignments. Genome Res. 2001;11:19521957.[Abstract/Free Full Text]
- Xiao HS, Huang QH, Zhang FX, et al. Identification of gene expression profile of dorsal root ganglion in the rat peripheral axotomy model of neuropathic pain. Proc Natl Acad Sci USA. 2002;99:83608365.[Abstract/Free Full Text]
- Schiller MR, Mains RE, Eipper BA. A neuroendocrine-specific protein localized to the endoplasmic reticulum by distal degradation. J Biol Chem. 1995;270:2612926138.[Abstract/Free Full Text]
- Bogumil R, Knipp M, Fundel SM, Vasak M. Characterization of dimethylargininase from bovine brain: evidence for a zinc binding site. Biochemistry. 1998;37:47914798.[CrossRef][Medline][Order article via Infotrieve]
- Knipp M, Charnock JM, Garner CD, Vasak M. Structural and functional characterization of the Zn(II) site in dimethylargininase-1 (DDAH-1) from bovine brain: Zn(II) release activates DDAH-1. J Biol Chem. 2001;276:4044940456.[Abstract/Free Full Text]
- Kimoto M, Tsuji H, Ogawa T, Sasaoka K. Detection of NG,NG-dimethylarginine dimethylaminohydrolase in the nitric oxide-generating systems of rats using monoclonal antibody. Arch Biochem Biophys. 1993;300:657662.[CrossRef][ISI][Medline][Order article via Infotrieve]
- MacAllister RJ, Parry H, Kimoto M, et al. Regulation of nitric oxide synthesis by dimethylarginine dimethylaminohydrolase. Br J Pharmacol. 1996;119:15331540.[ISI][Medline][Order article via Infotrieve]
- Niinobe M, Koyama K, Yoshikawa K. Cellular and subcellular localization of necdin in fetal and adult mouse brain. Dev Neurosci. 2000;22:310319.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Aizawa T, Maruyama K, Kondo H, Yoshikawa K. Expression of necdin, an embryonal carcinoma-derived nuclear protein, in developing mouse brain. Brain Res Dev Brain Res. 1992;68:265274.[Medline][Order article via Infotrieve]
- Stoilova D, Child A, Trifan OC, Crick RP, Coakes RL, Sarfarazi M. Localization of a locus (GLC1B) for adult-onset primary open angle glaucoma to the 2cen-q13 region. Genomics. 1996;36:142150.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Wirtz MK, Samples JR, Kramer PL, et al. Mapping a gene for adult-onset primary open-angle glaucoma to chromosome 3q. Am J Hum Genet. 1997;60:296304.[ISI][Medline][Order article via Infotrieve]
- Trifan OC, Traboulsi EI, Stoilova D, et al. A third locus (GLC1D) for adult-onset primary open-angle glaucoma maps to the 8q23 region. Am J Ophthalmol. 1998;126:1728.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Wirtz MK, Samples JR, Rust K, et al. GLC1F, a new primary open-angle glaucoma locus, maps to 7q35q36. Arch Ophthalmol. 1999;117:237241.[Abstract/Free Full Text]
- Assink JJ, Tijmes NT, ten Brink JB, et al. A gene for X-linked optic atrophy is closely linked to the Xp11.4-Xp11.2 region of the X chromosome. Am J Hum Genet. 1997;61:934939.[ISI][Medline][Order article via Infotrieve]
- Kerrison JB, Arnould VJ, Ferraz Sallum JM, et al. Genetic heterogeneity of dominant optic atrophy, Kjer type: Identification of a second locus on chromosome 18q12.2-12.3. Arch Ophthalmol. 1999;117:805810.[Abstract/Free Full Text]
- El-Shanti H, Lidral AC, Jarrah N, Druhan L, Ajlouni K. Homozygosity mapping identifies an additional locus for Wolfram syndrome on chromosome 4q. Am J Hum Genet. 2000;66:12291236.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Bomont P, Watanabe M, Gershoni-Barush R, et al. Homozygosity mapping of spinocerebellar ataxia with cerebellar atrophy and peripheral neuropathy to 9q33-34, and with hearing impairment and optic atrophy to 6p21-23. Eur J Hum Genet. 2000;8:986990.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Kremer H, Hamel BC, van den Helm B, et al. Localization of the gene (or genes) for a syndrome with X-linked mental retardation, ataxia, weakness, hearing impairment, loss of vision and a fatal course in early childhood. Hum Genet. 1996;98:513517.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Bishop JO, Smith GP. The determination of RNA homogeneity by molecular hybridization. Cell. 1974;3:341346.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Getz MJ, Birnie GD, Young BD, MacPhail E, Paul J. A kinetic estimation of base sequence complexity of nuclear poly(A)-containing RNA in mouse Friend cells. Cell. 1975;4:121129.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Bloomquist BT, Darlington DN, Mains RE, Eipper BA. RESP18, a novel endocrine secretory protein transcript, and four other transcripts are regulated in parallel with pro-opiomelanocortin in melanotropes. J Biol Chem. 1994;269:91139122.[Abstract/Free Full Text]
- Darlington DN, Schiller MR, Mains RE, Eipper BA. Expression of RESP18 in peptidergic and catecholaminergic neurons. J Histochem Cytochem. 1997;45:12651277.[Abstract/Free Full Text]
- Schiller MR, Mains RE, Eipper BA. A novel neuroendocrine intracellular signaling pathway. Mol Endocrinol. 1997;11:18461857.[Abstract/Free Full Text]
- Franco-Bourland RE, Guizar-Sahagun G, Garcia GA, et al. Retinal vulnerability to glutamate excitotoxicity in canine glaucoma: induction of neuronal nitric oxide synthase in retinal ganglion cells. Proc West Pharmacol Soc. 1998;41:201204.[Medline][Order article via Infotrieve]
- Morgan J, Caprioli J, Koseki Y. Nitric oxide mediates excitotoxic and anoxic damage in rat retinal ganglion cells cocultured with astroglia. Arch Ophthalmol. 1999;117:15241529.[Abstract/Free Full Text]
- Neufeld AH, Hernandez MR, Gonzalez M. Nitric oxide synthase in the human glaucomatous optic nerve head. Arch Ophthalmol. 1997;115:497503.[Abstract]
- Neufeld AH, Sawada A, Becker B. Inhibition of nitric-oxide synthase 2 by aminoguanidine provides neuroprotection of retinal ganglion cells in a rat model of chronic glaucoma. Proc Natl Acad Sci USA. 1999;96:99449948.[Abstract/Free Full Text]
- Smith MA, Vasak M, Knipp M, Castellani RJ, Perry G. Dimethylargininase, a nitric oxide regulatory protein, in Alzheimer disease. Free Radic Biol Med. 1998;25:898902.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Nakagomi S, Kiryu-Seo S, Kimoto M, Emson PC, Kiyama H. Dimethylarginine dimethylaminohydrolase (DDAH) as a nerve-injury-associated molecule: mRNA localization in the rat brain and its coincident up-regulation with neuronal NO synthase (nNOS) in axotomized motoneurons. Eur J Neurosci. 1999;11:21602166.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Surgucheva I, McMahan B, Ahmed F, Tomarev S, Wax MB, Surguchov A. Synucleins in glaucoma: implication of gamma-synuclein in glaucomatous alterations in the optic nerve. J Neurosci Res. 2002;68:97106.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Wang AG, Chen CH, Yang CW, et al. Change of gene expression profiles in the retina following optic nerve injury. Brain Res Mol Brain Res. 2002;101:8292.[Medline][Order article via Infotrieve]
- Szabo A, Dalmau J, Manley G, et al. HuD, a paraneoplastic encephalomyelitis antigen, contains RNA-binding domains and is homologous to Elav and Sex-lethal. Cell. 1991;67:325333.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Lee CH, Della NG, Chew CE, Zack DJ. Rin, a neuron-specific and calmodulin-binding small G-protein, and Rit define a novel subfamily of ras proteins. J Neurosci. 1996;16:67846794.[Abstract/Free Full Text]
- Wes PD, Yu M, Montell C. RIC, a calmodulin-binding Ras-like GTPase. EMBO J. 1996;15:58395848.[ISI][Medline][Order article via Infotrieve]
- Sano A, Mikami M, Nakamura M, Ueno S, Tanabe H, Kaneko S. Positional candidate approach for the gene responsible for benign adult familial myoclonic epilepsy. Epilepsia. 2002;43:2631.[CrossRef]
- Xiang M, Zhou L, Macke JP, et al. The Brn-3 family of POU-domain factors: primary structure, binding specificity, and expression in subsets of retinal ganglion cells and somatosensory neurons. J Neurosci. 1995;15:47624785.[Abstract]
- Gooley JJ, Lu J, Chou TC, Scammell TE, Saper CB. Melanopsin in cells of origin of the retinohypothalamic tract. Nat Neurosci. 2001;4:1165.[CrossRef][ISI][Medline][Order article via Infotrieve]
- Hattar S, Liao HW, Takao M, Berson DM, Yau KW. Melanopsin-containing retinal ganglion cells: architecture, projections, and intrinsic photosensitivity. Science. 2002;295:10651070.[Abstract/Free Full Text]
- Provencio I, Rodriguez IR, Jiang G, Hayes WP, Moreira EF, Rollag MD. A novel human opsin in the inner retina. J Neurosci. 2000;20:600605.[Abstract/Free Full Text]
- Blackshaw S, Kuo WP, Park PJ, et al. MicroSAGE is highly representative and reproducible but reveals major differences in gene expression among samples obtained from similar tissues. Genome Biol. 2003;4:R17.[CrossRef][Medline][Order article via Infotrieve]
This article has been cited by other articles:

|
 |

|
 |
 
I. Surgucheva, V. I. Shestopalov, and A. Surguchov
Effect of {gamma}-Synuclein Silencing on Apoptotic Pathways in Retinal Ganglion Cells
J. Biol. Chem.,
December 26, 2008;
283(52):
36377 - 36385.
[Abstract]
[Full Text]
[PDF]
|
 |
|