Whole genome assessment of the retinal response to diabetes reveals a progressive neurovascular inflammatory response
© Brucklacher et al; licensee BioMed Central Ltd. 2008
Received: 18 February 2008
Accepted: 13 June 2008
Published: 13 June 2008
Despite advances in the understanding of diabetic retinopathy, the nature and time course of molecular changes in the retina with diabetes are incompletely described. This study characterized the functional and molecular phenotype of the retina with increasing durations of diabetes.
Using the streptozotocin-induced rat model of diabetes, levels of retinal permeability, caspase activity, and gene expression were examined after 1 and 3 months of diabetes. Gene expression changes were identified by whole genome microarray and confirmed by qPCR in the same set of animals as used in the microarray analyses and subsequently validated in independent sets of animals. Increased levels of vascular permeability and caspase-3 activity were observed at 3 months of diabetes, but not 1 month. Significantly more and larger magnitude gene expression changes were observed after 3 months than after 1 month of diabetes. Quantitative PCR validation of selected genes related to inflammation, microvasculature and neuronal function confirmed gene expression changes in multiple independent sets of animals.
These changes in permeability, apoptosis, and gene expression provide further evidence of progressive retinal malfunction with increasing duration of diabetes. The specific gene expression changes confirmed in multiple sets of animals indicate that pro-inflammatory, anti-vascular barrier, and neurodegenerative changes occur in tandem with functional increases in apoptosis and vascular permeability. These responses are shared with the clinically documented inflammatory response in diabetic retinopathy suggesting that this model may be used to test anti-inflammatory therapeutics.
Diabetic retinopathy (DR) is one of the most debilitating diabetic complications. DR is best prevented by intensive glycemic control as demonstrated by the Diabetes Control and Complications Trial (DCCT) study . Established treatments of DR are limited to panretinal laser photocoagulation therapy and control of hypertension .
However, the ability to achieve intensive glycemic control through insulin is limited by hypoglycemia, and given that islet cell transplants are not widely available, pharmacotherapies for prevention and treatment of DR are needed. The PKC inhibitor ruboxistaurin has shown promise as a therapeutic agent , but this compound did not receive FDA approval. Therefore, continued development is needed for DR therapies.
The etiology and pathological processes of DR remain imperfectly understood with the major pathophysiolgical hypotheses being microvascular proliferation and/or permeability , advanced glycation endproducts (AGEs) , PKC activation , oxidative stress , inflammation , loss of insulin signaling, and neurodegeneration . Previously, we proposed that these different processes are not mutually exclusive, but rather could represent different components of a feed-forward cycle of nutrient overload leading to chronic inflammation, neurodegeneration, and compromised blood-retinal barrier function .
Previous work has provided important insights into the genomic and proteomic changes associated with DR [10–13]. This study sought to tie changes in retinal gene expression to functional changes in the retina as measured through vascular permeability and apoptosis. Clinical DR is a progressive disease; therefore the rodent model should also demonstrate progressive dysfunction with increasing duration of diabetes. This study accomplished two goals in examining aspects of the retinal response to diabetes: first, to demonstrate functional changes in retinal permeability and caspase-3 activity (as a surrogate for apoptosis) after 1 and 3 months of experimental diabetes, and second, to identify and confirm retinal gene expression alterations at 1 and 3 months of diabetes.
Animal Data. Mean weight and blood glucose at the time of sacrifice for the different sets of animals and their respective analyses. Mean ± S.E.M.
Blood Glucose (mg/dL)
411 ± 7
87 ± 3
Microarray & qPCR
317 ± 16
303 ± 8
617 ± 11
93 ± 4
Microarray & qPCR
303 ± 8
370 ± 21
391 ± 9
114 ± 12
264 ± 16
321 ± 16
587 ± 18
106 ± 3
319 ± 16
335 ± 10
574 ± 8
91 ± 13
364 ± 17
377 ± 33
379 ± 11
90 ± 3
277 ± 9
426 ± 23
545 ± 22
78 ± 2
341 ± 27
475 ± 24
379 ± 11
133 ± 12
296 ± 12
486 ± 21
592 ± 10
116 ± 4
325 ± 23
430 ± 37
Vascular permeability and apoptosis
Caspase and Permeability Data: Retinal vascular permeability and apoptosis was measured at 1 and 3 month time points.
Percent Mean Control
100 ± 8%
98 ± 8%
100 ± 3%
152 ± 8%*
100 ± 9%
100 ± 10%
100 ± 17%
125 ± 16%*
Gene names, aliases, and qPCR assay numbers.
Assay-on-Demand Catalogue #
Endothelin receptor B
Intercellular adhesion molecule 1
Monocarboxylate transporter 1
Natriuretic peptide precursor A
Natriuretic peptide receptor 3
Vascular cell adhesion molecule 1
Vascular endothelial growth factor A
Serping 1, C1NH, C1l
Chemokine (C-C motif) ligand 2
MCP1, GDCF2, MCAF, SCYA2, SMCCF
Chemokine, CC Motif, Receptor 5
CCCKR5, CMKBR5, CKR5
Chitinase 3-Like 1
Heat shock 27 kDa protein 1
HMN2B, HSP27, HSP28, HSP25
Janus Kinase 3
Interleukin 1 beta
BSF2, HGF, HSF, IFNB2, IL-6
Laminin, alpha 5
Lectin, galactoside-binding, soluble, 3
CBP35, GAL3, GALBP, MAC2
lectin, galactoside-binding, soluble, 3, binding protein
MAC2BP, Serum Protein 90K
Pigment epithelium-derived growth factor
Signal transducer and activator of transcription 3
Tissue inhibitor of metalloproteinase 1
Cholinergic receptor, nicotinic, alpha
BFNC, EBN, EBN1, NACRA4
Doublecortin-like kinase 1
Fasciculation and elongation zeta 2
Glutamate receptor, ionotropic, NMDA 2A
Isk-related voltage-gated K+ channel 2
LQT5, LQT6, MGC138292, MIRP1
Polycomb group ring finger 1
Proton-dependent high affinity oligopeptide transporter
GABA transporter 3
Zinc finger protein 219
qPCR confirmation and validation
Quantitative PCR analysis occurred in several steps. First, genes were confirmed for expression changes in the same set (Set 1, Table 1) of animals as used in the microarray analysis. Second, gene expression changes were validated in a second, independent set of animals (Set 2). Samples from the 3 month time point of Set 1 were exhausted before completion of the study and therefore another set (Set 3) of independent animals was generated for the 3 month time point. Genes were considered to be confirmed if the gene expression change was statistically significant at the appropriate time point in all of the animal sets tested.
The purposes of this study were to examine retinal permeability and apoptosis at 1 and 3 months of diabetes and to identify accompanying retinal mRNA expression changes. The time points of this study were chosen as models of the chronic progressive nature of complications such as DR and are consistent with previous studies [14, 19]. The findings of the study demonstrate that changes in retinal expression of inflammatory, microvascular, and neuronal transcripts occur in concert with functional changes in retinal permeability and apoptosis.
The increased vascular permeability reported here (52%) closely replicates our previous findings of increased vascular permeability (62%) at 3 months of diabetes . Similarly, we confirmed our previous finding of increased caspase-3 activity at 3 months of diabetes . We both confirm that these changes are evident at 3 months of diabetes and report that there are no changes in vascular permeability and caspase-3 activity at 1 month by these methods. Lack of significant change in permeability and apoptosis with the assays used at 1 month of diabetes provides the opportunity to assess those gene expression changes which occur before and/or concomitantly with functional alterations.
Thirty two genes significantly altered in the microarray analysis were subjected to qPCR confirmation. After qPCR confirmation of 26 of these genes in multiple sets of animals, these genes were placed in the microvascular, neuronal function and inflammatory gene sets because they represent central hypotheses of retinal cellular alterations with diabetes. While the interactions of the microvascular and neuronal cells along with inflammatory activation are not fully understood, the various retinal cells clearly function in a coordinated fashion with integrated metabolism and cell – cell contacts. We have previously postulated a "feed-forward" interaction concept in which defects in vascular, neural and inflammatory responses to diabetes perpetuate the neurovascular degeneration . These genes were tested across replicate sets of rats to determine those which were reproducibly altered. A major finding of this work is the reproducibility of these changes across independent experiments. Many more differentially regulated genes from the microarray analysis remain to be confirmed and expanded upon [as described elsewhere [17, 18]]. As mRNA levels do not always correspond to protein expression levels, future studies will be needed to address changes in protein levels and precise cellular localization of these changes.
Microvascular dysfunction is a hallmark of DR, with numerous reports of diabetes-induced changes in humans as well as in animal models. Several reproducibly altered microvascular transcripts were confirmed in this study (Figure 3). ICAM1 (intercellular adhesion molecule 1) is commonly associated with retinal disease states, as adhesion molecules are produced by vascular endothelial cells are induced during the proliferative stages of DR . ICAM1 suppression decreases leukostasis and vascular leakage in retinas of STZ-diabetic rats  and ICAM1 deficient STZ treated mice demonstrate reduced leukocyte adhesion and vascular permeability as compared to wild-type .
While the EDN1 (endothelin 1) alterations have been well described in DR, induction of EDN2 (endothelin 2), a potent vasoconstrictor, and the endothelin receptor EDNRB (endothelin receptor B) have been described after 6 months of diabetes and light injury [23, 24]. NPPA (natriuretic peptide precursor type A; also known as atrial natriuretic peptide, ANP) is an antigrowth factor in endothelial cells that counteracts the angiogenic and permeability actions of VEGF. NPPA is also protective against NMDA-induced retinal neurotoxicity . Our findings concur with a recent report of decreased NPPA mRNA and protein levels at 3 months but not 1 month of STZ-diabetes in rats .
Retinal mRNA levels of MCT1 (monocarboxylic acid transporter 1), transporter of lactate and pyruvate between astrocytes and neurons, were previously reported to be unchanged in the rat retina at 10 weeks of STZ-induced diabetes . However, MCT1 mRNA levels were consistently decreased at 3 months of diabetes in this study.
Induction of retinal VEGFA (vascular endothelial growth factor) protein with diabetes is well known, but previous studies of VEGF mRNA have produced conflicting results. For example, retinal VEGF mRNA levels have been reported to be unchanged at 3 months , increased at 6 months , and decreased at 6 months  duration of STZ-induced diabetes in the rat. We demonstrate a consistent decrease in VEGF mRNA expression in multiple independent sets of animals. The lack of concordance with previously described increases in VEGF protein suggests a post-translational mechanism for VEGF protein regulation. In total, these microvascular gene expression differences indicate a synergistic increase in microvascular dysfunction (ICAM1, EDN2, and EDNRB) and decrease in a potential compensatory mechanism (NPPA).
Altered inflammatory equilibrium has the potential to exacerbate or be induced by both the neurodegenerative and microvascular dysfunctions typical of DR . Clinically, elevated vitreal levels of pro-inflammatory cytokines have been described (e.g. [10, 13]). Our findings confirm the induction of many of these same factors in a rodent model and provide evidence of a number of novel inflammation-related changes (Figures 4 &5), some of which have also been described in Müller cells after 6 months of experimental diabetes . The majority of the genes can be described as pro-inflammatory (CCL2, CCR5, CD44, CHI3L1, JAK3, STAT3, TIMP1, LAMA5, LGALS3, and LGALS3BP) with some potentially anti-inflammatory changes (PEDF, HSPB1, C1-INH).
Pro-inflammatory changes included chemokines, receptors, and signaling molecules. All of these genes have multiple points of interaction with inflammatory and chemotactic cascades which themselves are often associated with breakdown of vascular barrier properties and neurodegeneration. Specifically, elevated CCL2 (chemokine, CC motif, ligand 2), (also known as Scya2; Small inducible Cytokine A2, and MCP-1; monocyte chemotactic protein 1) has been reported in the vitreous of patients with proliferative DR , with increased CCL2 correlating with increasing severity of proliferative DR . The high level of induction at 1 month (20.9 fold in set 1 and 5.4 fold in set 2) and even higher elevation at 3 months (84.1 fold in set 2 and 13.0 fold in set 3) demonstrates for the first time that CCL2 mRNA is increased in the retina of STZ-diabetic rats and suggests a very strong pro-inflammatory response that occurs within one month of diabetes and progresses with time.
CCR5 (chemokine, CC motif, receptor 5) a RANTES receptor facilitates leukocyte infiltration into the retina  and could work synergistically with elevated levels of RANTES present with diabetes . CD44 is a widely expressed glycoprotein that plays a direct role in leukocyte trafficking to the retina . The systemic plasma levels of the primary ligand for CD44, hyaluronan, are increased in diabetic patients . In a manner similar to CCR5, increased CD44 levels could work synergistically with elevated levels of its ligand to recruit leukocytes.
JAK3 (janus kinase 3) expression requires cytokine signaling and inhibition of JAK3 prevents STAT3 activation . STAT3 (signal tranducer and activator of transcription 3) has been localized to the inner nuclear layer and ganglion cells . TIMP-1 (tissue inhibitor of metallopeptidase1) expression is regulated through STAT3  and TIMP-1 is induced in the epiretinal membranes of Type 1 and Type 2 diabetic patients .
CHI3L1 (chitinase 3-like 1) is an inflammatory marker associated with acute and chronic inflammation, chemotaxis, and angiogenesis. CHI3L1 plasma levels in Type II diabetes patients are positively correlated with insulin resistance and endothelial dysfunction . LGALS3 (lectin, galactosidase-binding, soluble, 3), also known as GAL3 (galectin 3), plays many roles in inflammatory responses as an advanced glycation end-products (AGE) receptor . AGEs are elevated in diabetic tissues and sera . Removal of LGALS3 inhibits AGE-mediated retinal ischemia and suggests an important role for LGALS3 in AGE-related pathophysiology . LGALS3BP (lectin, galactosidase-binding, soluble, 3, binding protein), is a cell adhesion binding protein for LGALS3 and is elevated in patients with DR . There are no reported data describing the role of LAMA5 (Laminin, alpha-5) retina, however, laminin-LGALS3 interactions can promote leukocyte adhesion .
The three potentially anti-inflammatory gene expression changes confirmed (C1-INH, HSPB1, and PEDF) may represent a compensatory response to aspects of the pro-inflammatory response. C1-INH (complement component 1 inhibitor), also known as Serping1 [serine (or cysteine) peptidase inhibitor, clade G, member1], protein expression is increased in DR patients and C1-INH blocks the increased permeability caused by vitreal injection of carbonic anhydrase . HSPB1 (heat-shock 27-kD protein 1), also commonly known as Hsp27, inhibitis both caspase-dependent and -independent apoptosis  and is protective against oxidative stress breakdown of barrier integrity .
The observation that PEDF (pigment epithelium-derived factor) counteracts the effects of VEGF through antiangiogenic properties  has been challenged by recent report of elevated levels of PEDF in patients with proliferative DR than controls . In fact, increased vitreal levels of PEDF have been described in other DR clinical reports . We also observe an induction of PEDF and more studies will be needed to determine if PEDF induction is a compensatory response to other inflammatory and angiogenic changes or if it has a different purpose.
The majority of these inflammatory changes occur at 3 months of diabetes, the same time as the functional changes. However, the marked induction of CCL2 occurs before observed changes in permeability and apoptosis. This may represent an initial event which leads to the cascade of events eventually resulting in retinal dysfunction. Examination of the entire pattern of gene expression (Figure 7) raises the possibility that small magnitude changes could be working in concert to progressively cause later retinal dysfunction.
We have previously reviewed the neuronal components of DR pathophysiology and the role of neurodegeneration in vision loss . The data presented in this study establish transcriptional alterations in a set of genes related to neuronal function, all suggestive of neurodegeneration. DCAMKL1 (doublecortin-like kinase 1) is highly expressed in the ganglion cell layer  and functions with Doublecortin in microtubule stabilization. Further studies will be needed to determine if a decrease in DCAMKL1 could be associated with reductions or weakening of synapses observed in diabetic rat retinas .
PEPT2 (proton-coupled peptide transporter, Slc15a2) transports peptides across membranes in a low capacity/high affinity manner and is expressed in the inner nuclear layer Müller cells and in astrocytes of the retinal ganglion layer . GAT3 (GABA transporter 3, Slc6a11) is found on both neuronal and glial cells in the retina . The possible role of decreased GABA uptake in retinal pathology is unclear, but could induce electrophysiological changes in the diabetic retina .
KCNE2 (isk-related voltage-gated K+ channel 2) colocalizes with KCNQ channel members to regulate channel activity. While KCNE2 is uncharacterized functionally in the retina, it is highly enriched in rods . GRIN2A (ionotropic glutamate receptor 2a, NMDAR2a) receptors are poorly characterized in the retina, but could play a role in the altered dark adaptation as neurotransmission is increased in the dark .
PCGF1 (polycomb group ring finger 1), also known as nervous system polycomb 1 (NSPc1), is a nervous system transcription repressor regulated by PKC signaling . ZNF219 (zinc finger protein 219, or Zfp219) is an uncharacterized transcriptional repressor  with highest expression in the brain and nervous system.
While existing knowledge of these neuronal function-related genes in the retina is limited, the collective decreases observed indicate a compromised neuronal phenotype. Current work suggests that neurodegeneration in diabetic retinopathy includes a synaptic component that may lead to reduced neuronal function and signaling  further developing the concept of a compromised neuronal phenotype .
Lastly, to visualize the interactions between the different genes confirmed by qPCR, the Ingenuity Pathway Analysis (IPA) system was used to generate a network of gene interactions. Twenty-two of the genes confirmed by qPCR, along with increased caspase activity and decreased insulin levels created a network of interactions [See Additional file 4]. The additional genes and vascular permeability values were not able to be included in a rigorously generated network. This does not preclude them from having roles in the network. While these network analyses can only use known interactions, these increases in expression or activity (red fills) and decreases in expression or levels (green fills) could act synergistically, as we have previously proposed . Furthermore, when previously described roles of these transcripts in functions and diseases (inflammation, diabetes, non-proliferative DR, and neuropathy) are overlaid onto the network, a number of interactions are obvious (orange highlights). As we have discussed above, the combination of these changes could produce a neurovascular inflammatory response.
Each of the genes described has potential roles in the increased vascular permeability and apoptosis observed. A critical finding of this study is that a number of inflammatory genes (C1-INH, CCL2, LGALS3, LGALS3BP, and PEDF) are induced in a manner similar to clinical reports of patients with DR . These inductions occur in tandem with increases in vascular permeability and apoptosis that are analogous to the clinically observed changes in barrier properties and apoptosis. These findings suggest that the STZ-rat could be used to model aspects of clinical DR, especially inflammation. The STZ-rat model may provide a useful preclinical tool to investigate the specific roles of these genes in diabetes-induced retinal pathophysiology. This model may also serve as a preclinical model for studies of potential therapeutics to assess their efficacy in preventing or reversing deregulation of retinal gene expression and pathophysiology with diabetes.
The individual genes identified and confirmed here are points of interest for future studies focused on elucidating the role of each gene product in normal retinal physiology and DR. Together, the panel of genes altered with diabetes provides further evidence to support the hypothesis of a feed-forward cycle of chronic inflammation, neurodegeneration, and compromised blood-retinal barrier function in the pathophysiological development of DR.
Experiments were performed following the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research protocols. All rats were maintained by the Penn State JDRF Animal Core in accordance with the Institutional Animal Care and Use Committee guidelines under specific pathogen-free conditions and monitored by quarterly sentinel testing. Sprague Dawley male rats (Charles River Laboratories, Wilmington, MA) arrived at 100–125 grams. After one week and following an overnight fast, diabetes was induced by intraperitoneal injection of 65 mg/kg streptozotocin (Sigma-Aldrich, St. Louis, MO) in 10 mM sodium citrate pH 4.5 vehicle. Control rats were injected with an equal dose of vehicle only. Rats had free access to food and water, and were on a 12 hour light/dark cycle. Blood glucose level and body weight were measured 6 days post-STZ or vehicle injection, and biweekly throughout the experiment. Only rats with blood glucose levels >250 mg/dL at the time of the original test were included in the diabetic groups (Table 1, Online Appendix). No exogenous insulin was delivered. At the time of retinal harvest, rats were given a lethal dose of pentobarbital, 100 mg/kg, (Ovation Pharmaceuticals Inc., Deerfield, IL) by intraperitoneal injection. Retinas were excised and quickly frozen in liquid nitrogen.
Blood-retinal barrier permeability was measured in rats after 1 and 3 months after streptozotocin (STZ) treatment using a FITC-BSA modification of the method described by Xu et al.  using FITC-BSA. Under ketamine/xylazine anesthesia (67/6.7 mg/kg body weight, i.m.), animals received femoral vein injection of 100 mg/kg body weight of fluorescein isothiocyanate-bovine serum albumin (FITC-BSA) in sterile PBS (Sigma, St. Louis, MO). After 2 hours, the animals were anesthetized again with ketamine/xylazine. The abdominal cavity was opened and 1 ml blood was drawn from the inferior vena cava to obtain the FITC-BSA concentration in plasma. Subsequently, the heart of the animal was perfused for 2 min with citrate buffer (50 mM, pH 3.5, 37°C) containing 1% paraformaldehyde (Fisher, Pittsburgh, PA). After perfusion, retinas were harvested from both eyes and dried in a Savant Speed-Vac (Thermo Scientific, Waltham, MA). The FITC-BSA was extracted by incubating each retina in 200 μl PBS with 1% Triton X-100 (Sigma) and 0.1% sodium azide (Fisher) rocking at room temperature overnight. The extract was transferred to a filter tube and centrifuged at 5,000 × g for 3 hours. Plasma and retina extract samples were assayed in triplicate of 50 μl/well in a 96-well black/clear bottom plate (BD Biosciences, Franklin Lakes, NJ) for fluorescence with a SpectraMax Gemini EM fluorescent plate reader (Molecular Devices, Sunnyvale, CA) based on standard curves of FITC-BSA with excitation at 488 nm and emission at 520 nm. The auto-fluorescence background was subtracted with the retina obtained from the rat without FITC-BSA injection. The retinal permeability was calculated and expressed in μl plasma/g dry retina weight/hr circulation.
Caspase-3 activity was measured in retinal protein homogenates using the fluorometric CaspACE assay system (Promega). Rats were anesthetized with 100 mg/kg sodium pentobarbital. Retinas were excised, and placed in 60 μL of cold lysis buffer (25 mM Hepes, pH7.5, 5 mM MgCl2, 5 mM EDTA, 5 mM DTT, 2 mM PMSF, 10 μg/mL leupeptin, 1% NP40). The retina was gently sonicated followed by 30 minute incubation at 4°C and 20 min centrifugation at 16,000 g at 4°C. Caspase-3 activity was measured in the supernatants using the fluorometric CaspACE assay system (Promega) and protein content was measured by BioRad protein assay. CaspACE assay was performed according to the manufacturer's protocol using 50 μg retinal protein and 37°C incubations in a 96-well plate format.
Microarray analysis was performed in the Penn State College of Medicine Functional Genomics Core Facility according to standard procedures.
Total RNA was isolated with Tri-Reagent/BCP (Molecular Research Center, Cincinnati, OH) following standard methods  and quality and quantity was assessed using the RNA 6000 Nano LabChip with an Agilent 2100 Expert Bioanalyzer (Agilent, Palo Alto, CA).
1 μg RNA was transcribed to cRNA following the protocol of the Codelink iExpress cRNA Prep & Hyb Kit (GE Healthcare). After second strand synthesis and purification of dsDNA using Qiaquick spin columns (Qiagen, Valencia, CA), T7 reaction buffer, NTPs, 10 mM biotin-11-UTP, and T7 polymerase were added to the dsDNA and the reaction was incubated at 37°C for 14 hours. Biotin labeled cRNA was then purified using RNEasy columns (Qiagen) followed by quantitation. 10 μg from each sample was fragmented, denatured and then hybridized to Codelink rat whole genome microarray slides for 18 hours at 37°C. Slides were then incubated at room temperature with Alexa Fluor 647 labeled steptavidin for 30 minutes followed by washing.
Microarrays were scanned with an Axon 4000B scanner with GenePix4 v4.0 software at a 5 μm resolution at 635 nm with laser power at 100%, PMT voltage at 600 V, focus position 0 μm, and lines to average = 1. Images were then imported into CodeLink Expression Analysis Software v4.1 (GE Healthcare). Initial quality control (positive and negative controls), exclusion of manufacturing defects (MSR spots), background subtraction, and intra-array normalization was then performed and the results exported to GeneSpring GX 7.3 (Agilent Technologies).
For data analysis, data was imported into GeneSpring GX 7.3 (Agilent Technologies) and signal values less than 0.01 were set to 0.01, arrays were normalized to the 50th percentile, and individual genes normalized to the median. Values were then normalized on a per gene basis to the control group for each of the two time points (1 or 3 Month). Potential differential expression was determined with a one-way ANOVA (variances not assumed to be equal), p < 0.01 and filtered for 1.2 fold and greater differences in expression in accordance with standards for microarray analysis [59, 60]. The use of a combination of statistical and fold-change cutoffs as opposed to traditional multiple testing corrections (e.g. Bonferroni) produce gene lists with the lowest rate of type I and type II errors. 1.2 fold was chosen as the fold change cutoff, as this magnitude change is at the lower range of changes confirmable by qPCR. Lastly, probe sequences on the array were searched against rat genome sequences to eliminate any probes for sequences removed by standard genome processing and annotation. The full complement of microarray data has been deposited in the NIH/NLM Gene Expression Omnibus (GEO accession # GSE11733) .
Ontological analysis used Gene Ontology (GO) categories . Differentially expressed processes or functional categories were determined statistically, similarly to a previously described approach  using GeneSpring GX software. This analysis determined the number of genes in a category present on the array and the number of expression changes that would be part of that category by random chance given the number of differentially expressed genes.
Pathway and network analysis
Ingenuity Pathway Analysis (Ingenuity Systems, Redwood City CA) was used to create a network from the qPCR confirmed gene expression results. Insulin levels, caspase activity, and vascular permeability data were also included.
Clustering and statistical analyses
For statistical analysis of qPCR data, standard parametric t-tests (α < 0.05, two-tailed) were used. Each set (1 and 3 month) of animals was treated independently since the animals were generated and sacrificed separately. Principal components analysis clustering was performed using a mean centered and scaled algorithm (GeneSpring 7.3, Agilent Technologies).
Quantitative PCR analysis was performed using the 7900 HT Sequence Detection System (Applied Biosystems, Foster City, CA), 384-well optical plates, and Assay-On-Demand (Applied Biosystems) gene specific primers and probes. SDS 2.2.2 software and the 2-ΔΔCtCt analysis method  were used to quantitate relative amounts of product using β-actin as an endogenous control. β-actin levels were determined to be unchanged in an absolute quantitation experiment (data not shown). For a full listing of primer/probe sets see Table 2, Online Appendix. For statistical analysis of qPCR data, standard parametric t-tests (α < 0.05, two-tailed) were used. Each set (1 and 3 month) of animals was treated independently since the animals were generated and sacrificed separately.
This study was supported with research funding from the Juvenile Diabetes Research Foundation to WMF, DAA, KFL, AJB, SKB, and TWG, as well as Pennsylvania Tobacco Settlement Funds to WMF and TWG. TWG is the Jack and Nancy Turner Professor of Ophthalmology. This work was previously presented at Association for Research in Vision and Ophthalmology annual meeting.
The authors would like to thank Dennis Coles for assistance with figure preparation; J. Kyle Krady for samples; Wendy Dunton, Allison Collins, and Sara Mendoza De Granja in the Penn State JDRF Animal Models Core; Daniel J. Krissinger and Alan Kunselman for data analysis work, and Joseph F. Freeman for editing assistance.
- The effect of intensive treatment of diabetes on the development and progression of long-term complications in insulin-dependent diabetes mellitus. The Diabetes Control and Complications Trial Research Group. N Engl J Med. 1993, 329: 977-986. 10.1056/NEJM199309303291401.Google Scholar
- Mohamed Q, Gillies MC, Wong TY: Management of diabetic retinopathy: a systematic review. JAMA. 2007, 298: 902-916. 10.1001/jama.298.8.902.View ArticlePubMedGoogle Scholar
- Aiello LP, Davis MD, Girach A, Kles KA, Milton RC, Sheetz MJ, Vignati L, Zhi XE: Effect of ruboxistaurin on visual loss in patients with diabetic retinopathy. Ophthalmology. 2006, 113: 2221-2230. 10.1016/j.ophtha.2006.07.032.View ArticlePubMedGoogle Scholar
- Ng EW, Shima DT, Calias P, Cunningham ET, Guyer DR, Adamis AP: Pegaptanib, a targeted anti-VEGF aptamer for ocular vascular disease. Nat Rev Drug Discov. 2006, 5: 123-132. 10.1038/nrd1955.View ArticlePubMedGoogle Scholar
- Herold K, Moser B, Chen Y, Zeng S, Yan SF, Ramasamy R, Emond J, Clynes R, Schmidt AM: Receptor for advanced glycation end products (RAGE) in a dash to the rescue: inflammatory signals gone awry in the primal response to stress. J Leukoc Biol. 2007Google Scholar
- Aiello LP: The potential role of PKC beta in diabetic retinopathy and macular edema. Surv Ophthalmol. 2002, 47 Suppl 2: S263-S269. 10.1016/S0039-6257(02)00391-0.View ArticlePubMedGoogle Scholar
- Caldwell RB, Bartoli M, Behzadian MA, El-Remessy AE, Al-Shabrawey M, Platt DH, Liou GI, Caldwell RW: Vascular endothelial growth factor and diabetic retinopathy: role of oxidative stress. Curr Drug Targets. 2005, 6: 511-524. 10.2174/1389450054021981.View ArticlePubMedGoogle Scholar
- Wellen KE, Hotamisligil GS: Inflammation, stress, and diabetes. J Clin Invest. 2005, 115: 1111-1119.View ArticlePubMedPubMed CentralGoogle Scholar
- Antonetti DA, Barber AJ, Bronson SK, Freeman WM, Gardner TW, Jefferson LS, Kester M, Kimball SR, Krady JK, LaNoue KF, Norbury CC, Quinn PG, Sandirasegarane L, Simpson IA: Diabetic retinopathy: seeing beyond glucose-induced microvascular disease. Diabetes. 2006, 55: 2401-2411. 10.2337/db05-1635.View ArticlePubMedGoogle Scholar
- Kim SJ, Kim S, Park J, Lee HK, Park KS, Yu HG, Kim Y: Differential expression of vitreous proteins in proliferative diabetic retinopathy. Curr Eye Res. 2006, 31: 231-240. 10.1080/02713680600557030.View ArticlePubMedGoogle Scholar
- Joussen AM, Huang S, Poulaki V, Camphausen K, Beecken WD, Kirchhof B, Adamis AP: In vivo retinal gene expression in early diabetes. Invest Ophthalmol Vis Sci. 2001, 42: 3047-3057.PubMedGoogle Scholar
- Gerhardinger C, Costa MB, Coulombe MC, Toth I, Hoehn T, Grosu P: Expression of acute-phase response proteins in retinal Muller cells in diabetes. Invest Ophthalmol Vis Sci. 2005, 46: 349-357. 10.1167/iovs.04-0860.View ArticlePubMedGoogle Scholar
- Gao BB, Clermont A, Rook S, Fonda SJ, Srinivasan VJ, Wojtkowski M, Fujimoto JG, Avery RL, Arrigg PG, Bursell SE, Aiello LP, Feener EP: Extracellular carbonic anhydrase mediates hemorrhagic retinal and cerebral vascular permeability through prekallikrein activation. Nat Med. 2007, 13: 181-188. 10.1038/nm1534.View ArticlePubMedGoogle Scholar
- Krady JK, Basu A, Allen CM, Xu Y, LaNoue KF, Gardner TW, Levison SW: Minocycline reduces proinflammatory cytokine expression, microglial activation, and caspase-3 activation in a rodent model of diabetic retinopathy. Diabetes. 2005, 54: 1559-1565. 10.2337/diabetes.54.5.1559.View ArticlePubMedGoogle Scholar
- Edgar R, Domrachev M, Lash AE: Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002, 30: 207-210. 10.1093/nar/30.1.207.View ArticlePubMedPubMed CentralGoogle Scholar
- Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088.View ArticlePubMedGoogle Scholar
- VanGuilder HD, Brucklacher RM, Patel KM, Ellis RW, Freeman WM, Barber AJ: Diabetes downregulates presynaptic proteins and reduces basal synapsin I phosphorylation in rat retina. European Journal of Neuroscience. 2007, Submitted.Google Scholar
- Fort PE, Singh RJ, Losiewicz MK, Freeman WM, Jerfferson LS, Kimball SR, Gardner TW: Proteomic analysis reveals an anti-inflammatory component of diabetic retinopathy. International Ocular Inflammation Society. 2007Google Scholar
- Antonetti DA, Barber AJ, Khin S, Lieth E, Tarbell JM, Gardner TW: Vascular permeability in experimental diabetes is associated with reduced endothelial occludin content: vascular endothelial growth factor decreases occludin in retinal endothelial cells. Penn State Retina Research Group. Diabetes. 1998, 47: 1953-1959. 10.2337/diabetes.47.12.1953.View ArticlePubMedGoogle Scholar
- Limb GA, Chignell AH, Green W, LeRoy F, Dumonde DC: Distribution of TNF alpha and its reactive vascular adhesion molecules in fibrovascular membranes of proliferative diabetic retinopathy. Br J Ophthalmol. 1996, 80: 168-173. 10.1136/bjo.80.2.168.View ArticlePubMedPubMed CentralGoogle Scholar
- Miyamoto K, Khosrof S, Bursell SE, Rohan R, Murata T, Clermont AC, Aiello LP, Ogura Y, Adamis AP: Prevention of leukostasis and vascular leakage in streptozotocin-induced diabetic retinopathy via intercellular adhesion molecule-1 inhibition. Proc Natl Acad Sci U S A. 1999, 96: 10836-10841. 10.1073/pnas.96.19.10836.View ArticlePubMedPubMed CentralGoogle Scholar
- Joussen AM, Poulaki V, Le ML, Koizumi K, Esser C, Janicki H, Schraermeyer U, Kociok N, Fauser S, Kirchhof B, Kern TS, Adamis AP: A central role for inflammation in the pathogenesis of diabetic retinopathy. FASEB J. 2004, 18: 1450-1452.PubMedGoogle Scholar
- Rattner A, Nathans J: The genomic response to retinal disease and injury: evidence for endothelin signaling from photoreceptors to glia. J Neurosci. 2005, 25: 4540-4549. 10.1523/JNEUROSCI.0492-05.2005.View ArticlePubMedGoogle Scholar
- Chakrabarti S, Gan XT, Merry A, Karmazyn M, Sima AA: Augmented retinal endothelin-1, endothelin-3, endothelinA and endothelinB gene expression in chronic diabetes. Curr Eye Res. 1998, 17: 301-307. 10.1076/ceyr.17.3.301.5216.View ArticlePubMedGoogle Scholar
- Kuribayashi K, Kitaoka Y, Kumai T, Munemasa Y, Kitaoka Y, Isenoumi K, Motoki M, Kogo J, Hayashi Y, Kobayashi S, Ueno S: Neuroprotective effect of atrial natriuretic peptide against NMDA-induced neurotoxicity in the rat retina. Brain Res. 2006, 1071: 34-41. 10.1016/j.brainres.2005.11.068.View ArticlePubMedGoogle Scholar
- Rollin R, Mediero A, Fernandez-Cruz A, Fernandez-Durango R: Downregulation of the atrial natriuretic peptide/natriuretic peptide receptor-C system in the early stages of diabetic retinopathy in the rat. Mol Vis. 2005, 11: 216-224.PubMedGoogle Scholar
- Layton CJ, Chidlow G, Casson RJ, Wood JP, Graham M, Osborne NN: Monocarboxylate transporter expression remains unchanged during the development of diabetic retinal neuropathy in the rat. Invest Ophthalmol Vis Sci. 2005, 46: 2878-2885. 10.1167/iovs.04-1458.View ArticlePubMedGoogle Scholar
- Chung HK, Choi SM, Ahn BO, Kwak HH, Kim JH, Kim WB: Efficacy of troxerutin on streptozotocin-induced rat model in the early stage of diabetic retinopathy. Arzneimittelforschung. 2005, 55: 573-580.PubMedGoogle Scholar
- Hammes HP, Lin J, Bretzel RG, Brownlee M, Breier G: Upregulation of the vascular endothelial growth factor/vascular endothelial growth factor receptor system in experimental background diabetic retinopathy of the rat. Diabetes. 1998, 47: 401-406. 10.2337/diabetes.47.3.401.View ArticlePubMedGoogle Scholar
- Hernandez C, Segura RM, Fonollosa A, Carrasco E, Francisco G, Simo R: Interleukin-8, monocyte chemoattractant protein-1 and IL-10 in the vitreous fluid of patients with proliferative diabetic retinopathy. Diabet Med. 2005, 22: 719-722. 10.1111/j.1464-5491.2005.01538.x.View ArticlePubMedGoogle Scholar
- Matsumoto Y, Takahashi M, Ogata M: [Relationship between glycoxidation and cytokines in the vitreous of eyes with diabetic retinopathy]. Nippon Ganka Gakkai Zasshi. 2001, 105: 435-441.PubMedGoogle Scholar
- Crane IJ, Xu H, Wallace C, Manivannan A, Mack M, Liversidge J, Marquez G, Sharp PF, Forrester JV: Involvement of CCR5 in the passage of Th1-type cells across the blood-retina barrier in experimental autoimmune uveitis. J Leukoc Biol. 2006, 79: 435-443. 10.1189/jlb.0305130.View ArticlePubMedGoogle Scholar
- Meleth AD, Agron E, Chan CC, Reed GF, Arora K, Byrnes G, Csaky KG, Ferris FL, Chew EY: Serum inflammatory markers in diabetic retinopathy. Invest Ophthalmol Vis Sci. 2005, 46: 4295-4301. 10.1167/iovs.04-1057.View ArticlePubMedGoogle Scholar
- Xu H, Manivannan A, Liversidge J, Sharp PF, Forrester JV, Crane IJ: Involvement of CD44 in leukocyte trafficking at the blood-retinal barrier. J Leukoc Biol. 2002, 72: 1133-1141.PubMedGoogle Scholar
- Nieuwdorp M, Mooij HL, Kroon J, Atasever B, Spaan JA, Ince C, Holleman F, Diamant M, Heine RJ, Hoekstra JB, Kastelein JJ, Stroes ES, Vink H: Endothelial glycocalyx damage coincides with microalbuminuria in type 1 diabetes. Diabetes. 2006, 55: 1127-1132. 10.2337/diabetes.55.04.06.db05-1619.View ArticlePubMedGoogle Scholar
- O'Shea JJ, Pesu M, Borie DC, Changelian PS: A new modality for immunosuppression: targeting the JAK/STAT pathway. Nat Rev Drug Discov. 2004, 3: 555-564. 10.1038/nrd1441.View ArticlePubMedGoogle Scholar
- Zhang SS, Wei JY, Li C, Barnstable CJ, Fu XY: Expression and activation of STAT proteins during mouse retina development. Exp Eye Res. 2003, 76: 421-431. 10.1016/S0014-4835(03)00002-2.View ArticlePubMedGoogle Scholar
- Yamada E, Tobe T, Yamada H, Okamoto N, Zack DJ, Werb Z, Soloway PD, Campochiaro PA: TIMP-1 promotes VEGF-induced neovascularization in the retina. Histol Histopathol. 2001, 16: 87-97.PubMedGoogle Scholar
- Salzmann J, Limb GA, Khaw PT, Gregor ZJ, Webster L, Chignell AH, Charteris DG: Matrix metalloproteinases and their natural inhibitors in fibrovascular membranes of proliferative diabetic retinopathy. Br J Ophthalmol. 2000, 84: 1091-1096. 10.1136/bjo.84.10.1091.View ArticlePubMedPubMed CentralGoogle Scholar
- Rathcke CN, Vestergaard H: YKL-40, a new inflammatory marker with relation to insulin resistance and with a role in endothelial dysfunction and atherosclerosis. Inflamm Res. 2006, 55: 221-227. 10.1007/s00011-006-0076-y.View ArticlePubMedGoogle Scholar
- McFarlane S, Glenn JV, Lichanska AM, Simpson DA, Stitt AW: Characterisation of the advanced glycation endproduct receptor complex in the retinal pigment epithelium. Br J Ophthalmol. 2005, 89: 107-112. 10.1136/bjo.2004.045914.View ArticlePubMedPubMed CentralGoogle Scholar
- Stitt AW: The role of advanced glycation in the pathogenesis of diabetic retinopathy. Exp Mol Pathol. 2003, 75: 95-108. 10.1016/S0014-4800(03)00035-2.View ArticlePubMedGoogle Scholar
- Stitt AW, McGoldrick C, Rice-McCaldin A, McCance DR, Glenn JV, Hsu DK, Liu FT, Thorpe SR, Gardiner TA: Impaired retinal angiogenesis in diabetes: role of advanced glycation end products and galectin-3. Diabetes. 2005, 54: 785-794. 10.2337/diabetes.54.3.785.View ArticlePubMedGoogle Scholar
- Almkvist J, Karlsson A: Galectins as inflammatory mediators. Glycoconj J. 2004, 19: 575-581. 10.1023/B:GLYC.0000014088.21242.e0.View ArticlePubMedGoogle Scholar
- Whitlock NA, Lindsey K, Agarwal N, Crosson CE, Ma JX: Heat shock protein 27 delays Ca2+-induced cell death in a caspase-dependent and -independent manner in rat retinal ganglion cells. Invest Ophthalmol Vis Sci. 2005, 46: 1085-1091. 10.1167/iovs.04-0042.View ArticlePubMedGoogle Scholar
- Bailey TA, Kanuga N, Romero IA, Greenwood J, Luthert PJ, Cheetham ME: Oxidative stress affects the junctional integrity of retinal pigment epithelial cells. Invest Ophthalmol Vis Sci. 2004, 45: 675-684. 10.1167/iovs.03-0351.View ArticlePubMedGoogle Scholar
- Dawson DW, Volpert OV, Gillis P, Crawford SE, Xu H, Benedict W, Bouck NP: Pigment epithelium-derived factor: a potent inhibitor of angiogenesis. Science. 1999, 285: 245-248. 10.1126/science.285.5425.245.View ArticlePubMedGoogle Scholar
- Duh EJ, Yang HS, Haller JA, De JE, Humayun MS, Gehlbach P, Melia M, Pieramici D, Harlan JB, Campochiaro PA, Zack DJ: Vitreous levels of pigment epithelium-derived factor and vascular endothelial growth factor: implications for ocular angiogenesis. Am J Ophthalmol. 2004, 137: 668-674. 10.1016/S0002-9394(03)01422-3.PubMedGoogle Scholar
- Lin PT, Gleeson JG, Corbo JC, Flanagan L, Walsh CA: DCAMKL1 encodes a protein kinase with homology to doublecortin that regulates microtubule polymerization. J Neurosci. 2000, 20: 9152-9161.PubMedGoogle Scholar
- Berger UV, Hediger MA: Distribution of peptide transporter PEPT2 mRNA in the rat nervous system. Anat Embryol (Berl). 1999, 199: 439-449. 10.1007/s004290050242.View ArticleGoogle Scholar
- Hu M, Bruun A, Ehinger B: Expression of GABA transporter subtypes (GAT1, GAT3) in the adult rabbit retina. Acta Ophthalmol Scand. 1999, 77: 255-260. 10.1034/j.1600-0420.1999.770302.x.View ArticlePubMedGoogle Scholar
- Ishikawa A, Ishiguro S, Tamai M: Changes in GABA metabolism in streptozotocin-induced diabetic rat retinas. Curr Eye Res. 1996, 15: 63-71. 10.3109/02713689609017612.View ArticlePubMedGoogle Scholar
- Blackshaw S, Fraioli RE, Furukawa T, Cepko CL: Comprehensive analysis of photoreceptor gene expression and the identification of candidate retinal disease genes. Cell. 2001, 107: 579-589. 10.1016/S0092-8674(01)00574-8.View ArticlePubMedGoogle Scholar
- Xu Y, Ola MS, Berkich DA, Gardner TW, Barber AJ, Palmieri F, Hutson SM, LaNoue KF: Energy sources for glutamate neurotransmission in the retina: absence of the aspartate/glutamate carrier produces reliance on glycolysis in glia. J Neurochem. 2007, 101: 120-131. 10.1111/j.1471-4159.2006.04349.x.View ArticlePubMedGoogle Scholar
- Gong Y, Wang X, Liu J, Shi L, Yin B, Peng X, Qiang B, Yuan J: NSPc1, a mainly nuclear localized protein of novel PcG family members, has a transcription repression activity related to its PKC phosphorylation site at S183. FEBS Lett. 2005, 579: 115-121. 10.1016/j.febslet.2004.11.056.View ArticlePubMedGoogle Scholar
- Sakai T, Hino K, Wada S, Maeda H: Identification of the DNA binding specificity of the human ZNF219 protein and its function as a transcriptional repressor. DNA Res. 2003, 10: 155-165. 10.1093/dnares/10.4.155.View ArticlePubMedGoogle Scholar
- Xu Q, Qaum T, Adamis AP: Sensitive blood-retinal barrier breakdown quantitation using Evans blue. Invest Ophthalmol Vis Sci. 2001, 42: 789-794.PubMedGoogle Scholar
- Freeman WM, Nader MA, Nader SH, Robertson DJ, Gioia L, Mitchell SM, Daunais JB, Porrino LJ, Friedman DP, Vrana KE: Chronic cocaine-mediated changes in non-human primate nucleus accumbens gene expression. J Neurochem. 2001, 77: 542-549. 10.1046/j.1471-4159.2001.00252.x.View ArticlePubMedGoogle Scholar
- Allison DB, Cui X, Page GP, Sabripour M: Microarray data analysis: from disarray to consolidation and consensus. Nat Rev Genet. 2006, 7: 55-65. 10.1038/nrg1749.View ArticlePubMedGoogle Scholar
- Consortium MAQC: The MicroArray Quality Control (MAQC) project shows inter- and intraplatform reproducibility of gene expression measurements. Nat Biotechnol. 2006, 24: 1151-1161. 10.1038/nbt1239.View ArticleGoogle Scholar
- Harris MA, Clark J, Ireland A, Lomax J, Ashburner M, Foulger R, Eilbeck K, Lewis S, Marshall B, Mungall C, Richter J, Rubin GM, Blake JA, Bult C, Dolan M, Drabkin H, Eppig JT, Hill DP, Ni L, Ringwald M, Balakrishnan R, Cherry JM, Christie KR, Costanzo MC, Dwight SS, Engel S, Fisk DG, Hirschman JE, Hong EL, Nash RS, Sethuraman A, Theesfeld CL, Botstein D, Dolinski K, Feierbach B, Berardini T, Mundodi S, Rhee SY, Apweiler R, Barrell D, Camon E, Dimmer E, Lee V, Chisholm R, Gaudet P, Kibbe W, Kishore R, Schwarz EM, Sternberg P, Gwinn M, Hannick L, Wortman J, Berriman M, Wood V, de la CN, Tonellato P, Jaiswal P, Seigfried T, White R: The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. 2004, 32: D258-D261. 10.1093/nar/gkh066.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001, 25: 402-408. 10.1006/meth.2001.1262.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/1/26/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.