Abstract
Pupillary light reflex (PLR) and heart rate variability (HRV) parameters can be objective indicators of chronic rhinosinusitis (CRS) status from the viewpoint of autonomic nervous system activity. This study aimed to establish objective indicators for CRS using the 22-item Sino-Nasal Outcome Test (SNOT-22) and PLR/HRV parameters. Sixty-seven patients were prospectively and longitudinally followed up after surgical treatment. We investigated changes in SNOT-22 scores, representing CRS-specific quality of life (QOL). We prepared two models: linear regression model adjusting clinical factors as predictor variables (model 1) and linear mixed-effects model adjusting clinical factors and among-individual variability (model 2). We compared Akaike’s information criterion (AIC) values and regression coefficients. The model with lower AIC values was defined as the better-fit model. Model 2 showed lower AIC values in all parameters (better-fit model). Three parameters showed opposite results between the two models. The better-fit models showed significances in the five PLR parameters but not in any HRV parameters. Among these PLR parameters, constriction latency can be the most robust indicator because of the narrowest 95% confidence intervals. Adjusting the among-individual variability while investigating clinical potential of PLR/HRV parameters to reflect CRS-specific QOL can improve the model fit, thereby reaching robust conclusions from obtained data.
Similar content being viewed by others
Introduction
The nose possesses rich autonomic innervation of the nasal vasculature and nasal glands, contributing to the regulation of temperature and humidification of the air entering the lower airways1,2,3. Therefore, dysfunction of autonomic nervous system (ANS) activity may be associated with nasal symptoms and disorders. Chronic rhinosinusitis (CRS) is a representative upper-airway disorder characterized by a persisting inflammatory condition of the paranasal sinus. Some studies have investigated the association between CRS and ANS activity using tissue samples4,5 and subjective tools6.
One clinical study investigated the association between ANS dysfunction and CRS using the 31-item Composite Autonomic Symptom Score (COMPASS-31), a validated questionnaire for evaluating ANS dysfunction6. The results showed significant correlations between the 22-item Sino-Nasal Outcome Test (SNOT-22) and COMPASS-31 in participants with CRS with nasal polyps. Thus, these questionnaires can be used as subjective tools for assessing CRS-related ANS dysfunction, although they may be subject to bias. Therefore, the establishment of objective indicators reflecting CRS-related ANS activity can be helpful to obtain novel insights into the association between CRS and ANS activity. As the cardiac status and pupillary status are controlled by ANS activity7,8,9,10, recent studies investigated the association between ANS dysfunction and upper-airway disorders, such as obstructive sleep apnoea syndrome11,12,13 and allergic rhinitis14,15, using autonomic function tests for cardiac and/or pupillary status. However, to date, no clinical study has explored the association between ANS activity and CRS using objective indicators.
The pupillary light reflex (PLR) and/or heart rate variability (HRV) can be effective measures in investigating the association between ANS activity and physiological burden, such as that resulting from exercise16, foot bath17, and pain18. PLR is controlled by two kinds of muscles via the locus coeruleus containing both excitatory sympathetic premotor neurons and inhibitory parasympathetic premotor neurons, which is located in the brainstem: the iris sphincter muscle innervated by the parasympathetic nervous system (PNS) and the iris dilator muscle innervated by the sympathetic nervous system (SNS)7,10. HRV is the fluctuation in the length of the heartbeat R-R intervals, reflecting the involvement of ANS activity in cardiac regulation8,9. Thus, PLR mainly reflects the central autonomic regulation, while HRV reflects the peripheral autonomic regulation. Tests assessing these parameters might be effective in exploring the association between CRS and the central/peripheral ANS activity.
Deciding the study design and methods can be one of the critical issues in exploring the association between PLR/HRV parameters and CRS. SNOT-22 questionnaire scores represent CRS-specific quality of life (QOL)19 and have been used as the primary outcome in prospective20,21 and longitudinal22 clinical studies on CRS. As treatments including surgical intervention can result in long-term improvements in CRS-specific QOL22, there might be close associations between CRS-specific QOL and ANS activity. Therefore, a prospective longitudinal study investigating the association between SNOT-22 scores and PLR/HRV parameters in patients with CRS undergoing treatments, including surgical intervention, can contribute to exploring the relation between ANS activity and CRS. In a longitudinal study, data can be collected within a predefined period. Therefore, the data can show information about the degree of variations across time for each patient (within-individual variability) and among patients (between-individual variability). The use of linear mixed-effects models (LMMs) can be an appropriate statistical method for repeated data because LMMs possess the flexibility for assumptions and the advantage of adjusting individual variabilities23,24. Especially, as adjusting the individual variability can reduce false positives and false negatives in the obtained data24, the use of LMMs in a longitudinal study or repeated data can lead to robust conclusions without underestimating or overestimating the results.
The current study aimed to establish objective indicators of CRS status with the use of SNOT-22 and PLR/HRV parameters. For this purpose, we designed a 12-month prospective longitudinal study targeting patients with CRS who underwent treatments including surgical intervention and investigated the effects of adjusting the individual variability in results using two different statistical models.
Materials and methods
Patients
The Research Ethics Committee of Hyogo Prefectural Amagasaki General Medical Center approved this prospective longitudinal study.
Patients who underwent surgical intervention for medically refractory CRS at Hyogo Prefectural Amagasaki General Medical Center between August 2018 and April 2021 were enrolled. We performed this study in accordance with the tenets of the Declaration of Helsinki.
We collected sociodemographic and medical characteristics data through questionnaires, including sex, age, body mass index (BMI), smoking intensity according to the Brinkman index (number of cigarettes smoked per day multiplied by the number of smoking years), asthma, and diabetes. All patients signed an informed consent form after receiving information about the purpose and procedures of the study.
The inclusion criteria were age > 18 years and CRS that had been resistant to medical treatments for > 3 months. All patients underwent a biopsy when nasal polyps were observed endoscopically. The Japanese Epidemiological Survey of Refractory Eosinophilic Chronic Rhinosinusitis scoring system was preoperatively applied to the clinical data and biopsy samples to determine whether the patients had eCRS or neCRS25. The exclusion criteria were histories of glaucoma, arrhythmia, and cardiovascular disorders requiring a pacemaker. When enrolled patients did not visit postoperatively, these patients were excluded from the analysis.
BMI was classified into three categories26: (i) underweight, BMI < 18.5 kg/m2; (ii) healthy weight, 18.5 ≤ BMI < 25.0 kg/m2; and (iii) overweight and obese, 25 ≤ BMI kg/m2. Asthma status was classified as follows: (i) no asthma or (ii) inactive or active asthma. An inactive status was defined as having a history of asthma or asthma in remission, whereas an active status was defined as having current asthma requiring follow-up and/or treatments with specialists under the Japanese guidelines for adult asthma27. Diabetes status was classified as follows: (i) without diabetes or (ii) with diabetes requiring education for lifestyle changes, medication including oral hypoglycaemic agents, and insulin. CRS was divided into two phenotypes: (i) eCRS and (ii) neCRS.
Subjective and objective tools
All patients underwent subjective and objective measurements. As a subjective measurement tool, we used SNOT-22, which is a validated questionnaire that includes the 22 symptoms most related to CRS-specific QOL19. All symptoms were scored from 0 to 5 (0 = ‘no problem’, 1 = ‘very mild problem’, 2 = ‘mild or slight problem’, 3 = ‘moderate’, 4 = ‘severe problem’, and 5 = ‘problem as bad as it can be’). Higher scores represent worse CRS-specific QOL, and lower scores indicate improved QOL.
A hand-held infrared pupillometer (PLR-3000; NeurOptics, Irvine, CA, USA) and an HRV device (CheckMyHeart; Daily Care BioMedical, Taiwan) were used as objective measurement tools for the PLR and HRV tests, respectively. The details of the systems and tests have been previously described17,18.
In the PLR test, we assessed the following six parameters (Fig. 3A): (i) INIT (mm); (ii) ACV (mm/s); (iii) MCV (mm/s); (iv) constriction ratio, defined as the difference between the INIT and the minimum pupil size divided by the INIT (DELTA, %); (v) LAT (s); and (vi) ADV (mm/s). As PLR parameters can be affected by different light stimulus intensities (Fig. 3B)17,18, we used four different stimulus intensities using a flashlight with a fixed 800-ms pulse duration: (i) 10 μW, (ii) 50 μW, (iii) 121 μW, and (iv) 180 μW. The recording duration was 5 s. The following PLR parameters might be used as indicators of ANS activity: INIT as an indicator of sympathovagal balance28 and MCV10, DELTA10, and LAT29 as indicators of PNS activity.
In the HRV test, data recorded by the HRV device for 5 min (sampling rate: 250/s) were collected and analyzed using the software provided by the manufacturer (CheckMyHeart, Daily Care BioMedical). The software rejected the irregular R–R intervals automatically. The R–R intervals (resolution 4 ms) were obtained in the software. Cubical interpolation and resampling at 1 Hz were performed for the detrended time series. HRV analysis involves two different methods for assessing ANS activity: (i) time-domain analysis and (ii) frequency domain analysis. After the detrending procedures, the power spectral density of R-R intervals was calculated by an autoregressive method and a fast Fourier transform method. The order of autoregressive method was set as 16. We used three detrend parameters in the time-domain analysis: (i) standard deviation of the normal normal-to-normal (NN) intervals (SDNN, ms), (ii) root mean square of the successive NN interval differences (rMSSD, ms), and (iii) pNN50. We obtained five detrend parameters quantified using the autoregressive method for the frequency domain analysis: (i) low-frequency power (LF power, 0.04–0.15 Hz); (ii) high-frequency power (HF power, 0.15–0.4 Hz); (iii) normalized values of LF, (LF norm; LF/total power − very low frequency) × 100; (iv) normalized values of HF, (HF norm; HF/total power − very low frequency) × 100; and (v) LF/HF values. Previously, raw values such as LF and HF powers have been reported to have disadvantages due to large inter-variability and intra-variability, resulting in a long-tailed right-skewed exponential distribution8,17. Therefore, we excluded LF and HF powers. Finally, we adopted six HRV parameters for analysis: (i) SDNN, (ii) rMSSD, (iii) pNN50, (iv) LF norm, (v) HF norm, and (vi) LF/HF. The following parameters can be used as indicators of ANS activity: rMSSD, pNN50, and HF norm as indicators of PNS activity; SDNN and LF norm as an indicator of SNS and PNS activity modulation; and LF/HF as an indicator of sympathovagal balance8,9. However, these concepts have been recently criticized because of the non-linear interactions between SNS and PNS activity30. Thus, the interpretation of the HF norm, LF norm, and LF/HF as autonomic indicators remains debatable.
Treatments for CRS
The surgical strategies were determined on the basis of computed tomographic and endoscopic findings. All patients underwent endoscopic sinus surgery (ESS). The procedures consisted of either unilateral or bilateral maxillary antrostomy, ethmoidectomy, frontal sinusotomy, or sphenoidotomy, with septoplasty and inferior turbinate surgery as adjunctive procedures when needed. The patients underwent ESS under general anesthesia and received nasal packing (Merocel; Medtronic Xomed, Minneapolis, MN, USA) on the surgical side just before awakening to prevent postoperative bleeding. The nasal packing was removed 1 day after ESS. As a postoperative care procedure, all patients were recommended to perform saline irrigation during the follow-up period. Some patients, especially those with eCRS, can experience recurrence even after ESS31. When patients required additional treatments because of postoperative recurrence, a short course of systemic corticosteroid therapy (oral betamethasone: 1 mg/day for 1 week) or antibiotic therapy was administered. The number of the systemic corticosteroid medication courses was limited to two per month.
Study design
The PLR and HRV tests were performed in a seated position in a quiet consulting room kept at a comfortable temperature (21–25 °C). The patients were required to keep their position in a relaxed state during the tests. All patients received the short-term HRV test, followed by the PLR test.
We performed PLR and HRV tests at eight time points: (i) 1 week before ESS (defined as a reference), (ii) 1 week after ESS, (iii) 2 weeks after ESS, (iv) 1 month after ESS, (v) 2 months after ESS, (vi) 3 months after ESS, (vii) 6 months after ESS, and (viii) 12 months after ESS. The assessment time was between 8 am and 12 pm at all time points to minimize the effects of the circadian rhythm.
Statistical analysis
We expressed quantitative variables as means and standard deviations/95% CIs or as medians and interquartile ranges. We compared the SNOT-22 total symptom scores as non-parametric data across the time points using the Friedman test with Dunn’s multiple comparison test.
We investigated the association between SNOT-22 scores and PLR/HRV parameters using linear regression models. To reveal the effects of adjusting individual variability, we prepared two models: (i) linear regression model adjusting clinical factors as predictor variables (model 1) and (ii) LMM adjusting clinical factors and individual variability (model 2). The linear regression model consisted of fixed effects as predictor variables, whereas the LMM consisted of fixed and random effects. Random effects were subdivided into random intercepts and slopes, although the random slope models can be more prone to convergence problems than the random intercept models24. Therefore, we selected the random intercept model for the LMM. We introduced the total SNOT-22 scores, sex, age, smoking intensity, BMI categories, CRS phenotypes, measured eye side (only for PLR analyses), light stimulus intensity (only for PLR analyses), and the status of asthma and diabetes as fixed effects, whereas each patient ID was used as random effects for adjusting individual variability. PLR/HRV parameters were applied as response variables. The SNOT-22 total symptom scores, values of PLR and HRV parameters, light stimulus intensity, age, and smoking intensity were used as continuous variables, whereas sex, measured eye side, CRS phenotypes, and status of asthma and diabetes were used as binary variables. BMI was used as a categorical variable. To evaluate the effects of adjusting among-individual variability, we compared AIC values and regression coefficients of PLR/HRV parameters for every 1 unit increase in the SNOT-22 total symptom score. The model with the lower AIC values was defined as the better-fit model32. A two-sided p-value of < 0.05 was considered statistically significant. All statistical analyses were performed with R (version 3.6.1; R Foundation for Statistical Computing, Vienna, Austria; lme4 and lmerTest packages).
Ethics approval and consent to participate
The Research Ethics Committee of Hyogo Prefectural Amagasaki General Medical Center approved this study (approval number: 2-184). We obtained written informed consents from all subjects.
Results
Characteristics of patients
A total of 67 Japanese patients were enrolled in the current study. Table 1 shows the sociodemographic and medical characteristics of the patients.
Changes in total SNOT-22 scores over time
We compared the SNOT-22 total symptom scores across the eight time points (Fig. 1). The scores significantly changed during the follow-up period (p < 0.001). Compared with the preoperative values, post-hoc tests revealed a significant decrease in the scores from postoperative 1 to 12 months, indicating an improvement in CRS-specific QOL after treatments.
Association between SNOT-22 scores and PLR/HRV parameters
To establish PLR/HRV parameters as objective indicators reflecting CRS status, we investigated the association between SNOT-22 total scores and PLR/HRV parameters using two linear regression models. In the comparison of AIC values between the two models, model 2 showed lower Akaike’s information criterion (AIC) values for all parameters (Tables 2 and 3). Therefore, we defined model 2 as the better-fit model.
To evaluate the effects of adjusting individual variability with LMMs, we compared the regression coefficients of PLR/HRV parameters for every 1 unit change in the SNOT-22 total symptom scores (Tables 2 and 3). In PLR/HRV analyses, both models showed significances in five PLR parameters (see Models 2 in Table 2) but not in any HRV parameter (see Models 2 in Table 3). Initial pupil size (INIT) in PLR and percentage of adjacent normal-to-normal intervals that differ from each other by > 50 ms (pNN50) showed opposite results between the two models: model 1 showed significances but not model 2. In the constriction ratio (DELTA), model 1 showed a significance with negative values, while model 2 showed one with positive values. As shown in Fig. 2, decreased SNOT-22 total scores were associated with significant decreases in average constriction velocity (ACV), maximum constriction velocity (MCV), DELTA, and average dilation velocity (ADV), and prolonged constriction latency (LAT). Out of these PLR parameters, LAT showed the narrowest 95% confidence intervals (CIs).
Discussion
We assessed the association between PLR/HRV parameters and SNOT-22 total symptom scores using two statistical models in this prospective longitudinal study targeting patients with CRS. The results showed that model 2 had lower AIC values in all parameters. Three parameters showed opposite results between the two models: DELTA, INIT, and pNN50. Moreover, model 2 showed significant parameter changes with altered SNOT-22 scores in five PLR parameters but not HRV parameters. Among the PLR parameters, LAT showed the narrowest 95% CIs.
This study showed that the better-fit model of model 2 for all parameters and the opposite results between the two models in three parameters. These findings indicated that adjusting the among-individual variability with LMMs can improve the model fit, contributing to the reduction of false-positive ratios and the prevention of misleading conclusions. The findings were consistent with our previous study33.
Previously, we investigated whether PLR/HRV parameters can be altered by conditions that create a physiological burden, such as pain18 and foot bath17. Both PLR and HRV parameters showed significant alterations with altered subjective pain intensity scores in our previous study targeting clinical pain relief models18. However, our previous study investigating the effects of foot bath demonstrated significant alterations after foot bath in five out of seven PLR parameters but in only a few HRV parameters17. These findings indicated the importance of evaluating ANS activity targeting several organs and tissues. Our current study showed significant regression coefficients with altered SNOT-22 scores in five PLR parameters but not in any HRV parameters. The different responses might be due to the different sensitivities to altered CRS-specific QOL between PLR and HRV parameters. PLR parameters might be more effective than HRV parameters in assessing the association between CRS-specific QOL and ANS activity. Among the five PLR parameters, the narrowest 95% CIs was observed in LAT (see Fig. 2). The findings indicated that LAT is the most robust indicator reflecting the CRS-specific QOL.
Researchers may be interested in interpreting parameter alterations from the viewpoint of ANS activity. Although some HRV parameters remain debatable as autonomic roles (HF norm, LF norm, and LF/HF)30,34, this study showed that the altered ANS activity in the heart due to treatments can be less intense than that in the pupil. In PLR wave, the three phases are under the control of PNS and SNS activities (Fig. 3A): (i) phase 1 is mainly controlled by PNS activity; (ii) phase 2 is controlled by PNS and SNS activities; and (iii) phase 3 is predominantly controlled by SNS activity10. Therefore, parameters obtained from phase 1, such as ACV, MCV, DELTA, and LAT, can be used as indicators of PNS activity. For these parameters, a literature review reported MCV and relative constriction amplitude (i.e., DELTA) as indicators sensitive to PNS activity but questioned whether the same is true for LAT10. Previously, we investigated how PLR parameters can be altered by different light stimulus intensities in healthy individuals17. We observed shortened LAT and increased ACV, MCV, and DELTA with increased light stimulus intensity. As pupillary constriction can be a consequence of increased PNS outputs10, changes in LAT, ACV, MCV, and DELTA can be considered to indicate increased PNS activity17. As ADV is a parameter obtained from phase 2/3, it can be affected by SNS and PNS activities. However, the role of ADV as an indicator of ANS activity remains unclear. Taking these findings together, five PLR parameters can effectively reveal the association between CRS and ANS activity: ACV, MCV, DELTA, LAT, and ADV. The interpretation of altered PLR parameters from viewpoints of ANS activity requires attention as some PLR parameters remain unclear as autonomic indicators with reliability as well as HRV parameters. Prolonged LAT and decreased ACV, MCV, and DELTA due to decreased SNOT-22 scores, as observed in the current study, might be considered to indicate decreased PNS activity due to improvement of CRS-specific QOL.
This study had several limitations. First, the PLR alterations observed in the current study may have been affected by unconsidered factors such as psychological status35. As the five symptom domains, including psychological dysfunction, were defined using SNOT-2236, assessing the association between the PLR parameters and the domain might be effective in revealing the effects of psychological status on PLR parameters. In addition, the measurement was performed at a single time point before treatment. Preoperative measurements at several time points should be performed, and test–retest variability should be assessed. Second, this was designed as a one-group pre-treatment versus post-treatment longitudinal study. Therefore, it remains unclear whether the PLR parameter alterations observed in our current study can represent hyperactivity of the PNS in patients with CRS because we targeted a single CRS group. Further studies including several groups might overcome these issues. Third, we used SNOT-22 to evaluate the effects of CRS treatment. Although SNOT-22 can be an effective subjective tool for analyzing the CRS status, the scores represent CRS-specific QOL but not ANS dysfunction. Therefore, the use of ANS-focused questionnaires, such as the COMPASS-31, might provide different insights into the association between CRS-related PLR/HRV parameter alterations and ANS activity. Fourth, the current study data were obtained from a small sample size. Since eosinophilic CRS (eCRS) has close associations with type 2 immune responses31, and type 2 immune responses can be modulated by acetylcholine derived from PNS activity37, researchers might be interested in the differences in PLR/HRV parameter values between non-eCRS (neCRS) and eCRS. Future studies with larger samples for both neCRS and eCRS would enable assessing whether the CRS phenotype can affect ANS activity. Fifth, it was impossible to discuss the associations between SNS activity and CRS because PLR is predominantly controlled by PNS activity10. As INIT is an indicator of sympathovagal balance28, the lack of significant alterations in INIT observed in this study might reveal the decreased SNS and PNS activities with the same extent (non-linear interactions between SNS and PNS activity).
In summary, this study revealed the clinical potential of five PLR parameters as objective indicators of CRS-specific QOL, with LAT being the most robust. To obtain robust conclusions from PLR/HRV data obtained from a longitudinal study, not only clinical factors but also individual variability should be adjusted. As the autonomic roles in some parameters remain unclear, further studies are necessary to reveal the association between CRS and ANS activity.
Data availability
The datasets are available from the corresponding author upon request.
Abbreviations
- ACV:
-
Average constriction velocity
- ADV:
-
Average dilation velocity
- AIC:
-
Akaike’s information criterion
- ANS:
-
Autonomic nervous system
- BMI:
-
Body mass index
- COMPASS-31:
-
31-Item Composite Autonomic Symptom Score
- CRS:
-
Chronic rhinosinusitis
- DELTA:
-
Constriction ratio
- eCRS:
-
Eosinophilic chronic rhinosinusitis
- ESS:
-
Endoscopic sinus surgery
- HF:
-
High frequency
- HRV:
-
Heart rate variability
- INIT:
-
Initial pupil size
- LAT:
-
Constriction latency
- LF:
-
Low frequency
- LMM:
-
Linear mixed-effects model
- MCV:
-
Maximum constriction velocity
- neCRS:
-
Non-eosinophilic chronic rhinosinusitis
- norm:
-
Normalised
- PLR:
-
Pupillary light reflex
- pNN50:
-
Percentage of adjacent normal-to-normal intervals that differ from each other by > 50 ms
- PNS:
-
Parasympathetic nervous system
- QOL:
-
Quality of life
- rMSSD:
-
Root mean square of successive normal-to-normal interval differences
- SDNN:
-
Standard deviation of normal normal-to-normal intervals
- SNOT-22:
-
22-Item Sino-Nasal Outcome Test
- SNS:
-
Sympathetic nervous system
References
Naclerio, R. M., Bachert, C. & Baraniuk, J. N. Pathophysiology of nasal congestion. Int. J. Gen. Med. 3, 47–57 (2010).
de Assis EV, de Andrade IsidÃrio U, Feitosa AD et al. Autonomic nervous system and allergic diseases: Integrative literature review. Int. Arch. Med. 8(112), 1–6 (2015).
Sahin-Yilmaz, A. & Naclerio, R. M. Anatomy and physiology of the upper airway. Proc. Am. Thorac. Soc. 8, 31–39 (2011).
Cerejeira, R., Fernandes, S. & Pinto, M. C. Increased expression of alpha7nAChR in chronic rhinosinusitis: The intranasal cholinergic anti-inflammatory hypothesis. Auris Nasus Larynx 43, 176–181 (2016).
Do, H. B. et al. Decreased ciliary beat responsiveness to acetylcholine in the nasal polyp epithelium. Clin. Otolaryngol. 44, 356–365 (2019).
Chen, W. C. et al. The symptom burden of autonomic dysfunction is positively associated with chronic rhinosinusitis status. Rhinology 56, 227–233 (2018).
Szabadi, E. Functional organization of the sympathetic pathways controlling the pupil: Light-inhibited and light-stimulated pathways. Front. Neurol. 9, 1069 (2018).
Burr, R. L. Interpretation of normalized spectral heart rate variability indices in sleep research: A critical review. Sleep 30, 913–919 (2007).
Shaffer, F. & Ginsberg, J. P. An overview of heart rate variability metrics and norms. Front. Public Health 5, 258 (2017).
Wang, Y. et al. Parasympathetic nervous system dysfunction, as identified by pupil light reflex, and its possible connection to hearing impairment. PLoS One 11, e0153566 (2016).
Sequeira, V. C. C., Bandeira, P. M. & Azevedo, J. C. M. Heart rate variability in adults with obstructive sleep apnea: A systematic review. Sleep Sci. 12, 214–221 (2019).
Lombardi, C., Pengo, M. F. & Parati, G. Obstructive sleep apnea syndrome and autonomic dysfunction. Auton. Neurosci. 221, 102563 (2019).
Correia, F. J., Martins, L. E. B., Barreto, D. M. & Pithon, K. R. Repercussion of medium and long treatment period with continuous positive airways pressure therapy in heart rate variability of obstructive sleep apnea. Sleep Sci. 12, 110–115 (2019).
Ozsutcu, M. et al. Pupillometric assessment of autonomic nervous system in children with allergic rhinitis. Med. Princ. Pract. 22, 444–448 (2013).
Lan, M. Y., Lee, G. S., Shiao, A. S., Ko, J. H. & Shu, C. H. Heart rate variability analysis in patients with allergic rhinitis. Sci. World J. 2013, 947385 (2013).
Kaltsatou, A., Kouidi, E., Fotiou, D. & Deligiannis, P. The use of pupillometry in the assessment of cardiac autonomic function in elite different type trained athletes. Eur. J. Appl. Physiol. 111, 2079–2087 (2011).
Ishikawa, M. Effective combination of pupil light reflex and heart rate variability to assess foot bath effects on autonomic function in healthy adults. Biomed. Phys. Eng. Express 6, 015034 (2020).
Ishikawa, M. Clinical potential of pupillary light reflex and heart rate variability parameters as biomarkers for assessing pain relief effects on autonomic function: A prospective longitudinal study. Biomed. Phys. Eng. Express 6, 055003 (2020).
Hopkins, C., Gillett, S., Slack, R., Lund, V. J. & Browne, J. P. Psychometric validity of the 22-item Sinonasal Outcome Test. Clin. Otolaryngol. 34, 447–454 (2009).
DeConde, A. S. et al. Comparative effectiveness of medical and surgical therapy on olfaction in chronic rhinosinusitis: A prospective, multi-institutional study. Int. Forum Allergy Rhinol. 4, 725–733 (2014).
Hopkins, C. et al. Long-term outcomes from the English national comparative audit of surgery for nasal polyposis and chronic rhinosinusitis. Laryngoscope 119, 2459–2465 (2009).
DeConde, A. S. et al. Longitudinal improvement and stability of the SNOT-22 survey in the evaluation of surgical management for chronic rhinosinusitis. Int. Forum Allergy Rhinol. 5, 233–239 (2015).
Finucane, M. M., Samet, J. H. & Horton, N. J. Translational methods in biostatistics: Linear mixed effect regression models of alcohol consumption and HIV disease progression over time. Epidemiol. Perspect. Innov. 4, 8 (2007).
Harrison, X. A. et al. A brief introduction to mixed effects modelling and multi-model inference in ecology. PeerJ 6, e4794 (2018).
Tokunaga, T. et al. Novel scoring system and algorithm for classifying chronic rhinosinusitis: The JESREC Study. Allergy 70, 995–1003 (2015).
Weir, C. B. & Jan, A. BMI Classification Percentile and Cut Off Points. (BTI-StatPearls, 2021).
Ichinose, M. et al. Japanese guidelines for adult asthma 2017. Allegol. Int. 66, 163–189 (2017).
Fotiou, F., Fountoulakis, K. N., Goulas, A., Alexopoulos, L. & Palikaras, A. Automated standardized pupillometry with optical method for purposes of clinical practice and research. Clin. Physiol. 20, 336–347 (2000).
Pfeifer, M. A. et al. Quantitative evaluation of sympathetic and parasympathetic control of iris function. Diabetes Care 5, 518–528 (1982).
Billman, G. E. The LF/HF ratio does not accurately measure cardiac sympatho-vagal balance. Front. Physiol. 4, 26 (2013).
Fokkens, W. J. et al. European position paper on rhinosinusitis and nasal polyps 2020. Rhinology 58, 1–464 (2020).
Garson, G. D. Hierarchical Linear Modeling: Guide and Applications (SAGE Publications Inc., 2013).
Ishikawa, M. Clinical factors affecting pupillary light reflex parameters: A single-centre, cross-sectional study. Ophthalmic Physiol. Opt. 41, 952–960 (2021).
Houle, M. S. & Billman, G. E. Low-frequency component of the heart rate variability spectrum: A poor marker of sympathetic activity. Am. J. Physiol. 276, H215–H223 (1999).
Bradley, M. M., Miccoli, L., Escrig, M. A. & Lang, P. J. The pupil as a measure of emotional arousal and autonomic activation. Psychophysiology 45, 602–607 (2008).
DeConde, A. S., Bodner, T. E., Mace, J. C. & Smith, T. L. Response shift in quality of life after endoscopic sinus surgery for chronic rhinosinusitis. JAMA Otolaryngol. Head Neck Surg. 140, 712–719 (2014).
Bosmans, G. et al. Cholinergic modulation of type 2 immune responses. Front. Immunol. 8, 1873 (2017).
Acknowledgements
The authors thank Drs. Youhei Kumabe, Takeshi Morita, and Syuhei Ootsuki for the clinical care of the patients and Editage (www.editage.com) for English language editing. We obtained permission to use SNOT-22 tool for this study from Washington University, USA.
Funding
This research was supported by the Japanese Society of Allergology, Clinical Research Support Program.
Author information
Authors and Affiliations
Contributions
M.I.: Study conception and design, acquisition of data, analysis and interpretation of data, drafting of manuscript, and critical revision. H.H.: Acquisition of data, drafting of manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Hatsukawa, H., Ishikawa, M. Clinical potential of pupillary light reflex parameters as objective indicators reflecting chronic rhinosinusitis-specific quality of life: a 12-month prospective longitudinal study. Sci Rep 11, 21586 (2021). https://doi.org/10.1038/s41598-021-01153-1
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-021-01153-1
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.