Abstract Purpose Endotoxin-induced uveitis (EIU) is a well-established mouse model for studying human acute inflammatory uveitis. The purpose of this study is to investigate the genome-wide retinal transcriptome profile of EIU. Methods The anterior segment of the mice was examined with a slit-lamp, and clinical scores were evaluated simultaneously. The histological changes in the posterior segment of the eyes were evaluated with hematoxylin and eosin (H&E) staining. A high throughput RNA sequencing (RNA-seq) strategy using the Illumina Hiseq 2500 platform was applied to characterize the retinal transcriptome profile from lipopolysaccharide (LPS)-treated and untreated mice. The validation of the differentially expressed genes (DEGs) was analyzed with real-time PCR. Results At the 24th hour after challenge, the clinical score of the LPS group was significantly higher (3.83±0.75, mean ± standard deviation [SD]) than that of the control group (0.08±0.20, mean ± SD; p<0.001). The histological evaluation showed a large number of inflammatory cells infiltrated into the vitreous cavity in the LPS group compared with the control group. A total of 478 DEGs were identified with RNA-seq. Among these genes, 406 were upregulated and 72 were downregulated in the LPS group. Gene Ontology (GO) enrichment showed three significantly enriched upregulated terms. Twenty-one upregulated and seven downregulated pathways were remarkably enriched by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment. Eleven inflammatory response–, complement system–, fibrinolytic system–, and cell stress–related genes were validated to show similar results as the RNA-seq. Conclusions We first reported the retinal transcriptome profile of the EIU mouse with RNA-seq. The results indicate that the abnormal changes in the inflammatory response–, complement system–, fibrinolytic system–, and cell stress–related genes occurred concurrently in EIU. These genes may play an important role in the pathogenesis of EIU. This study will lead to a better understanding of the underlying mechanisms and shed light on discovering novel therapeutic targets for ocular inflammation. Introduction Uveitis is a common inflammatory ocular disease that compromises the vision of patients and even causes blindness within the working population worldwide [[32]1]. Several etiologies that could lead to uveitis have been proposed, including trauma, infection, and autoimmune disease [[33]2]. Uveitis can act either as an isolated intraocular inflammation or as a part of systemic autoimmune diseases [[34]3,[35]4]. Currently, the conventional treatment for uveitis includes topical or systemic administration of corticosteroids, non-steroidal anti-inflammatory drugs (NSAIDs), and immunosuppressive agents. Despite these available treatments, some patients do not respond well to these drugs; consequently, the inflammation becomes refractory or chronic. Even if the drugs are highly efficacious to patients, long-term application may bring about numerous adverse effects, including cataract, glaucoma, metabolic disorders, and even liver or kidney failure [[36]5-[37]7]. Therefore, identification of novel and effective therapeutic targets with fewer side effects for the treatment of uveitis is highly desired. Involving inflammation of the uveal tract, endotoxin-induced uveitis (EIU) is a widely used animal model that shares many clinical and histological characteristics with human acute infectious uveitis. EIU can be induced via a single systemic or local injection of lipopolysaccharide (LPS) [[38]7]. The process of inflammatory response induced by LPS is well understood. LPS stimulation activates Toll-like receptor 4 (TLR4) to initiate the cellular inflammatory pathways, including the transcription factors nuclear factor (NF)-κB and activator protein (AP)-1, and results in the release of cytokines and chemokines and increased expression of adhesion molecules [[39]8-[40]10]. Identification of the abnormally changed genes or molecular pathways may discover potential therapeutic targets, which can be facilitated with high-throughput systems analyses, particularly at the transcriptome levels. The transcriptome is an integrated set of transcripts at a specific stage or under a given physiologic condition [[41]11]. Transcriptome analysis is one of the approaches that aims to identify the genetic variations of a disease and investigate its pathogenesis. Next-generation sequencing (NGS) technology has launched a new era of tremendous potential and applications in genomic and transcriptomic analyses [[42]12,[43]13]. RNA-sequencing with NGS (RNA-seq) provides a comprehensive evaluation and quantification of all subtypes of RNA molecules expressed in a cell or tissue [[44]14]. RNA-seq technology has many advantages compared with earlier technologies. First, it is highly reproducible and has a much greater dynamic range than the microarray analysis. Second, RNA-seq can detect transcripts expressed at low levels [[45]15]. Within the constraints of the coverage depth, RNA-seq can measure the level of any transcripts present in the library and permit the identification of unannotated transcripts, different splice variants, and non-coding RNAs [[46]14,[47]16-[48]19]. Therefore, RNA-seq is expected to reveal a better representation of the transcriptome. Moreover, with a steady reduction in the costs of NGS, RNA-seq is now emerging as a method of choice for comprehensive transcriptome profiling. Accumulating evidence indicates that genome-wide characterization of the retinal transcriptome is essential to understanding the cell development, physiology, as well as disease of the eye [[49]20-[50]22]. Thus, RNA-seq technology may have the potential to reveal comprehensive and detailed information about the transcriptome changes during a disease. In the present study, we used NGS-based RNA-seq analysis to investigate the retinal transcriptome profile of EIU and control mice. We identified a series of differentially expressed genes. This study illustrates that RNA-seq offers a more complete, more accurate, and faster approach for comparative and comprehensive analysis of retinal transcriptomes in EIU and may open up new avenues for discovering new therapeutic targets for ocular inflammation. Methods Ethics statement The study was performed according to the ARVO Statement for the Use of Animals in Ophthalmic and Vision Research. The protocols were approved by the Ethics Committee of the First Affiliated Hospital of Chongqing Medical University. All surgeries were performed under anesthesia, and all efforts were made to minimize animal discomfort and stress. Animals and experimental procedures Female BALB/c mice at the age of 6 to 8 weeks were purchased from Jackson Laboratories (Bar Harbor, ME) and were housed under specific pathogen-free conditions at the Animal Care Service of Chongqing Medical University with a 12 h:12 h light-dark cycle. The mice were randomly divided into two groups: the LPS group and the control group. For the LPs group, to induce EIU, each mouse eye received a single intravitreal injection of 125 ng LPS prepared in 1 μl sterile PBS (140 mM NaCl, 2.7 mM KCl, 10 mM Na[2]HPO[4], 2 mM KH[2]PO[4], PH 7.4) as previously described [[51]23]. For the control group, each mouse eye was given an injection of 1 μl sterile PBS. Inflammation evaluation in the anterior chamber Inflammation in the anterior chamber was evaluated with a slit-lamp (Shangbang, Chongqing, China) observation at the 24th hour after LPS administration. The mice were euthanized by receiving an intraperitoneal injection with the mixture of 75 mg/kg (ketamine) and 13.6 mg/kg (xylazine). The clinical severity of ocular inflammation was graded by two independent observers according to criteria described previously [[52]23]. The clinical severity of ocular inflammation was graded according to the severity of iris hyperemia, exudate in anterior chamber, hypopyon and synechia by two independent observers according to the criteria described previously. Histological evaluation To evaluate the histological changes in the posterior segment, the eyes of the mice were then enucleated at the 24th hour after LPS injection and fixed in 4% paraformaldehyde (PFA) at 4 °C overnight. Eyeballs were embedded in paraffin. Serial 4 μm sections were cut through the cornea-optic nerve axis and stained with hematoxylin and eosin (H&E). Then the infiltrating cells in the posterior chamber and the retina were examined under a light microscope. Library preparation and mRNA sequencing The RNA-sequencing were performed by Novogene (Beijing, China). The library was prepared for sequencing using standard Illumina protocols. Total RNAs were isolated from the retinas of the LPS group and the control group using TRIzol reagent (Invitrogen, Carlsbad, CA) and treated with RNase-free DNase I (New England Biolabs, Ipswich, MA), to remove any contaminating genomic DNA. A total amount of 3 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using the NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (New England Biolabs) following the manufacturer’s instruction. mRNA extraction was performed using Dynabeads oligo (dT; Invitrogen Dynal). Double-stranded cDNAs were synthesized using Superscript II reverse transcriptase (Invitrogen) and random hexamer primers. The cDNAs were then fragmented with nebulization, and the standard Illumina protocol was followed thereafter to create the RNA-seq library. The library quality was assessed with an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA). The cDNA libraries were sequenced by the Illumina Hiseq 2500 device. For the data analysis, base calls were performed using CASAVA. Reads were aligned to the genome using the split read aligner [53]TopHat (v2.0.12) with the mismatches set at 2 and with the default parameters. Differentially expressed gene analysis Differentially expressed gene (DEG) analysis of the LPS group and the control group (each group has three biologic replicates per condition) was performed using the DESeq R package (1.10.1). DESeq provides statistical routines for determining differential expression in digital gene expression data using a model based on the negative binomial distribution. The p values were adjusted using the Benjamini and Hochberg approach for controlling the false discovery rate (FDR). Genes with an adjusted p value of less than 0.05 found by DESeq were assigned as differentially expressed. GO and KEGG pathways enrichment analysis of differentially expressed genes Gene Ontology (GO) enrichment analysis of the DEGs was performed by the GOseq R package, in which gene length bias was corrected. GO terms with a corrected p value of less than 0.05 were considered significantly enriched by differential expressed genes. Pathway enrichment was determined using the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway annotation. KEGG is a database resource for understanding high-level functions and utilities of the biologic system, such as the cell, the organism, and the ecosystem, from molecular-level information, especially large-scale molecular data sets generated by genome sequencing and other high-through put experimental technologies ([54]KEGG). We used KOBAS software (v 2.0) to test the statistical enrichment of differential expression genes in the KEGG pathways. Pathways were considered significantly enriched with an adjusted p value of less than 0.05. Real-time quantitative PCR Mice were euthanized at the 24th hour after LPS injection, and the retinas were freshly microdissected under a dissecting microscope to separate the neuroretina and the RPE-choroid complex. Total RNA of the neuroretina was then extracted using TRIzol reagent (Invitrogen) according to the manufacturer’s instructions. cDNA was synthesized by using the PrimeScript RT reagent kit (Takara Biotechnology, Dalian, China). Real-time PCR analysis of IRF8 was performed in a volume of 20 μl using SYBR Premix Ex Taq^TM II (Takara Biotechnology) with the ABI Prism 7500 system (Applied Biosystems, Foster City, CA). The conditions were 95 °C for 10 min, followed by 40 cycles of 15 s at 95 °C and 60 s at 60 °C. Each reaction was run in duplicate. Data were normalized to the expression of GAPDH. Relative quantification was achieved with the comparative 2^-ΔΔCt method as described previously [[55]23]. The sequences of the PCR primer pairs used in this study are shown in [56]Table 1. Table 1. Sequences of the primers for real-time PCR. Gene Accession number Forward Reverse Ccl5 __________________________________________________________________ [57]NM_013653.3 __________________________________________________________________ 5′-CTGCTGCTTTGCCTACCTCT-3′ __________________________________________________________________ 5′-ACACACTTGGCGGTTCCTT-3′ __________________________________________________________________ A2m __________________________________________________________________ [58]NM_175628.3 __________________________________________________________________ 5′-TCTGCCTACATCACCATTGC-3′ __________________________________________________________________ 5′-CTCTCCTCGCTGACTTCCAG-3′ __________________________________________________________________ Jak3 __________________________________________________________________ [59]NM_001190830.1 __________________________________________________________________ 5′-TGCCATCCTTGACTTACACG-3′ __________________________________________________________________ 5′-GGAACTCTCCCTGCTCCTTC-3′ __________________________________________________________________ Egr1 __________________________________________________________________ [60]NM_007913.5 __________________________________________________________________ 5′-AACACTTTGTGGCCTGAACC-3′ __________________________________________________________________ 5′-AGGCAGAGGAAGACGATGAA-3′ __________________________________________________________________ Atf3 __________________________________________________________________ [61]NM_007498.3 __________________________________________________________________ 5′-CCAGGTCTCTGCCTCAGAAG-3′ __________________________________________________________________ 5′-GATGGCGAATCTCAGCTCTT-3′ __________________________________________________________________ Cdkn1a __________________________________________________________________ [62]NM_001111099.2 __________________________________________________________________ 5′-CAAAGTGTGCCGTTGTCTCTT-3′ __________________________________________________________________ 5′-TCAAAGTTCCACCGTTCTCG-3′ __________________________________________________________________ Gadd45b __________________________________________________________________ [63]NM_008655.1 __________________________________________________________________ 5′-GTCGTTCTGCTGCGACAAT-3′ __________________________________________________________________ 5′-TGACAGTTCGTGACCAGGAG-3′ __________________________________________________________________ C1ra __________________________________________________________________ [64]NM_023143.3 __________________________________________________________________ 5′-CTGGAAGCTGCACTACACCA-3′ __________________________________________________________________ 5′-TGAAGTAATCCCGGAACTGG-3′ __________________________________________________________________ Sell __________________________________________________________________ [65]XM_006496718.3 __________________________________________________________________ 5′-CCAGTCCAAGTGTGCTTTCA-3′ __________________________________________________________________ 5′-TCTCTTGGCAGATTGGCTCT-3′ __________________________________________________________________ C3 __________________________________________________________________ [66]NM_009778.3 __________________________________________________________________ 5′-TGGTGGAGAAAGCAGTGATG-3′ __________________________________________________________________ 5′-ACGGGCAGTAGGTTGTTGTC-3′ __________________________________________________________________ F13a1 __________________________________________________________________ [67]NM_001166391.1 __________________________________________________________________ 5′-GCCCAATAACTCCAATGCTG-3′ __________________________________________________________________ 5′-CCCATCTCTCCTTGAACAGG-3′ __________________________________________________________________ GAPDH [68]XM_017321385.1 5′-GTATGACTCCACTCACGGCAAA-3′ 5′-GGTCTCGCTCCTGGAAGATG-3′ [69]Open in a new tab Statistical analysis Data are shown as mean±standard error of the mean (SEM). Real-time PCR data were analyzed with GraphPad Prism 5 software (GraphPad Software, Inc., San Diego, CA). The statistical significance of differences was determined with the independent-sample t test. A p value of less than 0.05 was considered statistically significant. Results Inflammation signs and clinical scores of the anterior chamber Inflammation in the anterior chamber was assessed with a slit-lamp at 24 hours after LPS injection. The clinical scores were evaluated based on the criteria described previously [[70]23] simultaneously with the slit-lamp. There was significant iris hyperemia and exudation into the anterior chamber, and hypopyon was observed in the LPS group. No inflammation was observed in the control group at the 24th hour. Evaluation of the clinical score showed that at the 24th hour, the clinical score for the mice in the LPS group was statistically significantly higher (3.83±0.75, mean ± SD) than that for the control group (0.08±0.20, mean ± SD; ***p<0.001; [71]Figure 1). Figure 1. [72]Figure 1 [73]Open in a new tab Inflammatory signs in the anterior chamber after LPS administration. The ocular inflammation was assessed with slit-lamp microscopy at the 24th hour after lipopolysaccharide (LPS) injection. Left panel: Representative images of the anterior chamber from the control- and LPS-treated groups at the 24th hour after the LPS injection. Right panel: The clinical score of the anterior chamber were evaluated at the 24th hour after LPS administration using the criteria of the endotoxin-induced uveitis (EIU) clinical score. The clinical score of ocular inflammation was assessed by two independent observers according to the criteria (n = 6). Histological changes in the posterior chamber To evaluate the histological changes, eyes were collected at the 24th hour after LPS injection. Histological examination showed that no inflammatory cells had infiltrated into the ciliary body or the vitreous cavity in the control group ([74]Figure 2A,C). However, severe intraocular inflammation was observed in the LPS group as evidenced by a large number of inflammatory cells infiltrated into the vitreous cavity and the ciliary body ([75]Figure 2B,D). Figure 2. [76]Figure 2 [77]Open in a new tab Histological evaluation of the eyes from the mice treated with or without LPS. Representative images (A, C) show infiltrating inflammatory cells in the posterior chamber of the eyes in the control group at the 24th hour. Representative images (B, D) show infiltrating inflammatory cells in the posterior chamber of the eyes in the lipopolysaccharide (LPS) group at the 24th hour after LPS injection. Magnification = 10X for A and B and 200X for C and D (n = 3). Hierarchical cluster of global gene expression To examine the global gene expression in the LPS-treated tissue and the control retinal tissue, mRNA was extracted from three individual biologic replicates, and RNA-seq was performed. To compare the difference in the global gene expression profiles of the different samples, average linkage hierarchical cluster analysis was performed ([78]Figure 3). The LPS-treated retina samples clearly separated from the control retina samples, suggesting that there is a distinct difference in gene expression between control and LPS-treated retinas. Figure 3. [79]Figure 3 [80]Open in a new tab Hierarchical cluster analysis of the retinas between the EIU mice and the control mice. The cluster analysis displays gene expression changes in the retinas of the lipopolysaccharide (LPS)-treated mice and the control mice. Gene expression profiles of the control mice that received sham injection with PBS in each eye (C1–3, n = 3) were compared with the endotoxin-induced uveitis (EIU) mice that received a single intravitreal injection with LPS in each eye (L1–3, n = 3). Red bars indicate increased expression, and blue bars indicate decreased expression. DEGs in the retinas of the EIU mice and the control mice To investigate the changes in gene expression after LPS administration, computational analysis of the DEGs between the LPS group and the control group was performed. Genes were considered to be differentially expressed with an adjusted p value of less than 0.05. A total of 478 genes exhibited altered expression in the LPS group versus the control group. Among these DEGs, 406 (85%) genes were upregulated in the LPS group ([81]Figure 4, red dots; Appendix 1), and 72 genes (15%) were downregulated ([82]Figure 4, green dots; Appendix 2). Figure 4. [83]Figure 4 [84]Open in a new tab Volcano plot of the upregulated and downregulated DEGs. The volcano plot shows the upregulated and downregulated differentially expressed genes (DEGs) of the retinas in the lipopolysaccharide (LPS) group and the control group. For each plot, the x-axis represents the log[2] fold change (FC), and the y-axis represents -log[10] (p values). Genes with an adjusted p value of less than 0.05 found with DESeq were assigned as differentially expressed. GO enrichment analysis of DEGs To investigate the GO, including the biological process, molecular function, and cellular location of the differentially expressed genes, GO enrichment analysis was applied to analyze all the DEGs in the LPS group compared with the control group. The top 30 most enriched upregulated and the top 30 most enriched downregulated GO terms are shown in [85]Figure 5. The adjusted p values of the three upregulated GO terms were less than 0.05; thus, these GO terms were considered significantly enriched GO terms, including antigen processing and presentation, immune response, and immune system process ([86]Table 2). These significantly upregulated GO terms were heavily associated with immune-related GO terms, indicating there was an increase in the immune response in the LPS-treated retinas. No significantly downregulated GO terms were enriched ([87]Figure 5). Figure 5. [88]Figure 5 [89]Open in a new tab The top 30 enriched upregulated and the top 30 enriched downregulated GO terms in the LPS group compared with the control group. The top 30 upregulated and the top 30 downregulated Gene Ontology (GO) terms in the lipopolysaccharide (LPS) group compared to the control group. Enriched GO terms with an adjusted p value of less than 0.05 were considered significantly enriched (*adjusted p<0.05). Table 2. The significantly enriched GO terms. Up/down-regulated GO accession Description Term type Adjusted p-value Upregulated GO Terms GO:0019882 __________________________________________________________________ antigen processing and presentation __________________________________________________________________ biologic process __________________________________________________________________ 3.38E-08 __________________________________________________________________ GO:0006955 __________________________________________________________________ immune response __________________________________________________________________ biologic process __________________________________________________________________ 3.11E-05 __________________________________________________________________ GO:0002376 immune system process biologic process 7.81E-05 [90]Open in a new tab KEGG pathway analysis of DEGs in the LPS group compared with the control group To identify the pathways involved in EIU, KEGG pathway enrichment analysis were applied. The top 20 upregulated pathways and the top 20 downregulated pathways are shown in [91]Figure 6. A total of 21 significantly enriched upregulated pathways were found, including phagosome, viral myocarditis, staphylococcus aureus infection, herpes simplex infection, HTLV-I infection, allograft rejection, graft-versus-host disease, type I diabetes mellitus, viral carcinogenesis, autoimmune thyroid disease, complement and coagulation cascades, antigen processing and presentation, Epstein-Barr virus infection, cell adhesion molecules (CAMs), measles, osteoclast differentiation, Chagas disease (American trypanosomiasis), legionellosis, pertussis, the p53 signaling pathway, and the tumor necrosis factor (TNF) signaling pathway. Six significantly downregulated pathways were enriched, including platelet activation, nitrogen metabolism, gastric acid secretion, bacterial invasion of epithelial cells, protein digestion and absorption, and glyoxylate and dicarboxylate metabolism ([92]Table 3). Figure 6. [93]Figure 6 [94]Open in a new tab The top 20 enriched upregulated and the top 20 enriched downregulated KEGG pathways in the LPS group compared with the control group. The top 20 upregulated and the top 20 downregulated Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways of the differentially expressed genes (DEGs) in the lipopolysaccharide (LPS) group compared with the control group. Enriched pathways with an adjusted p value of less than 0.05 were considered significantly enriched (*adjusted p<0.05). Table 3. The significantly enriched KEGG pathways. Up/down-regulated ID Pathway terms P-Value Corrected P-Value Upregulated pathways __________________________________________________________________ mmu04145 __________________________________________________________________ Phagosome __________________________________________________________________ 1.15E-07 __________________________________________________________________ 2.10E-05 __________________________________________________________________ mmu05416 __________________________________________________________________ Viral myocarditis __________________________________________________________________ 1.48E-06 __________________________________________________________________ 0.000135509 __________________________________________________________________ mmu05150 __________________________________________________________________ Staphylococcus aureus infection __________________________________________________________________ 4.24E-06 __________________________________________________________________ 0.000208759 __________________________________________________________________ mmu05168 __________________________________________________________________ Herpes simplex infection __________________________________________________________________ 5.03E-06 __________________________________________________________________ 0.000208759 __________________________________________________________________ mmu05166 __________________________________________________________________ HTLV-I infection __________________________________________________________________ 6.82E-06 __________________________________________________________________ 0.000208759 __________________________________________________________________ mmu05330 __________________________________________________________________ Allograft rejection __________________________________________________________________ 6.84E-06 __________________________________________________________________ 0.000208759 __________________________________________________________________ mmu05332 __________________________________________________________________ Graft-versus-host disease __________________________________________________________________ 9.27E-06 __________________________________________________________________ 0.000242254 __________________________________________________________________ mmu04940 __________________________________________________________________ Type I diabetes mellitus __________________________________________________________________ 1.88E-05 __________________________________________________________________ 0.000430529 __________________________________________________________________ mmu05203 __________________________________________________________________ Viral carcinogenesis __________________________________________________________________ 2.46E-05 __________________________________________________________________ 0.000500679 __________________________________________________________________ mmu05320 __________________________________________________________________ Autoimmune thyroid disease __________________________________________________________________ 5.16E-05 __________________________________________________________________ 0.000944675 __________________________________________________________________ mmu04610 __________________________________________________________________ Complement and coagulation cascades __________________________________________________________________ 0.0001248 __________________________________________________________________ 0.002076447 __________________________________________________________________ mmu04612 __________________________________________________________________ Antigen processing and presentation __________________________________________________________________ 0.0001865 __________________________________________________________________ 0.002621852 __________________________________________________________________ mmu05169 __________________________________________________________________ Epstein-Barr virus infection __________________________________________________________________ 0.0001981 __________________________________________________________________ 0.002621852 __________________________________________________________________ mmu04514 __________________________________________________________________ Cell adhesion molecules (CAMs) __________________________________________________________________ 0.0002006 __________________________________________________________________ 0.002621852 __________________________________________________________________ mmu05162 __________________________________________________________________ Measles __________________________________________________________________ 0.0004025 __________________________________________________________________ 0.004910452 __________________________________________________________________ mmu04380 __________________________________________________________________ Osteoclast differentiation __________________________________________________________________ 0.0005894 __________________________________________________________________ 0.006741672 __________________________________________________________________ mmu05142 __________________________________________________________________ Chagas disease (American trypanosomiasis) __________________________________________________________________ 0.0011778 __________________________________________________________________ 0.012678867 __________________________________________________________________ mmu05134 __________________________________________________________________ Legionellosis __________________________________________________________________ 0.0012619 __________________________________________________________________ 0.012828809 __________________________________________________________________ mmu05133 __________________________________________________________________ Pertussis __________________________________________________________________ 0.0014304 __________________________________________________________________ 0.013776561 __________________________________________________________________ mmu04115 __________________________________________________________________ p53 signaling pathway __________________________________________________________________ 0.003126 __________________________________________________________________ 0.028603336 __________________________________________________________________ mmu04668 __________________________________________________________________ TNF signaling pathway __________________________________________________________________ 0.0052882 __________________________________________________________________ 0.04608275 __________________________________________________________________ Downregulated pathways __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ __________________________________________________________________ mmu04611 __________________________________________________________________ Platelet activation __________________________________________________________________ 0.0012801 __________________________________________________________________ 0.032394895 __________________________________________________________________ mmu00910 __________________________________________________________________ Nitrogen metabolism __________________________________________________________________ 0.00197 __________________________________________________________________ 0.032394895 __________________________________________________________________ mmu04971 __________________________________________________________________ Gastric acid secretion __________________________________________________________________ 0.0024731 __________________________________________________________________ 0.032394895 __________________________________________________________________ mmu05100 __________________________________________________________________ Bacterial invasion of epithelial cells __________________________________________________________________ 0.002757 __________________________________________________________________ 0.032394895 __________________________________________________________________ mmu04974 __________________________________________________________________ Protein digestion and absorption __________________________________________________________________ 0.0040898 __________________________________________________________________ 0.033425935 __________________________________________________________________ mmu00630 Glyoxylate and dicarboxylate metabolism 0.0042671 0.033425935 [95]Open in a new tab Validation of the DEGs with real-time PCR To validate the DEGs from the RNA-seq findings, we prepared new mouse retinas in each group and performed real-time PCR. A total of 11 genes (log[2] fold change >2 and adjusted p<0.05) from the significantly enriched GO terms and KEGG pathways were chosen to be validated. The chosen genes are shown in [96]Table 4. The chosen genes that met the criteria were all upregulated in the LPS-treated group according to the RNA-seq results. The validated genes were chosen from the significantly enriched GO terms, including immune response and immune system process. These genes were also involved in 21 significantly upregulated enriched KEGG pathways, including phagosome, viral myocarditis, staphylococcus aureus infection, herpes simplex infection, HTLV-I infection, allograft rejection, graft-versus-host disease, type I diabetes mellitus, viral carcinogenesis, autoimmune thyroid disease, complement and coagulation cascades, antigen processing and presentation, Epstein-Barr virus infection, cell adhesion molecules (CAMs), measles, osteoclast differentiation, Chagas disease (American trypanosomiasis), legionellosis, pertussis, the p53 signaling pathway, and the TNF signaling pathway. No DEGs met the criteria of a log[2] fold-change greater than 2 and an adjusted p value of less than 0.05 in the significantly downregulated pathways. The 11 chosen genes can be divided into four groups: inflammatory response–related genes (Ccl5, JAK3, and Sell), complement system–related genes (C1ra and C3), fibrinolytic system–related genes (F13a1 and A2m), and cell stress–related genes (Cdkn1a, Gadd45b, Egr1, and Atf3). We found that the results obtained with real-time PCR were in consistent with those obtained with RNA-seq. All the validated genes were significantly upregulated in the LPS group compared with the control group (*p<0.05, **p<0.01,***p<0.001; [97]Figure 7). Table 4. Genes validated by real time PCR. Ensgene Symbol log2 Fold Change Adjusted p value Description ENSMUSG00000035042 __________________________________________________________________ Ccl5 __________________________________________________________________ 3.780941071 __________________________________________________________________ 9.15E-30 __________________________________________________________________ chemokine (C-C motif) ligand 5 __________________________________________________________________ ENSMUSG00000030111 __________________________________________________________________ A2m __________________________________________________________________ 2.741365525 __________________________________________________________________ 2.14E-17 __________________________________________________________________ alpha-2-macroglobulin __________________________________________________________________ ENSMUSG00000031805 __________________________________________________________________ Jak3 __________________________________________________________________ 2.525753144 __________________________________________________________________ 1.26E-11 __________________________________________________________________ Janus kinase 3 __________________________________________________________________ ENSMUSG00000038418 __________________________________________________________________ Egr1 __________________________________________________________________ 2.235443805 __________________________________________________________________ 3.43E-11 __________________________________________________________________ early growth response 1 __________________________________________________________________ ENSMUSG00000026628 __________________________________________________________________ Atf3 __________________________________________________________________ 2.122053519 __________________________________________________________________ 5.52E-08 __________________________________________________________________ activating transcription factor 3 __________________________________________________________________ ENSMUSG00000023067 __________________________________________________________________ Cdkn1a __________________________________________________________________ 2.19450178 __________________________________________________________________ 2.16E-07 __________________________________________________________________ cyclin-dependent kinase inhibitor 1A (P21) __________________________________________________________________ ENSMUSG00000015312 __________________________________________________________________ Gadd45b __________________________________________________________________ 2.179217262 __________________________________________________________________ 3.03E-07 __________________________________________________________________ growth arrest and DNA-damage-inducible 45 beta __________________________________________________________________ ENSMUSG00000055172 __________________________________________________________________ C1ra __________________________________________________________________ 2.171879166 __________________________________________________________________ 3.30E-07 __________________________________________________________________ complement component 1, r subcomponent A __________________________________________________________________ ENSMUSG00000026581 __________________________________________________________________ Sell __________________________________________________________________ 2.251491474 __________________________________________________________________ 1.72E-06 __________________________________________________________________ selectin, lymphocyte __________________________________________________________________ ENSMUSG00000024164 __________________________________________________________________ C3 __________________________________________________________________ 2.15560777 __________________________________________________________________ 2.24E-06 __________________________________________________________________ complement component 3 __________________________________________________________________ ENSMUSG00000039109 F13a1 2.01578419 5.52E-05 coagulation factor XIII, A1 subunit [98]Open in a new tab Figure 7. [99]Figure 7 [100]Open in a new tab Validation of DEGs using real-time PCR. The validated genes were chosen from the significantly enriched Gene Ontology (GO) terms and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway. The mRNA expression of inflammatory response–related genes (Ccl5, JAK3, and Sell; A), complement system–related genes (C1ra and C3; B), fibrinolytic system–related genes (F13a1 and A2m; C), and cell stress–related genes (Cdkn1a, Gadd45b, Egr1, and Atf3; D) were remarkably higher in the lipopolysaccharide (LPS) group compared with the control group (*p<0.05, **p<0.01, ***p<0.001, n = 6). Discussion In the present study, we identified a series of differentially expressed genes occurring in the retinal transcriptome of LPS-induced uveitis mice. We also attempted to identify relevant molecular pathways, with the goal of discovering potential therapeutic targets of uveitis. The results indicated that in addition to the inflammatory response, abnormal changes in complement system–, fibrinolytic system–, and cell stress–related genes may be involved in the pathogenesis of EIU. Treatment targeting these genes may assist in discovering novel therapeutic targets and finding candidate drugs. The top five genes with the greatest fold changes include Ccl5, A2m, Jak3, Egr1, and Atf3. Ccl5 is a key inflammatory mediator of uveitis in animal models [[101]24] and patients [[102]25]. This mediator attracts mononuclear cells to the inflammatory site [[103]26], allowing escalation of the inflammatory response in the ocular vasculature [[104]27]. α2-macroglobulin (A2m) is a specific cytokine carrier that binds inflammatory cytokines implicated in LPS-induced inflammation [[105]28] and has great potential to regulate cytokine homeostasis [[106]29]. Janus kinase 3 (Jak3) plays an important role in the activation of cytokines during immune and inflammation responses [[107]30]. Early growth response 1 (Egr1) is an important transcription factor for regulating the genes involved in inflammation, immunity, and coagulation [[108]31]. Activating transcription factor 3 (Atf3) is a counter-regulatory immune transcription factor induced by Toll-like receptor (TLR) signaling [[109]32]. Existing evidence showed these genes would be increased in response to LPS stimuli [[110]24,[111]33-[112]36]. In parallel with previous studies, the present RNA-seq results showed the mRNA expression of these genes was significantly increased in response to LPS. In general, the GO and KEGG pathway enrichment analysis of the RNA-seq data revealed that four types of genes were significantly upregulated, including the inflammatory response–related genes (Ccl5, JAK3, and Sell), complement system–related genes (C1ra and C3), fibrinolytic system–related genes (F13a1 and A2m), and cell stress–related genes (Cdkn1a, Gadd45b, Egr1, and Atf3). Consistent with the well-accepted notion, we showed a remarkable increase in the inflammatory response–related genes: (Ccl5, JAK3, and Sell). Notably, we also found some complement system– and fibrinolytic system–related genes were increased in response to LPS, indicating that not only inflammatory response occurs during uveitis, but also the complement and fibrinolytic system may be involved in the pathogenesis of EIU. In addition, previous studies showed that a drug targeted to the fibrinolytic system plays a role in the regulation of coagulation and inflammation and may also regulate the complement systems in infectious gastroduodenal disease and sepsis-induced liver injury in mice [[113]37,[114]38], indicating that the inflammation, complement, and fibrinolytic systems may be closely related to and affect each other in the infectious inflammatory responses. The present data suggested that abnormal changes in the inflammatory response–, complement system–, and fibrinolytic system–related genes happened concurrently during EIU and may be associated with the pathogenesis of LPS-induced uveitis. Additionally, candidate drugs that target the genes of these systems may provide better therapeutic strategies for treating patients who do not respond well to conventional drugs. However, the functional evidence of these DEGs needs to be proved in future work. This study has several limitations. First, we analyzed transcriptional changes in the retinas at the 24th hour after LPS injection. The retinal RNA profile presented at this time point could likely not only be reflective of retina cells but also show what is occurring within the infiltrated cells. Samples collected at a single time point may not be enough to reveal the transcriptional changes involved in the development of EIU. Thus, dynamic transcriptome studies of EIU at different time points would provide more comprehensive knowledge. Moreover, the marked influx of blood cells within the posterior segment would also be included in this transcriptome profile. The transcripts derived from these cells cannot be excluded in this study. Therefore, further studies of retina samples treated with perfusion should exclude the influence of these infiltrated cells. Taken together, the retinal transcriptomic analysis, which relied on RNA-seq, was an effective and comprehensive method for acquiring a global view of gene expression changes in the retina after the exposure to LPS. This study may provide new insights into the molecular mechanisms underlying LPS-induced ocular inflammation and may provide useful information toward discovering novel therapeutic targets for ocular inflammatory diseases. Acknowledgments