64
Views
15
CrossRef citations to date
0
Altmetric
Original Research

A single-center study evaluating the effect of the controlled adverse environment (CAESM) model on tear film stability

, , , , , & show all
Pages 1865-1872 | Published online: 13 Nov 2012

Abstract

Purpose

To investigate use of an improved ocular tear film analysis protocol (OPI 2.0) in the Controlled Adverse Environment (CAESM) model of dry eye disease, and to examine the utility of new metrics in the identification of subpopulations of dry eye patients.

Methods

Thirty-three dry eye subjects completed a single-center, single-visit, pilot CAE study. The primary endpoint was mean break-up area (MBA) as assessed by the OPI 2.0 system. Secondary endpoints included corneal fluorescein staining, tear film break-up time, and OPI 2.0 system measurements. Subjects were also asked to rate their ocular discomfort throughout the CAE. Dry eye endpoints were measured at baseline, immediately following a 90-minute CAE exposure, and again 30 minutes after exposure.

Results

The post-CAE measurements of MBA showed a statistically significant decrease from the baseline measurements. The decrease was relatively specific to those patients with moderate to severe dry eye, as measured by baseline MBA. Secondary endpoints including palpebral fissure size, corneal staining, and redness, also showed significant changes when pre- and post-CAE measurements were compared. A correlation analysis identified specific associations between MBA, blink rate, and palpebral fissure size. Comparison of MBA responses allowed us to identify subpopulations of subjects who exhibited different compensatory mechanisms in response to CAE challenge. Of note, none of the measures of tear film break-up time showed statistically significant changes or correlations in pre-, versus post-CAE measures.

Conclusion

This pilot study confirms that the tear film metric MBA can detect changes in the ocular surface induced by a CAE, and that these changes are correlated with other, established measures of dry eye disease. The observed decrease in MBA following CAE exposure demonstrates that compensatory mechanisms are initiated during the CAE exposure, and that this compensation may provide the means to identify and characterize clinically relevant subpopulations of dry eye patients.

Introduction

Dry eye disease is a term used to describe a collection of disorders with a shared diagnosis of tear film dysfunction that lead to decreased visual acuity, ocular pain, burning, and the potential for corneal scarring.Citation1 Prevalence is highest in older individuals and in women, as well those who have previously undergone laser vision correction.Citation1,Citation2 The most recent reports suggest that moderate to severe dry eye currently affects between 5 and 10 million Americans, with estimates of ten times that number worldwide.Citation2 Current treatments include artificial tears, tear plugs, or immunosuppressive drugs such as cyclosporine, but these treatments are often ineffective for many dry eye sufferers, and so, the unmet need for new therapeutics is significant.Citation2

Understanding of the pathophysiology of dry eye has advanced considerably in recent years, and what was once thought to be a condition due simply to insufficient tear production, is now recognized as a multifactorial collection of diseases.Citation3 This is a result of significant strides in ocular surface biology research, and improvements in clinical assessment techniques such as, corneal staining, blink pattern analysis, and various measures of tear film stability.Citation4,Citation5 In two recent papers we described the stepwise development of our improved method of tear film analysis, which we have designated the ocular tear film analysis protocol OPI 2.0 (Ora, Inc. Andover MA, USA; patent pending) system.Citation6,Citation7 This method combines a number of optimized parameters with automated video data capture and analysis, to generate a more objective, quantitative measure of tear film dynamics. These improvements have the potential to substantially enhance tear film metrics, and represent a key advance over previous methods of tear film analysis.

Tear film stability and blink behavior are inexorably linked. The tear film is established by the sweeping, squeegee-like action of the lids, and the rate and pattern of blinks is in turn, modulated by feedback from corneal sensory nerves.Citation3,Citation8 Efforts to measure properties of the tear film led to the development of the tear film break-up time (TFBUT), a method of analysis in which subjects are asked to refrain from blinking while an observer monitors the integrity of the tear film.Citation8 This “forced stare” approach allowed the first estimates of inherent tear film stability and provided the means for addressing the impact of different disorders, drug treatments, or environmental conditions, on the physiochemical properties (and therefore, functional attributes) of the tear film.Citation1,Citation4 The standard values developed for use with this method were >10 seconds for normal subjects, and <10 seconds for subjects with dry eye.Citation1,Citation4

More recently, studies have identified limitations to both the method used to measure TFBUT, as well as the metric itself. For example, it was found that reduction and standardization of the quantity of fluorescein used, led to a shift of the reference values for TFBUT (to a mean value of 7 seconds for normal subjects, and 2.5 seconds for those with dry eye).Citation5 Similarly, our goal in developing the OPI 2.0 system was to enhance and streamline dry eye metrics by measuring tear film stability under conditions of natural, rather than forced stare blinking.Citation6,Citation7 Furthermore, by capturing the natural dynamics of the tear film with automated methods, it is possible to expand upon previous studies of the interaction between blinking and tear film behavior.Citation8Citation13

Blinking is a reflex function regulated by the combination of autonomic inputs and sensory feedback about environmental conditions.Citation14 Blink rate is known to increase under adverse conditions such as those presented by the Controlled Adverse Environment (CAESM), a clinical model that reproduces a standard ocular challenge (by regulation of humidity, temperature, airflow, lighting conditions, and visual tasking) for the investigation of treatments for dry eye.Citation1,Citation4,Citation13,Citation15,Citation16 Our early assessments of both blink behavior and tear film stability led to our development of the Ocular Protection Index (OPI), a tool used in observational and clinical studies of dry eye, to evaluate the interaction between blinking and TFBUT.Citation9,Citation13Citation16 We defined OPI as the ratio of interblink interval (IBI) and TFBUT; lower values of the index and particularly values less than one, were associated with increased risk of keratitis (since on average, break-up of the tear film and subsequent corneal exposure would occur prior to the next blink).Citation9 Our more recent studies employed break-up and blink data captured with the OPI 2.0 system, to yield another estimate of corneal surface exposure.Citation7 This alternative index relied upon the mean break-up area (MBA) rather than a break-up time elicited under forced stare conditions.Citation6

The OPI 2.0 system was designed to evaluate ocular surface protection under a normal blink pattern and normal visual conditions, using fully automated software algorithms that provide real-time measurement of corneal exposure.Citation7 This system provides a simultaneous measurement of both break-up area and blink rate. Video capture allows calculation of a mean corneal surface exposure (we abbreviate this as mean break-up area, MBA) as well as additional metrics, such as palpebral fissure width. In an initial paper, we showed that video capture of MBA and blink rates was superior to the combination of forced stare TFBUT and IBI in distinguishing normal subjects from those with signs and symptoms of dry eye.Citation6 In a second paper, we demonstrated that our automated method provided accurate, reliable measures of corneal surface exposure, and again showed that these values could be used to distinguish between normal subjects and subjects with dry eye.Citation7

In this study, we explored the ability of the OPI 2.0 system to identify the changes and modifications of the tear film after exposing dry eye subjects to the CAE. A key aspect of the CAE is its utility in distinguishing subpopulations of dry eye patients. Subjects challenged by environmental changes (such as those presented by the CAE) normally respond with some degree of physiological compensation, and previous studies have shown that the ability of these mechanisms to adequately compensate for environmental challenges is reduced in those with dry eye.Citation17 Therefore, the nature and extent of the compensatory response will likely be a reflection of the underlying tear film pathology. These compensatory mechanisms, such as changes in blink rate or reflex tearing, are likely to have an impact on properties that are measured using the OPI 2.0 system. Our primary goal in developing new and refined metrics was their use as tools to identify and characterize patient subpopulations, particularly in multifactorial diseases such as dry eye. As the next step in the validation of OPI 2.0 System-based measures, we examined dry eye subjects before, immediately after, and 30 minutes after a 90-minute CAE challenge.

Methods

Inclusion criteria

Subjects were recruited from an existing database of dry eye patients, and a total of 33 subjects were enrolled. Criteria for inclusion in the database included a history of dry eye, use of artificial tears, and a Schirmer’s test scoreCitation4 of <5 mm in at least one eye. Enrolled subjects were at least 18 years of age, had a history of use, or desire to use an eye drop for dry eye symptoms within the past 6 months, and had a best corrected visual acuity of +0.7 or better in both eyes, as assessed by the Early Treatment of Diabetic Retinopathy Study scale. Subjects also had to satisfy each of the following criteria at baseline: a forced stare TFBUT score of <5 seconds in at least one eye; a deficient OPI (<1) during at least 30% of interblink intervals, as determined by a trained technician; and a total corneal fluorescein staining score of ≥3, based on the sum of the central, superior, and inferior regions of the cornea, as defined by the Ora scale (see below).

Exclusion criteria

Subjects were excluded from the study if they had clinically significant anterior blepharitis in the opinion of the investigator; were diagnosed with ongoing ocular infection (bacterial, viral, or fungal), or active ocular inflammation (eg, follicular conjunctivitis); wore contact lenses in the previous week; had used any eye drop in the 4 hours prior to the study; had previously undergone laser in situ keratomileusis surgery, or any other ocular surgery in the past year; were currently taking any topical ophthalmic prescription, or over-the-counter solutions, artificial tears, gels or scrubs that could not be discontinued for the duration of the trial; had used Restasis® in the previous 30 days; had a systemic disease, or uncontrolled medical condition that could interfere with study measurements or subject compliance; were currently pregnant or nursing; or had received an experimental drug or device within 30 days prior to the visit.

Study design

This was a single-center study, conducted in a single visit that included a 90-minute session in the CAE. Written informed consent was obtained prior to study procedures. Patient-reported and investigator-observed adverse events were captured and monitored for the duration of the study. A summary of the protocol, including the tests which were conducted, is shown in . The primary endpoint for this study was MBA prior to CAE exposure compared with MBA immediately and 30 minutes post-CAE exposure. Secondary endpoints included fluorescein staining, TFBUT, and redness prior to CAE compared with fluorescein staining, TFBUT, and redness immediately- and 30 minutes post-CAE exposure. The scoring of corneal staining was based on a 5-point Ora scale (0–4) where 0 = no staining and 4 = confluent fluorescein staining across the entire cornea. Redness was scored using a 5-point Ora scale (0–4), where 0 = no detectable vasodilation and 4 = broad, prominent ciliary and conjunctival vasodilation. Both scales have been used previously.Citation15,Citation16 Additional secondary endpoints included video-based measurements (including IBI and palpebral fissure width), as well as TFBUT, collected at all three time points.

Table 1 Overview of study protocol

Statistical analysis

Thirty-three subjects were enrolled in the study. Of the 66 eyes under study, 65 provided readable videos, and data from these 65 eyes comprised the complete analysis sample. All 33 subjects in the intent-to-treat population met the requirements of the protocol criteria. We also analyzed subgroups that were defined by their baseline MBA measure. The first of these were subjects with an initial MBA (MBAI) > 0.2, a value used to identify those with moderate to severe dry eye. This cutoff was based on previously published studies that showed subjects who reported dry eye signs and symptoms had a mean MBA score of 0.232.Citation5 Demographic variables (age, sex, duration of dry eye disease) were summarized by means and standard deviations. An additional measure, ocular discomfort, was provided by the subjects and was collected every 5 minutes during the 90-minute CAE, using a 5-point ordinal scale. Subjects scored their discomfort, and these measurements were used to calculate average discomfort during the CAE, and tearing. We defined tearing as the time at which discomfort either reduced or reached a plateau at a value less than the maximum, based upon previous studies demonstrating that this plateau is associated with a compensatory tearing response.Citation17

For those variables with long right-tailed distributions (including MBA, OPI2, blink rate, IBI, palpebral fissure, and break-up time) we used a gamma multiplicative model to obtain estimates for pre-CAE and post-CAE means, ratios of means, 95% confidence intervals, and P-values for tests of equality. Corneal fluorescein staining and conjunctival redness were analyzed using normal linear models. We used a generalized estimating equation to accommodate the within-subject correlation between eyes for both models. For the latter purpose, a sandwich variance was used in conjunction with a working independence correlation structure. All models were fit using the Genmod module of SAS version 9.2 (SAS Institute Inc, Cary, NC).Citation18

For variables measured pre- and post-CAE, bivariate relationships were obtained via correlations between change scores (eg, post-CAE MBA minus pre-CAE MBA versus pre-CAE blink rate minus post-CAE blink rate). These change scores were also correlated with variables collected in the CAE (average discomfort, and tearing [score = 1 if there was tearing within 90 minutes, score = 0 otherwise]), and demographic variables of age, duration of dry eye disease, and gender.

Results

In this single-visit study, we examined the effect of the CAE on MBA for subjects with confirmed dry eye disease. A total of 33 dry eye subjects completed the study and data from 65 qualified eyes were entered into the analysis algorithm. One video file was corrupted during the course of collection and could not be included in the analysis. No treatment was administered in the course of the study.

The demographics of the study populations were generally representative of the larger population of all dry eye patients.Citation1 The study population was 34% male, had a mean age of 64 years, and had experienced dry eye disease for an average of 13.4 years; a similar demographic profile also described all of the subpopulations analyzed in this study, as summarized in .

Table 2 Study populations: demographics and initial MBA values*

The primary endpoint was the change in MBA measured at baseline, immediately following CAE, and 30 minutes after exposure to the CAE, for subjects with MBAI measures of 0.2 or greater. The change in MBA for three subject populations (all subjects, MBAI scores < 0.2, MBAI scores > 0.2) is shown in . Subjects with higher MBAI are those with more severe dry eye,Citation7 and we found that in this population, as well as in the entire subject population, exposure to the CAESM caused a significant decrease in MBA. This figure also shows that for those subjects with less severe dry eye (MBAI < 0.2), there was a trend toward increased MBA which was not statistically significant.

Figure 1 A Comparison of endpoints before and after CAESM exposure.

Notes: Graphs show mean values (±SD) for the three endpoints in the total population, as well as for subjects with MBAI > 0.2 (severe dry eye), and MBAI < 0.2 (moderate dry eye). (A) Decrease in MBA was limited to subjects with MBAI > 0.2. (B) Compensatory increases in blink rates were primarily seen in subjects with MBAI <0.2. (C) All subjects experienced similar increases in corneal staining.
Abbreviations: CAESM, controlled adverse environment; MBA, mean break-up area; MBAI, initial mean break-up area.
Figure 1 A Comparison of endpoints before and after CAESM exposure.

show comparisons of two additional metrics in these three populations: blink rate and corneal staining. The change in blink rate that was observed over the course of the study is interesting in that while the population as a whole increased significantly, the group with MBAI > 0.2 did not exhibit this increase ( and ). Corneal fluorescein staining was noteworthy for two reasons. First, the difference between those with MBAI > 0.2 and those with MBAI < 0.2 was significant (P <0.001), consistent with our premise that those with higher MBA scores have more severe dry eye (). Second, all subjects showed significant increases in staining over the course of the visit, as expected for exposure to the CAE (, ). Of note, this increase in corneal staining occurred regardless of the change in MBA observed for all subjects.

Table 3 CAESM effects on measures of dry eye

Secondary endpoints included blink rate, palpebral fissure size, corneal staining, conjunctival redness, TFBUT, and IBI. The mean values for these parameters are shown in ; differences and P-values between pre-CAESM and post-CAESM values are shown in .

Table 4 Difference values for primary and secondary endpoints in different subpopulations

The major finding was the statistically significant decrease in the MBA immediately following CAE exposure. MBA values 30 minutes post-CAE were also significantly lower than baseline, although there was some recovery between the two post-CAE time points. As expected, both fluorescein staining and redness increased with CAE exposure ( and ). Blink rate and palpebral fissure width also decreased significantly across the population, at both post-CAESM times. In contrast, mean IBI and TFBUT changed only minimally, and not significantly. summarizes the difference scores for all endpoints, and shows the same pattern of significance for primary and secondary endpoints seen in mean value comparisons.

Based upon the pattern of decreased MBA and increased blink rates, we examined a second set of subpopulations that were defined using a higher threshold of MBAI (>0.5; n = 19). Even at this higher threshold, all subjects who displayed a decrease in MBA were included. These subjects had a higher initial corneal staining score, and exhibited a pronounced decrease in MBA in response to CAE exposure (, ). Analysis of these subpopulations (initial MBA > 0.2 or > 0.5) suggests that the MBA metric may provide the means to distinguish between individuals who respond in the CAESM with differing compensatory mechanisms. Comparison of subjects with either high or low MBA showed that there was a clear, statistically significant difference in the way these two groups responded to the CAESM. and show that the decrease in MBA observed in the full study population was entirely due to the decrease observed in the subpopulation with the higher MBAI, and that for the two lower MBAI groups, the mean value increased slightly over the course of the CAESM exposure (). A second distinction between those with low versus high initial MBA values is shown in and ; despite a decrease in MBA during the CAESM, subjects with high MBA values showed no change in blink rate, while those with low initial MBA values increased their blink rate in the CAESM by almost 50% (). Reduction in palpebral fissure accompanied this increase in blink rate. While it is reasonable to point out that the largest decrease in MBA may simply be due to a high MBAI, the differences in blink rates clearly represent a difference in compensatory responses.

Figure 2 Subpopulation comparisons.

Notes: Graphs show mean values (±SD) pre, immediately post and 30 min post CAESM for two alternative subpopulations of subjects with MBAI > 0.5 (n = 19) and MBAI < 0.5 (n = 46). (A) and (B) show that the decrease in MBA in the “severe” subpopulation was primarily seen in subjects with the highest MBAI (>0.5) who did not respond to the CAE with increased blink rates (seen in subjects with milder dry eye). Panel (C) shows that forced stare TFBUT was not altered in the total population, nor in these subgroups.
Abbreviations: MBA, mean break-up area; MBAI, initial mean break-up area; CAESM, controlled adverse environment; TFBUT, tear film break-up time.
Figure 2 Subpopulation comparisons.

TFBUT has been used as a standard metric in dry eye studies for many years, but in this study did not reflect the changes in dry eye signs and symptoms resulting from CAESM exposure. and and , show that none of the populations examined in this study showed a significant change in TFBUT, despite the fact that other metrics associated with dry eye, including corneal staining, ocular redness, and discomfort (not shown) increased over the time course of the CAESM.

We next did a correlation analysis to test whether specific endpoints might show an association with MBA changes; these data are summarized in for the population as a whole, and for the subpopulation of subjects with MBAI > 0.2. Of the endpoints analyzed, only tearing and blink rate (and the related parameter, IBI) were significantly correlated with MBA measures. In addition, all endpoints except IBI showed a positive correlation with MBA in the full study population. When this same analysis was applied to the high and low MBA subpopulations, the correlation between blink rate and MBA was restricted to the high MBA population, while subjects in the low MBA population exhibited an inverse correlation with these two parameters. This distinction indicates that the two subpopulations may be responding to the CAESM with different compensatory mechanisms. Collectively our results suggest that subjects in the low MBA group respond to the effects of the CAESM, at least in part, by increasing their blink rate, while the high MBA group does not. Despite this lack of change in compensatory metrics, subjects in the high MBA group exhibited a significant decrease in MBA and therefore are clearly responding to the CAESM in some way.

Table 5 Pearson correlation coefficients between MBA and secondary endpoint measures

Throughout the course of CAESM exposure, subjects were asked to rate their ocular discomfort (not shown). The values for all groups displayed a slow, consistent increase in scores over the time course of the CAESM, but none of the differences between group scores were statistically significant at any time point (that is, discomfort could not explain the differences in responses seen in the two groups).

Discussion

Dry eye disease is an exceedingly complex disease because of its variety of etiologies and the overlapping, interacting sensory elements and response mechanisms required for an optimally-tuned tear film. The identification of endpoints that are both meaningful and measurable has been problematic. The study presented here represents one step in the process of establishing robust quantitative metrics for clinical studies. Our goal was to identify measures that would be responsive in clinical models designed to replicate the disease process, and that could distinguish between the subpopulations that are characteristic of this disease. The data presented here demonstrated that the OPI 2.0 system and the measure of MBA can provide the assessment tools necessary to meet this goal.

The key finding of this study was that the metric of MBA was shown to identify subpopulations of dry eye subjects who displayed distinct compensatory mechanisms when challenged with an adverse environment. In contrast, the traditional metric of tear film stability, TFBUT, was not significantly altered by CAESM exposure. Patients also exhibited significant increases in corneal staining, ocular redness, and decreases in palpebral fissure width, which are all characteristic of dry eye disease. Thus, in a clinical model which reliably elicits signs and symptoms of dry eye disease, MBA provided a useful new metric that was superior to TFBUT.

Patients completing the CAESM exposure fell into two groups. The first group, which comprised about 70% of all subjects, was distinguished by a relatively stable MBA that was maintained in part, by an increase in blink rate and a decrease in palpebral fissure width. These changes indicated they were able to respond to the environmental challenge with these (and perhaps other) mechanisms in order to maintain a relatively constant corneal surface exposure, as measured by MBA. The remaining 30% of subjects began the study with levels of corneal staining similar to the low MBA group despite a baseline MBA that was 10-fold higher. Subjects in this group responded to the CAE by lowering their MBA more than 3-fold during the course of the CAE exposure. While we speculate that this response employed some combination of increased tearing, mucin secretion, or meibum expression, our study did not examine these specific tear film parameters. The net effect of their response however was evidenced by the fact that corneal staining between the two groups was comparable. Future studies will benefit from inclusion of techniques that can monitor changes in these tear film components during the course of CAE exposure.

Our study has direct impact on the design of therapeutic development strategies going forward. First, we provided direct evidence that break-up area, and not break-up time, is the more valuable parameter in studies of induced dry eye disease. Second, we established that the metrics of the OPI 2.0 system can distinguish subpopulations of subjects who are likely to require different therapeutic strategies for successful amelioration of their dry eye signs and symptoms. An important next step in the refinement of the OPI 2.0 system will be to associate the subgroups identified using metrics such as MBA, with clinically relevant patient populations.

Disclosure

The authors reported no conflicts of interest in this work.

References

  • The definition and classification of dry eye disease: report of the Definition and Classification Subcommittee of the International Dry Eye WorkShop2007Ocul Surf200752759217508116
  • PflugfelderSCPrevalence, burden, and pharmacoeconomics of dry eye diseaseAm J Manag Care2008143 SupplS102S10618452369
  • StevensonWChauhanSKDanaRDry eye disease: an immunemediated ocular surface disorderArch Ophthalmol20123019010022232476
  • OuslerGWGomesPJWelchDAbelsonMBMethodologies for the study of ocular surface diseaseOcul. Surf20053314315417131019
  • AbelsonMBOuslerGWNallyLAWelchDKrenzerKAlternative reference values for tear film break up time in normal and dry eye populationsSullivanDASternMETsubotaKLacrimal Gland, Tear Film, and Dry Eye Syndromes 3Part B: Basic Science and Clinical RelevanceNew YorkKluwer Academic/Plenum P200211211125
  • AbelsonRLaneKJAngjeliEJohnstonPOuslerGMontgomeryDMeasurement of ocular surface protection under natural blink conditionsClin Ophthalmol201151349135722034554
  • AbelsonRLaneKJRodriguezJValidation and verification of the OPI 2.0 SystemClin Opthalmol20126613622
  • Methodologies to diagnose and monitor dry eye disease: report of the Diagnostic Methodology Subcommittee of the International Dry Eye Workshop2007Ocul Surf20075210815217508118
  • OuslerGW3rdHagbergKWSchindelarMWelchDAbelsonMBThe Ocular Protection IndexCornea200827550951318520496
  • OuslerGWMichaelsonCChristensenMTAn evaluation of tear film breakup time extension and ocular protection index scores among three marketed lubricant eye dropsCornea200726894995217721294
  • SimmonsPAVehigeJGClinical performance of a mid-viscosity artificial tear for dry eye treatmentCornea200726329430217413956
  • TorkildsenGLOuslerGW3rdGomesPOcular comfort and drying effects of three topical antihistamine/mast cell stabilizers in adults with allergic conjunctivitis: a randomized, double-masked crossover studyClin. Ther20083071264127118691985
  • RolandoMAutoriSBadinoFBarabinoSProtecting the ocular surface and improving the quality of life of dry eye patients: a study of the efficacy of an HP-guar containing ocular lubricant in a population of dry eye patientsJ Ocul Pharmacol Ther200925327127819366323
  • CruzAAGarciaDMPintoCTCechettiSPSpontaneous eyeblink activityOcul Surf201191294121338567
  • SembaCPTorkildsenGLLonsdaleJDA phase 2 randomized, double-masked, placebo-controlled study of a novel integrin antagonist (SAR 1118) for the treatment of dry eyeAm J Ophthalmol201215361050106022330307
  • PataneMACohenAFromSTorkildsenGWelchDOuslerGW3rdOcular iontophoresis of EGP-437 (dexamethasone phosphate) in dry eye patients: results of a randomized clinical trialClin Ophthalmol2011563364321629568
  • OuslerGWAbelsonMBNallyLAEvaluation of the time to “natural compensation” in normal and dry eye subject populations during exposure to a controlled adverse environmentSullivanDASternMETsubotaKLacrimal Gland, Tear Film, and Dry Eye Syndromes 3, Part B: Basic Science and Clinical RelevanceNew YorkKluwer Academic/Plenum P200210571063
  • SAS Institute Inc2009SAS OnlineDoc® 9.2Cary, NCSAS Institute Inc