Purpose. To study the geographical distribution of corneal temperature (CT) and its influence on the intraocular pressure (IOP) of healthy human volunteers. Materials and Methods. Fifteen subjects (7 M, 8 F), 33.8±17.4 years old, were enrolled in this pilot, cross-sectional study. Measurements of CT were taken after one hour with closed eyelids (CET) or closed eyelids with a cooling mask (cm-CET) and compared to baseline. Results. If compared to baseline, after CET, average CT significantly increased by 0.56°C in the RE and by 0.48°C in the LE (p<0.001) and IOP concomitantly significantly increased by 1.13 mmHg and 1.46 mmHg, respectively, in each eye (p<0.001). After cm-CET, average CT significantly decreased by 0.11°C and 0.20°C, respectively, in the RE and LE (RE p=0.04; LE p=0.024), followed by a significant IOP decrease of 2.19 mmHg and 1.54 mmHg, respectively, in each eye (RE p<0.001; LE p=0.0019). Conclusion. Significant variations of CT occurred after CET and cm-CET and were directly correlated with significant differences of IOP. It can be speculated that both oxidative stress and sympathetic nerve fiber stimulation by temperature oscillations may affect the regulation of AH vortex flow and turnover, thus influencing IOP values.
1. Introduction
Temperature is one of the fundamental regulators of tissue metabolism [1, 2]. Interest in the temperature of the eye spans almost 130 years and the ability to measure the temperature of the eye, driven by prevailing technologies, has potential importance in both research and clinical situations, including the study of ocular physiology and pathology [3–9].
New infrared ocular thermographs allow a noncontact and nonintrusive characterization of the thermal profile across the ocular surface [10–17]. Applications have included dry eye, wearing contact lens, corneal sensitivity, and ophthalmic surgery [18–24]. In addition, some studies showed a correlation between ocular surface temperature and ocular blood flow. An increase of intraocular pressure (IOP) was found to be related to a contemporary decrease of ocular perfusion pressure and ocular temperature in monkeys [25]. A recent study on humans has also reported that eyes with ischemic central venous retinal occlusion (CRVO) have lower ocular surface temperatures than nonischemic ones [26]. Moreover, in carotid artery stenosis, the eye on the affected side has been found to have an impairment in retrobulbar hemodynamics along with a reduction in corneal temperature (CT) [26]. Thermography has also been applied to explore the role of vascular factors in the physiopathology of glaucoma and Galassi et al. have recently defined ocular surface temperature as a marker of impaired retrobulbar hemodynamics in patients with glaucoma [27].
However, to date, little work has been undertaken to determine the relationship between IOP and CT [3, 22, 28]. The variations in IOP following closed eyelid test (CET) both under normal conditions and after the administration of antioxidants have recently been investigated, leading to the conclusion that CET-induced ocular hypertension could be a response to mixed stress—oxidative and thermic—with degenerative effects on the trabecular meshwork (TM) [22, 29]. Moreover, given the known influence of ambient pressure and temperature on IOP, an underwater mask has been proposed as a provocative test in the diagnosis of primary open angle glaucoma (POAG) [30].
POAG is a multifactorial and not yet well-understood pathology [31, 32]. Ocular hypertension is one known critical risk factor for glaucomatous optic neuropathy. IOP is generated and maintained via the aqueous humor circulation system in the anterior chamber (AC) of the eye. The major factor controlling IOP is the dynamic balance between aqueous humor production in the ciliary body and its draining through so-called conventional—TM and Schlemm’s canal (SC)—and uveoscleral outflow pathways [33–36].
The goal of the study here reported was to map the CT in different areas of the cornea of healthy human volunteers and follow its variations after CET at room temperature or with a cooling mask (cm-CET) and correlate such variations to IOP values. The application of a physical model was then used to validate the experimental results.
2. Patients and Methods
This pilot, prospective, cross-sectional study was conducted following the tenets of the Declaration of Helsinki. Fifteen healthy Caucasian human volunteers (8 females and 7 males) were enrolled. The mean age was 33.8±17.4 years (range 18–72). After signing an informed consent, each study participant underwent a complete ophthalmological examination, including a medical history review, best-corrected visual acuity measurement (BCVA), slit-lamp biomicroscopy, dry eye tests, and dilated fundus examination. All subjects were free from ocular and systemic diseases with no history of previous ocular surgery. Patients enrolled in the study were accepted if they had no signs or symptoms of ocular dryness: Ocular Protection Index (OPI: calculated as the ratio between BUT and the blinking frequency per minute) score ≥ 1 and OSDI score ≤ 12. In addition, since corneal pachymetry may influence temperature fluctuations and their determination, patients’ central corneal thickness measured by ultrasonic pachymetry (Pacline compact multifeature pachymeter, Optikon 2000, Rome, Italy) had to be within normality limits, that is, 550±20 microns.
Further inclusion criteria were as follows:
IOP ≤ 21 mmHg, without any treatment.
Refraction values between −4 and +4 spheric diopters.
BCVA for far distance equal to 10/10.
Normal angle structure at gonioscopy.
Normal C/D ratio at slit-lamp examination.
Normal SAP (30-2 SITA standard program).
2.1. Corneal Temperature Measurement
Precise spatiotemporal measurements of CT were captured through the latest generation infrared thermometer (Sola Electro-Optics, Shanghai, China). Since it is known that there is an uneven distribution of temperature in the cornea, the average CT was calculated based on the recordings of CTs in five different areas of the cornea (nasal, temporal, superior, inferior, and central). All patients were acclimatized to the clinical environment for at least 15 min. The room temperature was specifically set and controlled at all times at 25°C (≈77°F), humidity was maintained at 42.0%, and the average indoor illumination was maintained at 300 lux, considered a standard indoor level of illumination. Any air drafts were avoided. The measurements were performed between 9:00 and 11:00 a.m. in a seated position and under the conditions described by Mori et al. [18]: the subject blinked normally, then closed both eyes for 5 s, and then kept the eyes open for more than 10 s. The thermography device was set up 20 cm in front of the eye and the head was held steady with a frame. During that time, the subject was asked not to blink. If the subject blinked, a new measurement was performed.
Three measurements were taken consecutively during a single session for each eye and the average value was recorded. The measurements were repeated after one hour of CET and after the same test following the application of a cooling mask on the volunteers’ eyelids. To avoid any operator related bias, all thermographic measurements were made by one single examiner.
2.2. IOP Measurements
IOP is reported as the mean value of three consecutive readings of each eye by Goldmann applanation tonometer (AT-900, Haag Streit Diagnostics, Switzerland) registered between 9:00 and 11:00 a.m. to minimize the effect of daily variations. To avoid interexaminer and intertonometer variances, all IOP measurements were taken by the same trained resident.
2.3. Physical Model
To simulate the temperature distribution of the human eye, we adopted the physical model developed by Karampatzakis and Samaras [37] that solves the Pennes bioheat transfer equation [38] coupled with the incompressible Navier-Stokes equation of fluid dynamics.
According to such a physical model, the eye can be modeled by seven regions with different thermal properties—the cornea, the anterior chamber, the trabecular meshwork, the iris, the lens, the vitreous humor, and the sclera. The basal metabolic heat generation, as well as the blood perfusion, mainly occurs in the iris and in the sclera. This model also includes secretory inflow, drainage, and circulation of the AH.
2.4. Statistical Analysis
The results obtained were statistically analyzed by Student’s t-test for paired samples. p values below 0.05 were considered significant and denoted by one (p∗<0.05) or two (p∗∗<0.01) asterisks (SPSS V.19, IBM SPSS Statistics, USA).
3. Results
We analyzed 30 eyes of 15 healthy volunteers. Figure 1 shows the distribution of basal values across the cornea. In each eye, the lowest temperature was observed at the temporal side, the highest temperature was observed at the nasal side, and intermediate values were observed along the corneal longitudinal axis (superior, central, and inferior). The same distribution was also maintained after CET and cm-CET.
(a) A graphic illustration of the distribution of temperature on the corneal surface. (b) and (c) Values of the temperature measured in the different areas of the cornea for the right (b) and left (c) eye with open eyelids (basal); after one hour with closed eyelids (CET) and after one hour with closed eyelids wearing a cooling mask (cm-CET).
CT mean basal values were 36.33±0.33°C in the right eye (RE) and 36.32±0.24°C in the left eye (LE). After CET, all the temperatures tended to increase with respect to the basal values (mean values: 36.89±0.20°C in the RE and 36.80±0.26°C in the LE) whereas after cm-CET all the temperatures showed a tendency to decrease (mean values: 36.22±0.39°C in the RE and 36.12±0.23°C in the LE). On average, after CET a highly significant increase of the corneal temperature of 0.56°C for the RE and 0.48°C for the LE (p<0.001 for both eyes) was observed and after cm-CET there was a significant decrease of 0.11°C for the RE and 0.20°C in the LE (RE p=0.04 and LE p=0.024 for both eyes) (Figure 2(a)).
Variation of the average corneal temperature (a) and the IOP (b) in the right (RE) and left (LE) eye with open eyelids (b); after one hour with closed eyelids (CET) and after one hour with closed eyelids wearing a cooling mask (cm-CET). p∗<0.05; p∗∗<0.01.
Correspondingly, we observed significant variations in IOP, which increased after CET by 1.13 mmHg in the RE and 1.46 mmHg in the LE (basal mean values: 12.20±1.72 mmHg in the RE and 12.87±3.6 mmHg in the LE; after CET mean values: 13.33±1.2 mmHg in the RE and 14.33±3.8 mmHg in the LE; p<0.001 for both eyes), whereas it significantly decreased after cm-CET by 2.19 mmHg in the RE and by 1.54 mmHg in the LE (basal mean values: 12.20±1.72 mmHg in the RE and 12.87±3.6 mmHg in the LE; after cm-CET mean values: 10.01±1.78 mmHg in the RE and 11.33±2.1 mmHg in the LE; RE p<0.001 and LE p=0.0019) (Figure 2(b)).
Figure 3 shows that there is indeed a direct correlation between average CT and IOP, with a correlation coefficient scoring 0.74 in the RE and 0.94 in the LE.
Graph of the relationship between average corneal temperature (y-axis) and IOP (x-axis) for the right eye (RE: closed, black circles) and left eye (LE: closed, grey circles). The correlation coefficient for a linear correlation is reported.
Finally, the application of the physical simulation model by Karampatzakis and Samaras [37] and Pennes bioheat transfer equation [38] (Figure 4) shows the following:
The mean temperature of the eye is influenced by the ambient temperature: it increases by about 1.0°C when the ambient temperature rises from 15°C to 30°C. The data that we obtained at an environmental temperature of 25°C were 36.3°C with open eyelids and 36.8°C after CET, in agreement with the theoretical prediction.
The evaporation rate of the thin tear film on the cornea decreases its temperature, similarly to sweat evaporation on the skin; therefore, blocking the evaporation increases the overall surface temperature, as less heat is being exchanged between the cornea and the environment. Mathematical simulations estimate such a temperature rise at about 0.4°C, in agreement with the experimental data after CET, when the evaporation was stopped for some time, and consequently the CT increased by 0.56°C (RE) and 0.48°C (LE).
Theoretical temperature values extrapolated from Karampatzakis’s physical model, relating the ambient temperature (x-axis) to the average corneal surface temperature (y-axis) with open eyelids (closed, grey circles, and dotted line), or after one hour of CET (closed, black circles, and continuous line). The closed squares represent the actual measured values at an ambient temperature of 25°C of averaged corneal temperatures for the right and left eyes with open eyelids (lower grey square), or after one hour of CET (upper black square).
4. Discussion
This study shows that, in normal eyes, following CET and cm-CET, IOP significantly changes in response to variations of CT. In agreement with previous studies [34–37], we also observed an overall increase in CT after CET. In addition, our study is the first to report an overall decrease of CT after application of a cooling mask on the lid surface (cm-CET). Moreover, a direct correlation was found between CT and IOP, with IOP values following the increase or the decrease of average surface CT values. To date, the physiological link between CT and IOP still remains somewhat inconclusive. Previous investigations demonstrated that eyelid closure, or the use of an underwater mask with open eyelids, can raise IOP because of an increase of local temperature [30]. Further studies confirmed the relationship between local temperature, variation in IOP, and anterior chamber oxidative stress following CET [29, 39–41].
As previously reported, our data confirm that CT varies in response to tear film evaporation rate, which in turn is influenced by environmental temperature and blinking rate [19, 42–44]. However, our results also show a peculiar pattern of CT distribution across the corneal surface, not in line with previous findings [12, 14, 15, 45]. In fact, in all experimental settings (basal, after CET, and cm-CET), we detected the coolest area at the corneal temporal region while the warmest area was at the nasal side. Along the corneal vertical axis, the temperatures showed intermediate values. To explain and validate our data, the physical model by Karampatzakis and Samaras [37] and Pennes bioheat transfer equation [38] was applied and allowed us to conclude that AH flow depends on CT distribution. Numerical simulations show that the circulation of AH follows the temperature difference between the temporal and nasal side of the cornea and that an increase of AH flow from stagnation to fast vortices correlates with the increasing asymmetry between temporal and nasal temperatures, thus influencing the balance between production and outflow, and finally IOP values. Accordingly, CT differences between nasal and temporal sides were higher after cm-CET (0.80 and 0.96 for the RE and LE, resp.) than after CET (0.52 and 0.44 for the RE and LE, resp.). This may indeed suggest a faster flow with lower temperatures, favoring a higher discharge of AH, and therefore a lower IOP.
IOP depends on the dynamics of AH turnover, which in turn is a balance between secretion and excretion. Three mechanisms are involved in AH formation by the ciliary body: diffusion, ultrafiltration, and active secretion, the last being the major contributor to its formation. AH can be considered analogous to a blood surrogate as it provides nutrition, removes excretory products from metabolism, transports neurotransmitters, stabilizes the ocular structure, and contributes to the regulation of the homeostasis of the surrounding ocular tissues [33–36]. AH, as part of a vascular circulatory loop, returns to the blood flow by passing through the TM, a uniquely modified vessel wall interposed between the anterior chamber and SC or through the uveoscleral pathway [34, 36].
Our results suggest that temperature may affect AH balance through the regulation of its secretion, excretion, and flow dynamics. Secretion is enhanced by increasing temperatures, because of increased blood flow in the anterior segment of the eye due to vasodilation and upregulation of metabolic processes in the ciliary body. The opposite effect is therefore expected after decreasing the temperature, since it is known that low temperatures stimulate the sympathetic system inducing vasoconstriction. Excretion is also expected to be regulated by temperature. At high temperatures, there is an increase of oxidative stress and inflammatory processes due to increased metabolic activities. These events may negatively affect the draining through the TM, since it has been shown that high temperature increases the cellularity of the anterior chamber angle structures [33, 34]. The opposite is expected with low temperatures.
The application of the physical model by Karampatzakis and Samaras [37] confirmed the direct correlation between IOP and CT and the influence of tear film evaporation and environmental temperature on CT. Moreover, and more interestingly, such numerical simulation suggested that AH responds to increasing surface temperature generating vortices that might contribute to the variations in IOP observed in our experiments.
Finally, the presence of several transient receptor potential (TRP) channel isotypes that are responsive to temperature variations in the three corneal layers and in corneal nerve fibers [46] may suggest an involvement of this receptor type also in controlling the aqueous humor dynamics, although no reports exist yet indicating such kind of interaction.
IOP is the only treatable risk factor in glaucoma, one of the world’s leading causes of blindness. The existing correlation between CT and IOP suggests that CT may affect short-term IOP control and its fluctuations. Such IOP fluctuations are related to faster glaucoma progression. Therefore, controlling IOP and its fluctuations is a main therapeutic goal in glaucoma treatment. Temperature and oxidative stress remain additional targets in the control of glaucoma progression [29, 31, 34, 36, 40]. It will be interesting to extend these findings studying the influence of temperature on IOP not only in the seated position but also in the supine position and on a larger population. Moreover, it could be interesting to see whether the same correlation persists in the presence of ocular hypertension. In such case, the use of a cooling mask during the night, when peaks of IOP are mostly expected with negative effects on POAG progression, might be advisable.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgment
English proofreading by Dr. Antony Bridgewood is greatly appreciated.
RingE. F.Progress in the measurement of human body temperature1998174192410.1109/51.687959JonesB. F.A reappraisal of the use of infrared thermal image analysis in medicine19981761019102710.1109/42.7466352-s2.0-0032216616PurslowC.WolffsohnJ.The relation between physical properties of the anterior eye and ocular surface temperature200784319720110.1097/opx.0b013e3180339f6e2-s2.0-34247607259SniegowskiM.ErlangerM.Velez-MontoyaR.OlsonJ. L.Difference in ocular surface temperature by infrared thermography in phakic and pseudophakic patients2015946146610.2147/OPTH.S69670AlioJ.PadronM.Normal variations in the thermographic pattern of the orbito-ocular region198251293982-s2.0-0020042859AlioJ.PadronM.Influence of age on the temperature of the anterior segment of the eye. Measurements by infrared thermometry198214315315910.1159/0002651872-s2.0-0020059278ShihS.-R.LiH.-Y.HsiaoY.-L.ChangT.-C.The application of temperature measurement of the eyes by digital infrared thermal imaging as a prognostic factor of methylprednisolone pulse therapy for Graves' ophthalmopathy2010885e154e15910.1111/j.1755-3768.2010.01941.x2-s2.0-77954634366ChangT.-C.HsiaoY.-L.LiaoS.-L.Application of digital infrared thermal imaging in determining inflammatory state and follow-up effect of methylprednisolone pulse therapy in patients with Graves' ophthalmopathy20082461454910.1007/s00417-007-0643-02-s2.0-36448953466AcharyaU. R.NgE. Y. K.YeeG. C.HuaT. J.KagathiM.Analysis of normal human eye with different age groups using infrared images200933320721310.1007/s10916-008-9181-52-s2.0-63149097867PurslowC.WolffsohnJ. S.Ocular surface temperature: a review200531311712310.1097/01.icl.0000141921.80061.172-s2.0-21244451409MorganP. B.SohM. P.EfronN.Corneal surface temperature decreases with age1999221111310.1016/S1367-0484(99)80025-32-s2.0-0006820928EfronN.YoungG.BrennanN. A.Ocular surface temperature1989899019062-s2.0-0024467803KlamannM. K. J.MaierA.-K. B.GonnermannJ.KleinJ. P.PleyerU.Measurement of dynamic ocular surface temperature in healthy subjects using a new thermography device201237867868310.3109/02713683.2012.6746102-s2.0-84863687392PattmöllerM.WangJ.PattmöllerJ.ZemovaE.EppigT.SeitzB.SzentmáryN.LangenbucherA.Interobserver and intraobserver reliability of corneal surface temperature measurements with the TG-1000 thermograph in normal eyes2015112974675110.1007/s00347-014-3210-12-s2.0-84922343312MoussaS.EppigT.PattmöllerJ.ZemovaE.SeitzB.LangenbucherA.SzentmáryN.Diurnal and zonal analysis of corneal surface temperature in young healthy adults201323564164610.5301/ejo.50002772-s2.0-84881528160MapstoneR.Measurement of corneal temperature196872237-IN2910.1016/s0014-4835(68)80073-92-s2.0-0014273662TanL.CaiZ.-Q.LaiN.-S.Accuracy and sensitivity of the dynamic ocular thermography and inter-subjects ocular surface temperature (OST) in Chinese young adults2009322788310.1016/j.clae.2008.09.0032-s2.0-60749108260MoriA.OguchiY.OkusawaY.OnoM.FujishimaH.TsubotaK.Use of high-speed, high-resolution thermography to evaluate the tear film layer1997124672973510.1016/s0002-9394(14)71689-72-s2.0-0030831118SuT. Y.HoW. T.LuC. Y.ChangS. W.ChiangH. K.Correlations among ocular surface temperature difference value, the tear meniscus height, Schirmer's test and fluorescein tear film break up time201599448248710.1136/bjophthalmol-2014-305183BuschschlüterS.KochC.von EickenJ.HöhH.Computation of the temperature rise at the corneal endothelium during cataract surgery by modeling of heat generation inside the anterior chamber201440102431244410.1016/j.ultrasmedbio.2014.05.0172-s2.0-84905636138MencucciR.MazzottaC.RossiF.PonchiettiC.PiniR.BaiocchiS.CaporossiA.MenchiniU.Riboflavin and ultraviolet A collagen crosslinking: in vivo thermographic analysis of the corneal surface20073361005100810.1016/j.jcrs.2007.03.0212-s2.0-34249063730GalassiF.GiambeneB.CorviA.FalaschiG.MenchiniU.Retrobulbar hemodynamics and corneal surface temperature in glaucoma surgery200828639940510.1007/s10792-007-9160-82-s2.0-57349197946SinghG.BhinderH. S.Comparison of noncontact infrared and remote sensor thermometry in normal and dry eye patients20051566686732-s2.0-30444444251MartinD. K.FattI.The presence of a contact lens induces a very small increase in the anterior corneal surface temperature19866455125182-s2.0-0022499424AukerC. R.ParverL. M.DoyleT.CarpenterD. O.Choroidal blood flow. I. Ocular tissue temperature as a measure of flow198210081323132610.1001/archopht.1982.010300403010202-s2.0-0019971158SodiA.GiambeneB.FalaschiG.CaputoR.InnocentiB.CorviA.MenchiniU.Ocular surface temperature in central retinal vein occlusion: preliminary data20071757557592-s2.0-38049017568GalassiF.GiambeneB.CorviA.FalaschiG.Evaluation of ocular surface temperature and retrobulbar haemodynamics by infrared thermography and colour Doppler imaging in patients with glaucoma200791787888110.1136/bjo.2007.1143972-s2.0-34347336544BillA.Autonomic nervous control of uveal blood flow19625670812-s2.0-0003035437PescosolidoN.MalagolaR.ScarsellaG.LenarduzziF.DapotoL.NebbiosoM.Oxidative stress in the closed-eyelid test: management of glaucoma20121610145314572-s2.0-84871884971PescosolidoN.La maschera subacquea e' veramente innocua? Effetti sulla pressione endoculare ed utilita' come test di provocazione per la diagnosi nel glaucoma cronico semplice19811079659666SaccàS. C.IzzottiA.RossiP.TraversoC.Glaucomatous outflow pathway and oxidative stress200784338939910.1016/j.exer.2006.10.0082-s2.0-33846810059QuigleyH.BromanA. T.The number of people with glaucoma worldwide in 2010 and 2020200690326226710.1136/bjo.2005.0812242-s2.0-33644655886SaccàS. C.CentofantiM.IzzottiA.New proteins as vascular biomarkers in primary open angle glaucomatous aqueous humor20125374242425310.1167/iovs.11-89022-s2.0-84866284869SaccàS. C.PullieroA.IzzottiA.The dysfunction of the trabecular meshwork during glaucoma course2015230351052510.1002/jcp.248262-s2.0-84920901897IzzottiA.BagnisA.SaccàS. C.The role of oxidative stress in glaucoma2006612210511410.1016/j.mrrev.2005.11.0012-s2.0-32344448577BagnisA.IzzottiA.CentofantiM.SaccàS. C.Aqueous humor oxidative stress proteomic levels in primary open angle glaucoma2012103556210.1016/j.exer.2012.07.0112-s2.0-84866098967KarampatzakisA.SamarasT.Numerical model of heat transfer in the human eye with consideration of fluid dynamics of the aqueous humour201055195653566510.1088/0031-9155/55/19/0032-s2.0-78149351847PennesH. H.Analysis of skin, muscle and brachial arterial blood temperatures in the resting normal human forearm19482153, article 3542-s2.0-77954247486PescosolidoN.BelcaroG.RuscianoD.SteigerwaltR. D.Jr.NebbiosoM.Retrospective study of glaucoma and closed-eyelid test: long-term outcomes in an Italian native population20135517PescosolidoN.‘The oxydative-stress’ role of free radicals in oxygen in hypertension which follow the ‘closed eyelid test’199211888678762-s2.0-0026672716BucciM. G.PescosolidoN.The closed eyelid test in the management of glaucoma1983528489FreemanR. D.FattI.Environmental influences on ocular temperature19731285966022-s2.0-0015906415TanJ.-H.NgE. Y. K.AcharyaU. R.CheeC.Infrared thermography on ocular surface temperature: a review20095249710810.1016/j.infrared.2009.05.0022-s2.0-67949111049MapstoneR.Determinants of corneal temperature1968521072974110.1136/bjo.52.10.7292-s2.0-0014346407NgE. Y. K.OoiE. H.Ocular surface temperature: a 3D FEM prediction using bioheat equation200737682983510.1016/j.compbiomed.2006.08.0232-s2.0-34047120180ReinachP. S.ChenW.MerglerS.Polymodal roles of transient receptor potential channels in the control of ocular function201521, article 51010.1186/s40662-015-0016-4