Keratoconus patients exhibit a distinct ocular surface immune cell and inflammatory profile

Quantitative evaluation of morphological and functional changes in meibomian glands and lipid layer thickness in patients with and without keratoconus

PURPOSE

To quantitatively evaluate the morphological parameters of meibomian glands (MGs) and lipid layer thickness (LLT) in patients with keratoconus (KC).

METHODS

In this prospective, cross-sectional study, 164 eyes of 164 keratoconus patients and 64 eyes of 64 age-matched control subjects were included. An advanced automatic MG analyzer was used to quantitatively measure the morphological and functional parameters of MGs. Morphological and functional parameters of MGs, LLT, and other ocular surface parameters were compared between the control and KC groups.

RESULTS

The mean meibomian gland diameter, length, square, and gland area ratio (GA) were all significantly decreased in the KC group (all P < 0.05), while no significant difference was observed in the gland tortuosity index (TI) and gland signal index (SI) between the KC and control groups (all P > 0.05). There was no significant difference in the number of total and incomplete blinking among patients with different stages of keratoconus (all P > 0.05). The gland diameter, square, and TI were all negatively associated with KC severity (all P < 0.05), while no significant difference was observed among all stages of KC in gland length, GA, and SI (all P > 0.05). Moreover, the LLTs were positively correlated with the gland diameter, square, GA, and TI and negatively correlated with anterior corneal curvature or KC severity (all P < 0.05).

CONCLUSIONS

Atrophic morphological changes in the meibomian glands were closely correlated with the severity of keratoconus. Gland diameter may be a sensitive functional morphology metric of meibomian glands in patients with keratoconus.

Inflammation and keratoconus: A comprehensive bidirectional Mendelian randomization analysis

An increasing body of evidence supports the involvement of inflammation and immune responses in the occurrence and development of keratoconus (KC). However, the causal relationship between inflammatory factors and KC remains unclear. We employed a 2-way Mendelian randomization (MR) approach to investigate the interaction between KC and inflammatory factors. Instrumental variables for 41 circulating inflammatory regulators and 12 matrix metalloproteinases (MMPs) were selected from genome-wide association studies of European ancestry. Summary statistics for KC were obtained from a genome-wide association study comprising 2116 cases and 24,626 controls of European ancestry. The primary analytical method for assessing causality was the inverse-variance weighted method. Two additional MR methods (MR-Egger and weighted median) were employed to complement the inverse-variance weighted results. In addition, several sensitivity analyses were conducted to evaluate heterogeneity, horizontal pleiotropy, and stability. Our findings indicated that genetically predicted higher levels of macrophage inflammatory protein-1β (odds ratio = 1.126, 95% confidence interval: 1.029-1.232, P = .01) and MMP-13 (odds ratio = 1.211, 95% confidence interval: 1.070-1.371, P = .003) were positively associated with an elevated risk of KC. Conversely, genetically predicted KC was associated with increased levels of interferon-gamma, interleukin-4, and MMP-1. Our current study provided suggestive evidence supporting causal associations of macrophage inflammatory protein-1β and MMP-13 with the risk of KC. In addition, KC appeared to affect the expression of interferon-gamma, interleukin-4, and MMP-1.

Investigation of the rhythmic recruitment of tear neutrophils to the ocular surface and their phenotypes

Hundreds of thousands of polymorphonuclear neutrophils (PMNs) are collected from the ocular surface upon waking, while few are harvested during daytime. This study aimed to investigate potential factors contributing to the circadian infiltration of tear PMNs, including changes in IL-8 and C5a in tears, and their phenotypes across different time points in a 24-h cycle. Tear PMNs were collected using a gentle eyewash after 2-h and 7-h of sleep (eye closure, EC) at night, after 2-h EC during the day, and towards the end of the afternoon. Significantly fewer cells were collected after 2-h EC during the day compared to 2-h EC at night. A positive correlation between IL-8 and PMN numbers existed, but not with C5a. Tear PMNs collected after 2-h EC at night were less degranulated and possessed a larger activation potential compared to 7-h EC. Tear PMNs from 7-h EC at night exhibited hyper-segmented nuclei and more NETosis compared to 2 h EC night, indicating an aged and activated phenotype. The diurnal-nocturnal recruitment pattern of tear PMNs may be driven by increased IL-8 in nighttime tears. Higher degranulation and NETs point to the significant activation of tear PMNs on the ocular surface during prolonged eye closure at night.

[Changes in the precorneal tear film at various stages of keratoconus]

The combination of keratoconus (KC) with signs of dry eye disease (DES) has been described in numerous scientific publications. At the same time, there is a relationship between KC stage and an increase in the severity of DES symptoms, however, there is still no common understanding of the severity of xerotic process depending on the clinical course of keratectasia.

PURPOSE

This study assesses the changes in the state of the precorneal tear film in KC relative to the stages of the disease.

MATERIAL AND METHODS

The study included 50 patients (100 eyes) with bilateral non-operated KC from subclinical to stage IV. The following methods were used for a comprehensive assessment of the condition of the tear film: biomicroscopy of the anterior segment of the eye using vital dyes, functional tests (Norn, Schirmer, Jones tests), tearscopy of the lipid layer of the tear film with software processing of the results, OCT-meniscometry, as well as filling out the OSDI questionnaire by patients.

RESULTS

The following significant changes were revealed as the KC stage progressed: an increase in the area of staining of the ocular surface with vital dyes, a decrease in the results of the Norn functional test and OCT-meniscometry, an increase in the areas of smaller thickness of lipids and areas of their complete absence according to tearscopy, as well as an increase in scores of the OSDI questionnaire.

CONCLUSION

A comprehensive assessment of the condition of the tear film in KC indicates the relationship and the strengthening of the signs of DES as keratectasia progresses, and can be considered as justification for the need to prescribe tear replacement and reparative therapy that improves the condition of the ocular surface and stabilizes the precorneal tear film.

Characteristics of Autophagy-Related Genes, Diagnostic Models, and Their Correlation with Immune Infiltration in Keratoconus

Purpose

Keratoconus (KTCN) is one of the most common degenerative keratopathies, significantly affecting vision and even leading to blindness. This study identifies potential biomarkers of KTCN based on the characterization of autophagy-related genes (ARGs) and the construction of a diagnostic model; and explores their relevance to immune infiltrating cells in KTCN.

Methods

Gene Expression Omnibus (GEO) data were downloaded and ARGs were acquired from GeneCards and Molecular Signatures Database (MSigDB). Autophagy-related differential expression genes (ARDEGs) were discovered through the integration of differentially expressed genes (DEGs) with ARGs, while hub genes of KTCN were discovered by protein-protein interaction (PPI) network analysis. The probable biological roles of these hub ARDEGs were examined using functional enrichment analysis, and a KTCN diagnostic model was generated using the least absolute shrinkage and selection operator (LASSO) regression analysis. We also employed the CIBERSORTx and ssGSEA algorithms to identify potential regulatory pathways to compare the abundance of immune cell infiltrates and their association with hub genes. Finally, the hub gene expression levels were confirmed using validation datasets as well as blood samples from KTCN and healthy individuals.

Results

In this study, we identified 12 hub ARDEGs, of which 9 genes were substantially distinct between KTCN patients and normal groups. The LASSO risk score was used to generate the nomogram, and the calibration curve evaluated the model's effective diagnostic performance (C index of 0.961). Patients with KTCN had greater percentages of M2 Macrophages and Gamma delta T cells, according to CIBERSORTx and ssGSEA. The outcomes of the bioinformatics analysis were supported by the DDIT3 and BINP3 expression levels in KTCN patients and healthy controls, according to the qRT-PCR data.

Conclusion

Five biomarkers (CFTR, PLIN2, DDIT3, BAG3, and BNIP3) and diagnostic models offer fresh perspectives on identifying and managing KTCN.

Screening and identification of genes related to ferroptosis in keratoconus

Corneal keratoconus (KC) is a dilated (ectatic) corneal disease characterized by a central thinning of the cornea, which causes protrusion into a conical shape that seriously affects vision. However, due to the complex etiology of keratoconus, its entire mechanism remains unclear and there is no mechanism-directed treatment method. Ferroptosis is a novel programmed cell death mechanism related to lipid peroxidation, stress, and amino acid metabolism, which plays a crucial role in various diseases. This study aimed to explore the relationship between keratoconus and ferroptosis, to provide new insights into the mechanism of keratoconus development, and potential treatment options based on further elucidation of this mechanism. The corresponding mRNA microarray expression matrix data of KC patients were obtained from GEO database (GSE204791). Weighted co-expression network analysis (WGCNA) and support vector machine recursive feature elimination (SVM-RFE) were selected to screen hub genes, which were overlapped with ferroptosis genes (FRGs) from FerrDb. GO and GSEA were performed to analyze differential pathways, ssGSEA was used to determine immune status, and then, feasible drugs were predicted by gene-drug network. Additionally, we predicted the miRNA and IncRNA of hub genes to identify the underlying mechanism of disease so as to predict treatment for the disease. The epithelial transcriptome from keratoconus tissue mRNA microarray data (GSE204791) was extracted for the main analysis, including eight epithelial cells and eight epithelial control cells. The differential genes that were overlapped by WGCAN, SVM-RFE and FRGs were mainly related to oxidative stress, immune regulation, cellular inflammation, and metal ion transport. Through further analysis, aldo-keto reductase family 1 member C3 (AKR1C3) was selected, and negatively correlated with mature CD56 natural killer (NK) cells and macrophages. Then, gene-drug interaction network analysis and miRNA prediction were performed through the website. It was concluded that four immune-related drugs (INDOMETHACIN, DAUNORUBICIN, DOXORUBICIN, DOCETAXEL) and a miRNA (has-miR-184) were screened to predict potential drugs and targets for disease treatment. To our knowledge, this was the first report of KC being associated with ferroptosis and prompted search for differential genes to predict drug targets of gene immunotherapy. Our findings provided insight and a solid basis for the analysis and treatment of KC.

Identifying immune cell infiltration and effective diagnostic biomarkers in Crohn's disease by bioinformatics analysis

Background

Crohn's disease (CD) has an increasing incidence and prevalence worldwide. It is currently believed that both the onset and progression of the disease are closely related to immune system imbalance and the infiltration of immune cells. The aim of this study was to investigate the molecular immune mechanisms associated with CD and its fibrosis through bioinformatics analysis.

Methods

Three datasets from the Gene Expression Omnibus data base (GEO) were downloaded for data analysis and validation. Single sample gene enrichment analysis (ssGSEA) was used to evaluate the infiltration of immune cells in CD samples. Immune cell types with significant differences were identified by Wilcoxon test and Least Absolute Shrinkage and Selection Operator (LASSO) regression analysis. Differentially expressed genes (DEGs) were screened and then subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional correlation analysis, as well as protein-protein interaction (PPI) network analysis. The cytoHubba program and the GSE75214 dataset were used to screen for hub genes and plot Receiver operating characteristic (ROC)curves to screen for possible biomarkers of CD based on diagnostic efficacy. The hub genes of CD were correlated with five significantly different immune cells. In addition, validation was performed by real time quantitative PCR (RT-qPCR) experiments in colonic tissue of CD intestinal fibrosis rats to further identify hub genes that are more related to CD intestinal fibrosis.

Results

The DEGs were analyzed separately by 10 algorithms and narrowed down to 9 DEGs after taking the intersection. 4 hub genes were further screened by the GSE75214 validation set, namely COL1A1, CXCL10, MMP2 and FGF2. COL1A1 has the highest specificity and sensitivity for the diagnosis of CD and is considered to have the potential to diagnose CD. Five immune cells with significant differences were screened between CD and health controls (HC). Through the correlation analysis between five kinds of immune cells and four biomarkers, it was found that CXCL10 was positively correlated with activated dendritic cells, effector memory CD8+ T cells. MMP2 was positively correlated with activated dendritic cells, gamma delta T cells (γδ T) and mast cells. MMP2 and COL1A1 were significantly increased in colon tissue of CD fibrosis rats.

Conclusion

MMP2, COL1A1, CXCL10 and FGF2 can be used as hub genes for CD. Among them, COL1A1 can be used as a biomarker for the diagnosis of CD. MMP2 and CXCL10 may be involved in the development and progression of CD by regulating activated dendritic cell, effector memory CD8+ T cell, γδ T cell and mast cell. In addition, MMP2 and COL1A1 may be more closely related to CD intestinal fibrosis.

Inflammatory profile of keratoconic corneal epithelium

BACKGROUND

Recent studies have presented inflammatory features on keratoconus (KC) and many inflammatory markers are described in the tears of patients with this disease. The KC pathogenesis is still unknown just like the correlation with inflammatory patterns. However, environmental and genetic issues may be part of the progress of KC. In addition, some systemic features, such as allergy and obesity, seem to be related to the progression of KC. Our purpose was to evaluate the neuropeptides vasoactive intestinal peptide (VIP), neuropeptide Y (NPY), chemokines ligand 2 (CCL-2) and 5 (CCL-5), and interleukins 6 (IL-6) and 8 (IL-8) on corneal epithelial cells and blood of patients with KC and in healthy controls. In addition, the neutrophil-to-lymphocyte ratio (NLR) was evaluated to predict inflammation.

METHODS

This including prospective observational study included 32 KC patients who underwent corneal crosslinking (CXL) and 32 control patients who underwent photorefractive keratectomy (PRK). Patients' corneal epithelial cells were removed surgically, and blood (buffy coat) was analyzed. Samples in triplicate were evaluated on rt-PCR for neuropeptides (VIP e NPY), interleukins (IL-6 e IL-8), and chemokines (CCL-2 and CCL-5).

RESULTS

Our study showed statistically higher CCL-5 and IL-8 on corneal epithelial cells in patients with KC. Blood cells were statistically higher in VIP and NPY in the KC group. Interleukin-8 on blood cells was statistically significant in KC'S group; for CCL-2 and CCL-5 they were statistically lower in patients with KC compared with controls. NLR showed no difference between the groups.

CONCLUSIONS

Our data support the findings of other studies that suggested altering KC status, such as inflammatory corneal disease. The presence of IL-8 in the cornea and blood samples of KC's group suggested systemic disease with a possible local or repercussion action. Further studies are warranted to elucidate KC pathogenesis and its correlation to systemic disease.

Sequence variants contributing to dysregulated inflammatory responses across keratoconic cone surface in adolescent patients with keratoconus

Background

Keratoconus (KTCN) is the most common corneal ectasia resulting in a conical shape of the cornea. Here, genomic variation in the corneal epithelium (CE) across the keratoconic cone surface in patients with KTCN and its relevance in the functioning of the immune system were assessed.

Methods

Samples from four unrelated adolescent patients with KTCN and two control individuals were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated towards the whole-genome sequencing (WGS) study embracing a total of 18 experimental samples. The coding and non-coding sequence variation, including structural variation, was assessed and then evaluated together with the previously reported transcriptomic outcomes for the same CE samples and full-thickness corneas.

Results

First, pathway enrichment analysis of genes with identified coding variants pointed to "Antigen presentation" and "Interferon alpha/beta signaling" as the most overrepresented pathways, indicating the involvement of inflammatory responses in KTCN. Both coding and non-coding sequence variants were found in genes (or in their close proximity) linked to the previously revealed KTCN-specific cellular components, namely, "Actin cytoskeleton", "Extracellular matrix", "Collagen-containing extracellular matrix", "Focal adhesion", "Hippo signaling pathway", and "Wnt signaling" pathways. No genomic heterogeneity across the corneal surface was found comparing the assessed topographic regions. Thirty-five chromosomal regions enriched in both coding and non-coding KTCN-specific sequence variants were revealed, with a most representative 5q locus previously recognized as involved in KTCN.

Conclusion

The identified genomic features indicate the involvement of innate and adaptive immune system responses in KTCN pathogenesis.

The Impaired Wound Healing Process Is a Major Factor in Remodeling of the Corneal Epithelium in Adult and Adolescent Patients With Keratoconus

Purpose

Keratoconus (KTCN) is the most common corneal ectasia, characterized by pathological cone formation. Here, to provide an insight into the remodeling of the corneal epithelium (CE) during the course of the disease, we evaluated topographic regions of the CE of adult and adolescent patients with KTCN.

Methods

The CE samples from 17 adult and 6 adolescent patients with KTCN, and 5 control CE samples were obtained during the CXL and PRK procedures, respectively. Three topographic regions, central, middle, and peripheral, were separated toward RNA sequencing and MALDI-TOF/TOF Tandem Mass Spectrometry. Data from transcriptomic and proteomic investigations were consolidated with the morphological and clinical findings.

Results

The critical elements of the wound healing process, epithelial-mesenchymal transition, cell-cell communications, and cell-extracellular matrix interactions were altered in the particular corneal topographic regions. Abnormalities in pathways of neutrophils degranulation, extracellular matrix processing, apical junctions, IL, and IFN signaling were revealed to cooperatively disorganize the epithelial healing. Deregulation of the epithelial healing, G2M checkpoints, apoptosis, and DNA repair pathways in the middle CE topographic region in KTCN explains the presence of morphological changes in the corresponding doughnut pattern (a thin cone center surrounded by a thickened annulus). Despite similar morphological characteristics of CE samples in adolescents and adults with KTCN, their transcriptomic features were different. Values of the posterior corneal elevation differentiated adults with KTCN from adolescents with KTCN and correlated with the expression of TCHP, SPATA13, CNOT3, WNK1, TGFB2, and KRT12 genes.

Conclusions

Identified molecular, morphological, and clinical features indicate the effect of impaired wound healing on corneal remodeling in KTCN CE.

Corneal Biomechanical Characteristics in Osteogenesis Imperfecta With Collagen Defect

Purpose

To identify the characteristic corneal biomechanical properties of osteogenesis imperfecta (OI), and to compare the corneal biomechanical properties between OI and keratoconus.

Methods

We included 46 eyes of 23 patients with OI, 188 eyes of 99 keratoconus patients, and 174 eyes of 92 normal controls to compare corneal biomechanical parameters between OI corneas, keratoconus, and normal controls by using Corneal Visualization Scheimpflug Technology (Corvis ST).

Results

Patients with OI had significantly higher Corvis biomechanical index (CBI) (P < 0.001), higher tomographic and biomechanical index (TBI) (P = 0.040), lower Corvis Biomechanical Factor (CBiF) (P = 0.034), and lower stiffness parameter at first applanation (SP-A1) (P < 0.001) compared with normal controls. In contrast, OI group showed lower CBI (P < 0.001), lower TBI (P < 0.001), higher CBiF (P < 0.001), and higher SP-A1 (P = 0.020) than keratoconus group. Notably, the stress-strain index (SSI) was not significantly different between the OI and normal controls (P = 1.000), whereas keratoconus showed the lowest SSI compared with OI group (P = 0.025) and normal controls (P < 0.001).

Conclusions

Although the corneal structures of OI patients are less stable and easier to deform as compared to those of the control group, there is no significant difference in material stiffness observed between the OI and normal controls. In contrast, the corneas of keratoconus showed not only lower structural stability and higher deformability but also lower material stiffness compared with those of OI cornea and normal controls.

Translational Relevance

The biomechanical alterations are different between OI corneas and keratoconus.

Non-allergic eye rubbing is a major behavioral risk factor for keratoconus

Since the environmental, behavioral, and socioeconomic factors in the etiology of keratoconus (KTCN) remain poorly understood, we characterized them as features influencing KTCN phenotype, and especially affecting the corneal epithelium (CE). In this case-control study, 118 KTCN patients and 73 controls were clinically examined and the Questionnaire covering the aforementioned aspects was completed and then statistically elaborated. Selected KTCN-specific findings were correlated with the outcomes of the RNA-seq assessment of the CE samples. Male sex, eye rubbing, time of using a computer after work, and dust in the working environment, were the substantial KTCN risk factors identified in multivariate analysis, with ORs of 8.66, 7.36, 2.35, and 5.25, respectively. Analyses for genes whose expression in the CE was correlated with the eye rubbing manner showed the enrichment in apoptosis (TP53, BCL2L1), chaperon-related (TLN1, CTDSP2, SRPRA), unfolded protein response (NFYA, TLN1, CTDSP2, SRPRA), cell adhesion (TGFBI, PTPN1, PDPK1), and cellular stress (TFDP1, SRPRA, CAPZB) pathways. Genes whose expression was extrapolated to the allergy status didn't contribute to IgE-related or other inflammatory pathways. Presented findings support the hypothesis of chronic mechanical corneal trauma in KTCN. Eye-rubbing causes CE damage and triggers cellular stress which through its influence on cell apoptosis, migration, and adhesion affects the KTCN phenotype.

Different Murine High-Risk Corneal Transplant Settings Vary Significantly in Their (Lymph)angiogenic and Inflammatory Cell Signatures

Purpose

Pathologic conditions in the cornea, such as transplant rejection or trauma, can lead to corneal neovascularization, creating a high-risk environment that may compromise subsequent transplantation. This study aimed to evaluate the impact of different types of corneal injury on hemangiogenesis (HA), lymphangiogenesis (LA) and immune cell pattern in the cornea.

Methods

We used five different corneal injury models, namely, incision injury, alkali burn, suture placement, and low-risk keratoplasty, as well as high-risk keratoplasty and naïve corneas as control. One week after incision and 2 weeks after all other different injuries, corneal HA and LA were quantified by morphometric analysis. In addition, immune cell patterns of the whole cornea and the recipient rim were analyzed by immunohistochemistry. Immune cells in the draining lymph nodes (dLNs) were quantified by flow cytometry.

Results

Different types of corneal injury caused significantly different HA and LA responses (both P < 0.0001). The infiltration of corneal macrophages, dendritic cells, neutrophils, major histocompatibility complex (MHC) II+ cells, CD4+ T cells, and CD8+ T cells varied significantly in different high-risk settings (all P < 0.0001). Both the expression of MHC II on macrophages (P = 0.0005) and the frequency of MHC II+ dendritic cells (P = 0.0014) in the draining lymph nodes were significantly different across the various high-risk scenarios.

Conclusions

Murine high-risk settings caused by different underlying pathologies vary significantly in their (lymph)angiogenic and inflammatory cell patterns. Therefore, anti(lymph)angiogenic or immunomodulatory strategies to prevent and/or treat immune responses after subsequent corneal transplantation may need to be customized according to their immune-vascular "signatures."

Ocular Surface Immune Cell Profiles in Contact Lens-Induced Discomfort

Purpose

Contact lens–induced discomfort (CLD) remains a primary factor in discontinuation or prevention of contact lens wear. Thus, we investigated the role of ocular surface immune cells in subjects with CLD.

Methods

Habitual contact lens (CL) wearers with CLD (n = 19; 38 eyes) and without CLD (n = 21; 42 eyes) as determined by the Contact Lens Dry Eye Questionnaire-8 was included in a trial. Enrolled subjects used either of the two types of CL (designated as CL-A or CL-D). Ocular surface cells from the bulbar conjunctiva were obtained by impression cytology. The collected cells were phenotyped using fluorochrome-conjugated antibodies specific for leukocytes (CD45⁺), neutrophils (CD66b^+,High,Low), macrophages (CD163⁺), T cells (CD3⁺CD4⁺, CD3⁺CD8⁺), natural killer (NK) cells (CD56^{+, High, Low}), natural killer T (NKT) cells (CD3⁺CD56⁺), and gamma delta T (γδT) cells (CD3⁺γδTCR⁺) by flow cytometry. Further, corneal dendritic cell density (cDCD) was also determined using in vivo confocal microscopy.

Results

Significantly higher proportions of CD45⁺ cells were observed in subjects with CLD compared to those without CLD. The percentages of CD66b^Total,Low, CD163⁺, pan T cells, CD4⁺T cells, CD8⁺T cells, CD56^{Total,High,Low} (NK) cells, and NKT cells, as well as the CD4/CD8 ratio, were significantly higher in CLD subjects. The proportion of T cells (CD4, CD8, CD4/CD8 ratio, NKT cells) and macrophages exhibited a direct association with discomfort score. The percentages of CD45⁺, CD66b^Total,Low, CD163⁺, CD3⁺, CD56^{Total,High,Low}, and NKT cells and cDCD were significantly higher in CLD subjects wearing CL-D. The percentages of CD66b^High, CD4⁺T cells, CD8⁺T cells, NKT cells, and CD4/CD8 ratio were significantly higher in CLD subjects wearing CL-A.

Conclusions

Increased proportions of ocular surface immune cells are observed in CLD, and the lens type could impact the immune cells associated with CLD.

Translational Relevance

The association between the proportion of altered ocular surface immune cell subsets and contact lens–induced discomfort underpins the importance of considering immune-related aspects during contact lens development and in the clinical management of ocular surface pain.