intravitreal injections of dexamethasone ( dex ) and triamcinolone ( taa ) are a mainstay in a clinical retinal practice to treat a wide range of pathologies , such as cystoid macular edema secondary to retinal vein occlusion , diabetic macular edema , choroidal neovascularization in macular degeneration , and noninfectious uveitis [ 16 ] . glucocorticoids have powerful anti - inflammatory action , but the exact underlying mechanism of this effect remains poorly defined . additionally , intravitreal administration of glucocorticoids may lead to side effects such as cataract progression , development of steroid - induced glaucoma , and central serous chorioretinopathy [ 7 , 8 ] . clinical decision making may be greatly enhanced with a better appreciation of the differential gene expression induced by intravitreal steroids . in a previous study using a mouse model , we demonstrated that intravitreal steroids affect the expression of several genetic pathways and alter the balance between neuroprotective and neurodegenerative processes in retina . in this sequel study we used the same model to examine the differential gene expression effects of intravitreal injections of dex and taa in rpe / choroid tissue in vivo , since they too are affected by both disease pathogenesis and intravitreal steroidal treatment at clinically relevant time points of 1 week and 1 month postinjection . the purpose of this study was to evaluate the gene expression effects that intravitreal dex and taa have on the rpe / choroid and the differences between them at clinically relevant time points . this study protocol was approved by the institutional animal care and use committee at the university of california davis before initiation . the study was conducted according to the approved protocol and in accordance with the arvo statement for the use of animals in ophthalmic and vision research . six- to eight - week - old male c57bl/6j mice ( jax - west , sacramento , ca ) were used for all experiments . mice were housed at 21c , under a 12-hour light/12-hour dark cycle , with food and water supplied ad libitum . one ml of ketamine ( 100  mg / ml ) , 0.1  ml of xylazine ( 100  mg / ml ) , and 8.9  ml of sterile water were mixed to prepare a final solution for injection . animals were anesthetized by injecting a 0.1  ml/10  g of body weight intraperitoneally . intravitreal injections were performed transconjunctivally in anesthetized c57bl/6j mice with hamilton 33  g needles delivering 1  l of solution . in group 1 ( n = 4 ) animals received balanced salt solution ( bss ) , group 2 ( n = 4 ) dexamethasone ( 10  g ) , and group 3 ( n = 4 ) taa ( 25  g ) . the solvent for dex and all mice included in this study were injected on the same day at the same time . at weeks 1 and 4 after the injections , the eyes were harvested for ocular dissection and total rna isolation . at days 7 and 30 postinjection mice from each group were sacrificed using a co2 chamber and the rpe / choroid tissue was harvested using published procedures . sample collections were done in the morning , in the period of 1  h. briefly , mouse eyes were removed with curved forceps , the eye glued to the petri dish with superglue , and a circumferential cut was made dorsal to limbus , removing the front of the eye . retinal tissue was removed after severing the optic nerve and the outside of the back of the eye was cleaned of muscles and connective tissue . dissected rpe / choroid was collected within one minute after the sacrifice of the animal , immediately placed in 300  l of rnalater , and stored at 20c . qiagen 's rneasy isolation kit was used to isolate total rna . as a first step tissue was placed in qiazol and pushed through syringe with 17  g needle several times so the pigmented cells of rpe would be released , while sclera stayed intact as a clear sheet and was removed from the tube with needle . since the choroid tends to be firmly attached to the inner sclera , it is likely that most of it stayed attached . rna samples were run on an agilent bioanalyzer microfluidics chip rna nano 6000 to assess quality and quantity . out of four samples , three samples having the highest quality ( rna integrity number ( rin )  7 value ) were labeled as probes for the affymetrix genechip microarrays ( affymetrix , santa clara , ca ) . rna samples were labeled using 0.100  ug of total rna , following the manufacturer 's protocol for genechip target labeling 3 ivt express kit ( affymetrix , santa clara , ca ) . after labeling , this array contains 45,101 probe sets , analyzing the expression of 39,000 transcripts and variants from over 34,000 well - characterized mouse genes . hybridization was performed in the uc davis genome center microarray core facility using a standard procedure ( affymetrix , santa clara , ca ) . microarray data were analyzed using genespring gx 12.5 microarray analysis software and ingenuity pathway analysis ( ipa ) software for biologically relevant changes of expression in genes and related pathways of the rpe / choroid . one - way anova was used in genespring to identify statistically significant genes at a significance level of p  0.01 . the benjamini - hochberg post hoc correction method was too conservative for microarray results with limited number of biological replicas such as this . no probe sets met this significance level of testing . therefore , a different set of stringency requirements was applied to reduce the risk of false positives . to identify biologically relevant gene expression changes for each of the time point / treatment conditions , volcano plot analysis   was performed , using a combination of unpaired t - test and asymptotic p value computation . the volcano plot is a scatter plot of the fold change versus the p value ( in log10p scale ) . it is commonly used to simultaneously depict the p value and the fold change for gene selection cutoff schemes in microarray data . it may be viewed as a summary of  statistical  significance and  biological  significance over a large number of genes . we employed the standard approach of using a p value ( p  0.05 ) as the primary criterion followed by fold change ( 1.5  fc  1.5 ) as the secondary criterion to select differentially expressed genes . this approach ensures control of false - positive error and preserves the desired biological significance . the data discussed in this publication have been deposited in ncbi 's gene expression omnibus [ 13 , 14 ] and are accessible through geo series accession number gse49872 ( http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=%20gse49872 ) . supplementary material ( available online at http://dx.doi.org/10.1155/2016/9856736 ) attached to this paper contains the complete list of the differentially expressed genes identified by anova at p  0.01 and volcano plot at ( 1.5  fc  1.5 ) p  0.05 . pathway analysis was performed using ingenuity pathway analysis ( ipa ) software ( qiagen , valencia , ca ) . the lists of candidate genes identified by genespring analysis were uploaded to the ipa site and ipa core analysis was applied to give us the list of activated pathways and gene networks with the highest degree of significance . quantitative pcr ( qpcr ) was used as a second method for validation of the results of microarray data analysis for several candidate genes ( ahrr , arntl , grik1 , lcn2 , per1 , and snap25 ) . all gene expression assays ( taqman , applied biosystems , foster city , ca ) were conducted in the real - time pcr research and diagnostic core facility at university of california davis . control genes ( b2 m , gapdh , and hprt1 ) were also included in order to validate the results of the qpcr . for the genes used for qpcr validation of the microarray data , threshold cycle ( ct ) values were used to determine the quantity of gene copies . for each gene , the fold change for each gene was determined compared to the mean of all 3 control genes used . raw data was imported as  .cel files and underwent rma normalization with using baseline values of median of all samples . next step was filtering by expression on the raw data for the signals between 20% and 100% , to eliminate too weak signals which were on the borderline of background noise . the signals for the probes were retained only if they were present in at least 1 out of 6 conditions within the range . the number of 36,136 probe sets passed out of total of 45,101 probe sets on the genechip . statistical analysis by one - way anova identified 1,868 out of 36,136 probe sets differentially expressed at a p value of p  0.05 and 430/36,136 at p  0.01 when the entire 18-chip sample set was considered . comparison of treatment versus control for each time point using an unpaired t - test identified a set of common genes similarly regulated by both steroids , as well as , a set of unique genes differentially regulated by each steroid . figures 1(a ) and 1(b ) represent venn diagrams detailing the numbers of genes in each category per each time point for fold change ( fc )  1.5 and p  0.05 . a list of common genes and unique genes for each steroid / time point / regulation group is presented in supplementary tables s1s6 . the full lists of differentially expressed genes by anova at p  0.01 and p  0.05 , as well as t - tests 1.5  fc  1.5 , p  0.05 , can be found in the supplementary tables s7 and s8 . biological networks and pathway analysis was performed utilizing ingenuity pathway analysis ( ipa ) software to identify integral biological pathways that are commonly or uniquely affected by the steroids . statistical analysis ( anova ) of the entire 18-chip sample set identified 430 probe sets differentially expressed at p  0.01 . from this data set the top candidate pathway affected was  circadian rhythm signaling figure 3 represents the canonical ipa pathway  circadian rhythm signaling .  in figure 4 the circadian rhythm signaling pathway is overlaid with expression data from each of the time point / steroid combinations . table 1 lists the genes that were identified by ipa analysis as key members of the pathway and their p values . mammalian circadian clock generates molecular circadian rhythm through coupled transcription / translation feedback loops which involves core clock genes : period ( per ) 1 and per 2 , cryptochrome ( cry ) 1 and cry 2 , clock , and aryl hydrocarbon receptor nuclear translocator - like ( arntl or bmal1 ) genes . recently , it has been shown that the transcription factor npas2 ( mop4 ) has been able to functionally substitute for clock in the mouse forebrain , thus adding the total number to 7 genes . the circadian rhythm generation is tissue specific , and changes of the rate of transcription of mrnas or the stability of per and cry proteins affect clock speed . prokineticin receptor 1 ( prokr1 ) , with an alternative name eg - vegfr1 , is a 7-transmembrane receptor protein from rhodopsin family , an antiapoptotic gene , involved in the signaling pathway in neovascularization . the glutamate receptors ( i.e. , grik1 , grm7 , and grin ) are the predominant excitatory neurotransmitter receptors in the mammalian brain . light sets off pacemakers mediated by glutamate ( glu ) , leading to glutamate receptor activation . of note is that most of the clock proteins have been reported to be localized in both the nucleus and cytoplasm [ 17 , 18 ] . at dex 1 week postinjection per1 , per3 , arntl , clock , cry , and grm7 were upregulated , while per2 , grik1 , grin , and prokr were downregulated ( figure 4(a ) ) . treatment with taa 1 week postinjection resulted in per1 , arntl , and clock upregulated , while grik1 , grm7 , prokr , and partially grin were downregulated ( figure 4(b ) ) . at dex 1 month postinjection most of these genes were downregulated , except arntl , which was still upregulated ( figure 4(c ) ) , while with taa 1 month postinjection most of the genes were downregulated except per2 , which was still upregulated ( figure 4(d ) ) . molecular timekeeping in cells ingenuity pathway analysis is a hypothesis generating tool that helps interpret the data in the context of biological processes , pathways , and networks . while each gene can be interpreted as a one - dimensional representation of data , we can think of a pathway ( linked list of interconnected genes ) as a two - dimensional representation of data , while gene networks represent a multidimensional representation of the data . figure 5 is the multidimensional representation of connections between circadian rhythm pathway with glucocorticoid , dopamine , calcium , gaba , glutamate , and vegf signaling pathways via common genes or molecules . because of that interdependence , regulatory changes in one pathway can reflect in many others , in this case , neurotransmitter signaling pathways , causing widespread consequences , some of which are discussed in section 4 . an alternative analysis by volcano plot in genespring ( 1.5  fc  1.5 , p  0.05 ) taking into account only genes that have the highest fold change was done for each steroid / time point to identify additional signaling pathways unique for each steroid / time point . in this section are highlighted changes in gene expression of genes of interest for ophthalmic diseases , neurodegeneration and neuroprotection , and related signal transduction pathways . the top 6 relevant canonical pathways for each time point with the significance of representation identified by ipa are presented in figures 6(a)6(d ) . the top 16 genes with altered expression in this group are all involved in visual perception and phototransduction . most of these genes were upregulated compared to controls ( snap25 , sag , rho , cnga1 , gnat1 , rcvrn , neurod1 , pde6 g , guca1a , cacna2d4 , pde6b , tulp1 , fscn2 , rdh12 , and olfm3 ) , and one ( grik1 ) was downregulated ( table s1 ) . the most upregulated gene in dex 1 week postinjection was snap25 with a fold change of 12.73 . this protein is regulated by optineurin , a pathogenic gene associated with primary open angle glaucoma ( poag ) . the most affected pathway for taa in week ca signaling has a role in neurotransmitter signaling , molecular transport , and vitamin and mineral metabolism . modulation of migration and ca signaling in retinal pigment epithelium cells has been observed as an initial step of proliferative vitreoretinopathy   and in angiotensin mediated ca signaling during pathologic neovascularization and inflammation . the genes involved in ca signaling that were affected were nachr , nicotinic acetylcholine receptor , which was downregulated , ccl28 , chemokine ( fc = 1.56 ) , trdn , triadin ( fc = 2.09 ) , asph , aspartate - beta - hydroxylase ( fc = 1.87 ) , ptger3 , prostaglandin e receptor 3 ( fc = 2.25 ) , sln , and sarcolipin ( fc = 2.56 ) . interestingly , the most upregulated gene was lcn2 , lipocalin 2 , an innate immunity protein , with a role in cellular response to oxidative stress , drug insult , il-1 response , and tnf - alpha response ( fc = 19.3 ) . the most downregulated ones were crystallins ( cryba1 , crygs and crybb2 , cryaa , cryba4 , crygb , and crygc ) , genes whose protein products play a role in development of the eye , in visual perception , and have a role as heat - shock proteins ( table s2 ) . the genes arntl , per1 , gnb3 , and arrdc2 were upregulated ( fc = 1.552.33 ) , while abra , ptger3 , dpp4 , chrna1 , ifi204 , murc , itsn1 , prokr1 , and ddit4l were downregulated ( fc = 2.42 to 1.56 ) . only three genes were found in common for both dex and taa : myh8 , tpm3 , identified as self - antigen in patients with behcet 's disease with posterior uveitis , and elov1 , with a role in fatty acid biosynthesis . ten genes had significantly altered expression following dex at 1 month postinjection ( table s4 ) . one of the upregulated genes was slc18a2 , a solute carrier with a role in neurotransmitter transport , as well as serotonin and dopamine receptor signaling [ 23 , 24 ] . the second gene was foxs1 , a member of the notch signaling pathway , which has a role in vasculogenesis , as well as insulin receptor and vegf receptor signaling . downregulated was nsg2 , neuron specific gene , part of the dopamine receptor signaling pathway . there were eight affected genes in this group ( table s5 ) , the most interesting being the downregulation of glutamate receptor grm7 . this receptor has a role in neuroprotection , synaptic transmission , and regulation of ca ion transport via voltage gated ca channel activity . as a general observation and contrary to the data of steroid effects in the retina , there were far fewer genes with an altered expression at 1 month than at 1 week postinjection . only one gene was common for both steroids after a t - test with fc = 1.5 , p < 0.05 ( table s6 ) . that gene was fibrillin2 ( fbn2 ) , an extracellular matrix constituent , with ca ion binding capabilities , mildly upregulated for both steroids ( fc = 1.51.6 ) . independent qpcr validation of representative differentially expressed genes was performed using commercial gene expression assays ( taqman , applied biosystems ) . the fold changes and their direction were confirmed for the chosen genes , thus corroborating the microarray data analysis results . rpe dysfunction plays a role in many ocular disorders including age - related macular degeneration , central serous chorioretinopathy , white - dot syndromes , and retinitis pigmentosa , which can lead to permanent visual distortion and blindness . one of the main functions of the rpe is regulation of disc shedding by photoreceptors , which is a rhythmic process essential for photoreceptor health and proper function . the mammalian retina has its own circadian clock mechanism , originally thought to be located in the photoreceptors [ 2931 ] . more recent studies have reported the presence of a circadian clock in rpe cells [ 3234 ] . it has been shown that signaling is mediated by the dopamine d1 receptor and gaba [ 35 , 36 ] . the mammalian retinal circadian clock mechanism is autonomous from the central master circadian clock , which in mammals resides in the hypothalamus . the cyclic pattern of expression of circadian proteins is dependent on a number of external cues or the primary and the strongest zeitgeber is light , which starts the process , and then there are many secondary ones , conveying the information to peripheral tissue . secondary ones start with serotonin and melatonin and then there are more downstream ones , such as insulin and glucocorticoids [ 38 , 39 ] . glucocorticoids have been demonstrated to be implicated in the entrainment of circadian rhythms [ 4042 ] . they were shown to be capable of overcoming other entrainment factors regulating the peripheral oscillators in the kidney and lung   and also act as strong entraining signals for peripheral circadian oscillators and even feedback on central oscillators . in our study this finding suggests that intravitreal injection of steroids has a significant effect on the rpe / choroid , causing dysregulation of the retinal circadian clock , which is observable at 1 week postinjection and continues through 1 month postinjection . disruption of circadian rhythm has been demonstrated to promote inflammation [ 4547 ] and has been implicated as a contributing factor to disease pathogenesis [ 48 , 49 ] . have shown in a diabetic mouse model that dysregulation of circadian rhythmicity resulted in increased retinal inflammation . since both inflammation and the physiological glucocorticoid metabolism are involved with circadian rhythm , it is not surprising that this is a major pathway influenced by steroid administration . although the clinical importance of this finding is not clear , we suggest that the effect of steroids on the expression of circadian rhythm genes may be part of their therapeutic anti - inflammatory action and is not necessarily an adverse effect . on the other hand , it has been reported that the release of glucocorticoid hormones happens in an hourly or ultradian rhythm in a pulsating manner and that levels of free glucocorticoid hormones are highly synchronized between blood , subcutaneous tissue , and the brain . by administering a bolus of steroids , this fine regulation is disrupted , and the loss of the ultradian rhythm may subsequently result in gene expression changes in downstream target tissues . our results demonstrated that the effect of both steroids on circadian rhythm gene expression was much more marked after 1 week than after 4 weeks . this may be reflecting the initial effect following the administration of a single intravitreal dose , which is subsequently waning . accordingly , it is possible that the use of a sustained - release implant , such as ozurdex or iluvien , may achieve a more stable concentration and lead to prolonged effects on rpe / choroid gene expression . these implants have already proven their longer - lasting therapeutic effect in a variety of diseases [ 2 , 4 , 5255 ] . the most affected pathway for dex week 1 postinjection is phototransduction pathway ( figure 6(a ) ) . why genes that belong to the phototransduction pathway change in the rpe / choroid was puzzling . we do not believe that retinal tissue contamination was the cause of this , since mouse retinal tissue is very loosely attached to rpe and upon dissection separates very easily and cleanly , so there is very little , if any , chance of contamination of rpe with retinal cells . to verify whether some of the rod and cone genes are actually expressed in rpe , we checked for the presence of several phototransduction genes listed in the tables from newman et al . ( table s7 ) and found several of them differentially expressed in rpe of amd versus normal ( gngt1 , guca , pde6 , both cone and rod specific , and rgs9 ) . our data suggests that their altered expression levels are a result of the dex treatment of the eye . we hypothesize that apparent upregulation of expression of these genes may have possibly resulted from accumulation of these mrnas due to slowed down phagocytic and autophagic processes that are dysregulated by steroid therapy . another one of the top five activated pathways in dex 1 week postinjection is gaba receptor signaling ( figure 6(a ) ) . the gaba transporter gat , located in the plasma membrane of nerve terminals and glial cells , plays an important role in the termination of synaptic transmission . this neurotransmitter has also been implicated in effecting the circadian rhythm [ 35 , 36 ] , and therefore its activation may reflect the significant alterations demonstrated in the gene expression of this pathway . since both dex and taa were shown to alter the gene expression of the circadian rhythm pathway , it is not surprising that they both altered the gaba receptor signaling pathway as well . the main signal transduction pathways that appear to be activated by steroids in the rpe are glutamate and ca signaling . light sets off pacemakers mediated by glutamate ( glu ) , leading to glutamate receptor activation [ 57 , 58 ] . a recent paper by suwanjang et al . shows how those two pathways are interconnected with glucocorticoid signaling : glucocorticoids reduce the intracellular ca concentration in neurons and astrocytes and protect neurons against glutamate toxicity . l - glutamate is the major excitatory neurotransmitter in the central nervous system that activates both ionotropic ( nmdar ) and metabotropic ( grm ) glutamate receptors . glutamanergic neurotransmission is involved in most aspects of normal neural function , such as intracellular signal transduction , axonal targeting , and synaptic clustering . our results show upregulation of grm7 with dex and downregulation with taa at 1 week . at 1 month grm7 , glutamate receptor 7 , with a role in neuroprotection , synaptic transmission , and regulation of ca ion transport via voltage gated ca channel activity , has been shown to confer susceptibility to many neuropathologic conditions such as age - related hearing impairment [ 57 , 60 , 61 ] , attention - deficit / hyperactivity disorder [ 58 , 62 ] , and schizophrenia [ 6365 ] . two genes that were consistently downregulated at both time points for taa were rbfox1 and neuronal mrna splicing factor fox-1 , which splices mrnas encoding proteins important in synaptic transmission and membrane excitation . a similar protein in the mouse may be involved in regulating the effects of dopamine in the basal ganglia . at dex 1-month postinjection time point , two members of dopaminergic neurotransmission were dysregulated , slc18a2 and nsg2 . our results from differential gene expression of over 34,000 well - characterized mouse genes of the rpe / choroid of c57bl/6j mice after intravitreal steroid injections at 1 week and 1 month postinjection suggest that dysregulation of glutamatergic and dopaminergic neurotransmission signaling may represent a novel target for therapeutic approaches in the treatment of steroid - induced glaucoma . in conclusion , in our previous study , we have reported that intravitreal dex and taa resulted in gene expression changes which potentially altered the balance between neuroprotective and neurodegenerative processes in the retina . in this study , using the same mouse model , we investigated the effects of these steroids on gene expression in the rpe / choroid and found that the major influenced pathways are circadian rhythm and several neurotransmission pathways . comparative analysis of differential gene expression that takes place in different tissues following the administration of dex and taa suggests tissue - specific changes . the activation of signal transduction pathways in the rpe / choroid may mirror or complement the changes in neuroprotection and neurodegeneration processes that occur in the retina following intravitreal steroid administration . the effect on circadian rhythm pathway the effect of both dex and taa was more significant at 1 week than at 4 weeks , which implies that a more controlled slow - release form of administration may achieve longer - lasting cellular and genetic effects . in conclusion , our findings have the potential to contribute to a more complete understanding of the mechanisms by which the steroids dex and taa regulate gene expression in the eye and to the development of improved therapeutics .
purpose . to identify retinal pigment epithelium ( rpe)/choroid genes and their relevant expression pathways affected by intravitreal injections of dexamethasone and triamcinolone acetonide in mice at clinically relevant time points for patient care 
 . methods . 
 differential gene expression of over 34,000 well - characterized mouse genes in the rpe / choroid of 6-week - old c57bl/6j mice was analyzed after intravitreal steroid injections at 1 week and 1 month postinjection , using affymetrix mouse genome 430 2.0 microarrays . 
 the data were analyzed using genespring gx 12.5 and ingenuity pathway analysis ( ipa ) microarray analysis software for biologically relevant changes . 
 results . both triamcinolone and dexamethasone caused differential activation of genes involved in 
 circadian rhythm signaling  pathway at both time points tested . 
 triamcinolone ( taa ) uniquely induced significant changes in gene expression in 
 calcium signaling  ( 1 week ) and  glutamate receptor signaling  pathways ( 1 month ) . 
 in contrast , dexamethasone ( dex ) affected the 
 gaba receptor signaling  ( 1 week ) and  serotonin receptor signaling  ( 1 month ) pathways . 
 understanding how intraocular steroids affect the gene expression of rpe / choroid is clinically relevant . conclusions . 
 this in vivo study has elucidated several genes and pathways that are potentially altering the circadian rhythms and several other neurotransmitter pathways in rpe / choroid during intravitreal steroid injections , which likely has consequences in the dysregulation of rpe function and neurodegeneration of the retina .