Abstract
Keratoconus is a corneal ectasia whose pathophysiological mechanisms, including biomolecular alterations and genetic influences, remain poorly understood. Recent studies have shown altered cytokine levels, increased proteinase activity, and other potential mediators in the tear film and corneal tissue, highlighting a possible involvement of inflammatory pathways in the pathophysiology of keratoconus. This observational study aims to characterize the tear proteome of keratoconus patients and compare it to a control group, reporting potential disease biomarkers in the tear film. 23 keratoconus patients were selected at the Cornea and External Diseases Outpatient Clinic of the Clinics Hospital of UNICAMP. The control group consisted of 17 age- and sex-matched participants. All study subjects underwent corneal tomography (Pentacam). Tear film samples were collected and sent for proteomic evaluation by mass spectrometry at the National Biosciences Laboratory (LNBio). After quantification, univariate and multivariate statistical analyses were performed. A total of 353 proteins were identified and quantified, of which 25 showed statistical differences in the univariate analysis (t-test), and 19 were selected in the multivariate analysis (PLS-DA). There was an overlap of 7 proteins identified in both uni- and multivariate analyses: chitinase-3-like protein 2, prosaposin, zymogen granule protein 16 homolog B, procollagen-lysine,2-oxoglutarate 5-dioxygenase 1, secretoglobin family 1D member 1, albumin, and Ig kappa chain V-I region. Thirty-seven proteins showed statistically significant variation between the keratoconus and control groups. Proteomic analysis revealed differentially expressed proteins in the tear film of keratoconus patients. We report the identified proteomic profile, which includes potential biomarkers that may help elucidate the disease’s pathophysiology.
Impact statement
This study provides a comprehensive proteomic characterization of the tear film in keratoconus, identifying differentially expressed proteins that reveal novel insights into the disease’s molecular mechanisms. By combining univariate and multivariate statistical analyses, including PLS-DA, we demonstrated a distinct proteomic signature that discriminates keratoconus from controls, implicating inflammatory, oxidative stress, and extracellular matrix remodeling pathways. These findings advance the field by expanding the catalog of candidate biomarkers and highlighting molecular pathways that may contribute to disease onset and progression. The data support the concept that keratoconus is not merely a structural ectasia but a complex disorder involving active biochemical dysregulation. This work provides a valuable proteomic resource for the ophthalmic research community and lays the groundwork for future translational studies aimed at early diagnosis and the development of targeted therapeutic strategies.
Introduction
Keratoconus is a corneal ectasia defined by progressive corneal thinning and protrusion, leading to irregular astigmatism and varying degrees of visual impairment. It typically manifests during puberty and progresses until the third or fourth decade of life. It is usually bilateral, often asymmetrical, though rarely presented unilaterally [1–3]. Keratoconus remains one of the leading indications for corneal transplantation worldwide [2] and has a significant financial impact as it affects economically active individuals [4].
A recent meta-analysis, which included 29 studies and a total population of over 50 million individuals from 15 different countries, found a global prevalence of 1.38 cases per 1,000 inhabitants [4]. Prevalence rates vary widely, reflecting differences in sample sizes, diagnostic methodologies, and the influences of genetic and ethnic factors in keratoconus presentation [4]. No population-based studies have been found on the prevalence of keratoconus in Brazil.
The main risk factors associated with keratoconus include eye rubbing, allergy, asthma, eczema, and a positive family history [3, 4], highlighting genetic and environmental influences. Approximately 8%–10% of keratoconus cases have a positive family history [2, 3]. Substantial evidence supports a genetic influence in keratoconus pathogenesis [2, 5, 6]. Twin studies provide insights into the relative contributions of genotype and environment to disease phenotypic expression. Keratoconus in monozygotic twins demonstrates high concordance with variable expressivity [1, 7, 8]. Studies in families of keratoconus patients have revealed topographic and tomographic abnormalities even in asymptomatic individuals [9, 10]. Indeed, a first-degree relative with keratoconus is considered a major risk factor for disease development [10]. Recent studies have identified over 20 genes implicated in keratoconus susceptibility [11]. Keratoconus is currently regarded as a multifactorial disease, where multiple genetic factors interact with environmental influences to determine its clinical presentation [12]. Recent proteomic work in offspring of patients with keratoconus further supports this concept, demonstrating early molecular alterations, particularly involving oxidative stress responses, cytoskeletal organization, and mechanotransduction, even before clinical or biomechanical abnormalities become detectable [13].
Keratoconus is a disease with a significant socio-economic impact due to its relatively high prevalence and its effect on an economically active age group. Despite being described nearly 300 years ago, the detailed pathophysiological, genetic, and environmental mechanisms involved in its onset and progression remain unclear. Current treatment is based on visual rehabilitation, ranging from optical aids such as glasses and rigid or scleral contact lenses to surgical interventions, including intrastromal rings and corneal transplantation [3]. Corneal crosslinking aims to halt disease progression by altering the structural properties of stromal collagen. The only preventive approach is controlling known risk factors, such as eye rubbing and ocular allergy. No pharmacological treatment is available to interfere with the altered metabolic pathways of the disease, nor is there a susceptibility test to identify individuals at risk of developing keratoconus or those already diagnosed at risk of progression [3].
Although keratoconus has been initially described as a noninflammatory disease due to the absence of clinical signs of inflammation, such as conjunctival hyperemia, corneal infiltrates, or anterior chamber reaction [1, 14], later studies have reported significant alterations in inflammatory mediators, including increased levels of pro-inflammatory cytokines, collagen degradation enzymes such as matrix metalloproteinases, indicating that keratoconic corneas exhibit some degree of inflammation [12, 15–18]. Oxidative stress markers and antioxidant systems are also dysregulated in keratoconus. Evidence suggests an increase in oxidative stress markers, particularly reactive oxygen species (ROS) and reactive nitrogen species (RNS), alongside a reduction in antioxidants such as ALDH/NADPH dehydrogenase, lactoferrin, transferrin, albumin, selenium, zinc, vitamin B12, and folic acid, among others [19, 20]. Large-scale proteomic studies have reinforced this inflammatory and oxidative-stress profile by identifying dysregulation of tear proteins involved in glycolytic pathways, reactive oxygen detoxification, and inflammatory regulation across different disease stages, including cystatin-S, lacritin, glutathione synthetase, and superoxide dismutase [21].
Proteomic analysis of human tissues and fluids has emerged as one of the most relevant recent approaches in biomarker research. The tear film has garnered increasing interest in recent years as a potential source of biomarkers for various diseases due to its accessibility, moderate complexity, and responsiveness to both ocular and systemic diseases [22]. Searching for biomarkers with high sensitivity and specificity for a given disease is crucial to improving diagnostic methods, identifying cellular and metabolic alterations that may elucidate underlying pathophysiological mechanisms, and providing potential therapeutic targets [23].
There is growing evidence that keratoconus exhibits a characteristic proteome [24]. Collagen and other structural proteins, such as lumican, keratan sulfate, and decorin, are decreased [24]. Conversely, there is an increased expression of degradative enzymes, including phosphatases, lipases, esterases, cathepsins, and matrix metalloproteinases, as well as elevated levels of pro-inflammatory proteins such as interleukin 1 (IL-1), interleukin 6 (IL-6), matrix metalloproteinase 9 (MMP-9), transforming growth factor beta (TGF-ß), and tumor necrosis factor alpha (TNF-α) [16, 17, 25]. Complementary evidence indicates that impaired epithelial wound healing, dysregulated epithelial-mesenchymal transition pathways, and altered cytokine signaling contribute to the characteristic topographic changes and epithelial remodeling observed in keratoconus [26].
In a previous study [27], we demonstrated that mass spectrometry-based proteomics performed on tear samples was able to differentiate three distinct diseases: keratoconus, pterygium, and dry eye secondary to graft-versus-host disease (GVHD). Each disease exhibited a characteristic proteomic profile, identifiable through multivariate statistical analysis methods. Furthermore, we identified the main differentially expressed proteins in each group compared to the control, which were reported as potential biomarkers for each disease. In our previous evaluation, the keratoconus group consisted of four samples, and the control group included six samples, which have been expanded in the present study.
Despite recent discoveries, the exact mechanisms initiating the cellular and molecular alterations that culminate in corneal degradation and shape distortion in keratoconus remain unknown. Furthermore, the interactions between genetic and environmental factors in modulating these alterations are yet to be fully elucidated. Tear proteomic analysis can be an important tool in the search for biomarkers in ocular diseases and has been used to investigate the pathophysiology of keratoconus.
This study aims to quantify and report differentially expressed proteins in the tear film of patients with keratoconus compared to a control group.
Materials and methods
This is a cross-sectional, observational, and non-interventional study. The research subjects were divided into two groups: 23 patients diagnosed with keratoconus under follow-up at the Cornea and External Diseases Outpatient Clinic of the Clinics Hospital of UNICAMP, and 17 control subjects. All research subjects underwent a comprehensive ophthalmological examination and corneal tomography using the OCULUS Pentacam® software 1.20r134. The diagnosis in the Keratoconus group was confirmed by identifying characteristic signs, including increased corneal curvature, stromal thinning, alterations in elevation maps, and irregular astigmatism. Corneal tomography was also used in the Control group to confirm the absence of the disease.
Patients diagnosed with keratoconus and under follow-up at the Keratoconus Outpatient Clinic of HC/UNICAMP who consented to participate in the study were included. For the control group, volunteers were recruited from hospital personnel, university students, and family members of patients from other ambulatories, all without clinical or tomographic signs of keratoconus, without a family history of the disease, and without any ocular surface pathology. Subjects were excluded if they were using anti-inflammatory, immunosuppressive, or immunomodulatory medications, either topical or systemic, or had active ocular inflammatory or infectious conditions at the time of data collection. Additionally, individuals with a history of ocular surface surgery, including laser refractive surgery, radial keratotomy, intrastromal ring implantation, corneal crosslinking, pterygium excision, cataract surgery, among others, were excluded, as well as previous trauma or signs of traumatic corneal scarring.
Tear samples were collected using microcapillary pipettes with atraumatic contact with the lower tear meniscus. The samples were then transferred to cryotubes and frozen at −80 °C. The samples (17 controls and 23 keratoconus cases) were later prepared and processed as previously described [27, 28].
Following data acquisition, processing was performed using MaxQuant software version 1.3.0.3 with the Andromeda algorithm against the UniProt Human Protein Database (downloaded in May 2019, containing 95,542 sequences and 38,078,700 residues). Bioinformatics analyses were conducted using Perseus software version 1.5.1.6. Logarithmic transformation was applied, and filters removed reverse sequences and proteins identified by only one modified peptide.
Mass spectrometry data underwent logarithmic transformation before statistical analysis. Univariate analyses were conducted using GraphPad Prism version 6.00. Measurements from keratoconus patient samples were compared with those of control samples using an unpaired t-test, both with and without correction for multiple analyses: false discovery rate 5% (FDR). Multivariate analyses were performed using the online platform MetaboAnalyst1 [29].
Results
A total of 40 individuals were evaluated, divided into a keratoconus patient group (n = 23) and a control group (n = 17). The main clinical characteristics of each group are shown in Table 1.
TABLE 1
| Variables | Keratoconus | Control | P-value |
|---|---|---|---|
| Sample size | 23 | 17 | |
| Sex (female/Male) | 10/13 | 11/6 | 0.184 |
| Age (years) | 21.65 | 24.47 | 0.075 |
| Kmax (diopters) | 57.7D | 44.9D | <0.001 |
| Pachymetry (µm) | 454 | 566 | <0.001 |
| Belin D (pentacam) | 9.56 | 0.94 | <0.001 |
| Corneal astigmatism (diopters) | 4.51 | 1.16 | <0.001 |
| Ocular allergy | 95.6% | 23.5% | <0.001 |
Clinical characteristics.
A total of 353 proteins were identified and quantified, of which 25 showed statistically significant differences in the univariate analysis using the t-test with p < 0.05 (Table 2; Figure 1), and 19 were selected in the multivariate partial least-squares discriminant analysis (PLS-DA) (Figure 2). Seven proteins were identified in both the uni- and multivariate analyses: chitinase-3-like protein 2, prosaposin, zymogen granule protein 16 homolog B, procollagen-lysine,2-oxoglutarate 5-dioxygenase 1, secretoglobin family 1D member 1, albumin, and Ig kappa chain V-I region. Thirty-seven proteins exhibited statistically significant variation between the keratoconus and control groups.
TABLE 2
| Gene | UniProt ID | Protein name | P-value | Control/KC ratio |
|---|---|---|---|---|
| DNAH5 | Q8TE73 | Dynein heavy chain 5, axonemal | 0.0011 | 0.54 |
| PEBP4 | Q96S96 | Phosphatidylethanolamine-binding protein 4 | 0.0017 | 0.55 |
| LSR | Q86X29 | Lipolysis-stimulated lipoprotein receptor | 0.0047 | 0.50 |
| CHI3L2 | Q15782 | Chitinase-3-like protein 2 | 0.0049 | 0.62 |
| CP | P00450 | Ceruloplasmin | 0.005 | 0.76 |
| PSAP | P07602 | Prosaposin | 0.005 | 0.74 |
| LPO | P22079 | Lactoperoxidase | 0.006 | 0.65 |
| ZG16B | Q96DA0 | Zymogen granule protein 16 homolog B | 0.0069 | 0.69 |
| SERPINA3 | P01011 | Alpha-1-antichymotrypsin | 0.0083 | 0.52 |
| SCGB1D1 | O95968 | Secretoglobin family 1D member 1 | 0.0088 | 0.49 |
| CST5 | P28325 | Cystatin-D | 0.0134 | 18.04 |
| PLA2G2A | P14555 | Phospholipase A2, membrane-associated | 0.014 | 0.60 |
| ANXA5 | P08758 | Annexin A5 | 0.0145 | 1.77 |
| MUCL1 | Q96DR8 | Mucin-like protein 1 | 0.0164 | 1.31 |
| PLOD1 | Q02809 | Procollagen-lysine,2-oxoglutarate 5-dioxygenase 1 | 0.0187 | 0.73 |
| DAG1 | Q14118 | Dystroglycan; Alpha-dystroglycan; Beta-dystroglycan | 0.0192 | 0.70 |
| GOLM1 | Q8NBJ4 | Golgi membrane protein 1 | 0.0237 | 0.72 |
| SCGB2A1 | O75556 | Mammaglobin-B | 0.0253 | 0.67 |
| GANAB | Q14697 | Neutral alpha-glucosidase AB | 0.0253 | 0.72 |
| Ig kappa | P01597 | Ig kappa chain V-I region | 0.0265 | 1.65 |
| SPRR3 | Q9UBC9 | Small proline-rich protein 3 | 0.0284 | 1.95 |
| CST2 | P09228 | Cystatin-SA | 0.031 | 12.11 |
| S100A7 | P31151 | Protein S100-A7 | 0.032 | 2.38 |
| ALB | P02768 | Serum albumin | 0.0359 | 1.79 |
| HYOU1 | Q9Y4L1 | Hypoxia up-regulated protein 1 | 0.0433 | 0.54 |
Differentially expressed proteins after t-test with P < 0.05.
FIGURE 1

Tear proteins (identified by gene name) with a p-value <0.05 after Student's t-test without FDR correction. After FDR correction, no protein reached p-value <0.05. Data distribution with standard deviation and mean are displayed. Blue: control. Red: keratoconus patients.
FIGURE 2
![Four-part figure showing data analysis related to keratoconus:A) Scatter plot of t[1] and t[2] scores with red and blue dots representing keratoconus and control groups, respectively, within an ellipse.B) Bar chart with green bars showing VIP[2] scores for various proteins, with error bars.C) Scatter plot correlating original and permuted Y-vectors, with green and blue points for R2Y(Cum) and Q2(Cum).D) Dendrogram clustering control (blue) and keratoconus (red) samples based on t[1], t[2].Legend included for color coding.](https://www.ebm-journal.org/files/Articles/10864/xml-images/ebm-250-10864-g002.webp)
PLS-DA and hierarchical clustering analysis. No data imputation was performed, and variables (proteins) with missing values were excluded. (A) PLS-DA score plots of keratoconus and control groups: 2 latent variables, R2X = 0.433, R2Y = 0.822, Q2 = 0.720. (B) Variable influence on this classification (VIP–variable importance in the prediction). (C) Permutation test of PLS-DA models, intercepts: R2=(0.0, 0.232), Q2=(0.0, −0.263). (D) Hierarchical clustering analysis based on the variables selected in the PLS-DA test, shown in Panel B.
Discussion
The present study evaluated the tear proteome of patients with keratoconus compared to a control group without the disease using mass spectrometry and reported differentially expressed proteins between the groups.
Proteomics studies have been used to identify biological markers across various medical fields, including ophthalmology. In 2020, our group published a study [27] on proteomic analysis of tears from patients with keratoconus, pterygium, and dry eye associated with GVHD, in which we reported the tear proteome for each disease and the differentially expressed proteins in each group compared to controls. Multivariate statistical analyses of tear proteome distinguished each distinct ocular disease by a characterizable proteomic profile. Recently, a meta-analysis of candidate proteins associated with keratoconus [27] analyzed 346 normal and 493 keratoconus eyes. Altered proteins involved in inflammatory pathways, extracellular matrix remodeling, and apoptosis were reported. The main proteins identified were MMP-9, IL-6, lysyl oxidase (LOX), TNF, and IL-1B.
Below is a discussion of the study’s main findings, based on the UniProt online database2 and literature review. Table 3 summarizes the biological functions of the main proteins.
TABLE 3
| Gene | Protein name | Main biological function |
|---|---|---|
| PEBP4 | Phosphatidylethanolamine-binding protein 4 | Modulates PI3K-AKT signaling; involved in cell survival |
| CHI3L2 | Chitinase-3-like protein 2 | Binds to carbohydrate structures, may play a role in immune response and tissue remodeling |
| CP | Ceruloplasmin | Ferroxidase activity; antioxidant defense; iron and copper metabolism |
| PSAP | Prosaposin | Lysosomal degradation of sphingolipids; neurotrophic activity |
| LPO | Lactoperoxidase | Antimicrobial defense; oxidative stress marker |
| ZG16B | Zymogen granule protein 16 homolog B | Secreted protein; proposed role in sustaining pro-inflammatory environments |
| SERPINA3 | Alpha-1-antichymotrypsin | Serine protease inhibitor; anti-inflammatory, antifibrotic, antioxidant roles |
| SCGB1D1 | Secretoglobin family 1D member 1 (Lipophilin A) | Modulates inflammation and tissue repair |
| SCGB2A1 | Mammaglobin B (Lipophilin C) | Secretory protein; androgen-binding; component of tear film protein complexes |
| CST5 | Cystatin D | Inhibitor of cysteine proteases; protects extracellular matrix from degradation |
| CST2 | Cystatin SA | Regulates proteolytic activity; contributes to tissue homeostasis |
| CST3 | Cystatin C | Key inhibitor of cysteine proteases; modulates extracellular proteolysis |
| PLA2G2A | Secretory phospholipase A2 | Hydrolyzes phospholipids; initiates arachidonic acid-mediated inflammatory pathways |
| ANXA5 | Annexin A5 | Epithelial repair; ECM structural organization |
| MUCL1 | Mucin-like 1 | Tear film lubrication and ocular surface protection |
| PLOD1 | Procollagen-lysine,2-oxoglutarate 5-dioxygenase 1 | Collagen fiber assembly and crosslinking |
| DAG1 | Dystroglycan 1 | Basement membrane organization; epithelial adhesion and signaling |
| SPRR3 | Small proline-rich protein 3 | Terminal epithelial differentiation; structural support |
| ALB | Albumin | Metal transport; antioxidant buffering |
| S100A7 | Protein S100-A7 (psoriasin) | Inflammation, angiogenesis, oxidative stress regulation |
| HYOU1 | Hypoxia up-regulated protein 1 | Endoplasmic reticulum stress response; cytoprotection under hypoxia/inflammation |
Main biological functions of differentially expressed proteins identified in the tear film of keratoconus patients.
The proteomic profile in the tear film of keratoconus patients shows dysregulation of proteins involved in inflammation, oxidative stress, tissue remodeling, and extracellular matrix homeostasis. These mechanisms are widely implicated in the disease’s pathophysiology. Rather than evaluating each differentially expressed protein in isolation, a pattern emerges in which multiple proteins converge on interconnected biological processes that may underlie the corneal thinning, epithelial instability, and stromal degradation seen in keratoconus.
Several upregulated proteins, including ceruloplasmin, lactoperoxidase, and prosaposin, suggest an enhanced inflammatory and oxidative environment on the ocular surface. Increased oxidative stress markers have been consistently reported in keratoconus corneas [19, 30]. Our findings reinforce this biochemical signature. The upregulation of enzymes such as phospholipase A2 further suggests a heightened pro-inflammatory state [31, 32], consistent with studies reporting increased levels of inflammatory mediators in dry eye disease and in subsets of keratoconus patients with allergic comorbidities. Albumin, however, showed decreased levels in the keratoconus group, consistent with previous studies [19, 33–35]. Its reduced presence in tears may increase susceptibility to oxidative stress [19]. S100 A7 protein is linked to ocular surface inflammation [36] and recurrent pterygium [37], but its levels were lower in the keratoconus group.
Proteins associated with tissue remodeling also showed dysregulation; higher levels of procollagen-lysine,2-oxoglutarate 5-dioxygenase 1, and dystroglycan 1 suggest altered extracellular matrix turnover and basement membrane dynamics [38], processes already implicated in epithelial fragility and stromal biomechanical weakening in keratoconus. Proteins with potential roles in maintaining epithelial integrity, such as Annexin A5 [39, 40], mucin-like protein 1 [41], and small proline-rich protein 3 [42], were reduced in keratoconus tears, consistent with previous reports of compromised epithelial barrier function in these patients. Cystatins (CST2, CST3, CST5), key inhibitors of cysteine proteases, were also decreased, potentially contributing to increased proteolytic activity and stromal degradation—an established pathogenic mechanism in keratoconus [43–45]. SerpinA3, which has been demonstrated to exhibit anti-inflammatory, anti-angiogenic, antioxidant, and anti-fibrotic activities [46], was upregulated in the keratoconus group, perhaps as a response to the proinflammatory environment. The behavior of secretoglobins [47], however, differed from prior studies that reported reduced levels [43, 48]; their elevation in our cohort may reflect heterogeneity among patient populations or differences in disease stage, severity, or environmental exposures. Secretoglobin 2A1 has previously been reported to be upregulated in patients with keratoconus [49].
Some proteins associated with known risk factors for keratoconus were upregulated in the disease group. Hypoxia up-regulated protein 1, which is correlated with vernal keratoconjunctivitis [50], may be associated with the high prevalence of ocular allergy in this population. Saposins are associated with the sphingolipid metabolic pathway, which has been linked to eye-rubbing behavior [51].
Some proteins, such as chitinase-3-like protein 2 and zymogen granule protein 16 homolog B, remain poorly characterized in ocular physiology [49, 52] but were consistently upregulated in our analysis. Their repeated identification across proteomic studies suggests they may represent underexplored components of the inflammatory and remodeling responses in keratoconus.
Together, these differentially expressed proteins support a multifactorial disease model in which chronic inflammation, oxidative imbalance, epithelial barrier disruption, and extracellular matrix instability act synergistically to promote corneal thinning. Continued investigation into these pathways may help refine biomarkers for diagnosis and progression monitoring, as well as identify new therapeutic targets.
The PLS-DA model constructed to discriminate between keratoconus and control tear proteomes demonstrated strong explanatory and predictive performance. Using two latent variables, the model achieved R2X = 0.433, R2Y = 0.822, and Q2 = 0.720. The proportion of explained variance in the predictor matrix (R2X) indicates that the latent variables efficiently captured the most relevant structure within the proteomic dataset. The high R2Y value demonstrates that the model accounts for the majority of variance associated with group separation, supporting that the discrimination observed is largely driven by biological differences rather than random noise. Moreover, the Q2 value of 0.720 confirms the model’s high predictive capacity, suggesting that it can reliably classify new samples according to disease status.
Model validity was further supported by permutation testing, which yielded intercepts of R2 = 0.232 and Q2 = −0.263. The low R2 and negative Q2 intercepts from randomized models indicate that the original model’s performance was not due to overfitting, but rather reflects genuine structure–response relationships within the data. Together, these findings provide strong evidence that the proteomic profiles of keratoconus and control groups are distinctly segregated in the multivariate space, underscoring the robustness and biological relevance of the discriminatory model.
The main limitation of this study is that the findings from the univariate analysis did not remain statistically significant after FDR correction. This is due to the small sample size, which reduces the study’s power. However, the literature shows that proteomics studies using mass spectrometry rarely have larger sample sizes due to such analyses’ high cost and time-consuming nature. It is also important to make the current findings available, including raw proteomics data, as these can be used in meta-analyses and bioinformatics tools to contribute to a larger database.
Another limitation is that no filtering of differentially expressed proteins was performed based on the percentage of missing values among subjects or fold-change thresholds.
Our study demonstrated altered levels of several proteins related to inflammatory pathways and oxidative stress, which agrees with recent findings in the literature. Building on our previous study, we conducted a second proteomic evaluation of tear samples from keratoconus patients, with a sample size four times larger, providing further data for future biomarker analyses and aiding the scientific community in unraveling the intricate pathophysiological mechanisms involved in keratoconus. This study represents a significant step forward in elucidating the complex pathophysiological mechanisms underlying keratoconus. Identifying differentially expressed proteins in the tear film might serve as insight for potential biomarkers. By expanding the understanding of molecular alterations associated with the disease, our findings reinforce the role of inflammation and oxidative stress in keratoconus progression and pave the way for future translational research aimed at improving early diagnosis, risk stratification, and the development of targeted therapeutic interventions. Given the socio-economic burden of keratoconus and the current lack of disease-modifying treatments, identifying specific proteomic signatures in tear fluid holds immense potential for transforming clinical management in corneal ectatic diseases.
Conclusion
Proteomic analysis revealed differentially expressed proteins in the tear film of keratoconus patients. We report the identified proteomic profile, which includes potential biomarkers that may help elucidate the disease’s pathophysiology, providing data for future studies.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics statement
The studies involving humans were approved by Research Ethics Committee of UNICAMP. The studies were conducted in accordance with the local legislation and institutional requirements. The participants provided their written informed consent to participate in this study.
Author contributions
DA did data collection; RD and DA performed experiments; MRA, DA, and MA performed data analysis and interpretation; DA, MRA, and MA drafted the article. All authors contributed to the article and approved the submitted version.
Funding
The author(s) declared that financial support was received for this work and/or its publication. The study was funded by FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo), Brazil, JP1 2014/19138-5, JP2 2021/04045-5.
Acknowledgments
We acknowledge the Mass Spectrometry Laboratory at the Brazilian Biosciences National Laboratory (LNBio), Brazilian Center for Research in Energy and Materials (CNPEM), Campinas, Brazil, for their support with the mass spectrometry analysis.
Conflict of interest
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Generative AI statement
The author(s) declared that generative AI was not used in the creation of this manuscript.
Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.ebm-journal.org/articles/10.3389/ebm.2025.10864/full#supplementary-material
References
1.
Rabinowitz YS . Keratoconus Surv Ophthalmol (1998) 42(4):297–320.
2.
Mas Tur V MacGregor C Jayaswal R O’Brart D Maycock N . A review of keratoconus: diagnosis, pathophysiology, and genetics. Surv Ophthalmol (2017) 62(6):770–83. 10.1016/j.survophthal.2017.06.009
3.
Gomes JAP Rapuano CJ Belin MW Ambrósio R , Group of Panelists for the Global Delphi Panel of Keratoconus and Ectatic Diseases. Global consensus on keratoconus diagnosis. Cornea (2015) 34(12):e38–9. 10.1097/ICO.0000000000000623
4.
Hashemi H Heydarian S Hooshmand E Saatchi M Yekta A Aghamirsalim M et al The prevalence and risk factors for keratoconus: a systematic review and meta-analysis. Cornea (2020) 39(2):263–70. 10.1097/ICO.0000000000002150
5.
Sugar J Macsai MS . What causes keratoconus?Cornea (2012) 31(6):716–9. 10.1097/ICO.0b013e31823f8c72
6.
Abu-Amero KK Al-Muammar AM Kondkar AA . Genetics of keratoconus: where do we stand?J Ophthalmol (2014) 2014:641708. 10.1155/2014/641708
7.
Loukovitis E Sfakianakis K Syrmakesi P Tsotridou E Orfanidou M Bakaloudi DR et al Genetic aspects of keratoconus: a literature review exploring potential genetic contributions and possible genetic relationships with comorbidities. Ophthalmol Ther (2018) 7(2):263–92. 10.1007/s40123-018-0144-8
8.
Tuft SJ Hassan H George S Frazer DG Willoughby CE Liskova P . Keratoconus in 18 pairs of twins. Acta Ophthalmol (2012) 90(6):e482–e486. 10.1111/j.1755-3768.2012.02448.x
9.
Steele TM Fabinyi DC Couper TA Loughnan MS . Prevalence of orbscan II corneal abnormalities in relatives of patients with keratoconus. Clin Exp Ophthalmol (2008) 36(9):824–30. 10.1111/j.1442-9071.2009.01908.x
10.
Lapeyre G Fournie P Vernet R Roseng S Malecaze F Bouzigon E et al Keratoconus prevalence in families: a French study. Cornea (2020) 39(12):1473–9. 10.1097/ICO.0000000000002546
11.
Bykhovskaya Y Rabinowitz YS . Update on the genetics of keratoconus. Exp Eye Res (2021) 202(October 2020):108398. 10.1016/j.exer.2020.108398
12.
Santodomingo-Rubido J Carracedo G Suzaki A Villa-Collar C Vincent SJ Wolffsohn JS . Keratoconus: an updated review. Contact Lens and Anterior Eye (2022) 45(3):101559. 10.1016/j.clae.2021.101559
13.
López-López M Regueiro U Bravo S Pena C Pastoriza Y Conde-Amboage M et al Tear proteomic analysis from offspring of keratoconus patients: new insights into corneal biomechanical weakness and disease risk stages. Invest Ophthalmol Vis Sci (2025) 66(5):41. 10.1167/iovs.66.5.41
14.
Krachmer JH Feder RS Belin MW . Keratoconus and related noninflammatory corneal thinning disorders. Surv Ophthalmol (1984) 28(4):293–322. 10.1016/0039-6257(84)90094-8
15.
Wisse RPL Kuiper JJW Gans R Imhof S Radstake TRDJ Van Der Lelij A . Cytokine expression in keratoconus and its corneal microenvironment: a systematic review. Ocul Surf (2015) 13(4):272–83. 10.1016/j.jtos.2015.04.006
16.
Jun AS Cope L Speck C Feng X Lee S Meng H et al Subnormal cytokine profile in the tear fluid of keratoconus patients. PLoS One (2011) 6(1):1–8. 10.1371/journal.pone.0016437
17.
Lema I Sobrino T Durán JA Brea D Díez-Feijoo E . Subclinical keratoconus and inflammatory molecules from tears. Br J Ophthalmol (2009) 93(6):820–4. 10.1136/bjo.2008.144253
18.
Balasubramanian SA Mohan S Pye DC Willcox MDP . Proteases, proteolysis and inflammatory molecules in the tears of people with keratoconus. Acta Ophthalmol (2012) 90(4):303–9. 10.1111/j.1755-3768.2011.02369.x
19.
Navel V Malecaze J Pereira B Baker JS Malecaze F Sapin V et al Oxidative and antioxidative stress markers in keratoconus: a systematic review and meta-analysis. Acta Ophthalmol (2021) 99(6):e777–94. 10.1111/aos.14714
20.
Monteiro de Barros MR Chakravarti S . Pathogenesis of keratoconus: NRF2-antioxidant, extracellular matrix and cellular dysfunctions. Exp Eye Res (2022) 219(March):109062. 10.1016/j.exer.2022.109062
21.
Kannan R Shetty R Panigrahi T Koh SK Khamar P Deshpande V et al Untargeted tear proteomics in a large South-Asian cohort reveals inflammatory signaling, ECM remodeling, and altered metabolism in Keratoconus. Invest Ophthalmol Vis Sci (2025) 66(2):60. 10.1167/iovs.66.2.60
22.
Ponzini E Santambrogio C De Palma A Mauri P Tavazzi S Grandori R . Mass spectrometry-based tear proteomics for noninvasive biomarker discovery. Mass Spectrom Rev (2022) 41(5):842–60. 10.1002/mas.21691
23.
Khanna RK Catanese S Emond P Corcia P Blasco H Pisella PJ . Metabolomics and lipidomics approaches in human tears: a systematic review. Surv Ophthalmol (2022) 67(4):1229–43. 10.1016/j.survophthal.2022.01.010
24.
Ferrari G Rama P . The keratoconus enigma: a review with emphasis on pathogenesis. Ocul Surf (2020) 18:363–73. 10.1016/j.jtos.2020.03.006
25.
Lema I Durán JA . Inflammatory molecules in the tears of patients with keratoconus. Ophthalmology (2005) 112(4):654–9. 10.1016/j.ophtha.2004.11.050
26.
Jaskiewicz K Maleszka-Kurpiel M Matuszewska E Kabza M Rydzanicz M Malinowski R et al The impaired wound healing process is a major factor in remodeling of the corneal epithelium in adult and adolescent patients with Keratoconus. Invest Ophthalmol Vis Sci (2023) 64(2):22. 10.1167/iovs.64.2.22
27.
de Almeida Borges D Alborghetti MR Franco Paes Leme A Ramos DR Duarte B Veiga M et al Tear proteomic profile in three distinct ocular surface diseases: keratoconus, pterygium, and dry eye related to graft-versus-host disease. Clin Proteomics (2020) 17(1):1–16. 10.1186/s12014-020-09307-5
28.
Winck FV Ribeiro ACP Domingues RR Ling LY Riaño-Pachón DM Rivera C et al Insights into immune responses in oral cancer through proteomic analysis of saliva and salivary extracellular vesicles. Sci Rep (2015) 5(November):1–13.
29.
Xia J Broadhurst DI Wilson M Wishart DS . Translational biomarker discovery in clinical metabolomics: an introductory tutorial. Metabolomics (2013) 9(2):280–99. 10.1007/s11306-012-0482-9
30.
Arnal E Peris-Martínez C Menezo JL Johnsen-Soriano S Romero FJ . Oxidative stress in keratoconus?Invest Ophthalmol Vis Sci (2011) 52(12):8592–7. 10.1167/iovs.11-7732
31.
Wei Y Epstein SP Fukuoka S Birmingham NP Li XM Asbell PA . sPLA2-IIa amplifies ocular surface inflammation in the experimental dry eye (DE) BALB/c mouse model. Invest Ophthalmol Vis Sci (2011) 52(7):4780–8. 10.1167/iovs.10-6350
32.
Soria J Acera A Merayo-Lloves J Durán JA González N Rodriguez S et al Tear proteome analysis in ocular surface diseases using label-free LC-MS/MS and multiplexed-microarray biomarker validation. Sci Rep (2017) 7(1):1–15. 10.1038/s41598-017-17536-2
33.
Balasubramanian SA Pye DC Willcox MDP . Levels of lactoferrin, secretory IgA and serum albumin in the tear film of people with keratoconus. Exp Eye Res (2012) 96(1):132–7. 10.1016/j.exer.2011.12.010
34.
Yenihayat F Altıntaş Ö Kasap M Akpınar G Güzel N Çelik OS . Comparative proteome analysis of the tear samples in patients with low-grade keratoconus. Int Ophthalmol (2018) 38(5):1895–905. 10.1007/s10792-017-0672-6
35.
López-López M Regueiro U Bravo SB Chantada-Vázquez M del P Varela-Fernández R et al Tear proteomics in keratoconus: a quantitative SWATH-MS analysis. Invest Ophthalmol Vis Sci (2021) 62(10):30. 10.1167/iovs.62.10.30
36.
Tong L Lan W Lim RR Chaurasia SS . S100A proteins as molecular targets in the ocular surface inflammatory diseases. Ocul Surf (2014) 12:23–31. 10.1016/j.jtos.2013.10.001
37.
Zhang Y Liu F . Elevation of S100 calcium-binding protein A7 in recurrent pterygium. Exp Ther Med (2019) 18:3147–52. 10.3892/etm.2019.7922
38.
Gupta V Kawahara G Gundry SR Chen AT Lencer WI Zhou Y et al The zebrafish dag1 mutant: a novel genetic model for dystroglycanopathies. Hum Mol Genet (2011) 20(9):1712–25. 10.1093/hmg/ddr047
39.
André da Silva R Moraes de Paiva Roda V Philipe de Souza Ferreira L Oliani SM Paula Girol A Gil CD Annexins as potential targets in ocular diseases, 27. Drug Discovery Today. Elsevier Ltd (2022).
40.
Watanabe M Kondo S Mizuno K Yano W Nakao H Hattori Y et al Promotion of corneal epithelial wound healing in vitro and in vivo by annexin A5. Invest Ophthalmol Vis Sci (2006) 47(5):1862–8. 10.1167/iovs.05-0591
41.
Zhou L Beuerman RW . Tear analysis in ocular surface diseases. Prog Retin Eye Res (2012) 31(6):527–50. 10.1016/j.preteyeres.2012.06.002
42.
Lupasco T He Z Cassagne M Sagnial T Brion L Fournié P et al Corneal epithelium in keratoconus underexpresses active NRF2 and a subset of oxidative stress-related genes. PLoS One (2022) 17(10 October):e0273807. 10.1371/journal.pone.0273807
43.
Acera A Vecino E Rodríguez-Agirretxe I Aloria K Arizmendi JM Morales C et al Changes in tear protein profile in keratoconus disease. Eye (2011) 25(9):1225–33. 10.1038/eye.2011.105
44.
Recalde JI Duran JA Rodriguez-Agirretxe I Soria J Sanchez-Tena MA Pereiro X et al Changes in tear biomarker levels in keratoconus after corneal collagen crosslinking. Mol Vis (2019) 25:12–21.
45.
Balasubramanian SA Wasinger VC Pye DC Willcox MDP . Preliminary identification of differentially expressed tear proteins in keratoconus. Mol Vis (2013) 19:2124–34.
46.
Kontoh-Twumasi R Budkin S Edupuganti N Vashishtha A Sharma S . Role of serine protease inhibitors A1 and A3 in ocular pathologies. Invest Ophthalmol Vis Sci (2024) 65(2):16. 10.1167/iovs.65.2.16
47.
Jackson BC Thompson DC Wright MW Mcandrews M Bernard A Nebert DW et al Update of the human secretoglobin (SCGB) gene superfamily and an example of “evolutionary bloom” of androgen-binding protein genes within the mouse Scgb gene superfamily. Hum Genomics (2011) 5(6):691–702. 10.1186/1479-7364-5-6-691
48.
Soria J Villarrubia A Merayo-Lloves J Elortza F Azkargorta M Alvarez de Toledo J et al Label-free LC-MS/MS quantitative analysis of aqueous humor from keratoconic and normal eyes. Mol Vis (2015) 21:451–60.
49.
López-López M Regueiro U Bravo SB Chantada-Vázquez M del P Pena C et al Shotgun proteomics for the identification and profiling of the tear proteome of keratoconus patients. Invest Ophthalmol Vis Sci (2022) 63(5):12. 10.1167/iovs.63.5.12
50.
Leonardi A Donato A Rosani U Di Stefano A Cavarzeran F Brun P . Endoplasmic reticulum stress and unfolded protein response in vernal Keratoconjunctivitis. Invest Ophthalmol Vis Sci (2024) 65(4):23. 10.1167/iovs.65.4.23
51.
Jaskiewicz K Maleszka-Kurpiel M Michalski A Ploski R Rydzanicz M Gajecka M . Non-allergic eye rubbing is a major behavioral risk factor for keratoconus. PLoS One (2023) 18(4 APRIL):1–21. 10.1371/journal.pone.0284454
52.
Perumal N Funke S Pfeiffer N Grus FH . Proteomics analysis of human tears from aqueous-deficient and evaporative dry eye patients. Sci Rep (2016) 6(March):1–12. 10.1038/srep29629
Summary
Keywords
biomarkers, cornea, keratoconus, proteomics, tear film
Citation
de Almeida Borges D, Alborghetti MR, Domingues RR, Paes Leme AF and Alves M (2026) Tear proteomic analysis in keratoconus patients and potential biomarkers: a case-control study. Exp. Biol. Med. 250:10864. doi: 10.3389/ebm.2025.10864
Received
15 October 2025
Revised
10 December 2025
Accepted
22 December 2025
Published
12 January 2026
Volume
250 - 2026
Updates
Copyright
© 2026 de Almeida Borges, Alborghetti, Domingues, Paes Leme and Alves.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Daniel de Almeida Borges, dborges@outlook.com.br
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.