|Molecular Vision 2007;
Received 12 January 2007 | Accepted 30 May 2007 | Published 18 June 2007
A microarray analysis of retinal transcripts that are controlled by image contrast in mice
Section for Neurobiology of the Eye, University Eye Hospital Tuebingen, Tuebingen, Germany
Correspondence to: Marita Pauline Feldkaemper, Section for Neurobiology of the Eye, University Eye Hospital Tuebingen, Calwerstr. 7/1, 72076 Tuebingen, Germany; Phone: (+49) 7071/29-87424; FAX: (+49) 7071/29-5196; email: email@example.com
Purpose: The development of myopia is controlled by still largely unknown retinal signals. The aim of this study was to investigate the changes in retinal mRNA expression after different periods of visual deprivation in mice, while controlling for retinal illuminance.
Methods: Each group consisted of three male C57BL/6 mice. Treatment periods were 30 min, 4 h, and 6+6 h. High spatial frequencies were filtered from the retinal image by frosted diffusers over one eye while the fellow eyes were covered by clear neutral density (ND) filters that exhibited similar light attenuating properties (0.1 log units) as the diffusers. For the final 30 min of the respective treatment period mice were individually placed in a clear Perspex cylinder that was positioned in the center of a rotating (60 degrees) large drum. The inside of the drum was covered with a 0.1 cyc/degree vertical square wave grating. This visual environment was chosen to standardize illuminances and contrasts seen by the mice. Labeled cRNA was prepared and hybridized to Affymetrix GeneChip® Mouse Genome 430 2.0 arrays. Alterations in mRNA expression levels of candidate genes with potential biological relevance were confirmed by semi-quantitative real-time reverse transcription polymerase chain reaction (RT-PCR).
Results: In all groups, Egr-1 mRNA expression was reduced in diffuser-treated eyes. Furthermore, the degradation of the spatial frequency spectrum also changed the cFos mRNA level, with reduced expression after 4 h of diffuser treatment. Other interesting candidates were Akt2, which was up-regulated after 30 min of deprivation and Mapk8ip3, a neuron specific JNK binding and scaffolding protein that was temporally regulated in the diffuser-treated eyes only.
Conclusions: The microarray analysis demonstrated a pattern of differential transcriptional changes, even though differences in the retinal images were restricted to spatial features. The candidate genes may provide further insight into the biochemical short-term changes following retinal image degradation in mice. Because deprivation of spatial vision leads to increased eye growth and myopia in both animals and humans, it is believed some of the identified genes play a role in myopia development.
Myopia is increasing in prevalence world-wide , but information remains fragmented about the regulation of ocular growth and the development of refractive errors. In animal models, it has been shown that the alteration of retinal images leads to altered gene expression and a change in eye growth patterns. Myopia can be artificially induced by placing negative lenses in front of the eye or by reducing the retinal image quality by a diffuser [2,3], whereas hyperopia can be generated by positive lens wear . These eye growth responses were shown in several animal models, including chicks , tree shrews , marmosets  and rhesus monkeys . In mice, a shift toward myopia can be induced by form deprivation [8-10]. During negative spectacle lens wear, axial eye growth rates are accelerated until a sharp retinal image is achieved. Depriving the retina of high contrast and high spatial frequencies by diffusers also induces axial eye growth, but a sharp retinal image cannot be restored in this case. A local mechanism within the retina is involved in eye growth regulation since neither accommodation nor contributions from the brain are necessary .
Little is known about genes and proteins whose expression is susceptible to altered visual stimulation. Among the known genes and substances are glucagon [12,13], retinoic acid , vasointestinal peptide (Vip) [15-17], sonic hedgehog (Shh) [18,19], and the transcription factor Egr-1 [20,21]. The identification of novel retinal genes that are influenced by altered visual conditions is important as they may lead to new targets for a pharmacological therapy of myopia. In the present study, DNA-microarrays, followed by semi-quantitative real-time RT-PCR, were used to screen for differentially regulated genes and subsequent validation.
Black wildtype C57BL/6 male mice were raised on a 12 h: 12 h light-dark cycle with light onset at 8:00 a.m. The mice were reared in the local animal facility and given free access to water and food. For the GeneChip® experiments the mice (3 per group) were studied at postnatal (P) days 30 and 32, while the mice used for the validation of the results by semi-quantitative real-time RT-PCR (6 animals per group) were studied at P29 to P36. The experimental treatment was in accordance with the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research and was approved by the university commission for animal welfare (reference AK 3/05).
Diffuser and neutral density filter design
Diffusers were made from transparencies that were sanded down with emery paper to give a frosted appearance. Neutral density (ND) filters of 0.1 log unit light attenuation (Kodak Wratten Gelatin Filter No. 96; Rochester, NY) were used to match image illumination [21,22]. Velcro rings with an inner diameter of 6 mm and a total diameter of 10 mm were fitted with the frosted transparency or ND filter. The contralateral eyes served as controls to minimize effects due to differences between individual animals in factors such as hormone level and immunological state.
Velcro rings were attached to the periorbital fur the day before the experiments while animals were under diethyl ether anesthesia. For anesthesia, animals were put in a diethyl ether containing glass jar supplied with a grid at the bottom. The glass jar was covered by a lid until the animals were unconscious (which takes about 20-30 s). The complementary Velcro rings carrying the diffusers and ND filters were attached on the experimental day. Mice were kept under cool white light of approximately 120 lux (Lumilux 30W/840; Osram, Munich, Germany) during the treatment period and were sacrificed after a total of six hours in light (Figure 1). At the end of the respective treatment period, each mouse was separately placed on a stationary platform in the center of the rotating drum with a diameter of 60 cm for 30 min. Illumination in the drum was about 400 lux. The stripe pattern had a spatial frequency of 0.1 cycles per degree. The angular velocity was about 60 degrees per second, which was much higher than that used by Prusky et al. . Previous studies in our lab, using the drum, have shown that an angular velocity of 50 to 60 deg/s was optimal for studying spatial vision in C57BL/6 mice in our optomotor set-up .
The whole body optomotor grating acuity of C57BL/6 mice was recently measured in the same drum and was found to be limited to approximately 0.3-0.4 cycles/degree. The stripe contrast measured at 400 lux was approximately 90% . Previous experiments have shown that the image contrast with the diffusers used in the current study is reduced to about 60-70% of the initial value . The illumination before and during the experiment was far below the illumination that was previously used to induce retinal degeneration in mice . To insure that gene expression changes were a result of the treatment and not the hormone status, age, light, or time of day, we only used male litter mates and sacrificed the mice after they had been exposed to 6 h of light. This was usually between 2 and 3 p.m. in the afternoon (Figure 1).
Tissue preparation and RNA isolation
All mice were sacrificed by an overdose of diethyl ether. Their eyes were enucleated and immediately transferred to a petri dish containing chilled Ringer's solution until preparation. The eyes were subsequently placed on a filter paper, perforated by a canula and opened with scissors, cutting around the iris. The lens was removed, and the retina was extracted, snap frozen in liquid nitrogen, and stored at -70 °C. The preparation was done using a stereo boom microscope (10X magnification) and RPE cells were carefully removed. The total length of time between death/enucleation and tissue freezing was about 5-10 min. Each retina was treated and analyzed as a separate sample. The retinas were homogenized for 1 min with speed increasing from 11,000 to 20,000 rpm (Diax 900 Homogenizer; Heidolph, Kelheim, Germany). Total RNA was isolated with a kit (RNeasy Mini Kit; Qiagen, Hilden, Germany) according to the manufacturer's instructions. RNA quality and quantity were determined by the Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA), which yielded RNA integrity numbers (RIN) from 8.4 to 9.3.
The GeneChip® mouse genome 430 2.0 Array from Affymetrix (Santa Clara, CA) allows the comprehensive analysis of genome-wide expression on a single array. We used 45,000 probe sets to analyze the expression level of over 39,000 transcripts and variants from over 34,000 well-characterized mouse genes.
Microarray analysis was performed by the Affymetrix Resource Facility at the University of Tuebingen as follows. Target labeling for the expression analysis was performed for 1 μg total RNA per sample according to the manufacturer's protocol (Affymetrix). The GeneChips® were automatically stained and washed in a fluidics station as recommended by the manufacturer. The scanning and the analysis were done using the Affymetrix Microarray Suite Software (v. 5.0).
Affymetrix data analysis and statistics
The relative abundance of individual genes is based on the signal intensities of the corresponding probe sets that were analyzed by ArrayAssist 4.0 (Stratagene, La Jolla, CA). The microarray data were not Benjamini/Hochberg corrected since no significant differences in gene expression were obtained when this procedure was performed. Normalization and confirmation is covered in the Discussion section. Data from the experiments involving 30 min, 4 h and 6+6 h treatments were individually normalized using robust multi-array analysis with correction for GC content (gcRMA). The strategy of this normalization method is the calculation of a background adjustment step that ignores the mismatch (MM) intensities but incorporates sequence information from the probes (GC content) . Afterwards the data were analyzed group-wise by paired t-tests to find genes that were differentially expressed in the diffuser- and filter-treated eyes of the same animal. The p-value was below 0.05 and a minimum fold-change (FC) was 1.5. mRNA level changes observed throughout the diffuser treatment were evaluated by normalizing together raw data from all three microarray experiments (30 min, 4 h, 6+6 h). One kind of background adjustment was calculated for all the 18 arrays. The quality of all microarray experiments was assessed by GAPDH 3'/5' and Actin 3'/5' ratios as well as background, scalar factors and present calls (Table 1) .
Network generation and pathway analysis
The networks and functional analyses were generated through the use of Ingenuity Pathways Analysis (Ingenuity® Systems).
A data set containing gene identifiers and corresponding fold change values was uploaded into the application. Each gene identifier was mapped to its corresponding gene object in the Ingenuity Pathways Knowledge Base. These genes, called focus genes, were overlaid onto a global molecular network developed from information contained in the Ingenuity Pathways Knowledge Base. Networks of these focus genes were then algorithmically generated based on their connectivity.
Semi-quantitative real-time polymerase chain reaction
All RNA samples used for semi-quantitative RT-PCR were treated with RNase-free DNase I (Roche, Mannheim, Germany). Gel electrophoresis was used to check quality of RNA. Concentration and purity were determined by spectrophotometry at 260 and 280 nm (average ratio 1.95±0.08). Of each sample, 1 μg total RNA was reverse transcribed (M-MLV reverse transcriptase; Promega, Madison, WI) using a combination of 50 ng random hexamers and 500 ng oligo (dT)15 primers in a total volume of 20 μl.
Gene sequences chosen for validation were obtained from the National Center for Biotechnology Information NCBI. Primers (Table 2) were preferentially designed to bind within the same coding region as the Affymetrix probes using a commercial program (Prime; GeniusNet Husar; KYE Systems, Heidelberg, Germany) and Primer Premier 5 (Premier Biosoft International; Palo Alto, CA). They were ordered from a commercial synthesis service (VBC Genomics, Vienna, Austria).
PCR was performed in a thermocycler (iCycler; Bio-Rad, Hercules, CA) using a commercial fluorescence detection kit (QuantiTect SYBR Green PCR kit; Qiagen). After an initial heat activation (15 min at 95 °C) 40 cycles of 15 s at 94 °C, 30 s at 59 °C and 45 s at 72 °C were run. Cyclic fluorescence measurements were taken at the end of the annealing phase. The volume of a single reaction added up to 15 μl containing 2 ng template and a final primer concentration of 0.6 μM each. All PCR products were verified by automated sequencing or by restriction enzyme digestion (Usp36). Both the acidic ribosomal protein (Arp) and β-actin were chosen as internal reference genes and used for subsequent analyses.
Semi-quantitative real-time reverse transcriptase polymerase chain reaction data analysis and statistics
Measurements were performed in triplicates and the CT (Threshold Cycle) means were used for further data analyses by following a procedure described in reference .
Briefly, the efficiencies (E) of the individual primer sets during PCR were determined by dilution series and following equation: E=10(-1/slope).
Statistical data analysis was performed after conversion of the raw data into normalized expression (NE). Therefore, the CT means were transformed into relative quantities (RQ), and corresponding normalization factors (NF) were determined.
For the RQs, all CT values obtained for a single gene were considered (e.g. control, 30 min, 4 h, and 6+6 h group). The sample showing the highest expression level (e.g. the lowest CT) was set to 1 (reference sample, RS). All other RQs were calculated by following formula:
Relative Quantity (Sample, S)=Efficiency(CT (RS)-CT (S)).
The normalization factor represents the geometric mean (root) of the reference gene RQs:
The NE represents the ratio of the relative quantity of a sample and its corresponding normalization factor: NE=RQ/NF
NE differences between diffuser-treated and ND filter-covered fellow eyes were calculated and statistically analyzed by un-paired Student's t-tests against null.
GeneChip® analysis was performed on three animals per treatment group. mRNA expression levels in the diffuser-covered eyes were calculated relative to those in the ND filter-treated fellow eyes. After gcRMA normalization, expression of 16 genes was found to be affected after 30 min, with 13 genes being up-regulated. The 4 h treatment resulted in 27 differentially expressed genes, 23 of which were down-regulated. After 6+6 h, the number of up- and down regulated genes was balanced, with 10 genes being up-regulated and 11 being down-regulated.
The genes were arranged in a list according to their major function as determined by Ingenuity Pathways Analysis (Table 3).
To identify diffuser specific genes that were temporally regulated, the raw data of all three microarray experiments were normalized together and expression changes were determined by pair-wise comparisons. The p-value was set to 0.01 and a minimum fold change of 2.0 was established. The expression of 249 genes was significantly different after 6+6 h versus 30 min of treatment, and that of 191 genes was significantly different after 4 h versus 30 min of treatment but the expression of only 49 genes was significantly different after 6+6 h versus 4 h of treatment. Six genes were chosen to confirm the temporal expression analysis: fibroblast growth factor 13 (Fgf13) mitogen-activated protein kinase 8 interacting protein 3 (Mapk8ip3), trafficking protein, kinesin binding 2 (Trak2), Jagged 1 (Jag1), serine/threonine kinase 3 (Stk3), and kruppel-like factor 9 (Klf9).
Semi-quantitative real-time reverse transcriptase polymerase chain reaction
mRNA expression levels were evaluated by real-time RT-PCR. Six animals were used in each group. Genes chosen for validation met certain criteria. First, the signal intensity in either group (diffuser or ND filter) exceeded 65, and second, the Affymetrix probe sets covered a coding region (as determined by the University of California Santa Cruz genome browser). The genes shown in Table 3 were identified after gcRMA normalization of the data set. A previous analysis of the 30 min experiment using MAS5 normalization resulted in a higher number of differentially regulated genes and the confirmation of three genes (Akt2, Krt2-6b, Myb) from that list. The issue of normalization is dealt with in the Discussion section. Accordingly, eight potentially interesting genes were selected for validation: early growth response 1 (Egr-1), protein kinase B (Akt2), keratin complex 2, basic, gene 6b (Krt2-6b), myeloblastosis oncogene (Myb), FBJ osteosarcoma oncogene (cFos), ubiquitin specific peptidase 36 (Usp36), zinc finger, MYM-type 4 (Zmym4), and dedicator of cytokinesis 4 (Dock4). The mRNA expression levels of these genes were measured in diffuser- and ND filter-treated eyes by RT-PCR (Figure 2).
The differential expression of most genes could be confirmed by semi-quantitative real-time RT-PCR. However, the differential expression of Dock4 (p=0.07) and Myb (p=0.09) just missed significance, and the differential expression of Usp36 and Zmym4 also could not be confirmed. The correlation (Table 4) between the results obtained by the GeneChip® experiment and by real-time PCR was good (R=0.86; p<0.01).
Expression of Egr-1 and cFos mRNA as a function of time
The impact of diffuser and ND filter treatment on Egr-1 and cFos mRNA expression was analyzed over time and also compared to a group of untreated control animals (Figure 3).
The Egr-1 mRNA level in the eyes of treated animals was remarkably reduced when compared to that in the control group. Real-time RT-PCR data were analyzed statistically by ANOVA followed by Dunnett's test (p<0.001). After 30 min, the average Egr-1 mRNA level dropped to about 26% of the control group level. After 4 h of treatment the mean mRNA expression level was about 14% of the initial value, and at the end of the 6+6 h treatment, mean mRNA expression level was approximately 10% of the baseline level. The reduction of image intensity therefore had a higher impact on Egr-1 expression level than did the changes in image contrast and spatial frequency content that were induced by diffuser wear. Nevertheless, the mRNA expression differences (Figure 2) between the diffuser- and ND filter-treated eyes persisted for 6+6 h.
cFos mRNA expression was also found to be decreased by the treatment in comparison to the control group. The real-time RT-PCR data were analyzed statistically by ANOVA followed by Dunnett's test (C versus ND Filter 24 h: p<0.05; C versus all other groups: p<0.001). It was found that after 30 min, the average cFos mRNA level dropped to 35% of the baseline level but increased subsequently. After 4 h the mean mRNA expression level was approximately 44% of the initial value, and after 6+6 h the mean mRNA expression level increased by 17% reaching 61% of the control group level. As in the case of Egr-1, the reduction of image intensity had a high impact on cFos expression, but differences between diffuser- and ND filter-treated eyes were still apparent after 6+6 h (Figure 2).
Baseline expression levels of Akt2, Krt2-6b, Myb, Zmym4, Usp36, and Dock4
Baseline expression level of the genes chosen for validation was evaluated by measuring target gene expression in an untreated control group. The control group was sacrificed after 6 h of light exposure to circumvent diurnal influences on mRNA levels (Figure 1A). Real-time RT-PCR data (Figure 4) were analyzed statistically by ANOVA followed by Tukey-Kramer HSD tests.
The mRNA expression of Akt2 and Dock4 increased in both diffuser- and ND filter-treated eyes, however, the increase was more pronounced in the diffuser-treated eyes. In contrast, Krt2-6b expression was substantially high in untreated animals, while it was lower in the treated animals. Myb expression was hardly detectable in the control group and therefore seems to be induced by the treatment conditions. Usp36 mRNA level seemed to be higher in diffuser- and ND filter-treated eyes, although not statistically significant. Zmym4 mRNA expression was not influenced by either treatment condition.
Semi-quantitative real-time reverse transcriptase polymerase chain reaction
Six potentially interesting genes were selected to validate the diffuser specific expression changes over time: Fgf13, Mapk8ip3, Trak2, Jag1, Stk3, and Klf9 mRNA. Their levels were measured in diffuser-treated eyes only. Baseline expression levels were determined in untreated control animals that spent 6 h in light, including 30 min in the rotating drum. The statistical analysis was done by ANOVA followed by pair-wise Tukey-Kramer HSD post hoc tests (Figure 5).
The microarray results were fully confirmed by real-time RT-PCR for Fgf13 (Figure 5A) and Mapk8ip3 (Figure 5B), in part for Stk3 (Figure 5E), and Trak2 (Figure 5C), but not for Klf9 (Figure 5F), and Jag1 (Figure 5D). Possible reasons for this are addressed in the Discussion section. For the purpose of clarity, not all significant changes are shown in Figure 5B but were as follows: the Mapk8ip3 mRNA level was significantly decreased in the treatment groups when compared to the control group (control versus 30 min: p<0.01; control versus 4 h and 6+6 h: p<0.001).
Molecular networks comprising the differentially regulated genes were generated through the use of Ingenuity Pathways Analysis. The relatively low number of differentially expressed genes could not be consigned to a specific pathway. After 30 min and 4 h of deprivation, most of the genes that changed were involved in cellular development, whereas after 1 day of treatment, they were mostly involved in DNA replication, recombination, nucleic acid metabolism, and small molecule biochemistry.
Placing diffusers over the eyes can artificially induce myopia in animal model. The identification of retinal genes and proteins that are induced or repressed by such conditions can provide new candidates for pharmacological intervention of myopia also in humans. The present study demonstrated that short as well as longer periods of retinal image degradation cause significant changes in the expression level of some genes.
Egr-1 mRNA expression
The influence of form deprivation on Egr-1 mRNA expression  could be confirmed and strengthened in the present study, since Egr-1 expression was decreased in the diffuser-treated eyes at all time points even though the intensity of the retinal images was matched in both eyes. Moreover, Egr-1 expression in the retina was strongly regulated by light. In other species it was shown that in a subset of cells, Egr-1 was not regulated by the light intensity but instead regulated by the sign of defocus (i.e. the glucagonergic amacrine cells in the chicken retina and the GAD65 cells in the macaque retina). In these cells, Egr-1 might induce growth signals independently from the lighting situation. We suggest that there exists a focus-sensitive subpopulation of Egr-1 expressing cells in the murine retina as well.
The early growth response protein 1 (Egr-1)  was first identified as an immediate early gene responsive to growth factors and various differentiation signals and later confirmed as a transcriptional regulatory protein. It is induced in the absence of de novo protein synthesis by mitogens, developmental or differentiation cues, tissue or radiation injury or signals that cause neuronal excitation . The zinc-finger protein Egr-1 is located in the nucleus [32-34] and has numerous target genes , among which are PDGF-A  and PDGF-B , bFGF  and TGF-β1 [37,39]. A potential connection between Egr-1 and myopia was first described in the chick, where it was found that the expression of ZENK correlates with the sign of defocus imposed by lenses in a subset of amacrine cells (AC), specifically the glucagon AC . In the macaque retina it was shown that focus-sensitive immunoreactivity for Egr-1 is induced in a subpopulation of GABAergic amacrine cells (GAD65-immunoreactive cells) .
cFos mRNA expression
The decreased expression of cFos mRNA as a result of reduced image quality was determined by microarray analysis after 4 h of treatment. Furthermore, the analysis by real-time RT-PCR also revealed a down-regulation of cFos mRNA in the diffuser-treated eyes after 30 min and even after 6+6 h. The level of cFos mRNA, in contrast to that of Egr-1 mRNA, increased slightly over time in both diffuser and in ND filter-treated eyes, implying that the two genes are independently regulated when the image contrast is reduced by diffusers. The independent regulation of those two immediate early genes has been previously shown in the chick retina .
cFos is an immediate early gene that belongs to the activator protein (AP)-1 transcription factor family [40-42]. In the retina, Fos-like protein expression has been found to be induced by light in amacrine cells and ganglion cells of dark-adapted rabbits . The activation of cFos in response to light onset, and its circadian regulation, have been described in many studies [44-47]. The role of cFos in myopia research has been previously investigated in the chick retina, where it was found that switching from diffuse blur to focused vision induced cFos protein expression in an amacrine cell subpopulation . Thus changes in cFos gene expression might indicate the activation of retinal interneurons or circuits that mediate the growth responses to well-focused images. Furthermore, the transcriptional activity of c-Jun and cFos can be inhibited by retinoic acid receptors in response to their ligands . Conversely, AP-1 represses the transactivation of retinoid receptors [50,51]. Retinoic acid has been previously implicated in the regulation of eye growth. The synthesis of retinoic acid is regulated in a focus-dependent manner in the choroid and the retina [52-54] and the expression of retinal retinoic acid receptors was found to be increased under form-deprivation conditions .
Akt2 (also called protein kinase B β, PKBβ) is a serine/threonine protein kinase and a downstream effector of the phosphatidylinositol 3-kinase (PI3K). Many physiological effectors and pharmaceuticals, including insulin, insulin-like growth factor, vascular endothelial growth factor, nerve growth factor, carbachol, forskolin and vanadate, are capable of inducing Akt kinase activity, primarily in a PI3K-dependent manner . However, it has been found that Akt is also activated independently of the PI3K by heat shock and hyperosmolarity  as well as oxidative stress and chemical stressors . Akt was first implicated in signal transduction by the demonstration that its kinase activity is induced by growth factors such as platelet derived growth factor (PDGF) and basic fibroblast growth factor (bFGF) [59,60]. Furthermore, it has been shown that Akt2 detaches from the inner surface of the plasma membrane, where it is initially activated, and translocates to the nucleus within 30 min of its activation by growth factors . In the nucleus, Akt isoforms have been hypothesized to phosphorylate and modulate the activity of transcription factors .
In mammals, Akt2 is expressed in most tissues and organs including the retina, but especially in insulin-responsive tissues [62,63]. Previous experiments of our group have shown that insulin might act as a growth stimulator (unpublished data by MF obtained from studies in the chick). Akt2 therefore may be a new and interesting candidate for involvement in signal transduction during myopia development.
Krt2-6b, Myb, and Dock4 expression
Fewer details are known about the other genes that were found to be regulated by deprivation: Krt2-6b (also called mK6b), the expression of which was reduced after 30 min of diffuser treatment and even more after neutral density filter treatment compared to untreated control eyes, encodes an intermediate filament featuring a complex expression pattern. Krt2-6b is constitutively expressed in a variety of internal stratified epithelia but the expression can also be induced by injury and other acute challenges. mRNA is induced in the epidermis as early as 1 h following acute injury or topical application of phorbol esters  or retinoic acid .
Myb (also called c-myb) expression was induced by diffuser treatment and to an even higher extent by filter wear. Myb is a myeloblastosis oncogene acting as a transcriptional transactivator. Gene knockout experiments have shown that genes of the Myb family play an essential role in development . Myb is expressed in ganglion, amacrine, horizontal, and photoreceptor cells of adult mice, suggesting that it might play a role in the physiology of retinal cells . Myb family target genes accomplish diverse functions in cell death, cell adhesion, transcription and signal transduction.
Dock4 is a potent rac activator and an unconventional guanine exchange factor for the Rho family of guanosine triphosphatases (Rho GEF GTPases), as a protein interacting with harmonin, which is expressed in the inner and outer photoreceptor segments. Mutations in the actin bundling and PDZ domain-containing protein harmonin are the causes of Usher syndrome type 1C (USH1C), a syndrome of congenital deafness and progressive blindness .
Data normalization and statistical analysis of microarray experiments
Initially, the data of the 30 min experiment were normalized using MAS5 (Microarray Suite 5.0), resulting in 169 differentially expressed genes (FC 1.5, p <0.01). Since gcRMA normalization is currently recommended as the method of choice by some authors , it was subsequently applied to normalize the data of all experiments. Based on this method, the obtained list contained only 16 genes. The same was true for the 4 h and the 6+6 h experiments, where 173 and 164 differentially expressed genes were obtained when MAS5 normalization was conducted, respectively. In a comparison of the MAS5 and gcRMA generated gene lists, only cFos was found in both lists obtained for the 4 h experiment. The other lists showed no overlap. This result clearly shows that differences in the normalization method may have a dramatic influence on the outcome of microarray studies. However, no clear consensus exists as to which method is best under a given set of circumstances . A major difference between both normalization strategies is the correction for hybridization to the mismatch (MM) probes for that particular probe set as applied by MAS5 while gcRMA alternatively calculates a background adjustment step that ignores the MM intensities but incorporates sequence information from the probes (GC content) .
Although it is obvious that the mode of normalization has a considerable effect, the genes chosen from either list could be confirmed in our case. Both approaches give therefore only an incomplete, and perhaps partially flawed, picture.
It is also important to mention that even though the data were not Benjamini/Hochberg corrected we were nevertheless able to confirm some of the differentially expressed genes which were detected with less rigorous statistical testing. This may suggest that a more pragmatic approach with less statistical rigor, but more repeated experiments, may represent a more successful screening procedure.
The low number of differentially expressed genes might mirror the fact that a limited number of quite specific changes are expected in the retina during eye growth regulation. It seems likely that the circuits that control eye growth have to operate without compromising other retinal and ocular functions, and that they are therefore localized to a sub-set of retinal neurons.
Another reason for the low number of differentially expressed transcripts might also be that, unless a gene is expressed by just one type of cell, the data could conceal opposite changes in expression in two or more cell types. Thus, for example, "zero" difference could result from equal and opposite changes in two cell types, and "increase" could result from a larger increase in one cell type and a lesser decrease in another type of cell.
Confirmation of microarray results
Most of the changes in gene expression detected by microarray and chosen for further study were validated by RT-PCR. Nevertheless, some of the microarray results, particularly the changes in expression over time, were not confirmed. One reason might be that comparing the independent microarray experiments is more complicated, since the preparations and all following procedures of the different treatment periods were done with some time delay. As previously discussed, normalization is a critical point in microarray analysis, and normalizing the data of all 18 GeneChips® together may carry a risk of potential bias.
Another reason for the difficulties of validation could be that some of the differences between expression in treated and control eyes were minute and therefore difficult to confirm.
The identification of genes that are temporally regulated by diffuser treatment would be of interest, even though the confirmation of the results proved to be difficult. We think that the results obtained by real-time RT-PCR are more reliable, considering the issues of normalization.
Unfortunately, there is little information about the genes for which the time kinetics was examined: Longer treatment periods (4 h, 6+6 h diffuser wear) increased Fgf13 mRNA levels in our study. Its exact function is unknown but Fgf13 was also found to be up-regulated after light-induced retinal damage .
Mapk8ip3 (also called Jip3 or JSAP1) mRNA levels were strongly reduced in diffuser-treated eyes. Mapk8ip3 is a scaffold protein and brings together consecutive members of the c-Jun NH2-terminal kinase (JNK) signaling cascade . It is expressed exclusively in neurons and binds preferentially to JNK3 and also to the MAPK kinase 7 .
Conclusion and outlook
The microarray technique is a powerful technology, offering the possibility to assay thousands of genes in one reaction. Each method has its limitation, and it will therefore never be possible to identify all gene expression changes using only one method. There is evidence from the literature that the extent of changes in retinal mRNA levels during deprivation is sometimes only small and therefore in the range of the detection limit of the used methods (for example Northern blots or real-time PCR). We set the cut-off level for the discovery of differentially expressed genes to more than 1.5 fold, which also included the risk of missing some differentially expressed genes.
It is important to keep in mind that changes in the mRNA level are not always translated into changes in protein content and that post-translational modifications of proteins may also be involved in eye growth regulation. Therefore, a histological study of Akt2 and cFos protein expression in the mouse retina following deprivation seems to be interesting and is planned. It might also be interesting to investigate the expression of the Jun proteins (c-Jun, JunB and JunD) since they dimerize with Fos family members (cFos, FosB, Fra-1, and Fra-2) to form transcriptionally active complexes . Moreover, the mRNA expression pattern of some candidate genes identified in this screening could also be studied in the chick, which is easier to handle and responds more reliably with development of myopia when diffusers are attached.
We thank Dr. M. Bonin from the Microarray Facility Tuebingen and M. Papke from the group of Dr. B. Wissinger for sequence analysis of the PCR products. This study was supported by the German Research Council (DFG Fe 450/1-1 and DFG Scha 518/13-1) and the Kerstan Foundation. The authors also thank Prof. Dr. E. Zrenner for support.
1. Dirani M, Chamberlain M, Garoufalis P, Chen C, Guymer RH, Baird PN. Refractive errors in twin studies. Twin Res Hum Genet 2006; 9:566-72.
2. Wallman J, Turkel J, Trachtman J. Extreme myopia produced by modest change in early visual experience. Science 1978; 201:1249-51.
3. Hodos W, Kuenzel WJ. Retinal-image degradation produces ocular enlargement in chicks. Invest Ophthalmol Vis Sci 1984; 25:652-9.
4. Schaeffel F, Howland HC. Visual optics in normal and ametropic chickens. Clinical Vision Sciences 1988; 3:83-98.
5. Sherman SM, Norton TT, Casagrande VA. Myopia in the lid-sutured tree shrew (Tupaia glis). Brain Res 1977; 124:154-7.
6. Troilo D, Judge SJ. Ocular development and visual deprivation myopia in the common marmoset (Callithrix jacchus). Vision Res 1993; 33:1311-24.
7. Smith EL 3rd, Harwerth RS, Crawford ML, von Noorden GK. Observations on the effects of form deprivation on the refractive status of the monkey. Invest Ophthalmol Vis Sci 1987; 28:1236-45.
8. Schaeffel F, Burkhardt E, Howland HC, Williams RW. Measurement of refractive state and deprivation myopia in two strains of mice. Optom Vis Sci 2004; 81:99-110.
9. Faulkner AE, Kim MK, Iuvone PM, Pardue MT. Head-mounted goggles for murine form deprivation myopia. J Neurosci Methods 2007; 161:96-100.
10. Tejedor J, de la Villa P. Refractive changes induced by form deprivation in the mouse eye. Invest Ophthalmol Vis Sci 2003; 44:32-6.
11. Wallman J, Winawer J. Homeostasis of eye growth and the question of myopia. Neuron 2004; 43:447-68.
12. Buck C, Schaeffel F, Simon P, Feldkaemper M. Effects of positive and negative lens treatment on retinal and choroidal glucagon and glucagon receptor mRNA levels in the chicken. Invest Ophthalmol Vis Sci 2004; 45:402-9.
13. Feldkaemper MP, Schaeffel F. Changes in retinal glucagon content during myopia development and effects of glucagon on contrast sensitivity. Invest Ophthalmol Vis Sci 2001; 42:S59.
14. Bitzer M, Feldkaemper M, Schaeffel F. Identification of retinoic acid synthesizing enzymes in the chicken retina and their possible roles during myopia development. Invest Ophthalmol Vis Sci 1999; 40:S848.
15. Tkatchenko AV, Walsh PA, Tkatchenko TV, Gustincich S, Raviola E. Form deprivation modulates retinal neurogenesis in primate experimental myopia. Proc Natl Acad Sci U S A 2006; 103:4681-6.
16. Seltner RL, Stell WK. The effect of vasoactive intestinal peptide on development of form deprivation myopia in the chick: a pharmacological and immunocytochemical study. Vision Res 1995; 35:1265-70.
17. Stone RA, Laties AM, Raviola E, Wiesel TN. Increase in retinal vasoactive intestinal polypeptide after eyelid fusion in primates. Proc Natl Acad Sci U S A 1988; 85:257-60.
18. Escano MF, Fujii S, Sekiya Y, Yamamoto M, Negi A. Expression of Sonic hedgehog and retinal opsin genes in experimentally-induced myopic chick eyes. Exp Eye Res 2000; 71:459-67.
19. Akamatsu S, Fujii S, Escano MF, Ishibashi K, Sekiya Y, Yamamoto M. Altered expression of genes in experimentally induced myopic chick eyes. Jpn J Ophthalmol 2001; 45:137-43.
20. Fischer AJ, McGuire JJ, Schaeffel F, Stell WK. Light- and focus-dependent expression of the transcription factor ZENK in the chick retina. Nat Neurosci 1999; 2:706-12.
21. Brand C, Burkhardt E, Schaeffel F, Choi JW, Feldkaemper MP. Regulation of Egr-1, VIP, and Shh mRNA and Egr-1 protein in the mouse retina by light and image quality. Mol Vis 2005; 11:309-20 <http://www.molvis.org/molvis/v11/a36/>.
22. Zhong X, Ge J, Smith EL 3rd, Stell WK. Image defocus modulates activity of bipolar and amacrine cells in macaque retina. Invest Ophthalmol Vis Sci 2004; 45:2065-74.
23. Prusky GT, Alam NM, Beekman S, Douglas RM. Rapid quantification of adult and developing mouse spatial vision using a virtual optomotor system. Invest Ophthalmol Vis Sci 2004; 45:4611-6.
24. Schmucker C, Seeliger M, Humphries P, Biel M, Schaeffel F. Grating acuity at different luminances in wild-type mice and in mice lacking rod or cone function. Invest Ophthalmol Vis Sci 2005; 46:398-407.
25. Bartmann M, Schaeffel F. A simple mechanism for emmetropization without cues from accommodation or colour. Vision Res 1994; 34:873-6.
26. Grimm C, Wenzel A, Hafezi F, Reme CE. Gene expression in the mouse retina: the effect of damaging light. Mol Vis 2000; 6:252-60 <http://www.molvis.org/molvis/v6/a34/>.
27. Zhijin Wu and Rafael Irizarry, and Robert Gentleman, and Francisco Martinez Murillo and Forrest Spencer, 2004. A Model Based Background Adjustment for Oligonucleotide Expression Arrays. Johns Hopkins University Dept. of Biostatistics Working Paper Series 1001, Berkeley Electronic Press.
28. Tumor Analysis Best Practices Working Group. Expression profiling--best practices for data generation and interpretation in clinical trials. Nat Rev Genet 2004; 5:229-37.
29. Schippert R, Brand C, Schaeffel F, Feldkaemper MP. Changes in scleral MMP-2, TIMP-2 and TGFbeta-2 mRNA expression after imposed myopic and hyperopic defocus in chickens. Exp Eye Res 2006; 82:710-9.
30. Sukhatme VP, Cao XM, Chang LC, Tsai-Morris CH, Stamenkovich D, Ferreira PC, Cohen DR, Edwards SA, Shows TB, Curran T. A, Le Beau MM, Adamson ED. zinc finger-encoding gene coregulated with c-fos during growth and differentiation, and after cellular depolarization. Cell 1988; 53:37-43.
31. Gashler A, Sukhatme VP. Early growth response protein 1 (Egr-1): prototype of a zinc-finger family of transcription factors. Prog Nucleic Acid Res Mol Biol 1995; 50:191-224.
32. Cao XM, Koski RA, Gashler A, McKiernan M, Morris CF, Gaffney R, Hay RV, Sukhatme VP. Identification and characterization of the Egr-1 gene product, a DNA-binding zinc finger protein induced by differentiation and growth signals. Mol Cell Biol 1990; 10:1931-9.
33. Day ML, Fahrner TJ, Aykent S, Milbrandt J. The zinc finger protein NGFI-A exists in both nuclear and cytoplasmic forms in nerve growth factor-stimulated PC12 cells. J Biol Chem 1990; 265:15253-60.
34. Waters CM, Hancock DC, Evan GI. Identification and characterisation of the egr-1 gene product as an inducible, short-lived, nuclear phosphoprotein. Oncogene 1990; 5:669-74.
35. Fu M, Zhu X, Zhang J, Liang J, Lin Y, Zhao L, Ehrengruber MU, Chen YE. Egr-1 target genes in human endothelial cells identified by microarray analysis. Gene 2003; 315:33-41.
36. Wang ZY, Deuel TF. An S1 nuclease-sensitive homopurine/homopyrimidine domain in the PDGF A-chain promoter contains a novel binding site for the growth factor-inducible protein EGR-1. Biochem Biophys Res Commun 1992; 188:433-9.
37. Khachigian LM, Lindner V, Williams AJ, Collins T. Egr-1-induced endothelial gene expression: a common theme in vascular injury. Science 1996; 271:1427-31.
38. Hu RM, Levin ER. Astrocyte growth is regulated by neuropeptides through Tis 8 and basic fibroblast growth factor. J Clin Invest 1994; 93:1820-7.
39. Kim SJ, Glick A, Sporn MB, Roberts AB. Characterization of the promoter region of the human transforming growth factor-beta 1 gene. J Biol Chem 1989; 264:402-8.
40. Abate C, Curran T. Encounters with Fos and Jun on the road to AP-1. Semin Cancer Biol 1990; 1:19-26.
41. Pennypacker KR. AP-1 transcription factor complexes in CNS disorders and development. J Fla Med Assoc 1995; 82:551-4.
42. Curran T, Morgan JI. Fos: an immediate-early transcription factor in neurons. J Neurobiol 1995; 26:403-12.
43. Sagar SM, Sharp FR. Light induces a Fos-like nuclear antigen in retinal neurons. Brain Res Mol Brain Res 1990; 7:17-21.
44. Gudehithlu KP, Neff NH, Hadjiconstantinou M. c-fos and NGFI-A mRNA of rat retina: evidence for light-induced augmentation and a role for cholinergic and glutamate receptors. Brain Res 1993; 631:77-82.
45. Yoshida K, Kawamura K, Imaki J. Differential expression of c-fos mRNA in rat retinal cells: regulation by light/dark cycle. Neuron 1993; 10:1049-54.
46. Guido ME, de Arriba Zerpa GA, Bussolino DF, Caputto BL. Immediate early gene c-fos regulates the synthesis of phospholipids but not of gangliosides. J Neurosci Res 1996; 43:93-8.
47. Bussolino DF, de Arriba Zerpa GA, Grabois VR, Conde CB, Guido ME, Caputto BL. Light affects c-fos expression and phospholipid synthesis in both retinal ganglion cells and photoreceptor cells in an opposite way for each cell type. Brain Res Mol Brain Res 1998; 58:10-5.
48. McGuire JJ, Stell WK. Fos-labelling of retinal neurons activated by complex visual stimuli: Relevance to form-deprivation myopia and ocular growth control. Society for Neuroscience Abstracts. 1998; 24:2091.
49. Pfahl M. Nuclear receptor/AP-1 interaction. Endocr Rev 1993; 14:651-8.
50. Yang-Yen HF, Zhang XK, Graupner G, Tzukerman M, Sakamoto B, Karin M, Pfahl M. Antagonism between retinoic acid receptors and AP-1: implications for tumor promotion and inflammation. New Biol 1991; 3:1206-19.
51. Schule R, Rangarajan P, Yang N, Kliewer S, Ransone LJ, Bolado J, Verma IM, Evans RM. Retinoic acid is a negative regulator of AP-1-responsive genes. Proc Natl Acad Sci U S A 1991; 88:6092-6.
52. Mertz JR, Wallman J. Choroidal retinoic acid synthesis: a possible mediator between refractive error and compensatory eye growth. Exp Eye Res 2000; 70:519-27.
53. Seko Y, Shimokawa H, Tokoro T. In vivo and in vitro association of retinoic acid with form-deprivation myopia in the chick. Exp Eye Res 1996; 63:443-52.
54. Seko Y, Shimizu M, Tokoro T. Retinoic acid increases in the retina of the chick with form deprivation myopia. Ophthalmic Res 1998; 30:361-7.
55. Feldkaemper MP, Wang HY, Schaeffel F. Changes in retinal and choroidal gene expression during development of refractive errors in chicks. Invest Ophthalmol Vis Sci 2000; 41:1623-8.
56. Datta SR, Brunet A, Greenberg ME. Cellular survival: a play in three Akts. Genes Dev 1999; 13:2905-27.
57. Konishi H, Matsuzaki H, Tanaka M, Ono Y, Tokunaga C, Kuroda S, Kikkawa U. Activation of RAC-protein kinase by heat shock and hyperosmolarity stress through a pathway independent of phosphatidylinositol 3-kinase. Proc Natl Acad Sci U S A 1996; 93:7639-43.
58. Konishi H, Matsuzaki H, Tanaka M, Takemura Y, Kuroda S, Ono Y, Kikkawa U. Activation of protein kinase B (Akt/RAC-protein kinase) by cellular stress and its association with heat shock protein Hsp27. FEBS Lett 1997; 410:493-8.
59. Burgering BM, Coffer PJ. Protein kinase B (c-Akt) in phosphatidylinositol-3-OH kinase signal transduction. Nature 1995; 376:599-602.
60. Franke TF, Yang SI, Chan TO, Datta K, Kazlauskas A, Morrison DK, Kaplan DR, Tsichlis PN. The protein kinase encoded by the Akt proto-oncogene is a target of the PDGF-activated phosphatidylinositol 3-kinase. Cell 1995; 81:727-36.
61. Meier R, Alessi DR, Cron P, Andjelkovic M, Hemmings BA. Mitogenic activation, phosphorylation, and nuclear translocation of protein kinase Bbeta. J Biol Chem 1997; 272:30491-7.
62. Altomare DA, Guo K, Cheng JQ, Sonoda G, Walsh K, Testa JR. Cloning, chromosomal localization and expression analysis of the mouse Akt2 oncogene. Oncogene 1995; 11:1055-60.
63. Altomare DA, Lyons GE, Mitsuuchi Y, Cheng JQ, Testa JR. Akt2 mRNA is highly expressed in embryonic brown fat and the AKT2 kinase is activated by insulin. Oncogene 1998; 16:2407-11.
64. Takahashi K, Yan B, Yamanishi K, Imamura S, Coulombe PA. The two functional keratin 6 genes of mouse are differentially regulated and evolved independently from their human orthologs. Genomics 1998; 53:170-83.
65. Rothnagel JA, Seki T, Ogo M, Longley MA, Wojcik SM, Bundman DS, Bickenbach JR, Roop DR. The mouse keratin 6 isoforms are differentially expressed in the hair follicle, footpad, tongue and activated epidermis. Differentiation 1999; 65:119-30.
66. Oh IH, Reddy EP. The myb gene family in cell growth, differentiation and apoptosis. Oncogene 1999; 18:3017-33.
67. Lee E, Chung YH, Park JY, Mun GH, Oh CS, Oh SH, Shin DH. The distribution of c-myb immunoreactivities in the adult mouse retina. Neurosci Lett 2004; 366:297-301.
68. Yan D, Li F, Hall ML, Sage C, Hu WH, Giallourakis C, Upadhyay G, Ouyang XM, Du LL, Bethea JR, Chen ZY, Yajnik V, Liu XZ. An isoform of GTPase regulator DOCK4 localizes to the stereocilia in the inner ear and binds to harmonin (USH1C). J Mol Biol 2006; 357:755-64.
69. Qin LX, Beyer RP, Hudson FN, Linford NJ, Morris DE, Kerr KF. Evaluation of methods for oligonucleotide array data via quantitative real-time PCR. BMC Bioinformatics 2006; 7:23.
70. Shippy R, Fulmer-Smentek S, Jensen RV, Jones WD, Wolber PK, Johnson CD, Pine PS, Boysen C, Guo X, Chudin E, Sun YA, Willey JC, Thierry-Mieg J, Thierry-Mieg D, Setterquist RA, Wilson M, Lucas AB, Novoradovskaya N, Papallo A, Turpaz Y, Baker SC, Warrington JA, Shi L, Herman D. Using RNA sample titrations to assess microarray platform performance and normalization techniques. Nat Biotechnol 2006; 24:1123-31.
71. Huang H, Frank MB, Dozmorov I, Cao W, Cadwell C, Knowlton N, Centola M, Anderson RE. Identification of mouse retinal genes differentially regulated by dim and bright cyclic light rearing. Exp Eye Res 2005; 80:727-39.
72. Raivich G, Behrens A. Role of the AP-1 transcription factor c-Jun in developing, adult and injured brain. Prog Neurobiol 2006; 78:347-63.
73. Kelkar N, Delmotte MH, Weston CR, Barrett T, Sheppard BJ, Flavell RA, Davis RJ. Morphogenesis of the telencephalic commissure requires scaffold protein JNK-interacting protein 3 (JIP3). Proc Natl Acad Sci U S A 2003; 100:9843-8.
74. Angel P, Karin M. The role of Jun, Fos and the AP-1 complex in cell-proliferation and transformation. Biochim Biophys Acta 1991; 1072:129-57.