Variability of autonomic nerve activity in dry eye with decreased tear stability

Abstract

The autonomic nervous system plays a crucial role in the maintenance of homeostasis. Neurogenic inflammation due to decreased stability of tear film may induce changes in autonomic nerve activity, which could be associated with symptom expression.This study aimed to measure biological parameters that represent autonomic nerve activity in dry eye (DE)s caused by tear film dysfunction and investigate their relationship with symptom intensity. This prospective, cross-sectional, comparative study evaluated 34 eyes of 34 participants (mean age: 52.5 ± 13.4 years; range: 20–81 years) without keratoconjunctival damage. Nineteen eyes in the DE group showed DE symptoms and tear break-up time (TBUT) of ≤5 seconds (short TBUT DE); the 15 eyes in the non-DE group showed no DE symptoms. Autonomic nerve activity was measured for 10 minutes—starting and ending 5 minutes before and after instilling ophthalmic solution—and evaluated using the low-frequency component (LF) to the high-frequency component (HF) ratio of heart rate variability (autonomic balance). The pre-ophthalmic solution administration LF/HF ratio was not significantly different (P = 0.59) between the two groups, however, the standard deviation of the LF/HF ratio (LF/HF-SD) tended to be higher in the DE group than that in the non-DE group (P = 0.086). The DE symptom intensity was significantly related to LF/HF-SD (P = 0.005), which significantly decreased after ophthalmic solution administration in the DE group (P = 0.04). The large fluctuations in autonomic balance may be key for the understanding of the mechanism underlying DE symptoms.

Introduction

There is a discrepancy between the symptoms and clinical findings of dry eye (DE) [1,2]. DE is classified into aqueous-deficient DE and evaporative DE [3,4]. However, the concept that decreased tear film stability is of essence to DE is gaining prominence, especially in Asia [5]. Tear film stability is determined by several factors, including tear volume, tear quality, keratoconjunctival epithelial condition, and eyelid condition. DE with decreased tear film stability and symptoms, with or without increased tear evaporation rate, was described as DE with short tear film break-up time (TBUT) in a Japanese epidemiological study [6]. This type of DE accounts for 80% of all DE cases in Japan; however, its pathogenesis remains largely unknown [6]. The symptoms of DE include eye pain and eye discomfort, but there is no correlation between the severity of symptoms and clinical test findings in DE [1,2]. This discrepancy is especially evident in short TBUT DE [7–9], which has relatively severe symptoms despite having few corneal epithelium clinical test findings [6–9]. As a result, short TBUT DE is integral in discussions concerning DE symptoms.

The Dry Eye Workshop II identified neuropathy as one of the etiologies of DE disease [10]. Neuropathic pain due to neurogenic inflammation caused by damage or irritation to sensory nerve endings on the ocular surface has been implicated as an important cause of DE symptoms [11–14]. Several studies on corneal sensory nerves related to neuropathic pain have focused on transient receptor potential vanilloid 1 (TRPV1), which is associated with pain sensation, and transient receptor potential melastatin 8 (TRPM8), which is associated with temperature change sensation on the corneal surface [15]. In rodent DE models, Bereiter et al. [16] observed high activity in TRPV1, while Kovács et al. [17] observed high activity in TRPM8. In addition, Masuoka et al. [18] showed that both TRPM8 and TRPV1 are expressed in the cold-sensitive corneal nerve endings of individuals with DE and that TRPV1 activation is involved in pain perception. These findings suggest that the aberrant activity of peripheral neurosensitivity may be linked to DE symptoms.

Additionally, higher brain functions may be involved in the onset and progression of DE symptoms [19–21]. We previously used menthol to investigate cold sensitivity and its association with short TBUT DE and found no difference in cold sensitivity between individuals with short TBUT DE and normal eyes. However, there was a difference in whether the cold stimulus was perceived as pleasant or unpleasant [22]; that is, the qualitative perception of cold stimuli in patients with DE may be different from that of patients with normal eyes, while the quantitative perception of cold stimuli may be the same in both DE and normal patients.

The autonomic nervous system plays a crucial role in the maintenance of homeostasis. In the field of DE, the parasympathetic nervous system, which is activated by abdominal breathing, is important to increases tear volume [23]. According to an animal study on lacrimal innervation, the parasympathetic nervous system is more important than the sympathetic nervous system for tear secretion [24]. Neurogenic inflammation due to decreased stability of tear film may induce changes in autonomic nerve activity, associated with symptom expression. Since parasympathetic activity is relatively reduced under stress, it might manifest as discomfort observed in short TBUT DE. Measurement of autonomic nerve activity may be a useful index for objective evaluation of unpleasant emotions associated with DE. In this study, we aimed to investigate the characteristics of autonomic nerve activity and its relationship with symptom intensity in short TBUT DE.

Materials and methods

The study protocol was reviewed and approved by the Ethics Committee of the Institutional Review Board of Ito Clinic, Saitama, Japan (registration number IRIN2021-11). All procedures were performed according to the ethical standards of the responsible committee on human experimentation (institutional and national) and the Helsinki Declaration of 1964, as revised in 2013. Informed consent was obtained from all participants. This study was registered in the University Hospital Medical Information Network (registration number UMIN000045019).

Participants

This prospective, cross-sectional, comparative clinical study evaluated the eyes of keratoconjunctival damage-free volunteers (vital staining score <3) at Wada Eye Clinic. Participants with DE symptoms and TBUTs of ≤5 seconds (short TBUT DE) were designated as the DE group, while those without DE symptoms were designated as the non-DE group (the TBUT value was not considered a criterion of non-DE) by an ophthalmologist. The Schirmer test values were not considered for the diagnosis of short TBUT DE. Exclusion criteria were previous ocular surgery and/or trauma within the last 12 months; anatomic abnormalities in the cornea, conjunctivitis, and/or eyelids; glaucoma; current contact lens use; and systemic diseases that may affect autonomic nerve activity such as heart disease excluding hypertension, diabetes mellitus, Parkinson’s disease, pulmonary disease, and Sjögren’s syndrome.

Sample size determination

The sample size was determined based on a previous study [25]. With a mean difference of 0.51 between the two groups, a standard deviation of 0.5, an alpha error of 0.05, and a power of 0.8, the required sample size was 16 in each group. Forty cases were registered for consideration as ineligible cases.

DE questionnaire

A Japanese version of the Ocular Surface Disease Index (J-OSDI) [26] questionnaire was completed by participants. The J-OSDI, which was translated and culturally adapted from the Ocular Surface Disease Index into Japanese, has good internal consistency, test-retest reliability, and discriminant validity by known-group (DE vs. non-DE) comparisons. J-OSDI scores of ≥13 points were regarded as positive for DE symptoms and those of <13 points as negative.

Stress check questionnaire

Participants were then given the Brief Job Stress Questionnaire (BJSQ) as a second questionnaire, which includes questions concerning participants’ health during the past month [27]. The BJQS consists of 29 questions used to asses 6 parameters (liveliness, frustration, fatigue, anxiety, depression, and physical complaints). Participants choose one of four responses: “almost never,” “sometimes,” “often,” and “always.” Scores were then calculated using a conversion table, with participants achieving a score of 6–30 points. The lower the score, the higher the stress the individual is experiencing; the higher the score, the lower the stress. A score of ≤17 points was considered to be indicative of in mental and physical health distress in everyday life.

Evaluation of the tear function and ocular surface

We performed ocular surface examinations, including TBUT measurement, keratoconjunctival vital staining, and the Schirmer test, before measuring autonomic nerve activity. The TBUT was measured after administering 2 μL preservative-free 1% sodium fluorescein into the conjunctival sac using a micropipette. Keratoconjunctival epithelial staining was evaluated after TBUT measurement. Overall epithelial damage was scored on a scale of 0–9 points, as described previously [28]. The Schirmer test was performed last during the eye surface evaluations.

Autonomic nerve activity measurement

The autonomic nerve activity was measured > 10 minutes after performing the Schirmer test using the Silmee™ Bar type Lite (TDK, Tokyo, Japan). This biosensor automatically calculates heartbeat intervals, pulse wave intervals, and autonomic nerve activity by measuring and analyzing electrocardiogram pulse wave, acceleration, and skin temperature. Only data for heart rate, RR interval, low-frequency component (LF), high-frequency component (HF), and the LF/HF ratio as autonomic parameters were extracted to avoid excessive amounts of data.

The device was attached approximately 3 cm below the center of both clavicles, a position where the electrocardiographic potential can be measured at a position close to the heart, and the influence of the movement of the upper arm and chest muscles was limited (Fig 1). Autonomic nerve activity was measured with the participant in a rested sitting position and blinking naturally for 10 minutes. We instructed the participants not to sleep, to sit quietly, and to stare blankly into the distance. Measurement was started 5 minutes before ophthalmic solution administration of 0.3% hyaluronic acid (0.1% Hyalein®; Santen Pharmaceutical Co., Ltd., Osaka, Japan) and ended 5 minutes after ophthalmic solution administration. The temperature of the room was maintained at 23°C –25°C during examinations, with 50%–60% humidity.

[Figure omitted. See PDF.]

This biosensor automatically calculates heartbeat intervals, pulse wave intervals, and autonomic nerve activity by measuring and analyzing electrocardiograms, pulse waves, acceleration, and/or skin temperature.

Autonomic nerve function was assessed by frequency analysis of cardiac beat movements (fluctuations in heart rate [RR] intervals) and quantification of sympathetic and parasympathetic activities [29,30]. This technique partitions the total variance (the “power”) of a continuous series of beats into its frequency components, typically identifying two main peaks: low frequency (0.04–0.15 Hz) and high frequency (0.15–0.4 Hz) [31]. The HF peak is widely believed to reflect cardiac parasympathetic nerve activity, while the LF is often assumed to have a dominant sympathetic component [29,30,32]. Based on these assumptions, it was proposed that the LF/HF ratio could be used to quantify the degree of sympathovagal balance [29–31,33,34], Heart rate variability metrics obtained in this study were RR, heart rate interval (RR intervals), and the LF/HF ratio.

Temperature sensation and the comfort level to ophthalmic solution

We evaluated the participants’ temperature sensation and comfort level to the ophthalmic solution at room temperature. The temperature sensation and comfort level perceived immediately after administration of the ophthalmic solution were rated using a visual analog scale: a temperature sensation with a score of −5 indicated “most intense cooling sensation imaginable,” and a score of +5 indicated “most intense warming sensation imaginable”; a comfort level with a score of −5 indicated “most intense unpleasant sensation imaginable,” and a score of +5 indicated “most intense pleasant sensation imaginable.”

Statistical analyses

Data are presented as mean ± standard deviation, where applicable. Baseline J-OSDI scores, tear function, and autonomic nerve activity were compared between both groups using the Mann-Whitney test, since it was difficult to regard the normal distribution as a result of the normality test. Heart rate variability parameters pre- and post-ophthalmic solution administration were compared in each group using the Wilcoxon test. The temperature sensation and the comfort level to the ophthalmic solution were also compared between both groups using the Mann-Whitney test. Multiple regression analysis was performed to determine factors affecting J-OSDI scores, with the pre-ophthalmic solution administration parameters of DE and heart rate variability in the overall cohort as the dependent variables.

All statistical analyses were performed using SPSS Statistics for Windows, version 17.0, (SPSS Inc., Chicago, IL, USA). Statistical significance was set at a P-value of <0.05.

Results

Demographics of the study population

Of the 40 eligible eyes of 40 volunteers with no keratoconjunctival damage (vital staining score <3), 34 eyes of 34 participants (11 males, 23 females; mean age: 52.5 ± 13.4 years; range: 20–81 years) were included in this study (S1 Table). There were six participants with DE symptoms and TBUTs of >5 seconds who were excluded from the study. Table 1 presents the demographic characteristics of the study participants. The DE and non-DE groups comprised 19 eyes of 19 participants (4 males, 15 females; mean age: 51.4 ± 12.5 years; age range: 20–74 years) and 15 eyes of 15 participants (7 males, 8 females; mean age: 53.9 ± 14.7 years; age range: 20–84 years), respectively. Regarding the presence or absence of systemic disease, hypertension was observed in 3 participants each in the DE and non-DE groups, and breast cancer in remission was seen in 1 particpant in the DE group.

[Figure omitted. See PDF.]

The TBUT value was significantly lower (P < 0.05), Schirmer test value tended to be lower (P = 0.05), and the J-OSDI scores were significantly higher in the DE group than those in the non-DE group (P < 0.05). There were no significant differences in the keratoconjunctival staining scores between the two groups (P > 0.05).

Temperature sensation and the comfort level to ophthalmic solution

There were no significant differences in temperature sensitivity to the ophthalmic solution at room temperature (cold stimulus) or stress levels between the DE and non-DE groups (P > 0.05; Table 1).

Autonomic nerve activity

Table 2 shows the results of each heart rate variability metric in the DE and non-DE groups. No significant differences were observed in the pre-ophthalmic solution administration LF, HF, and LF/HF ratio between the DE and non-DE groups (P > 0.05). There were no significant changes in the LF, HF, and LF/HF ratio pre- and post-ophthalmic solution administration in the DE group (P > 0.05; Fig 2A). Similar results were obtained in the non-DE group (P > 0.05).

[Figure omitted. See PDF.]

(A) The LF/HF ratio pre- and post-instillation of ophthalmic solution in the DE and non-DE groups. There are no significant changes in either the DE or non-DE groups (P > 0.05). (B) The standard deviation of the LF/HF ratio pre- and post-instillation of ophthalmic solution in the DE and non-DE groups. The standard deviation of the LF/HF ratio significantly decreased post-instillation in the DE group (P < 0.05), whereas it was low pre-instillation and increased post-instillation in the non-DE group. (C) Typical DE and non-DE cases show fluctuation of the LF/HF ratio over time. The DE case presented is of a 39-year-old woman with a J-OSDI score of 18.8 points. The non-DE case is of a 21-year-old woman with a J-OSDI score of 9.1 points. The fluctuation in the LF/HF ratio of the DE case is larger than of the non-DE case, and it decreases post-instillation. DE, dry eye; HF, high-frequency component; J-OSDI, Japanese version of the Ocular Surface Disease Index Questionnaire; LF, low-frequency component.

[Figure omitted. See PDF.]

The pre-ophthalmic solution administration standard deviation of the LF/HF ratio (LF/HF-SD), which represents the fluctuation in the LH/FH ratio over time, tended to be higher in the DE group than that in the non-DE group (P = 0.086). It significantly decreased post-ophthalmic solution administration in the DE group (P = 0.04), whereas it increased in the non-DE group (P = 0.05; Fig 2B). Fig 2C shows the fluctuation of the LF/HF ratio over time in typical DE and non-DE cases.

Relationship between autonomic nerve activity and DE symptom intensity

Table 3 shows the results of multiple regression analysis for determining factors affecting the J-OSDI scores with independent variables, such as the Schirmer test values, TBUTs, keratoconjunctival scores, and LF/HF-SD. The LF/HF-SD value was a significant factor affecting the J-OSDI scores (P < 0.05), whereas the other variables were not.

[Figure omitted. See PDF.]

Discussion

Heart rate variability analysis has been widely used to assess autonomic nervous function [25,29–35]. The HF is believed to reflect cardiac parasympathetic nerve activity, while LF mainly reflects sympathetic activity [29–31,33]. The LF/HF ratio has gained wide acceptance as a tool to assess cardiovascular autonomic regulation, where increases in the LF/HF ratio are assumed to reflect a shift to “sympathetic dominance” and decreases in this index to correspond to a “parasympathetic dominance” [31]. Positive and negative emotional states result in differential heart rate variability responses [36–40].

Based on the hypothesis that the symptoms of short TBUT DE may be induced by the central nervous system, which is associated with the autonomic nervous system, we measured autonomic nervous activity and investigated 1) whether short TBUT DE has specific characteristics with regards to autonomic activity and 2) if it is associated with clinical symptoms.

Short TBUT DE was characterized by substantial changes in autonomic balance over time. In addition, DE symptom intensity was related to fluctuations in autonomic balance over time. Accordingly, the onset of DE symptoms may be related to fluctuations in the autonomic balance. We initially expected sympathetic nerve activity to be high under the stress induced by unpleasant emotions caused by eye discomfort and/or pain in DE. However, the subjective stress levels observed in the DE group were not high and were similar to those in the non-DE group. There was no difference in the LF/HF ratio between the DE and non-DE groups, indicating that DE symptoms were not linked to either subjective stress or the magnitude of autonomic nervous activity. Instead, the LF/HF ratio fluctuations over time were larger in the DE group than those in the non-DE group, and DE symptom intensity was related to fluctuations in the LF/HF ratio, while DE parameters such as keratoconjunctival epithelial damage, tear stability, and amount of tear secretion, were not. The result that there was no difference in temperature sensation to a cold stimulus (instillation of eye drops) with or without DE indicates that peripheral nerve sensitivity is not involved in the induction of DE symptoms. It is speculated that changes in the central nervous system, rather than changes in the peripheral nervous system, are involved.

Fluctuations in autonomic activity can play an important role in disease pathogenesis. Living organisms are constantly subjected to changes and stimuli in their external environment and are equipped with mechanisms to maintain optimal functioning through regulating tissues and organs. The regulation of the autonomic nervous system is one such mechanism. However, DE, as a pathological condition, might cause the autonomic nervous system to initiate unstable physiological functioning. Fluctuations in autonomic activity are involved in the etiology of several diseases [41,42]. For example, a comparative study of carotid endarterectomy and carotid artery stenting, which are both treatments for hypertrophic myocardial infarction, showed that long-term autonomic fluctuations due to parasympathetic muscle tone hypertonia occur perioperatively to carotid endarterectomy, and increase the incidence of myocardial infarction [41]. A study of autonomic activity in patients with Brugada syndrome, a heart disease that causes idiopathic arrhythmias (ventricular fibrillation), showed that the circadian variation of the LF and LF/HF ratio were low in patients with the syndrome and that the LF/HF-SD and LF could be useful parameters for the diagnosis of Brugada syndrome [25]. Similarly, our study shows that such fluctuations may be associated with DE symptoms.

Fluctuations in the LH/HF ratio decreased post-ophthalmic solution administration in the DE group, but not in the non-DE group. However, there was no significant difference between the two groups in terms of the comfort level experienced by participants when they were exposed to the ophthalmic solution, suggesting that the LF/HF-SD was unrelated to comfort when exposed to a stimulus. These findings are contradictory. Temporarily improving DE symptoms by ophthalmic solution administration and feeling comfortable with ophthalmic solution administration may not be related.

Our study has a few limitations. First, we did not confirm if the participants had any autonomic dysfunctions other than DE. Menopausal disorders, depression, and some other diseases are known to impair the balance of the autonomic nervous system in a broad sense; this could have been confirmed by performing an autonomic nerve test, such as a head-up tilt test or assessing changes in blood pressure due to changes in posture. Strictly speaking, it cannot be ruled out that these factors may have influenced the current results. However, since mental stress was similar between the DE and non-DE groups, it is unlikely that mental stress affected changes in autonomic nerve activity. Second, participants with decreased TBUT without DE symptoms were included in the non-DE group. Patients with decreased TBUT included those who remained asymptomatic (non-DE) and those who developed symptoms (DE). Thus, in this study, we compared the DE and non-DE groups, focusing on the presence or absence of symptoms. In fact, in Asia, DE diagnosis is defined as a case in which tear film instability and DE symptoms are present [5]. However, tear film instability without DE symptoms is not considered as confirmed DE diagnosis. Third, the sex of the participants in this study did not match between the DE and non-DE groups, possibly influencing the results. In addition, lifestyle habits such as smoking [42] and sleeping [43], which can be confounding factors and may affect the results, were not evaluated. In addition, since the autonomic nervous system is affected by environmental factors (e.g., temperature, humidity, and room illuminance) and/or body conditions (e.g., aging, drowsiness, and fatigue), it was difficult to assess autonomic nerve activity that reflected only DE.

Conclusion

In conclusion, the fluctuation in the autonomic balance was large and related to DE symptom intensity in short TBUT DE. This fluctuation was decreased by instilling hyaluronic acid ophthalmic solution. Fluctuations in autonomic activity may be an important key to understanding the mechanism of DE symptoms.

Supporting information

S1 Table. Raw data.

https://doi.org/10.1371/journal.pone.0276945.s001

Acknowledgments

Ethics

The study protocol was reviewed and approved by the Ethics Committee of the Institutional Review Board of Ito Clinic, Saitama, Japan (registration number IRIN2021-11). All procedures were performed according to the ethical standards of the responsible committee on human experimentation (institutional and national) and the Helsinki Declaration of 1964, as revised in 2013. Written informed consent was obtained from all participants. This study was registered in the University Hospital Medical Information Network (registration number UMIN000045019).

Citation: Kaido M, Arita R, Mitsukura Y, Ishida R, Tsubota K (2022) Variability of autonomic nerve activity in dry eye with decreased tear stability. PLoS ONE 17(11): e0276945. https://doi.org/10.1371/journal.pone.0276945

About the Authors:

Minako Kaido

Roles: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Project administration, Resources, Software, Validation, Visualization, Writing – original draft

E-mail: [email protected]

Affiliations Wada Eye Clinic, Chiba, Japan, Tsubota Laboratory, Inc., Tokyo, Japan

https://orcid.org/0000-0001-5604-5847

Reiko Arita

Roles: Project administration, Supervision, Writing – review & editing

Affiliations Tsubota Laboratory, Inc., Tokyo, Japan, Itoh Clinic, Saitama, Japan

Yasue Mitsukura

Roles: Supervision, Validation, Writing – review & editing

Affiliation: Faculty of Science and Technology, Keio University, Kanagawa, Japan

https://orcid.org/0000-0001-5575-8589

Reiko Ishida

Roles: Data curation

Affiliation: Ishida Eye Clinic, Shizuoka, Japan

Kazuo Tsubota

Roles: Supervision, Writing – review & editing

Affiliation: Tsubota Laboratory, Inc., Tokyo, Japan

https://orcid.org/0000-0002-8874-7111

References

1. Nichols KK, Nichols JJ, Mitchell GL. The lack of association between signs and symptoms in patients with dry eye disease. Cornea. 2004;23(8):762–770. pmid:15502475

2. Hua R, Yao K, Hu Y, Chen L. Discrepancy between subjectively reported symptoms and objectively measured clinical findings in dry eye: a population based analysis. BMJ Open. 2014;4:e005296. pmid:25168038

3. Lemp MA, Crews LA, Bron AJ, Foulks GN, Sullivan BD. Distribution of aqueous-deficient and evaporative dry eye in a clinic-based patient cohort: a retrospective study. Cornea. 2012 May;31(5):472–8. pmid:22378109.

4. Asiedu K, Dzasimatu SK, Kyei S. Clinical subtypes of dry eye in youthful clinical sample in Ghana. Cont Lens Anterior Eye. 2019 Apr;42(2):206–211. Epub 2018 Oct 15. pmid:30337142.

5. Tsubota K, Yokoi N, Shimazaki J, Dogru M, Kojima T, Yamada M, et al. New perspectives on dry eye definition and diagnosis: a consensus report by the Asia Dry Eye Society. Ocul Surf. 2017;15(1):65–76. pmid:27725302

6. Yokoi N, Uchino M, Uchino Y, Dogru M, Kawashima M, Komuro A, et al. Importance of tear film instability in dry eye disease in office workers using visual display terminals: the Osaka study. Am J Ophthalmol. 2015;159(4):748–754. pmid:25555800

7. Toda I, Shimazaki J, Tsubota K. Dry eye with only decreased tear break-up time is sometimes associated with allergic conjunctivitis. Ophthalmology. 1995;102(2):302–309. pmid:7862418

8. Tsubota K. Short tear film breakup time-type dry eye. Invest Ophthalmol Vis Sci. 2018;59(14):DES64–DES70. pmid:30481808

9. Tsubota K, Yokoi N, Watanabe H, Dogru M, Kojima T, Yamada M. A new perspective on dry eye classification: proposal by the Asia Dry Eye Society. Eye Contact Lens. 2020;46(1):S2–S13. pmid:31425351

10. Craig JP, Nichols KK, Akpek EK, Caffery B, Dua HS, Joo CK, et al. TFOS DEWS II definition and classification report. Ocul Surf. 2017;15(3):276–783. pmid:28736335

11. Sabatino F, Di Zazzo A, De Simone L, Bonini S. The intriguing role of neuropeptides at the ocular surface. Ocul Surf. 2017;15(1):2–14. pmid:27840126

12. Asiedu K, Markoulli M, Bonini S, Bron AJ, Dogru M, Kwai N, et al. Tear film and ocular surface neuropeptides: Characteristics, synthesis, signaling and implications for ocular surface and systemic diseases. Exp Eye Res. 2022;218:108973. pmid:35149082

13. Yang LWY, Mehta JS, Liu YC. Corneal neuromediator profiles following laser refractive surgery. Neural Regen Res. 2021;16(11):2177–2183. pmid:33818490

14. Hwang DD, Lee SJ, Kim JH, Lee SM. The role of neuropeptides in pathogenesis of dry dye. J Clin Med. 2021;10(18):4248. pmid:34575359

15. Asiedu K. Role of ocular surface neurobiology in neuronal-mediated inflammation in dry eye disease. Neuropeptides. 2022;95:102266. pmid:35728484

16. Bereiter DA, Rahman M, Thompson R, Stephenson P, Saito H. TRPV1 and TRPM8 channels and nocifensive behavior in a rat model for dry eye. Invest Ophthalmol Vis Sci. 2018;59(8):3739–3746. pmid:30046815

17. Kovács I, Luna C, Quirce S, Mizerska K, Callejo G, Riestra A, et al. Abnormal activity of corneal cold thermoreceptors underlies the unpleasant sensations in dry eye disease. Pain. 2016(2);157:399–417. pmid:26675826

18. Masuoka T, Yamashita Y, Nakano K, Takechi K, Niimura T, Tawa M, et al. Chronic tear deficiency sensitizes transient receptor potential vanilloid 1-mediated responses in corneal sensory nerves. Front Cell Neurosci. 2020;14:598678. pmid:33424555

19. Yu K, Guo Y, Ge QM, Su T, Shi WQ, Zhang LJ, et al. Altered spontaneous activity in the frontal gyrus in dry eye: a resting-state functional MRI study. Sci Rep. 2021;11(1):12943. pmid:34155240

20. Sun T, Shu HY, Wu JL, Su T, Liu YJ, Zhang LJ, et al. Investigation of changes in the activity and function of dry eye-associated brain regions using the amplitude of low-frequency fluctuations method. Biosci Rep. 2022;42(1):BSR20210941. pmid:34981112

21. Liang RB, Liu LQ, Shi WQ, Sun T, Ge QM, Li QY, et al. Abnormal fractional amplitude of low frequency fluctuation changes in patients with dry eye disease: a functional magnetic resonance imaging study. Front Hum Neurosci. 2022;16:900409. pmid:35693538

22. Kaido M, Inoue S, Kawashima M, Ishida R, Nakamura S, Tsubota K. Role of transient receptor potential melastatin 8 activity in menthol-induced cold sensitivity and its qualitative perception in dry eye. Ocul Surf. 2021;19:307–312. pmid:33127598

23. Sano K, Kawashima M, Ikeura K, Arita R, Tsubota K. Abdominal breathing increases tear secretion in healthy women. Ocul Surf. 2015;13(1):82–87. pmid:25557347

24. Jin K, Imada T, Hisamura R, Ito M, Toriumi H, Tanaka KF. Identification of lacrimal gland postganglionic innervation and its regulation of tear secretion. Am J Pathol 2020;190(5):1068–1079. pmid:32084368

25. Tokuyama T, Nakano Y, Awazu A, Uchimura-Makita Y, Fujiwra M, Watanabe Y, et al. Deterioration of the circadian variation of heart rate variability in Brugada syndrome may contribute to the pathogenesis of ventricular fibrillation. J Cardiol. 2014;64(2):133–138. pmid:24495503

26. Midorikawa-Inomata A, Inomata T, Nojiri S, Nakamura M, Iwagami M, Fujimoto K, et al. Reliability and validity of the Japanese version of the Ocular Surface Disease Index for dry eye disease. BMJ Open. 2019;9:e033940. pmid:31772113

27. Shimomitsu T. Final development of the brief job stress questionnaire mainly used for assessment of the individuals in the ministry of labor sponsored grant for the prevention of work-related illness. The 1999 report. 2000.

28. van Bijsterveld OP. Diagnostic tests in the Sicca syndrome. Arch Ophthalmol. 1969;82(1):10–14. pmid:4183019

29. Akselrod S, Gordon D, Ubel FA, Shannon DC, Berger AC, Cohen RJ. Power spectrum analysis of heart rate fluctuation: a quantitative probe of beat-to-beat cardiovascular control. Science. 1981;213(4504):220–222. pmid:6166045

30. Pomeranz B, Macaulay RJ, Caudill MA, Kutz I, Adam D, Gordon D, et al. Assessment of autonomic function in humans by heart rate spectral analysis. Am J Physiol. 1985;248(1 Pt 2):H151–H153. pmid:3970172

31. Hayano J, Sakakibara Y, Yamada A, Yamada M, Mukai S, Fujinami T, et al. Accuracy of assessment of cardiac vagal tone by heart rate variability in normal subjects. Am J Cardiol. 1991;67(2):199–204. pmid:1987723

32. Malik M, Camm AJ, Bigger JT, Breithardt G, Cerutti S, Cohen RJ. Heart rate variability: standards of measurement, physiological interpretation and clinical use. Task Force of the European Society of Cardiology and the North American Society of pacing and electrophysiology. Eur Heart J. 1996;17:354–381.

33. Hayano J, Taylor JA, Yamada A, Mukai S, Hori R, Asakawa T, et al. Continuous assessment of hemodynamic control by complex demodulation of cardiovascular variability. Am J Physiol 1993;264(4 Pt 2):H1229–H1238. pmid:8476100

34. Pagani M, Lombardi F, Guzzetti S, Rimoldi O, Furlan R, Pizzinelli P, et al. Power spectral analysis of heart rate and arterial pressure variabilities as a marker of sympatho-vagal interaction in man and conscious dog. Circ Res. 1986;59(2):178–193. pmid:2874900

35. Billman GE. Heart rate variability–a historical perspective. Front Physiol. 2011;2:86. pmid:22144961

36. McCraty R, Atkinson M, Tiller WA, Rein G, Watkins AD. The effects of emotions on short-term power spectrum analysis of heart rate variability. Am J Cardiol 1995;76(14):1089–1093. pmid:7484873

37. Kop WJ, Synowski SJ, Newell ME, Schmidt LA, Waldstein SR, Fox NA. Autonomic nervous system reactivity to positive and negative mood induction: the role of acute psychological responses and frontal electrocortical activity. Biol Psychol. 2011;86(3):230–238. pmid:21182891

38. Mauss IB, Robinson MD. Measures of emotion: a review. Cogn Emot. 2009;23(2):209–237. pmid:19809584

39. von Rosenberg W, Chanwimalueang T, Adjei T, Jaffer U, Goverdovsky V, Mandic DP. Resolving ambiguities in the LF/HF ratio: LF-HF scatter plots for the categorization of mental and physical stress from HRV. Front Physiol 2017;8:360. pmid:28659811

40. Ghiasi S, Greco A, Barbieri R, Scilingo EP, Valenza G. Assessing autonomic function from electrodermal activity and heart rate variability during cold-pressor test and emotional challenge. Sci Rep. 2020;10(1):5406. pmid:32214158

41. Amino M, Nakano M, Komatsu T, Yoshizawa R, Kunugita F, Kiyono K, et al. Prolonged autonomic fluctuation derived from parasympathetic hypertonia after carotid endarterectomy but not stenting. J Stroke Cerebrovasc Dis. 2019;28(1):10–20. pmid:30366866

42. Moss SE, Klein R, Klein BE. Prevalence of and risk factors for dry eye syndrome. Arch Ophthalmol. 2000;118(9):1264–1268. pmid:10980773

43. Zhang Y, Chen H, Wu X. Prevalence and risk factors associated with dry eye syndrome among senior high school students in a county of Shandong Province, China. Ophthal Epidemiol. 2012;19(4):226–230.