- Research article
- Open Access
Whole genome assessment of the retinal response to diabetes reveals a progressive neurovascular inflammatory response
BMC Medical Genomics volume 1, Article number: 26 (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.
At harvest, as predicted by earlier assessments, diabetic rats were hyperglycemic and underweight in comparison to control rats (Table 1). These data and the methods described are representative of experiments performed by the Penn State JDRF Animal Models Core since 2003 .
Vascular permeability and apoptosis
Retinal vascular permeability, as measured by FITC-BSA incorporation into nonvascular retinal tissue, was unchanged at 1 month of diabetes but was significantly increased (52%, p < 0.05) after 3 months of diabetes (Table 2). Similarly, caspase-3 activity was unchanged after 1 month of diabetes but was increased significantly after 3 months of diabetes (24%, p < 0.05) (Table 2).
22,578 probes had detectable signal at both the 1 and 3 month time points. Expression values were generally consistent across rats and arrays in the study. Diabetes-associated changes were of a higher number and magnitude at 3 months as compared to 1 month post-STZ injection (Figure 1). The mean fold changes were 1.3 ± 0.6 fold for 1-month changes (124 genes) and 1.7 ± 0.4 fold for the 3 month changes (2310 genes). The full set of gene expression data has been deposited in the NIH/NLM Gene Expression Omnibus (GEO accession # GSE11733) . To determine processes or functional groupings which are altered with diabetes, ontological analysis was performed using Gene Ontology (GO) categories . Ontology results identified a number of physiological processes and molecular functions altered in the diabetic animals (Figure 2). Of special note was the significant induction of the gene ontologies for the physiological processes of cytokine production/immune response at 3 months of diabetes. Also regulated were ontologies related to cytoskeleton/crystallin proteins and GABA receptor activity.
Based on statistical analysis of microarray data and available gene information, including previous reports, 32 genes were chosen for qPCR confirmation and validation (Table 3). These 32 genes represent both genes not previously associated with DR as well as those described with DR or other retinal pathologies, including clinical reports . These genes were placed post hoc in three categories to aid in organization and presentation. These groups are not intended to be definitive or exclusionary nor do they preclude genes from having multiple functions across the different categories and these groupings do not represent canonical gene ontology classes. Nine genes were placed in the neuronal function group, 8 in the microvascular, and 15 in the inflammatory-related group. The 32 genes examined by qPCR are only a subset of the changes observed in the microarray analysis since qPCR confirmation of microarray changes (2401 genes) was not logistically feasible. Other genes/processes which were examined further are described in other reports [17, 18].
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.
Microvascular-related genes demonstrated highly reproducible changes across three independent sets of rats with 5 genes significantly altered at the 3 month time point (EDN2, EDNRB, ICAM1, MCT1, NPPA). For a complete listing of full gene names and synonyms see Table 3. VEGFA reproducibly decreased at 1 month (Figure 3). Two microvascular-related genes (NPR3, VCAM1) did not reach a statistically significant level of change in all sets of rats [See Additional file 1]. Of the 15 inflammatory genes tested, 12 genes (C1-INH, CCR5, CD44, CHI3L1, HSPB1, JAK3, LGALS3, LGALS3BP, LAMA5, PEDF, STAT3, TIMP1) were reproducibly upregulated at 3 months (Figures 4 &5). CCL2 (MCP-1) was increased at both 1 and 3 month time points (Figure 5). Of particular note were commonalities with a recent seminal report of vitreous protein levels in proliferative DR patients (Figure 5) . The magnitude of induction in the inflammatory genes was much larger than other groups with inductions as high as 80 fold. Two inflammatory transcripts were not confirmed as gene expression changes [IL1β, IL6 [See Additional file 2]. All 7 confirmed neuronal function-related genes (DCAMKL1, GAT3, GRIN2A, KCNE2, PCGF1, PEPT2, ZNF219) revealed reproducible decreases at 3 months of STZ-induced diabetes (Figure 6). Two genes (CHRNA4, FEZ2) did not change in a statistically significant manner in all sets of rats [See Additional file 3]. In general, expression ratios were similar between microarray and qPCR measurements. An exact technical comparison is not possible as, in most cases, different animals (Set 1vs Set 2 or 3) and different numbers (N) of animals used in the microarray and qPCR analyses.
To visualize the entire pattern of gene expression measured by qPCR in the different sets of animals the qPCR gene expression data was used to generate a Principle Components Analysis (PCA) plot (Figure 7). This allows for visualization of the relationship between sets of animals at different time points and treatments. The first component provides separation between control and diabetic sets and accounts for 75.6% of the study variance providing further evidence of the reproducibility of the diabetic effect. Furthermore, the 1 month animal sets represent a midpoint between control and 3 month diabetic animals. This demonstrates that while most of these genes are not statistically altered at 1 month, there is a gene expression pattern developing prior to onset of significant increases in vascular permeability and apoptosis. The second component accounted for a much smaller percentage of the study variance (8.5%) and represents the variability between the independent sets of animals. Therefore, these gene expression changes are a robust effect compared to inter-experiment variability.
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.
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.
Mohamed Q, Gillies MC, Wong TY: Management of diabetic retinopathy: a systematic review. JAMA. 2007, 298: 902-916. 10.1001/jama.298.8.902.
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.
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.
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. 2007
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.
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.
Wellen KE, Hotamisligil GS: Inflammation, stress, and diabetes. J Clin Invest. 2005, 115: 1111-1119.
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.
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.
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.
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.
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.
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.
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.
Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088.
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.
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. 2007
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Almkvist J, Karlsson A: Galectins as inflammatory mediators. Glycoconj J. 2004, 19: 575-581. 10.1023/B:GLYC.0000014088.21242.e0.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Xu Q, Qaum T, Adamis AP: Sensitive blood-retinal barrier breakdown quantitation using Evans blue. Invest Ophthalmol Vis Sci. 2001, 42: 789-794.
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.
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.
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.
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.
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.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1755-8794/1/26/prepub
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.
RMB, KMP, and HDVB carried out the microarray and qPCR assays, HDVG and AJB generated samples and contributed to the writing of the manuscript, DAA and C-ML performed the vascular permeability studies, KFL performed the caspase assays, TWG contributed to the experimental design and manuscript preparation, SKB directed the generation of samples, and WMF performed the bioinformatic analyses and prepared the manuscript. All authors read and approved the final manuscript.
Electronic supplementary material
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Brucklacher, R.M., Patel, K.M., VanGuilder, H.D. et al. Whole genome assessment of the retinal response to diabetes reveals a progressive neurovascular inflammatory response. BMC Med Genomics 1, 26 (2008). https://doi.org/10.1186/1755-8794-1-26
- Diabetic Retinopathy
- Vascular Permeability
- Gene Expression Change
- Ingenuity Pathway Analysis
- Proliferative Diabetic Retinopathy