PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 17478.
Published online 2017 December 12. doi:  10.1038/s41598-017-17536-2
PMCID: PMC5727318

Tear proteome analysis in ocular surface diseases using label-free LC-MS/MS and multiplexed-microarray biomarker validation

Abstract

We analyzed the tear film proteome of patients with dry eye (DE), meibomian gland dysfunction (MGD), and normal volunteers (CT). Tear samples were collected from 70 individuals. Of these, 37 samples were analyzed using spectral-counting-based LC-MS/MS label-free quantitation, and 33 samples were evaluated in the validation of candidate biomarkers employing customized antibody microarray assays. Comparative analysis of tear protein profiles revealed differences in the expression levels of 26 proteins, including protein S100A6, annexin A1, cystatin-S, thioredoxin, phospholipase A2, antileukoproteinase, and lactoperoxidase. Antibody microarray validation of CST4, S100A6, and MMP9 confirmed the accuracy of previously reported ELISA assays, with an area under ROC curve (AUC) of 87.5%. Clinical endpoint analysis showed a good correlation between biomarker concentrations and clinical parameters. In conclusion, different sets of proteins differentiate between the groups. Apolipoprotein D, S100A6, S100A8, and ceruloplasmin discriminate best between the DE and CT groups. The differences between antileukoproteinase, phospholipase A2, and lactoperoxidase levels allow the distinction between MGD and DE, and the changes in the levels of annexin A1, clusterin, and alpha-1-acid glycoprotein 1, between MGD and CT groups. The functional network analysis revealed the main biological processes that should be examined to identify new candidate biomarkers and therapeutic targets.

Introduction

DE and MGD diseases involve the disruption of the lacrimal functional unit, resulting in symptoms of discomfort and visual disturbance and tear film instability. These illnesses may exist independently as either symptomatic or asymptomatic disorder; however, they are frequently found in the same patient. The progression of these diseases typically leads to alterations in the tear film and tear hyperosmolarity1 and the secretion of inflammatory mediators into the tears, initiated by cytokine release and metalloproteinase activation2. Ocular dryness disorders also promote squamous metaplasia (SM), mainly as a consequence of inflammation3, whose severity can be used for grading different diseases affecting the ocular surface4,5.

Tear film proteomics offers powerful analytical tools for studying the proteins involved in ocular diseases. Using these tools, we might be able to establish the precise functions of these proteins in the underlying pathophysiological processes and provide diagnostics biomarkers. The studies of differential protein expression in complex biofluids such as tear film require rapid, highly reproducible, and accurate quantification. These requirements can be satisfied by mass spectrometry-based analysis6. Label-based mass spectrometry techniques have been employed to analyze the tear film, mainly by using isobaric tags for relative and absolute quantitation (iTRAQ)710. However, label-free shotgun proteomics techniques are becoming increasingly popular since they are faster and supply cleaner and simpler results than the label-based techniques6. One of the most used label-free methods is spectral counting of identified proteins. This method is based on the fact that the number of spectra corresponding to peptides originating from a given protein shows a good linear correlation with the abundance of that protein1114. Nevertheless, count ratios can be biased because of particular physiochemical properties of peptides or sampling issues in Data-Dependent Acquisition (DDA) shotgun procedure that may affect their detection by MS. As not all peptides can be perfectly ionized, appropriate measurement corrections should be applied to obtain accurate data15,16.

To date, several proteomic techniques have been applied to study variations in tear proteome caused by pathological conditions and/or medical treatments. Many MS-based proteomic studies of tear film have been carried out using iTRAQ710, label-free14,1719, SELDI-TOF20, or MALDI-TOF2123 analysis. However, the results have been, at times, conflicting. Some of the studies have demonstrated changes in the expression of several proteins. Overexpression has been reported for protein ANXA114,2325, S100A97,14,24,25, S100A87,14,20,2325, and S100A47,24 and downregulation, for prolactin-inducible protein (PIP)7,8,2426, lactotransferrin (LT)7,8,14,24,27, proline-rich protein 4 (PRP4)8,14,20,25,26, lysozyme (LYZ)7,8,14,20,25,27, and cystatin-S (CST4)8,14,24,25, among others. Nevertheless, the changes in the expression of some proteins remain controversial. For example, several studies20,27 have demonstrated upregulation of lipocalin-1 (LCN1), but its downregulation has been shown in others7,8,24,28. Similarly, conflicting results have been reported for other proteins such as secretoglobin 2A28,26,28 and serum albumin14,20,2628. These discrepancies might have been caused by several factors such as different sample collection methods or sensitivity and dynamic range of the employed quantitation techniques29,30. Therefore, more extensive tear proteome studies focusing on systems biology, supported by orthogonal validation assays, are needed to extract meta-information common to all these studies. The results of such studies should reduce the technological and experimental bias and improve our knowledge of the role of relevant proteins in the pathophysiology of ocular surface-related disorders.

We have previously characterized tear film proteomes and identified biomarkers for the aqueous-deficient dry eye (ADDE) and MGD, the main cause of evaporative DE (EDE) condition31, using 2D gel-based proteomics24.

The purpose of the present study was to identify proteins that can be used for discrimination between DE and MGD pathologies. To achieve this, we performed a quantitative differential study of tear protein expression using LC-MS/MS spectral-counting quantitative proteomics. Three candidate biomarkers were validated employing an orthogonal technique (customized microarray assays) using independent DE and CT tear samples.

Materials and Methods

Patients

A retrospective case-controlled study was carried out, in which 70 patients were enrolled. The experimental design consisted of a discovery phase (including 7 DE patients, 12 chronic MGD patients, and 18 CT individuals for biomarker candidate selection), and a validation phase, for which 24 DE patients and 9 CT individuals were recruited (Fig. 1).

Figure 1
Workflow overview. A general workflow outlining the steps used in the study for the identification of proteins deregulated in DE and MGD tear samples, using LC-MS/MS label-free quantitative proteomics, candidate biomarkers validation, and clinical correlation ...

This research was conducted by medically qualified personnel after approval by the Cruces Hospital and the Principado de Asturias Hospital Ethics Committees. Approval was obtained in strict accordance with the tenets of the Declaration of Helsinki on Biomedical Research Involving Human Subjects. Patients were recruited at the Instituto Clínico Quirúrgico de Oftalmología (Bilbao, Bizkaia, Spain) and the Instituto Universitario Fernández-Vega (Oviedo, Asturias, Spain). Before tear collection, the signed informed consent was obtained from all patients, once the nature and possible consequences of the study had been explained. The diagnosis was based on clinical examination including the Schirmer I test with anesthesia to measure the basal secretion, slit-lamp examination of the lid margin and Meibomian glands, fluorescein staining results according to Oxford scale, and subjective symptoms. Each patient answered the OSDI questionnaire, which included some statements about the quality of their vision and wellbeing.

Patients were classified as having DE if they had dry eye symptoms, dynamics abnormalities in Schirmer I test (≤5 mm/5 min), characteristic fluorescein staining of the corneal epithelium, and a tear breakup time (TBUT) ≥5 s. MGD was diagnosed if the patients had eyelid inflammation, Schirmer I test results >5 mm/5 min, TBUT <5 s and alterations in Meibomian glands. Healthy subjects in the CT group were not suffering from any ocular disease (no allergic or atopic history). They presented Schirmer I test values >5 mm/5 min, TBUT >5 s, no corneal fluorescein staining or sensations of discomfort, and no evident eyelid inflammation. The individuals with (or history of) any systemic or ocular disorder or condition (including ocular surgery, trauma, and disease) and patients with Sjögren’s syndrome were excluded. Contact lens users were also excluded to avoid any possible interference with the interpretation of the results.

Clinical evaluation of patients and sample collection

The order of tests and sample collection was always the same. First, the OSDI questionnaire was conducted to assess the symptoms of ocular irritation and their effect on the vision-related function. Second, one drop of topical anesthetic was applied to the ocular surface and the Schirmer I and fluorescein tests were performed. The tear samples and conjunctival impression cytology samples were collected one day later to avoid the interference between clinical tests and the proteomic study.

All tear samples were collected using calibrated 10-µl glass microcapillary tubes (BLAUBRAND intraMark, Wertheim, Germany). Tear samples were obtained from the inferior temporal tear meniscus, minimizing irritation of the ocular surface or lid margin, without anesthesia. After collection, the tear samples were placed in Eppendorf tubes and stored at −80 °C until analysis. Protein quantification was performed using the EZQ Protein Quantification Kit (Invitrogen Dynal AS, Oslo, Norway).

To evaluate the degree of damage of the conjunctival epithelium (by examining the extent of SM), conjunctival impression cytology (CIC) samples were obtained for PAS-hematoxylin staining. To achieve this, 5 × 5-mm membranes of cellulose acetate (HAWP304, Millipore, Bedford, MA, USA) were applied to the bulbar superior conjunctival epithelium following topical anesthesia (Colircusi double anesthetic, Alcon Cusí, Barcelona, Spain). The samples were immediately fixed in 96% ethanol and stained with PAS-hematoxylin, as previously described32. The samples were later examined under a light microscope to evaluate the grade of SM33. The cytoplasmic and nuclear areas of non-secretory cells, cytoplasmic alterations and staining, nuclear alterations, the nuclear to cytoplasmic area ratio (N:C ratio), and the number of goblet cells were examined.

Tear protein digestion, LC-MS/MS, and protein identification

To digest the proteins, RapiGest SF (Waters, Milford, MA, USA) was added to 4 µg of total tear protein at a concentration of 0.25% (w/v), and the mixture was heated at 85 °C for 15 min with agitation. Proteins were reduced (5 mM dithiothreitol, 60 °C, 25 min), alkylated (15 mM iodoacetamide, room temperature, 30 min), and digested with trypsin (0.2 µg, 37 °C, overnight; Roche Diagnostics, Penzberg, Germany). Finally, RapiGest SF was inactivated following the manufacturer’s instructions.

LC-MS/MS analysis was performed in a nanoACQUITY UPLC system interfaced to a SYNAPT HDMS mass spectrometer (Waters Corporation, Milford, MA). An aliquot of 1 µg of protein from each sample was loaded onto a Symmetry 300 C18, 180 µm × 20 mm precolumn (Waters Corporation, Milford, MA, USA) and washed with 0.1% formic acid for 3 min at a flow rate of 5 µl/min. The precolumn was connected to a BEH130 C18, 75 µm × 200 mm, 1.7 µm (Waters Corporation, Milford, MA, USA), equilibrated in 3% acetonitrile with 0.1% formic acid. Peptides were directly eluted with a 120-min linear gradient of 3–60% acetonitrile onto a homemade nanoESI emitter. Data-dependent MS/MS acquisitions were performed on the 3 most intense precursors per scan, with charge states of 2, 3, or 4 over a survey m/z range of 400–1990 and a dynamic peak exclusion of 30 s. Collision energies were varied as a function of the m/z and charge state of each peptide. [Glu1]-fibrinopeptide B (Sigma-Aldrich, St. Louis, MO, USA) at a concentration of 100 fmol/ml was sprayed through the NanoLockSpray source and sampled every 30 s. The obtained spectra were processed using VEMS34 and searched using MASCOT version 2.2.03 (Matrix Science, London, UK) against the UniProtKB/Swiss-Prot (version 2015_12) database, using only Homo sapiens entries. For protein identification, the following parameters were adopted: carbamidomethylation of cysteines as fixed modification, oxidation of methionines as variable modification, peptide mass tolerance of 50 ppm, fragment mass tolerance of 0.1 Da, and 1 missed cleavage.

Quantitative protein analysis using spectral counting (APEX quantitation)

Protein quantitation was performed using label-free spectral counting, employing APEX Quantitative Proteomics Tool v.1.1 as described previously35. MASCOT output files were converted to the pepXML format using Trans-Proteomic Pipeline software v. 3.2.2 (Institute for Systems Biology, Seattle, WA)36. The resulting pepXML files were analyzed using PeptideProphet37. Subsequently, ProteinProphet38 was employed to assemble peptides into protein identification groups. Only proteins with ProteinProphet probability  0.95 and at least two identified peptides with a PeptideProphet probability  0.95 were reported. A probability-based penalty for peptide detection was applied, given that some peptides are detected by MS more readily than others. APEX quantitation analysis involves three steps. In step 1, a classification model is built, based on physiochemical properties of the peptide sequence (peptide mass, length, amino acid composition, properties related to charge, hydrophobicity, and amino acid frequencies within the secondary peptide structures). The model is used to estimate MS-detectability (O(i)) for any given protein. Step 2 uses spectral count information to obtain absolute abundances of proteins in each MS/MS experiment, corrected by applying an MS-detectability correction factor. Step 3 allows statistical significance analysis of differential expression in distinct biological samples. The protein false discovery rate was set at 1%; only proteins identified at a 99% confidence level were used for spectral counting.

Statistics

Data matrix obtained from APEX quantitation was used for multivariate statistical analysis. Before statistical analysis, protein expression data were filtered by considering only the proteins identified in at least 50% of biological samples. K-nearest neighbor data imputation was performed to complete the data matrix and manage missing values imposed by intrinsic sampling issues in DDA shotgun approaches. Data were then normalized using total spectral count normalization to reduce run-to-run variation39. Once the data matrix had been prepared for analysis, feature subset subtraction was performed by stepwise discriminant analysis with a p-value cutoff of 0.01, to obtain the most significant attributes defining each group. This process extracts the information that contributes most to the discrimination between groups. To reduce the dimensionality of the dataset for 2D-visualization including the possible overlap between the groups, sample clustering was performed (using an exploratory data technique, the canonical discriminant analysis). Normalization, imputation, feature subset, and clustering analysis were carried out using R-Package v2.15.240.

Functional network analysis

Once the most significant tear proteins had been identified by LC-MS/MS analysis, the list of the most significantly upregulated or downregulated proteins was loaded into Cytoscape v. 2.8.341 using Reactome Functional Interaction Network Cytoscape plugin42, to find network patterns related to diseases. The Cytoscape plugin accesses the Reactome Functional Interaction Network database. This database extends curated pathways with non-curated sources of information, including protein–protein interactions, gene coexpression, protein domain interaction, Gene Ontology (GO) annotations, and text-mined protein interactions. The process obtains a biological network based on the experimentally found proteins and the directly related proteins required to connect them in a network. With this information, functional interaction sub-networks (modules) were constructed. Then, the analysis of network modules of highly interacting groups of proteins and the functional enrichment and pathway analysis were performed. To select the biological processes and/or pathways that are predominant in each of the modules, we filtered the functional enrichment results using a false discovery rate cutoff value of 0.01. Finally, the Hubba Cytoscape Plugin43 was applied to the functional interaction subnetwork to find the most important nodes/hubs within the network. Topological analyses were performed using the Edge Percolated Component ranking method to evaluate node essentiality.

Validation of candidate biomarkers using antibody microarray assays

Three proteins with significant associations in mass spectrometry and network analysis (S100A6, CST4, and MMP9) were evaluated using customized sandwich-like antibody microarrays. The process of customized microarray preparation for quantitative tear biomarker analysis included several steps. These were i) identification of compatible antibody pairs for each biomarker and selection of high-affinity pairs by surface plasmon resonance (SPR) and sandwich ELISA assays, ii) evaluation of cross-reactivity between components in the multiplexed immunoassay, iii) determination of detection limits for each biomarker, iv) integration into the microarray platform and technical validation, and v) measurement of biomarker concentration in the selected tear samples (for validation purposes).

Briefly, mouse monoclonal antibodies against the selected proteins were spotted on 3D-epoxy-activated glass slides and cyclo-olefin polymer slides (PolyAn GmbH, Berlin, Germany) with the non-contact spotter sciFLEXARRAYER S11 (Scienion AG, Berlin, Germany), in a 12-subarray format. Each antibody was spotted 20 times to reduce the technical variability. The nozzle used was PDC90, and a voltage of 82 V was applied, with a pulse width of 46 µs, loading a final volume of 1.5 nl/spot. The arrays were blocked with TBS-T supplemented with BSA (1%) for 1 hour. Antigens S100A6, CST4, and MMP9 were incubated for 1 h to obtain a standard curve from 0.195 to 166 ng/ml. The reaction volume was 70 µl/well for all the steps in the immunoassay. Tear samples were diluted (1/30) in PBS for microarray analysis. Subsequently, the samples were incubated for 1 h with rabbit detection antibodies. Finally, after washing the slides with TBS-T, the secondary Alexa Fluor 647-labeled anti-rabbit antibodies were added and incubated for 1 h. Fluorescence of the spots was measured using a Tecan LS Reloaded Microarray Scanner (Tecan Deutschland GmbH, Crailsheim, Germany) at 633 nm, and protein concentration was determined based on standard curve intensity values.

After testing for normal distribution, significant differences between the groups were determined using Kruskal–Wallis non-parametric test. Once the concentration of each protein had been obtained, ROC curve analysis was performed, using logistic regression. The classification algorithm was trained using the biomarker panel, and to assess learner performance, a 10-fold random sampling was carried out with 70% of the samples as a training set. The remaining 30% of the samples constituted a test set. Statistical analysis was carried out using the Orange Canvas statistics package, version 2.6 (http://orange.biolab.si/). Finally, the Pearson correlation analysis was performed to assess correlations between clinical parameters and protein levels. The statistical analysis was conducted using R statistical program40.

Data availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Results

Patients

The patients recruited for the study were grouped into two sub-study groups, for discovery and validation. In the first group (37 volunteers), tear proteomics and CIC analyses were conducted. Demographical and clinical data for the discovery phase are presented in Table 1.

Table 1
Demographic and clinical data of patients included in the discovery phase of the study.

We followed the recommendations for tests appropriate for diagnosis and assessment of the target ocular surface diseases studied here, DE and MGD44. Various diagnostic tests were used, depending on the objective. Thus, in the first part of the study (biomarker candidate selection), we analyzed patients with DE and MGD. We used the symptoms as a screening tool (OSDI), tear volume assessment to ascertain the subtype classification of DE (Schirmer test), and damage to the ocular surface (impression cytology) as a standard test to study squamous metaplasia and goblet cell density of the conjunctiva. In the second part of the study (the validation of candidate biomarkers), we assessed the DE patients and CT subjects; we employed the recommended diagnostic tests for DE44. Following the clinical protocol proposed for DE diagnosis, a broad range of diagnostic tests was performed, including OSDI (ocular surface symptoms), TBUT (tear film stability), Schirmer I test (tear volume), and fluorescein staining (damage to the ocular surface). This was also done to determine the correlation matrix for biomarker concentrations and clinical parameters.

PAS-hematoxylin impression cytology revealed that in control samples, 72% of individuals (13/18) presented grade 0 SM (G0), whereas the remaining 28% (5/18) showed grades 1–2 (Table 1). In the MGD group, 17% of patients (2/12) had SM of grade 0, whereas mild grades (G1 and G2) were seen in 50% (6/12), and moderate (G3) to severe (G4) grades in 33% of patients (4/12). In the DE group, a mild grade of SM was found in 57% of patients (4/7), while the remaining 43% of DE patients exhibited moderate (G3) and severe (G4) grades. Schirmer test results showed significant differences (p-value < 0.05) between the groups, with values of 1.7, 8.0, and 13.8 for the DE, MGD, and CT groups, respectively. Mean age differences between the groups were not statistically significant (p-value > 0.05).

In the second part of the study, tear samples from 33 patients were analyzed on microarrays for validation purposes. To evaluate the correlation between biomarkers and clinical parameters, OSDI, TBUT, Schirmer I test, and fluorescein staining results were obtained for each patient in this phase of the study. All these parameters differed significantly between DE and CT individuals. The mean values for each parameter are presented in Table 2.

Table 2
Demographic and clinical data of patients included in the validation phase of the study.

Quantitative protein analysis using spectral counting

Quantitative LC-MS/MS spectra analysis of tear proteome was performed computationally using the adjusted spectral counting method.

The analysis identified 603 significantly distinct proteins, by filtering the data to a minimum of two peptides identified with a PeptideProphet probability of 0.9 and a ProteinProphet probability of 0.95. However, as only the proteins identified in at least 50% of the samples were considered, this number was reduced to 135 proteins.

Protein expression data for those 135 proteins were used in further analysis. Model-based data imputation was performed to complete the data matrix and manage missing values. To establish which proteins were characteristic for each pathological condition, a subtraction of expression data, using stepwise discriminant analysis modeling, was performed. As a result, we obtained a list of 26 specific proteins that undergo expression changes, achieving an effective discrimination between the groups (Table 3).

Table 3
Proteins identified by LC-MS/MS with the most significant changes in their expression levels in dry eye (DE), Meibomian gland dysfunction (MGD), and control (CT) groups.

The clearest discrimination between the groups could be obtained by examining the expression of the following proteins (overexpressed in the DE group): complement C3 (C3), protein S100-A6 (S100A6), protein S100-A8 (S100A8), ceruloplasmin (CP), apolipoprotein D (APOD), alpha-1-acid glycoprotein 2 (ORM2), thioredoxin (TXN), Ig gamma-1 (IGHG1), membrane-associated phospholipase A2, (PLA2G2A), alpha-1-antitrypsin (SERPINA1), and antileukoproteinase (SLPI). Of these proteins, TXN, IGHG1, PLA2G2A, SERPINA1, and SLPI showed the lowest levels of expression in the MGD group, while the remaining proteins had the lowest expression in the control group. Similarly, in the MGD group, we found overexpression of the following proteins: annexin-A1 (ANXA1), clusterin (CLU), alpha-1-acid glycoprotein 1 (ORM1), and lactoperoxidase (LPO). LPO was the most downregulated protein in the DE group. Finally, this analysis also revealed, both in the DE and MGD groups, reduced expression of several proteins. These were proline-rich protein 1 (PROL1), secretoglobin family 1D member 1 (SCGB1D1), mammaglobin-B (SCGB2A1), lipocalin-1 (LCN1), cystatin-S (CST4), prolactin-inducible protein (PIP), dermcidin (DCD), submaxillary gland androgen-regulated protein 3B (PROL3), galectin-3-binding protein (LGALS3BP), extracellular glycoprotein lacritin (LACRT), and proline-rich protein 4 (PRR4). The lowest levels of expression of these proteins were found in the DE group. For each protein, Table 3 lists fold-changes and the groups with its highest and lowest expression values.

Using data mining techniques, we also evaluated the spatial separation that these proteins are capable of producing in the three tested groups. We used the expression data of the 26 proteins in the canonical discriminant analysis. We found a good separation between the experimental groups, reflecting the discriminative power presented by this list of proteins. However, it should be noted that there was a larger spatial separation between the DE and control group than between MGD and control (Fig. 2).

Figure 2
Canonical Discriminant Analysis results showing the separation between the samples as a function of the values for the 26 discriminant proteins obtained by stepwise discriminant analysis. Each of the points represents a sample from each group. Good separation ...

Functional interaction networks

The Cytoscape Reactome Functional Interaction Network plugin was used to find network patterns by accessing the Reactome Functional Interaction Network database. The disease-related network consisted of 19 of the candidate proteins associated with the diseases, found in this study, and the minimum number of interacting proteins necessary to interconnect them in a network (Fig. 3). Once the network had been constructed, functional clustering was performed using the same Cytoscape plugin. Cluster construction was based on proximity between the nodes. Five topologically well-differentiated modules were obtained (modules 1 to 5), comprising 11, 11, 15, 4, and 4 proteins, respectively. Functional enrichment analyses were then performed for each module to determine the biological processes involved in the pathologies, in relation to the structure of the network. With a false discovery rate cutoff of 0.01, we found statistically significant biological processes in 4 of the 5 modules. As illustrated in Table 4, module 1 of the network mainly represents the processes such as inflammatory response, regulation of chemokine production, response to lipopolysaccharide, cytokine-mediated signaling pathway, apoptotic signaling, acute-phase response, cell redox homeostasis, innate and humoral immune response, regulation of JNK cascade, and cellular defense response. Module 2 was principally formed by proteins involved in the extracellular matrix organization, collagen catabolic process, response to cytokine stimulus, angiogenesis, aging, response to hypoxia, regulation of cell proliferation, and neutrophil chemotaxis. Module 3 was made up of proteins significantly related to the defense response to bacteria and positive regulation of the inflammatory response. Module 4 was related to signal transduction, regulation of apoptosis, and fibroblast proliferation processes. Module 5 did not show significant relationships with any biological processes.

Figure 3
The functional interaction network obtained using the Functional Interaction Reactome Cytoscape plugin41. Proteins in the interaction network are represented as nodes (red dashed-line circles), while the interaction between any two proteins is represented ...
Table 4
Involvement of each module in different biological processes, showing the proteins implicated in each process42.

Finally, to find the most significant nodes in the network, the BisoGenet plugin was applied to the network. This plugin generates numerical and graphical output, making it easy to find the key nodes in the complex networks. This analysis can provide biologically meaningful node identification and functional classification by identifying the major nodes, the so-called “hub nodes,” which have many incoming and/or outgoing connections within a network. As a result, we found that the proteins MMP9, JUN, RELA, STAT3, TLR4, ESR1, GNAI3, ANXA1, and TIMP1 were the hub nodes of the functional interaction network. These proteins are involved in inflammation, adherens junction process, signaling pathways, cell–cell communication, and signal transduction.

Validation of candidate biomarkers using microarray assays

To validate further the relative protein expression data of tear proteins biomarkers obtained via MS-based spectral-counting techniques, we employed an orthogonal technique based on antibody microarrays, using an independent group of DE and control samples. The S100A6 and CST4 proteins, significantly deregulated according to the results of the LC-MS/MS analysis, were selected for validation. The MMP9 protein, inferred from network analysis, was also used in this validation. The antibody microarray assays confirmed significant changes in CST4 and S100A6 levels seen in LC-MS/MS assays. They also corroborated the network analysis inference of the implication of MMP9 in the DE disease (Fig. 4). The results showed that the tear concentrations of S100A6 (p-value = 0.0008) and MMP9 (p-value = 0.0097) were significantly higher in the DE group in comparison with the control group, confirming the mass spectrometry findings. However, the CST4 level was significantly lower in the DE group than in the CT group (p-value = 0.0028) (Fig. 5). The ratios between the DE and CT group for S100A6, MMP9, and CST4 were 2.67, 2.48, and −1.59, respectively.

Figure 4
Antibody microarrays. (A) Spotting pattern for the 12-subarray format. (B) Representative image of the arrays showing the distribution of the standard calibration curve and samples. Each microarray slide contains a standard calibration curve (left column) ...
Figure 5
Comparison of the concentrations of candidate biomarkers in control (CT) and dry eye (DE) groups, measured using the customized antibody microarrays. Concentration is expressed in ng/mL. Green circles represent the mean concentration for the group. ( ...

In addition, we used these three markers to generate multivariate predictive models. To this end, we performed training employing random sampling with 70% of the data and a logistic regression algorithm. Subsequently, the classifier method was validated using the remaining 30% of data. The results showed an AUC value of 87.5% (Fig. 5D).

Correlation of validated tear biomarkers with clinical endpoints

The data obtained for proteins S100A6, CST4, and MMP9 were used to assess the correlation between the biomarkers and clinical parameters. A post-hoc Pearson correlation analysis was performed to determine whether the individual protein levels were associated with the tested clinical parameters (Schirmer I test, TBUT, OSDI, fluorescein staining, and age).

S100A6 and MMP9 showed significant positive correlation with OSDI and fluorescein staining and negative correlation with TBUT and Schirmer I test (Table 5). Similarly, S100A6 and MMP9 biomarkers positively correlated with each other. CST4 protein negatively correlated with fluorescein staining. There was no correlation between the biomarkers and the age of patients.

Table 5
Pearson correlation matrix for biomarker concentrations and clinical parameters.

Discussion

In our MS-based proteomic study, we first compared tear samples of DE and MGD patients with control tear samples. The following proteins were mostly overexpressed in DE: C3, S100A6, S100A8, CP, APOD, ORM2, TXN, IGHG1, PLA2G2A, SERPINA1, and SLPI. These proteins are involved in biological processes related to the defense response, inflammatory process, and response to wounding. Among these proteins, TXN, IGHG1, PLA2G2A, SERPINA1, and SLPI were most downregulated in MGD.

The most representative overexpressed proteins in the MGD group were ANXA1, CLU, ORM1, and LPO; the LPO protein had its lowest expression level in the DE group. These proteins are implicated mainly in apoptosis, oxidative stress, immune response, and keratinocyte differentiation, which represent some of the principal processes involved in ocular surface diseases45. A panel of TXN, IGHG1, PLA2G2A, SERPINA1, SLPI, and LPO proteins could unequivocally discriminate between DE and MGD as, according to our observations, their regulatory trends in these disease groups are opposite.

TXN is implicated in immune response, cell proliferation, cell–cell signaling, and oxidation–reduction processes46. IGHG1 is associated with immune response47, and PLA2G2A is involved in the defense and inflammatory responses and lipid catabolic processes48. SERPINA1 is involved in acute-phase response, proteolysis and platelet activation49, and SLPI, in the defense response50. The increased expression of these proteins in DE and the reduction of their levels in MGD suggest an enhanced activation of defense response, inflammatory, proteolytic, or cell–cell signaling processes triggered by a decrease in tear production.

In contrast, LPO is associated with the response to oxidative stress51. The direct relationship between hyperosmolarity and oxidative stress has been described by Zhang et al.52. Gaffney et al.53 have predicted substantially higher levels of hyperosmolarity in evaporative DE than in aqueous-deficient DE condition, which might explain the high expression levels of LPO protein in MGD group in comparison with DE. This observation and the spatial separation of MGD and DE groups in the canonical discriminant analysis, seen here suggest a different activation of osmolarity compensatory mechanisms, depending on the typology of dry eye.

In summary, the changes in regulatory processes in which these proteins are involved (such as endopeptidase activity, response to bacteria, acute-phase response, cell–cell signaling, or response to oxidative stress) might be the key biological events which distinguish DE from MGD. Validation of those candidate protein biomarkers in new tear samples will be a target of further studies.

In contrast, the expression of another group of proteins was significantly downregulated in both pathologies: PROL1, SCGB1D1, SCGB2A1, LCN1, CST4, PIP, DCD, PROL3, LGALS3BP, LACRT, and PRR4. Most of these are proteins secreted by the lacrimal and Meibomian glands54,55. Taking into account the main functions of these proteins, this result suggests an association of the secretory deficiency of aqueous-deficient DE or obstructive MGD with a decreased antibacterial activity, which could increase a risk of ocular infection. Some of these proteins, such as CST4, LCN1, SCGB2A1, and PIP, have also been downregulated in a previous study reported by our group24. Clustering analysis revealed that this protein group, characterized by a reduction in their expression levels, could be used to discriminate between DE and MGD.

We also found a very significant separation between the control and DE groups, with MGD occupying an intermediate position between these two groups (as presented in Fig. 2). Interestingly, a similar behavior was observed by our research group in a recent study56, in which we compared the conjunctival epithelium proteomes of DE, MGD, and CT individuals. In that work, we observed, based on a group of 17 proteins, a clear separation between DE and CT groups, with MGD in an intermediate position, closest to the CT group. These independent results from tear and epithelium analysis support the hypothesis that tear proteome alterations accompany the changes in conjunctival epithelium.

Apart from the unclear etiology of the studied pathologies, the heterogeneity of the reported results constitutes another key challenge in proteomic studies of complex diseases. Several factors can affect the results obtained for the tear proteome, including the tear collection method and sample storage29, analytical methods30, and differences in the patient inclusion and exclusion criteria. To validate a robust biomarker, it is of paramount importance to establish its reproducibility using several different technological approaches. To achieve this, we compared our results with those obtained in other proteomic tear studies. Tear samples used in these studies, representing the same pathologies, have been analyzed using 2D-PAGE24,28,57, SELDI-TOF profiling20, MS-based quantification approaches79,14,17,19,23,25,26,28,5860, or ELISA assays6163. The comparison revealed discrepancies between the results of different studies, probably caused by some of the previously mentioned factors. However, the proteins whose expression behavior is similar in different studies acquire greater weight in terms of their robustness as biomarkers. Thus, despite the differences in the analytical technologies employed, we found a group of proteins whose expression was consistently upregulated in the DE group, in comparison with the controls, in the cited studies (S100A8, S100A6, C3, CP, ORM2, IGHG1, SERPINA1, ANXA1, CLU, and ORM179,14,20,23,24,26,28,5863). Similarly, we found another group of proteins with reduced expression levels in the pathological groups, including PIP, proline-rich family (PROL1, PROL3, PRR4), SCGB1D1, CST4, LACRT, and LCN179,14,17,20,2426,28,57,60. The proteins whose expression behavior is consistent in different studies are listed in Table 6.

Table 6
Proteins with expression profiles similar to those in other studies.

The proteins S100A6, ANXA1, and CST4 are particularly noteworthy as molecular markers. Apart from their consistent expression behavior in various studies, their validity has been already confirmed using absolute quantification in ELISA assays with independent samples24. That validation study has been conducted using 44 individuals for a biomarker discovery pipeline and 100 for validation. The study has shown a correct assignment of 97.9% in sample classification for a panel of five markers (S100A6, ANXA1, CST4, PLAA, and ANXA11), with a precision of 93.5% for the first three markers. These results confirm the validity of these biomarkers.

Our study also found novel markers associated with ocular surface disorders: APOD, TXN, PLA2G2A, SLPI, LPO, LGALS3BP, and DCD.

Label-free approaches are increasingly considered advantageous and reliable since no additional chemical manipulation is needed. These proteomic methods have gained popularity due to their compatibility with high-throughput systems, their speed, and good reproducibility in the complex peptide mixture analyses6. Spectral-counting-based quantification is more reproducible and has a larger dynamic range than other label-free quantification methods64. However, in spectral-counting proteomics, physicochemical properties of the peptides can affect MS detection. This might result in imprecise results due to errors in the number of repeat observations of MS/MS spectra of the peptides. Therefore, we applied a peptide-specific normalization factor, which depends on peptide detection probability and improves the accuracy and reproducibility of the analysis65,66.

Complex diseases constitute intricate systems of altered biological processes. The study of systems biology using functional interaction networks reveals the changes triggered by various stressors (evaporation, osmolarity, hyposecretion) and the molecular machinery that produces the diseases. Clusters of highly interconnected nodes can be seen within the functional interaction networks; these are often the protein complexes involved in pathways with a high probability of being regulated by the same mechanisms67. Our analysis of the functional network revealed five modules implicated in biological processes such as defense, inflammatory and acute-phase responses, response to wounding, response to lipopolysaccharide, T cell secretory granule organization, apoptosis, oxidative stress, extracellular matrix organization, signal transduction, keratinocyte differentiation, and fibroblast proliferation. Similar results of DE tear proteomics analyses have been reported in other studies24,31,45. In addition, we identified the most essential nodes as a function of their topological characteristics within the network and their interconnectivity with other nodes68. These essential central nodes in the network may provide a detailed insight into the functions involved and their relationships. Our results indicate that DE and MGD share a functional interaction network in which the proteins JUN, RELA, STAT3, TLR4, ESR1, MMP9, GNAI3, ANXA1, and TIMP1 behave as principal nodes. These proteins are implicated in biological processes such as adherens junction processes, Wnt signaling pathways, cell–cell communication, and signal transduction. These processes appear to be associated with the very nature of the studied pathologies45.

Using a second, independent analytical method, we confirmed the results of the mass spectrometry analysis for two proteins (S100A6 and CST4). We also included a third potential biomarker inferred from the network analysis (MMP9), a principal node, which has often been reported as related to DE disease6973.

To select the candidate biomarkers to validate in our study, we considered statistical significance (FDR), upregulation and downregulation fold-change, biomarkers representing different physiological processes (biological significance), number of replications (number of studies detecting a significant change), and the value of the protein in the previously reported panels of biomarkers or clinical trials. The S100A6 and CST4 proteins selected in our MS/MS study comply with the established criteria. One of these biomarkers was upregulated (S100A6), and one, downregulated (CST4). They are involved in different physiological processes (S100A6 participates in calcium binding/epithelial integrity and growth and CST4, in cysteine protease inhibitor/antimicrobial activity).

S100A6 has the most significant FDR. Although this protein showed only an intermediate fold-change value in this study, in our previous report using 2D-PAGE24, its fold-change was high (8.5). It has a suitable number of replicates (number of studies reporting similar changes): three in unbiased mass spectrometry screenings (US) and one validation in candidate immuno-detection study (CA)74. This protein has been included in a panel of DE biomarkers24. In addition, S100A6 has been used for monitoring the response of patients to changing glaucoma treatment from preserved latanoprost to preservative-free tafluprost. The increase in the level of tear S100A6 in patients treated with preserved latanoprost was reduced a year after switching to preservative-free tafluprost75, suggesting the potential application of this protein not only as a diagnostic but also as a prognostic biomarker (as in cases of topical drug-induced DE)76. In this study, CST4 showed a significant (but intermediate-level) FDR and fold-change. However, in our previous analysis using 2D-PAGE, we have found a 4.6-fold change for this protein24. In addition, CST4 has a considerable number of replicates: five in unbiased mass spectrometry screenings (US) and one validation in a candidate immuno-detection (CA)74. It has been included in a panel of DE biomarkers24. Finally, the two biomarkers, when used together as panel, have shown correct assignment (CA) of 86% in the diagnosis and classification of patients (unpublished data).

The MMP9 protein has been found in candidate immuno-detection and activity detection (CA) studies but not in unbiased mass spectrometry screening74. We selected this biomarker in our network analysis because of the nine proteins identified as central nodes, MMP9 was the protein best represented in modules (Table 4). Accordingly, the top proteins in the central nodes, in decreasing order of interaction networks (5, 4, 3, and 2), were: MMP9 (module 2) in 5 biological processes (BP); RELA and TRL4 (both in module 1) in 4 BP; ANXA 1 (module 4) and JUN and TIMP1 (module 2) in 3 BP; and STAT 3 (module 1) in 2 BP. Moreover, in module 2, MMP9 was the main player in the nodes with the most significant FDR.

The microarrays were specifically designed for 1-microliter tear sample analysis. The standardization of the miniaturized immunoassays had included several important steps. We selected the high-affinity compatible antibody pairs (using SPR and sandwich ELISA assays) and chose appropriate slide surfaces for the linkage of capture antibodies. We experimentally determined the incubation times and washing and blocking buffers. We also performed cross-reactivity studies to assure specificity of the assays and adjusted the calibration curves to clinically relevant concentration ranges (unpublished data).

The results of microarray analysis corroborated the proteomic results obtained using mass spectrometry and network analysis inference. They were comparable with the results of a study previously published by our group, in which these biomarkers were examined using commercially available ELISA kits24. The ROC curve values in the two studies are similar. Here, the results showed a high degree of correlation between the biomarkers and the clinical data, which is essential when considering a biomarker as surrogate endpoint77. However, there was no correlation with the age of patients, indicating that the biomarkers are valid for any age range.

Taking into account the scarcity of tear samples (only a few microliters can be obtained from each patient), the microarray technology is a promising alternative to other proteomic techniques for simultaneous protein/peptide quantification. Another important advantage of this technique is that a large number of proteins can be examined at the same time. Multifactorial pathologies such as DE are complex and involve many proteins and biological processes; a simultaneous analysis of multiple proteins in the same sample can provide a broad overview of the changes in protein expression during disease development.

We are currently evaluating the selected biomarkers (using multiplexed quantification techniques) in new tear samples obtained from different pathologies.

In summary, we performed a study of differential protein expression in DE and MGD disorders, using quantitative proteomics based on label-free MS. We studied the behavior of these two pathologies, from a biological perspective, using functional interaction networks. The analysis revealed nine central nodes, represented by the MMP9, JUN, RELA, STAT3, TLR4, ESR1, GNAI3, ANXA1, and TIMP1 proteins. These essential nodes might determine the biological functions to be studied in the future. Such studies should help to understand the diverse pathologies and accelerate the development of new effective therapeutic agents and strategies. Our comparison with the results obtained using different proteomic technologies corroborated our candidates for stable and robust molecular markers of the studied diseases. Among these, the S100A6 and CST4 proteins are some of the most noteworthy. They are now confirmed by orthogonal validation using ELISA24 and the customized microarray assays employed here. Thus, they are validated using two different technologies, in different tear samples and in independent assays. Finally, we explored the correlation of selected biomarkers with clinical parameters, obtaining good correlation levels and confirming the validity of these proteins as tear biomarkers for DE. However, additional coordinated interlaboratory validation should be conducted to test the potential tear biomarkers in different populations.

Acknowledgements

Mass spectrometry analysis was performed at the Proteomics Core Facility SGIKer (member of ProteoRed-ISCIII) at the University of the Basque Country. This work was partially supported by the Spanish Ministry of Science and Innovation through the INNPACTO subprogram, grant IPT-2011-1429-010000, and the Centre for the Development of Industrial Technology (CDTI), NEOTEC Program, grant IDI-20080118. The funding organizations had no role in the design or conduct of this research.

Author Contributions

Author Contributions

T.S., J.S., and A.A. designed the study, interpreted the data, and wrote the manuscript; J.S., A.A., T.S. and S.S., conceived and designed the experiments; J.S., N.G., S.R., N.B., and W.S., carried out the experiments; J.S. performed all data analysis; J.M. and J.D recruited subjects, and provided essential samples; T.S., S.S., F.B., and H.P. provided supervision, contributed essential materials; all authors reviewed the manuscript.

Notes

Competing Interests

The authors declare that the results reported in this work may be of commercial interest to Bioftalmik Applied Research.

Footnotes

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Craig JP, Tomlinson A. Importance of the lipid layer in human tear film stability and evaporation. Optom. Vis. Sci. 1997;74:8–13. doi: 10.1097/00006324-199701000-00014. [PubMed] [Cross Ref]
2. Gulati A, Sacchetti M, Bonini S, Dana R. Chemokine receptor CCR5 expression in conjunctival epithelium of patients with dry eye syndrome. Arch. Ophthalmol. 2006;124:710–6. doi: 10.1001/archopht.124.5.710. [PubMed] [Cross Ref]
3. De Paiva CS, et al. Dry eye-induced conjunctival epithelial squamous metaplasia is modulated by interferon-gamma. Invest. Ophthalmol. Vis. Sci. 2007;48:2553–60. doi: 10.1167/iovs.07-0069. [PubMed] [Cross Ref]
4. Tseng SC. Staging of conjunctival squamous metaplasia by impression cytology. Ophthalmology. 1985;92:728–33. doi: 10.1016/S0161-6420(85)33967-2. [PubMed] [Cross Ref]
5. Baudouin C, et al. Role of hyperosmolarity in the pathogenesis and management of dry eye disease: Proceedings of the ocean group meeting. Ocul. Surf. 2013;11:246–258. doi: 10.1016/j.jtos.2013.07.003. [PubMed] [Cross Ref]
6. Zhu W, Smith JW, Huang C-M. Mass spectrometry-based label-free quantitative proteomics. J. Biomed. Biotechnol. 2010;2010:840518. [PMC free article] [PubMed]
7. Zhou L, et al. Identification of tear fluid biomarkers in dry eye syndrome using iTRAQ quantitative proteomics. J. Proteome Res. 2009;8:4889–905. doi: 10.1021/pr900686s. [PubMed] [Cross Ref]
8. Srinivasan S, Thangavelu M, Zhang L, Green KB, Nichols KK. iTRAQ quantitative proteomics in the analysis of tears in dry eye patients. Invest. Ophthalmol. Vis. Sci. 2012;53:5052–9. doi: 10.1167/iovs.11-9022. [PMC free article] [PubMed] [Cross Ref]
9. Tong L, Zhou L, Beuerman RW, Zhao SZ, Li XR. Association of tear proteins with Meibomian gland disease and dry eye symptoms. Br. J. Ophthalmol. 2011;95:848–52. doi: 10.1136/bjo.2010.185256. [PubMed] [Cross Ref]
10. Csősz É, et al. Quantitative analysis of proteins in the tear fluid of patients with diabetic retinopathy. J. Proteomics. 2012;75:2196–204. doi: 10.1016/j.jprot.2012.01.019. [PubMed] [Cross Ref]
11. Motoyama A, Yates JR. Multidimensional LC separations in shotgun proteomics. Anal. Chem. 2008;80:7187–93. doi: 10.1021/ac8013669. [PubMed] [Cross Ref]
12. Domon B, Aebersold R. Mass spectrometry and protein analysis. Science. 2006;312:212–7. doi: 10.1126/science.1124619. [PubMed] [Cross Ref]
13. Chen EI, Yates JR. Cancer proteomics by quantitative shotgun proteomics. Mol. Oncol. 2007;1:144–59. doi: 10.1016/j.molonc.2007.05.001. [PMC free article] [PubMed] [Cross Ref]
14. 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:29629. doi: 10.1038/srep29629. [PMC free article] [PubMed] [Cross Ref]
15. Veenstra TD. Global and targeted quantitative proteomics for biomarker discovery. J. Chromatogr. B. Analyt. Technol. Biomed. Life Sci. 2007;847:3–11. doi: 10.1016/j.jchromb.2006.09.004. [PubMed] [Cross Ref]
16. Patel VJ, et al. A comparison of labeling and label-free mass spectrometry-based proteomics approaches. J. Proteome Res. 2009;8:3752–9. doi: 10.1021/pr900080y. [PubMed] [Cross Ref]
17. Cappadona S, et al. Improved label-free LC-MS analysis by wavelet-based noise rejection. J. Biomed. Biotechnol. 2010;2010:131505. doi: 10.1155/2010/131505. [PMC free article] [PubMed] [Cross Ref]
18. Acera A, et al. Changes in tear protein profile in keratoconus disease. Eye. 2011;25:1225–33. doi: 10.1038/eye.2011.105. [PMC free article] [PubMed] [Cross Ref]
19. Li B, et al. Tear proteomic analysis of patients with type 2 diabetes and dry eye syndrome by two-dimensional nano-liquid chromatography coupled with tandem mass spectrometry. Invest. Ophthalmol. Vis. Sci. 2014;55:177–86. doi: 10.1167/iovs.13-12080. [PubMed] [Cross Ref]
20. Grus FH, et al. SELDI-TOF-MS ProteinChip array profiling of tears from patients with dry eye. Invest. Ophthalmol. Vis. Sci. 2005;46:863–76. doi: 10.1167/iovs.04-0448. [PubMed] [Cross Ref]
21. González, N. et al. Human tear peptide/protein profiling study of Dry eye and Meibomian gland dysfunction by SPE-MALDI -TOF Mass Spectrometry analyses. J. Proteomics (2013).
22. Funke S, Azimi D, Wolters D, Grus FH, Pfeiffer N. Longitudinal analysis of taurine induced effects on the tear proteome of contact lens wearers and dry eye patients using a RP-RP-Capillary-HPLC-MALDI TOF/TOF MS approach. J. Proteomics. 2012;75:3177–90. doi: 10.1016/j.jprot.2012.03.018. [PubMed] [Cross Ref]
23. Matheis N, et al. Proteomics Differentiate Between Thyroid-Associated Orbitopathy and Dry Eye Syndrome. Invest. Ophthalmol. Vis. Sci. 2015;56:2649–56. doi: 10.1167/iovs.15-16699. [PubMed] [Cross Ref]
24. Soria J, et al. Tear proteome and protein network analyses reveal a novel pentamarker panel for tear film characterization in dry eye and meibomian gland dysfunction. J. Proteomics. 2013;78:94–112. doi: 10.1016/j.jprot.2012.11.017. [PubMed] [Cross Ref]
25. Li B, et al. Tear proteomic analysis of Sjögren syndrome patients with dry eye syndrome by two-dimensional-nano-liquid chromatography coupled with tandem mass spectrometry. Sci. Rep. 2014;4:5772. doi: 10.1038/srep05772. [PMC free article] [PubMed] [Cross Ref]
26. Nichols JJ, Green-Church KB. Mass spectrometry-based proteomic analyses in contact lens-related dry eye. Cornea. 2009;28:1109–17. doi: 10.1097/ICO.0b013e3181a2ad81. [PubMed] [Cross Ref]
27. Grus FH, Sabuncuo P, Augustin AJ. Analysis of tear protein patterns of dry-eye patients using fluorescent staining dyes and two-dimensional quantification algorithms. Electrophoresis. 2001;22:1845–50. doi: 10.1002/1522-2683(200105)22:9<1845::AID-ELPS1845>3.0.CO;2-N. [PubMed] [Cross Ref]
28. Versura P, et al. Tear proteomics in evaporative dry eye disease. Eye. 2010;24:1396–402. doi: 10.1038/eye.2010.7. [PubMed] [Cross Ref]
29. Fullard RJ, Snyder C. Protein levels in nonstimulated and stimulated tears of normal human subjects. Invest. Ophthalmol. Vis. Sci. 1990;31:1119–26. [PubMed]
30. Green-Church KB, Nichols KK, Kleinholz NM, Zhang L, Nichols JJ. Investigation of the human tear film proteome using multiple proteomic approaches. Mol. Vis. 2008;14:456–70. [PMC free article] [PubMed]
31. Nelson JD, et al. The international workshop on meibomian gland dysfunction: report of the definition and classification subcommittee. Invest. Ophthalmol. Vis. Sci. 2011;52:1930–7. doi: 10.1167/iovs.10-6997b. [PMC free article] [PubMed] [Cross Ref]
32. Rivas L, Oroza MA, Perez-Esteban A, Murube-del-Castillo J. Topographical distribution of ocular surface cells by the use of impression cytology. Acta Ophthalmol. 1991;69:371–6. doi: 10.1111/j.1755-3768.1991.tb04830.x. [PubMed] [Cross Ref]
33. Nelson JD. Impression cytology. Cornea. 1988;7:71–81. doi: 10.1097/00003226-198801000-00012. [PubMed] [Cross Ref]
34. Matthiesen, R., Trelle, M. B., Højrup, P., Bunkenborg, J. & Jensen, O. N. VEMS 3.0: algorithms and computational tools for tandem mass spectrometry based identification of post-translational modifications in proteins. J. Proteome Res. 4, 2338–47 (2005). [PubMed]
35. Braisted JC, et al. The APEX Quantitative Proteomics Tool: generating protein quantitation estimates from LC-MS/MS proteomics results. BMC Bioinformatics. 2008;9:529. doi: 10.1186/1471-2105-9-529. [PMC free article] [PubMed] [Cross Ref]
36. Deutsch EW, et al. A guided tour of the Trans-Proteomic Pipeline. Proteomics. 2010;10:1150–9. doi: 10.1002/pmic.200900375. [PMC free article] [PubMed] [Cross Ref]
37. Keller A, Nesvizhskii AI, Kolker E, Aebersold R. Empirical statistical model to estimate the accuracy of peptide identifications made by MS/MS and database search. Anal. Chem. 2002;74:5383–92. doi: 10.1021/ac025747h. [PubMed] [Cross Ref]
38. Nesvizhskii AI, Keller A, Kolker E, Aebersold R. A statistical model for identifying proteins by tandem mass spectrometry. Anal. Chem. 2003;75:4646–58. doi: 10.1021/ac0341261. [PubMed] [Cross Ref]
39. Gokce E, Shuford CM, Franck WL, Dean RA, Muddiman DC. Evaluation of normalization methods on GeLC-MS/MS label-free spectral counting data to correct for variation during proteomic workflows. J. Am. Soc. Mass Spectrom. 2011;22:2199–208. doi: 10.1007/s13361-011-0237-2. [PubMed] [Cross Ref]
40. Team, R. D. C. R: A Language and Environment for Statistical Computing. (2011).
41. Shannon P, et al. Cytoscape: A software Environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–2504. doi: 10.1101/gr.1239303. [PubMed] [Cross Ref]
42. Wu G, Feng X, Stein L. A human functional protein interaction network and its application to cancer data analysis. Genome Biol. 2010;11:R53. doi: 10.1186/gb-2010-11-5-r53. [PMC free article] [PubMed] [Cross Ref]
43. Lin C-Y, et al. Hubba: hub objects analyzer–a framework of interactome hubs identification for network biology. Nucleic Acids Res. 2008;36:W438–43. doi: 10.1093/nar/gkn257. [PMC free article] [PubMed] [Cross Ref]
44. Wolffsohn, J. S. et al. TFOS DEWS II Diagnostic Methodology report. Ocul. Surf. 10.1016/j.jtos.2017.05.001 (2017). [PubMed]
45. Lemp M, et al. The definition and classification of dry eye disease: report of the definition and classification of the Dry Eye WorkShop (2007) Ocul. Surf. 2007;5:75–92. doi: 10.1016/S1542-0124(12)70081-2. [PubMed] [Cross Ref]
46. Jee D, Park SH, Kim MS, Kim EC. Antioxidant and inflammatory cytokine in tears of patients with dry eye syndrome treated with preservative-free versus preserved eye drops. Investig. Ophthalmol. Vis. Sci. 2014;55:5081–5089. [PubMed]
47. Li, X. et al. The presence of IGHG1 in human pancreatic carcinomas is associated with immune evasion mechanisms. Pancreas, 10.1097/MPA.0b013e318213d51b (2011). [PubMed]
48. Cormier R. PLA2G2A (phospholipase A2, group IIA (platelets, synovial fluid) Atlas Genet Cytogenet Oncol Haematol. 2011;15:198–207.
49. Boehm N, et al. New insights into autoantibody profiles from immune privileged sites in the eye: A glaucoma study. Brain. Behav. Immun. 2012;26:96–102. doi: 10.1016/j.bbi.2011.07.241. [PubMed] [Cross Ref]
50. McDermott AM. Antimicrobial compounds in tears. Experimental Eye Research. 2013;117:53–61. doi: 10.1016/j.exer.2013.07.014. [PMC free article] [PubMed] [Cross Ref]
51. Alves, M. et al. Dry Eye And Oxidative Stress: The Role Of Peroxidases And Hidrogen Peroxide In Rat Lacrimal Gland. Invest. Ophthalmol. Vis. Sci. 53 (2012).
52. Zhang Z, Dmitrieva NI, Park J-H, Levine RL, Burg MB. High urea and NaCl carbonylate proteins in renal cells in culture and in vivo, and high urea causes 8-oxoguanine lesions in their DNA. Proc. Natl. Acad. Sci. USA. 2004;101:9491–6. doi: 10.1073/pnas.0402961101. [PubMed] [Cross Ref]
53. Gaffney EA, Tiffany JM, Yokoi N, Bron AJ. A mass and solute balance model for tear volume and osmolarity in the normal and the dry eye. Prog. Retin. Eye Res. 2010;29:59–78. doi: 10.1016/j.preteyeres.2009.11.002. [PubMed] [Cross Ref]
54. Tsai PS, et al. Proteomic analysis of human meibomian gland secretions. Br. J. Ophthalmol. 2006;90:372–7. doi: 10.1136/bjo.2005.080846. [PMC free article] [PubMed] [Cross Ref]
55. Liu S, et al. Changes in gene expression in human meibomian gland dysfunction. Invest. Ophthalmol. Vis. Sci. 2011;52:2727–40. doi: 10.1167/iovs.10-6482. [PMC free article] [PubMed] [Cross Ref]
56. Soria, J. et al. The analysis of human conjunctival epithelium proteome in ocular surface diseases using impression cytology and 2D-DIGE. Exp. Eye Res., 10.1016/j.exer.2017.03.006 (2017). [PubMed]
57. Koo B-S, Lee D-Y, Ha H-S, Kim J-C, Kim C-W. Comparative analysis of the tear protein expression in blepharitis patients using two-dimensional electrophoresis. J. Proteome Res. 2005;4:719–24. doi: 10.1021/pr0498133. [PubMed] [Cross Ref]
58. Li, B. et al. Tear proteomic analysis of patients with type 2 diabetes and dry eye syndrome by two-dimensional-nano-liquid chromatography coupled with tandem mass spectrometry. Invest. Ophthalmol. Vis. Sci., 10.1167/iovs.13-12080 (2013). [PubMed]
59. Zhou, L. et al. Proteomic Analysis Revealed the Altered Tear Protein Profile in a Rabbit Model of Sjögren’s Syndrome-Associated Dry Eye. Proteomics10.1002/pmic.201200230 (2013). [PMC free article] [PubMed]
60. Boehm N, et al. Alterations in the tear proteome of dry eye patients–a matter of the clinical phenotype. Invest. Ophthalmol. Vis. Sci. 2013;54:2385–92. doi: 10.1167/iovs.11-8751. [PubMed] [Cross Ref]
61. Mackie IA, Seal DV. Diagnostic implications of tear protein profiles. Br. J. Ophthalmol. 1984;68:321–4. doi: 10.1136/bjo.68.5.321. [PMC free article] [PubMed] [Cross Ref]
62. Mackie IA, Seal DV. Confirmatory tests for the dry eye of Sjögren’s syndrome. Scand. J. Rheumatol. Suppl. 1986;61:220–3. [PubMed]
63. Seal DV, et al. Bacteriology and tear protein profiles of the dry eye. Br. J. Ophthalmol. 1986;70:122–5. doi: 10.1136/bjo.70.2.122. [PMC free article] [PubMed] [Cross Ref]
64. Zybailov B, Coleman MK, Florens L, Washburn MP. Correlation of relative abundance ratios derived from peptide ion chromatograms and spectrum counting for quantitative proteomic analysis using stable isotope labeling. Anal. Chem. 2005;77:6218–24. doi: 10.1021/ac050846r. [PubMed] [Cross Ref]
65. Lu P, Vogel C, Wang R, Yao X, Marcotte EM. Absolute protein expression profiling estimates the relative contributions of transcriptional and translational regulation. Nat. Biotechnol. 2007;25:117–24. doi: 10.1038/nbt1270. [PubMed] [Cross Ref]
66. Vogel C, Marcotte EM. Label-free protein quantitation using weighted spectral counting. Methods Mol. Biol. 2012;893:321–41. doi: 10.1007/978-1-61779-885-6_20. [PMC free article] [PubMed] [Cross Ref]
67. Dezso Z, Oltvai ZN, Barabási A-L. Bioinformatics analysis of experimentally determined protein complexes in the yeast Saccharomyces cerevisiae. Genome Res. 2003;13:2450–4. doi: 10.1101/gr.1073603. [PubMed] [Cross Ref]
68. Ekman D, Light S, Björklund AK, Elofsson A. What properties characterize the hub proteins of the protein-protein interaction network of Saccharomyces cerevisiae? Genome Biol. 2006;7:R45. doi: 10.1186/gb-2006-7-6-r45. [PMC free article] [PubMed] [Cross Ref]
69. Uchino Y, et al. Alteration of galectin-3 in tears of patients with dry eye disease. Am. J. Ophthalmol. 2015;159:1027–1035.e3. doi: 10.1016/j.ajo.2015.02.008. [PMC free article] [PubMed] [Cross Ref]
70. Aragona P, et al. Matrix metalloproteinase 9 and transglutaminase 2 expression at the ocular surface in patients with different forms of dry eye disease. Ophthalmology. 2015;122:62–71. doi: 10.1016/j.ophtha.2014.07.048. [PubMed] [Cross Ref]
71. López-Miguel, A. et al. Clinical and Molecular Inflammatory Response in Sjögren Syndrome-Associated Dry Eye Patients Under Desiccating Stress. Am. J. Ophthalmol. 161, 133-41-2 (2016). [PubMed]
72. Messmer, E. M., von Lindenfels, V., Garbe, A. & Kampik, A. Matrix Metalloproteinase 9 Testing in Dry Eye Disease Using a Commercially Available Point-of-Care Immunoassay. Ophthalmology10.1016/j.ophtha.2016.07.028 (2016). [PubMed]
73. Chan TCY, Ye C, Chan KP, Chu KO, Jhanji V. Evaluation of point-of-care test for elevated tear matrix metalloproteinase 9 in post-LASIK dry eyes. Br. J. Ophthalmol. 2016;100:1188–91. doi: 10.1136/bjophthalmol-2015-307607. [PubMed] [Cross Ref]
74. Willcox, M. D. P. et al. TFOS DEWS II Tear Film Report. Ocul. Surf., 10.1016/j.jtos.2017.03.006 (2017). [PubMed]
75. Uusitalo, H., Aapola, U., Jylhä, A., Näattinen, J. & Beuerman, R. Improved tear fluid proteome and dynamics after switch from preserved latanoprost to preservative free tafluprost. A 1-year follow-up study. Investig Ophthalmol Vis Sci56 (2015).
76. Gomes, J. A. P. et al. TFOS DEWS II iatrogenic report. Ocul. Surf., 10.1016/j.jtos.2017.05.004 (2017). [PubMed]
77. Firestein GS. A biomarker by any other name. Nat. Clin. Pract. Rheumatol. 2006;2:635. doi: 10.1038/ncprheum0347. [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group