|Year : 2020 | Volume
| Issue : 1 | Page : 39-46
Analysis of differentially expressed genes in bacterial and fungal keratitis
Rui Tian, He Zou, Lufei Wang, Lu Liu, Meijiao Song, Hui Zhang
Department of Ophthalmology, The Second Hospital of Jilin University, Changchun, Jilin Province, China
|Date of Submission||09-Jan-2019|
|Date of Acceptance||12-Jul-2019|
|Date of Web Publication||19-Dec-2019|
Dr. Hui Zhang
Department of Ophthalmology, The Second Hospital of Jilin University, No. 218, Ziqiang Street, Nan Guan District, Changchun, 130000, Jilin Province
Source of Support: None, Conflict of Interest: None
Purpose: This study was aimed at identifying differentially expressed genes (DEGs) in bacterial and fungal keratitis. The candidate genes can be selected and quantified to distinguish between causative agents of infectious keratitis to improve therapeutic outcomes. Methods: The expression profile of bacterial or fungal infection, and normal corneal tissues were downloaded from the Gene Expression Omnibus. The limma package in R was used to screen DEGs in bacterial and fungal keratitis. The Co-Express tool was used to calculate correlation coefficients of co-expressed genes. The "Advanced network merge" function of Cytoscape tool was applied to obtain a fusional co-expression network based on bacterial and fungal keratitis DEGs. Finally, functional enrichment analysis by DAVID software and KEGG analysis by KOBAS of DEGs in fusion network were performed. Results: In total, 451 DEGs in bacterial keratitis and 353 DEGs in fungal keratitis were screened, among which 148 DEGs were found only in bacterial keratitis and 50 DEGs only in fungal keratitis. Besides, 117 co-expressed gene pairs were identified among bacterial keratitis DEGs and 87 pairs among fungal keratitis DEGs. In total, nine biological pathways and seven KEGG pathways were screened by analyzing DEGs in the fusional co-expression network. Conclusion: TLR4 is the representative DEG specific to bacterial keratitis, and SOD2 is the representative DEG specific to fungal keratitis, both of which are promising candidate genes to distinguish between bacterial and fungal keratitis.
Keywords: Bacterial keratitis, co-expression network, differentially expressed genes (DEGs), fungal keratitis
|How to cite this article:|
Tian R, Zou H, Wang L, Liu L, Song M, Zhang H. Analysis of differentially expressed genes in bacterial and fungal keratitis. Indian J Ophthalmol 2020;68:39-46
|How to cite this URL:|
Tian R, Zou H, Wang L, Liu L, Song M, Zhang H. Analysis of differentially expressed genes in bacterial and fungal keratitis. Indian J Ophthalmol [serial online] 2020 [cited 2020 Aug 7];68:39-46. Available from: http://www.ijo.in/text.asp?2020/68/1/39/273262
Diseases of the cornea are a major cause of blindness worldwide.  The etiology of corneal blindness encompasses a wide variety of inflammatory and infectious eye diseases that ultimately cause functional blindness. [1,2] Keratitis is a type of corneal inflammation resulting in vision loss. It typically arises due to noninfectious causes such as eye trauma but can manifest as a result of microbial infection by pathogens such as fungi, bacteria, viruses or amebae.  Until now, infectious keratitis remains one of the main causes of corneal blindness and poses a diagnostic dilemma due to its varied presentation and visual morbidity. ,
Currently, bacterial keratitis and fungal keratitis are the most common corneal infectious diseases posing a risk to patient vision. , The major causative pathogens for bacterial keratitis are Staphylococcus aureus and Pseudomonas aeruginosa.  Bacterial keratitis frequently leads to severe visual impairment from corneal ulceration, perforation, and scarring.  Following an infection, topical antimicrobial therapy is crucial for managing symptoms. [9,10] Risk factors of fungal keratitis include ocular trauma, topical steroid use, ocular surface disease, and contact lens use.  Aspergillus spp., Fusarium spp., Candida spp., are the major causative pathogens of fungal keratitis among many.  Fungal keratitis commonly leads to poor visual acuity,  and is typically managed by polyenes and azoles. 
Each case of infectious keratitis must be confirmed by evaluating corneal infiltrate cultures. , Clinically, corneal ulcers are often treated empirically without the use of microbiological analysis due to urgent requests for treatment to achieve optimal therapeutic outcomes.  In order to rely on empirical treatment, the clinician must distinguish between infectious agents based on clinical history, symptoms and characteristics. This method remains highly subjective and risky as incorrect identification of the pathogen facilitates further development of the corneal infection, ultimately leading to a worsened therapeutic outcome. Hence, it is necessary to identify and develop novel approaches to quickly recognize or identify bacterial versus fungal keratitis.
For the specific treatment of infectious keratitis, it is important to reveal the functional and molecular aspects of the disease to develop a possible treatment strategy. In recent years, the analysis of differentially expressed genes (DEGs) in disease has attracted a lot of attention, and may be a promising approach to develop more efficient treatments for keratitis. In this study, we aimed to screen the DEGs in bacterial and fungal keratitis by comparing total gene expression levels in infected versus healthy corneal tissues. This strategy will allow for the identification of candidate genes that can be therapeutically targeted to treat keratitis originating from different infectious agents.
| Methods|| |
Affymetrix microarray data
The transcription profile of GSE58291 was downloaded from the Gene Expression Omnibus (GEO; http://www.ncbi.nlm.nih.gov/geo/), a public functional genomics data repository that archives and freely distributes high-throughput molecular abundance data at the National Center for Biotechnology Information. In total, 30 corneal tissue samples were acquired. Among them, three samples showed empty expression profile data, numbered as GSM1406007 (fungal infection), GSM1406009 (bacterial infection), and GSM1406015 (normal control). Hence, 27 tissue samples (12 normal corneas, 7 bacteria-infected corneas, and 8 fungi-infected corneas) were reserved for bioinformatic analysis. Detailed information of the 27 samples is listed in [Table 1]. The causative organisms for bacterially infected corneas included Streptococcus pneumonia (n = 6) and Pseudomonas. aeruginosa (n = 1). The causative organisms for fungal keratitis were Fusarium sp. (n = 5), Aspergillus sp. (n = 2, A. flavus and A. terreus) and Lasiodiplodia sp. (n = 1). Platform information was GPL10558 Illumina Human HT-12 V4.0 expression beadchip. Platform annotation information of the chip expression profiles was also downloaded.
Data preparation and differential gene expression analysis
The raw expression profile data in text format were mapped to the corresponding gene names using the GPL10558 Illumina HumanHT-12 V4.0 expression beadchip platform. The average expression value was calculated as the single expression value of this gene, when multiple probes matched to the same gene. Then, logarithm to the base 2 (log2) of expression values was calculated to acquire approximately normally distributed gene expression data, which were continuously subjected to median normalization. [17,18] According to sample infection types, comparisons were performed in the bacterial infection versus normal control group, as well as in the fungal infection versus normal control group. Here, samples of normal cornea tissues were classified as the normal control group. The limma  package in R was used to screen DEGs by analyzing the gene expression data of corneal tissues from the above three groups. The Bonferroni's method  in multi-test package was applied to adjust raw P values for false discovery rate (FDR).  FDR <0.05 and the absolute value of log2FC >1 were used as cut-off criteria.
Comparisons of gene expression profiles
Gene expression profiles are species-specific, suggesting that gene expression is significantly altered in diseased tissues.  According to the expression profile of screened DEGs in bacteria versus fungal and normal samples, we extracted the expression value of DEGs in each sample from the downloaded expression value files. Then, the pheatmap package in R was used to generate expression values by biclustering , based on Euclidean distance.  The results are shown as a heatmap.
Calculations of co-expression correlation coefficient among DEGs
Although there are approximately 25,000 genes in the human genome, only a fraction of these genes are expressed simultaneously in a single cell or specific tissues during a specific developmental stage.  There are many methods to identify whether co-expression exists between two genes, iamong which the most common method is to use Pearson's correlation coefficient.  To obtain DEGs with correlations, the CoExpress tool  (http://www.bioinformatics.lu/CoExpress/) was used to calculate correlation coefficients among co-expressed DEGs in the bacterial versus normal group or the fungal versus normal group. Finally, gene pairs with the absolute value of correlation coefficients >0.9 were retained.
Difference between bacterial versus normal DEGs and fungus versus normal DEGs
Through comparison of the gene expression profile between bacterial and normal groups, we acquired the screened DEGs, which are referred to as DEGs1. Similarly, through comparison between fungal and normal groups, we acquired more screened DEGs, which are referred to as DEGs2. To compare the differences between DEGs1 and DEGs2, a Venn diagram was used. 
The fusion of co-expression network
Based on DEGs1 and DEGs2, we acquired two corresponding co-expression networks by gene co-expression network analysis. Then, the DEGs1-based co-expression network was merged with the DEGs2-based co-expression network using the Cytoscape tool "Advanced network merge" to obtain a unique network. 
Functional enrichment analysis of DEGs in fusion co-expression network
Currently, there are multiple tools for gene function enrichment analysis, among which DAVID has been widely used.  Using DAVID software, the biological pathways significantly enriched by DEGs in the fusion co-expression network were identified. A P value less than 0.05 was used as a screening threshold.
Pathway analysis of DEGs in fusion co-expression network
Continuously, the pathway annotations and enrichment analysis were completed using KOBAS  based on algorithm of accumulative hypergeometric distribution. A P value less than 0.05 was used as a screening threshold.
| Results|| |
Data preprocessing and DEGs' identification
To remove system errors under sequencing, the data were preprocessed. Through data preparation described in the methods section, we obtained the normalized gene expression data. After preprocessing, the medians of expression values in all the samples were relatively linear, suggesting that the expression data were well-normalized [Figure 1]. A total of 451 DEGs were obtained from the comparison between the bacterial and normal groups and 353 DEGs from the comparison between the fungal and normal groups.
|Figure 1: The boxplot of expression profiling at prestandardization (a) and poststandardization (b). Boxes with white, light gray, and dark gray represent normal cornea, bacterial-infected cornea, and fungus-infected cornea samples, respectively. Y-axis represents gene expression value|
Click here to view
We observed that the screened DEGs could significantly distinguish bacteria/fungus-infected from normal corneal samples [Figure 2]. 2a and b. These results indicated that significant sample differences existed among screened DEGs between bacteria-infected and normal groups, as well as between fungus-infected and normal groups.
|Figure 2: The heap map of screened DEGs. (a) Heap map of DEGs screened between bacterial-infected and normal groups. (b) Heap map of DEGs screened between fungus-infected and normal groups. Red color represents high expression, while blue color represents low expression. Color changes from blue to red indicate the corresponding expression value change from lower to higher|
Click here to view
Difference between DEGs between bacterial versus normal and fungus versus normal groups
To compare the difference in DEGs in bacterial versus normal and that in fungal versus normal, a Venn diagram was constructed. We observed that the number of overlapped DEGs was 303, which accounted for 67.18% (303/451) among bacterial versus normal DEGs and 85.84% (303/353) among fungal versus normal DEGs, respectively [Figure 3]. There were 148 DEGs specific to bacterial keratitis, such as CD34 (CD34 molecule, P = 1.40E-09, low expression), HK2 (hexokinase 2, P = 6.79E-07, overexpression), and TLR4 (toll-like receptor 4, P = 2.35E-09, overexpression) and 50 specific DEGs in fungal keratitis, such as ADH7 (alcohol dehydrogenase 7, P = 3.85E-04, low expression), ASGR1 (asialoglycoprotein receptor 1, P = 5.49E-08, overexpression), and SOD2 (superoxide dismutase 2, P = 7.89E-05, overexpression) [Table 2]. These results suggest that there were a large number of DEGs identified both in bacteria-infected and fungus-infected corneas, exhibiting tremendous similarities in the above two keratopathies. In addition, the same DEG showed homodromous up- or down-regulation in bacteria and fungus-infected cornea samples, also exhibiting complete uniformity in the above two keratopathies.
|Figure 3: Venn diagram of DEG sets between bacterial vs. normal (left) and fungus vs. normal groups (right)|
Click here to view
Calculations of co-expression correlation coefficient among DEGs
The number of co-expressed gene pairs was 117 pairs among bacterial versus normal DEGs and 87 pairs among fungal versus normal DEGs, respectively. The co-expression networks were visualized using Cytoscape tool to obtain the corresponding network graphs.
The fusion of co-expression networks
After fusion of bacterial versus normal and fungal versus normal coexpression networks, a novel fusion coexpression network was generated [Figure 4]. This fusion co-expression network included 79 DEG nodes and 190 connecting edges. Among these 79 DEG nodes, 19 were unique to bacterial versus normal DEGs, 5 were unique to fungal versus normal DEGs, and 55 were present in both comparisons.
|Figure 4: The fusion coexpression network merged from bacterial vs. normal and fungus vs. normal coexpression networks. Dark gray and light gray represent bacterial vs. normal DEGs and fungus vs. normal DEGs, respectively. Triangle and inverted triangle represent up- and down-regulation DEGs, respectively. White rhombus represents DEGs identified both in bacterial vs. normal DEGs and fungus vs. normal DEGs|
Click here to view
Function enrichment analysis of DEGs in fusion co-expression networks
Through analysis of DEGs in fusion co-expression networks using DAVID, we searched nine biological pathways in total that were significantly differentially regulated [Supplemental Table 1] [Additional file 1]. Among these nine biological pathways, the immune response was the most significant pathway. Notably, the other eight biological pathways were mainly associated with the immune system.
Pathway analysis of DEGs in fusion co-expression network
In total, DEGs in fusion co-expression network were involved seven KEGG pathways [Supplemental Table 2] [Additional file 2], among which antigen processing and presentation (hsa04612) was the most striking. Specifically, there were five DEGs identified both between the bacterial versus normal groups and the fungal and normal groups, including HLA-DRB3, IFI30, HLA-DPA1, HLA-DMB, and HLA-DMA, involved in antigen processing and presentation pathways. Among them, HLA-DRB3, HLA-DPA1, HLA-DMB, and HLA-DMA were also involved in the immune response.
| Discussion|| |
In this study, we primarily screened the DEGs in bacterial keratitis and fungal keratitis through analyzing the gene expression profiles of corneal tissues. In total, there were 451 DEGs identified from bacterial keratitis versus normal corneal tissues and 353 DEGs identified from fungal keratitis versus normal corneal tissues. The number of overlapping DEGs between bacterial keratitis and fungal keratitis was 303, which accounted for a larger proportion in corresponding keratitis. In addition, through co-expression network analysis, 117 co-expressed gene pairs were identified in bacterial keratitis DEGs and 87 pairs in fungal keratitis DEGs. After constructing the fusional co-expression network based on bacterial and fungal keratitis co-expression DEGs, nine biological pathways by function enrichment analysis and seven KEGG pathways by KEGG analysis were identified as significant.
Toll-like receptor 4 (TLR4) is a crucial pattern recognition molecule that participates in the innate immune response to lipopolysaccharide, a vital component of Gram-negative bacteria.  It is reported that TLR4 mRNA levels were significantly upregulated in bacterial (P. aeruginosa) infected mouse cornea tissue.  In accordance with this study, our data confirmed that TLR4 levels were significantly increased by approximately 5-fold in human cornea tissues with bacterial keratitis compared to those in cornea tissues from healthy donors.  The deficiency of TLR4 in mouse could result in increased polymorphonuclear neutrophil infiltration and proinflammatory cytokine production, as well as decreased β-defensin-2 and inducible nitric oxide synthase production in mouse with P. aeruginosa infection of the cornea.  Yan et al. reported that TLR4 found in corneal macrophages could regulate P. aeruginosa keratitis by signaling through myeloid differentiation factor 88 (MyD88)-dependent and -independent pathways.  In addition, TLR4 was also reported to regulate fungal keratitis such as fusarium keratitis  and A. fumigatus keratitis,  although from our analysis, there were no significant differences in TLR4 expression between corneal tissues infected with a fungal pathogen versus normal tissues. This suggests that TLR4 is a candidate target gene to distinguish bacterial keratitis from fungal keratitis. A promising strategy for the diagnosis of infectious keratitis may be developed based on TLR4 expression.
Among the 50 non-overlapping DEGs in fungal keratitis, SOD2 levels were found to be significantly increased by about 9-fold in human corneal tissues with fungal keratitis compared to those in normal human corneal tissues. It has been previously reported that SOD2 expression is significantly increased by about 2-fold in mouse corneas with fungal (Candida albicans) keratitis compared to that in healthy mouse corneas.  We speculated that the differences in SOD2 fold-change between mouse and human could be due to species-specific differences and fungal species. Moreover, SOD2 was a pivotal DEG node in the fusional co-expression network, which was derived from fungus keratitis DEGs, and co-expression is associated with MYOC, KRT6B, and CSF1R. Meanwhile, SOD2 was identified to participate in responses to wounding and oxidation reduction pathways. Although there is currently no report describing a role for SOD2 in keratitis, SOD2 still remains a potential candidate gene to distinguish between bacterial and fungal keratitis due to its specific association with fungal keratitis.
Furthermore, through functional analysis of DEGs in fusion co-expression networks, we identified nine biological pathways such as pro-inflammatory and anti-inflammatory responses, antigen processing and presentation,  and wounding responses.  The majority of the identified pathways are associated with the immune system. Through KEGG analysis of DEGs in fusion co-expression network, we identified seven KEGG pathways, of which antigen processing and presentation  intestinal immune network for IgA production, autoimmune thyroid disease, and viral myocarditis  are more associated with the immune system. This suggests that perturbations in the immune system induced by pathogen exposure in the cornea leads to the malignant advance of infectious keratitis. Based on our findings, we speculate that strategies aimed at controlling inflammation are a compensatory therapy to alleviate the pain experienced by patients with keratitis excepting to eliminate pathogens, which requires further investigation of the identified immune-related DEGs in infectious keratitis.
In this study we identified novel DEGs associated with bacterial or fungal keratitis; however, our study does have limitations. Firstly, we were limited in the clinical materials including verification of our samples. Secondly, due to the limitations of obtaining human tissue samples, the sample size remains small. Future studies will be required to corroborate our findings using larger sample sizes. Finally, the expression levels of the candidate genes may be affected by the nature of the pathogen, the stage of the disease, or the genetic background of the host. Thus, studies with large sample sizes are warranted to validate our findings in the near future.
| Conclusion|| |
In summary, our work screened 451 DEGs in corneas with bacterial keratitis and 353 DEGs in corneas with fungal keratitis, in which 148 DEGs were found to be specific to bacterial keratitis and 50 DEGs specific in fungal keratitis. TLR4 was an upregulated gene specific in bacterial keratitis and SOD2 was an upregulated gene specific to fungal keratitis. Both genes are promising candidate targets to distinguish bacterial and fungal keratitis.
Financial support and sponsorship
Conflicts of interest
There are no conflicts of interest.
| References|| |
Whitcher JP, Srinivasan M, Upadhyay MP. Corneal blindness: A global perspective. Bull World Health Organ 2001;79:214-21.
Cao J, Yang Y, Yang W, Wu R, Xiao X, Yuan J, et al
. Prevalence of infectious keratitis in Central China. BMC Ophthalmol 2014;14:43.
Estimated Burden of Keratitis - United States, 2010. MMWR Rep 2014;63:1027-30.
Solanki S, Rathi M, Khanduja S, Dhull CS, Sachdeva S, Phogat J. Recent trends: Medical management of infectious keratitis. Oman J Ophthalmol 2015;8:83-5. [Full text]
Abdullah AM, Nadia AK, Archana T, Wali UK. Bacterial keratitis: perspective on epidemiology, clinico-pathogenesis, diagnosis and treatment. Sultan Qaboos Univ Med J 2009;9:184-95.
Yuan X, Mitchell BM, Wilhelmus KR. Gene profiling and signaling pathways of Candida albicans keratitis. Mol Vis 2008;14:1792-8.
Henri S, Kaye SB, Timothy N, Amanda H, Stephen T, Parry CM. An in vitro
investigation of synergy or antagonism between antimicrobial combinations against isolates from bacterial keratitis. Invest Ophthalmol Vis Sci 2010;51:4151-5.
Palioura S, Henry CR, Amescua G, Alfonso EC. Role of steroids in the treatment of bacterial keratitis. Clin Ophthalmol 2016;10:179.
Bower KS, Kowalski RP, Gordon YJ. Fluoroquinolones in the treatment of bacterial keratitis. Am J Ophthalmol 1996;121:712.
Moshirfar M, Mirzaian G, Feiz V, Kang PC. Fourth-generation fluoroquinolone-resistant bacterial keratitis after refractive surgery. J Cataract Refract Surg 2006;32:515-8.
Jurkunas U, Behlau I, Colby K. Fungal keratitis: Changing pathogens and risk factors. Cornea 2009;28:638.
Kalyana CS, Jayasudha R, Ranjith K, Dutta A, Pinna NK, Mande SS, et al
. Alterations in the gut bacterial microbiome in fungal Keratitis patients. PLoS One. 2018;13:e0199640.
Venkatesh P, Krishnan T, Mascarenhas J, Srinivasan M, Oldenburg CE, Toutain-Kidd CM, et al
. Predictors of outcome in fungal keratitis. Eye 2012;26:1226-31.
Chin Fen N, Mark D, Sharon CA, Kay S, Kong DC. Clinical utility of caspofungin eye drops in fungal keratitis. Int J Antimicrob Agents 2014;44:96-104.
Ibrahim YW, Boase DL, Cree IA. Epidemiological characteristics, predisposing factors and microbiological profiles of infectious corneal ulcers: The Portsmouth corneal ulcer study. Br J Ophthalmol 2009;93:1319.
Cyril D, Porco TC, Lietman TM, Venkatesh P, Lalitha P, Mano Ranjan D, et al
. The clinical differentiation of bacterial and fungal keratitis: A photographic survey. Invest Ophthalmol Vis Sci 2012;53:1787.
Troyanskaya O, Cantor M, Sherlock G, Brown P, Hastie T, Tibshirani R, et al
. Missing value estimation methods for DNA microarrays. Bioinformatics 2001;17:520-5.
Fujita A, Sato JR, Rodrigues LDO, Ferreira CE, Sogayar CM. Evaluating different methods of microarray data normalization. BMC Bioinform 2006;7:469.
Smyth GK. Limma: Linear models for microarray data. In: Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, R Irizarry WH. New York: Springer; 2005. p. 397-420.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: A practical and powerful approach to multiple testing. J Roy Stat Soc Ser B 1995;57:289-300.
Benjamini Y. Discovering the false discovery rate. J Roy Stat Soc Ser B 2010;72:405-16.
Ester M, Kriegel H-P, Sander J, Xu X. A density-based algorithm for discovering clusters in large spatial databases with noise. In: Simoudis E, Han J, Fayyad UM, editors. Proceedings of the Second International Conference on Knowledge Discovery and Data Mining (KDD-96). AAAI Press; 1996. p. 226-231. ISBN 1-57735-004-9.
Szekely GJ, Rizzo ML. Hierarchical clustering via joint between-within distances: Extending Ward's minimum variance method. J Classification 2005;22:151-83.
Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Section 16.4. Hierarchical Clustering by Phylogenetic Trees; Numerical Recipes: The Art of Scientific Computing, 2007. p. 868-82.
Deza EDMM. Encyclopedia of Distances. Springer; 2009. p. 94.
Ross D, Scherf U, Mb, Perou C, Rees C, Spellman P, Iyer V, et al
. Systematic variation in gene expression patterns in human cancer cell lines (see comments). Nat Genet 2000;24:227.
Takeshi O, Kengo K, Kenta N, Masayuki S, Shinpei H, Motoshi S, et al
. ATTED-II: A database of co-expressed genes and cis elements for identifying co-regulated gene groups in Arabidopsis. Nucl Acids Res 2007;35(Database issue):D863-9.
Nazarov PV, Muller A, Khutko V, Vallar L. Co-expression analysis of large microarray data sets using CoExpress software tool. Conference: WCSB 2010.
Ruskey FW, M. Venn diagram survey. Elect J Combinatorics 2005. p. 5.
Yang L, Lee YF, Ng MK. SNP and gene networks construction and analysis from classification of copy number variations data. BMC Bioinform 2011;12:S4-S.
Da Wei H, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucl Acids Res 2009;37:1.
Jianmin W, Xizeng M, Tao C, Jingchu L, Liping W. KOBAS server: A web-based platform for automated annotation and pathway identification. Nucl Acids Res 2006;34(Web Server issue):W720-W4.
Hara Y, Shiraishi A, Ohashi Y. Hypoxia-altered signaling pathways of toll-like receptor 4 (TLR4) in human corneal epithelial cells. Mol Vis 2009;15:2515-20.
Huang X, Du WS, Barrett RP, Hazlett LD. TLR4 is required for host resistance in Pseudomonas aeruginosa keratitis. Invest Ophthalmol Vis Sci 2006;47:4910.
Yan S, Mausita K, Sanhita R, Ramadan RT, Williams SR, Scott H, et al
. TLR4 and TLR5 on corneal macrophages regulate Pseudomonas aeruginosa keratitis by signaling through MyD88-dependent and -independent pathways. J Immunol 2010;185:4272-83.
Tarabishy AB, Bishr A, Yan S, Yoshifumi I, Mukherjee PK, Lass JH, et al
. MyD88 regulation of Fusarium keratitis is dependent on TLR4 and IL-1R1 but not TLR2. J Immunol 2008;181:593-600.
Gao X, Zhao G, Cui L, Jing L, Nan J, Qian W, et al
. LOX-1 and TLR4 affect each other and regulate the generation of ROS in A. fumigatus keratitis. Int Immunopharmacol 2016;40:392-9.
Hua X, Chi W, Su L, Li J, Zhang ZD, Yuan XY. ROS-induced oxidative injury involved in pathogenesis of fungal keratitis via p38 MAPK activation. Sci Rep 2017;7:10421.
Francesca S, Kyla O, David MD, Seipp RP, Jennifer H, Rayshad G, et al
. Epigenetic enhancement of antigen processing and presentation promotes immune recognition of tumors. Cancer Res 2008;68:9601-7.
Silvio E, Mario P, Lattorff HMG. Dynamics of immune system gene expression upon bacterial challenge and wounding in a social insect (Bombus terrestris). PLoS One 2011;6:e18126.
Maisch B, Ristiæ AD, Hufnagel G, Pankuweit S. Pathophysiology of viral myocarditis: The role of humoral immune response. Cardiovasc Pathol 2002;11:112-22.
[Figure 1], [Figure 2], [Figure 3], [Figure 4]
[Table 1], [Table 2]