Purpose
Dry eye disease (DED) is a common ocular surface disorder associated with ocular discomfort, visual disturbance, and reduced quality of life. Vitamin D has anti-inflammatory and immunomodulatory properties, and deficiency may predispose individuals to ocular surface disease. This study evaluated the association between vitamin D deficiency and incident DED.
Design
Retrospective cohort study.
Subjects
Adults aged ≥18 years identified from the TriNetX U.S. Collaborative Network. After 1:1 propensity score matching, 6,047,502 patients with vitamin D deficiency were compared with 6,047,502 matched controls without vitamin D deficiency.
Methods
De-identified electronic health records with up to 20 years of follow-up were analyzed. Vitamin D deficiency was defined using ICD-10 code E55. Controls were patients undergoing routine medical examinations without abnormal findings (ICD-10: Z00.00). Propensity score matching was performed for age, sex, and relevant comorbidities. Time-to-event analyses were conducted using Cox proportional hazards models. Kaplan–Meier curves and log-rank tests compared cumulative incidence between cohorts. Risk ratios (RRs) and odds ratios (ORs) were calculated from end-of-follow-up incidence proportions.
Main Outcome Measures
Incident dry eye disease, defined by ICD-10 code H04.12.
Results
DED developed in 196,639 patients (3.3%) with vitamin D deficiency and 160,141 patients (2.7%) without deficiency. Vitamin D deficiency was associated with a significantly increased risk of DED (hazard ratio [HR], 1.286; 95% CI, 1.277–1.294; p <.001). Kaplan–Meier curves demonstrated a higher cumulative incidence in the vitamin D–deficient cohort (log-rank χ² = 5590.618; p <.001). The risk ratio was 1.242 (95% CI, 1.234-1.250; p <.001), and the odds ratio was 1.250 (95% CI, 1.242-1.259; p <.001). Limitations include reliance on diagnostic coding and lack of serum vitamin D level data.
Conclusions
Vitamin D deficiency was significantly associated with an increased risk of developing DED, with a 28.6% higher hazard compared with matched controls. These findings suggest that vitamin D deficiency may represent a modifiable risk factor for DED and support further prospective studies to evaluate the role of vitamin D supplementation in prevention and management.
INTRODUCTION
Dry eye disease (DED) is a heterogenous ocular surface condition characterized by a loss of homeostasis and associated ocular symptoms. Tear film instability, hyperosmolarity, ocular surface inflammation and damage and neurosensory abnormalities play etiological roles. A meta‐analysis of United States (US) studies estimated a DED prevalence of approximately 8.1% (95% confidence interval (CI) 4.9-13.1%) in the adult population, acknowledging that different diagnostic definitions (symptoms and/or signs) were used across studies. Despite differences in prevalence reports, DED is widely acknowledged as a common condition that drives a large proportion of eye clinic visits and imposes a substantial burden on patients and the healthcare system. ,, Functionally, ocular pain symptoms coupled with fluctuating vision can impact the ability to perform activities of daily living and significantly reduce quality of life. , Current DED therapies target multiple facets of disease, including tear supplementation with ocular lubricants, anti-inflammatory treatment with limited-duration topical corticosteroids or immunomodulators such as cyclosporine, and Meibomian gland dysfunction (MGD)-directed therapies including lid/physical measures and oral macrolide or tetracycline antibiotics; dietary omega-3/polyunsaturated fatty acids have also been investigated as adjunctive options. However, these therapies vary in effectiveness and are often associated with burdensome side effects. The chronic burden of DED and the limited efficacy of existing treatments highlight the need for effective preventive and new therapeutic strategies.
Vitamin D is a multifunctional hormone whose main physiological role is to maintain calcium and phosphate balance, mineral metabolism, and skeletal health. It also has important anti-inflammatory and immune-modulating effects, mediated through regulation of gene transcription, and it can reduce inflammatory signaling and influence programmed cell death. Given that DED is multifactorial with inflammation as an important component, vitamin D may serve as a potential mediator in its pathogenesis.
Several studies have explored the relationship between vitamin D and DED with mixed findings. Some have found that vitamin D deficiency impacts symptoms and signs of DED. In a prospective clinical study of 34 individuals with serum vitamin D deficiency and 21 controls with normal vitamin D levels, the case group demonstrated significantly lower tear film break-up time (5.18 ± 2.15 vs 7.36 ± 3.10 seconds; p =.01) and Schirmer test scores (12.18 ± 6.44 vs 18.57 ± 8.99 mm; p =.007) compared to controls, supporting an association between vitamin D deficiency and impaired tear film stability and production. Vitamin D deficiency has also been implicated in DED symptoms. A 2020 systematic review and meta-analysis included 10 clinical studies (18,919 participants). Pooling 5 studies with DED as the exposure, heterogeneously defined in each study by a mix of symptoms and/or signs, and 5 studies with serum vitamin D as the exposure, the authors found that patients with DED had mean serum 25-hydroxyvitamin D levels 3.99 ng/mL lower than controls. In subgroup meta-analysis restricted to 5 studies that used vitamin D deficiency as the exposure, individuals with vitamin D deficiency had OSDI scores 10.70 points higher and Schirmer test values 6.38 mm/5 min lower than controls, while tear break-up time was similar between the groups, indicating that vitamin D deficiency is associated with a more severe DED phenotype across symptoms and signs. On the other hand, other studies did not find relationships between vitamin D and DED. A cross-sectional study of 740 Korean adults found no significant relationship between serum vitamin D levels and OSDI-measured dry eye symptoms. A similar pattern was observed in a study of 247 US male Veterans, where serum vitamin D levels showed no relationship with tear parameters. Taken together, these studies suggest that vitamin D deficiency may influence certain aspects of DED, although results across the literature remain inconsistent.
Vitamin D supplementation has been investigated as a treatment for dry eye disease, with predominantly beneficial effects, particularly among individuals with vitamin D deficiency. In a 2024 systematic review and meta-analysis of eight studies including 439 participants, vitamin D supplementation was associated with improved symptoms on the OSDI, with a standardized mean difference (SMD) of −1.10 (95% confidence interval −1.45 to −0.74), and improved clinical signs, including increased tear production on the Schirmer tear secretion test (SMD 1.43, 95% confidence interval 0.81 to 2.05) and improved tear film stability measured by tear film break-up time (SMD 1.19, 95% confidence interval 0.83 to 1.55), supporting a potential therapeutic role for supplementation. Similar findings were noted in smaller studies. In a retrospective observational study, 105 patients with DED refractory to conventional treatment and vitamin D deficiency or insufficiency (mean serum 25-hydroxyvitamin D level 10.52 ± 4.61 ng/mL) received a single intramuscular injection of 200,000 international units of cholecalciferol. Several tear film parameters improved, at least transiently, including tear stability (3.16 ± 2.27 → 5.19 ± 2.34 seconds at 6 weeks, p =.001) and production (6.69 ± 3.92 → 8.92 ± 7.60 millimeters at 6 weeks ( p =.015). Taken together, these findings suggest that vitamin D supplementation may serve as a useful adjunct therapy for DED, but larger, high-quality randomized trials are still needed to define optimal dosing, durability of effect, and the patient subgroups most likely to benefit.
Overall, although significant associations between vitamin D deficiency and DED have been reported, the literature is limited by small sample sizes, cross-sectional designs, and methodological inconsistencies across studies. Large-scale, longitudinal analyses are thus needed to clarify the strength and direction of this relationship. This study aims to address these gaps through a retrospective cohort study using a large, federated database to evaluate the association between vitamin D deficiency and incident DED.
METHODS
This retrospective cohort study utilized data from the TriNetX US Collaborative Network, a federated research platform that aggregates de-identified electronic health records from more than 180 healthcare organizations across the United States. The platform is fully compliant with the Health Insurance Portability and Accountability Act (HIPAA), with all patient-level data de-identified in accordance with Section 164.514(a) of the HIPAA Privacy Rule. TriNetX’s integrated analytic tools, including incidence estimation, outcomes analyses, and propensity score matching, allow patient-level analyses while reporting only aggregate-level results, thereby preserving data confidentiality. Because this study relied exclusively on de-identified data, it did not constitute human subjects research, and review by an Institutional Review Board was not required. The study adhered to the ethical principles outlined in the Declaration of Helsinki and its subsequent amendments and was conducted in accordance with the STROBE (STrengthening the Reporting of Observational Studies in Epidemiology) guidelines.
The TriNetX database was queried for records spanning September 2005 through September 2025. Eligible participants were adults aged 18 years or older at the time of their most recent qualifying clinical encounter. Two cohorts were defined using International Classification of Diseases, Tenth Revision (ICD-10) codes. Within the TriNetX platform, historical ICD-9 diagnoses are algorithmically mapped to their ICD-10-CM equivalents, allowing cohort definition using ICD-10 codes over the entire study period. Cohort 1 consisted of patients with a documented diagnosis of Vitamin D deficiency (ICD-10: E55), while Cohort 2 comprised patients who had undergone routine adult medical examinations without abnormal findings (ICD-10: Z00.00) and who had no prior history of vitamin D deficiency. The index date was defined as the first recorded encounter for vitamin D deficiency in Cohort 1 and the first documented routine adult examination in Cohort 2. Encounters occurring before September 2005 were excluded. Outcomes were assessed beginning one day after the index event, and all subsequent events recorded through September 2025 were included in the observation period.
To reduce baseline imbalances and minimize confounding, propensity score matching was performed at a 1:1 ratio using the nearest-neighbor method without replacement. Matching was restricted to 5 prespecified variables: age at index and sex in the demographics category; ischemic heart disease (ICD-10: I20–I25) and diabetes mellitus (ICD-10: E08–E13) in the diagnosis category; and body mass index (BMI) in the laboratory category. Other characteristics were captured and reported in the analysis, including race and ethnicity, as well as clinical diagnoses relevant to ocular health such as menopausal and perimenopausal disorders (ICD-10: N95), essential (primary) hypertension (ICD-10: I10), nicotine dependence (ICD-10: F17), other hypothyroidism (ICD-10: E03), thyrotoxicosis/hyperthyroidism (ICD-10: E05), depressive episode (ICD-10: F32), and generalized anxiety disorder (ICD-10: F41.1). Covariate balance after propensity score matching was assessed using SMDs, with an absolute SMD <0.1 considered indicative of ideal balance and values <0.2 regarded as acceptable for residual imbalance. Follow-up duration for both cohorts was reported as the mean (± standard deviation) and the median (with interquartile range), expressed in days.
The primary outcome was the diagnosis of Dry Eye Disease (ICD10: H04.12). Statistical analyses were conducted using the TriNetX analytics platform. Hazard ratios (HRs) with 95% CIs were estimated using Cox proportional hazards regression models, and time-to-event outcomes were summarized with Kaplan–Meier cumulative incidence curves and compared using log-rank test. Differences in incidence proportions between cohorts were assessed with chi-square tests, and effect sizes were expressed as risk ratios (RRs) and odds ratios (ORs) with corresponding 95% confidence intervals (CIs). A 2-sided P value <.05 was considered statistically significant.
RESULTS
Cohort Characteristics
Cohort selection is summarized in Figure 1 . Within the TriNetX US Collaborative Network, Cohort 1 included adult patients (≥18 years old) with a documented diagnosis of vitamin D deficiency (ICD-10: E55), totaling 6,634,934 individuals prior to propensity score matching. Cohort 2 included 11,069,982 adults with no history of vitamin D deficiency who had undergone routine adult medical examinations without abnormal findings (ICD-10: Z00.00). To align the study population with the 20-year observation window, 22,595 patients in Cohort 1 and 276,769 patients in Cohort 2 whose index event occurred before September 2005 were excluded. After propensity score matching, 6,047,502 patients were included in each cohort, yielding a combined matched study population of 12,095,004 individuals. To restrict analyses to incident DED (ICD-10: H04.12) and preserve temporal ordering between exposure and outcome, patients with a prior diagnosis of DED on or before the index date were excluded from the outcome analysis, resulting in the removal of 156,098 patients from Cohort 1 and 88,290 patients from Cohort 2. The final analytic sample, therefore, comprised 5,891,404 patients in the vitamin D–deficient cohort and 5,959,212 patients in the control cohort (total N =11,850,616).
Cohort selection flowchart. The final analytic sample comprised 5,891,404 patients in the vitamin D–deficient cohort and 5,959,212 patients in the control cohort (total N = 11,850,616). Abbreviation: BMI: Body Mass Index.
Baseline characteristics were well balanced between the cohorts, with standardized mean differences (SMDs) <0.1 for most covariates; residual imbalances were observed for Black race (SMD = 0.122), BMI (SMD = 0.158), and not Hispanic or Latino ethnicity (SMD = 0.196), all of which remained below the acceptable 0.2 threshold, whereas unknown ethnicity showed a modest residual imbalance slightly above this cutoff (SMD = 0.238) ( Table 1 ). In the post-match cohorts, the mean follow-up times were comparable, with 1,647 days (SD ±1,417) in the vitamin D deficiency group and 1,690 days (SD ±1,504) in controls. The corresponding median follow-up times were 1,261 days (IQR 2,059) and 1,294 days (IQR 2,201), respectively.
Table 1
Pre and Post Matching Cohort Characteristics
| Characteristic | Vitamin D Deficiency (Pre-Match) N = 6,634,934 | No Deficiency (Pre-Match) N = 11,069,982 | SMD | Vitamin D Deficiency (Post-Match) N = 6,047,502 | No Deficiency (Post-Match) N = 6,047,502 | SMD |
|---|---|---|---|---|---|---|
| Age at Index (Mean ± SD) | 51.8 ± 19.9 | 43.8 ± 19.7 | 0.403 | 50.3 ± 19.7 | 50.8 ± 19.5 | 0.025 |
| Female (%) | 4,455,822 (67.4%) | 5,669,254 (52.5%) | 0.307 | 3,941,883 (65.2%) | 3,875,091 (64.1%) | 0.023 |
| Black or African American (%) | 1,033,092 (15.6%) | 1,265,073 (11.7%) | 0.114 | 952,227 (15.7%) | 698,462 (11.5%) | 0.122 |
| White (%) | 4,458,999 (67.4%) | 7,372,492 (68.3%) | 0.019 | 4,047,080 (66.9%) | 4,246,512 (70.2%) | 0.071 |
| American Indian or Alaska Native (%) | 31,661 (0.5%) | 36,074 (0.3%) | 0.023 | 29,268 (0.5%) | 19,248 (0.3%) | 0.026 |
| Native Hawaiian or Other Pacific Islander (%) | 19,026 (0.3%) | 36,567 (0.3%) | 0.009 | 17,380 (0.3%) | 20,185 (0.3%) | 0.008 |
| Asian (%) | 297,824 (4.5%) | 462,698 (4.3%) | 0.011 | 277,524 (4.6%) | 245,936 (4.1%) | 0.026 |
| Other Race (%) | 265,861 (4.0%) | 449,365 (4.2%) | 0.007 | 249,356 (4.1%) | 229,076 (3.8%) | 0.017 |
| Unknown Race (%) | 505,876 (7.7%) | 1,170,944 (10.8%) | 0.111 | 474,667 (7.8%) | 588,083 (9.7%) | 0.066 |
| Not Hispanic or Latino (%) | 4,761,297 (72.0%) | 6,623,357 (61.4%) | 0.227 | 4,341,180 (71.8%) | 3,786,885 (62.6%) | 0.196 |
| Hispanic or Latino (%) | 541,872 (8.2%) | 860,643 (8.0%) | 0.008 | 507,095 (8.4%) | 441,629 (7.3%) | 0.040 |
| Unknown Ethnicity (%) | 1,309,170 (19.8%) | 3,309,213 (30.7%) | 0.252 | 1,199,227 (19.8%) | 1,818,988 (30.1%) | 0.238 |
| Menopausal and other perimenopausal disorders (%) | 25,274 (0.4%) | 15,888 (0.1%) | 0.046 | 19,314 (0.3%) | 15,316 (0.3%) | 0.012 |
| Ischemic heart diseases (%) | 119,543 (1.8%) | 102,934 (1.0%) | 0.073 | 93,545 (1.5%) | 96,676 (1.6%) | 0.004 |
| Essential (primary) hypertension (%) | 558,985 (8.5%) | 438,827 (4.1%) | 0.182 | 443,001 (7.3%) | 376,709 (6.2%) | 0.044 |
| Diabetes mellitus (%) | 327,451 (5.0%) | 226,341 (2.1%) | 0.155 | 238,875 (3.9%) | 223,842 (3.7%) | 0.013 |
| Nicotine dependence (%) | 87,591 (1.3%) | 86,184 (0.8%) | 0.051 | 73,909 (1.2%) | 68,219 (1.1%) | 0.009 |
| Hypothyroidism (%) | 154,025 (2.3%) | 94,069 (0.9%) | 0.116 | 122,194 (2.0%) | 83,195 (1.4%) | 0.050 |
| Thyrotoxicosis [hyperthyroidism] (%) | 18,499 (0.3%) | 8,140 (0.1%) | 0.049 | 15,477 (0.3%) | 6,750 (0.1%) | 0.034 |
| Depressive episode (%) | 156,749 (2.4%) | 99,211 (0.9%) | 0.114 | 131,633 (2.2%) | 76,745 (1.3%) | 0.070 |
| Generalized anxiety disorder (%) | 59,706 (0.9%) | 45,444 (0.4%) | 0.059 | 52,364 (0.9%) | 30,653 (0.5%) | 0.043 |
| BMI | 30.5 ± 8.4 | 29.0 ± 7.3 | 0.183 | 30.4 ± 8.5 | 29.1 ± 7.3 | 0.158 |
Stay updated, free articles. Join our Telegram channel
Full access? Get Clinical Tree