Effects of phytoncide on immune cells and psychological stress of gynecological cancer survivors: randomized controlled trials

Article information

J Exerc Rehabil Vol. 19, No. 3, 170-180, June, 2023
Publication date (electronic) : 2023 June 28
doi : https://doi.org/10.12965/jer.2346150.075
1Department of Physical Education, Chungnam National University, Daejeon, Korea
2Research Institute of Sports and Industry Science, Hanseo University, Seosan, Korea
3Exercise Immunology Center, Seoul Songdo Hospital, Seoul, Korea
*Corresponding author: Yong-Seok Jee, https://orcid.org/0000-0001-6797-0843, Research Institute of Sports and Industry Science, Hanseo University, Hanseo 1-ro, Haemi-myeon, Seosan 31962, Korea, Email: jeeys@hanseo.ac.kr
*Co-correspondence author: Sang-Kyun Park, https://orcid.org/0000-0002-8140-5194, Department of Physical Education, Chungnam National University, Daehak-ro, Yuseong-gu, Daejeon 34134, Korea, Email: ttl0033@gmail.com
Received 2023 April 03; Accepted 2023 April 23.

Abstract

This study aimed to identify neuroendocrine substances and natural killer (NK) cells, including their subsets and receptors, to determine whether phytoncides scented in an urban hospital could help improve stress in cancer survivors. Fifty-five gynecological cancer survivors were assigned to either the control group (CG, n=28) or phytoncide group (PTG, n=27). The PTG meditated by lying down in a space scented with phytoncide for 1 hr a day, 5 days a week for 8 weeks. Stress levels in both groups were high prior to the experiment and decreased by 9.31%±45.98% (P=0.003) only in the PTG after the experiment. Although the parasympathetic nerve activity of the PTG increased, the epinephrine and cortisol levels were significantly (P<0.001) decreased by 5.29%±25.5% and by 24.94%±11.62%. Moreover, the PTG showed a significant increase in the levels of NK cell subsets after 8 weeks, whereas there was no improvement in the CG. In conclusion, phytoncide fragrance reduces stress, increases the number of NK cells and their family even in a nonforest environment, and improves innate immunocytes in gynecological cancer survivors; parasympathetic nerve activity and cortisol hormones play critical roles in this process. That is, a phytoncide essential oil helps to stimulate changes in immunocytes’ mobility by affecting the human nervous and endocrine systems, thereby providing relief for psychological stress among cancer survivors who previously had cancer cells.

INTRODUCTION

Although cancer survivorship is clearly increasing, there appears to be no improvement in the rate among older women over 70 years of age (Lundberg et al., 2022). The incidence of gynecological cancers has gradually increased (Benedet, 2000; Zhang et al., 2021). Cancer survivors experience considerable stress due to health concerns such as cancer cell metastasis and the sense of isolation they have experienced in the past (Tang et al., 2016). Long-term stress may lead to the rapid growth of cancer cells or metastasis of new cancer cells in patients who underwent treatment previously (Neeman and Ben-Eliyahu, 2013). Therefore, it is necessary to develop an intervention program that can help reduce the stress of survivors and an immune-enhancing therapy program capable of blocking the metastasis of new cancer cells.

Several studies have shown that the fragrance of phytoncides reduces stress and promotes immunocyte growth, including natural killer (NK) cells (Akram et al., 2020; Li et al., 2007; Vivier et al., 2008). Phytoncides are volatile substances that relieve stress and exert antibacterial effects when inhaled from the forests (Jo et al., 2021; Li et al., 2006; Li et al., 2007; Li et al., 2008), and they are most released from conifers, particularly pines, nut pins, and cypress trees. Essential oils from the Hinoki cypress and pine trees are primarily used in commercial products containing various terpenes (Jo et al., 2021; Lee et al., 2015), of which monoterpenes have adrenal cortex stimulation, antisepsis, sterilization, and antiviral effects. Tannin, the second most abundant component of phytoncides, has detoxification, hemostatic, and anti-inflammatory effects (Akram et al., 2020).

People have long believed that living in the woods improves the health of patients with cancer in Korea; however, the reason for this remains unclear. Although some researchers have argued that phytoncides activate the NK cells of patients with cancer to enhance immunocytes, others have suggested that phytoncides through the forests help reduce stress and the negative effects of chemotherapy (Tønnesen et al., 1987). Bodner et al. (1998) reported that when the level of cortisol increases, that of NK cell decreases, whereas when the level of cortisol decreases, that of NK cell increases. Tønnesen et al. (1987) reported that adrenaline produced an immediate increase in NK cell activity and a selective increase in circulating NK cells. Salleh (2008) also reported that NK cells respond sensitively to changes in cortisol or epinephrine in the blood.

The diverse experimental settings make it challenging to draw definitive conclusions from the findings of multiple previous studies regarding the effects of phytoncide fragrance on stress hormone responses and adaptations. That is, as reported in previous studies, stress hormones affect NK cells, and it can be said that phytoncide fragrance affects NK cells because it affects stress hormones. However, research conducted using phytoncide fragrance made from essential oils in urban areas is extremely rare, as research has only been conducted in external forest environments. Specifically, the mechanism by which the inhalation of the scent of phytoncides oil affects the endocrine system, the nervous system, and NK cells and subsets remains unclear. Thus, this study investigated the effects of phytoncide oil fragrance on the stress index, neuroendocrine concentration, and NK cell family in gynecological cancer survivors.

MATERIALS AND METHODS

Study design and participants

This prospective single-blinded randomized controlled trial was conducted in accordance with the Declaration of Helsinki and approved by the ethics committee (HS22-11-01). The study was registered with the Korean Clinical Research Information Service (KCT0008283). All the participants listened to the experimental procedures and provided informed consent before the study, and all data were acquired from November 2, 2022, to January 3, 2023. The participants were between 61 and 79 years old and gynecological cancer (cervical, ovarian, or breast cancers) survivors who were diagnosed between January 2016 and December 2020 and underwent surgery, radiation therapy, or chemotherapy, with a survival rate of less than 5 years. They were subjected to continuous management at Seoul Songdo Hospital, Korea. Patients with serious inflammation, osteoporosis, or psychological disorders were excluded. Although this study aimed to recruit 60 participants considering the dropout rate, 63 participants initially participated. Only 62 participants were included in the random classification because one person encountered the coronavirus disease. Participants were randomized into groups using the Research Randomizer program (Urbaniak and Plous, 2023); each participant was asked to draw a numbered ticket from a jar, where the number ranged from 1 to 62. Thirty-one participants each were assigned to the control group (CG) or phytoncide group (PTG). Some participants in both groups dropped out during the experiment (Fig. 1).

Fig. 1

Participants’ allocation (consolidated standards for reporting of trials flow diagram).

Intervention

Phytoncides oil (Chamaecyparis obtusa) was selected for olfactory stimulation (Raha et al., 2019), and an oil burner (Good Earth Lab, Seoul, Korea) made of cypress wood was used. The main oil components were terpenes, primarily monoterpenes and sesquiterpene (Lee et al., 2010). Three to five drops of phytoncides oil were dropped onto the burner to allow the wood to absorb and emit a subtle scent. The temperature and humidity of the closed laboratory space (approximately 66 m2) were maintained at 25°C–26°C and 50%–60%, respectively, to maintain the scent diffusion. The PTG meditated by lying down in a space scented with phytoncide for 1 hr a day, 5 days a week, for 8 weeks. The CG was allowed to meditate while lying down for the same time and period in a different place without phytoncide scent. The participants were classified according to their community areas to prevent communication between the groups.

Stress scale

This study used a 10-item version of the perceived stress scale (PSS) to identify stress levels (Cohen et al., 1983). PSS consisted of six positive items and four negative items. Each item was rated based on the past month on a Likert scale, which composed ‘never’ (0) to ‘very often’ (4) (Table 1). In scoring the measures, the positive items were reverse-scored and then aggregated (Golden-Kreutz, 2004). The score ranges 0–13, 14–26, and 27–40 were considered low, moderate, and high perceived stress, respectively. The reliability coefficient, assessed using Cronbach alpha, was 0.913. We conducted a structural model analysis to confirm if the potential variables on the stress scale explain the observed variables.

Perceived stress scale for gynecological cancer survivor

Body composition and controlled variables measures

The participants were asked to limit their food and caffeine intake before the test, urinate and not exercise the day before the test. All participants removed the metals attached to their bodies, and the body composition of each participant was measured using an InBody 770 (Biospace Co. Ltd., Seoul, Korea) after their height was measured. Additionally, daily diet by CAN-Pro 5.0 (The Korean Nutrition Society, 2023) and physical activity level by the International Physical Activity Questionnaire (Booth, 2000) were recorded. Calorie intake and expenditure through physical activity were checked daily, recorded on a computer log, and averaged at the end of the month.

Autonomic nervous system measures

All participants were asked to abstain from working out, ingest food 3 hr, and take a break 30 min before the test to obtain a measurement of heart rate variability (HRV). The uBioMacpa (Biosensecreative Co. Ltd., Seoul, Korea) was used for HRV test. Frequency analysis of HRV provided data on the following HRV components: high frequency and low frequency (Lane et al., 2009). Since it was clear that high frequency is an activity indicator of the parasympathetic nervous system (PNS), the prevailing view was that low frequency mainly reflects sympathetic nervous system (SNS) activity (White and Raven, 2014).

Stress hormones and immunocytes measures

The participants had fasted for at least 10 hr and 20 mL of blood was taken from the antecubital vein. It was centrifuged at 3,000 rpm for 15 min and stored in a −80°C freezer until analysis. Whole blood samples were used for automated differential blood cell count and fluorescence-activated cell sorting (FACS) analyses. Peripheral blood samples were collected from heparinized collection tubes and ethylenediaminetetraacetic acid (EDTA) tubes. Cortisol levels were measured using an enzyme-linked immunospecific assay (ELISA) kit (DRG Instruments GmbH, Hamburg, Germany). Epinephrine levels were analyzed using a high-performance liquid chromatography with a plasma catecholamine ELISA kit (Labor Diagnostika Nord, Nordhorn, Germany). EDTA-anticoagulant blood samples were tested using Sysmex XN-550 (Sysmex Corp., Kobe, Japan). Particularly, the samples were stored at room temperature for 4 hr before the test. A melted whole-blood technique with a maximum of 8-color staining was used for flow cytometry and antibody staining. Notably, 50 μL of blood was with anti-human antibodies against anti-CD56, anti-CD3, anti-CD314 (NKG2D), anti-CD158b (KIR2DL3) (BioLegend, San Diego, CA, USA), and isotype controls (Martinez-Borra and Khakoo, 2008). Moreover, T-cell immunoreceptors with immunoglobulin (Ig) and immunoreceptor tyrosine-based inhibitory motif (TIGIT) and lymphocyte function-associated antigen-1 (LFA-1) were analyzed. After incubation for 30 min at room temperature in the dark, the erythrocyte was lysed by adding 500 μL of FACS lysing solution to each test tube for 15 min. Upon completion of staining, the cells were analyzed using FACS Canto II (BD Biosciences, San Jose, CA, USA) and the FlowJo program (Treestar, Ashland, OR, USA) and the results are presented as percentages.

Statistical analyses

Utilizing G·Power program, the sample size was obtained based on an a priori effect size (0.25) of f 2 (V), α error probability (0.05), power (0.95), two groups, and two measurements (Faul et al., 2007), and the total sample size was 54. The software used in this study for analyses were SPSS and AMOS (ver. 25; IBM Corp., Armonk, NY, USA) and GraphPad Prism 9.5.1 (Graph Pad Software, San Diego, CA, USA). Normal data distribution was estimated using D’Agostino and Pearson omnibus tests. The stress scale showed scores for 10 questions (A to J) for each participant in each group using a heatmap. Confirmatory factor analysis was conducted to verify the validity of the observed variables constituting the latent variables of the stress questionnaire. This study analyzed the prescores of the variables among the groups using the Mann–Whitney U-test. Two-way analysis of variance and Wilcoxon signed-rank tests were used to compare the results between groups and times. A delta percentage (%) analysis was also conducted, and effect size (η2) was computed according to Cohen d (Cohen, 1992). All data are presented as mean±standard deviation and analyzed based on P≤0.05.

RESULTS

Demographic and clinical characteristics

As shown in Table 2, there were no significant differences in age, body composition components, demographic properties, surgical stages, and therapeutic methods between the groups before the experiment.

Demographic and clinical characteristics of the participants

Specifically, the calorie intake of CG at week 4 was 1,608.36± 176.31 kcal/day, whereas the calorie intake of PTG was 1,588.85± 167.51 kcal/day, showing no statistically significant difference (Z= −0.338; P=0.736, η2=0.003). Even at week 8, the calorie intake of CG was 1,752.29±249.99 kcal/day, whereas the calorie intake of PTG was 1,725.41±141.00 kcal/day, showing no significant difference (Z=−1.207; P=0.228, η2=0.004) between the groups. The calorie consumption of CG was 534.21±62.79 kcal/day, whereas that of PTG was 575.37±10.04 kcal/day at week 4, showing no statistically significant difference (Z=−1.672; P=0.095, η2= 0.103). At week 8, the average calorie consumption of CG was 583.82±33.43 kcal/day, whereas that of PTG was 583.30±32.42 kcal/day, showing no significant difference (Z=−0.110; P=0.913, η2=0.007) between the groups.

Stress scale and implication of phytoncide fragrance on stress level

In the fit of the assessing model, the χ2 (df=35) value was 170.58 (P<0.001). Since the χ2 value was greatly affected by the sample size and complexity of the model, other goodness-of-fit indices were identified. The results showed that the comparative fit index=0.905, Tucker–Lewis index=0.914, and root mean square error of approximation=0.036 were good, and the test model was found to be suitable (Fig. 2). All factor loadings were found to have a significant positive effect in the confirmatory factor analysis except for Q3 and Q5 (Table 3).

Fig. 2

Confirmatory factor analysis model.

Confirmatory factor analysis result on stress scale

Fig. 3 shows the stress levels of participants before and after the experiment based on the stress scale. Before the experiment, the total stress levels in CG and PTG were 29.14±7.46 and 28.96± 8.77, respectively; however, that of CG increased to 29.36±4.53 (Z=−0.048, P=0.962), while that of PTG decreased to 22.70± 3.10 (Z=−2.921, P=0.003) after 8 weeks (Fig. 3E). The Δ% total stress levels of CG was 10.89±49.59%, while that of PTG was −9.31%±45.98%, indicating a significant difference between the groups (Z=−2.568; P=0.010, η2=0.210) (Fig. 3F).

Fig. 3

Differences and changes in stress level. In the figure, panels A and B represent stress scores analyzed before and after intervention in CG. Panels C and D represent stress scores analyzed before and after phytoncide intervention in PTG. Panel E is the comparison result of stress scores before and after intervention in CG and PTG, and panel F is the comparison result of the change rate of total stress levels in the two groups. CG, control group; PTG, phytoncide group. **P<0.01.

Implication of phytoncide fragrance on autonomic nervous activity

As shown in Fig. 4A, the SNS of CG was 7.16±0.85 before the experiment and increased to 7.38±1.00 (Z=−1.166, P=0.243) after 8 weeks. The SNS of PTG was 7.34±0.62 before intervention and decreased to 6.07±0.82 (Z=−3.965, P=0.001) after 8 weeks. The Δ% of CG was 3.82%±14.7%, while that of PTG was −16.75%±13.25% (P=0.001, η2=0.359). These results showed significant differences and changes in the group (P=0.002, η2= 0.166), time (P=0.001, η2=0.205), and interaction (P=0.001, η2=0.339).

Fig. 4

Differences and changes in sympathetic nerve activity and parasympathetic nerve activity. Panels A and B represent nerve activities analyzed before and after intervention in CG and PTG. CG, control group; PTG, phytoncide group. ***P<0.001.

The PNS of CG was 3.72±0.49 and decreased to 3.39±0.37 (Z=−3.737, P=0.001) after the experiment, as shown in Fig. 4B. Meanwhile, the PNS of PTG increased from 3.56±0.45 to 4.10± 0.47 (Z=−3.245, P=0.001) after 8 weeks. The Δ% of CG was −8.31%±10.05%, while that of PTG was 16.93%±20.42% (P= 0.001, η2=0.392). These results showed significant differences and changes in the group (P=0.004, η2=0.143) and interaction (P=0.001, η2=0.379).

The implication of phytoncide fragrance on stress hormones

As shown in Fig. 5A, the epinephrine of CG was 37.90±8.84 pg/mL before the experiment and increased to 48.68±13.25 pg/mL (Z=−3.131, P=0.002) after it. Conversely, the epinephrine of PTG was 39.03±8.53 pg/mL before intervention and decreased to 35.74±8.35 pg/mL (Z=−1.412, P=0.158) after 8 weeks. The Δ% of CG was 34.31%±46.46%, while that of PTG was −5.29%± 25.5% (P=0.001, η2=0.223). These results showed significant differences and changes in the group (P=0.005, η2=0.139), time (P=0.040, η2=0.077), and interaction (P=0.001, η2=0.228).

Fig. 5

Differences and changes in epinephrine and cortisol hormones. Panels A and B represent two endocrine concentrations analyzed before and after intervention in CG and PTG. CG, control group; PTG, phytoncide group. **P<0.01. ***P<0.001.

The cortisol of CG was 16.18±2.63 μg/dL before the experiment and increased to 17.36±3.19 μg/dL (Z=−1.139, P=0.255) after 8 weeks, as shown in Fig. 5B. In contrast, the cortisol of PTG was 15.96±2.22 μg/dL before intervention and decreased to 11.87± 2.01 μg/dL (Z=−4.474, P=0.001) after the experiment. The Δ% of CG was 9.0%±23.62%, while that of PTG was −24.94%± 11.62% (P=0.001, η2=0.460). These results showed significant differences and changes in the group (P=0.001, η2=0.330), time (P=0.001, η2=0.193), and interaction (P=0.001, η2=0.441).

The implication of phytoncide fragrance on immunocytes

Table 4 presents NK cells and their subsets showed significant positive changes except for CD56bright NK cell. Among them, KIR2DL3 showed an increased tendency in CG and a decreased change in PTG. The remaining NK cells, subsets, and receptors showed decreased changes in the CG and increased changes in the PTG. Regarding the rate of change, the total NK cell of CG was 0.61%±37.1%, while that of PTG was 23.32%±32.94% (P= 0.019, η2=0.098). Moreover, the CD56dim NK Δ% in CG was −2.43%±3.94%, while that of PTG was 1.45%±4.02% (P=0.001, η2=0.198). The NKG2D Δ% in CG was −1.97%±4.22%, while that of PTG was 3.02%±3.24% (P=0.001, η2=0.311). The KIR2DL3 Δ% in CG and PTG were 8.36%±21.82% and −7.56%± 19.87%, respectively, indicating a significant difference between groups (P=0.006, η2=0.131). The TIGIT Δ% in CG and PTG were −4.33%±23.6% and 23.26%±25.69%, respectively (P= 0.001, η2=0.245). The perforin Δ% in CG and PTG were −6.17%± 7.30% and 3.21%±4.84%, demonstrating a significant difference between groups (P=0.001, η2=0.372). Finally, the LFA-1+NK Δ% in CG and PTG were −3.56%±5.68% and 6.08%±10.02%, representing a significant difference between groups (P=0.001, η2=0.269).

Changes and differences of NK cells’ family and related receptors

DISCUSSION

This study confirmed that gynecological cancer survivors had moderate-to-high stress levels that decreased after phytoncide intervention. In addition, the study discovered that cancer survivors who were exposed to phytoncide fragrance experienced a beneficial impact on their stress hormones, as it suppressed the SNS and activated the PNS. Furthermore, a decrease in stress hormones was found to enhance NK cells and their subsets after being treated with phytoncide.

A fragrance is a volatile chemical compound that is perceived through olfaction (Ko et al., 2021; Woo and Lee, 2020; Zarzo, 2007). Olfactory receptors detect molecules and produce emotional responses that carry signals to the brain (Royet and Plailly, 2004). The main components of phytoncides are terpenes, which have rejuvenating effects in humans (Houdkova et al., 2020). Ikei et al. (2015) reported that essential oils extracted from phytoncides could be used to treat stress. Similarly, this study found that the stress level of the PTG decreased by approximately 9% after phytoncide intervention; this stress reduction was due to changes in the levels of neuroendocrine substances. Among the results of this study, the SNS activity was significantly decreased in the PTG (~−17%), while the PNS activity was significantly increased in the PTG (~17%). It is estimated that the effect of the phytoncide scent on stress is positive when it is absorbed into the body (Li et al., 2008). Phytoncide reduces stress by stimulating the olfactory bulb, which activates the autonomic nerve system (Angelucci et al., 2014). However, few studies have examined which hormones are affected by the autonomic nerve system and immunocytes (Segerstrom and Miller, 2004). Therefore, this study aimed to fill this gap, and it was found that phytoncide scent suppressed the SNS while activating the PNS. According to the results of this study, epinephrine and cortisol levels tended to decrease after applying phytoncide fragrance, though cortisol levels reduced more than epinephrine levels. Notably, the epinephrine concentration Δ% of PTG decreased by ~5%, while that of CG increased by ~34% after 8 weeks. The cortisol concentration Δ% of the PTG decreased by ~25%, while that of CG increased by ~9% after 8 weeks. Based on the effects, the phytoncide scent caused a positive change in the cortisol hormone in gynecological cancer survivors.

NK cells and their subsets were altered by stress hormones, as verified in this study, suggesting that the reduction of epinephrine and cortisol activates NK cells (Bodner et al., 1998; Tønnesen et al., 1987). Li et al. (2006) stated that the target cells are killed when phytoncides are used to treat an NK cell line. They also suggested that phytoncides increase NK cell activity by inhibiting cortisol. Moreover, NK cells activated by increased phytoncide levels are related to increased perforin, granzyme B, and granulin (Li et al., 2006). Thus, NK cells were influenced by their subsets. Specifically, the expression intensity of CD56 further divided NK cells into two subsets (CD56dim and CD56bright NK cells) (Phan et al., 2017). While CD56dim exhibited cytotoxicity, CD56bright showed less cytotoxicity but regulated immunocytes by secreting various cytokines (Raulet, 2003). Among these two subsets, this study found that CD56dim NK, rather than CD56bright NK, had a greater effect on the phytoncide intervention.

Typically, NKG2D+NK receptor detects UL16-binding proteins and major histocompatibility complex (MHC) class-I-chain-related proteins, which are intracellular molecules whose expression is increased during cancer (Martinez-Borra and Khakoo, 2008). In contrast, the inhibitory receptors of NK cells measure the presence or absence of intracellular molecules that are constitutively expressed on the surfaces of target cells (Long, 2008). Specifically, human KIR is an MHC class I-specific inhibitory receptor expressed in NK cells (Di Santo, 2008). This study showed that NKG2D+NK increased, while KIR2DL3+NK of the PTG decreased. These results suggest that phytoncides affect immunocyte receptors, even if provided as an aroma compound, as the blood flow of the chemical function in the tissues is improved. In addition, TIGIT, perforin, and LFA-1+NK cell levels decreased in the CG, while those increased in the PTG. Specifically, genetic ablation of TIGIT has been reported to enhance the NK cell-killing function (Urlaub et al., 2017). One of the most important weapons of killer T cells is perforin, a protein that pierces the outer membrane of target cells and causes cell necrosis (Trapani and Smyth, 2002). Moreover, LFA-1 was critical for the efficient activation and cytotoxicity of NK cells (Urlaub et al., 2017). As known, the subsets of NK cells are more mature and have higher cytotoxic potential to activate LFA-1 more.

This study observed that NK cells, including their subsets and receptors, were highly positively enhanced in gynecological cancer survivors. These results suggest that when the phytoncide scent is transmitted through olfactory cells, it stimulates the PNS and reduces stress hormones, thereby enhancing the number of NK family cells. Previous studies have reported that forest bathing enhances NK cell activity and expression of intracellular anticancer proteins (Akram et al., 2020; Eagle and Trowsdale, 2007; Gatti et al., 1987; Prager and Watzl, 2019). In addition, various chemicals emitted from forests have been reported to reduce psychological stress by activating the PNS while inhibiting the SNS (Chae et al., 2021; Lee et al., 2017; Li, 2010). This study also found that phytoncide essential oil was similar to the effect of phytoncide exposure in forests. Here, we speculate that the inhalation of phytoncide scent can prevent the decline of immunocytes and activate the function of innate immunocytes that destroy cancer cells. That is, a phytoncide fragrance from essential oil has been found to affect the human nervous and endocrine systems, which in turn affects the mobility of immunocytes and provides relief from psychological anxiety and stress for cancer survivors who previously had cancer cells. In summary, the fragrance of essential oil has been confirmed to significantly affect immune cell mobility and provide psychological stability. However, this study had the limitations of small sample size and limited diversity of participants in terms of demographic background. Consequently, we encourage future researchers to include a larger and more diverse pool of participants and observe a greater number of immunocytes.

ACKNOWLEDGMENTS

The authors wish to thank the participants of this study. The authors also would like to thank Doctor Jong-Kyun Lee for the clinical evaluation of survivors and for registration of clinical data. This work was supported by the Ministry of Education of the Republic of Korea and the National Research Foundation of Korea (NRF-2022S1A5A2A01048017).

Notes

CONFLICT OF INTEREST

No potential conflict of interest relevant to this article was reported.

References

Akram, M., Riaz, M., Munir, N., Rasul, A., Daniyal, M., Ali Shah, S.M., Shariati, M.A., Shaheen, G., Akhtar, N., Parveen, F., Akhter, N., Owais Ghauri, A., Chishti, A.W., Usman Sarwar, M., & Said Khan, F. Progress and prospects in the management of bacterial infections and developments in Phytotherapeutic modalities. Clin Exp Pharmacol Physiol, (2020). 47, 1107–1119.
Angelucci, F.L., Silva, V.V., Dal Pizzol, C., Spir, L.G., Praes, C.E., & Maibach, H. Physiological effect of olfactory stimuli inhalation in humans: an overview. Int J Cosmet Sci, (2014). 36, 117–123.
Benedet, J.L. Progress in gynecologic cancer detection and treatment. Int J Gynaecol Obstet, (2000). 70, 135–147.
Bodner, G., Ho, A., & Kreek, M.J. Effect of endogenous cortisol levels on natural killer cell activity in healthy humans. Brain Behav Immun, (1998). 12, 285–296.
Booth, M. Assessment of physical activity: an international perspective. Res Q Exerc Sport, (2000). 71, Suppl 2. 114–20.
Chae, Y., Lee, S., Jo, Y., Kang, S., Park, S., & Kang, H. The effects of forest therapy on immune function. Int J Environ Res Public Health, (2021). 18, 8440.
Cohen, J. A power primer. Psychol Bull, (1992). 112, 155–159.
Cohen, S., Kamarck, T., & Mermelstein, R. A global measure of perceived stress. J Health Soc Behav, (1983). 24, 385–396.
Di Santo, J.P. Natural killer cells: diversity in search of a niche. Nat Immunol, (2008). 9, 473–475.
Eagle, R.A., & Trowsdale, J. Promiscuity and the single receptor: NKG2D. Nat Rev Immunol, (2007). 7, 737–744.
Faul, F., Erdfelder, E., Lang, A.G., & Buchner, A. G*Power 3: a flexible statistical power analysis program for the social, behavioral, and biomedical sciences. Behav Res Methods, (2007). 39, 175–191.
Gatti, G., Cavallo, R., Sartori, M.L., del Ponte, D., Masera, R., Salvadori, A., Carignola, R., & Angeli, A. Inhibition by cortisol of human natural killer (NK) cell activity. J Steroid Biochem, (1987). 26, 49–58.
Golden-Kreutz, D.M., Browne, M.W., Frierson, G.M., & Andersen, B.L. Assessing stress in cancer patients: a second-order factor analysis model for the perceived stress scale. Assessment, (2004). 11, 216–223.
Houdkova, M., Albarico, G., Doskocil, I., Tauchen, J., Urbanova, K., Tulin, E.E., & Kokoska, L. Vapors of volatile plant-derived products significantly affect the results of antimicrobial, antioxidative and cytotoxicity microplate-based assays. Molecules, (2020). 25, 6004.
Ikei, H., Song, C., & Miyazaki, Y. Physiological effect of olfactory stimulation by Hinoki cypress (Chamaecyparis obtusa) leaf oil. J Physiol Anthropol, (2015). 34, 44.
Jo, H., Cha, B., Kim, H., Brito, S., Kwak, B.M., Kim, S.T., Bin, B.H., & Lee, M.G. α-Pinene enhances the anticancer activity of natural killer cells via ERK/AKT pathway. Int J Mol Sci, (2021). 22, 656.
Ko, L.W., Su, C.H., Yang, M.H., Liu, S.Y., & Su, T.P. A pilot study on essential oil aroma stimulation for enhancing slow-wave EEG in sleeping brain. Sci Rep, (2021). 11, 1078.
Lane, R.D., McRae, K., Reiman, E.M., Chen, K., Ahern, G.L., & Thayer, J.F. Neural correlates of heart rate variability during emotion. Neuroimage, (2009). 44, 213–222.
Lee, G.S., Hong, E.J., Gwak, K.S., Park, M.J., Choi, K.C., Choi, I.G., Jang, J.W., & Jeung, E.B. The essential oils of Chamaecyparis obtusa promote hair growth through the induction of vascular endothelial growth factor gene. Fitoterapia, (2010). 81, 17–24.
Lee, I., Choi, H., Bang, K.S., Kim, S., Song, M., & Lee, B. Effects of forest therapy on depressive symptoms among adults: a systematic review. Int J Environ Res Public Health, (2017). 14, 321.
Lee, S.H., Do, H.S., & Min, K.J. Effects of essential oil from Hinoki cypress, Chamaecyparis obtusa, on physiology and behavior of flies. PLoS One, (2015). 10, e0143450.
Li, Q. Effect of forest bathing trips on human immune function. Environ Health Prev Med, (2010). 15, 9–17.
Li, Q., Morimoto, K., Kobayashi, M., Inagaki, H., Katsumata, M., Hirata, Y., Hirata, K., Suzuki, H., Li, Y.J., Wakayama, Y., Kawada, T., Park, B.J., Ohira, T., Matsui, N., Kagawa, T., Miyazaki, Y., & Krensky, A.M. Visiting a forest, but not a city, increases human natural killer activity and expression of anti-cancer proteins. Int J Immunopathol Pharmacol, (2008). 21, 117–127.
Li, Q., Morimoto, K., Nakadai, A., Inagaki, H., Katsumata, M., Shimizu, T., Hirata, Y., Hirata, K., Suzuki, H., Miyazaki, Y., Kagawa, T., Koyama, Y., Ohira, T., Takayama, N., Krensky, A.M., & Kawada, T. Forest bathing enhances human natural killer activity and expression of anti-cancer proteins. Int J Immunopathol Pharmacol, (2007). 20, (2 Suppl 2)3–8.
Li, Q., Nakadai, A., Matsushima, H., Miyazaki, Y., Krensky, A.M., Kawada, T., & Morimoto, K. Phytoncides (wood essential oils) induce human natural killer cell activity. Immunopharmacol Immunotoxicol, (2006). 28, 319–333.
Long, E.O. Negative signaling by inhibitory receptors: the NK cell paradigm. Immunol Rev, (2008). 224, 70–84.
Lundberg, F.E., Kroman, N., Lambe, M., Andersson, T.M., Engholm, G., Johannesen, T.B., Virtanen, A., Pettersson, D., Ólafsdóttir, E.J., Birgisson, H., Lambert, P.C., Mørch, L.S., & Johansson, A.L.V. Age-specific survival trends and life-years lost in women with breast cancer 1990–2016: the NORDCAN survival studies. Acta Oncol, (2022). 61, 1481–1489.
Martinez-Borra, J., & Khakoo, S.I. Speed and selection in the evolution of killer-cell immunoglobulin-like receptors. Int J Immunogenet, (2008). 35, 89–96.
Neeman, E., & Ben-Eliyahu, S. Surgery and stress promote cancer metastasis: new outlooks on perioperative mediating mechanisms and immune involvement. Brain Behav Immun, (2013). 30, (Suppl)S32–S40.
Phan, M.T., Chun, S., Kim, S.H., Ali, A.K., Lee, S.H., Kim, S., Kim, S.H., & Cho, D. Natural killer cell subsets and receptor expression in peripheral blood mononuclear cells of a healthy Korean population: reference range, influence of age and sex, and correlation between NK cell receptors and cytotoxicity. Hum Immunol, (2017). 78, 103–112.
Prager, I., & Watzl, C. Mechanisms of natural killer cell-mediated cellular cytotoxicity. J Leukoc Biol, (2019). 105, 1319–1329.
Raha, S., Kim, S.M., Lee, H.J., Lee, S.J., Heo, J.D., Venkatarame Gowda Saralamma, V., Ha, S.E., Kim, E.H., Mun, S.P., & Kim, G.S. Essential oil from Korean Chamaecyparis obtusa leaf ameliorates respiratory activity in Sprague-Dawley rats and exhibits protection from NF-κB-induced inflammation in WI38 fibroblast cells. Int J Mol Med, (2019). 43, 393–403.
Raulet, D.H. Roles of the NKG2D immunoreceptor and its ligands. Nat Rev Immunol, (2003). 3, 781–790.
Royet, J.P., & Plailly, J. Lateralization of olfactory processes. Chem Senses, (2004). 29, 731–745.
Salleh, M.R. Life event, stress and illness. Malays J Med Sci, (2008). 15, 9–18.
Segerstrom, S.C., & Miller, G.E. Psychological stress and the human immune system: a meta-analytic study of 30 years of inquiry. Psychol Bull, (2004). 130, 601–630.
Tang, P.L., Mayer, D.K., Chou, F.H., & Hsiao, K.Y. The experience of cancer stigma in Taiwan: a qualitative study of female cancer patients. Arch Psychiatr Nurs, (2016). 30, 204–209.
The Korean Nutrition Society. CAN-Pro 5.0 [Internet], Seoul (Korea): The Korean Nutrition Society; (2023). [cited 2023 Jan 3]. Available from: https://www.kns.or.kr/can/CanPro5_result.asp .
Tønnesen, E., Christensen, N.J., & Brinkløv, M.M. Natural killer cell activity during cortisol and adrenaline infusion in healthy volunteers. Eur J Clin Invest, (1987). 17, 497–503.
Trapani, J.A., & Smyth, M.J. Functional significance of the perforin/granzyme cell death pathway. Nat Rev Immunol, (2002). 2, 735–747.
Urbaniak, G.C., & Plous, S. Research Randomizer (Version 4.0) [Internet], Social Psychology Network; (c1997–2023). [cited 2023 Jan 3]. Available from: http://www.randomizer.org/ .
Urlaub, D., Höfer, K., Müller, M.L., & Watzl, C. LFA-1 activation in NK cells and their subsets: influence of receptors, maturation, and cytokine stimulation. J Immunol, (2017). 198, 1944–1951.
Vivier, E., Tomasello, E., Baratin, M., Walzer, T., & Ugolini, S. Functions of natural killer cells. Nat Immunol, (2008). 9, 503–510.
White, D.W., & Raven, P.B. Autonomic neural control of heart rate during dynamic exercise: revisited. J Physiol, (2014). 592, 2491–2500.
Woo, J., & Lee, C.J. Sleep-enhancing effects of phytoncide via behavioral, electrophysiological, and molecular modeling approaches. Exp Neurobiol, (2020). 29, 120–129.
Zarzo, M. The sense of smell: molecular basis of odorant recognition. Biol Rev Camb Philos Soc, (2007). 82, 455–479.
Zhang, M., Cheng, S., Jin, Y., Zhao, Y., & Wang, Y. Roles of CA125 in diagnosis, prediction, and oncogenesis of ovarian cancer. Biochim Biophys Acta Rev Cancer, (2021). 1875, 188503.

Article information Continued

Fig. 1

Participants’ allocation (consolidated standards for reporting of trials flow diagram).

Fig. 2

Confirmatory factor analysis model.

Fig. 3

Differences and changes in stress level. In the figure, panels A and B represent stress scores analyzed before and after intervention in CG. Panels C and D represent stress scores analyzed before and after phytoncide intervention in PTG. Panel E is the comparison result of stress scores before and after intervention in CG and PTG, and panel F is the comparison result of the change rate of total stress levels in the two groups. CG, control group; PTG, phytoncide group. **P<0.01.

Fig. 4

Differences and changes in sympathetic nerve activity and parasympathetic nerve activity. Panels A and B represent nerve activities analyzed before and after intervention in CG and PTG. CG, control group; PTG, phytoncide group. ***P<0.001.

Fig. 5

Differences and changes in epinephrine and cortisol hormones. Panels A and B represent two endocrine concentrations analyzed before and after intervention in CG and PTG. CG, control group; PTG, phytoncide group. **P<0.01. ***P<0.001.

Table 1

Perceived stress scale for gynecological cancer survivor

Items Never Almost never Some times Fairly often Very often
Q1 How often have you been upset because of something that happened unexpectedly? 0 1 2 3 4
Q2 How often have you felt that you were unable to control the important things in your life? 0 1 2 3 4
Q3 How often have you felt nervous and stressed? 0 1 2 3 4
Q4 How often have you felt confident about your ability to handle your personal problems? 0 1 2 3 4
Q5 How often have you felt that things were going your way? 0 1 2 3 4
Q6 How often have you found that you could not cope with all the things that you had to do? 0 1 2 3 4
Q7 How often have you been able to control irritations in your life? 0 1 2 3 4
Q8 How often have you felt that you were on top of things? 0 1 2 3 4
Q9 How often have you been angered because of things that happened that were outside of your control? 0 1 2 3 4
Q10 How often have you felt difficulties were piling up so high that you could not overcome them? 0 1 2 3 4

Questions on this scale ask about the cancer survivors’ feelings and thoughts in the last month. For each occasion or situation, indicate how often you felt or thought a certain way. Some questions are similar but have differences, so each should be treated as a separate question.

Table 2

Demographic and clinical characteristics of the participants

Variable Groups Z P-value η2
CG (n=28) PTG (n=27) Total (n=55)
Age (yr) 69.04±5.12 68.67±5.94 68.85±5.49 −0.051 0.960 0.001
Height (cm) 162.71±6.75 163.26±7.16 162.98±6.90 −0.076 0.939 0.002
Body weight (kg) 68.01±11.08 69.14±8.78 68.57±9.94 −0.304 0.761 0.003
Body mass index (kg/m2) 25.60±3.37 26.04±3.80 25.82±3.56 −0.712 0.477 0.004
Muscle mass (kg) 27.51±5.53 28.70±5.79 28.09±5.64 −0.820 0.412 0.011
Fat mass (kg) 18.84±5.17 19.44±4.67 19.14±4.90 −0.668 0.504 0.004
Percent fat (%) 27.85±6.55 27.80±6.90 27.83±6.66 −0.228 0.819 0.001
Married 1.14±0.36 1.19±0.40 1.16±0.37 −0.420 0.674 0.003
Employed 1.68±0.48 1.85±0.36 1.76±0.43 −1.498 0.134 0.042
Medications’ number 2.29±1.18 2.37±0.88 2.33±1.04 −0.402 0.688 0.002
Supplements’ number 1.82±0.55 1.93±0.62 1.87±0.58 −0.631 0.528 0.008
Comorbidities’ number 2.82±0.86 2.89±0.89 2.85±0.87 −0.400 0.689 0.002
Cancer types 1.71±0.76 1.56±0.70 1.64±0.73 −0.779 0.436 0.012
Cancer diagnosed years 4.07±0.48 4.03±0.53 4.05±0.50 −0.712 0.476 0.002
Lymph node metastasis 1.79±0.42 1.70±0.47 1.75±0.44 −0.692 0.489 0.009
Organ metastasis 1.93±0.26 1.85±0.36 1.89±0.31 −0.904 0.366 0.015
Surgical stages 1.21±0.42 1.15±0.36 1.18±0.39 −0.630 0.529 0.007
Stages 0 or FIGO I 3.50±1.87 4.50±2.45 4.07±2.20 −0.780 0.436 0.054
Stages I or FIGO II 8.00±4.47 6.50±3.61 7.33±4.10 −0.880 0.379 0.034
Stages II–III or FIGO IIIA 2.00±1.00 3.00±1.58 2.63±1.41 −0.911 0.362 0.135
Stages IV or FIGO IIIB 2.50±1.29 1.50±0.71 2.17±1.17 −0.953 0.340 0.195
Radiation therapy 1.29±0.46 1.41±0.50 1.35±0.48 −0.940 0.347 0.016
Chemotherapy 1.11±0.31 1.15±0.36 1.13±0.34 −0.452 0.651 0.004

Values are presented mean±standard deviation.

CG, control group; PTG, phytoncide group; FIGO, International Federation of Gynecology and Obstetrics.

Table 3

Confirmatory factor analysis result on stress scale

Factors Question B SE β CR P-value
Stress Q1 1 - 0.427 - -
Stress Q2 1.335 0.629 0.684 2.124 0.034
Stress Q3 1.053 0.566 0.510 1.861 0.063
Stress Q4 1.976 0.851 0.910 2.322 0.020
Stress Q5 0.910 0.525 0.449 1.735 0.083
Stress Q6 1.913 0.829 0.889 2.309 0.021
Stress Q7 1.117 0.573 0.560 1.957 0.050
Stress Q8 1.397 0.668 0.657 2.091 0.037
Stress Q9 1.188 0.597 0.176 1.991 0.046
Stress Q10 2.024 0.873 0.903 2.318 0.020

Q, question; SE, standard error; CR, concept reliability.

Table 4

Changes and differences of NK cells’ family and related receptors

Variable Times (T) Groups (G) P-value η2


CG PTG G T G*T
Total NK (cells/μL) Pre 524.18±240.83 520.74±227.17 0.331 0.415 0.009 0.122
Post 486.93±186.01 590.41±140.71

Total NK (%) Pre 25.35±10.45 24.40±6.87 0.270 0.454 0.001 0.253
Post 22.71±7.93 28.18±6.11

CD56bright NK (%) Pre 2.22±1.12 2.12±0.84 0.295 0.305 0.283 0.022
Post 1.95±0.96 14.28±60.90

CD56dim NK (%) Pre 87.99±4.11 88.33±4.90 0.119 0.349 0.001 0.197
Post 85.87±5.62 89.57±5.38

NKG2D (%) Pre 88.50±3.72 88.54±4.46 0.053 0.327 0.001 0.316
Post 86.71±4.43 91.21±5.50

KIR2DL3 (%) Pre 38.99±10.26 39.99±10.42 0.328 0.393 0.005 0.140
Post 41.29±9.39 35.79±7.16

TIGIT (%) Pre 44.64±13.64 44.21±10.87 0.025 0.113 0.001 0.254
Post 40.56±8.83 53.27±11.89

Perforin (%) Pre 92.13±3.49 91.23±3.87 0.001 0.058 0.001 0.370
Post 86.37±6.54 94.01±2.79

LFA-1+NK (%) Pre 83.69±6.86 83.08±7.97 0.110 0.338 0.001 0.265
Post 80.94±9.80 87.39±2.34

Values are presented mean±standard deviation.

NK, natural killer; CG, control group; PTG, phytoncide inhales group; NKG2D, natural-killer group 2 member D; KIR2DL3, killer cell immunoglobulin (Ig) like receptor, 2 Ig domains and long cytoplasmic tail 3; TIGIT, T cell immunoreceptor with Ig and ITIM (immunoreceptor tyrosine-based inhibitory motif) domains; LFA-1, lymphocyte function-associated antigen 1.