Estudio Feromonas Mujeres Perfume
Estudio Feromonas Mujeres Perfume
Correspondence
ktouhara@g.ecc.u-tokyo.ac.jp
In brief
Human physiology; Natural sciences;
Chemistry; Biological sciences
Highlights
• Axillary odor is more pleasant, fragrant, and less vinegary
during the ovulatory phase
iScience ll
OPEN ACCESS
Article
Human ovulatory phase-increasing
odors cause positive emotions
and stress-suppressive effects in males
Nozomi Ohgi,1,3 Mika Shirasu,1,3 Yusuke Ogura,1 Yukei Hirasawa,1 Masako Okamoto,1 Rieko Kawamura,1
Hirosato Takikawa,1 and Kazushige Touhara1,2,4,*
1Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Tokyo, Japan
2International Research Center for Neurointelligence (WPI-IRCN), Institutes for Advanced Study, The University of Tokyo, Tokyo, Japan
3These authors contributed equally
4Lead contact
*Correspondence: ktouhara@g.ecc.u-tokyo.ac.jp
https://doi.org/10.1016/j.isci.2025.113087
SUMMARY
In many animals, olfactory cues are key for recognizing their sex, kinship, and reproductive status. In humans,
it has been proposed that males find female body odors attractive during the ovulatory phase, coinciding with
peak fertility. However, the molecular and physiological foundations of this attractiveness remain unclear. In
this study, we combined sensory evaluation with gas chromatography-mass spectrometry for chemical anal
ysis and identified volatile compounds that fluctuate across menstrual cycle phases. We identified three
ovulatory phase-increasing compounds that mitigate the basal unpleasant axillary odor, resulting in the
ovulatory-phase body odor being perceived as the most pleasant by males. Additionally, these compounds
alleviated hostility and stress induced by the basal axillary odor, leading to relaxation in males and an
enhanced positive impression of female facial images. Our results suggest that several ovulatory phase-
increasing compounds may contribute to the positive effects of ovulatory female odor on males, as proposed
in previous research.
iScience ▪▪, 113087, ▪▪, 2025 © 2025 The Authors. Published by Elsevier Inc. 1
This is an open access article under the CC BY-NC license (http://creativecommons.org/licenses/by-nc/4.0/).
Please cite this article in press as: Ohgi et al., Human ovulatory phase-increasing odors cause positive emotions and stress-suppressive effects in
males, iScience (2025), https://doi.org/10.1016/j.isci.2025.113087
ll iScience
OPEN ACCESS Article
Figure 1. Sensory evaluation of female axil
lary odor during the menstrual cycle
Pleasantness (A) and intensity (B) ratings of 84
axillary odor samples collected from 21 females
during four phases (M: menstrual, F: follicular, O:
ovulatory, L: luteal) by male raters (n = 21).
Pleasantness was rated on a scale from -4
(extremely unpleasant) to +4 (extremely pleasant),
and intensity was rated on a scale from 0 (not
perceived) to 6 (extremely intense). Error bars
show standard errors. The effect of menstrual
phase on rating values was examined using a
linear mixed model (LMM), with rating values as
the dependent variable and menstrual cycle phase
as the explanatory variable. Statistical signifi
cance of the differences between the O phase and
each of the other phases was evaluated by testing
the significance of the beta values (see STAR
Methods for details). *p < 0.05 and **p < 0.01.
See also Figure S1 and Table S1.
RESULTS were excluded because they may not accurately assess odor
quality (STAR Methods). The selection rate of 15 descriptors,
Axillary odor during the ovulatory phase is more based on quantitative descriptive analysis (QDA, see STAR
pleasant than that during the menstrual phase Methods), chosen by male raters (n = 16) to describe the axillary
To assess the intensity and pleasantness of axillary odor across samples is presented in Table 1. Among the 15 descriptors, with
different phases of the menstrual cycle, axillary odors were a focus on those with higher frequency of selection, ‘‘vinegary
collected on gauze from 21 female donors at each of the four odor’’ was checked 24–28 times in M, F, and L-phases, and 17
stages of their menstrual cycle (i.e., M[enstrual], F[ollicular], O times in the O-phase, with a lower number of selections in the
[vulatory], L[uteal]-phases; a total of 84 odor samples) O-phase. In contrast, ‘‘fragrant odor’’ was examined 37–40 times
(Figure S1A). Each phase of the menstrual cycle was determined in the M-, F-, and L-phases, and 53 times in the O-phase, with a
by the luteinizing hormone (LH) surge and body temperature, and higher number of selections in the O-phase. All these results sug
the amounts of estradiol and progesterone in saliva were evalu gest that axillary odors during the O-phase are relatively
ated for some donors (Figure S1B). Given that previous reports pleasant.
suggest differences in body odor under stress from normal
body odor,14,15 we used the state-trait anxiety inventory (STAI) Chemical composition analysis of axillary odor identified
score to assess the anxiety state of the donors at the time of increasing or decreasing compounds for each phase
sample collection. No significant difference in anxiety state, as To explore the chemical profiles of female axillary odors during
indicated by the STAI score, was observed among the four the M-, F-, O-, and L-phases, we analyzed volatiles emitted
phases (see STAR Methods). Out of the 84 odor samples, each from the axilla during the menstrual cycle in 21 female donors,
male rater (n = 21) provided odor ratings (intensity and pleasant the same donors as in the previous section, using GC/MS.
ness) on samples from 10 to 11 donors (40–44 samples in total). Initially, we compared four different adsorbents, including sterile
Linear mixed model (LMM) analyses revealed that axillary odor gauze, polydimethylsiloxane (PDMS) membranes, Monotrap,
during the O-phase was rated significantly more pleasant than and Twister, for a female (Table S2 and Method for detailed sam
those during M and L-phases (M vs. O: t(846) = 3.1, p = 0.006, ple preparation and analytical methods for each adsorbent). We
L vs. O: t(846) = 2.4, p = 0.048, Figure 1 and Table S1). There found that the PDMS membrane adsorbed the largest number of
was also a tendency for pleasantness to be higher compared compounds (Figure S2) and that its close fit to the armpit allowed
to the F-phase (F vs. O: t(846) = 2.2, p = 0.077, Figure 1 and for the efficient adsorption of low-volatile compounds such as
Table S1). When excluding raters with lower olfactory ability long-chain and unsaturated fatty acids secreted from the axillary
(Methods), the O-phase showed significantly higher pleasant surface (see Tables S3 and S4). Therefore, we used the PDMS
ness compared to the M, F, and L-phases (M vs. O: t(636) = membrane to sample axillary odors.
2.8, p = 0.017, F vs. O: t(636) = 2.5, p = 0.043, L vs. O: t(636) = Subsequently, 84 axillary odor samples (four-phase samples
2.4, p = 0.045, Table S1). These findings suggest that the from 21 female donors) collected using the PDMS membrane
O-phase tends to have a pleasant odor, with this tendency being were subjected to analysis using GC/MS combined with a ther
more pronounced among raters with higher olfactory abilities. mal desorption system (Figure S3). Among numerous peaks
No statistically significant differences in intensity were observed identified following GC/MS analyses, we focused on compounds
among the phases. that were structurally identified using an authentic sample or that
Next, the axillary odor quality at each of the four phases was existed in more than five donors, at least in one phase. The re
investigated using the check-all-that-apply (CATA) sensory sulting 98 peaks were further analyzed (average occurrence
profiling method. Raters with low odor identification quality rate in all phases of all compounds = 71%, see STAR
iScience ll
Article OPEN ACCESS
Table 1. Sensory evaluation of the odor quality of female axillary odor during the menstrual cycle
Male raters, n = 16
Number of selections Selection rate (%) Cochran’s
Descriptor Q test Effect
Descriptor selection rate (%) M F O L M F O L (p value) size (η2Q)
Pleasant Fragrant 25% 37 37 53 40 22% 22% 32% 24% 0.044 0.016
Woody/Pencil-like 11% 16 23 15 17 23% 32% 21% 24% 0.41 0.0056
Sweet 6% 11 8 12 9 28% 20% 30% 23% 0.73 0.0026
Milky 5% 11 7 11 8 30% 19% 30% 22% 0.65 0.0032
Citrus 4% 5 6 13 6 17% 20% 43% 20% 0.14 0.011
Grassy/Green 4% 6 4 10 6 23% 15% 38% 23% 0.31 0.0070
Earthy 3% 4 7 4 2 24% 41% 24% 12% 0.22 0.0086
Yogurt 2% 4 0 8 3 27% 0% 53% 20% 0.027 0.018
Unpleasant Sweaty 29% 51 55 45 46 26% 28% 23% 23% 0.40 0.0058
Vinegary 14% 26 24 17 28 27% 25% 18% 29% 0.17 0.010
Sebaceous 13% 21 23 21 25 23% 26% 23% 28% 0.82 0.0019
Damp clothing 8% 15 13 12 14 28% 24% 22% 26% 0.90 0.0019
Dusty 6% 7 10 11 11 18% 26% 28% 28% 0.72 0.0027
Musty 5% 13 10 3 7 39% 30% 9% 21% 0.048 0.016
Stinky 5% 9 8 7 8 28% 25% 22% 25% 0.96 0.00061
Evaluation of the sensory profile of female axillary odors at each phase of the menstrual cycle using the check-all-that-apply (CATA) method (male
raters, n = 16). Descriptor selection rate indicates ‘‘number of descriptor selections/number of raters (16 males) × number of axillary odor samples
evaluated by each male rater (40 or 44 samples).’’ Number of selections indicates the ‘‘number of descriptor selections for each phase samples.’’ Se
lection rate indicates ‘‘number of selections/number of total selection from M to L phase.’’ Cochran’s Q-test was used to evaluate the differences in
axillary odors between the four different phases.
Methods). To mitigate the effect of inter-individual data vari p values were 1,3-butanediol (ID41) (Z = 2.6, p = 0.005, r = 0.56),
ability, the peak area values were standardized using a robust octadecane (ID#45) (Z = 2.7, p = 0.0032, r = 0.60), and
standardization method. To discern changes in the secretion of 1-pentadecanol (ID#72) (Z = 2.6, p = 0.0043, r = 0.58)
the 98 compounds during the menstrual cycle, Ward’s method (Figures 2 and 3; Tables S3 and S4). In the L-phase, four com
cluster analysis was performed using the mean of the robust Z pounds significantly increased, and the identified compounds
scores for each time period. The 98 compounds were arranged were hydrogen sulfide (ID#1) (Z = 2.1, p = 0.02, r = 0.45),
so that compounds with similar changes between the four 2-pyrrolidone (ID#58) (Z = 2.3, p = 0.01, r = 0.51), and hexahy
phases were located in close proximity (Figure 2). This analysis dro-2H-azepin-2-one (ID#68) (Z = 1.9, p = 0.029, r = 0.41)
reveals several groups of compounds that tend to decrease dur (Tables S3 and S4).
ing the M-phase (top cluster) or increase during the O- and
F-phases (top or middle cluster, respectively). Men perceive the odor of ovulatory phase-increasing
To identify increasing or decreasing compounds at each compounds as more pleasant and comforting
phase, we searched for compounds that significantly increased In conjunction with the sensory evaluation (refer to Figure 1) and
above the median of all phases for each phase. A one-sample GC/MS analysis (refer to Figures 2 and 3), (E)-geranylacetone,
Wilcoxon signed-rank test (one-tailed) was performed to deter tetradecanoic acid, and (Z)-9-hexadecenoic acid emerged as
mine whether the amount of the compound in each phase was potential compounds contributing to the reduction of the un
greater than 0 (median of the robust Z score). Out of the 98 com pleasantness of axillary odor during the O-phase. To probe
pounds, seven significantly increased during the O-phase how the increase in these three compounds during the
(Figures 2 and 3; Tables S3 and S4). Three compounds were O-phase influences the quality of axillary odor for men, we exam
identified as (E)-geranylacetone (ID#49) (Z = 2.1, p = 0.016, r = ined changes in odor quality when the odor-O (odor-O: three
0.47), tetradecanoic acid (ID#86) (Z = 1.7, p = 0.043, r = 0.38), O-phase-increasing compounds) was introduced to a basal axil
and (Z)-9-hexadecenoic acid (ID#94) (Z = 2.1, p = 0.019, r = lary odor, closely mimicking axillary odor16 (odor-A: a basal axil
0.45) (Figure S4). In the M-phase, four compounds significantly lary odor), consistently present in the female axilla throughout
increased, and the identified compounds were octanal (ID#15) the menstrual cycle with little or no difference (refer to Star
(Z = 1.9, p = 0.031, r = 0.41) and 1-tetradecanol (ID#65) (Z = Methods and Table S6). Sensory evaluations for four odor stimuli
1.7, p = 0.045, r = 0.37) (Tables S3 and S4). The remaining two (Ctrl, no odor; odor-O, three O-phase-increasing compounds;
compounds with small p value (like ID#2) were unknown (not odor-A, basal axillary odor; odor-A+O, a mixture of odor-A and
identified) but representative m/z was obtained (Table S10). In odor-O; Table S6) were conducted on 72 male raters. The eval
the F-phase, 26 compounds significantly increased, and of the uation terms included the 15 descriptors in Table 1, along with
15 compounds that could be identified, the three with the lowest intensity, pleasantness, preference, feminine/masculine, relax,
ll iScience
OPEN ACCESS Article
iScience ll
Article OPEN ACCESS
and arousal/excitement. Figure 4 illustrates the results for 13 tering other participants, sat calmly in a chair, and were initially
terms where significant differences were observed in the sensory asked to perform a baseline multiple mood scale (MMS,
evaluation (all results are available in Table S7). Odor-O was Table S8) (a type of psychological test) assessment and saliva
weakly perceptible in intensity (t(64) = 5.4, q < 0.0005, r = collection (base condition). In session 1–3, the raters then wore
0.56), and the overall odor quality was deemed pleasant with a headset with a microphone, in which the experimenter
fragrant (t(64) = 3.9, q < 0.005, r = 0.44), citrus (t(64) = 3.1, applied an odor to the microphone portion without informing
q < 0.05, r = 0.36), and relaxing (t(60) = 3.6, q < 0.005, r = the raters, and rated female face images (impression rating)
0.42). The addition of odor-O to odor-A led to significantly higher and completed MMS, followed by saliva collection immediately
ratings for pleasantness (t(71) = 10.1, q < 0.0005, r = 0.77) and after removing the headset. The 19 male rater had counter-
preference (t(71) = 9.5, q < 0.0005, r = 0.75), a more feminine rat balanced exposure to Ctrl or odor-O in session 1 and 2. Of
ing for feminine/masculine (t(69) = 4.9, q < 0.0005, r = 0.51), and the 19 male raters, 9 were exposed to odor-A, and 10 were
lower ratings for intensity (t(71) = 4.0, q < 0.0005, r = 0.43) (odor-A exposed to odor-A+O in session 3.
+O in Figure 4A). Fragrant (t(71) = 2.9, q < 0.05, r = 0.32), citrus (t The analysis of MMS in sessions 1 and 2 (Figure 5A) revealed
(71) = 2.6, q < 0.05, r = 0.29), woody (t(71) = 3.7, q < 0.005, r = that ‘‘hostility’’ was significantly lower in the odor-O condition
0.40), sweet (t(71) = 2.3, q < 0.05, r = 0.27), and relaxing compared to the Ctrl condition (no odor) (q = 3.7, p = 0.05,
(t(65) = 3.6, q < 0.005, r = 0.41) received significantly higher rat d = 0.66) (Figure 5B, top panel, second from the left; Table
ings, whereas sweaty (t(71) = 5.0, q < 0.0005, r = 0.51), damp S14). The impression rating tasks during sessions showed a sig
clothing (t(71) = 2.9, q < 0.05, r = 0.33), sebaceous (t(71) = 2.3, nificant increase in ‘‘boredom’’ (q = 4.2, p = 0.025, d = 0.90) and a
q < 0.05, r = 0.26), and stinky (t(71) = 4.8, q < 0.0005, r = 0.50) decrease in ‘‘concentration’’ (q = 4.5, p = 0.016, d = 0.84) in the
received significantly lower ratings in odor-A+O (Table S12). Ctrl condition compared to the base condition (before session 1),
This observation aligns with the finding that male raters associ with no significant difference observed between the odor-O and
ated fragrant and citrus scents when rating the ovulatory female base conditions (Figure 5B, top panel; Table S14). These findings
axillary odor (Table 1). suggest that odor-O has a positive psychological effect that
We posited a hypothesis that the three O-phase-increasing might decrease the feeling of ‘‘hostility.’’
compounds (odor-O) are pivotal female axillary odors evoking In session 3, a comparison of the scales for each odor condi
pleasant and positive emotional changes in males. Thus, we tion and the corresponding control condition demonstrated
tested the effects of adding each of the odor-O compounds a significant decrease in ‘‘relaxation’’ in the odor-A condition
to odor-A for 37 male raters in terms of pleasantness. A single (t(6) = 2.5, p = 0.047, d = 0.69). However, this decrease was miti
odor-O compound significantly reduced the unpleasantness of gated by adding odor-O into odor-A (odor-A+O condition)
odor-A, although the combined effect of the odor-O mix was (Figure 5B, middle and lower panel, forth column from the right;
more pronounced than that of a single compound (Figure 4B, Table S14). These results suggest that odor-O has a positive
A + O vs. A + G: t(36) = 2.1, p = 0.016, d = 0.48; A + O vs. A + psychological effect that increases the feeling of relaxation
T: t(36) = 1.8, p = 0.023, d = 0.45; and A + O vs. A + H: t(36) = when added to the basal axillary odor.
2.1, p = 0.012, d = 0.48, odor-G, T, and H are (E)-geranylace To delve into the relaxation effect of odor-O, we focused on
tone, tetradecanoic acid, and (Z)-9-hexadecenoic acid, α-amylase in saliva, known to be elevated in malodor stress.17
respectively; Tables S6, S12, and Figure S5). To evaluate We observed a significant negative correlation between
whether the same unpleasantness-suppression effect could α-amylase and the ‘‘relaxation’’ score in MMS at all time points
be observed in compounds that increase during phases other (p = 0.0005) (base, session 1 to 3, n = 19 × 4, Figure 5C), indi
than the O-phase, a similar experiment was conducted using cating that the amount of α-amylase in saliva increases as the
the three compounds that increase during the F-phase ‘‘relaxation’’ score decreases. Comparing α-amylase levels in
(odor-F). Odor-F also demonstrated the effect of reducing the the base, Ctrl (no odor), and odor-O conditions revealed no sig
unpleasantness of the basal axillary odor; however, the effect nificant differences (Figure 5D, left). Comparing the Ctrl(A) (Ctrl
was not as significant as that of odor-O (Figure 4C, A + O vs. data from raters who smelled odor-A in session 3, n = 7) and
A + F: q = 6.2, p < 0.005, d = 0.70). odor-A conditions showed a significant increase in α-amylase
under odor-A (t(6) = 2.9, p = 0.027, d = 0.94, Figure 5D, middle,
Psychological and physiological changes induced by the Table S14), whereas this increase was suppressed by adding
O-phase-increasing compounds on men odor-O (odor-A+O condition) (Ctrl data from raters who smelled
To investigate the psychological and physiological effects of A + O in session 3, n = 10). These results demonstrate that
the O-phase-increasing odor, we conducted an experiment odor-O has the effect of suppressing the increase in salivary
with the timeline depicted in Figure 5A. In essence, male raters amylase caused by a baseline axillary odor, correlating with
(n = 19) entered a private experimental room to avoid encoun the level of relaxation.
Figure 2. The secretion levels of the 98 compounds in each phase of the menstrual cycle
Cluster analysis with the Ward method was performed using the mean of the robust Z scores for each time period, and the 98 compounds were arranged such
that compounds showing similar changes between the four phases were located in close proximity to each other. The box extends from the 25th to the 75th
percentiles, and individual scores are plotted as points superimposed on the graph. The pink arrows indicate the three compounds that increased during the
ovulatory phase, and the yellow arrows indicate the three compounds that increased during the follicular phase used in the following experiments.
See also Tables S2, S3, S4, and S5; Figures S2 and S3.
ll iScience
OPEN ACCESS Article
Figure 3. Variability in volatile compounds that significantly increased during ovulatory or follicular phases
Line graphs show the changes in the amounts of the three O-phase-increasing compounds (#49, #86, and #94) and the three identified compounds (#41, #45, and
#72) during the menstrual cycle. The chemical structure of each compound is shown in the graph. The y axis shows the robust Z score of the peak area of each
compound obtained using GC/MS analysis. The box extends from the 25th to the 75th percentiles; individual scores are plotted as points superimposed on the
graph, and the gray lines indicate the connection of each donor’s score. Statistical analyses were performed using a one-sample Wilcoxon signed-rank test (one-
tailed) to determine whether the amount of each compound was greater than 0 (median of robust Z score). The significance level was set at p < 0.05. *p < 0.05, and
**p < 0.01.
See also Figure S4.
Effects of the O-phase-increasing compounds on male the average ratings for these categories (Figure 5E, red line) and
evaluation of female faces employed it as the ‘‘attractiveness’’ in the subsequent analysis.
In the same experimental timeline (Figure 5A), we examined the Face images were divided into two groups based on attractive
impact of the O-phase-increasing odor on the impressions male ness under Ctrl conditions: a low-rated group (4 images) and a
raters had for female faces. Images of eight female faces were high-rated group (4 images). It was observed that the attractive
prepared (50% smiling, Figure 5E), and raters were instructed ness of faces in the low-rated group decreased with odor-A and
to rate their impression of each face in four categories (cate significantly increased with odor-A+O compared to the no-odor
gories: beautiful, elegant and intellectual, want to spend time condition (F (1, 8) = 9.2, p = 0.016, η2 = 0.004, Figure 5F, left,
with, want to keep gazing at). Because all four categories of 8 im Table S14). When the differences between odor-A and no-odor
ages exhibited similar rating trends under Ctrl conditions (odor-A - Ctrl(A) in Figure 5F) and odor-A+O and no-odor (A +
(Figure 5E, purple, blue, green, and orange lines), we calculated O - Ctrl(A + O) in Figure 5F) were calculated and compared, it
iScience ll
Article OPEN ACCESS
Figure 4. Three ovulatory phase-increasing compounds positively affected the unpleasant quality of basal axillary odor
(A) Results of sensory evaluation for four odor stimuli (Ctrl: no odor; odor-O: O-phase-increasing three compounds; odor-A: a basal axillary odor; odor-A+O:
mixture of odor-A, and odor-O, Table S6). The y axis indicates the rating score (male, n = 61–72, mean ± SEM). Paired t test with FDR for all descriptors. *q < 0.05,
**q < 0.01, ***q < 0.005, and ****q < 0.001, #q < 0.0005. Intensity (0 = not perceived to +6 = extremely intense), pleasantness (− 3 = very unpleasant to +3 = very
pleasant), preference (− 3 = extremely dislike to +3 = extremely like), feminine/masculine (− 3 = very masculine to +3 = very feminine), relax (0 = not perceived
to +6 = extremely perceived), arousal/excitement (0 = not perceived to +6 = extremely perceived), and the top 15 descriptors in Table 1 (0 = not perceived to +6 =
extremely perceived). See also Tables S6 and S7.
(B) Evaluation of the intensity and pleasantness of each odor-O component y axis indicates the rating score (male, n = 37, mean ± SEM). (odor-G,-T, and-H: (E)-
geranylacetone, tetradecanoic acid, and (Z)-9-hexadecenoic acid, respectively; odor-A+G, A + T, A + H: mixtures of odor-A and odor-G, T, and H, respectively).
The odor-A results are shown for reference. Dunnett’s multiple comparison test was used for comparisons. *p < 0.05 and **p < 0.01. See also Figure S5.
(C) Evaluation of the intensity and pleasantness of each stimulus. y axis indicates the rating score (male, n = 37, mean ± SEM). (odor-F: mixture of three increasing
F-phase-compounds; odor A + F: mixture of odor-A and odor-F). The odor-F results are shown for reference. Tukey’s multiple-comparison test was used.
***p < 0.005 and ****p < 0.001.
ll iScience
OPEN ACCESS Article
iScience ll
Article OPEN ACCESS
was evident that the attractiveness of low-rated images was palmitoleic acid, it is nearly odorless. It is derived from palmitic
significantly higher under odor-A+O than under odor-A (F (1, acid and is presumed to be a precursor to (E)-2-nonenal, a repre
16) = 15, p = 0.0013, η2 = 0.1, Figure 5F, left). In the faces of sentative odor of ‘‘aging odor,’’ which is generated through its
the high-rated group, there were no significant differences be breakdown by skin bacteria.34,35 Although the three compounds
tween the odor stimuli condition and the no-odor state aforementioned increased during the ovulatory phase, it is
(Figure 5F, right). These results suggest that when odor-O was important to note that this study was designed as an exploratory
added to the basal axillary odor, a positive behavioral effect investigation. To minimize type II error, multiple testing correc
was observed in terms of impression ratings. tion was not applied, which unfortunately increased susceptibil
ity to type I error.
DISCUSSION Changes in axillary volatile compounds within the menstrual
cycle may be attributed to fluctuations in hormone concentra
This study investigated the role of olfactory cues in human inter tions, including estradiol and progesterone, within the cycle. In
sexual communication, with a focus on the ovulatory phase of particular, estradiol receptors are expressed in many tissues
the menstrual cycle, and identified ovulatory phase-increasing and are believed to be involved in systemic metabolism and
odorants that elicited positive emotional, physiological, and sweat gland secretory activity. Studies have revealed that amino
behavioral effects, facilitating male-female interaction. acids and their derivatives in blood tend to increase during the
Key points to consider when analyzing body odors are compre menstrual phase and decrease during the luteal phase, whereas
hensiveness and reproducibility in a well-controlled condition. phospholipids and other substances in blood increase during the
Conventional research on axillary odor often relies on a single follicular phase and decrease during the luteal phase.36 In this
collection or analytical method, resulting in the detection of a study, we demonstrate that long-chain fatty acids, including tet
limited number of compounds.18–21 In this study, we compared radecanoic acid and (Z)-9-hexadecenoic acid, increase from the
four collection adsorbents and GC/MS pretreatment methods. follicular phase to around ovulation and decrease during the
Based on the diversity of collected compounds and consistency luteal phase. This possibly reflects the degradation of blood-
with previous findings,22 we adopted a method combining derived lipids on the skin surface, which increase during the
PDMS and a thermal desorption system with GC/MS analysis. follicular phase. Moreover, several compounds also increased
This method proved to be the most efficient in detecting the during the follicular phase, and some of these tended to further
largest number of compounds among those tested and could increase during ovulation, suggesting the strong effects of estra
reproducibly detect compounds from an average of more than diol on the variability of amounts of each compound during the
70% of the donors. However, since competitive effects are known menstrual cycle.
to occur with adsorbents, the composition of the compounds ad This study revealed that the odor of the three O-phase-
sorbed on PDMS cannot be said to be identical to that of the orig increasing compounds induces a sense of comfort and relaxa
inal axillary odor. Through careful optimization of the analytical tion in men. Furthermore, when the odor of the basal axillary
method allowed for the identification of compounds during the odor was combined with the three O-phase increasing com
menstrual cycle that affect male-female interaction. pounds, it not only enhanced the attractiveness of low-attrac
The current knowledge on the identified three O-phase- tive-level female faces but also suppressed the increase in sali
increasing compounds is as follows: (1) (E)-geranylacetone: it vary amylase, an indicator of stress-related bad odor, in
is produced by oxidation of squalene and emits a floral or green comparison to the odor of the basal axillary scent alone. It has
odor quality and is among the prominently recognized sub been reported that pleasant essential oil (like rose) enhance
stances in humans.23–27 (2) Tetradecanoic acid: it has an facial attractiveness,37 but it is noteworthy that such effects
oxidized sebum odor quality and is detected in skin, saliva,28 may also be present in components of body odor that fluctuate
breast milk, and amniotic fluid. Its primary production pathway with the menstrual cycle. If other compounds that increase dur
is believed to be triglyceride hydrolysis by skin bacteria lipase,29 ing ovulation that could not be identified in this study can be clar
or it may be derived from blood owing to its ubiquitous presence ified, the function of ovulatory axillary odor will be elucidated in
in cells and plasma.30,31 Human amniotic fluid and breast milk more detail. Previous reports have suggested that when others
contain this odorant to which newborns exhibit a suckling smell body odor collected in high-stress situations, they experi
response.32,33 (3) (Z)-9-hexadecenoic acid: commonly termed ence anxiety from ambiguous facial expressions,38 and anxiety
(B) The y axis indicates the score for each mood on the MSS for each condition by male raters. The first row shows base (basal condition) vs. Ctrl (no odor) vs.
odor-O (three ovulatory-phase increasing compounds) (n = 19), the second row shows Ctrl(A) vs. odor-A (basal axillary odor) (n = 9), and the third row shows Ctrl
(A + O) vs. odor-A+O (mixture of A and odor-O) (n = 10). *p < 0.05, **p < 0.01, and ***p < 0.005.
(C) Salivary α-amylase and ‘‘relaxation’’ scores in MMS at all time points (baseline and sessions 1–3, n = 76) are negatively correlated by simple linear regression
analysis. This indicates that the amount of α-amylase in the saliva increased with lower relaxation in this experiment. Dotted lines are 95% confidence intervals of
the best fit lines.
(D) The y axis shows the amount of salivary α-amylase under each condition. α-amylase significantly increased under odor-A stimulation (n = 7) but not under
odor-A+O (n = 10). *p < 0.05.
(E) Ratings for images of 8 females by male raters (n = 19) under Ctrl conditions for four impression categories (purple, blue, green, and orange lines). The average
values of these categories (red line) are indicated (mean ± SE).
(F) Average rating of a low-rated group (8 images) and a high-rated group (8 images) under each condition and the differences between odor and Ctrl conditions
(mean ± SE). *p < 0.05 and **p < 0.01.
ll iScience
OPEN ACCESS Article
and nervousness increase,39 whereas when others smell body Material availability
odor collected in a happy state, a facial expression and percep This study did not generate new unique reagents.
tual processing style indicative of happiness were shown to be
induced.40 Additionally, women who smelled their partner’s Data and code availability
body odor experienced less stress.41,42 Body odors from new • GC/MS raw data supporting the current study will not be deposited in a
borns has been reported to activate reward-related cerebral public repository because we are still analyzing the data for projects not
areas in their mothers43 and mother’s odor did not cause typical reported in this study; however, these data will be available from the
brain responses to fearful stimuli in their child.44 These findings corresponding author upon reasonable request.
It is known that the composition of axillary odor is partially influ • METHOD DETAILS
enced by genetic factors,50 and that the perception of body odor ○ Menstrual cycle assessment
○ Stress related emotional state assessment in axillary odor sampling
is also affected by genetic background51,52 and cultural
○ Collection of axillary odor for sensory evaluation
context.53,54 Therefore, by conducting sensory evaluations and ○ Procedure for sensory evaluation of female axillary odor
chemical analysis of axillary odor across diverse ethnic groups ○ Odor questionnaire for axillary odor
beyond the Japanese population, while accounting for the ef ○ Generation of descriptive terms
tion, it should be possible to perform a detailed investigation ○ Collection of axillary odors to compare the four phases of menstru
iScience ll
Article OPEN ACCESS
SUPPLEMENTAL INFORMATION 15. Mujica-Parodi, L.R., Strey, H.H., Frederick, B., Savoy, R., Cox, D., Bota
nov, Y., Tolkunov, D., Rubin, D., and Weber, J. (2009). Chemosensory
Supplemental information can be found online at https://doi.org/10.1016/j.isci. Cues to Conspecific Emotional Stress Activate Amygdala in Humans.
2025.113087. PLoS One 4, e6415. https://doi.org/10.1371/journal.pone.0006415.
16. Hara, T., Kyuka, A., and Shimizu, H. (2015). Butane-2,3-dione: the Key
Received: June 4, 2024 Contributor to Axillary and Foot Odor Associated with an Acidic Note.
Revised: September 6, 2024 Chem. Biodivers. 12, 248–258. https://doi.org/10.1002/cbdv.201400272.
Accepted: July 8, 2025
17. Hirasawa, Y., Shirasu, M., Okamoto, M., and Touhara, K. (2019). Subjec
tive unpleasantness of malodors induces a stress response. Psychoneur
REFERENCES oendocrinology 106, 206–215. https://doi.org/10.1016/j.psyneuen.2019.
03.018.
1. O’Connell, R.J., Singer, A.G., Macrides, F., Pfaffmann, C., and Agosta, W.
18. Zeng, X.N., Leyden, J.J., Spielman, A.I., and Preti, G. (1996). Analysis of
C. (1978). Responses of the male golden hamster to mixtures of odorants
characteristic human female axillary odors: Qualitative comparison to
identified from vaginal discharge. Behav. Biol. 24, 244–255. https://doi.
males. J. Chem. Ecol. 22, 237–257. https://doi.org/10.1007/bf02055096.
org/10.1016/S0091-6773(78)93110-3.
19. Natsch, A., Derrer, S., Flachsmann, F., and Schmid, J. (2006). A broad di
2. Cerda-Molina, A.L., Hernández-López, L., Rojas-Maya, S., Murcia-Mejı́a, versity of volatile carboxylic acids, released by a bacterial aminoacylase
C., and Mondragón-Ceballos, R. (2006). Male-induced sociosexual from axilla secretions, as candicate molecules for the determination of hu
behavior by vaginal secretions in Macaca arctoides. Int. J. Primatol. 27, man-body odor type. Chem. Biodivers. 3, 1–20. https://doi.org/10.1002/
791–807. https://doi.org/10.1007/s10764-006-9045-0. cbdv.200690015.
3. Clarke, P.M.R., Barrett, L., and Henzi, S.P. (2009). What role do olfactory 20. Gallagher, M., Wysocki, C.J., Leyden, J.J., Spielman, A.I., Sun, X., and
cues play in chacma baboon mating? Am. J. Primatol. 71, 493–502. Preti, G. (2008). Analyses of volatile organic compounds from human
https://doi.org/10.1002/ajp.20678. skin. Br. J. Dermatol. 159, 780–791. https://doi.org/10.1111/j.1365-
4. Ziegler, T.E., Schultz-Darken, N.J., Scott, J.J., Snowdon, C.T., and Ferris, 2133.2008.08748.x.
C.F. (2005). Neuroendocrine response to female ovulatory odors depends 21. Curran, A.M., Rabin, S.I., Prada, P.A., and Furton, K.G. (2005). Compari
upon social condition in male common marmosets, Callithrix jacchus. son of the volatile organic compounds present in human odor using
Horm. Behav. 47, 56–64. https://doi.org/10.1016/j.yhbeh.2004.08.009. SPME-GC/MS. J. Chem. Ecol. 31, 1607–1619. https://doi.org/10.1007/
5. Doty, R.L., Ford, M., Preti, G., and Huggins, G.R. (1975). Changes in the s10886-005-5801-4.
intensity and pleasantness of human vaginal odors during the menstrual 22. Dormont, L., Bessière, J.M., and Cohuet, A. (2013). Human Skin Volatiles:
cycle. Science 190, 1316–1318. https://doi.org/10.1126/science. A Review. J. Chem. Ecol. 39, 569–578. https://doi.org/10.1007/s10886-
1239080. 013-0286-z.
6. Singh, D., and Bronstad, P.M. (2001). Female body odour is a potential cue 23. Leal, H.M., Hwang, J.K., Tan, K., and Leal, W.S. (2017). Attraction of Culex
to ovulation. Proc. Biol. Sci. 268, 797–801. https://doi.org/10.1098/rspb. mosquitoes to aldehydes from human emanations. Sci. Rep. 7, 17965.
2001.1589. https://doi.org/10.1038/s41598-017-18406-7.
�ek, J., Dvo�ráková, R., Barto�
7. Havlı́c s, L., and Flegr, J. (2006). Non-adver 24. Zhao, Z., Zung, J.L., Hinze, A., Kriete, A.L., Iqbal, A., Younger, M.A., Mat
tized does not mean concealed: Body odour changes across the human thews, B.J., Merhof, D., Thiberge, S., Ignell, R., et al. (2022). Mosquito
menstrual cycle. Ethology 112, 81–90. https://doi.org/10.1111/j.1439- brains encode unique features of human odour to drive host seeking. Na
0310.2006.01125.x. ture 605, 706–712. https://doi.org/10.1038/s41586-022-04675-4.
8. Gildersleeve, K.A., Haselton, M.G., Larson, C.M., and Pillsworth, E.G. �ska, N., Flynn, C., Ratcliffe, N., Belluomo, I., Myridakis, A., Gould,
25. Drabin
(2012). Body odor attractiveness as a cue of impending ovulation in O., Fois, M., Smart, A., Devine, T., and Costello, B.D. (2021). A literature
women: Evidence from a study using hormone-confirmed ovulation. survey of all volatiles from healthy human breath and bodily fluids: the hu
Horm. Behav. 61, 157–166. https://doi.org/10.1016/j.yhbeh.2011.11.005. man volatilome. J. Breath Res. 15, 034001. https://doi.org/10.1088/1752-
9. Haselton, M.G., and Gildersleeve, K. (2011). Can men detect ovulation? 7163/abf1d0.
Curr. Dir. Psychol. Sci. 20, 87–92. https://doi.org/10.1177/09637214 26. Syed, Z., and Leal, W.S. (2009). Acute olfactory response of Culex mosqui
11402668. toes to a human- and bird-derived attractant. Proc. Natl. Acad. Sci. USA
10. Gangestad, S.W., and Thornhill, R. (2008). Human oestrus. Proc. Biol. Sci. 106, 18803–18808. https://doi.org/10.1073/pnas.0906932106.
275, 991–1000. https://doi.org/10.1098/rspb.2007.1425. 27. Owsienko, D., Goppelt, L., Hierl, K., Schäfer, L., Croy, I., and Loos, H.M.
�ek, J., Cobey, K.D., Barrett, L., Klapilová, K., and Roberts, S.C.
11. Havlic (2024). Body odor samples from infants and post-pubertal children differ
(2015). The spandrels of Santa Barbara? A new perspective on the peri- in their volatile profiles. Commun. Chem. 7, 53. https://doi.org/10.1038/
ovulation paradigm. Behav. Ecol. 26, 1249–1260. https://doi.org/10. s42004-024-01131-4.
1093/beheco/arv064. 28. de Lacy Costello, B., Amann, A., Al-Kateb, H., Flynn, C., Filipiak, W.,
12. Kuukasjärvi, S., Eriksson, C.P., Koskela, E., Mappes, T., Nissinen, K., and Khalid, T., Osborne, D., and Ratcliffe, N.M. (2014). A review of the volatiles
Rantala, M.J. (2004). Attractiveness of women’s body odors over the men from the healthy human body. J. Breath Res. 8, 014001. https://doi.org/10.
strual cycle: The role of oral contraceptives and receiver sex. Behav. Ecol. 1088/1752-7155/8/1/014001.
15, 579–584. https://doi.org/10.1093/beheco/arh050. 29. Nicolaides, N. (1974). Skin Lipids: Their Biochemical Uniqueness. Science
13. Zetzsche, M., Weiß, B.M., Kücklich, M., Stern, J., Birkemeyer, C., Widdig, 186, 19–26. https://doi.org/10.1126/science.186.4158.19.
A., and Penke, L. (2024). Combined perceptual and chemical analyses 30. Aparicio, E., Martı́n-Grau, C., Hernández-Martinez, C., Voltas, N., Canals,
show no compelling evidence for ovulatory cycle shifts in women’s axillary J., and Arija, V. (2021). Changes in fatty acid levels (saturated, monounsat
odour. Proc. Biol. Sci. 291, 20232712. https://doi.org/10.1098/rspb. urated and polyunsaturated) during pregnancy. BMC Pregnancy Child
2023.2712. birth 21, 778. https://doi.org/10.1186/s12884-021-04251-0.
14. de Groot, J.H.B., Kirk, P.A., and Gottfried, J.A. (2021). Titrating the Smell of 31. Aktas, M., Elmastas, M., Ozcicek, F., and Yilmaz, N. (2016). Erythrocyte
Fear: Initial Evidence for Dose-Invariant Behavioral, Physiological, and Membrane Fatty Acid Composition in Premenopausal Patients with Iron
Neural Responses. Psychol. Sci. 32, 558–572. https://doi.org/10.1177/ Deficiency Anemia. J. Oleo Sci. 65, 225–231. https://doi.org/10.5650/
0956797620970548. jos.ess15211.
ll iScience
OPEN ACCESS Article
32. Contreras, C.M., Gutiérrez-Garcı́a, A.G., Mendoza-López, R., Rodrı́guez- creases during the fertile phase of the menstrual cycle. Proc. Biol. Sci.
Landa, J.F., Bernal-Morales, B., and Dı́az-Marte, C. (2013). Amniotic fluid 271, S270–S272. https://doi.org/10.1098/rsbl.2004.0174.
elicits appetitive responses in human newborns: Fatty acids and appetitive 49. Burriss, R.P., Troscianko, J., Lovell, P.G., Fulford, A.J.C., Stevens, M.,
responses. Dev. Psychobiol. 55, 221–231. https://doi.org/10.1002/dev. Quigley, R., Payne, J., Saxton, T.K., and Rowland, H.M. (2015). Changes
21012. in women’s facial skin color over the ovulatory cycle are not detectable
33. Gutiérrez-Garcı́a, A.G., Contreras, C.M., and Dı́az-Marte, C. (2017). Myris by the human visual system. PLoS One 10, e0130093. https://doi.org/
tic acid in amniotic fluid produces appetitive responses in human new 10.1371/journal.pone.0130093.
borns. Early Hum. Dev. 115, 32–37. https://doi.org/10.1016/j.earlhum 50. Toyoda, Y., Sakurai, A., Mitani, Y., Nakashima, M., Yoshiura, K.i., Naka
dev.2017.08.009. gawa, H., Sakai, Y., Ota, I., Lezhava, A., Hayashizaki, Y., et al. (2009).
34. Haze, S., Gozu, Y., Nakamura, S., Kohno, Y., Sawano, K., Ohta, H., and Earwax, osmidrosis, and breast cancer: why does one SNP (538G> A) in
Yamazaki, K. (2001). 2-Nonenal newly found in human body odor tends the human ABC transporter ABCC11 gene determine earwax type? FA
to increase with aging. J. Invest. Dermatol. 116, 520–524. https://doi. SEB J. 23, 2001–2013. https://doi.org/10.1096/fj.09-129098.
org/10.1046/j.0022-202x.2001.01287.x. 51. Keller, A., Zhuang, H., Chi, Q., Vosshall, L.B., and Matsunami, H. (2007).
35. Porter, N.A. (1986). Mechanisms for the Autoxidation of Polyunsaturated Genetic variation in a human odorant receptor alters odour perception.
Lipids. Acc. Chem. Res. 19, 262–268. https://doi.org/10.1021/ Nature 449, 468–472. https://doi.org/10.1038/nature06162.
ar00129a001. 52. Menashe, I., Abaffy, T., Hasin, Y., Goshen, S., Yahalom, V., Luetje, C.W.,
36. Draper, C.F., Duisters, K., Weger, B., Chakrabarti, A., Harms, A.C., and Lancet, D. (2007). Genetic elucidation of human hyperosmia to isova
Brennan, L., Hankemeier, T., Goulet, L., Konz, T., Martin, F.P., et al. leric acid. PLoS Biol. 5, e284. https://doi.org/10.1371/journal.pbio.
(2018). Menstrual cycle rhythmicity: metabolic patterns in healthy women. 0050284.
Sci. Rep. 8, 14568. https://doi.org/10.1038/s41598-018-32647-0. 53. Ferdenzi, C., Razafindrazaka, H., Baldovini, N., Poupon, D., Pierron, D.,
37. Seubert, J., Gregory, K.M., Chamberland, J., Dessirier, J.M., and Lund and Bensafi, M. (2019). Influence of gender and culture on the perception
ström, J.N. (2014). Odor valence linearly modulates attractiveness, but of acidic compounds of human body odor. Physiol. Behav. 210, 112561.
not age assessment, of invariant facial features in a memory-based rating https://doi.org/10.1016/j.physbeh.2019.112561.
task. PLoS One 9, e98347. https://doi.org/10.1371/journal.pone.0098347. 54. Schleidt, M., Hold, B., and Attili, G. (1981). A cross-cultural study on the
38. Zhou, W., and Chen, D. (2009). Fear-Related Chemosignals Modulate attitude towards personal odors. J. Chem. Ecol. 7, 19–31. https://doi.
Recognition of Fear in Ambiguous Facial Expressions. Psychol. Sci. 20, org/10.1007/bf00988632.
177–183. https://doi.org/10.1111/j.1467-9280.2009.02263.x. 55. Spielberger, C., Gorsuch, R., Lushene, R., Vagg, P.R., and Jacobs, G.
39. Albrecht, J., Demmel, M., Schöpf, V., Kleemann, A.M., Kopietz, R., May, (1983). Manual for the State-Trait Anxiety Inventory (Form Y1 – Y2), IV
J., Schreder, T., Zernecke, R., Brückmann, H., and Wiesmann, M. (Consulting Psychologists Press).
(2011). Smelling chemosensory signals of males in anxious versus non 56. Rossi, V., and Pourtois, G. (2012). Transient state-dependent fluctuations
anxious condition increases state anxiety of female subjects. Chem. in anxiety measured using STAI, POMS, PANAS or VAS: a comparative re
Senses 36, 19–27. https://doi.org/10.1093/chemse/bjq087. view. Anxiety Stress Coping 25, 603–645. https://doi.org/10.1080/
40. de Groot, J.H.B., Smeets, M.A.M., Rowson, M.J., Bulsing, P.J., Blonk, C. 10615806.2011.582948.
G., Wilkinson, J.E., and Semin, G.R. (2015). A sniff of happiness. Psychol. 57. Adams, J., Williams, A., Lancaster, B., and Foley, M. (2007). Advantages
Sci. 26, 684–700. https://doi.org/10.1177/0956797614566318. and uses of check-all-that-apply response compared to traditional scaling
of attributes for salty snacks. In 7th Pangborn Sensory Science Sympo
41. Granqvist, P., Vestbrant, K., Döllinger, L., Liuzza, M.T., Olsson, M.J.,
sium 12–16 August Minneapolis USA.
Blomkvist, A., and Lundström, J.N. (2019). The scent of security: Odor
of romantic partner alters subjective discomfort and autonomic stress re 58. Nishida, K., Kobayashi, M., Ogihara, H., Takeo, T., Kitano, M., and Take
sponses in an adult attachment-dependent manner. Physiol. Behav. 198, uchi, K. (2010). Clinical Usefulness of Smell Identification Test Card :
144–150. https://doi.org/10.1016/j.physbeh.2018.08.024. Open Essence. Nihon Jibiinkoka Gakkai Kaiho 113, 751–757. https://
doi.org/10.3950/jibiinkoka.113.751.
42. Hofer, M.K., Collins, H.K., Whillans, A.V., and Chen, F.S. (2018). Olfactory
cues from romantic partners and strangers influence women’s responses 59. Suzuki, M., Hashimoto, M., Yoshioka, M., Murakami, M., Kawasaki, K.,
to stress. J. Pers. Soc. Psychol. 114, 1–9. https://doi.org/10.1037/ and Urashima, M. (2011). The odor stick identification test for Japanese
pspa0000110. differentiates Parkinson’s disease from multiple system atrophy and pro
gressive supra nuclear palsy. BMC Neurol. 11, 157. https://doi.org/10.
43. Lundström, J.N., Mathe, A., Schaal, B., Frasnelli, J., Nitzsche, K., Gerber,
1186/1471-2377-11-157.
J., and Hummel, T. (2013). Maternal status regulates cortical responses to
the body odor of newborns. Front. Psychol. 4, 597. https://doi.org/10. 60. Fujio, H., Inokuchi, G., Fukuda, Y., Kuroki, S., Kokan, N., and Nibu, K.
3389/fpsyg.2013.00597. (2018). A New Screening Tool, An Odor Identification Test Card to Detect
Age–related Changes in Olfaction Hisami. Nippon Jibiinkoka Gakkai Kaiho
44. Schaal, B., Saxton, T.K., Loos, H., Soussignan, R., and Durand, K. (2020).
121, 38–43. https://doi.org/10.3950/jibiinkoka.121.38.
Olfaction scaffolds the developing human from neonate to adolescent and
beyond. Philos. Trans. R. Soc. Lond. B Biol. Sci. 375, 20190261. https:// 61. Stevens, D., Cornmell, R., Taylor, D., Grimshaw, S.G., Riazanskaia, S., Ar
doi.org/10.1098/rstb.2019.0261. nold, D.S., Fernstad, S.J., Smith, A.M., Heaney, L.M., Reynolds, J.C., et al.
(2015). Spatial variations in the microbial community structure and diver
45. Pipitone, R.N., and Gallup, G.G. (2008). Women’s voice attractiveness
sity of the human foot is associated with the production of odorous vola
varies across the menstrual cycle. Evol. Hum. Behav. 29, 268–274.
tiles. FEMS Microbiol. Ecol. 91, 1–11. https://doi.org/10.1093/fem
https://doi.org/10.1016/j.evolhumbehav.2008.02.001.
sec/fiu018.
46. Bryant, G.A., and Haselton, M.G. (2009). Vocal cues of ovulation in human 62. Riazanskaia, S., Blackburn, G., Harker, M., Taylor, D., and Thomas, C.L.P.
females. Biol. Lett. 5, 12–15. https://doi.org/10.1098/rsbl.2008.0507. (2008). The analytical utility of thermally desorbed polydimethylsilicone
47. Pavela Banai, I. (2017). Voice in different phases of menstrual cycle among membranes for in-vivo sampling of volatile organic compounds in and
naturally cycling women and users of hormonal contraceptives. PLoS One on human skin. Analyst 133, 1020–1027. https://doi.org/10.1039/
12, e0183462. https://doi.org/10.1371/journal.pone.0183462. b802515k.
48. Roberts, S.C., Havlicek, J., Flegr, J., Hruskova, M., Little, A.C., Jones, B. 63. Soini, H.A., Bruce, K.E., Wiesler, D., David, F., Sandra, P., and Novotny, M.
C., Perrett, D.I., and Petrie, M. (2004). Female facial attractiveness in V. (2005). Stir Bar Sorptive Extraction: A New Quantitative and
iScience ll
Article OPEN ACCESS
Comprehensive Sampling Technique for Determination of Chemical Signal 66. Terasaki, M., Klshimoto, Y., and Koga, A. (1992). Construction of a multiple
Profiles from Biological Media. J. Chem. Ecol. 31, 377–392. https://doi. mood scale. Japanese J. Psychol. 62, 350–356. https://doi.org/10.4992/
org/10.1007/s10886-005-1347-8. jjpsy.62.350.
67. Ogawa, T., Monchi, R., Kikuya, M., and Suzuki, N. (2000). Development of
64. Williams, C.M., and Mander, L.N. (2001). Chromatography with silver ni
the General Affect Scaled. Shinrigaku Kenkyu 71, 241–246. https://doi.
trate. Tetrahedron 57, 425–447. https://doi.org/10.1016/S0040-4020(00)
org/10.4992/jjpsy.71.241.
00927-3.
68. Nagy, T., van Lien, R., Willemsen, G., Proctor, G., Efting, M., Fülöp, M.,
65. Leleu, A., Demily, C., Franck, N., Durand, K., Schaal, B., and Baudouin, J. Bárdos, G., Veerman, E.C.I., and Bosch, J.A. (2015). A fluid response:
Y. (2015). The Odor Context Facilitates the Perception of Low-Intensity Alpha-amylase reactions to acute laboratory stress are related to sample
Facial Expressions of Emotion. PLoS One 10, e0138656. https://doi.org/ timing and saliva flow rate. Biol. Psychol. 109, 111–119. https://doi.org/10.
10.1371/journal.pone.0138656. 1016/j.biopsycho.2015.04.012.
ll iScience
OPEN ACCESS Article
STAR★METHODS
Ethics
The study received approval from the Ethics Committee of the University of Tokyo and adhered to the principles of the Declaration of
Helsinki. Informed written consent was obtained from all participants before conducting the experiments.
iScience ll
Article OPEN ACCESS
Odor raters
Sensory evaluation of female axillary odor
Twenty-one raters (males: n = 21, age: mean ± SD = 24 ± 2.9 years) were recruited for evaluating axillary odor samples. Raters were
heterosexual, single, free of respiratory, brain and nervous system, mental, and chronic diseases, not taking anabolic steroids, and
nonsmokers. Raters were informed that they would be evaluating various odors, including female-related odors. In this study, the
influence of sex and gender was not evaluated.
Odor evaluation using O-phase-increasing odor
Seventy-two single men (age: mean ± SD = 29 ± 5.4 years, the score of Open Essence Kit (OPE score): 10 ± 1.2) were recruited as
odor raters. All raters were heterosexual, did not suffer from respiratory, brain and nervous system, mental, and chronic diseases, did
not take anabolic steroids, and were nonsmokers. The ability to identify odors was measured using the Open Essence Kit, and all
raters scored 8 or higher, indicating no problem in identifying odor quality.
In this study, the influence of sex and gender was not evaluated.
METHOD DETAILS
ll iScience
OPEN ACCESS Article
indicates the presence of anxiety symptoms. The mean STAI-S scores for each phase were as follows: M-phase = 32.6 ± 5.48 points,
F-phase = 33.7 ± 9.16 points, O-phase = 34.9 ± 10.3 points, and L-phase = 34.9 ± 7.72 points. No significant difference in anxiety state
was observed among the four phases (p = 0.423, Friedman’s test).
(1) First, the sensory terms were generated via a pilot test comprising an open-ended questionnaire that was administered to 15
raters (male: n = 8, female: n = 7). They evaluated 17 samples of axillary odor during the menstrual cycle (menstrual phase: n =
4, follicular phase: n = 4, ovulatory phase: n = 6, luteal phase: n = 3). The raters sniffed the axillary odor samples in Teflon
squeeze bottles (Big Boy, Japan) and immediately described the quality of the odor after sniffing each sample. The list of de
scriptors consisted of 103 terms. Of the 103 terms, those used only once were omitted.
(2) The next procedure was conducted by five experts (male: n = 2, female: n = 3), involved in sensory evaluation and analytical
chemistry, to summarize the terms of the descriptor list. When several terms suggested a synonym, the most common word
was selected. After eliminating redundancies and controversial terms, the descriptor list had a total of 45 terms.
(3) To ensure that the 45 terms were sufficient to evaluate axillary odors, we conducted a pilot evaluation. Five raters assessed 21
odor samples (menstrual phase: n = 4, follicular phase: n = 4, ovulatory phase: n = 5, luteal phase: n = 8) and rated the 45 terms
on a seven-point scale (0 = not perceived to +6 = extremely perceived). We retained characteristic descriptors that were likely
to be important to discriminate between samples, as well as descriptors that were commonly perceived in most samples and
were considered important in assessing axillary odors. Furthermore, we re-confirmed the definitions of terms that were
believed to have multiple definitions and could cause individual perception-related differences. Terms with high-frequency
correlations were consolidated. Thus, the final selection resulted in 20 descriptors.
iScience ll
Article OPEN ACCESS
GCMS analysis
Four different sampling techniques and associated GC/MS analysis
To determine the appropriate sampling method for the analysis of volatile organic compounds (VOCs) in axillary odors, we compared
four different sampling methods. These included (1) sterile 100% cotton gauze (Oosaki Medical Corp., Japan), which is widely used
for body odor collection18,19; (2) PDMS membranes (Cat. No.: 751-624-16, Goodfellow Cambridge Limited) that are suitable for
capturing body odor VOCs due to their high adsorption capacity and flexibility61,62; (3) Monotrap (GL Sciences, Japan); and (4)
Twister (GERSTEL Inc., USA) adsorbents, which are widely used for odor analysis and are easy to handle.63
To compare each method, axillary odors were obtained from a female (laboratory member who met the inclusion/exclusion criteria
in the method of ‘‘female axillary odor donors,’’ age = 27) axilla using all four methods on the same day. The PDMS membranes were
first placed on the axillae surface, covered with PET film cut from Flek-Sampler (OMI ODOR-AIR SERVICE Corporation, Japan), and
secured to the skin using waterproof films (New Tac Kasei, Japan) for 2 h (refer to ‘‘collection of axillary odors to compare the four
phases of menstruation’’). Subsequently, sterile gauze (5 cm × 5 cm) was placed on the surface of the axillae and held in place for 2 h
using baby skin surgical tape. Monotrap or Twister was placed on the surface of the axillae and covered with a plate (TissueTec Cry
omold No. 2, Sakura Finetech Inc Japan.) and secured with waterproof film for 2 h. Following removal, PDMS membranes, Monotrap,
and Twister samples were transferred to a 2-mL glass vial (Thermo Fisher Scientific Inc., USA) and immediately frozen at − 30◦ C. The
gauze samples were immediately frozen at − 80◦ C in a 50-mL glass vial (NICHIDENRIKA-GLASS CO., LTD., Japan). While the sam
ples were attached to the axilla, the donor was instructed to change into clothes (UNIQLO, Japan) that were washed with fragrance-
free detergent.
We analyzed the axillary odor samples obtained from the four sampling methods using the appropriate equipment and conditions
for each sample. The sterile gauze sample was analyzed by GCMS-QP2010 Ultra (Shimadzu, Japan) (SH-Stabilwax column of
30 m × 0.53 mm i.d. with a film thickness of 1 μm) combined with the large volume static headspace system (7650HS-CTS, Entech,
USA). The gauze samples were placed in 500-mL inert glass vials (EN39-75500W, Entech, USA) and equilibrated in the oven for
20 min (80◦ C) before extracting the headspace into the capillary traps. The headspace (200 cc) was directed to the capillary traps,
followed by backflushing directly to the GC analytical column. The column temperature was programmed to rise from 40◦ C (5 min
hold) to 230◦ C (66 min hold) at 10◦ C/min. The interface temperature was maintained at 230◦ C, and the ion source temperature at
230◦ C. Mass spectra were obtained in full-scan mode (range: m/z 29–400) at 70ev electron impact.
PDMS membrane, Twister, and Monotrap samples were analyzed using GCMS-TQ8030 (Shimadzu, Japan) (Stabilwax column of
60 m × 0.32 mm i.d. with 0.5-μm film thickness) combined with a thermal desorption system OPTIC4 (GL Science, Japan). The con
ditions of thermal desorption controlled by OPTIC4 are as follows. The vaporization chamber increased from 35◦ C to 230◦ C at 5◦ C/s,
septum purge flow 5mL/min, sample sweep time 30 s, column flow 2mL/min, transfer time 2 min, split flow 10mL/min, cryotrap low
temperature − 150◦ C (130 s hold), cryotrap high temperature 250◦ C, cryotrap heat ramp rate 60◦ C/s, and the GCMS column temper
ature was programmed to rise from 50◦ C (2.5 min hold) followed by 10◦ C/min to 150◦ C, then 5◦ C/min to 230◦ C (61.5 min hold). The
interface temperature was maintained at 230◦ C, and the temperature of the ion source was 230◦ C. Mass spectra were obtained in
full-scan mode (range: m/z 29–400) at 70ev electron impact. Volatile compounds were detected using the automatic peak detector
(GCMS solution software; Shimadzu, Japan), and individual peaks were identified using the NIST 17 (National Institute of Standards
and Technology, USA) and the Wiley Registry of Mass Spectral Data, 9th edition (Wiley, USA). Visual inspection and fragment match
ing were used to confirm peaks with >70% similarity in the library search. Moreover, about half of the compounds were confirmed by
the MS and retention time of purchased authenticated compounds with >95% similarity. A peak retention index (RI) was subse
quently calculated to compare the analysis results for each condition. If a peak detected in the axillary odor sample was also detected
in the control sample (experimental room), only the peak detected with a peak area value more than twice that of the control was
adopted as the body odor compounds and the total number of such peaks were calculated (Figure S2). The threshold of more
than twice was set to avoid picking up noise while ensuring that no potential candidates were excluded.
ll iScience
OPEN ACCESS Article
We collected 84 axillary odor samples using PDMS membranes from November 2017 to March 2018, and all GC/MS data were
obtained from January 2018 to April 2018 using the aforementioned methods. Volatile compounds were detected using the auto
matic peak detector (GCMS solution software; Shimadzu, Kyoto, Japan), and individual peaks were verified using NIST 17 (National
Institute of Standards and Technology, USA) and the Wiley Registry of Mass Spectral Data, 9th edition (Wiley, USA). Compounds with
fatty acid characteristics, even if they did not match the database exactly, were labeled as ‘‘long-chain fatty acid’’ (Table S4) ‘‘Un
known compound’’ refers to a compound with an RI or mass spectrum that matched between samples but could not be identified.
Subsequently, the peaks were aligned, and the area of each peak was calculated using the GCMS solution software (Shimadzu,
Japan) by automatic integration. The compounds information is shown in Tables S10 and S11. For statistical analysis, 98 compounds
that could be identified by authentic compounds or had a minimum occurrence rate of 24% (equivalent to 5 donors) or more in the
phase in which they were detected in the largest number of people (average occurrence rate in all phases of all compounds = 71%)
were used. Of the 98 compounds, of which 52 compounds were identified by authentic compounds, 28 compounds were detected in
the control sample.
Because the samples were carefully stored and analyzed in a short period of time, there was little difference in the peak areas of the
internal standard toluene-D8 between axillary samples, so absolute peak areas were used in subsequent statistical analyses. When
compounds were not detected by GCMS solution (LOD&LOQ, a signal-to-noise ratio of 3 was considered), we set 0 for peak area.
iScience ll
Article OPEN ACCESS
Statistical analyses of 98 axillary odor compounds in each phase of the menstrual cycle
A priori power analysis using G*Power 3.1 indicated that 22 donors were enough in the condition of effect size d = 0.57, α level = 0.05,
Power (1-β err prob) = 0.8 in Wilcoxon signed-rank test (one sample case). We converted the peak area of 98 compounds into a robust
Z score for standardization. Robust z-scores were calculated as follows: robust z = (x – median)/normalized interquartile range
(NIQR), where ‘‘x’’ is the respective peak area value, and median and NIQR were those of the respective female donors. The cluster
ll iScience
OPEN ACCESS Article
analysis used in Figure 2 is ward method with squared euclidean distance. IBM SPSS version 27 (IBM Co., USA) was used for these
analyses. The p-value was calculated using a one-sample Wilcoxon signed-rank test (one-tailed) to determine whether the amount of
the compound in each phase was greater than 0 (median of the robust Z score), and the significance level was set at p < 0.05
(Tables S3 and S4). To test whether it is possible to discriminate between the ovulatory phase and other phases using the 98 com
pounds that were used in this analysis, conditional logistic regression on the corresponding data for 21 donors was performed using
the COXREG command of SPSS ver. 27 (Conditional logistic regression using COXREG (ibm.com)). First, the problem of multicolli
nearity was encountered, so the Spearman’s coefficient for the 98 compounds was calculated. Based on > 0.7, 13 variables (#5, 6,
14, 19, 24, 56, 83, 89, 90, 91, 92, 93, 94) were removed. A logistic regression analysis was conducted using forward selection (likeli
hood ratio) for the remaining 85 variables (Table S5).