Next Article in Journal
Identification and Validation of Metabolic Markers for Adulteration Detection of Edible Oils Using Metabolic Networks
Previous Article in Journal
Improved Dried Blood Spot-Based Metabolomics: A Targeted, Broad-Spectrum, Single-Injection Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Chemical Fingerprints of Emotional Body Odor

1
Unilever R&D Vlaardingen, Olivier van Noortlaan 120, 3133 AT Vlaardingen, The Netherlands
2
Faculty of Social and Behavioural Sciences, Utrecht University, Heidelberglaan 1, 3584 CS Utrecht, The Netherlands
3
William James Center for Research, ISPA-Instituto Universitário, 1149-041 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Metabolites 2020, 10(3), 84; https://doi.org/10.3390/metabo10030084
Submission received: 2 February 2020 / Revised: 22 February 2020 / Accepted: 24 February 2020 / Published: 28 February 2020

Abstract

:
Chemical communication is common among animals. In humans, the chemical basis of social communication has remained a black box, despite psychological and neural research showing distinctive physiological, behavioral, and neural consequences of body odors emitted during emotional states like fear and happiness. We used a multidisciplinary approach to examine whether molecular cues could be associated with an emotional state in the emitter. Our research revealed that the volatile molecules transmitting different emotions to perceivers also have objectively different chemical properties. Chemical analysis of underarm sweat collected from the same donors in fearful, happy, and emotionally neutral states was conducted using untargeted two-dimensional (GC×GC) coupled with time of flight (ToF) MS-based profiling. Based on the multivariate statistical analyses, we find that the pattern of chemical volatiles (N = 1655 peaks) associated with fearful state is clearly different from that associated with (pleasant) neutral state. Happy sweat is also significantly different from the other states, chemically, but shows a bipolar pattern of overlap with fearful as well as neutral state. Candidate chemical classes associated with emotional and neutral sweat have been identified, specifically, linear aldehydes, ketones, esters, and cyclic molecules (5 rings). This research constitutes a first step toward identifying the chemical fingerprints of emotion.

1. Introduction

The type of information human chemosignals (body odors) carry has been the subject of an increasing number of systematic investigations. The accumulating evidence from this research has highlighted the significant role that the neglected medium of olfactory communication (chemosignals) and the sensory system (the sense of smell) play in communication (e.g., [1,2,3]). This research has spanned the effects of body odors produced under different conditions of the donors, ranging from illness (e.g., [4]) to emotional states (e.g., fear: [5]). Despite the rapidly accumulating research on the effect of the information contained in body odors and its range, there is to our knowledge, no study that has analyzed the chemical composition of body odors and identified the distinctive chemical differences in the features of, for instance volatiles produced in a fearful versus happy state. This is precisely the focus of the current report. Below we first provide a very brief overview of the research on the information transmitted via human body odors. We then turn to the nature of the chemical composition of human odors and its characteristics that is the focus of the empirical work reported here.
“Communication occurs when a signal produced by an individual cause a reaction in another organism where both the signal and the reaction have been designed for these purposes” (Scott-Philips, 2015, p. 30; see [6]). Body odors have been found to serve precisely this signal function, and can range from bodily fluids such as sweat, to urine, tears, blood, and mother’s milk. The accumulating research to date has shown that they can signal information about the states and traits of individuals, such as emotional state [7,8,9,10], illness [4], gender [8,11,12], sexual preference [13,14], sexual arousal [15], kin [16], mate compatibility [17], menstrual phase [18] albeit some controversy on this subject [19,20], and personality traits such as extraversion [21]. The effects of human body odors do not appear to adhere to the strict definition of pheromones, as being single molecular compounds that activate very specific reactions (e.g., stereotyped behaviors, or developmental processes [22]). Therefore, the type of communication carried by human odors is increasingly referred to as chemosignaling [10,23].
Considerable progress has been made in our understanding of how emotions are chemosignaled, namely how emotion is transferred from one agent to another by means of body odors. Usually, these examinations used an experimental paradigm in which donors are exposed to emotion either inducing scenarios (e.g., film clips inducing fear, happiness, disgust states) or situations (e.g., stressful exams). During these sessions, the excreted odor is sampled from the armpit region. The collected samples are frozen and retained until the experimental session are thawed, and a recipient agent is then exposed to the odor. The experimental paradigms used a number of dependent measures depending on the type of emotional condition under which the donor produced the armpit odor and the predicted neural, cognitive, and behavioral consequences for the recipient. For example, in the case of fear odor one general question has been its consequences for increased vigilance (e.g., [24,25]), which has been measured by the caution manifested in a word choice task [24], the accuracy on an easy visual search task [25], the speed with which facial expressions were classified [26]. Another frequently used indicator of armpit sweat effects is facial muscle activity measurement (electromyography) indicative of facial emotion expressions [1,8,26,27,28]. Other assessments have focused on effects of chemical signals in relation to real-life settings such as influencing social judgment [29], and clinical performance of dental students [30]. Similarly, recordings of brain activity have been assessed with findings indicating higher event-related potentials at various time points (e.g., [31,32]), increased amygdala activity [33,34], and recruiting a wider network of brain regions in relation to empathy [35].
While it is clear that human chemosignals emitted by the sender are carriers of the information that is transmitted and received by the olfactory sensory system, the nature of the carriers remains a mystery. Body odor emanating from the underarm consists of a complex mixture of hundreds of chemical volatiles that is perceived holistically [36]. It is possible that a specific emotion signal is conveyed by a distinctive pattern of volatiles within an entire “package” of a volatile complex [36]. Thus, the dynamic nature of emotional states could mean that a distinctive set of volatiles are emitted or modified, generating in the brain a representation of that odor as being distinctively different from other such odors. An example not related to body odor would be the smell of coffee or wine, which are smelled by many people as distinctive and unique to the beverages that generate them, even though these smells emerge from chemical cocktails made up out of hundreds of volatile compounds [37]. We surmise that the emotion signal in body odor is conveyed by a pattern of chemical volatiles embedded in the larger mixture of volatiles, which in their combination may be regarded as characteristic of the emotional state experienced by the sender. It is unlikely that emotional specificity in body odor depends on a single (magical) molecule [2,22,23,38]. So rather than a single “fear” or “happiness” molecule, the more prevalent view about the nature of odor-print or odor signature is that there is a larger bouquet of chemicals shared across individuals that differ in terms of relative amounts [36]. This is why we introduced the concept of the emotional odor signals as patterns of chemical volatiles rather than single key molecules.
The current study was designed to examine the feasibility of emotion-specific odor signals or fingerprints in underarm sweat. First, we examined whether the patterns of volatile chemicals associated with underarm sweat samples from fearful, happy, and neutral donors, respectively, are distinctly different. Second, based on a first discrimination between chemical patterns, we tentatively identified the molecular classes associated with emotional and control sweat. An overview of volatile compounds emitted from the skin by Dormont et al. [39] lists “carboxylic acids of various chain lengths and derivative esters, aldehydes, alkanes, short chain alcohols, and some ketones” [39] with mainly alkanes and C6–C11 carboxylic acids being associated with the axillae. In another review, de Lacy Costello et al. [40] identified 532 named volatile compounds analyzed from skin secretions, although not in relation to the axillae specifically. The study reported here does not pretend to elucidate the definitive fingerprints of fear or happy emotion—the very existence of which might even be the topic of debate—yet it aims to progress our understanding of the chemical carriers of emotion communication by focusing on the primary candidate chemical classes, which may include any of the above.
We sampled underarm sweat from 24 male donors (for full details please see Methods) in a within-subjects design in which each donor participated in three emotion induction sessions (fear, happy, (pleasant) neutral) with a one-week interval between sessions during which the donors followed a protocolled regimen to keep contamination of sweat odor from odors related to other sources to a minimum. The emotion induction sessions involved watching pre-validated emotional movie clips for approximately half an hour. Following previous research [1,8,26,27,28], and for reasons explained in the Methods section, the order of sessions was not random, always starting with fear and ending with neutral. In the happy condition, as opposed to the fear and neutral condition, donors sat together in groups of 2, 3, or 4 to boost the positive experience. This was done because the co-presence of others heightens the positive emotional experience [41,42,43,44]. Self-report rating of emotion states was collected to check the effectiveness of the manipulation. The sweat was collected onto textile pads attached to the left and right underarm of pre-cleaned t-shirts. Textile pads were removed from the t-shirts and kept in sterile vacuum packs at temperatures of maximally −24 °C until the time of analysis.
Chemical analyses were conducted at Unilever R&D Vlaardingen, the Netherlands. All 144 textile pads were subjected to HeadSpace Sorptive Extraction combined with Thermal Desorption, followed by separation and detection using GC×GCMS. The Statistical Compare software feature within ChromaTOF® was subsequently used to perform peak alignment between all samples. The resulting mass/retention time (MS/RT) peak table (3796 peaks) was exported and used for further multivariate statistical analysis. After data clean-up, 1655 volatiles were retained and subjected to principal component analysis (PCA), followed by partial least squares (using non-normalized, Pareto scaling, 7-fold cross-validation, and permutation testing) to identify statistically significant peaks testing three models: Fear vs. Neutral, Happy vs. Neutral, and Happy vs. Fear. The analysis focuses on commonalities of combined variables within a group and how they differ between the groups. Box plots were generated by using a statistical program (JMP®) for identifying and selecting volatiles showing significant differences.
We deliberately chose a conservative analysis that reveals only the strongest discriminators. As we did not know in advance which chemical features are associated with emotion signals, we believe than a comprehensive untargeted screening approach is appropriate.
Based on visual detection of a pattern of results corresponding to the Happy condition suggesting subclusters in the chemical data, cross-validated regularized logistic regression was conducted on the self-reported emotional states in all emotion conditions to explore whether these could account for these chemical subclusters. In other words, does the position donors occupy in the PCA score plot when in the Happy condition relate to their scores on the various self-reported emotional states?

2. Results

2.1. Emotion Manipulation

A total of 144 sweat pads were collected from the left and right underarm of 24 non-smoking Caucasian men (Mage = 22.42, SD = 3.54; range 19–34) during the three conditions Fear, Happy, and Neutral.
Self-report data on emotional state were first subjected to a nonparametric Friedman test (degrees of freedom: 2, n = 24) to check for overall effects of condition (fear, happy, neutral) on self-report variables. Follow-up Wilcoxon signed-ranks tests were used to distinguish between specific pairs of conditions (fear-happy, fear-neutral, happy-neutral). Table 1 displays the results on effects of emotion condition on self-reported emotional state, demonstrating that emotion induction via movies was successful in eliciting the three target states (see also Figure 1): Happiness, χ²(2, n = 24) = 40.77, p < 0.001; Fear, χ²(2, n = 24) = 38,00 p < 0.001; and Calmness, χ²(2, n = 24) = 33.11, p < 0.001. Following [45] we refer to calmness as an indicator of a pleasant neutral state.

2.2. Chemical Analyses

To identify GC-MS peaks that have systematically higher or lower concentrations in one condition group than in the other, partial least squares discriminant analysis (PLS-DA) was performed on two condition classes, respectively. All three significant PLS-DA models, namely Neutral vs. Fear, Neutral vs Happy, and Fear vs. Happy sweat (Figure 2) were statistically significant. However, the strongest model was the Neutral vs. Fear (Q2 = 0.85), followed by Neutral vs. Happy (Q2 = 0.6). The Fear versus Happy model was weaker but still significant (Q2 = 0.33). From each model, the most discriminative peaks were identified, amounting to 94 peaks in total. They are listed in Table A1 (Appendix A). To visualize the relationships of the three condition groups to each other, principal component analysis (PCA) was performed based on these 94 peaks (Figure 3). Individual box plots for each of the peaks across the emotion conditions can be found in Supplementary Figure S1. From this analysis, it became evident, that Neutral and Fear are distinct classes while the happy condition showed a bipolar pattern, with a subcluster of data overlapping with the “fear” cluster, and another subcluster with the “neutral” cluster. An explanation for the two subclusters was explored via post-hoc testing below.
Following the multivariate statistical analysis, the peaks found to be discriminative in the PLS-DA models were further analyzed using univariate statistics. For this, receiver operating characteristics (ROC) curves were calculated to measure the ability of each peak to discriminate between two types of emotion. In Figure A1 (Appendix A) two examples are depicted showing that peaks 2209 and 234 clearly discriminated Fear from Neutral. In contrast, the discriminative ability of these peaks was less when comparing Happy and Neutral, or Fear and Happy.
In addition, box plots comparing the emotion conditions on 94 statistically significant peaks were inspected, revealing clear differences between the conditions for 84 peaks. The remaining ten peaks showed only minor differences. From the 84 peaks, 36 peaks had increased intensities for both the Fear and Happy conditions when compared to Neutral (Figure 4A); 23 peaks had decreased intensities for both Fear and Happy conditions when compared to Neutral (Figure 4B); five peaks had increased intensities upon Fear when compared to Happy and Neutral (Figure 4C); 16 and 4 peaks were more increased (Figure 4D) and decreased (Figure 4E), respectively, in Happy when compared to Fear and Neutral. It is worth mentioning that most peak changes were not due to increased or decreased total sweat production, because there was (i) no significant difference in the total sweat production between the conditions (Figure A2, Appendix A) and (ii) no correlation between the peak intensities and the total intensities. Only for about 31 out of 94 peaks, the condition-dependent changes of the peak were related to increased total sweat production.
The data on 94 significant compounds were thereupon reviewed to investigate whether a logical chemical pattern could be detected to explain the chemical pattern of data. Two main types of chemical compound classes were identified: class 1 compounds (aldehydes/ketones) and class 2 compounds (esters and cyclic molecules (5 rings)). There was also a third class with 22 “happy compounds” that were mostly changed in the Happy when compared to Fear and Neutral. However, most of these compounds could not be identified or were identified as non-biological compounds, such as caprolactam. The Class 1 compounds (aldehydes/ketones) were high in the fear sweat and low in the neutral sweat condition and Class 2 compounds (esters and cyclic molecules) were high in the neutral condition and low in the fear sweat condition. The sweat of the subjects in the happy condition was either similar to the fear sweat or to the neutral sweat.
Further analysis on the subjects (n = 24) based on bipolar response in the happy condition resulted in two types of responses: response group 1 (n = 12) which we will refer to as sub-cluster HF (Happy overlapping with Fear) as it refers to individuals whose chemical patterns when in the happy condition overlap with the fear condition chemical patterns, and response group 2 (n = 8) which we will refer to as subcluster HN (Happy overlapping with Neutral) as their chemical patterns overlap with the neutral condition chemical patterns. This finding seems to suggest that the presumably emotional dimension that seems to be captured mostly by PC 1 and to some extent PC 2 (t[1] versus t[2] respectively in Figure 3) may be characterized by an increase in the linear aldehydes/ketones and a decrease in the esters and cyclic molecules (5 rings). That is, the high t[1](neutral both groups + subcluster HN) sweat in Figure 3 was high in esters and 5 rings, but low in linear aldehydes/ketones, while the low t[1] (fear both groups + subcluster HF) sweat was low in esters and 5 rings, but high in linear aldehydes/ketones (Table 2). More details about the members of the two chemical classes can be found in Table 3.

2.3. Exploring the Pattern of Bipolarity

Additional analyses were run to explore whether the bipolar pattern in the chemical data from the happy condition, referred to as HF and HN sub-clusters, can be explained by differences in emotional profile based on self-report, or on other factors. Figure 1A shows that profiles for self-reported emotional states for happy and neutral in particular, overlap, and hence it is possible that those individuals whose emotional state when happy overlapped to some extent with neutral, may have shown a correlated pattern for their chemical sweat odor profiles. We thereupon investigated whether the sub-clusters HF and HN differed in self-reported happiness. Happiness was virtually indistinguishable between the chemical sub-clusters, and very high on both (MHF = 3.79, SDHF = 0.42, MHN = 3.75, SDHN = 0.46, t(20) = 0.18, p = 0.856), so the difference in chemical configuration could not be explained by bimodality in happiness experienced by the donors. We thereupon investigated whether the bimodal pattern could be explained by systematic differences in emotional states other than happiness (see Table 1, Figure 1A). This turned out not to be the case. For these analyses and statistical approach we refer to Supplementary material S2A,B: Additional testing to explore bipolar pattern in happy chemical clusters.

3. Discussion

This study was conducted to explore the possibility of identifying chemical fingerprints of distinct emotional states in sweat odor sampled from the axillae of male donors in fear, happy, and neutral emotional states. Support was obtained for the notion of signal specificity in axillary sweat odor, as clusters of chemical data in a two-dimensional space spanning PC1 and PC2 obtained from PCA conducted on a matrix of 1655 chemical volatile peaks showed clearly distinctive clusters significantly separating fingerprints from fear versus neutral, as well as happy versus neutral, and happy versus fear. The distinction between the clusters showed up most clearly for the contrast fear versus neutral, while the pattern corresponding to the happy condition is best described as consisting of two subclusters, with one such cluster overlapping in chemical space with fear (HF), the other with neutral (HN). In general, the patterns are best described as Fear and HF being high in intensity of linear aldehydes and ketones, and low on esters and cyclic structures, and Neutral and HN showing the opposite pattern. Referring to the literature, the overview of chemical volatiles analyzed from the body compiled by de Lacy Costello et al. [40] as displayed in Figure 2 of their review [40] shows a large variety of chemical classes listing ten classes of chemical compounds for skin that include the four identified in our study, alongside others such as nitrogen- and sulfur-containing classes. Our analyses suggest that differences across emotion versus neutral sweat are most likely to be found in intensity differences of the aforementioned four chemical classes. Previously, using data mining techniques on a database of over 100 trace gas species, Williams et al. [48] had found differences in airborne chemicals emitted into the air of a cinema by audiences watching movies. Specific film events (i.e., “suspense” versus “comedy”) were associated with different levels of chemicals (such as acetone and isoprene) being emitted into the air. Furthermore, a recent publication by Preti et al. [49] comparing VOC composition of exhaled breath from individuals after validated stress induction to breath collected after a non-stressed control condition using GC/MS analytical techniques revealed three candidate biomarkers of stress: dimethyl sulfide and again, acetone, isoprene. It is not immediately clear how these results relate to our findings, as this would entail a comparison of stress-related volatiles in the oral cavity versus underarm, which are different environments. However, the results from both abovementioned studies support the notion of some specificity of volatile biomarkers relating to emotional state.
Previous reports on chemical analysis from odor arising from the axillary have also listed long-chain acids such as (E)-2-methyl-3-hexenoic-acid [50,51] as strong contributors to malodor. This compound, along with other typical compounds such hexanoic acid, butanoic acid, valeric acid, and 3-hydroxy-3-methylhexanoic acid (3H3MH), was picked up in our analyses as well but did not appear among the set of 94 compounds for which significant differences in peak intensity between the emotion conditions were found. This suggests that compounds previously identified as “smelly” are not necessarily also involved in chemosignalling of emotion.
In order to better understand the subclusters within the happy condition we examined the possibility of a relation to self-reported emotional state. However, we did not find support for a correspondence between these chemical clustering in space of the happy donors and their self-reported happiness or other emotional states. One could have expected a relation with arousal to emerge, as it is possible that some males experienced enhanced arousal in both the fear and happy states. This could have come about if the happy movie clips made them feel excited and invigorated in a way that corresponded to the arousal state resulting from the fear induction. Likewise, in other males the happy clips may have lowered the arousal levels comparable to the neutral condition because the happy clips had a soothing effect instead. However, this pattern did not come out. Self-reported arousal may not have been able to pick the proposed systematic relationship up. It would have been more appropriate to gage arousal with physiological endpoints such as heart rate and skin conductance.
Finally, we considered whether differences in movement across the conditions could account for the differences in patterns of volatiles. However, well-known indicators of energy metabolism like branched amino acids [52] and metabolites like acetoacetate and beta-hydroxybutyrate [53] are not volatile enough to detect with the used GC-MS technique.
Therefore, except for acetone which did not show up as different across conditions, neither acetoacetate nor beta-hydroxybutyrate were detected without further derivatization agents. However, assuming that the significant linear aldehydes and ketones we detected were correlated with movement, then the observed pattern across conditions does not seem to be entirely logical. While we encountered low levels of linear aldehydes and ketones during the emotionally neutral condition (as expected), these levels were high in the fear condition when participants did not move much. As noted earlier, in the happy condition 14 individuals showed a high response in terms of emitted volatiles whereas eight showed a low response. In conclusion, it is unclear when considering the present data whether differences in movement across the conditions caused the observed patterns of volatile emissions. Moreover, we would like to point out that rather than considering movement as a confound, we would consider movement (or lack thereof) as intrinsic to, or part of, the emotional states of interest. For future studies, it would be informative to include movement measures to explore the correlation between volatiles and movement in relation to emotional state.
There are a few limitations to our research: we used a fixed order of emotion induction for reasons that goes against the common practice of counterbalancing conditions. These reasons have been explained in the Methods section and are related to the inherent difficulty of inducing emotions in a laboratory setting knowing they are also influenced by pre-existing affective states [45]. Consequently, inducing fear during the first visit when cortisol levels are known to be high [54], served to maximize the chances of success at inducing fear. The self-reported levels of emotion from the donors in our study corroborated effective induction of fear, happiness, and a calm “pleasant neutral” state by following this fixed order.
To address the issue of fixed order, in a recently published study by our team [55] fear and pleasant-neutral were counterbalanced and sessions were held on the same day to simplify logistics. We found non-significant order effects on physiology and subjective ratings of fear. Counterbalancing also did not abolish the typical fear responses shown by receivers in a subsequent receiver study [56]. As this finding emerged after we had conducted our study, we recommend the use of counterbalanced designs in future research; however, in the present research, we wanted to use a cautious approach that mimicked prior research.
Another limitation may apply to the happy emotion induction condition. Here, happiness was induced while 2–4 donors were sitting together as opposed to being alone as was the case in the fear and neutral emotion induction conditions. The decision to do this was based on the desire to maximize the chances of success of happiness induction as the social presence of others can boost the positive experience, as pointed out earlier. The self-reported happiness/positive emotion ratings demonstrate that the happiness induction was indeed effective. However, the possibility remains that chemical volatiles associated with the body odor emanated during that session may have reflected “social” processes rather than “happy.” Data collected independent of the present study to address this issue are reported as Supplementary material S3. We conclude from the results of this study that effective emotion condition can emerge independently from the social presence or others suggesting that it is unlikely that this factor influenced the chemical composition of sweat as much as the experienced emotion.
Finally, and referring to Lenochova, Roberts, and Havlicek [57] who showed no differences in perceptual evaluations of sweat collected fresh versus after 4 months of freezing, we do not expect an effect of freezer storage time on sweat composition. In the present study, sweat was stored in the freezer at −24 °C directly after collection; all samples were collected over a one month period and shipped together to the analytical facility where they were stored at −80 °C and analyzed shortly after arrival. Thus, we conclude that the difference in composition of volatiles between emotional sweat clusters reflects the differences in emotional state induced in and reported by the donors as intended.
In previous research by our group, the emotional signal in the sweat volatiles was validated by the finding that receivers exposed to donors’ odor displayed a simulacrum in their facial expressions of the emotion experienced by the donors [1,8,25,26,27,28] The correspondence between donor and recipient could not be tested in the present study, because all sweat that was collected was used for the purpose of chemical analysis. Nevertheless, the evidence for such a correspondence has been documented in numerous studies. Notably, as a first step, differences in sweat composition correlating to emotional versus neutral state in sweat donors was confirmed.
Likewise, the biological activity of the extracted chemical characteristics for the different body orders can only be established if it were possible to artificially compose odors with the same properties that were extracted and show that these reproduce the reactions found earlier for the biological material collected (“Koch’s postulates”, see [22]). This is an avenue for future research.
As far as we know, this was the first study to investigate a systematic relation between emotional states and sweat odor emanated during such states. Our results underscore the likelihood of chemical fingerprints of emotion in sweat volatiles related to the emotional state. Specifically these differences emerge in the chemical compound classes of linear aldehydes, ketones, esters, and 5 rings as candidate key compounds. The distinctive fingerprints reported here support the possibility of emotion-specific signals excreted by a sender enabling a social communication with a receiver. Previous research by our group suggests that this chemical communication can elicit emotional contagion between the sender and receiver putting both individuals on the same page. These results again demonstrate the remarkable social functionality (see also Stevenson [58]) of the sense of smell, being not just limited to animals, but extending to humans [59].
In view of the exploratory nature of this research, we used a comprehensive, nontargeted screening approach. Using this conservative method, only the strongest discriminators emerged. As a next step, we call on others to conduct separate studies to replicate these results and validate our variables in a concerted effort of unravelling the signal of happiness and fear.

4. Materials and Methods

4.1. Emotion Induction and Sweat Collection from Donors

Twenty-four non-smoking Caucasian men (Mage = 22.52 years, SD = 3.50) were recruited at Utrecht University in 2014 to participate in this study. Participants were included if they reported experiencing fear during horror movies and not scoring extremely high on neuroticism or psychoticism. All men provided written informed consent to follow a regimen to prevent odor contamination of underarm sweat which prescribed adherence to a diet, the use of fragrance-free personal care products, and restriction of certain behaviors for 2 days prior to each sweat sampling session, and to donate sweat in three consecutive sessions (fear, happy, neutral) separated by one week. Order of emotion-induction conditions was fixed, and not counterbalanced because of known difficulties with the induction of emotions like fear in a laboratory setting. Emotion experiences elicited in the lab are not influenced by pre-existing affective states [45]. In previous experiments (e.g., [54]), we found that when participants enter the lab for the first time their salivary cortisol levels (indicator of psychological stress) were high (e.g., de Groot et al. [54], Figure 3C). As it can take more than 20–30 min for this level to decline to a “relaxed” state, being stressed on the first visit might confound the induction of happiness or neutral/calm emotional states, while at the same time, the fear induction would benefit from it. So, to maximize the chance of the fear induction being successful, we decided to administer it always first in the context of pre-existing high cortisol in initial donor sessions. For the inverse reason, the neutral condition was always last. At this point, we expected donors to have become sufficiently familiar with the procedure and the experimenter and that they would be fairly relaxed when they arrived. At the same time, because the first (fear induction when being alone) and the second session (happiness together with others) had been so different we also expected them to be still curious about what would happen in the third session to counteract the possible boredom.
During the second emotion induction, donors sat together in groups varying from 2–4 because the co-presence of others is contagious and heightens a positive experience. As stated above, inducing emotions in the lab can be difficult [45]. Rather than making an effort to keep all inductions similar on other domains (number of people being tested in a condition), we made efforts to design the induction so as to maximize the chances that the emotion induction would be successful.

4.1.1. Emotion Induction

Each of the emotion induction procedures involved watching a series of movie-clips for approximately 30 min. In the fear condition participants watched the original version of nine English-spoken horror film clips [60]; film database codes: 7, 16, 28, 32, 38, 46, 50, 55, 66; see also de Groot et al. 1. In the happy condition, participants watched four scenes carefully selected scenes (“Bear Necessities”, The Jungle Book; a short movie of Kurt Kuenne Validation; a comical opera scene from the Intouchables; and a televised prank on a company called Mobistar from the television programme Basta). In the neutral condition participants watched film clips that were effective in maintaining a calm (pleasant-neutral) state 1, such as a scene from the TV show “Rail Away.” The neutral video induced a pleasant, relaxing form of neutrality rather than plain neutrality, which could result in negative feeling states such as boredom or irritation [45]. The movie clips used in this study have been made publicly available on: https://osf.io/v9uf5/ (see Materials > Film clips.zip).

4.1.2. Manipulation Check

First, participants dichotomously rated their experience of 20 different emotional states (yes/no). If they answered “no,” participants were scored a 0 on that state. If they answered “yes,” participants additionally rated their experience of the respective state on a separate 5-point scale ranging from 1 (a little bit) to 5 (extreme). The yes/no queries were used to discourage participants from treating a unipolar measure of a state as a bipolar measure of one emotion (e.g., sadness) and its polar opposite (e.g., happiness) [61]. The questionnaire included all 16 items from Russell’s [46] circumplex model (e.g., happy, sad) and four additional discrete emotion terms, namely surprise, anger, fear, and disgust. Given the happy, fear, and neutral emotion induction condition, the discrete emotions happiness and fear, but also calmness, were marked as “target states” and these were analyzed separately. Furthermore, based on the circumplex model of emotion [46], different arrangements of feeling states were averaged to create indicators of high arousal (alert, excited, elated, happy, tensed, nervous, stressed, upset), low arousal (sad, depressed, lethargic, fatigued, calm, relaxed, serene, contented), positive valence (alert, excited, elated, happy, contented, serene, relaxed, calm), and negative valence (sad, depressed, lethargic, fatigued, upset, stressed, nervous, tense).
Arousal and valence were also measured directly via the affect grid [47] —a 7 × 7 matrix, where valence (negative-positive) was represented horizontally, and arousal (high-low) vertically. Participants had to place an “X” in one of the squares, yielding two scores (valence, arousal) ranging from 1 to 7.

4.1.3. Sweat Donation and Storage

Participants wore pre-cleaned t-shirts with absorbent pads attached via stainless steel poppets, 4 for each armpit. Following each session, pads from both the left and right underarm (n = 144 in total: 24 × 3 conditions × 2 axillae) were removed from the t-shirts using tweezers and scissors precleaned in a 96% ethanol solution. The pads were placed in labelled bags and vacuum-sealed, and bags were thereupon immediately stored in a freezer at −28 °C (Utrecht University, Utrecht, the Netherlands). After completion of the collection the samples were sent on dry ice by DHL same day delivery to URDV where they were stored at −80 °C.

4.2. Sweat Analysis Procedure with GC×GC-TOFMS

4.2.1. Instrumentation

A Pegasus® 4D GC×GC-TOFMS system (LECO) equipped with a liquid nitrogen cryogenic quad jet modulator and secondary oven was used for the sample analysis. The thermal desorption (TD) was performed with a thermal desorption unit (TDU) equipped with a MPS 2 auto-sampler and a liquid nitrogen (LN2) cooled CIS 4 programmed temperature vaporization (PTV) inlet (all from Gerstel) installed on an Agilent 6890 gas chromatograph (for further reference to methods used see Harker et al. [62]).

4.2.2. Sample Preparation/Extraction

Pads were analyzed using the HeadSpace Sorptive Extraction (HSSE) technique. For the extraction, 250 mL Scott-Duran GLS80 bottles were used, with screw caps with PFTE septum. A small magnet was attached to the outside wall of the bottle, just below the neck of the bottle, using adhesive tape. Then the pad was transferred into the bottle using a clean pair of tweezers and 5 µL internal standard solution was added using a 5 µL transfer pipettor at the bottom of the bottle. Finally, a PDMS stir bar (1 cm × 1 mm; Gerstel) was attached with forceps to the inside bottle wall (attached to the magnet) and the bottle was closed with the screw cap with PFTE septum. Extraction was done in an oven at 60 °C for 2 h. After the extraction the stir bar was removed with forceps and placed in a clean glass thermal desorption liner. The glass liner was placed in the Twister tray of the MPS sampler, ready for analysis.
Reconditioning of stir bars was done after use by soaking in ultrapure water and methanol 1:1 (v/v) for 1–2 h each; stir bars were then removed from the solvent and dried on a clean surface at room temperature for 1 h. Finally, the stir bars were thermally conditioned in the glass thermal desorption liners for 30 min at 250 °C in a flow of helium using the Gerstel tube conditioner (TC) and stored in the MPS-2 Twister tray. Used bottles and caps were rinsed with water and methanol and finally thermally cleaned in an oven at 75 °C for approximately 8 h.
The internal standard solution contained 30.4 µg/mL five times deuterium labelled benzaldehyde (Aldrich) using methanol as the solvent. The peak area of the internal standard was only used to monitor the efficiency of the extraction and the stability of the instrument.
Per day/sequence all pads from two participants were measured. The participants were randomly selected. All six pads from one participant were measured after each other and the left and right axilla of one treatment were analyzed together (L&R, R&L, L&R). Finally the order of the treatment was varied (H, N,F & N,H,F or F,N,H & H, F,N). Each sequence was started with a blank tube and at the end of the sequence a defrost cycle of the modulator was applied.

4.2.3. Thermal Desorption (TD) – GC×GC TOF MS

The stir bars were thermally desorbed in splitless mode by programming the TDU from 35 °C (held for 1 min) to 250 °C (held for 5 min) at 720 °C min−1 with 75 mL min−1 desorption flow. Desorbed compounds were focused at −120 °C on a liner filled with deactivated quartzwool in the LN2 cooled PTV inlet for subsequent GC×GC–TOF MS analysis. After desorption, the PTV inlet was programmed from −120 °C to 250 °C (held for 5 min) at 12 °C/sec to inject trapped compounds onto the analytical column. The injection was performed in split mode with a split ratio of 1:50. The first dimension (1D) column was a 30 m × 0.25 mm, df = 0.5 µm non-polar Rxi-5MS (Restek Corporation) and the second dimension (2D) column was 2 m × 0.15 mm, df = 0.15 µm mid-polar RTX-200 (Restek Corporation). Deactivated universal presstight connectors (Restek Corporation) were used for connecting the 1D and 2D columns. Helium was used as the carrier gas at a constant flow of 1.5 mL/min. The 1D GC oven temperature was initially set at 40 °C and held for 2 min, before being increased at 4 °C/min to 250 °C where it was held for 15 min. The temperature offset for the secondary oven was +10 °C relative to the 1D GC oven. The modulation period (PM) was 3 s, with a hot pulse time of 0.9 s and a cold pulse time of 0.60 s on each jet. The modulator temperature offset was +15 °C relative to the secondary oven. The transfer line was held at 275 °C. The TOFMS was operated in electron ionization (EI) mode at 70 eV, with an acquisition mass range of 35–300 amu, an acquisition rate of 150 Hz, and a detector voltage of 1800 V. The ion source was heated to 200 °C. Mass calibration and tuning was performed using perfluorotributylamine (PFTBA). An acquisition delay of 3 min was applied.

4.2.4. Data Processing

Data processing was performed using ChromaTOF® (version 4.50; LECO, St Joseph, MI, USA). This software was used for peak finding, with each peak defined as being a unique feature corresponding to a unique deconvoluted combination of retention times and mass spectrum, mass deconvolution, integration, and library searching. The baseline offset was set to 1, just above the noise level and automatic smoothing was applied. The first dimension peak width was set to 10 s, while the second dimension peak width was set to 0.1 s. A 65% mass spectral match was required to combine the subpeaks with each other and their corresponding base peak. A minimum signal-to-noise ratio (S/N) of 100 was used for the base peak detection with a minimum of 2 apexing masses and sub-peak detection was cut-off at a minimum S/N of 10. Traditional, not adaptive, integration was used. Peak identification was performed by a forward search to the 2011 National Institute of Standards and Technology (NIST) mass spectral library database with a minimum similarity match of >80%. The Statistical Compare software feature within ChromaTOF® was subsequently used to perform the peak alignment between samples using a mass spectral match threshold of 60%. Samples were input into the Statistical Compare and separated into different classes: Neutral (n = 48), Fear (n = 47, Happy (n = 47), Room (n = 24; denoted as Reference), and Clean Pads (n = 4). One pad of the Happy and Fear classes were missing, respectively. For the QC samples (clean pads spiked with deuterium labeled internal standards) a separate SC table was created. During the alignment, peak researching was performed with a minimum S/N cut-off of 50 in order to search for the peaks not found during the initial peak-finding step. A maximum retention time difference of 3 s (i.e., 1 modulation period) in the 1D and 0.05 s in the 2D was permitted during alignment to allow for retention time deviations between samples. During alignment, analytes were only retained if found in at least three samples OR if found in 50% of the samples within a class. The resulting peak table (3796 peaks) was exported and used for further multivariate analysis. Figure A3 (Appendix A) shows an example of a typical comprehensive GC×GC TOF MS contour plot of a test person.
The statistical analysis was performed in SIMCA (version 15.02.5959, Sartorius Stedim Biotech, Umëa, Sweden). Initially, principal component analysis (PCA) was performed on all 3796 peaks and pads (141 condition samples +24 reference samples + 4 clean samples = 169) to distinguish between reference and clean pads from the sweat pads from the subjects. Samples from left and right underarms of the same individual were treated as different samples. A visual inspection of PLS-DA score plots on all left panels of Figure 2 by donor ID and arm (L vs. R) distribution showed that L and R samples from the same donor were not identical but also that variation between underarms was smaller than between donors. The peak areas were scaled to Pareto prior to analysis [63]. This scaling method attenuates large peaks and enhances small peaks, and compromises between mean centering and unit variance scaling. This method was chosen to focus on the main effects. The scores plot (Figure A4; Appendix A) shows a clear separation between the control pads and the pads from the subjects, except for the one reference pad that was contaminated with sweat from a subject. Various (tentatively assigned) volatiles were detected in the control samples including for example acetic acid butyl ester, 1,1,3-trimethylcyclohexane, 1,1-dimethylehtyl-2-methyl-propionic acid, tricane, diethyltridecane, dodecane heptanone, hexadecane, 5,5-diethyltridecane, and ethylacetate. PCA analysis excluding the reference and clean samples identified one pad in the Happy group as outlier.
For the comparison between the condition classes the data were further cleaned up. For this, only peak intensities above an (arbitrary) threshold of 2400 were considered to focus on higher abundant peaks (23% of the dataset). Moreover, we excluded peaks (n = 2139) for which the median concentrations calculated over all classes were zero (defined as <2400) to focus on peaks that are representative for the conditions. Then, all zeros were removed, to prevent the interference of zero values, because they were not true values but rather indicated an absence of detection, small signals, or noise. In addition, the condition-matching pads of the missing pads (n = 6) were excluded to have for each subject and each arm the same number of condition pads. In total, 135 pads including 45 pads per condition class and 1655 peaks were taken to further analysis.
To identify the peaks that discriminate between the condition classes, partial least squares discriminant analysis (PLS-DA) was performed. This supervised method allows for finding the maximum differences between the conditions (rather than the maximum variance as in case of PCA). In total, three PLS-DA models were calculated, namely Neutral versus Happy, Neutral versus Fear, Happy versus Fear. The data were not normalized and Pareto scaling was applied prior to analysis. All models were validated using seven-fold cross-validation and permutation testing with 900 models of randomly permuted data [64]. Q2-values were calculated to assess the predictive performance of the models. A Q2 value of 1 indicates maximum predictive power, whereas Q2 values close to or below 0 indicate a lack of predictive power.
To visualize the final variances of all three conditions, PCA was performed on 94 peaks that were in total found to be discriminative in the three PLS-DA models. The data were not normalized and Pareto scaling was applied prior to analysis.
Univariate statistical analysis was performed in JMP® Pro 12.0.1. (SAS Institute Inc. 2013, 198 Cary, NC, USA). Receiver operating characteristic (ROC) curves were plotted to additionally confirm the discriminative power of each individual discriminative peak identified from the PLS-DA models Neutral-Fear, Neutral-Happy, Fear-Happy. The ROC curve displays the true positive rate (TPR; probability of correct classification, sensitivity) versus the false positive rate (FPR; probability of wrong classification, 1-specificity) based on the ranked peak areas against the binary classifier system (condition I e.g., Neutral versus condition II e.g., Fear). The condition class with the higher concentrations was defined as the positive level. Peak areas from condition I are perfectly discriminated from condition II, if they are systematically higher (or lower) in condition I than in condition II. In this case, ROC curve passes through the upper left corner and the area under the ROC curve (AUROC) is maximal. The steepness of the ROC curve indicates the relation between the true positive rate and the false positive rate. A peak is maximally discriminative when the true positive rate is maximum and the false positive rate is minimum.
Box plots were made on the discriminative peaks from the PLS-DA models and on the sum of all peak areas to visualize differences between the three condition groups. The box plot displays the maximum, 75 quartile, median, 25 quartile, and minimum values. In addition, the peak areas of the individual peaks were plotted against the total intensities to assess whether condition-dependent changes of the peaks are related to the total sweat production.

4.3. Ethics Statement

The study protocol was approved by the Wageningen University ethics committee in 2014 (U-IEC dossier 14/23; approval date 30 September 2014). All methods were performed in accordance with national and international guidelines and regulations, and with the Helsinki Declaration of 1975, as revised in 2013 [65]. Written informed consent was obtained from the sweat donors prior to starting the study. Sweat donors received a financial remuneration of €50, for participating in the study.

Supplementary Materials

The following are available online at: Movie clips used in this study: https://osf.io/v9uf5/ (see Materials > Film clips.zip); raw data etc., will be made available upon publication. The following are available online at https://www.mdpi.com/2218-1989/10/3/84/s1. Figure S1: Individual box plots for each peak across the emotion condition fear, happy, and emotionally neural, S2 A+B: Additional testing to explore bipolar pattern in happy chemical cluster: results + statistical analyses, S3: Study to address the social factor during happiness induction, Figure S3: Scatter plot of self-reported emotions in happiness condition (top panel) and neutral (bottom panel) on PC1 and PC2 comparing scores from four studies (n = 84) in which multiple men sat together and happiness condition was always second (pink) vs. a study (n = 24) in which men were sitting alone and condition was counterbalanced (with happiness condition being first or second; turquoise).

Author Contributions

Conceptualization, M.A.M.S., J.H.K., C.B., J.d.G. and G.R.S.; data curation, E.A.E.R., D.M.J., and J.H.K.; formal analysis, E.A.E.R., D.M.J., E.v.V., J.H.K., M.B.E., B.M., and J.d.G.; investigation, E.A.E.R., D.M.J., M.B.E., J.d.G., and G.R.S.; methodology, M.A.M.S., E.A.E.R., D.M.J., E.v.V., J.H.K., M.B.E., B.M., J.d.G., and G.R.S.; project administration, C.B. and I.G.; supervision, M.A.M.S., C.B., and I.G.; validation, E.A.E.R., D.M.J., J.H.K., B.M., and J.d.G.; visualization, E.A.E.R., D.M.J., J.H.K., and J.d.G.; writing—original draft, M.A.M.S., E.A.E.R., D.M.J., J.H.K., B.M., J.d.G., and G.R.S.; writing—review and editing, E.A.E.R., D.M.J., J.H.K., C.B., I.G., B.M., J.d.G., and G.R.S. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by Unilever R&D, Vlaardingen.

Acknowledgments

The authors would like to express their thanks to Joy Wilkinson for supporting and initiating this project; Patricia Bulsing and Vanessa Martens for management of the project; Sonja Kaal for her help with statistical analyses; Svetlana Riazanskaia, David Smith, Daniel C’Ailceta for feedback and advice on chemical analytical work; Anna Thomas and Timo Giesbrecht for their helpful comments and support throughout the project.

Conflicts of Interest

This study was funded by Unilever. Unilever is Fast-Moving Consumer Goods company that sells Beauty and Personal Care products. All authors declare no potential (financial nor non-financial) conflicts of interest. Unilever authors were involved in chemical analyses of human sweat and relevant statistical analyses. The design of the study, sweat collection, and statistical analyses related to human self-report data was led by the authors from Utrecht University.

Appendix A. Additional Data Figures and Tables

Figure A1. Receiver operating characteristic (ROC) curves of peak 2209 (left panel) and 234 (right panel) for the classifications Happy–Neutral (blue), Fear–Neutral (red), and Fear–Happy (green). For peak 2209 Neutral was the positive level for Happy–Neutral and Fear–Neutral, and Happy was the positive level for Fear–Happy. For 234 fear was the positive level for all three classifications.
Figure A1. Receiver operating characteristic (ROC) curves of peak 2209 (left panel) and 234 (right panel) for the classifications Happy–Neutral (blue), Fear–Neutral (red), and Fear–Happy (green). For peak 2209 Neutral was the positive level for Happy–Neutral and Fear–Neutral, and Happy was the positive level for Fear–Happy. For 234 fear was the positive level for all three classifications.
Metabolites 10 00084 g0a1
Figure A2. Box plot of three condition classes based on the sum of 1655 peak areas per sweat sample.
Figure A2. Box plot of three condition classes based on the sum of 1655 peak areas per sweat sample.
Metabolites 10 00084 g0a2
Figure A3. Example of a typical comprehensive GC×GC TOF MS contour plot of a test person, with x-axis denoting first dimension rt time in min (on non-polar column in Table A1), and y-axis denoting second dimension rt time is seconds (mid polar column in Table A1).
Figure A3. Example of a typical comprehensive GC×GC TOF MS contour plot of a test person, with x-axis denoting first dimension rt time in min (on non-polar column in Table A1), and y-axis denoting second dimension rt time is seconds (mid polar column in Table A1).
Metabolites 10 00084 g0a3
Figure A4. PCA performed on 3796 RT/MS peaks and 169 pads using Pareto scaling. The scores plot (upper panel) displays the variation of the pads classified as clean (from clean pads), reference (from pads deposited in the room), and the condition samples (Neutral, Happy, Fear from the subjects). The loading plot (lower panel) displays the variation of the analytes with tentative assignments.
Figure A4. PCA performed on 3796 RT/MS peaks and 169 pads using Pareto scaling. The scores plot (upper panel) displays the variation of the pads classified as clean (from clean pads), reference (from pads deposited in the room), and the condition samples (Neutral, Happy, Fear from the subjects). The loading plot (lower panel) displays the variation of the analytes with tentative assignments.
Metabolites 10 00084 g0a4
Table A1. Discriminative peaks from three PLS-Da models with peak in order of first dimension retention time on non-polar column (RT Dim 1) with addition of second dimension retention time on mid-polar column (RT Dim 2).
Table A1. Discriminative peaks from three PLS-Da models with peak in order of first dimension retention time on non-polar column (RT Dim 1) with addition of second dimension retention time on mid-polar column (RT Dim 2).
PeakRT Dim 1RT Dim 2MassDiscriminationUp or Down
37213.0651.5037343F H vs ND
61252.9861.5444887F H vs ND
75266.8961.6986684F H vs ND
82278.8811.5221687F H vs ND
94290.8671.549787F H vs ND
122311.8421.7020761F H vs ND
146347.7982.0713369F H vs NU
161368.0591.6856959F H vs ND
177381.9862.0108559F H vs ND
199410.8521.8462445F H vs NU
211422.8492.5979243F H vs ND
221431.6982.1741655F H vs NU
234446.541.9437456F H vs NU
235445.6112.111545F H vs ND
251470.6511.8519543F H vs ND
299536.5541.8113545F H vs ND
456713.8311.9689861F H vs ND
505773.2872.3509883F H vs NU
531794.2622.8830396F H vs ND
560824.2261.7422857F H vs NU
562823.5321.91573103H vs F NU
638878.1612.0662257F H vs NU
662896.752.418681F H vs NU
747967.7542.2181146F H vs NU
760979.7952.229891H vs F NU
7831001.010.42477985F H vs NU
7941007.012.3344870F H vs NU
79910102.3005699F H vs NU
8411038.332.2012959F H vs NU
9161084.91.5754293H vs F ND
9201084.971.6566881F H vs NU
9571111.882.0459257F H vs NU
9651118.61.60533267F H vs ND
9791123.872.3847981F H vs NU
10211146.782.9112787F H vs ND
10371150.082.7332199F H vs NU
11341225.920.28506685F H vs NU
11661240.732.2519113F H vs NU
11851252.261.75714207H vs F NU
12061263.012.21667101F H vs ND
12961336.582.0116657F H vs NU
13131354.662.3280281F H vs NU
13671411.521.8379859F H vs NU
13871422.391.8460459F H vs NU
14051433.160.157468113H vs F NU
14371458.511.54951163F H vs NU
14491462.461.976758F H vs NU
14501462.462.1945127F H vs NU
14571465.462.42029131F H vs NU
14971507.462.1505859H vs F NU
15211525.411.82177207H vs F NU
15531549.822.1710559F H vs ND
15621555.151.69778281H vs F NU
15871570.332.2686381F H vs NU
15991581.971.5498573H vs F NU
16161591.931.7407159F H vs NU
16561628.350.058274184F H vs NU
16721639.321.8459259F H vs NU
16731638.782.08221133H vs F NU
17351687.011.571671F vs H NU
17651711.211.7153969F vs H NU
17761717.151.8235659F H vs NU
17801723.151.93293114F H vs NU
18101741.732.2816259F H vs ND
18421762.11.84815109H vs F NU
18501767.882.3002887F H vs ND
18661777.232.7036287F H vs ND
19111804.521.5426257H vs F NU
20111848.741.9391895F H vs NU
20191848.071.7377897H vs F ND
20501867.031.6533492H vs F NU
20731875.982.19588177F H vs NU
20881884.551.736870F vs H NU
21071893.491.543970F vs H NU
21271902.931.9109772F H vs NU
21491913.331.9694158H vs F NU
22091944.881.68335191F H vs ND
22621974.911.7981859H vs F NU
22771980.841.7227557F H vs NU
23482013.62.3127287F H vs ND
23602022.792.48686100F H vs ND
24122046.831.84665113F H vs NU
24902091.721.88072114F H vs NU
25292112.931.62908119H vs F ND
25782139.681.87183150F vs H NU
26082156.881.7681170H vs F ND
26252165.991.5824391F H vs NU
27562220.551.8675758F H vs NU
31002391.312.21692219F H vs NU
31902457.321.60417151H vs F NU
32312493.191.8400569F vs H NU
32342496.291.60178151H vs F NU
34372739.031.59881151H vs F NU
34892830.621.62092149H vs F NU
Note: “Mass” corresponds to the unique ion in the deconvoluted EI spectrum determined by the LECO ChromaToF software. “Discrimination” indicates significance from either of Happy (H) and Fear (F), or both, over Neutral (N) associated with donor peaks being up (U) or down (D) for Happy/Fearful.

References

  1. De Groot, J.H.B.; Smeets, M.A.M.; Rowson, M.J.; Bulsing, P.J.; Blonk, C.G.; Wilkinson, J.E.; Semin, G.R. A sniff of happiness. Psychol. Sci. 2015, 26, 684–700. [Google Scholar] [CrossRef] [PubMed]
  2. Pause, B.M. Processing of body odor signals by the human brain. Chemosens. Percept. 2012, 5, 55–63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Semin, G.R.; de Groot, J.H.B. The chemical bases of human sociality. Trends Cogn. Sci. 2013, 17, 427–429. [Google Scholar] [CrossRef] [PubMed]
  4. Olsson, M.J.; Lundström, J.N.; Kimball, B.A.; Gordon, A.R.; Karshikoff, B.; Hosseini, N.; Sorjonen, K.; Olgart Höglund, C.; Solares, C.; Soop, A.; et al. The scent of disease. Psychol. Sci. 2014, 25, 817–823. [Google Scholar] [CrossRef] [PubMed]
  5. De Groot, J.H.B.; Smeets, M.A.M. Human fear chemosignaling: Evidence from a meta-analysis. Chem. Senses 2017, 42, 663–673. [Google Scholar] [CrossRef]
  6. Laidre, M.E.; Johnstone, R.A. Animal signals. Curr. Biol. 2013, 23, R829–R833. [Google Scholar] [CrossRef] [Green Version]
  7. Chen, D.; Haviland-Jones, J. Human olfactory communication of emotion. Percept. Mot. Skills 2000, 91, 771–781. [Google Scholar] [CrossRef]
  8. De Groot, J.H.B.; Semin, G.R.; Smeets, M.A.M. Chemical communication of fear: A case of male-female asymmetry. J. Exp. Psychol. Gen. 2014, 143, 1515–1525. [Google Scholar] [CrossRef]
  9. Mutic, S.; Parma, V.; Brünner, Y.F.; Freiherr, J. You smell dangerous: Communicating fight responses through human chemosignals of aggression. Chem. Senses 2016, 41, 35–43. [Google Scholar] [CrossRef] [Green Version]
  10. Zhou, W.; Chen, D. Fear-related chemosignals modulate recognition of fear in ambiguous facial expressions. Psychol. Sci. 2009, 20, 177–183. [Google Scholar] [CrossRef]
  11. Bensafi, M.; Brown, W.M.; Tsutsui, T.; Mainland, J.D.; Johnson, B.N.; Bremner, E.A.; Young, N.; Mauss, I.; Ray, B.; Gross, J.; et al. Sex-steroid derived compounds induce sex-specific effects on autonomic nervous system function in humans. Behav. Neurosci. 2003, 117, 1125–1134. [Google Scholar] [CrossRef] [PubMed]
  12. Parma, V.; Tirindelli, R.; Bisazza, A.; Massaccesi, S.; Castiello, U. Subliminally perceived odours modulate female intrasexual competition: An eye movement study. PLoS ONE 2012, 7, e30645. [Google Scholar] [CrossRef] [PubMed]
  13. Martins, Y.; Preti, G.; Crabtree, C.R.; Runyan, T.; Vainius, A.A.; Wysocki, C.J. Preference for human body odors is influenced by gender and sexual orientation. Psychol. Sci. 2005, 16, 694–701. [Google Scholar] [CrossRef] [PubMed]
  14. Lübke, K.T.; Hoenen, M.; Pause, B.M. Differential processing of social chemosignals obtained from potential partners in regards to gender and sexual orientation. Behav. Brain Res. 2012, 228, 375–387. [Google Scholar] [CrossRef] [PubMed]
  15. Gelstein, S.; Yeshurun, Y.; Rozenkrantz, L.; Shushan, S.; Frumin, I.; Roth, Y.; Sobel, N. Human tears contain a chemosignal. Science 2011, 331, 226–231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Porter, R.H.; Schaal, B. Olfaction and development of social preferences in neonatal organisms. In Handbook of Olfaction and Gustation; Doty, R.L., Ed.; Marcel Dekker: New York, NY, USA, 1995; pp. 299–321. [Google Scholar]
  17. Havlicek, J.; Roberts, S.C.; Flegr, J. Women’s preference for dominant male odour: Effects of menstrual cycle and relationship status. Biol. Lett. 2005, 1, 256–259. [Google Scholar] [CrossRef] [Green Version]
  18. McClintock, M.K. Menstrual Synchrony and Suppression. Nature 1971, 229, 244–245. [Google Scholar] [CrossRef]
  19. Zhengwei, Y.; Schank, J.C. Women do not synchronize their menstrual cycles. Hum. Nat. 2006, 17, 433–447. [Google Scholar]
  20. Ziomkiewicz, A. Menstrual synchrony: Fact or artifact? Hum. Nat. 2006, 17, 419–432. [Google Scholar] [CrossRef]
  21. Sorokowska, A.; Sorokowski, P.; Szmajke, A. Does personality smell? Accuracy of personality assessments based on body odour. Eur. J. Pers. 2012, 26, 496–503. [Google Scholar] [CrossRef]
  22. Wyatt, T.D. The search for human pheromones: The lost decades and the necessity of returning to first principles. Proc. R. Soc. B Biol. Sci. 2015, 282, 20142994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Doty, R.L. The Great Pheromone Myth; Johns Hopkins University Press: Baltimore, MD, USA, 2010. [Google Scholar]
  24. Chen, D.; Katdare, A.; Lucas, N. Chemosignals of fear enhance cognitive performance in humans. Chem. Senses 2006, 31, 415–423. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. De Groot, J.H.B.; Smeets, M.A.M.; Kaldewaij, A.; Duijndam, M.J.A.; Semin, G.R. Chemosignals communicate human emotions. Psychol. Sci. 2012, 23, 1417–1424. [Google Scholar] [CrossRef] [PubMed]
  26. Kamiloglu, R.G.; Smeets, M.A.M.; de Groot, J.H.B.; Semin, G.R. Fear odor facilitates the detection of fear expressions over other negative expressions. Chem. Senses 2018, 43, 419–426. [Google Scholar] [CrossRef] [PubMed]
  27. De Groot, J.H.B.; Semin, G.R.; Smeets, M.A.M. I can see, hear, and smell your fear: Comparing olfactory and audiovisual media in fear communication. J. Exp. Psychol. Gen. 2014, 143, 825–834. [Google Scholar] [CrossRef] [PubMed]
  28. De Groot, J.H.B.; van Houtum, L.A.E.M.; Gortemaker, I.; Ye, Y.; Chen, W.; Zhou, W.; Smeets, M.A.M. Beyond the west: Chemosignaling of emotions transcends ethno-cultural boundaries. Psychoneuroendocrinology 2018, 98, 177–185. [Google Scholar] [CrossRef] [PubMed]
  29. Dalton, P.; Mauté, C.; Jaén, C.; Wilson, T. Chemosignals of stress influence social judgments. PLoS ONE 2013, 8, e77144. [Google Scholar] [CrossRef] [Green Version]
  30. Singh, P.B.; Young, A.; Lind, S.; Leegaard, M.C.; Capuozzo, A.; Parma, V. Smelling anxiety chemosignals impairs clinical performance of dental students. Chem. Senses 2018, 43, 411–417. [Google Scholar] [CrossRef] [Green Version]
  31. Adolph, D.; Meister, L.; Pause, B.M. Context counts! Social anxiety modulates the processing of fearful faces in the context of chemosensory anxiety signals. Front. Hum. Neurosci. 2013, 7, 1–17. [Google Scholar] [CrossRef] [Green Version]
  32. Pause, B.M.; Lübke, K.; Laudien, J.H.; Ferstl, R. Intensified neuronal investment in the processing of chemosensory anxiety signals in non-socially anxious and socially anxious individuals. PLoS ONE 2010, 5, e10342. [Google Scholar] [CrossRef] [Green Version]
  33. Mujica-Parodi, L.R.; Strey, H.H.; Frederick, B.; Savoy, R.; Cox, D.; Botanov, Y.; Tolkunov, D.; Rubin, D.; Weber, J. Chemosensory cues to conspecific emotional stress activate amygdala in humans. PLoS ONE 2009, 4, e6415. [Google Scholar] [CrossRef] [PubMed]
  34. Maier, A.; Scheele, D.; Spengler, F.B.; Menba, T.; Mohr, F.; Güntürkün, O.; Stoffel-Wagner, B.; Kinfe, T.M.; Maier, W.; Khalsa, S.S.; et al. Oxytocin reduces a chemosensory-induced stress bias in social perception. Neuropsychopharmacol. 2019, 44, 281–288. [Google Scholar] [CrossRef] [PubMed]
  35. Prehn-Kristensen, A.; Wiesner, C.; Bergmann, T.O.; Wolff, S.; Jansen, O.; Mehdorn, H.M.; Ferstl, R.; Pause, B.M. Induction of empathy by the smell of anxiety. PLoS ONE 2009, 4, e5987. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. De Groot, J.H.B.; Semin, G.R.; Smeets, M.A.M. On the communicative function of body odors: A theoretical integration and review. Perspect. Psychol. Sci. 2017, 12, 306–324. [Google Scholar] [CrossRef] [PubMed]
  37. Goyert, H.F.; Frank, M.E.; Gent, J.F.; Hettinger, T.P. Characteristic component odors emerge from mixtures after selective adaptation. Brain Res. Bull. 2007, 72, 1–9. [Google Scholar] [CrossRef] [Green Version]
  38. Wysocki, C.J.; Preti, G. Facts, fallacies, fears, and frustrations with human pheromones. Anat. Rec. Part A Discov. Mol. Cell. Evol. Biol. 2004, 281, 1201–1211. [Google Scholar] [CrossRef]
  39. Dormont, L.; Bessière, J.M.; Cohuet, A. Human skin volatiles: A review. J. Chem. Ecol. 2013, 39, 569–578. [Google Scholar] [CrossRef]
  40. De Lacy Costello, B.; Amann, A.; Al-Kateb, H.; Flynn, C.; Filipiak, W.; Khalid, T.; Osborne, D.; Ratcliffe, N.M. A review of the volatiles from the healthy human body. J. Breath Res. 2014, 8, 014001. [Google Scholar] [CrossRef]
  41. Fuller, R.G.; Sheehy-Skeffington, A. Effects of group laughter on responses to humourous material, a replication and extension. Psychol. Rep. 1974, 35, 531–534. [Google Scholar] [CrossRef]
  42. Devereux, P.G.; Ginsburg, G.P. Sociality effects on the production of laughter. J. Gen. Psychol. 2001, 128, 227–240. [Google Scholar] [CrossRef]
  43. Provine, R.R. Contagious laughter: Laughter is a sufficient stimulus for laughs and smiles. Bull. Psychon. Soc. 1992, 30, 1–4. [Google Scholar] [CrossRef]
  44. Gervais, M.; Wilson, D.S. The evolution and functions of laughter and humor: A synthetic approach. Q. Rev. Biol. 2019, 80, 395–430. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Rottenberg, J.; Ray, R.D.; Gross, J.J. Emotion elicitation using films. In Handbook of Emotion Elicitation and Assessment; Coan, J.A., Allen, J.J.B., Eds.; Oxford University Press: Oxford, UK, 2007; pp. 9–28. [Google Scholar]
  46. Russell, J.A. A circumplex model of affect. J. Pers. Soc. Psychol. 1980, 39, 1161–1178. [Google Scholar] [CrossRef]
  47. Russell, J.A.; Weiss, A.; Mendelsohn, G.A. Affect grid: A single-item scale of pleasure and arousal. J. Pers. Soc. Psychol. 1989, 57, 493–502. [Google Scholar] [CrossRef]
  48. Williams, J.; Stönner, C.; Wicker, J.; Krauter, N.; Derstroff, B.; Bourtsoukidis, E.; Klüpfel, T.; Kramer, S. Cinema audiences reproducibly vary the chemical composition of air during films, by broadcasting scene specific emissions on breath. Sci. Rep. 2016, 6, 25464. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  49. Preti, G.; Dalton, P.; Maute, C. Analytical identification of stress odors from human breath. HDIAC J. 2019, 6, 10–15. [Google Scholar]
  50. Zeng, X.; Leyden, J.J.; Lawley, H.J.; Sawano, K.; Nohara, I.; Preti, G. Analysis of characteristic odors from human male axillae. J. Chem. Ecol. 1991, 17, 1469–1492. [Google Scholar] [CrossRef]
  51. Pandey, S.K.; Kim, K.H. Human body-odor components and their determination. Trends Anal. Chem. 2011, 30, 784–796. [Google Scholar] [CrossRef]
  52. Shimomura, Y.; Kitaura, Y. Physiological and pathological roles of branched-chain amino acids in the regulation of protein and energy metabolism and neurological functions. Pharmacol. Res. 2018, 133, 215–217. [Google Scholar] [CrossRef]
  53. Hagenfeldt, L.; Wahren, J. Human Forearm Muscle Metabolism during Exercise III Uptake, release and oxidation of β-hydroxybutyrate and observations on the β-hydroxybutyrate/acetoacetate ratio. Scand. J. Clin. Lab. Invest. 1968, 21, 314–320. [Google Scholar] [CrossRef]
  54. De Groot, J.H.B.; Smeets, M.A.M.; Semin, G.R. Rapid stress drives chemical transfer of fear from sender to receiver. PLoS ONE 2015, 10, e0118211. [Google Scholar] [CrossRef] [PubMed]
  55. De Groot, J.H.B.; Kirk, P.A.; Gottfried, J.A. Encoding fear intensity in human sweat. Philos. Trans. R. Soc. B Biol. Sci. 2020. [Google Scholar] [CrossRef]
  56. De Groot, J.H.B.; Kirk, P.A.; Gottfried, J.A. Titrating the smell of fear: Dose-invariant responses of human chemosignals on behavior, physiology, and the brain. Prep.
  57. Lenochova, P.; Roberts, S.C.; Havlicek, J. Methods of human body odor sampling: The effect of freezing. Chem. Senses 2009, 34, 127–138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Stevenson, R.J. An initial evaluation of the functions of human olfaction. Chem. Senses 2010, 35, 3–20. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. McGann, J.P. Poor human olfaction is a 19th-century myth. Science 2017, 356, eaam7263. [Google Scholar] [CrossRef] [Green Version]
  60. Schaefer, A.; Nils, F.; Sanchez, X.; Philippot, P. Assessing the effectiveness of a large database of emotion-eliciting films: A new tool for emotion researchers. Cogn. Emot. 2010, 24, 1153–1172. [Google Scholar] [CrossRef]
  61. Larsen, J.T.; McGraw, A.P.; Cacioppo, J.T. Can people feel happy and sad at the same time? J. Pers. Soc. Psychol. 2001, 81, 684–696. [Google Scholar] [CrossRef]
  62. Harker, M.; Carvell, A.-M.; Marti, V.P.J.; Riazanskaia, S.; Kelso, H.; Taylor, D.; Grimshaw, S.; Arnold, D.S.; Zillmer, R.; Shaw, J.; et al. Functional characterisation of a SNP in the ABCC11 allele—Effects on axillary skin metabolism, odour generation and associated behaviours. J. Dermatol. Sci. 2014, 73, 23–30. [Google Scholar] [CrossRef]
  63. Van den Berg, R.A.; Hoefsloot, H.C.; Westerhuis, J.A.; Smilde, A.K.; van der Werf, M.J. Centering, scaling, and transformations: Improving the biological information content of metabolomics data. BMC Genom. 2006, 7, 142. [Google Scholar] [CrossRef] [Green Version]
  64. Van Velzen, E.J.J.; Westerhuis, J.A.; van Duynhoven, J.P.M.; van Dorsten, F.A.; Hoefsloot, H.C.J.; Jacobs, D.M.; Smit, S.; Draijer, R.; Kroner, C.I.; Smilde, A.K. Multilevel data analysis of a crossover designed human nutritional intervention study. J. Proteome Res. 2008, 7, 4483–4491. [Google Scholar] [CrossRef] [PubMed]
  65. Association, W.M. World Medical Association Declaration of Helsinki: Ethical Principles for Medical Research Involving Human Subjects. JAMA 2013, 310, 2191–2194. [Google Scholar]
Figure 1. Effects of emotion induction (fear, neutral, happy) on self-reported emotional state: (A) Spider plot showing participants’ emotion profiles on affective circumplex [46]; (B) experienced core affect (valence, arousal) on affect grid [47]; (C) experienced discrete emotions. Error bars: ± 1 standard error of the mean (SEM).
Figure 1. Effects of emotion induction (fear, neutral, happy) on self-reported emotional state: (A) Spider plot showing participants’ emotion profiles on affective circumplex [46]; (B) experienced core affect (valence, arousal) on affect grid [47]; (C) experienced discrete emotions. Error bars: ± 1 standard error of the mean (SEM).
Metabolites 10 00084 g001
Figure 2. Partial least squares discriminant analysis (PLS-DA) score plots (left panels) and loading plots (right panels) based on 1655 peak areas and on two condition classes at a time: Neutral versus Fear (upper panels), Neutral versus Happy (middle panels), Fear versus Happy (lower panel). The first two principal components (PCs) t[1] and t[2] are displayed. The loading plots show weighted coefficients of the principal components. The data were not normalized and Pareto-scaled prior to analysis. Q2 values > 0.2 were statistically significant.
Figure 2. Partial least squares discriminant analysis (PLS-DA) score plots (left panels) and loading plots (right panels) based on 1655 peak areas and on two condition classes at a time: Neutral versus Fear (upper panels), Neutral versus Happy (middle panels), Fear versus Happy (lower panel). The first two principal components (PCs) t[1] and t[2] are displayed. The loading plots show weighted coefficients of the principal components. The data were not normalized and Pareto-scaled prior to analysis. Q2 values > 0.2 were statistically significant.
Metabolites 10 00084 g002
Figure 3. Principal component analysis (PCA) score plot based on 135 sweat samples and 94 peak areas that were significantly discriminating the condition groups in the PLS-DA models. The principal components t[1], t[2] are displayed, explaining R2X[1]= 0.437 and R2X[2] = 0.141 of the variances, respectively. The data were not normalized and Pareto scaling was applied.
Figure 3. Principal component analysis (PCA) score plot based on 135 sweat samples and 94 peak areas that were significantly discriminating the condition groups in the PLS-DA models. The principal components t[1], t[2] are displayed, explaining R2X[1]= 0.437 and R2X[2] = 0.141 of the variances, respectively. The data were not normalized and Pareto scaling was applied.
Metabolites 10 00084 g003
Figure 4. Box plots showing individual peaks per condition class: Panels (AE) show the boxplots of the peaks 146, 1850, 2578, 2019, and 1842, respectively, per condition for each of the 24 panelists as observed in underarm sweat both left and right (so 48 datapoints for each of the three conditions in each graph).
Figure 4. Box plots showing individual peaks per condition class: Panels (AE) show the boxplots of the peaks 146, 1850, 2578, 2019, and 1842, respectively, per condition for each of the 24 panelists as observed in underarm sweat both left and right (so 48 datapoints for each of the three conditions in each graph).
Metabolites 10 00084 g004
Table 1. Results of emotion manipulation on 20 discrete emotion states, three of which are the target states, rated from 0 to 5, and separately on valence, and arousal (affect grid).
Table 1. Results of emotion manipulation on 20 discrete emotion states, three of which are the target states, rated from 0 to 5, and separately on valence, and arousal (affect grid).
FeelingFriedman TestComparison Condition (Intended Emotions)Order
Fear vs. HappyFear vs. NeutralHappy vs. Neutral
Target State
Fearc2 = 38.00, p < 0.001Z = 3.9, p < 0.001Z = 3.9, p < 0.001Z = 0.0, p > 0.999F>H=N
Happinessc2 = 40.77, p < 0.001Z = −4.5, p < 0.001Z = −2.4, p = 0.017Z = 4.2, p < 0.001H>N>F
Calmnessc2 = 33.11, p < 0.001Z = −3.9, p < 0.001Z = −4.0, p < 0.001Z = −1.7, p = 0.096N=H>F
Affect Grid
Valencec2 = 34.63, p < 0.001Z = −4.2, p < 0.001Z = −3.2, p< 0.001Z = 3.8, p < 0.001H>N>F
Arousalc2 = 31.73, p < 0.001Z = 2.6, p < 0.001Z = 4.3, p< 0.001Z = 3.2, p < 0.001F>H>N
Other States
Tensec2 = 41.38, p < 0.001Z = 4.1, p < 0.001Z = 4.1, p < 0.001Z = −1.0, p = 0.317F>N=H
Nervousc2 = 40.00, p < 0.001Z = 4.0, p < 0.001Z = 4.0, p < 0.001Z = 0.0, p > 0.999F>H=N
Relaxedc2 = 39.74, p < 0.001Z = −4.3, p < 0.001Z = −4.2, p < 0.001Z = 1.5, p = 0.137H=N>F
Alertc2 = 39.00, p < 0.001Z = 4.1, p < 0.001Z = 4.2, p < 0.001Z = 1.8, p = 0.061F>H=N
Elatedc2 = 36.90, p < 0.001Z = −4.2, p < 0.001Z = −1.5, p = 0.141Z = 4.0, p < 0.001H>N=F
Serenec2 = 30.90, p < 0.001Z = −3.8, p < 0.001Z = −4.1, p < 0.001Z = −1.2, p = 0.213N=H>F
Contentc2 = 29.54, p < 0.001Z = −4.0, p < 0.001Z = −2.0, p = 0.050Z = 3.6, p < 0.001H>N>F
Disgustc2 = 29.39, p < 0.001Z = 3.4, p < 0.001Z = 3.4, p < 0.001Z = −1.0, p = 0.317F>N=H
Surprisedc2 = 27.49, p < 0.001Z = 3.3, p < 0.001Z = 3.7, p < 0.001Z = 2.6, p = 0.011F>H>N
Fatiguedc2 = 20.35, p < 0.001Z = 0.3, p = 0.752Z = −3.5, p < 0.001Z = −3.2, p < 0.001N > H=F
Stressedc2 = 17.43, p < 0.001Z = 2.7, p = 0.007Z = 2.7, p = 0.007Z = −1.0, p = 0.317F>N=H
Upsetc2 = 16.29, p < 0.001Z = 2.9, p = 0.004Z = 2.9, p = 0.004Z = −1.4, p = 0.157F>N=H
Lethargicc2 = 11.35, p = 0.003Z = −0.5, p = 0.593Z = −1.9, p = 0.062Z = −2.2, p = 0.030N>H=F
Excitedc2 = 11.15, p = 0.004Z = 0.0, p = 0.972Z = 3.0, p = 0.002Z = 2.9, p = 0.004F=H>N
Depressedc2 = 4.00, p = 0.135Z = 1.4, p = 0.157Z = 0.1, p = 0.916Z = 1.8, p = 0.068F=H=N
Sadc2 = 2.00, p = 0.368Z = 1.0, p = 0.317Z = −0.6, p = 0.564Z = −1.4, p = 0.157F=H=N
Angryc2 = 1.40, p = 0.497Z = 0.3, p = 0.785Z = 1.4, p = 0.157Z = 0.5, p = 0.655F=H=N
Note: Left column: results from non-parametric Friedman test, degrees of freedom: 2, n = 24. Middle columns: follow-up Wilcoxon signed-ranks tests for specific comparisons. Right column: expected ordering of conditions on a particular self-report variable, with signs (=, <, >) indicating relative effect direction. Manipulation checks in bold.
Table 2. Response pattern (high vs. low) of the two chemical compound classes by emotion condition (neutral, fear, and happy) with happy showing a bipolar pattern with individuals in subcluster HF overlapping with Fear pattern; in subcluster HN overlapping with Neutral. n = number of donors.
Table 2. Response pattern (high vs. low) of the two chemical compound classes by emotion condition (neutral, fear, and happy) with happy showing a bipolar pattern with individuals in subcluster HF overlapping with Fear pattern; in subcluster HN overlapping with Neutral. n = number of donors.
Chemical Class Condition
NeutralFearHappySubcluster
Esters and 5-rings High LowLow HF (n = 14)
HighHN (n = 8)
Linear aldehydes/ketones Low HighHigh HF (n = 14)
LowHN (n = 8)
Table 3. Significant compounds by class 1 and class 2 compounds using Tukey–Kramer, all pairs testing with p < 0.05.
Table 3. Significant compounds by class 1 and class 2 compounds using Tukey–Kramer, all pairs testing with p < 0.05.
Chemical Class
-CLASS 1CLASS 2
Linear Aldehydes / Ketones / AlcoholsEsters and Some Ring Structures *
1Peak 146: 3-penten-2-onePeak 37: ethyl acetate
2Peak 234: hexanalPeak 61: 4-methyl-1,3-dioxolane
3Peak 505: 2-heptenalPeak 122: propyl acetate
4Peak 560: 1-octen-3-olPeak 161: 1-ethoxy-2-propanol
5Peak 638: octanalPeak 251: butyl acetate
6Peak 662: 2,4-heptadienalPeak 531: 2-cyclopenten-1-one, 3-methyl
7Peak 794: 2-octenalPeak 2360: 1,6-dioxacyclododecane-7,12-dione
8(Peak 957: nonanal)
9Peak 920: 2-hexylfuran
10Peak 979: octadienal
11Peak 1134: γ-heptyl-lacton
12(Peak 1296: decanal)
13Peak 1313: nonadienal
14Peak 1449: 2-Me-decanal
15Peak 1587: decadienal
16Peak 2127: 3-tridecanone
Note: Compounds in parentheses are peaks that were identified with minor differences. * Not all ring structures are in this class.

Share and Cite

MDPI and ACS Style

Smeets, M.A.M.; Rosing, E.A.E.; Jacobs, D.M.; van Velzen, E.; Koek, J.H.; Blonk, C.; Gortemaker, I.; Eidhof, M.B.; Markovitch, B.; de Groot, J.; et al. Chemical Fingerprints of Emotional Body Odor. Metabolites 2020, 10, 84. https://doi.org/10.3390/metabo10030084

AMA Style

Smeets MAM, Rosing EAE, Jacobs DM, van Velzen E, Koek JH, Blonk C, Gortemaker I, Eidhof MB, Markovitch B, de Groot J, et al. Chemical Fingerprints of Emotional Body Odor. Metabolites. 2020; 10(3):84. https://doi.org/10.3390/metabo10030084

Chicago/Turabian Style

Smeets, Monique A.M., Egge A.E. Rosing, Doris M. Jacobs, Ewoud van Velzen, Jean H. Koek, Cor Blonk, Ilse Gortemaker, Marloes B. Eidhof, Benyamin Markovitch, Jasper de Groot, and et al. 2020. "Chemical Fingerprints of Emotional Body Odor" Metabolites 10, no. 3: 84. https://doi.org/10.3390/metabo10030084

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop