Thyroid gland dysfunction and vitamin D receptor gene polymorphism in keratoconus

Objectives To detect the serum level of thyroid hormones, vitamin D and vitamin D receptors (VDR) polymorphism in keratoconus (KC) patients and to identify the association between vitamin D deficiency and thyroid dysfunction in KC. Methods This cross sectional study included 177 KC patients with no thyroid disorders compared to 85 healthy controls with normal corneal tomography. Measurements of thyroid stimulating hormone (TSH), free triiodothyronine (FT3), free tetraiodothyronine (FT4) and serum 25-OH vitamin D were done using Enzyme linked immusoassay (ELISA test). VDR polymorphisms were tested including [Taq I (rs731236), Apa I (rs7975232) and Bsm I (rs1544410)] using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP). Results An increase in frequency of thyroid disorders (P = 0.04), decrease in serum 25(OH) vitamin D level (P < 0.001), Taq 1 and tt genotype (P < 0.001) were significantly distributed in KC patients. A significantly higher serum 25(OH) vitamin D level was reported in TT genotype, while insufficient level was more common in Tt genotype (P < 0.001). A deficient serum 25(OH) vitamin D level was predominant in tt genotype (P < 0.001). A 95% confidence interval was in TSH (1.603, 2.946), FT4 (24.145, 77.06), hypothyroidism (1.062, 67.63), insufficient (2.936, 11.643) and deficient vitamin D (5.283, 28.704) and all were significant risk factors for KC with (P < 0.05). Conclusions Both thyroid disorders and low vitamin D are potential factors for KC development. Studying VDR at the molecular level provides interesting avenues for future research toward the identification of new KC cases.


INTRODUCTION
Keratoconus (KC) is a progressive ectatic corneal disease, characterized by conical protrusion of the cornea with progressive thinning resulting in myopia, irregular astigmatism, and associated with severe visual impairment [1]. The exact pathophysiology of KC is not fully explained, many theories had been suggested based on immunological, genetic, and environmental factors [2][3][4][5][6]. Various inflammatory mediators (Cytokines) were linked with KC suggesting a pivotal role of inflammation in the pathogenesis, of the disease [7]. The prevalence of thyroid gland dysfunction was reported in KC patients. Thyroxin was found to cause biochemical changes in the corneal stroma as a result of thyroxin-receptor interaction [8][9][10][11]. Vitamin D receptors (VDR) are found in the cornea, lens, and retinal pigment epithelium [12,13]. The effect of vitamin D on KC was evaluated by Shivakumar et al. [14] in their in vitro study which showed that vitamin D enhances VDR and activates autophagic lysosomal clearance in oxidatively damaged human corneal epithelial cells. An existing relationship between vitamin D deficiency and thyroid gland disorders had been reported in the pathogenesis of keratoconus [15,16]. Many studies had attempted to determine whether low vitamin D levels were associated with the presence and severity of KC [12] and whether they are caused by an isolated inflammatory action [5] or due to VDR and thyroid receptor interaction [17,18].
The aim of this study was to detect the serum level of thyroid hormones, serum vitamin D and VDR polymorphism in KC patients and to identify any association between vitamin D deficiency and thyroid disorders in those patients

SUBJECTS AND METHODS
This prospective, observational, cross sectional study was conducted at Mansoura ophthalmic center, faculty of medicine, Mansoura University in the period from March 2021 to September 2021. It included 177 keratoconus patients without any previously diagnosed thyroid disorders in comparison with 85 healthy controls with normal cornea and without any systemic diseases. This sample of individuals was selected from population in Dakahlia Governorate in Delta, Lower Egypt. The control group was selected from candidates of refractive surgery with normal topography that were free from post-Lasik ectasia for at least six months postoperatively. This study adhered to the tenets of Declaration of Helsinki and was approved by Mansoura faculty of medicine Institutional Review Board (code No R.21.01.1157.R1) and was registered on www.clinicaltrials.gov (NCT05073601). Each patient provided a written informed consent.
Diagnosis of keratoconus was done by a cornea specialist depending on slit lamp biomicroscopy findings of localized corneal thinning and ectasia which was confirmed by Pentacam (Oculus Optikgeräte GmbH, Wetzlar, G-ermany) examination. For each patient flat (K1), steep (K2), and maximum (Kmax) simulated keratometric readings, and corneal pachymetry were recorded. Patients who could not provide the informed consent or the necessary samples for any reason, patients who had any ophthalmic pathology other than KC, and those with systemic disorders were excluded from the study.

Sample collection
A 5 ml venous blood was withdrawn from all subjects included in the study; 2 ml on EDTA for genetic study and 3 ml on a plain tube for thyroid hormones and serum 25-(OH) vitamin D measurement. The blood samples for hormonal assessment were centrifuged then serum samples were stored at −20°C till the analysis time
For all PCR reactions, 50 µL containing 10 µg of genomic DNA, 25 pmoL of each primer, 25 µL of Dream Taq Green PCR Master Mix (including 2.5 units of Taq Polymerase), and 18 µL of double distilled water were used. PCR was performed by an initial denaturation cycle at 95°C for 3 min, annealing for 30 cycles at 95°C for 30 s, primer specific temperatures used for 30 s (TaqI 61°C, ApaI 72°C and 58.5°C for BsmI) and a temperature of 72°C for 30 s, and the final step of extension at 72°C for 5 min [21].
The PCR products were then digested by restriction enzymes (Taq I, Apa I, and Bsm I) supplied by Thermo Scientific. The enzyme-specific restriction reactions were performed according to the insert supplied with the kits. Then, the separation of DNA fragments was done using 2% agarose gel electrophoresis and visualized under UV light ( Fig. 1A-C).

Statistical analysis
Sample size is calculated using online sample size calculator available at (https://www.openepi.com/SampleSize/SSMean.htm). The collected data was analyzed using the Statistical Package of Social Science (SPSS) program for Windows (SPSS, Inc, Chicago, IL) version 22. Kolmogorov-Smirnov test was used to test data normality. Quantitative data were presented as median (range) as they were not normally distributed while qualitative data were described using number and percent. Mann-Whitney test was used to compare two medians. Chi-square test or Fisher's exact test were used to test the association between categorical variables. Chi-squared test was used to determine the deviations from Hardy-Weinberg equilibrium (HWE) expectations. A logistic regression model using a backward stepwise method was used for the significant parameters. Odds ratio (OR) and 95% confidence interval (CI) were calculated. The results were considered to be significant when the probability of error is less than 5% (p < 0.05).

RESULTS
This cross-sectional study included 177 keratoconus patients compared to 85 age and sex-matched healthy subjects with normal cornea which were enrolled as a control group. Demographic and topographic parameters of the included subjects showed a significant increase in K1, K2 and K max in keratoconus group compared to controls, and a significant reduction in pachymetry in keratoconus group compared to controls. Laboratory investigations showed a significant increase in TSH level and FT4 (P = 0.001 and <0.001) respectively and the frequency of thyroid disorders was statistically significant higher in KC group. There was a significant decrease in serum 25(OH) vitamin D level and significant increase in vitamin D insufficiency and deficiency in KC group compared to control group (Table 1).
Applying Hardy Weinberg equation revealed that TaqI, ApaI and BsmI genotypes in two groups are in HWE. As regards the distribution of VDR genotypes Taq 1 was significantly associated with KC. The tt gene variant alleles were of higher statistically significant distribution in KC patients when compared to controls indicating that tt genotype has higher risk to develop KC (P = 0.001, OR = 3.459). On the other hand, there was no statistically significant association between the APA 1 or BSM1 and KC ( Table 2).
Comparison of vitamin D level among the studied VDR genotypes in KC patients reported a significantly higher serum 25(OH) vitamin D level in TT genotype, while insufficient level was more common in Tt genotype. In addition, a deficient serum 25(OH) vitamin D level was predominant in tt genotype (Table 3).
Univariate analysis was done to study risk factors for KC development. The following parameters: TSH, T4, Hypothyroidism, Hyperthyroidism, vitamin D level, insufficient and deficient vitamin   D, and tt genotype of Taq1 were found to be statistically significant (P < 0.05). Logistic regression analysis was conducted in KC group, using the above significant covariates, TSH, T4, hypothyroidism, insufficient and deficient vitamin D were predicted to be risk factors for the disease. On the other hand, sufficient vitamin D was a significant protective factor (Table 4).

DISCUSSION
This study aimed to identify any association between vitamin D deficiency and thyroid disorders in Kc patients as both of them are potential factors that were recently linked to KC development [8,22,23].
Regarding thyroid dysfunction, this study reported a statistically significant higher incidence of thyroid dysfunction with more predominance of hypothyroidism in KC patients in comparison to healthy controls. This was in agreement with previous studies which reported that hypothyroidism may play a role in the exaggeration or even initiation of KC [8,[24][25][26]. Others have demonstrated that thyroxin in tear fluid may have a possible role in KC development [24]. The relationship between KC and thyroid disorders is not yet well settled. Karabulut et al. [25] showed that patients with Grave's disease have significantly altered corneal biomechanical properties, while Gatzioufas et al. [11,26] found that hypothyroidism may induce corneal topographical and biomechanical changes or even exacerbate KC resulting in acute corneal hydrops. Further studies are needed to determine whether the increase or the decrease in thyroid hormones levels is more associated with KC development and progression.
Previous studies [15,16] reported the involvement of Vitamin D and thyroid dysfunctions in KC development as they both work on receptor belonging to steroid/thyroid receptor family. The present study revealed a significant decrease in serum vitamin D level in KC patients than the control group. In addition, there was a significant increase in both vitamin D insufficiency and vitamin D deficiency (53.7% and 35.6% in KC group compared to 36.5% and 12.9% in control group respectively). This finding goes in   [18] who demonstrated that inflammatory, apoptotic, or oxidative mechanisms may cause ocular surface affection due to vitamin D deficiency. Moreover, researchers found that decreased vitamin D level significantly increased non progressive KC probability by 1.23 and progressive KC probability by 1.29 more than the control group. They concluded that evaluating serum vitamin D levels in KC patients at the onset and the follow-up examinations may predict the disease course [23].
In this study, we also evaluated the association between the most common VDR polymorphisms: TaqI (rs731236), BsmI (rs1544410) and ApaI (rs7975232) with KC. To the best of our knowledge this is the first research to study this association. It was found that Taq1 polymorphism was associated with an increased risk of KC. The tt genotype had a higher risk to develop KC as t allele was more prevalent in KC patients (56.5%) while T allele was shown to have a protective effect against the disease (60% in control group compared to 43.5% in KC group). On the other hand, APA 1 and BSM1 were silent polymorphisms not associated with an increased risk of KC. Assessment of vitamin D level among studied VDR genotypes in KC patients revealed a significantly higher serum vitamin D level in TT genotype and insufficient serum vitamin D level in Tt genotype, while deficient serum vitamin D level was predominant in tt genotype.
In our study, Logistic regression analysis for prediction of keratoconus using the multivariate analysis revealed that TSH, T4, hypothyroidism, insufficient and deficient vitamin D levels were significantly associated with keratoconus.
One of the Limitations of our study is that we didn't correlate the vitamin D level to keratoconus grade and the state of progression. Future studies are needed to assess this relationship and to confirm the effect of vitamin D supplementation on keratoconus stability. In addition, Assessment of parathyroid hormone is also needed to study its effect on thyroid hormones and vitamin D which may have a direct role in KC development.

SUMMARY TABLE
What was known before • The exact pathophysiology of keratoconus is not well understood. Thyroid disorders and vitamin D may play a role in disease development What this study adds • Thyroid dysfunction and vitamin D deficiency are significantly associated with keratoconus. vitamin D receptor polymorphism is associated with keratoconus development. Taq1 gene and its tt alleles were reported as high risk factors for development of keratoconus.

DATA AVAILABILITY
Derived data supporting the findings of this study are available from the corresponding author on request.