|
|
||||||||
1From the Institute of Ophthalmology, Molecular Genetics, University College, London, United Kingdom; the 2Department of Visual Neuroscience, Imperial College, London, United Kingdom; and the 3Macular Research Centre, University of British Columbia, Vancouver, British Columbia, Canada.
| Abstract |
|---|
|
|
|---|
METHODS. Two-millimeter foveomacular and four-millimeter macular punches from human donor eyes, in addition to sections of midperipheral retina, were used to study differential gene expression. Multiple microarray experiments were combined with quantitative PCR and bioinformatic analyses. In the present study, the expression of both known and previously unidentified retinal genes was determined.
RESULTS. Several macula enriched transcripts were revealed. Nuclear pore complex interacting protein (NPIP) and eukaryotic translation initiation factor 2
kinase (GCN2) were expressed at levels approaching that of red/green cone opsin in the macula. The protein products of several genes highlighted using these expression analyses were also localized in the retina. Both NPIP and histone deacetylase 9 (HDAC9) proteins were detected in cone photoreceptor outer segments.
CONCLUSIONS. Characterizing macula enriched transcripts is an important stepping-stone in understanding the molecular basis for visual acuity in the retina. The approach also provides excellent candidates for diseases that affect the macula and fovea such as age-related macular degeneration (AMD). Indeed, several of these transcripts, such as NPIP and GCN2, have genomic loci that are consistent with being candidate genes for AMD.
Although many expressed sequence tags (ESTs) have been identified from central retinal cDNA libraries (www.ncbi.nlm.nih.gov/dbEST/ National Center for Biotechnology Information [NCBI], Bethesda, MD), most are unlikely to vary significantly in their expression throughout the retina. Despite this, Northern blot analysis has shown that there are significant differences in gene expression between the central and peripheral retina.2 More recently, serial analysis of gene expression (SAGE) comparing the perimacular region (6 mm) with the peripheral retina has revealed several genes that are overrepresented in the central retina (www.ncbi.nlm.nih.gov/geo/ NCBI), HMAC2_SAGE).3 However, a greater number of gene tags have been found to be differentially expressed between two separate peripheral retinas (HPR1/2_SAGE).
Age-related macular degeneration (AMD) is a multifactorial disease that causes loss of central vision. It is the most common cause of blindness in the developed world and is responsible for more than 50% of blind or partial sight registrations in the UK.4 Approximately 10% of people older than 65 and 30% older than 75 show retinal changes consistent with AMD. Moreover, the prevalence of AMD is increasing.5 Twin and sibling studies suggest that the genetic component of AMD is at least one quarter,6 7 8 9 with smoking being the only well established environmental risk factor. Several whole genome scans for linkage in families with AMD have defined four main chromosomal regions of susceptibility to AMD: 1q13, 6p, 10q26, and 16p12. Recently, complement factor H (CFH/HF1) has been implicated as the susceptibility gene at 1q13, and some genetic variants are likely to double the risk of having AMD.10 11 12 13 Variants in two other regulatory components of the alternative complement pathway BF and C2 have also been associated with increased or decreased risk of AMD.14 At the 10q26 locus HTRA1, a serine protease activated by cellular stress has been associated with major susceptibility to AMD independent of CFH.15 16
As AMD susceptibility is due to both environmental and hereditary factors, gene expression analyses may help to understand the molecular context of the disease. Microarray-initiated experiments have been used in the retina to identify both disease-causing (autosomal dominant retinitis pigmentosa17 ) and major regulatory (homeobox CRX18 19 ) genes.
In the present study, we identified retinal genes that are consistently overexpressed in the macula, some of which are also abundant transcripts. An important observation was that two of the protein products of these transcripts are enriched in cone photoreceptor outer segments. One of these genes, nuclear pore complex interacting protein (NPIP), lies at an as yet undefined susceptibility locus for AMD (16p12).
| Methods |
|---|
|
|
|---|
Northern Blot Analysis
Probes were designed from human rod and long/medium-wave cone opsin of lengths 345 and 424 bp respectively (primers available on request). These probes were amplified individually by RT-PCR from retinal RNA (SuperScript; Invitrogen). The probes were labeled with
32-P-dCTP using a labeling kit (Prime-a-Gene; Promega, Madison, WI) and cleaned up with a nucleotide removal kit (QIAquick; Qiagen, Crawley, CA). For each probe 5 µg of total RNA, from peripheral retina and 4 mm macula, were electrophoresed in adjacent lanes of a formaldehyde agarose gel and transferred to a nylon membrane. The probes were hybridized (ULTRAhyb; Ambion) according to the manufacturers protocol.
Microarrays
Microarray slides were designed at the Gene Array Facility, University of British Columbia (Vancouver, BC, Canada). The set consisted of 13,899 70-mers from a probe set (Operon ver. 1.1; Qiagen), and these were spotted robotically in duplicate onto glass slides, together with control probes contained in the set. Each slide was subject to strict quality control. Targets were prepared with 7.5 µg total RNA, with either Cy3 or Cy5 dUTP and reverse transcriptase (Superscript II; Invitrogen) and were cleaned up with a PCR purification kit (QiaQuick; Qiagen) in light-shielded conditions and precipitated overnight with ethanol and sodium acetate. Each target was resuspended in 50 µL hybridization solution containing: 25 µL formamide, 12.5 µL 20x SSC, 0.5 µL 10% SDS, 5 µL of 2 g/L BSA, 5 µL of 5 g/L yeast tRNA and 2 µL of 10 g/L salmon testes DNA. The slides were denatured at 95°C just before use. The Cy3 and Cy5 target were mixed together and hybridized to the slide under a coverslip in a moisture-rich environment overnight at 42°C. The coverslip was removed by floating the slide in 0.2x SSC. The slide was then washed three times for 5 minutes with 0.1x SSC and 0.1% SDS and three times for 5 minutes in 0.1x SSC. The slide was centrifuged in a tube at 2000 rpm for 5 minutes to dry before imaging (ChipReader; Virtek Vision Corp, Waterloo, ON, Canada).
Grids were overlaid onto images of each slide by computer (ImageGene software; BioDiscovery, Los Angeles, CA), which was also used to quantify each spot in the two channels corresponding to the Cy3 and Cy5 fluorescence images. The individual spot intensities were subjected to an in-house background correction algorithm that uses the lowest 10% of intensity values for each subgrid as background. In microarray analysis, typically less than half of the probes light up with a particular target. The data were parsed by using custom-written algorithms to associate each spot with its corresponding probe identity. Data were then calibrated by using several different statistical methods, including log-ratios and variance stabilization (www.dkfz-heidelberg.de/mga/whuber/ German Cancer Research Center, Heidelberg, Germany).20 Data were analyzed by writing scripts in the R statistical computing language (www.r-project.org). The processed data from each slide were inserted into an SQL database to facilitate data mining.
Quantitative PCR
First, intron-spanning primers (Supplementary Data, online at http://www.iovs.org/cgi/content/full/48/12/5388/DC1) were designed within two adjacent exons of rhodopsin. Two micrograms of pooled peripheral retinal and pooled foveomacular RNA (as measured by spectrophotometry) was reverse transcribed with random decamers (RetroScript; Ambion). Peripheral retinal cDNA was amplified with the intron-spanning rhodopsin primer pair to check for genomic DNA contamination.
Primer pairs (Supplementary Data) were generated to regions of the cDNAs chosen for qPCR analysis, including four internal controls (18S rRNA, acidic ribosomal phosphoprotein, β-actin, and ubiquitin C), by using the GCG-Wisconsin package (Accelrys, Campbell, CA). Peripheral retinal cDNA was used to find optimal primer pairs to amplify each chosen cDNA specifically, with SYBR-green PCR master mix (Applied Biosystems, Inc.; Foster City, CA). Products were checked on agarose. Optimized reactions were melted (model 7700; ABI) to check for a single product and to discern its melting temperature. To double check the specificity of the primers, four PCR products (RHO, OPN1LW, 6q, and BHM2) were subjected to direct sequencing.
Quantitative PCR reactions were performed in triplicate for each gene with peripheral, foveomacular, and no-template control samples (model 7700; ABI). An additional PCR step was inserted at the appropriate melting temperature at which to measure the fluorescence of the product. Cyclical fluorescence values were recorded and preprocessed with the system software (ABI). These data were then analyzed on a custom-written spreadsheet (Excel; Microsoft, Redmond, WA; DART-PCR21 ) to extrapolate the initial fluorescence values and normalize them to the internal controls.
Immunocytochemistry
A rabbit polyclonal antibody to the C terminus of HDAC9 (NP_848512, 1069 amino acids) was obtained from the laboratory of Arthur Zelent (Institute of Cancer Research). It was raised to the peptide DVEQPFAQEDSRTAG, from amino acid 1046-1060 that detects the full-length protein.22 For NPIP, a 15-amino-acid synthetic peptide, CSLPFHPQRMIISRN, corresponding to the C terminus was used for production of rabbit polyclonal antisera (Sigma-Aldrich, St. Louis, MO).
For immunocytochemistry, tissue culture cells were washed, fixed, and permeabilized using 3.7% paraformaldehyde and 0.1% Triton X-100 in eight-well chamber glass slides (Nunc, Naperville, IL) before blocking with normal donkey serum (NDS 1:10) in BSA 3%, to counter nonspecific binding of the primary antibody. All incubation steps were performed at room temperature, and appropriate washing steps were included. 4',6-Diamidino-2-phenylindole (DAPI) was used to counterstain nuclei. The retinal cryosections were processed under similar conditions except that Triton was not used. Control experiments were performed in parallel with omission of the primary antibody only. A confocal laser scanning microscope (LSM510; Carl Zeiss Meditec) was used to image the slides.
Cloning
Primers were designed (Supplementary Data) to add restriction sites to each end of NPIP to facilitate subcloning into the GFP fusion expression vector (pEGFP-C2; BD-Clontech, Palo Alto, CA) that fuses the insert to the C terminus of GFP. A 5' EcoRI site and 3' SalI site were chosen, as they are compatible with in-frame subcloning. NPIP was amplified from brain cDNA, and the products were gel purified and cloned into a vector (pGEM-T Easy; Promega). DNA extracted from the amplified vector was sequenced bidirectionally before proceeding to the subcloning step.
| Results |
|---|
|
|
|---|
Genes Enriched in Macula and Fovea
Six duplicate microarray experiments were performed to find differentially expressed genes in the human macula (Table 1) . Expression data (Fig. 1) were plotted as the difference in expression between channels: M = h(Cy5) – h(Cy3) against the rank of the mean expression
![]() |
|
|
Data were mined by querying an SQL database for probes that showed consistent differential expression (Table 2) , according to strict selection criteria based on diagnostic plots and analyses. First, genes that did not exhibit an absolute fluorescence of at least 200 above background in at least two channels were discarded, because this was less than twice the highest background value. Second, genes that showed variability of one expression unit (twofold) or greater were excluded. All central retina enriched genes included had expression ratios of >3 (eightfold increase) and also had higher expression ratios in the foveomacula-versus-periphery compared with the macula-versus-periphery experiments. Peripheral enriched genes were included, for comparison, with an expression of –2 (fourfold decrease) or less. These strict selection criteria yielded several candidate differentially expressed genes for further analysis. Many more genes were found to be differentially expressed in macula compared with peripheral retina across the whole data set (Supplementary Data).
|
The data were analyzed using the DART-PCR system.21 DART-PCR uses an amplification plot method that is more conservative in calculating the expression change (in x-fold) than the common 2–
Ct estimate of reaction efficiency. The expression level of each gene and the ratio in foveomacula:peripheral retina was calculated (Table 3) and plotted on a log base 2 scale (Fig. 2) , where +1 is a twofold increase and –1 is a twofold decrease in expression. The fluorescence due to SYBR green is also proportional to the length of DNA fragment. By correcting for amplicon length, the percentage expression of selected genes relative to either rod or cone opsin were also calculated (Table 3) . Reaction efficiency, as measured by the amplification plot method, typically varied by up to 10% between any two given probes, so that the values of relative expression are estimates.
|
|
kinase (GCN2) showed high expression levels, which approached that of red/green cone opsin in the foveomacula. GCN2 phosphorylates the
subunit of eukaryotic translation initiation factor 2 and is thought to be involved in the cellular response to stress.24 HDAC9 was more than threefold enriched and SRGAP3 was approximately twofold enriched in foveomacula. SRGAP3 is thought to play an important role in neuronal migration and axonal connectivity.25 S-antigen (rod arrestin), the only peripheral retina enriched gene that fitted the strict microarray selection criteria, was over five times enhanced in the periphery according to the qPCR data, with expression approaching that of rhodopsin.
HDAC9 in Macula and Photoreceptors
HDACs are chromatin-remodeling factors that act as transcriptional repressors in response to signals such as cellular stress. HDAC9 was originally identified while investigating cDNAs that encode large proteins expressed in the human brain.26 HDAC9s huge genomic size (7p21.1, 500 kb) and the degree to which it is regulated are unprecedented. This finding may indicate a wider role for HDAC9 than just histone modification.22
HDAC9 mRNA was eightfold enriched in the foveomacula according to the microarray data and was the most highly enriched gene after cone opsin according to the qPCR data. Full-length HDAC9 was localized to the nucleus of human neuroblastoma (SK-N-SH) cells (Fig. 3) using immunofluorescence with an antibody to the C terminus of HDAC9. Full-length HDAC9 has been localized to the nucleus by transiently transfecting simian fibroblast (COS-7) cells with an HDAC9-FLAG construct, followed by immunofluorescence with an anti-FLAG antibody.22 Of interest, HDAC9 was also detected in association with filamentous structures in the SK-N-SH cells (Fig. 3) .
|
NPIP in Cone Outer Segments
NPIP is a member of a primate-specific gene family that has undergone recent adaptive evolution.27 The protein contains one predicted transmembrane domain and is located on the nuclear envelope in cultured cells.
NPIP was eightfold overexpressed in the foveomacula compared with the peripheral retina according to the microarray data. It was the second most highly enriched gene after cone opsin according to the qPCR data. NPIP was very highly expressed in foveomacula at approximately 76% of the level of cone opsin. Endogenous NPIP was localized to the nuclear envelope region of adherent SK-N-SH cells (Supplementary Data), by using immunofluorescence with an antibody raised to the C terminus, as determined by nuclear counterstain. These data for the endogenous protein were in accordance with a previous report that localized NPIP-GFP to this subcellular location by transiently transfecting COS-7 cells.27
We localized NPIP in the retina with the anti-NPIP antibody on cryosections of human fovea (Fig. 4) . This revealed specific staining on the sheaths of the cone outer segments. There was apparent colocalization with peanut agglutinin (PNA), which stains cone photoreceptor sheaths.
|
| Discussion |
|---|
|
|
|---|
From diagnostic analysis of the microarray data, variance stabilization20 was chosen as the most appropriate method with which to perform calibration. Variance stabilization is a recent normalization algorithm that is designed to produce a constant signal-to-noise ratio over the entire dataset. This method allowed us to distinguish meaningful data from noise by using human macular mRNA. Furthermore, although many genes were found to be enriched in the macula according to the microarrays, only the few that fit strict data-mining criteria were studied further. Similarly, using a sample pooling approach together with an amplification plot method of efficiency calculation,21 we were able to obtain reasonable estimates of macular gene expression relative to the opsins.
The gene products of NPIP, GCN2, SRGAP3, and HDAC9 were previously uncharacterized in the retina. Our data provide evidence that these transcripts are macula-enriched and that NPIP and GCN2 represent a significant proportion of the macular transcriptome. Only SRGAP3 is as yet associated with a monogenic disease phenotype, in which abnormal development of neuronal structures leads to mental retardation.25 However, NPIP (16p13-11), GCN2 (15q15) and SRGAP3 (3p25) are all located within possible regions of linkage suggested by whole-genome screens of age-related maculopathy.28 29 30
We localized full-length HDAC9 to the nucleus of cultured cells, as would be expected by its HDAC activity and as previously reported.22 However, it was also localized to filamentous structures. Of note, HDAC6 has been found to act as a molecular motor in the cytoskeleton.31 In the human fovea, we found that HDAC9 localized to photoreceptors, especially the outer segments. Transcriptional repression (by histone deacetylation) may have a protective role in neuronal ageing and degeneration by preventing a stress response that leads to apoptosis of neuronal cells.32 It is also possible that HDAC9 is involved in maintaining structural integrity in the photoreceptor via cytoskeletal elements and/or protein trafficking, perhaps via the ubiquitin-proteasome degradation system.33 This possibility may be especially important in the extremely cell dense macula.
NPIP is a member of a gene family originating from human chromosome 16 (LCR16a) that has been duplicated very recently and evolved rapidly.27 It is an extreme example of positive selection, with up to 43% divergence of amino acids between primates, a hallmark of adaptive evolution. BLAST analysis of the reference sequence against dbEST or searching for motifs with PROSITE confirmed that NPIP is primate specific. Strikingly, among mammals it is only primate retinas that possess a macula. (Birds, such as hawks, have a very different retinal structure and can be bifoveate.) Apart from approximately 25% similarity with chicken caldesmon (NP_989489), and a portion of an archaeal coiled coil domain (NP_613407), there are no significant similarities of NPIP with any other known protein. Johnson et al.27 described eight exons for NPIP, determined by genomic sequence analysis and reported a single putative transmembrane domain (from amino acid 54-72) surrounded by
-coil.
This is the first report to our knowledge of tissue localization of NPIP. It is of particular interest that NPIP is localized to cone photoreceptor outer segments in humans and that NPIP appears to be a significant part of the macular transcriptome. The longer novel isoform NPIP1 and the reference isoform AAD3439427 localize to the nuclear envelope, whereas the shorter novel isoform NPIP2 localizes within the nucleus of cultured cells. This suggests that there is a membrane localization signal at the N terminus of the longer isoforms. It may be this same localization signal that directs NPIP to the outer segments of human cone photoreceptors although this, along with ultrastructural localization, remains to be seen.
Although the role of NPIP in cone photoreceptors remains unclear, it is of interest that other proteins thought to be associated with the nuclear pore also have a function in photoreceptors. For example, the Drosophila eyes closed gene (Eyc, an analogue of p47) is involved in both nuclear envelope and photoreceptor membrane reorganization.34 There is also some evidence that RanBP2, a large protein that binds the nuclear pore-associated GTPase Ran, is highly expressed in cone photoreceptors of the vertebrate retina.35 Indeed, RanBP2 was thought to be a chaperone for red/green cone opsin.
Autosomal recessive retinitis pigmentosa (arRP) has been linked to a chromosomal region that overlaps with NPIP (16p12.3-p12.1).36 However, one of the two families showed physical signs frequently seen in patients with BBS (Bardet-Biedl syndrome). The other family reported night blindness from the second decade, which is characteristic of rod-cone dystrophy, although no retinal electrophysiology was published. NPIP is located within a region of suggested linkage (16p12) obtained from three separate genome-wide scans of age-related maculopathy28 29 30 and an additional genome-wide scan for predisposition to Drusen formation.37
Characterizing macula enriched transcripts is an important stepping-stone in understanding the molecular basis for visual acuity in the retina. This approach also provides excellent candidates for diseases that affect the macula and fovea such as AMD. Work to define roles for NPIP and the other genes highlighted in this study is ongoing.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication March 24, 2007; revised July 5, 2007; accepted October 8, 2007.
Disclosure: D.M. Hornan, None; S.N. Peirson , None; A.J. Hardcastle, None; R.S. Molday, None; M.E. Cheetham, None; A.R. Webster, 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: Daniel M. Hornan, Molecular Genetics, Institute of Ophthalmology, 11-43 Bath St., London EC1V 9EL, UK; d.hornan{at}doctors.org.uk.
| References |
|---|
|
|
|---|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |