Affectionate touch and diurnal oxytocin levels: An ecological momentary assessment study

  1. Ekaterina Schneider
  2. Dora Hopf
  3. Corina Aguilar-Raab
  4. Dirk Scheele
  5. Andreas B Neubauer
  6. Uta Sailer
  7. René Hurlemann
  8. Monika Eckstein  Is a corresponding author
  9. Beate Ditzen  Is a corresponding author
  1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Germany
  2. Heidelberg University, Germany
  3. Department of Social Neuroscience, Faculty of Psychology, Ruhr University Bochum, Germany
  4. Department for Education and Human Development, DIPF|Leibniz Institute for Research and Information in Education, Germany
  5. Center for Research on Individual Development and Adaptive Education of Children at Risk, Germany
  6. Department of Behavioural Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Norway
  7. Department of Psychiatry, University of Oldenburg, Germany

Abstract

Background:

Affectionate touch, which is vital for mental and physical health, was restricted during the Covid-19 pandemic. This study investigated the association between momentary affectionate touch and subjective well-being, as well as salivary oxytocin and cortisol in everyday life during the pandemic.

Methods:

In the first step, we measured anxiety and depression symptoms, loneliness and attitudes toward social touch in a large cross-sectional online survey (N = 1050). From this sample, N = 247 participants completed ecological momentary assessments over 2 days with six daily assessments by answering smartphone-based questions on affectionate touch and momentary mental state, and providing concomitant saliva samples for cortisol and oxytocin assessment.

Results:

Multilevel models showed that on a within-person level, affectionate touch was associated with decreased self-reported anxiety, general burden, stress, and increased oxytocin levels. On a between-person level, affectionate touch was associated with decreased cortisol levels and higher happiness. Moreover, individuals with a positive attitude toward social touch experiencing loneliness reported more mental health problems.

Conclusions:

Our results suggest that affectionate touch is linked to higher endogenous oxytocin in times of pandemic and lockdown and might buffer stress on a subjective and hormonal level. These findings might have implications for preventing mental burden during social contact restrictions.

Funding:

The study was funded by the German Research Foundation, the German Psychological Society, and German Academic Exchange Service.

Editor's evaluation

This important study combines a large cross-sectional survey with detailed ecological momentary assessment to examine the relationship between affectionate touch and well-being during the first wave of the COVID-19 pandemic in Germany. The manuscript reports valuable and solid findings extending previous research in this domain. Specifically, the combination of ecologically momentary assessment data with repeated measurements of salivary cortisol and oxytocin adds to the current understanding of how affectionate touch relates to psychological burden and affect. Due to its correlational nature, the causality of effects remains speculative and needs to be addressed by future work.

https://doi.org/10.7554/eLife.81241.sa0

Introduction

Social integration and close social contact have been shown to improve mental and physical health as well as increase longevity (Holt-Lunstad, 2018). This effect has been suggested to be mediated through physical proximity and affectionate touch, with touch serving as a social safety signal (Eckstein et al., 2020). Affectionate touch has been associated with beneficial effects on human development and psychological well-being throughout the lifespan (Atzil et al., 2018; Cascio et al., 2019). Touch activates reward-related brain regions (Kreuder et al., 2017) and reduces stress-induced cortisol (Ditzen et al., 2007) and pain (Kreuder et al., 2019). On a neuroendocrine level, the stress-buffering effects of affectionate touch on subjective measures and activity of the hypothalamic–pituitary–adrenal (HPA) axis have been hypothesized to be mediated by the neuropeptide hormone oxytocin (Eckstein et al., 2020).

The outbreak of the Covid-19 pandemic was a continuous stressor with major health and societal consequences (Fancourt et al., 2021; Mata et al., 2021; Pierce et al., 2020). Immediate restrictions and physical distancing were necessary measures to control the spread of the virus. The resulting physical isolation has been linked to higher self-reported loneliness, especially as a result of the first lockdown (Fancourt et al., 2021; Mata et al., 2021; Pierce et al., 2020) and in individuals with previous higher loneliness (Bu et al., 2020). In general, loneliness and social isolation have been associated with poorer mental and physical health as well as increased mortality (Lee et al., 2021; Leigh-Hunt et al., 2017). Thus, it is not surprising that several recent studies emphasize the potential impact of loneliness during the Covid-19 pandemic on mental health (Brooks et al., 2020; Campion et al., 2020; Rozenkrantz et al., 2020). Large population-based studies suggest that levels of mental distress with clinical significance increased from 18.9% in 2018–19 to 27.3% during the pandemic (Pierce et al., 2020), and the increase was most prominent in the initial phase of the lockdown (Chandola et al., 2022). On the other hand, perceived social support as well as frequent social contact during the pandemic were associated with lower depression scores (Sommerlad et al., 2021).

The request to minimize social contact and increase physical distance during lockdown consequently led to less physical contact and lower frequency of interpersonal touch, potentially increasing the feeling of longing for touch. Moreover, higher longing for touch was associated with prolonged and more severe Covid-19 restrictions (Meijer et al., 2022). Literature on touch deprivation suggests that a lack of touch is associated with lower levels of general well-being and an increased risk of mental health problems (Banerjee et al., 2021). A recent study by von Mohr and colleagues showed that self-reported deprivation of intimate touch (but not other types such as friendly or professional touch) during the Covid-19 lockdown was associated with higher loneliness scores. In addition, they found that intimate touch deprivation was associated with higher anxiety levels; however, this association was no longer significant when accounting for loneliness (von Mohr et al., 2021). The authors suggested that the lack of intimate touch may increase anxiety in individuals with higher loneliness. Burleson and colleagues reported that reduced affectionate touch was associated with more psychological distress, especially for those participants, who typically use touch for affect regulation (Burleson et al., 2022).

Initial laboratory research has demonstrated that receiving touch such as a massage can have beneficial effects evident in reduced self-reported anxiety and stress levels (Kirschner and Kirschner, 2019), as well as decreased cortisol (Maratos et al., 2017) and increased oxytocin (Morhenn et al., 2012) concentrations. Similarly, a more recent study found a significant increase in plasma oxytocin and corresponding neural responses after a foot massage. Interestingly, basal oxytocin concentrations, as well as oxytocin increase after the massage were associated with more positive attitudes toward social touch (Li et al., 2019). Moreover, touching a dog as compared to merely observing it was associated with not only decreased self-reported stress, but also increased self-reported happiness (Sokal et al., 2021). Based on these findings, we hypothesized that a positive attitude toward touch and increased loneliness would be associated with higher anxiety and depression symptoms during the lockdown. On a momentary level, we expected that affectionate touch would be associated with decreased subjective anxiety, distress, and decreased HPA axis activity (cortisol levels), as well as with higher endogenous oxytocin levels. Furthermore, we expected that the link between subjective anxiety and distress with affectionate touch would be mediated by elevated oxytocin levels. To the best of our knowledge, there has not yet been a study investigating the associations of affectionate touch with mental health and neuroendocrine variables during the Covid-19 lockdown. We addressed this gap using ecological momentary assessment (EMA) of both repeated psychological and endocrine measures in a large sample with frequent repeated everyday life measuring.

Methods

For this study, ethical approval was granted from the ethics committee of the Heidelberg University Medical Faculty (approval no. S-214/2020), and the study was registered online at https://drks.de/search/en/trial/DRKS00021671. All participants provided written informed consent. We used the disclosure of interest form of the International Committee of Medical Journal Editors (ICMJE) to report no conflicts of interest. STROBE protocol was used to standardize reporting.

Study design and population

In a large online survey launched in April 2020, structural social factors, such as housing situation, anxiety, and depressive symptoms, as well as subjective psychosocial burden, loneliness, and the perception of touch during the physical distancing measures, were assessed (Hopf et al., 2022). Study participants were recruited via local newspapers, radio programs, and social media. In an attempt to more actively involve the study participants in the research (collecting data, carrying out measurements in open formats, reporting unexpected results, i.e. citizen science approach), all participants (N = 1050) who had completed the online survey were invited to take part in a 2-day psychobiological EMA. Participants were given standardized instructions via phone on how to use their smartphones to collect momentary subjective data, as well as saliva samples via a passive drool method at six time points per day over the course of two consecutive days (i.e. in total, each individual provided 12 saliva samples). They received the collecting devices via mail along with the informed consent documents to sign. Sampling times on each day were adapted to the individual wake-up time and were taken directly after awakening, 30 min after, 45 min after, 2½ hr after, 8 hr after, and directly before going to sleep. To reduce potential missing values, minimize irregularities, and increase adherence, the data sampling was monitored by study members.

Measures

Hospital Anxiety and Depression Scale (HADS)

General psychological distress was assessed using the total score of the Hospital Anxiety and Depression Scale (HADS) (Hinz and Brähler, 2011). Sum scores were calculated for anxiety and depression subscales as well as for the total score. The internal consistency of the global HADS score in our data was high (HADS total score: Cronbach’s α = 0.89; HADS Anxiety subscale: Cronbach’s α = 0.82; HADS Depression subscale: Cronbach’s α = 0.82).

UCLA Loneliness Scale

Loneliness was measured using the 20-item UCLA Loneliness Scale (Döring and Bortz, 1993). Participants rated how often they felt in a certain way during the past 2 wk, with higher scores indicating higher levels of loneliness. The sum scores were used for statistical analyses. In our sample, the scale showed a high internal consistency (Cronbach’s α = 0.91).

Social Touch Questionnaire (STQ)

To measure attitudes toward social touch, we used the Social Touch Questionnaire (STQ) (Wilhelm et al., 2001), assessing different aspects of social touch such as touch involving family and friends vs. touch involving strangers, touch occurring in different settings, as well as touch with sexual vs. without sexual connotation. Internal consistency in our data was high with Cronbach’s α = 0.84. Low values of STQ indicate a high liking of social touch, whereas high values indicate a high aversion to social touch. To interpret the results more intuitively, individual scores of the STQ were inverted (i.e. high STQ values indicate a more positive attitude towards touch).

Ecological momentary assessment

Momentary levels of well-being (anxiety, stress, general and Covid-19 related burden, as well as happiness levels) were assessed through single items (‘Please indicate how you feel at the moment …’) using visual analog scales from 0 (not at all) to 100 (very much). Affectionate touch was assessed with the question ‘Since the last time point, did you experience touch, hugs, kisses, cuddles, etc.?’ and an additional visual analog scale for the intensity rating of the experienced touch from 0 (low intensity) to 100 (high intensity).

Neuroendocrine measures

On seeing the prompt on their smartphones, participants self-sampled their saliva into Salicaps (small plastic tubes) via passive drool and stored each sample immediately after collection in their home freezers. At the end of data collection, the study team personally visited to collect the samples on dry ice. The saliva samples were stored at –80°C until analyses at the Institute of Medical Psychology’s biochemical lab at Heidelberg University Hospital.

For the analyses of endogenous oxytocin concentrations, saliva samples were thawed and centrifuged at 4°C at 1.500 × g for 15 min and subsequently analyzed without extraction (50% of the samples in duplicates) following the protocol of oxytocin enzyme-linked immunosorbent assay from Enzo Life Sciences (ELISA; ENZO Life Sciences, Switzerland). The detection limit was 15 pg/ml, and the variation coefficient for intra- and inter-assay precision was 6.12 and 11.13%, respectively. For cortisol analyses, 20% of the samples were analyzed in duplicates and an ELISA from Demeditec Diagnostics (Demeditec Diagnostics, Germany) was used with a reported detection limit of 0.019 ng/ml. Intra- and inter-assay variations in our sample were 2.95 and 7.51%, respectively.

Statistical analyses

For data processing, IBM SPSS version 27 was used. Statistical analyses were conducted using R studio (R version 4.1.1) and Mplus (version 8.6). We analyzed the relationship of attitudes toward social touch (STQ) and loneliness (UCLA Loneliness) with anxiety and depression symptoms (HADS total) controlling for age, sex, and presence of mental disorder using multiple regression analyses. No violations of general assumptions for multiple regression (linearity, homoscedasticity, normality, and independence of errors) were detected. The total score of HADS, as well as HADS Anxiety and HADS Depression subscales, were included as dependent variables, whereas STQ, UCLA Loneliness, as well as the interaction variable (STQ × UCLA Loneliness) were entered as independent variables into the model. STQ and UCLA Loneliness scores were centered around their respective means. Missing data were deleted listwise.

To test whether affectionate touch was associated with well-being and neuroendocrine markers in everyday life, we conducted multiple hierarchical linear models. To separate within-person and between-person effects, self-reported touch (yes/no) and the intensity of touch were centered around each person’s mean and the person’s mean was centered on the grand mean. First, we included momentary affectionate touch (yes/no) controlling for age, sex, and day as independent variables to predict individual momentary self-reported anxiety, stress, general and Covid-19 related burden, as well as happiness levels in separate models. Subsequently, we analyzed whether the intensity of experienced touch was associated with these momentary psychological states following the same analytical approach. For models including cortisol and oxytocin measures as dependent variables, we additionally controlled for body mass index (BMI) and several potential confounders: momentary food and drink intake, alcohol, caffeine, and cigarette consumption, as well as physical activity, sleep duration, and quality, problems falling asleep, intake of sleeping pills, forced awakening, and brushing teeth. Furthermore, we controlled for assessment time points by including time (coded from 0 to 3 for the assessment time points 3–6) to control for linear diurnal changes after the awakening response (Ning and Luo, 2017). Additionally, for these models, we conducted random slope models and compared the fit of these models to models without random slope for the focal predictor (touch; touch intensity) using likelihood ratio tests. For the models on affectionate touch as a binary variable (yes/no), we report random intercept and random slopes models in the ‘Results’ section since these showed a statistically better fit compared to random intercept and fixed slopes models. However, for the dimensional intensity of affectionate touch, we report random intercept and fixed slopes models since the random slope models did not yield a better model fit. Before analyses, cortisol and oxytocin levels were log-transformed (natural logarithm) to normalize the distribution. Cortisol and oxytocin awakening response was calculated using the formula for calculating the area under the curve concerning increase (Pruessner et al., 2003).

Results

Sample characteristics

From April to August 2020, 1483 participants filled out the online survey, of whom 433 were excluded from data analysis (see Figure 1). A total of 1050 participants (n = 815 women, n = 227 men, n = 4 non-binary, n = 4 no information on gender) were included in the analyses. Participants’ age ranged from 18 to 81 y, with a mean age of 36.34 (SD = 14.77). 20.2% (n = 212) indicated that they suffered from a diagnosed mental disorder. The most frequent single diagnosis was depression (35%) followed by anxiety disorders (10%). Of those with at least one diagnosis, 27.5% indicated having multiple diagnoses.

Flowchart of the recruitment process.

Figure 1 depicts the recruitment stages of both the online and the ecologically momentary assessments (EMA) study. Participants were recruited between April 1 and July 30, 2020, via online media and local newspapers. Inclusion criteria: fluency in German, minimum age of 18 y, and willingness to participate voluntarily. In total, 1483 individuals agreed to participate, of whom 1050 participants filled out the online questionnaires of interest. Out of the 472 participants who were interested in the EMA study, 247 finished the assessments.

After completion of the online survey, 472 individuals indicated that they were interested in the EMA study, of whom 257 confirmed their participation after receiving detailed information. Ten participants withdrew from the study due to personal reasons, resulting in a total of 247 participants (n = 173 women, n = 74 men) completing the 2 d EMA. The mean age of the sample was 32.02 y (SD = 13.12) ranging from 18 to 78 y (for more details on sample characteristics, please see Table 1).

Table 1
Sample characteristics of online survey and ecological momentary assessment.
Sample characteristics of online survey participants
Men (n = 227)Women (n = 815)Non-binary (n = 4)Missing (n = 4)
M (SD)M (SD)M (SD)M (SD)
Age (years)34.67 (15.18)36.74 (14.58)45.50 (24.73)40.25 (15.39)
General psychological distress*10.22 (6.75)13.18 (7.49)21.00 (5.89)19.50 (14.66)
Anxiety5.15 (3.68)6.95 (4.12)10.50 (3.11)9.25 (7.68)
Depression 5.07 (3.55)6.22 (4.10)10.50 (3.32)10.25 (8.18)
Loneliness §37.18 (10.15)39.33 (10.95)53.00 (13.24)47.50 (19.50)
Attitude toward social touch 33.58 (10.18)34.76 (12.27)38.50 (23.39)40.75 (16.92)
Sample characteristics of ecological momentary assessment participants
Men (n = 74)Women (n = 173)
MSDMSD
Age (years)30.9913.6233.0512.41
Cortisol (ng/ml) **8.402.028.682.31
Oxytocin (pg/ml)**176.12106.18164.5496.74
Covid-19-related burden**36.9824.6141.7823.49
General burden39.9623.8947.2021.78
Stress levels††29.4915.8835.6217.00
Anxiety levels††18.3915.9824.1420.08
Happiness levels††71.1317.0967.8718.42
Intensity of affectionate touch††65.2120.0056.5723.13
  1. Table depicts means (M) and standard deviations (SD). Number of participants indicated as (n).

  2. *

    Hospital Anxiety and Depression Scale (HADS total score).

  3. HADS Anxiety subscale.

  4. HADS Depression subscale.

  5. §

    University of California, Los Angeles Loneliness Scale (UCLA Loneliness).

  6. Social Touch Questionnaire (STQ).

  7. **

    Out of 2964 possible data points, n = 2724 remained for analysis after excluding outliers, samples that were not stored as instructed or below detection limit, sampling problems.

  8. ††

    Momentary self-reported state.

Attitude toward touch and its association with anxiety, depression, and loneliness

On average, participants’ HADS total scores were M = 12.58 (SD = 7.49, range = 0–37). 39.7% of the sample had values above the cut-off score (>13) compared to a reference sample (Hinz and Brähler, 2011). The average HADS Anxiety subscale score was M = 6.58 (SD = 4.19, range = 0–20), whereas the HADS Depression subscale score was M = 5.99 (SD = 4.04, range = 0–21) with values exceeding the cut-off scores (>8) in 29.8 and 24.5% of cases, respectively. The results of multiple regression analyses showed significant main effects of sex (β = 0.111; t(1031) = 4.655, p<0.001), presence of diagnosed mental disorder (β = 0.151; t(1031) = 5.882, p<0.001), UCLA Loneliness (β = 0.548; t(1031) = 20.403, p<0.001), STQ (β = -0.052; t(1031) = −2.083, p=0.038), as well as a significant interaction of UCLA Loneliness × STQ (β = 0.052; t(1031) = 2.104, p=0.036) on the outcome variable total HADS score. Thus, anxiety and depression symptoms were higher in women, individuals with a mental disorder and participants with higher loneliness; and lower in participants with a more positive attitude toward touch. In contrast, although the moderation effects were small, they indicate that the association of loneliness with anxiety and depression symptoms was more pronounced for individuals with a more positive attitude toward social touch. The model tested here was significant overall (F(6,1031) = 125.1, p<0.001) with an R² of 0.421.

Next, we analyzed the association of the UCLA Loneliness × STQ interaction with the subscales of the HADS by following the same analytical approach. Here, we found that the outcome variable HADS Anxiety was also significantly and positively associated with female sex (β = 0.142; t(1031) = 5.369, p<0.001), presence of mental disorder (β = 0.180; t(1031) = 6.3, p<0.001), and UCLA Loneliness (β = 0.404; t(1031) = 13.54, p<0.001). Furthermore, we observed a significant interaction of UCLA Loneliness × STQ (β = 0.079; t(1031) = 2.907, p=0.004). However, the subscale HADS Depression showed only a significant association with sex (β = 0.060; t(1031) = 2.593, p=0.010), presence of mental disorder (β = 0.097; t(1031) = 3.838, p<0.001), and UCLA Loneliness (β = 0.603; t(1031) = 22.974, p<0.001). The UCLA Loneliness × STQ interaction was not significant (p=0.539). Both models with the outcome variable HADS Anxiety as well as with HADS Depression were overall significant (F(6,1031) = 69.25, p<0.001; with an R² of 0.287 and F(6,1031) = 139.1, p<0.001; with an R² of 0.447, respectively).

Affectionate touch, anxiety, oxytocin, and stress-related outcomes on a momentary level

Descriptive statistics of outcomes of interest are displayed in Table 1. In addition, an explorative graphical illustration of daily profiles of oxytocin and cortisol shows their variation throughout the day (Figure 2). The patterns of daily profiles did not appear to differ based on participants’ relationship status (single vs. in a relationship) or living situation (alone vs. with others) (see Figure 2—figure supplement 1). A positive correlation between the two assessment days was found for individual (ln-transformed) mean values of oxytocin (r(227) = 0.850, p<0.001), as well as cortisol (r(243) = 0.571, p<0.001) levels. Additionally, we found a significant negative correlation between mean cortisol and oxytocin awakening response (r(181) = –0.195, p=0.008).

Figure 2 with 1 supplement see all
Diurnal oxytocin and cortisol trajectories.

Panels (A) and (B) illustrate the daily oxytocin (pg/ml) and cortisol (ng/ml) trajectories across 2 d and all participants. Gray area indicates cortisol and oxytocin awakening response. Error bars indicate 95% confidence intervals.

Results from separate random intercept and random slopes multilevel analyses showed that on a momentary (within-person) level, presence of affectionate touch was significantly and negatively associated with stress (b = −4.187; t(793) = −2.100; p=0.036), but not with general burden, anxiety, happiness, cortisol, or with oxytocin levels (see Tables 2 and 3, respectively). The negative association with Covid-19-related burden did not reach statistical significance (b = −2.660; t(792) = −1.867; p=0.062).

Table 2
Results of the associations between affectionate touch and self-reported psychological affective states.
‍(A) Random intercept and random slopes models
EffectsGeneral burdenCovid-19 burdenStressAnxietyHappiness
Fixed effectsWithin-person
 Intercept47.734 (4.668); p<0.00143.145 (5.009); p<0.00134.431 (3.596); p<0.00121.096 (3.958); p<0.00171.080 (3.778); p<0.001
 ‍Touch*0.462 (1.561); p=0.767–2.660 (1.424); p=0.062–4.187 (1.994); p=0.036–0.217 (1.510); p=0.8861.557 (1.599); p=0.331
‍Between-person
 Touch*–5.560 (3.791); p=0.144–7.478 (4.068); p=0.067–7.534 (2.907); p=0.010–1.483 (3.210); p=0.64512.420 (3.068); p<0.001
Covariates
 ‍Age–0.186 (0.122); p=0.128–0.158 (0.131); p=0.228–0.136 (0.094); p=0.150–0.090 (0.103); p=0.3880.051 (0.099); p=0.605
 ‍Sex4.709 (3.325); p=0.1583.761 (3.565); p=0.2934.986 (2.524); p=0.0506.568 (2.814); p=0.021–3.318 (2.681); p=0.217
 Day–2.196 (0.864); p=0.011–2.602 (0.981); p=0.008–4.189 (1.228); p<0.001–3.145 (0.875); p<0.0011.637 (0.901); p=0.070
‍Random effects (SD)
 ‍Intercept21.48722.79413.69417.64716.597
 ‍Touch*9.1410.6428.0947.9068.945
 Residual12.82114.86418.74113.07613.478
(B) Random intercept and fixed slopes models
EffectsGeneral burdenCovid-19 burdenStressAnxietyHappiness
Fixed effects
Within-person
 Intercept44.439 (6.130); p<0.00139.748 (6.114); p<0.00132.966 (4.609); p<0.00124.466 (5.277); p<0.00171.437 (4.626); p<0.001
 ‍Touch intensity–0.077 (0.028); p=0.008–0.068 (0.036); p=0.058–0.148 (0.044); p<0.001–0.065 (0.029); p=0.0260.085 (0.030); p=0.005
‍Between-person
 ‍‍Touch intensity–0.121 (0.086); p=0.163–0.138 (0.087); p=0.115–0.223 (0.067); p=0.001–0.102 (0.074); p=0.1710.314 (0.066); p<0.001
Covariates
 Age–0.031 (0.161); p=0.847–0.082 (0.161); p=0.610–0.090 (0.125); p=0.475–0.158 (0.143); p=0.2700.007 (0.122); p=0.952
 ‍Sex2.647 (4.253); p=0.5354.143 (4.214); p=0.3272.650 (3.069); p=3895.932 (3.645); p=0.106–0.030 (3.171); p=0.993
 ‍‍Day–3.788 (1.056); p<0.001–4.930 (1.302); p<0.001–4.791 (1.583); p=0.003–3.718 (1.078); p<0.0013.695 (1.112); p=0.001
Random effects (SD)
 Intercept22.46721.60113.20618.62215.842
 ‍Residual11.81214.62618.30012.12112.591
  1. Table depicts coefficients (standard errors in parentheses) and p-values of associations between (A) the presence and (B) intensity of affectionate touch and psychological variables. Number of observations = 593–1023, Number of participants 162–227.

  2. *

    0 = no, 1 = yes.

  3. 0 = male, 1 = female.

Table 3
Results of the associations between affectionate touch and hormonal levels.
(A) Random intercept and random slopes models(B) Random intercept and fixed slopes models
EffectsCortisolOxytocinEffectsCortisolOxytocin
Fixed effectsWithin-personFixed effectsWithin-person
Intercept2.941 (0.165); p<0.0014.973 (0.405); p<0.001Intercept2.744 (0.370); p<0.0014.657 (0.798); p<0.001
‍Touch*–0.019 (0.060); p=0.756–0.030 (0.076); p=0.688Touch intensity–0.001 (0.001); p=0.3670.006 (0.002); p=0.003
‍Between-personBetween-person
Touch*–0.121 (0.057); p=0.036–0.145 (0.147); p=0.329Touch intensity–0.001 (0.002); p=0.5040.002 (0.003); p=0.489
‍CovariatesCovariates
‍Age–0.001 (0.002); p=0.633–0.013 (0.005); p=0.011Age–0.001 (0.003); p=0.776–0.019 (0.007); p=0.007
‍Sex–0.020 (0.044); p=0.647–0.129 (0.123); p=0.293Sex–0.001 (0.063); p=0.989–0.167 (0.151); p=0.272
‍Day–0.050 (0.035); p=0.155–0.010 (0.059); p=0.863Day–0.032 (0.051); p=0.5280.077 (0.086); p=0.370
‍Time-fall –0.476 (0.023); p<0.0010.051 (0.035); p=0.145Time-fall –0.471 (0.032); p<0.0010.025 (0.049); p=0.611
‍Body mass index–0.014 (0.005); p=0.0090.012 (0.015); p=0.422Body mass index0.005 (0.009); p=0.5850.014 (0.022); p=0.528
‍Eating*–0.077 (0.075); p=0.3020.032 (0.117); p=0.788Eating*–0.125 (0.099); p=0.2100.099 (0.149); p=0.506
‍Drinking*–0.007 (0.081); p=0.9330.045 (0.126); p=0.722‍Drinking*–0.011 (0.108); p=0.919–0.039 (0.161); p=0.810
‍Caffeine*0.127 (0.043); p=0.004–0.099 (0.070); p=0.158‍Caffeine*0.131 (0.064); p=0.0430.157 (0.102); p=0.125
‍Alcohol*–0.030 (0.063); p=0.628–0.174 (0.096); p=0.070‍Alcohol*–0.012 (0.072); p=0.865–0.154 (0.112); p=0.171
‍Cigarettes*0.104 (0.072); p=0.1530.053 (0.139); p=0.705Cigarettes*0.089 (0.098); p=0.368–0.054 (0.179); p=0.766
‍Physical activity*0.039 (0.040); p=0.3320.108 (0.065); p=0.098‍Physical activity*–0.072 (0.056); p=0.1990.014 (0.086); p=0.873
‍Sleep duration §–0.003 (0.004); p=0.371–0.009 (0.008); p=0.272‍Sleep duration §–0.029 (0.029); p=0.3210.046 (0.060); p=0.438
‍Sleep quality –0.001 (0.001); p=0.289–0.000 (0.002); p=0.912‍Sleep quality –0.001 (0.001); p=0.477–0.000 (0.002); p=0.960
Problem falling asleep*0.031 (0.053); p=0.552–0.245 (0.101); p=0.015Problem falling asleep*0.030 (0.081); p=0.709–0.133 (0.148); p=0.370
Sleeping pills*0.032 (0.109); p=0.7720.031 (0.254); p=0.904Sleeping pills*–0.203 (0.290); p=0.4870.771 (0.453); p=0.093
Forced awake*–0.008 (0.041); p=0.8380.012 (0.088); p=0.894Forced awake*0.016 (0.058); p=0.7880.075 (0.117); p=0.520
Brushing teeth*0.045 (0.036); p=0.2170.096 (0.056); p=0.090Brushing teeth*0.029 (0.051); p=0.5680.016 (0.077); p=0.839
Random effects (SD)Random effects (SD)
‍Intercept0.1490.542Intercept0.1790.546
‍Touch*0.2730.166Touch intensity
‍Residual0.3340.477Residual0.3390.461
  1. Table depicts unstandardized coefficients (standard errors in parentheses) and p-values of hormonal associations with (A) the presence and (B) intensity of affectionate touch. Number of observations = 251–545. Number of participants = 88–152.

  2. *

    0 = no, 1 = yes.

  3. 0 = male, 1 = female.

  4. 0 = time point 1–3, 1 = time point 4, 2 = time point 5, 3 = time point 6.

  5. §

    In hours.

  6. 1 = very bad, 101 = very good.

On a between-person level, affectionate touch was significantly associated with lower cortisol (b = −0.121; t(128) = −2.118; p=0.036) (see Table 3), stress (b = −7.534; t(223) = −2.592; p=0.010), as well as with higher happiness (b = 12.420; t(223) = 4.049; p<0.001) levels (see Table 2), but not with general burden or anxiety. The negative association with Covid-19-related burden did not reach statistical significance (b = −7.478; t(223) = −1.838; p=0.067) (see Figure 3).

Associations between occurrence of affectionate touch and psychological and hormonal state.

Panels (A) to (G) illustrate violin plots with density distributions of subjective ratings of general and Covid-19-related burden, stress, anxiety, happiness, cortisol, and oxytocin, depending on whether touch occurred or not. Each dot represents one assessment. Central dots (black) represent each mean. Black lines represent the standard deviations. * indicates statistically significant results (p<0.05). +indicates a statistical trend (p<0.1). ♠ indicates statistically significant between-person effect. ♣ indicates statistically significant within-person effect.

We also analyzed the intensity of experienced affectionate touch as a predictor for psychological and hormonal outcomes, separating within-person and between-person effects. We found that within a person there were significant and negative associations of the intensity of momentary affectionate touch with anxiety (b = −0.065; t(430) = −2.232; p=0.026), stress (b = −0.148; t(430) = −3.363; p<0.001), general burden (b = −0.077; t(430) = –2.687; p=0.008), and positive associations with momentary happiness (b = 0.085; t(430) = 2.795; p=0.005). The negative association with Covid-19-related burden, however, did not reach statistical significance (b = −0.068; t(428) = −1.900; p=0.058) (see Table 2).

Momentary oxytocin levels were significantly higher with more intensive affectionate touch (b = 0.006; t(149) = 3.058; p=0.002) and cortisol levels were descriptively slightly lower; however, this effect was not statistically significant (see Table 3).

Furthermore, on the between-person level, higher intensity of affectionate touch was significantly associated with less stress (b = −0.223; t(158) = −3.318; p=0.001) and greater happiness (b = 0.314; t(159) = 4.764; p<0.001) (see Figure 4), but not with anxiety, general burden, Covid-19-related burden or hormonal levels. No statistically significant sex effects emerged in any of these analyses.

Associations between momentary affectionate touch intensity and psychological and hormonal state.

Panels (A) to (G) illustrate the results of random intercept and fixed slopes models depicting associations of momentary intensity of touch with self-reported general burden and Covid-19 related burden, stress, anxiety, happiness, cortisol, and oxytocin. Gray lines indicate the overall predicted slope, whereas the blue lines indicate the individual’s predicted slopes with their minimum and maximum predicted values as endpoints. Gray areas depict the 95% confidence band. * indicates statistically significant results (p<0.05).+indicates a statistical trend (p<0.1). ♠ indicates statistically significant between-person effect. ♣ indicates statistically significant within-person effect.

In a final set of analyses, we conducted multilevel structural equation models to investigate whether there was evidence for oxytocin mediating the effects of affectionate touch on cortisol and/or self-report outcomes. None of the indirect effects of the presence of affectionate touch (p>0.773) or intensity of affectionate touch (p>0.194) on the within-person level were statistically significant. Furthermore, on the within-person level, the correlations of momentary oxytocin with cortisol (r = −0.016, p=0.521), Covid-19-related burden (r = −0.030, p=0.210), stress (r = −0.020, p=0.442), anxiety (r = −0.013, p=0.660), and happiness (r = 0.031, p=0.328) were not statistically significant.

Discussion

This study investigated the associations of affectionate touch with self-reported mental health and mood, as well as with momentary endogenous oxytocin and cortisol levels during the first Covid-19 lockdown in the spring of 2020.

In our online survey data, we found significant main effects of sex, psychopathology, and loneliness on psychological distress (HADS total score) and, more specifically, on anxiety (HADS Anxiety) and depressive (HADS Depression) symptoms. Individuals reported higher levels of depression and anxiety, especially if they were female or burdened by mental illness or loneliness. These data are in line with previous studies (Fancourt et al., 2021; Mata et al., 2021; Pierce et al., 2020). Interestingly, analyses showed that the attitude toward touch significantly moderated the association between loneliness and the HADS total score, as well as the HADS Anxiety subscale. Thus, individuals with a positive attitude and affect toward social touch experiencing loneliness showed higher distress and anxiety in times of Covid-19-related lockdown. Whether these moderation effects are apparent outside of pandemic-caused physical restrictions is unknown and should be addressed in future studies. These findings support our hypothesis that touch deprivation and loneliness could be related to anxiety symptoms. However, it is further important to note that about 20% of the participants reported having at least one psychiatric diagnosis. In comparison, the pre-pandemic 12 mo prevalence in the general population in Germany is about 28% (Jacobi et al., 2014). Thus, our sample seems to be slightly less burdened compared to the general population, which partly limits the generalizability of the results.

Results of the psychobiological EMA study in a large sample of a broad age range showed that the presence of affectionate touch was negatively associated with stress and cortisol levels and positively linked with happiness. Moreover, the more intensely affectionate touch was experienced, the lower were subsequent subjective anxiety, general burden and stress levels. Higher intensity of affectionate touch was associated with elevated oxytocin and self-reported happiness. Note that in this data assessment in everyday life affectionate touch could not be experimentally manipulated and these correlative results should be interpreted with caution. However, the results could be interpreted that affectionate touch during the Covid-19 pandemic buffers anxiety and stress and downregulates the HPA response, particularly cortisol. At the same time, the intensity of affectionate touch was associated with increased endogenous oxytocin levels and subjective happiness. Interestingly, a recent study demonstrated that foot massage was rated as more pleasurable and rewarding and was associated with a higher increase of oxytocin after the massage administered by hand as compared to machine-administered massage, although the intensity of the massage was rated similarly (Li et al., 2019). Furthermore, the oxytocin system and its potential stress-ameliorating effects seem to be triggered by meaningful and intense touch in particular (Eckstein et al., 2020).

To our knowledge, the present data provide the first empirical evidence to suggest that affectionate touch is related to reduced anxiety, stress, and general burden, as well as stress-responsive cortisol levels and at the same time is linked to higher endogenous oxytocin levels and well-being in an ecologically valid everyday life setting.

Up to now, no systematic data on daily oxytocin profiles and momentary oxytocin levels in a large sample of men and women from varying age groups have been available. Based on single peripheral oxytocin measures from relatively small samples so far (see Valstad et al., 2017, for an overview), it was the object of debate whether peripheral oxytocin levels might be interpretive of emotional functioning or correspond with stressful experiences (Engel et al., 2019). Moreover, in the last decade, there has been an extensive discussion about the reliability and validity of peripheral oxytocin measures (Martins et al., 2020; Szeto et al., 2011; Tabak et al., 2023). There are several methodological issues and challenges associated with measuring oxytocin in blood plasma and saliva samples (Tabak et al., 2023). For example, studies have shown that oxytocin concentrations after sample extraction are much lower compared to unextracted oxytocin measurements (Szeto et al., 2011). Additionally, the correlations between extracted and unextracted oxytocin levels as well as between saliva and plasma oxytocin concentrations have been inconsistent across studies (e.g. Hoffman et al., 2012; Martins et al., 2020; Nagahashi-Araki et al., 2022; Szeto et al., 2011). These inconsistencies might be due to numerous reasons including differences in the study populations, methods of sample processing, and analyses. In particular, using different assay types (e.g. radioimmunoassay vs. enzyme immunoassay), as well as sample preparation (extraction vs. non-extraction), may contribute to these inconsistencies and make it difficult to compare results between studies (Tabak et al., 2023). Since unextracted samples were used in this study, the concentrations probably represent both free and bound oxytocin (MacLean et al., 2019), thereby potentially limiting the comparability with studies using extracted samples.

Another important issue is the intraindividual stability of oxytocin over time (Feldman et al., 2013; Martins et al., 2020; Schneiderman et al., 2012). A recent study reports no correlation of single oxytocin measures between several assessments, indicating that single measures of oxytocin might not be reliable to represent oxytocin baseline levels (Martins et al., 2020). In our sample, we found a significant positive correlation of mean oxytocin values between the two assessment days. As the fluctuations of oxytocin throughout the day were apparent in our study, correlating mean values of six oxytocin measures over the day might represent the individual baseline oxytocin levels better than single measures (this issue has been also discussed in Tabak et al., 2023). However, the difference between our data and the findings reported by Martins et al., 2020 might also reflect methodological differences as previous studies with unextracted samples also reported positive correlations over time (Feldman et al., 2013; Schneiderman et al., 2012).

In our study, we found a positive association between momentary oxytocin levels and the intensity of affectionate touch on a within-person level. Our findings are in line with previous studies showing an increase in salivary oxytocin after self-touch (de Jong et al., 2015), standardized touch (Portnova et al., 2020), and massage (Li et al., 2019; Moussa et al., 2021). Notably, in recent animal research, it has been shown that the density of oxytocin neurons in the brains of male mice increased after social isolation and that oxytocin neurons are involved in regulating social craving (Musardo et al., 2022). Moreover, social touch has been associated with central nervous system oxytocin activation and secretion into the periphery (Tang et al., 2020). Our data are in line with this research; however, we did not find support for the hypothesis that peripheral oxytocin directly mediated the effects of affective touch on momentary subjective distress and well-being on a statistical level. Central nervous oxytocin dynamics (receptor sensitivity, real-time levels, local gene expression, or methylation) and their interaction with the HPA axis cannot be measured in the living human brain presently (Quintana et al., 2019). Thus, although speculative, our results might suggest that central nervous system oxytocin mechanisms as triggered by touch can modulate endocrine outcomes (peripheral oxytocin, cortisol) and subjective distress, as well as well-being. Furthermore, affectionate touch might influence different outcome levels in parallel (not mediated) processes, or the effects of oxytocin on perceived distress, well-being, and cortisol could unfold across a longer time span, instead of on a moment-to-moment basis. Alternatively, it is also possible that initially higher levels of participants’ well-being (including lower stress, anxiety, or burden) might have increased feelings of closeness and therefore promoted affectionate touch. Also, other aspects that accompany physical contact such as eye contact, compliments, or affective closeness with loved ones could have contributed to the beneficial effects.

The study has some limitations that need to be addressed. The assessment of individual depression and anxiety levels was based on the self-reports using the HADS. Although this instrument has been validated and repeatedly used in clinical practice (Hinz and Brähler, 2011), it does not replace clinical interviews and might be influenced by self-report bias. Furthermore, in the EMA measures we used single items to minimize the drop-out rate during the study, but these items might not comprehensively reflect the individual’s experience of burden, stress, anxiety, or happiness. In addition, touch from strangers was restricted during the pandemic; thus, affectionate touch experiences were probably mostly from family contacts. Therefore, we cannot draw differential conclusions about varying contexts of touch. While in general, being touched by strangers may be rated as less pleasant, during times of a pandemic it is also associated with a higher risk for infection. In contrast, during the lockdown, touch at home may be experienced either as harmless and pleasant (Sorokowska et al., 2021) or as too close to feel comfortable during times of limited distraction and constant and close physical contact with family members. The latter is particularly relevant when it comes to the association of touch (yes/no) with cortisol and oxytocin. While touch, per se, seemed to be associated with reduced cortisol levels in the present sample, oxytocin secretion appears to be related to the intensity of touch. However, since this is a cross-sectional study, we here interpret associations rather than causal effects. The Covid-19-related lockdown provided us with a social situation to study the effects of touch between family/household members (romantic couples, parent–child dyads, etc.) in a relatively controlled setting. While the situation was quite specific and limited the generalizability of the results to everyday life in pre- or post-pandemic conditions, the risk of viral infection was not the only concern among the population. Participants reported significant concerns about being isolated from others and how long it might take for them to get back to normal (Hopf et al., 2021). These concerns related, at least in part, to the fear of loneliness, defined as a perceived lack of social connection and the distress this causes (Bekhet et al., 2008). Thus, these results obtained in the general population during pandemic-related restrictions can be partially generalized to other situations, such as a lack of social contacts due to migration, physical illnesses/disabilities, or other reasons.

Conclusion

Taken together, the present findings provide support for the links of affectionate touch with more positive mental health outcomes during times of prolonged stress. Notably, the above associations with lower anxiety, better mood, and reduced cortisol levels in everyday life during the Covid-19 lockdown showed that more intense affectionate touch is related to higher salivary oxytocin levels on a moment-to-moment basis. This suggests that endogenous oxytocin might be stimulated through targeted behavior (e.g. social touch), which could have implications for prevention and interventions for individuals who are particularly vulnerable during times of stress and social isolation.

Data availability

The datasets analyzed and presented in this manuscript are openly available online (https://doi.org/10.11588/data/WFNWJT).

The following data sets were generated
    1. Schneider E
    2. Hopf D
    3. Ditzen B
    (2023) heiDATA
    Affectionate touch and diurnal oxytocin levels: An ecological momentary assessment study [Research Data].
    https://doi.org/10.11588/data/WFNWJT

References

    1. Döring N
    2. Bortz J
    (1993)
    Psychometrische Einsamkeitsforschung: Deutsche Neukonstruktion der UCLA loneliness scale
    Diagnostica 39:224–239.
    1. Kirschner M
    2. Kirschner R
    (2019)
    Hand massage reduces perceived stress, anxiety and fatigue
    International Journal of Innovative Studies in Medical Sciences 3:1–6.
    1. Morhenn V
    2. Beavin LE
    3. Zak PJ
    (2012)
    Massage increases oxytocin and reduces Adrenocorticotropin hormone in humans
    Alternative Therapies in Health and Medicine 18:11–18.

Decision letter

  1. Matthias Gamer
    Reviewing Editor; University of Würzburg, Germany
  2. Christian Büchel
    Senior Editor; University Medical Center Hamburg-Eppendorf, Germany
  3. Matthias Gamer
    Reviewer; University of Würzburg, Germany
  4. Elizabeth Anne Broadbent
    Reviewer
  5. Keith M Kendrick
    Reviewer; University of Electronic Science and Technology of China, China

In the interests of transparency, eLife publishes the most substantive revision requests and the accompanying author responses.

Decision letter after peer review:

Thank you for submitting your article "Affectionate touch and diurnal oxytocin levels: An ecological momentary assessment study" for consideration by eLife. Your article has been reviewed by 3 peer reviewers, and the evaluation has been overseen by a Reviewing Editor and Christian Büchel as the Senior Editor. The following individuals involved in review of your submission have agreed to reveal their identity: Matthias Gamer (Reviewer #1); Elizabeth Anne Broadbent (Reviewer #2); Keith M Kendrick (Reviewer #3).

The reviewers have discussed their reviews with one another, and the Reviewing Editor has drafted this to help you prepare a revised submission.

Essential revisions:

1) The interpretation and conclusions of this study should be toned down as this purely correlative approach does not allow for drawing causal inferences about the relationship between affectionate touch, anxiety, stress, and hormone levels.

2) It should be acknowledged that different types of affectionate touch and contexts might differentially affect wellbeing and hormone levels. Furthermore, other social and physical factors might have been affected by the restrictions during the first wave of the COVID-19 pandemic in Germany. Such potential confounds should be discussed in more detail with respect to the generalizability of the current findings to other contexts.

3) Given ongoing discussions in the field how saliva oxytocin concentrations should be measured and interpreted, some more details on the current procedure should be added and discussed with respect to reported concerns in the literature.

4) The statistical analyses should be better justified and revised in accordance with the suggestions of Reviewer #1.

Reviewer #1 (Recommendations for the authors):

Schneider et al. investigated the association of affectionate touch and subjective well-being while also considering salivary cortisol and oxytocin as objective measures of stress responsiveness or buffering, respectively. The study was carried out during the first wave of the COVID-19 pandemic in Germany that was characterized by numerous restrictions of the social life to minimize the risk of infections. Following a cross-sectional online survey with more than 1,000 participants, nearly 250 subjects provided ecologically momentary assessment (EMA) data over two consecutive days including repeated measurements of salivary cortisol and oxytocin. Multilevel models revealed that the occurrence of affectionate touch was related to lower stress and increased happiness as well as lower cortisol levels across individuals. Moreover, momentary intensity of affectionate touch within individuals correlated negatively with general burden, stress, and anxiety and positively with happiness and oxytocin levels. There was no evidence for oxytocin mediating the effects of affectionate touch on objective (cortisol levels) or subjective stress responses (self-report ratings).

The current article has several strengths but also some weaknesses that are outlined below.

Strengths:

1. The current study including detailed EMA data as well as neuroendocrine variables nicely extends previous cross-sectional examinations that solely relied on self-report measures. This approach using large sample sizes allows for a comprehensive analysis of how affectionate touch is related to subjective well-being as well as salivary cortisol and oxytocin levels. Thereby, it seems possible not only to identify correlations between psychological constructs but also to address potential effects on a biological level.

2. The analyses are mostly adequate (see comments below) and consider a large amount of potential confounds.

3. The hypotheses, methods, and results are presented clearly and comprehensively, and I enjoyed reading the article.

Weaknesses:

1. Due to the correlational approach of the current study, explicit conclusions about the causality of effects are problematic. While it seems plausible that "affectionate touch might act as a buffer against anxiety and stress and might downregulate HPA axis outcomes" (line 305), there are also numerous alternative explanations for such relationship. For example, low levels of stress (e.g., caused by some unknown variable) might enhance social closeness and thereby modulate the occurrence and experience of affectionate touch as well as other psychological factors (e.g., anxiety or happiness). Alternatively, the mere presence of a loved one could be responsible for all the observed effects while at the same time also increasing the likelihood of affectionate touch. To consider such explanations, the discussion should be toned down.

2. The current study was conducted in a very specific context – the first wave of a pandemic that was associated with significant societal restrictions and uncertainties. While I fully agree that the current study is very relevant to understanding how affectionate touch might mitigate the deleterious mental health effects of such a situation, I was wondering to what extent the current findings could be generalized to other contexts. I would like to see at least a brief discussion of this issue.

3. While I appreciate that the authors carefully controlled for numerous potential confounds, I did not understand all relevant analysis steps. For example, why are random intercept and random intercept + random slopes models presented for all analyses? To my opinion, it would make more sense to only report the model that fits the data better (taking model complexity into account). Relatedly, the calculation of the time-rise (change from time point 1 to 3) and time-fall covariates (change from time point 3 to 6) seemed rather arbitrary to me. Figure 2 seems to indicate that time point 2 seems more appropriate than time point 3.

– In the abstract, it should be clarified whether the results describe between-person or within-person effects.

– In the section on "sample characteristics", it is mentioned that a rather large proportion of the sample indicated to suffer from a diagnosed mental disorder (>20%). I was wondering how representative such sample is?

– Table 1: The n's of the online survey participants (< 300) seem very low given that 1,050 participants provided data.

– Line 186: Covid-19 burden also failed to show a significant effect.

– Figure 3: Please use consistent terminology – "general burden" instead of "self-reported burden"

– Line 231: "general burden" instead of only "burden".

– Table 4 and 5: No model is marked with an asterisk, but the meaning of this symbol is still mentioned in the notes.

– Lines 312-315: Please double-check the results description. Significant effects were only obtained for general but not for Covid-19 related burden.

– Line 375: Typo "to report n conflicts of interest"?

– Line 391: What is meant by the sentence "To reduce potential noise in the data sampling was constantly monitored by study members."

– Line 432: The version of R is more relevant than the version of R Studio but I suppose that 4.1.1 already refers to the version of R.

Reviewer #2 (Recommendations for the authors):

The authors were seeking to investigate whether affectionate touch was associated with mental health during the pandemic using self-report measures and hormone levels.

A major strength of the paper was the large sample size. A weakness is that the study is observational so causation cannot be inferred. The authors achieved their aims and results support their conclusions. The paper is important as it shows the association between affectionate touch and lower distress during social isolation.

Reviewer #3 (Recommendations for the authors):

The research described in this paper mainly used online surveys which attempted to establish whether attitudes to and receipt of affectionate touch from others might act as a buffer to reduce stress in challenging circumstances with reduced social contact such as during the COVID-19 lockdowns.

The main strength of the study is in its large scale and systematic and rigorous analysis of the measures collected online and also the detailed analysis of behavioral and endocrine measures in a smaller cohort of subjects at multiple time points over two days. The main weakness resides in insufficient measures being included to clearly establish the importance of different types of affectionate touch, and their context, and whether other social (i.e. social interactions not involving touch) or physical factors (such as activity/exercise) likely to be influenced during the lockdown period showed similar associations with stress, anxiety and endocrine measures.

Overall, the findings reported by the authors suggest firstly that individuals with a positive attitude to social touch are more likely to exhibit stress and loneliness problems in circumstances where social contact is reduced, such as during the COVID-19 lockdowns. Secondly, findings provide some initial support for the potential importance of receiving affectionate touch in helping to reduce stress and anxiety and influence concentrations of hormones associated with their regulation (i.e. cortisol and oxytocin). However, further more detailed and controlled experiments will be needed to establish fully the specificity of these potentially beneficial effects of affectionate touch, the most important aspects and contexts of affectionate touch which are potentially influencing wellbeing and the mechanisms involved.

Touch, particularly in the form of massage, has well established anxiolytic and other effects on well-being, although the current study has extended this to demonstrate that the amount and intensity of interpersonal affectionate received in everyday contexts may also be important. Giving and receiving warm hugs and cuddles may be more important to our general mental wellbeing than we realize.

In terms of the measures recorded while the focus is on affectionate touch do the authors also have any data on other social (face to face conversations and other mutual activities not involving intimate touch) and physical activity likely to have some effect on cortisol and oxytocin concentrations?

While online information collection clearly has limitations concerning the amount of different variables which can be included I feel that given the central importance of the question concerning the role of receipt of affective touch many relevant details appear not to have been either included or considered. For example, were all subjects in a relationship and co-habiting/sleeping together? Was only one individual per household recruited? Given that a lot of affective touch can also be received from children and even pets, was this also considered? If subjects included a range of individuals both living alone and cohabiting I would have anticipated this being included as a variable, especially given the focus on COVID where individuals living alone might be expected to experience the greatest reduction in affective touch relative to those co-habiting/living in a family environment. As far as I can see there is no indication as to whether the subjects included in the study considered that they did receive less affectionate touch during the sampling period relative to before it?

The peak of oxytocin concentrations and reduced cortisol ones specifically at awakening is interesting but I wonder whether this simply reflects couples experiencing morning cuddles in a warm bed etc or a specific effect of awakening?

The reported saliva oxytocin concentrations are high and this is an ongoing issue in the field. I assume that the samples were run unextracted and also not in duplicate? While there is a clear recommendation now from both researchers and assay manufacturers to use extraction for plasma samples there is less consensus concerning saliva ones. With the widely used ELISA ENZO assay used in the current paper findings clearly shown that in extracted saliva samples generally result in baselines of <10pg/ml (see for example Fujii et al. 2016 Sci Reps 6:38662; Le et al. 2022 Psychother and Psychosom 91(5):335-347), while for unextracted ones this figure is much higher (up to 200 pg/ml), as in the current study. With plasma samples correlations between extracted and unextracted samples are weak, however a recent study comparing baseline and stimulated concentrations in saliva concentrations in post-partum women did report a significant low correlation (0.32-0.38) although profiles of release were different and concentrations were 10-fold higher in unextracted samples (Nagahasi-Araki et al. 2022, BMC Pregnancy and Childbirth 22:711). It is likely that higher concentrations in unextracted samples are contributed to by both the assay type (RIA vs ELISA) and whether the antibody used is sensitive to oxytocin fragments as well as oxytocin itself. Using a gold standard RIA, for example, as opposed to an ELISA, papers have consistently reported similar low oxytocin concentrations in both extracted plasma samples and unextracted saliva samples (generally below 5pg/ml – de Jong et al. 2015 Psychoneuroendocrinology 62:381-388; Martins et al., 2020 eLife 9:e62456), whereas with the ENZO ELISA extracted plasma and extracted saliva samples yield similar concentrations (Le et al., 2022). Indeed, there is a clear expectation that oxytocin concentrations should be similar in plasma and saliva since there is not reason to expect radically higher concentrations in saliva. The authors really need to at least acknowledge this issue and include it as a potential limitation. A recent review has, for example, highlighted the need for papers to make it very clear whether extracted or unextracted values are being reported (Tabak et al. 2023. Mol Psychiatry 28:127-140).

The authors do not mention a relevant study on humans reporting that foot massage both increases plasma oxytocin concentrations and corresponding neural responses in the orbitofrontal cortex and superior temporal sulcus but not somatosensory cortex. The study also reports that greater basal and manual touch-evoked oxytocin concentrations were associated with more positive attitudes to social touch (also using the STQ).(Li et al. 2019 Psychoneuroendocrinology 101:193-203).

While the authors have mainly decided to focus on presenting their findings and avoided entering into mechanistic explanations of how increased saliva oxytocin concentrations might influence levels of stress and anxiety etc I think it would help if there was a little more explanation of how saliva changes might influence the brain and endocrine system. In particular, the authors may wish to discuss the conclusion in another eLife paper that salivary and plasma oxytocin concentrations are not reliable trait markers for the oxytocin system in humans (Martins et al. 2020, eLife 9:62456).

In the discussion the authors state that affectionate touch reduced the COVID-19 related burden (line 313) but as far as I can see there was no significant effect found for this?

A small point but in the discussion of the study on rats by Tang et al. on line 326 animals are described as receiving "affectionate touch" from each other whereas perhaps simply receiving "social touch" might be more appropriate in this case.

https://doi.org/10.7554/eLife.81241.sa1

Author response

Essential revisions:

1) The interpretation and conclusions of this study should be toned down as this purely correlative approach does not allow for drawing causal inferences about the relationship between affectionate touch, anxiety, stress, and hormone levels.

We thank the editor for this comment, and we agree that our study design does not allow causal conclusions. We, therefore, toned down the interpretation of our results and provided alternative possible explanations for the correlational results we find (p. 20-22).

2) It should be acknowledged that different types of affectionate touch and contexts might differentially affect wellbeing and hormone levels. Furthermore, other social and physical factors might have been affected by the restrictions during the first wave of the COVID-19 pandemic in Germany. Such potential confounds should be discussed in more detail with respect to the generalizability of the current findings to other contexts.

We agree with the editor and now more precisely discuss the potential effects of different kinds of touch on health-related outcomes and the generalizability of our results to non-pandemic contexts (p. 18- 22)

3) Given ongoing discussions in the field how saliva oxytocin concentrations should be measured and interpreted, some more details on the current procedure should be added and discussed with respect to reported concerns in the literature.

Measuring oxytocin concentrations in saliva and blood samples is indeed a highly discussed topic. We now address this topic in our manuscript, taking into account the most recently published review articles and discussing the descriptive rather than mechanistic findings we report here (p. 19-20)

4) The statistical analyses should be better justified and revised in accordance with the suggestions of Reviewer #1.

As suggested, we have thoroughly revised and adjusted the statistical approach where necessary (p. 9).

Reviewer #1 (Recommendations for the authors):

Schneider et al. investigated the association of affectionate touch and subjective well-being while also considering salivary cortisol and oxytocin as objective measures of stress responsiveness or buffering, respectively. The study was carried out during the first wave of the COVID-19 pandemic in Germany that was characterized by numerous restrictions of the social life to minimize the risk of infections. Following a cross-sectional online survey with more than 1,000 participants, nearly 250 subjects provided ecologically momentary assessment (EMA) data over two consecutive days including repeated measurements of salivary cortisol and oxytocin. Multilevel models revealed that the occurrence of affectionate touch was related to lower stress and increased happiness as well as lower cortisol levels across individuals. Moreover, momentary intensity of affectionate touch within individuals correlated negatively with general burden, stress, and anxiety and positively with happiness and oxytocin levels. There was no evidence for oxytocin mediating the effects of affectionate touch on objective (cortisol levels) or subjective stress responses (self-report ratings).

The current article has several strengths but also some weaknesses that are outlined below.

Strengths:

1. The current study including detailed EMA data as well as neuroendocrine variables nicely extends previous cross-sectional examinations that solely relied on self-report measures. This approach using large sample sizes allows for a comprehensive analysis of how affectionate touch is related to subjective well-being as well as salivary cortisol and oxytocin levels. Thereby, it seems possible not only to identify correlations between psychological constructs but also to address potential effects on a biological level.

2. The analyses are mostly adequate (see comments below) and consider a large amount of potential confounds.

3. The hypotheses, methods, and results are presented clearly and comprehensively, and I enjoyed reading the article.

Weaknesses:

1. Due to the correlational approach of the current study, explicit conclusions about the causality of effects are problematic. While it seems plausible that "affectionate touch might act as a buffer against anxiety and stress and might downregulate HPA axis outcomes" (line 305), there are also numerous alternative explanations for such relationship. For example, low levels of stress (e.g., caused by some unknown variable) might enhance social closeness and thereby modulate the occurrence and experience of affectionate touch as well as other psychological factors (e.g., anxiety or happiness). Alternatively, the mere presence of a loved one could be responsible for all the observed effects while at the same time also increasing the likelihood of affectionate touch. To consider such explanations, the discussion should be toned down.

We want to thank the reviewer for this constructive comment. We agree that with this nonexperimental design causal interpretations are critical, and the data must be interpreted cautiously. Since in this data assessment in everyday life, affectionate touch could not be experimentally manipulated and the results are correlational, we agree that causal conclusions are difficult to draw. We have therefore rephrased our conclusions (p. 18):

“Note that in this data assessment in everyday life, affectionate touch could not be experimentally manipulated and these correlative results should be interpreted with caution…”

As suggested, we have included alternative explanations in the Discussion section (p. 21):

” Alternatively, it is also possible that initially higher levels of participants’ well-being (including lower stress, anxiety, or burden) might have increased feelings of closeness and therefore promoted affectionate touch. Also, other aspects that accompany physical contact, such as eye contact, compliments, or affective closeness with loved ones could have contributed to the beneficial effects.”

2. The current study was conducted in a very specific context – the first wave of a pandemic that was associated with significant societal restrictions and uncertainties. While I fully agree that the current study is very relevant to understanding how affectionate touch might mitigate the deleterious mental health effects of such a situation, I was wondering to what extent the current findings could be generalized to other contexts. I would like to see at least a brief discussion of this issue.

We thank the reviewer for bringing up this very important aspect of our data. On the one hand, the pandemic lockdown strengthens the impact of our results due to the clinically relevant and ecologically valid stress situation. On the other hand, it is (luckily and hopefully) a very uncommon and rare occurrence, which clearly raises the question of how to interpret the data in light of other situations. We now address this question in our manuscript (p. 22):

“While the situation was quite specific and limited the generalizability of the results to everyday life in pre- or post-pandemic conditions, the risk of viral infection was not the only concern among the population. Participants reported significant concerns about being isolated from others and how long it might take for them to get back to normal (Hopf et al., 2021). These concerns are related, at least in part, to the fear of loneliness, defined as a perceived lack of social connection and the distress this causes (Bekhet et al., 2008). Thus, these results obtained in the general population during pandemic-related restrictions can be partially transferred to other situations, such as a lack of social contacts due to migration, physical illnesses/disabilities, or other reasons.”

3. While I appreciate that the authors carefully controlled for numerous potential confounds, I did not understand all relevant analysis steps. For example, why are random intercept and random intercept + random slopes models presented for all analyses? To my opinion, it would make more sense to only report the model that fits the data better (taking model complexity into account). Relatedly, the calculation of the time-rise (change from time point 1 to 3) and time-fall covariates (change from time point 3 to 6) seemed rather arbitrary to me. Figure 2 seems to indicate that time point 2 seems more appropriate than time point 3.

We appreciate this reviewer’s comment. For the sake of transparency, we initially reported all random intercept and fixed slopes as well as random intercepts and random slopes models in the previous version of this manuscript. However, we agree with the reviewer that it makes sense to optimize and reduce the number of reported models. To keep the report consistent, we decided to report random intercept and random slopes models on affectionate touch (yes/no), since four of the seven models showed a statistically better fit when the slopes were set at random. For models with the intensity of touch as the outcome variable, we now report random intercepts and fixed slopes models (none of these models showed a better model fit when slopes were set at random). The following information is now in the statistical analyses section (p. 9):

“For the models on affectionate touch as a binary variable (yes/no), we report random intercept and random slopes models in the Results section, since these showed a statistically better fit as compared to random intercept and fixed slopes models. However, for the dimensional intensity of affectionate touch, we report random intercept and fixed slopes models, since the random slope models did not yield a better model fit.”

Regarding the time-rise (change from time point 1 to 3) and time-fall (change from time point 3 to 6), we initially made this decision based on the theoretical assumption that the first three time points (all scheduled within 45 minutes after awakening) represent the typical cortisol awakening response (CAR), which represents a distinct process from the diurnal change throughout the rest of the day (we relied on tutorials on how to calculate and interpret CAR (Stalder et al., 2016)). We had therefore made the a-priori decision for a separate representation of the cortisol awakening response in our statistical models. However, after re-evaluating our decision in response to this comment, we realized that the first two time points are irrelevant to our research question because affectionate touch was not assessed at the first two time points of a day (only at time points 3-6). Consequently, there is no need to control for the initial rise in cortisol from time point 1 to 3 in our models. We have therefore adjusted our models and removed the covariate representing time-rise from our multilevel models. These models now contain only a time variable representing the linear decrease in cortisol levels throughout the day after the cortisol awakening response.

– In the abstract, it should be clarified whether the results describe between-person or within-person effects.

We thank the reviewer for this comment and have included information on within-person and between-person effects in the abstract.

– In the section on "sample characteristics", it is mentioned that a rather large proportion of the sample indicated to suffer from a diagnosed mental disorder (>20%). I was wondering how representative such sample is?

We agree that these are high numbers of self-reported mental impairment, indeed. However, the pre-pandemic 12-month prevalence of any psychiatric disease has been estimated to be about 28% (Jacobi et al., 2014). Thus, the observed prevalence rate of about 20% is, actually, even lower than the rate reported for the general population in Germany. As we believe that this is an important point, we now comment on this data in the Discussion section (p.17):

"However, it is further important to note that about 20% of the participants reported having at least one psychiatric diagnosis. In comparison, the pre-pandemic 12-month prevalence in the general population in Germany is about 28 % (Jacobi et al., 2014). Thus, our sample seems to be slightly less burdened compared to the general population, which partly limits the generalizability of the results."

– Table 1: The n's of the online survey participants (< 300) seem very low given that 1,050 participants provided data.

The reviewer might have misunderstood the first column as the total sample: n = 227, which refers to the number of male participants, whereas n = 815 refers to the number of female participants. We now have changed the presentation of the numbers in table 1.

– Line 186: Covid-19 burden also failed to show a significant effect.

We apologize for this mistake and changed it in the text.

– Figure 3: Please use consistent terminology – "general burden" instead of "self-reported burden"

Thanks for pointing out the inconsistent wording, we have changed it in the figure.

– Line 231: "general burden" instead of only "burden".

Done.

– Table 4 and 5: No model is marked with an asterisk, but the meaning of this symbol is still mentioned in the notes.

We have removed the asterisk from the table notes.

– Lines 312-315: Please double-check the results description. Significant effects were only obtained for general but not for Covid-19 related burden.

Thanks for pointing this out, we have corrected it in the text.

– Line 375: Typo "to report n conflicts of interest"?

We corrected the typo.

– Line 391: What is meant by the sentence "To reduce potential noise in the data sampling was constantly monitored by study members."

We changed the wording and made additions to the monitoring description to make it clearer: “To reduce potential missing values, minimize irregularities and to increase adherence, the data sampling was constantly monitored by study members.”

– Line 432: The version of R is more relevant than the version of R Studio but I suppose that 4.1.1 already refers to the version of R.

We agree with the reviewer’s comment on R, we, in fact, used R version 4.1.1 and changed it within the text: “Statistical analyses were conducted using R studio (R version 4.1.1).”

Reviewer #2 (Recommendations for the authors):

The authors were seeking to investigate whether affectionate touch was associated with mental health during the pandemic using self-report measures and hormone levels.

A major strength of the paper was the large sample size. A weakness is that the study is observational so causation cannot be inferred. The authors achieved their aims and results support their conclusions. The paper is important as it shows the association between affectionate touch and lower distress during social isolation.

We thank the reviewer for this positive feedback!

Reviewer #3 (Recommendations for the authors):

The research described in this paper mainly used online surveys which attempted to establish whether attitudes to and receipt of affectionate touch from others might act as a buffer to reduce stress in challenging circumstances with reduced social contact such as during the COVID-19 lockdowns.

The main strength of the study is in its large scale and systematic and rigorous analysis of the measures collected online and also the detailed analysis of behavioral and endocrine measures in a smaller cohort of subjects at multiple time points over two days. The main weakness resides in insufficient measures being included to clearly establish the importance of different types of affectionate touch, and their context, and whether other social (i.e. social interactions not involving touch) or physical factors (such as activity/exercise) likely to be influenced during the lockdown period showed similar associations with stress, anxiety and endocrine measures.

Overall, the findings reported by the authors suggest firstly that individuals with a positive attitude to social touch are more likely to exhibit stress and loneliness problems in circumstances where social contact is reduced, such as during the COVID-19 lockdowns. Secondly, findings provide some initial support for the potential importance of receiving affectionate touch in helping to reduce stress and anxiety and influence concentrations of hormones associated with their regulation (i.e. cortisol and oxytocin). However, further more detailed and controlled experiments will be needed to establish fully the specificity of these potentially beneficial effects of affectionate touch, the most important aspects and contexts of affectionate touch which are potentially influencing wellbeing and the mechanisms involved.

Touch, particularly in the form of massage, has well established anxiolytic and other effects on well-being, although the current study has extended this to demonstrate that the amount and intensity of interpersonal affectionate received in everyday contexts may also be important. Giving and receiving warm hugs and cuddles may be more important to our general mental wellbeing than we realize.

We appreciate the reviewer’s thorough assessment of our work and the constructive comments.

In terms of the measures recorded while the focus is on affectionate touch do the authors also have any data on other social (face to face conversations and other mutual activities not involving intimate touch) and physical activity likely to have some effect on cortisol and oxytocin concentrations?

We thank the reviewer for this question. We included the variable of physical activity in our initial models to control for its potential influence on hormonal outcomes and did not find any significant associations. Indeed, in this paper we focused on the associations of affectionate touch in times of physical contact restrictions. To avoid α error inflation, originally, we did not conduct exploratory analyses. However, in response to the reviewer’s comment, we ran additional analyses on a specific other social activity. We used conversations as a predictor, as we think this could serve as an interesting comparison to touch. We analyzed whether a conversation (yes/no) was associated with hormonal levels and did not find significant associations, neither with cortisol nor oxytocin levels. In the following and equivalent to the touch data, we analyzed the associations of hormonal levels with the intensity of conversations and found a significant positive link between the intensity of conversation and oxytocin levels on within-person level, but not with cortisol concentrations. However, it must be noted that both variables, conversations, and affectionate touch, are highly correlated. If the editor and reviewer recommend us to report these findings, we can add the analyses in supplementary files to the manuscript.

While online information collection clearly has limitations concerning the amount of different variables which can be included I feel that given the central importance of the question concerning the role of receipt of affective touch many relevant details appear not to have been either included or considered. For example, were all subjects in a relationship and co-habiting/sleeping together? Was only one individual per household recruited? Given that a lot of affective touch can also be received from children and even pets, was this also considered? If subjects included a range of individuals both living alone and cohabiting I would have anticipated this being included as a variable, especially given the focus on COVID where individuals living alone might be expected to experience the greatest reduction in affective touch relative to those co-habiting/living in a family environment. As far as I can see there is no indication as to whether the subjects included in the study considered that they did receive less affectionate touch during the sampling period relative to before it?

Although we asked participants whether they received affectionate touch including hugs, kisses, cuddles, caresses, etc., it is not possible to extract the source of touch from this data.

Regarding relationship and cohabitation status, about 29 % of the subjects were single and 21% were living alone. We did not get into more detail within this paper as the focus of this study was to analyze within-subject associations of affectionate touch and hormonal and self-report measures. Nevertheless, we agree with the reviewer that living arrangements and relationship status might make a difference, especially during Covid-19 lockdowns, where social touch outside the own household was not permitted. In response to the reviewer’s comment, we considered the variables in our analyses. The group of individuals living with others reported more affectionate touch (514 times) than individuals living alone (42 times). Similarly, those in a relationship reported more affectionate touch (525 times) than singles (30 times). Within-subject association between touch, cortisol or oxytocin did not change overall, when relationship status or living situation were included. However, the between-subject association between affectionate touch and cortisol was mainly driven by the group with romantic partner. These results are linked with the reviewer’s earlier comment and suggests that those in a romantic relationship exchanged other types and more intense and erotic touch than singles – which again is in line with the association of touch intensity and oxytocin found here.

We unfortunately have only indirect information on whether exchanged touch had changed with the pandemic: via the online assessment participants were asked whether the quality of social contacts had changed since the beginning of the pandemic. Within the EMA sample, about 71% of the participants reported that the perceived quality of their social contacts was lower compared to before the pandemic. Based on this item, we cannot conclude whether the affectionate touch was received less often, we can only speculate that the reported reduction in the quality of social contact could possibly represent the reduction of affectionate touch.

The peak of oxytocin concentrations and reduced cortisol ones specifically at awakening is interesting but I wonder whether this simply reflects couples experiencing morning cuddles in a warm bed etc or a specific effect of awakening?

We thank the reviewer for bringing up this important question. To test this hypothesis, we now graphed oxytocin levels throughout the day, but separated the data based on the participants’ living situation relationship status (single vs. in a relationship) or (living alone vs. living with others). Of note, the patterns look almost identical (Figure 2—figure supplement 1). Even though we cannot rule out that there might be some other (unknown) factors responsible for the oxytocin “awakening response”, it seems that it is not exclusively due to the presence of a loved one.

We now also mention this in the manuscript (p. 13) and add these graphs as a supplement to our manuscript:

“The patterns of daily profiles did not appear to differ based on participants’ relationship status (single vs. in a relationship) or living situation (alone vs. with others).”

The reported saliva oxytocin concentrations are high and this is an ongoing issue in the field. I assume that the samples were run unextracted and also not in duplicate? While there is a clear recommendation now from both researchers and assay manufacturers to use extraction for plasma samples there is less consensus concerning saliva ones. With the widely used ELISA ENZO assay used in the current paper findings clearly shown that in extracted saliva samples generally result in baselines of <10pg/ml (see for example Fujii et al. 2016 Sci Reps 6:38662; Le et al. 2022 Psychother and Psychosom 91(5):335-347), while for unextracted ones this figure is much higher (up to 200 pg/ml), as in the current study. With plasma samples correlations between extracted and unextracted samples are weak, however a recent study comparing baseline and stimulated concentrations in saliva concentrations in post-partum women did report a significant low correlation (0.32-0.38) although profiles of release were different and concentrations were 10-fold higher in unextracted samples (Nagahasi-Araki et al. 2022, BMC Pregnancy and Childbirth 22:711). It is likely that higher concentrations in unextracted samples are contributed to by both the assay type (RIA vs ELISA) and whether the antibody used is sensitive to oxytocin fragments as well as oxytocin itself. Using a gold standard RIA, for example, as opposed to an ELISA, papers have consistently reported similar low oxytocin concentrations in both extracted plasma samples and unextracted saliva samples (generally below 5pg/ml – de Jong et al. 2015 Psychoneuroendocrinology 62:381-388; Martins et al., 2020 eLife 9:e62456), whereas with the ENZO ELISA extracted plasma and extracted saliva samples yield similar concentrations (Le et al., 2022). Indeed, there is a clear expectation that oxytocin concentrations should be similar in plasma and saliva since there is not reason to expect radically higher concentrations in saliva. The authors really need to at least acknowledge this issue and include it as a potential limitation. A recent review has, for example, highlighted the need for papers to make it very clear whether extracted or unextracted values are being reported (Tabak et al. 2023. Mol Psychiatry 28:127-140).

We thank the reviewer for these thoughtful comments. We agree that the information about whether the saliva samples were extracted prior to analyses is important and included this information in the methods section. We further discuss the methodological issues and debates regarding peripheral measures of oxytocin and include the recent reviews of Tabak et al., 2023 and Martins et al., 2020 that were mentioned by the reviewer (p. 19-20):

“Moreover, in the last decade, there has been an extensive discussion about the reliability and validity of peripheral oxytocin measures (Martins et al., 2020, Szeto et al., 2011, Tabak et al., 2023). There are several methodological issues and challenges associated with measuring oxytocin in blood plasma and saliva samples (Tabak et al., 2023). For example, studies have shown that oxytocin concentrations after sample extraction are much lower compared to unextracted oxytocin measurements (Szeto et al., 2011). Additionally, the correlations between extracted and unextracted oxytocin levels as well as between saliva and plasma oxytocin concentrations have been inconsistent across studies (e.g. Hoffman et al., 2012, Martins et al., 2020, Nagahashi-Araki et al., 2022, Szeto et al., 2011). These inconsistencies might be due to numerous reasons including differences in the study population, methods of sample processing and analyses. In particular, using different assay types (e.g. radioimmunoassay vs. enzyme immunoassay), as well as sample preparation (extraction vs. non-extraction), may contribute to these inconsistencies and make it difficult to compare results between studies (Tabak et al., 2023). Since unextracted samples were used in this study, the concentrations probably represent both free and bound oxytocin (MacLean et al., 2019) thereby potentially limiting the comparability with studies using extracted samples. Another important issue is the intraindividual stability of oxytocin over time (Feldman et al., 2013, Martins et al., 2020, Schneiderman et al., 2012). A recent study reports no correlation of single oxytocin measures between several assessments indicating that single measures of oxytocin might be not reliable to represent oxytocin baseline levels (Martins et al., 2020). In our sample, we found a significant positive correlation of mean oxytocin values between the two assessment days. As fluctuations of oxytocin throughout the day were apparent in our study, correlating mean values of six oxytocin measures over the day might represent the individual baseline oxytocin levels better than single measures (this has also been discussed in Tabak et al., 2023). However, the difference between our data and the findings reported by Martins et al., 2020 might also reflect methodological differences, as previous studies with unextracted samples also reported positive correlations over time (Feldman et al., 2013, Schneiderman et al., 2012).”

The authors do not mention a relevant study on humans reporting that foot massage both increases plasma oxytocin concentrations and corresponding neural responses in the orbitofrontal cortex and superior temporal sulcus but not somatosensory cortex. The study also reports that greater basal and manual touch-evoked oxytocin concentrations were associated with more positive attitudes to social touch (also using the STQ).(Li et al. 2019 Psychoneuroendocrinology 101:193-203).

We thank the reviewer for bringing this relevant study to our attention and have now included this paper in our introduction (p. 4):

“Similarly, a more recent study found a significant increase in plasma oxytocin and corresponding neural responses after a foot massage. Interestingly, basal oxytocin concentrations, as well as oxytocin increase after the massage, were associated with more positive attitudes towards social touch.”

As well as in the Discussion section (p. 18):

“Interestingly, a recent study demonstrated that foot massage was rated as more pleasurable and rewarding and was associated with a higher increase of oxytocin after the massage administered by hand, as compared to machine-administered massage, although the intensity of the massage was rated similarly (Li et al., 2019).”

While the authors have mainly decided to focus on presenting their findings and avoided entering into mechanistic explanations of how increased saliva oxytocin concentrations might influence levels of stress and anxiety etc I think it would help if there was a little more explanation of how saliva changes might influence the brain and endocrine system. In particular, the authors may wish to discuss the conclusion in another eLife paper that salivary and plasma oxytocin concentrations are not reliable trait markers for the oxytocin system in humans (Martins et al. 2020, eLife 9:62456).

We agree with the reviewer that this is an additional important point of debate. We are thankful to you for bringing up this paper. One of the arguments, that oxytocin measures in the periphery are not reliable, was that the authors did not find any correlation between oxytocin levels from different time points. We tested how the mean values from day 1 were correlated with day 2 and found a high correlation for both oxytocin and cortisol between days. We have now included that information in the manuscript (p. 13):

“A positive correlation between the two assessment days was found for individual (ln-transformed) mean values of oxytocin (r(227)=.850, p<.001) as well as cortisol (r(243)=.571, p<.001) levels.” We further discuss this paper, along with other papers regarding methodological issues of oxytocin assessments, in our Discussion section (text in the comment above).

We think that our data suggest a specific diurnal pattern in salivary oxytocin, which might have not been captured in previous studies. However, please note that we interpret salivary oxytocin levels as a consequence, rather than a driving factor of central nervous system phenomena, such as anxiety, stress or HPA axis activation. Recently, there has been data on touch and central nervous activation of oxytocin in female rats (Tang et al., 2020). We assume that such dynamics might predict peripheral oxytocin increases. We discuss this in more detail in the Discussion section (p. 21).

In the discussion the authors state that affectionate touch reduced the COVID-19 related burden (line 313) but as far as I can see there was no significant effect found for this?

We agree with the reviewer and corrected this mistake by changing the word “Covid-19 related” into “general”.

A small point but in the discussion of the study on rats by Tang et al. on line 326 animals are described as receiving "affectionate touch" from each other whereas perhaps simply receiving "social touch" might be more appropriate in this case.

We agree that in this case “social touch” is more appropriate and exchanged the words in the manuscript.

Additionally added references:

Bekhet, A. K., Zauszniewski, J. A., & Nakhla, W. E. (2008). Loneliness: a concept analysis. Nursing forum,

Feldman, R., Gordon, I., Influs, M., Gutbir, T., & Ebstein, R. P. (2013). Parental oxytocin and early caregiving jointly shape children’s oxytocin response and social reciprocity. Neuropsychopharmacology, 38(7), 1154-1162.

Hoffman, E. R., Brownley, K. A., Hamer, R. M., & Bulik, C. M. (2012). Plasma, salivary, and urinary oxytocin in anorexia nervosa: a pilot study. Eating Behaviors, 13(3), 256-259.

Hopf, D., Schneider, E., Eckstein, M., Aguilar-Raab, C., & Ditzen, B. (2021). COVID-und Social Distancing bezogene Sorgen und ihre Beziehung zu psychischer und körperlicher Erkrankung. PPmPPsychotherapie· Psychosomatik· Medizinische Psychologie, 71(02), 57-60.

Jacobi, F., Höfler, M., Strehle, J., Mack, S., Gerschler, A., Scholl, L., Busch, M. A., Maske, U., Hapke, U., & Gaebel, W. (2014). Psychische störungen in der allgemeinbevölkerung. Nervenarzt, 85(1), 77-87.

Li, Q., Becker, B., Wernicke, J., Chen, Y., Zhang, Y., Li, R., Le, J., Kou, J., Zhao, W., & Kendrick, K. M. (2019). Foot massage evokes oxytocin release and activation of orbitofrontal cortex and superior temporal sulcus. Psychoneuroendocrinology, 101, 193-203.

MacLean, E. L., Wilson, S. R., Martin, W. L., Davis, J. M., Nazarloo, H. P., & Carter, C. S. (2019). Challenges for measuring oxytocin: The blind men and the elephant? Psychoneuroendocrinology, 107, 225-231.

Martins, D., Gabay, A. S., Mehta, M., & Paloyelis, Y. (2020). Salivary and plasmatic oxytocin are not reliable trait markers of the physiology of the oxytocin system in humans. ELife, 9, 1-19. https://doi.org/https://doi.org/10.7554/eLife.62456

Musardo, S., Contestabile, A., Knoop, M., Baud, O., & Bellone, C. (2022). Oxytocin neurons mediate the effect of social isolation via the VTA circuits. ELife, 11.

Nagahashi-Araki, M., Tasaka, M., Takamura, T., Eto, H., Sasaki, N., Fujita, W., Miyazaki, A., Morifuji, K., Honda, N., & Miyamura, T. (2022). Endogenous oxytocin levels in extracted saliva elevates during breastfeeding correlated with lower postpartum anxiety in primiparous mothers. BMC Pregnancy and Childbirth, 22(1), 1-10.

Quintana, D. S., Rokicki, J., van der Meer, D., Alnæs, D., Kaufmann, T., Córdova-Palomera, A., Dieset, I., Andreassen, O. A., & Westlye, L. T. (2019). Oxytocin pathway gene networks in the human brain. Nature communications, 10(1), 1-12.

Schneiderman, I., Zagoory-Sharon, O., Leckman, J. F., & Feldman, R. (2012). Oxytocin during the initial stages of romantic attachment: relations to couples’ interactive reciprocity. Psychoneuroendocrinology, 37(8), 1277-1285.

Stalder, T., Kirschbaum, C., Kudielka, B. M., Adam, E. K., Pruessner, J. C., Wüst, S., Dockray, S., Smyth, N., Evans, P., & Hellhammer, D. H. (2016). Assessment of the cortisol awakening response: Expert consensus guidelines. Psychoneuroendocrinology, 63, 414-432.

Szeto, A., McCabe, P. M., Nation, D. A., Tabak, B. A., Rossetti, M. A., McCullough, M. E., Schneiderman, N., & Mendez, A. J. (2011). Evaluation of enzyme immunoassay and radioimmunoassay methods for the measurement of plasma oxytocin. Psychosomatic medicine, 73(5), 393-400.

Tabak, B. A., Leng, G., Szeto, A., Parker, K. J., Verbalis, J. G., Ziegler, T. E., Lee, M. R., Neumann, I. D., & Mendez, A. J. (2023). Advances in human oxytocin measurement: Challenges and proposed solutions. Molecular psychiatry, 28(1), 127-140.

Tang, Y., Benusiglio, D., Lefevre, A., Hilfiger, L., Althammer, F., Bludau, A., Hagiwara, D., Baudon, A., Darbon, P., Schimmer, J., Kirchner, M. K., Roy, R. K., Wang, S., Eliava, M., Wagner, S., Oberhuber, M., Conzelmann, K. K., Schwarz, M., Stern, J. E.,... Grinevich, V. (2020). Social touch promotes interfemale communication via activation of parvocellular oxytocin neurons. Nature neuroscience, 23(9), 1125-1137.

https://doi.org/10.7554/eLife.81241.sa2

Article and author information

Author details

  1. Ekaterina Schneider

    1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Heidelberg, Germany
    2. Heidelberg University, Heidelberg, Germany
    Contribution
    Data curation, Formal analysis, Investigation, Visualization, Methodology, Writing – original draft
    Competing interests
    No competing interests declared
  2. Dora Hopf

    1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Heidelberg, Germany
    2. Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Data curation, Formal analysis, Funding acquisition, Investigation, Methodology, Project administration, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9476-0478
  3. Corina Aguilar-Raab

    1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Heidelberg, Germany
    2. Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  4. Dirk Scheele

    Department of Social Neuroscience, Faculty of Psychology, Ruhr University Bochum, Bochum, Germany
    Contribution
    Conceptualization, Resources, Investigation, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7613-0376
  5. Andreas B Neubauer

    1. Department for Education and Human Development, DIPF|Leibniz Institute for Research and Information in Education, Frankfurt, Germany
    2. Center for Research on Individual Development and Adaptive Education of Children at Risk, Frankfurt, Germany
    Contribution
    Formal analysis, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  6. Uta Sailer

    Department of Behavioural Medicine, Faculty of Medicine, Institute of Basic Medical Sciences, University of Oslo, Oslo, Norway
    Contribution
    Formal analysis, Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-9728-8738
  7. René Hurlemann

    Department of Psychiatry, University of Oldenburg, Bad ZwischenahnOldenburg, Germany
    Contribution
    Methodology, Writing – review and editing
    Competing interests
    No competing interests declared
  8. Monika Eckstein

    1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Heidelberg, Germany
    2. Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Methodology, Project administration, Writing – review and editing
    For correspondence
    monika.eckstein@med.uni-heidelberg.de
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1846-4992
  9. Beate Ditzen

    1. Institute of Medical Psychology, Center for Psychosocial Medicine, Heidelberg University Hospital, Heidelberg, Germany
    2. Heidelberg University, Heidelberg, Germany
    Contribution
    Conceptualization, Resources, Software, Formal analysis, Supervision, Funding acquisition, Validation, Investigation, Methodology, Project administration, Writing – review and editing
    For correspondence
    beate.ditzen@med.uni-heidelberg.de
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-5853-4572

Funding

German Research Foundation (SFB 1158)

  • Beate Ditzen

German Psychological Society (Corona Scholarship)

  • Ekaterina Schneider

German Research Foundation (Young Seed Grant)

  • Monika Eckstein

German Academic Exchange Service

  • Monika Eckstein

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The study team is grateful for the help with data collection received from C Gäbel, M Fischer, L-M Müller, N Stockburger, L Fischer, R Dahlke, F Frech, J Gerber, and J Zimmer. For the publication fee we acknowledge financial support by Deutsche Forschungsgemeinschaft within the funding programme „Open Access Publikationskosten“ as well as by Heidelberg University.

Ethics

For this study ethical approval was granted from the ethics committee of the Heidelberg University Medical Faculty (approval no. S-214/2020) and the study was registered online at https://drks.de/search/en/trial/DRKS00021671. All participants provided written informed consent.

Senior Editor

  1. Christian Büchel, University Medical Center Hamburg-Eppendorf, Germany

Reviewing Editor

  1. Matthias Gamer, University of Würzburg, Germany

Reviewers

  1. Matthias Gamer, University of Würzburg, Germany
  2. Elizabeth Anne Broadbent
  3. Keith M Kendrick, University of Electronic Science and Technology of China, China

Version history

  1. Preprint posted: March 28, 2022 (view preprint)
  2. Received: June 20, 2022
  3. Accepted: May 12, 2023
  4. Version of Record published: May 30, 2023 (version 1)

Copyright

© 2023, Schneider et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 378
    Page views
  • 51
    Downloads
  • 0
    Citations

Article citation count generated by polling the highest count across the following sources: Crossref, PubMed Central, Scopus.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Ekaterina Schneider
  2. Dora Hopf
  3. Corina Aguilar-Raab
  4. Dirk Scheele
  5. Andreas B Neubauer
  6. Uta Sailer
  7. René Hurlemann
  8. Monika Eckstein
  9. Beate Ditzen
(2023)
Affectionate touch and diurnal oxytocin levels: An ecological momentary assessment study
eLife 12:e81241.
https://doi.org/10.7554/eLife.81241

Further reading

    1. Epidemiology and Global Health
    2. Genetics and Genomics
    Arturo Torres Ortiz, Michelle Kendall ... Louis Grandjean
    Research Article

    Accurate inference of who infected whom in an infectious disease outbreak is critical for the delivery of effective infection prevention and control. The increased resolution of pathogen whole-genome sequencing has significantly improved our ability to infer transmission events. Despite this, transmission inference often remains limited by the lack of genomic variation between the source case and infected contacts. Although within-host genetic diversity is common among a wide variety of pathogens, conventional whole-genome sequencing phylogenetic approaches exclusively use consensus sequences, which consider only the most prevalent nucleotide at each position and therefore fail to capture low frequency variation within samples. We hypothesized that including within-sample variation in a phylogenetic model would help to identify who infected whom in instances in which this was previously impossible. Using whole-genome sequences from SARS-CoV-2 multi-institutional outbreaks as an example, we show how within-sample diversity is partially maintained among repeated serial samples from the same host, it can transmitted between those cases with known epidemiological links, and how this improves phylogenetic inference and our understanding of who infected whom. Our technique is applicable to other infectious diseases and has immediate clinical utility in infection prevention and control.

    1. Epidemiology and Global Health
    2. Genetics and Genomics
    Jenny van Dongen, Gonneke Willemsen ... Michael C Neale
    Research Article Updated

    Background:

    Smoking-associated DNA methylation levels identified through epigenome-wide association studies (EWASs) are generally ascribed to smoking-reactive mechanisms, but the contribution of a shared genetic predisposition to smoking and DNA methylation levels is typically not accounted for.

    Methods:

    We exploited a strong within-family design, that is, the discordant monozygotic twin design, to study reactiveness of DNA methylation in blood cells to smoking and reversibility of methylation patterns upon quitting smoking. Illumina HumanMethylation450 BeadChip data were available for 769 monozygotic twin pairs (mean age = 36 years, range = 18–78, 70% female), including pairs discordant or concordant for current or former smoking.

    Results:

    In pairs discordant for current smoking, 13 differentially methylated CpGs were found between current smoking twins and their genetically identical co-twin who never smoked. Top sites include multiple CpGs in CACNA1D and GNG12, which encode subunits of a calcium voltage-gated channel and G protein, respectively. These proteins interact with the nicotinic acetylcholine receptor, suggesting that methylation levels at these CpGs might be reactive to nicotine exposure. All 13 CpGs have been previously associated with smoking in unrelated individuals and data from monozygotic pairs discordant for former smoking indicated that methylation patterns are to a large extent reversible upon smoking cessation. We further showed that differences in smoking level exposure for monozygotic twins who are both current smokers but differ in the number of cigarettes they smoke are reflected in their DNA methylation profiles.

    Conclusions:

    In conclusion, by analysing data from monozygotic twins, we robustly demonstrate that DNA methylation level in human blood cells is reactive to cigarette smoking.

    Funding:

    We acknowledge funding from the National Institute on Drug Abuse grant DA049867, the Netherlands Organization for Scientific Research (NWO): Biobanking and Biomolecular Research Infrastructure (BBMRI-NL, NWO 184.033.111) and the BBRMI-NL-financed BIOS Consortium (NWO 184.021.007), NWO Large Scale infrastructures X-Omics (184.034.019), Genotype/phenotype database for behaviour genetic and genetic epidemiological studies (ZonMw Middelgroot 911-09-032); Netherlands Twin Registry Repository: researching the interplay between genome and environment (NWO-Groot 480-15-001/674); the Avera Institute, Sioux Falls (USA), and the National Institutes of Health (NIH R01 HD042157-01A1, MH081802, Grand Opportunity grants 1RC2 MH089951 and 1RC2 MH089995); epigenetic data were generated at the Human Genomics Facility (HuGe-F) at ErasmusMC Rotterdam. Cotinine assaying was sponsored by the Neuroscience Campus Amsterdam. DIB acknowledges the Royal Netherlands Academy of Science Professor Award (PAH/6635).