Skip to main content
  • Research article
  • Open access
  • Published:

Endocannabinoids and related lipids in blood plasma following touch massage: a randomised, crossover study



The endocannabinoid system is involved in the regulation of stress and anxiety. In a recent study, it was reported that short-term changes in mood produced by a pleasant ambience were correlated with changes in the levels of plasma endocannabinoids and related N-acylethanolamines (Schrieks et al. PLoS One 10: e0126421, 2015). In the present study, we investigated whether stress reduction by touch massage (TM) affects blood plasma levels of endocannabinoids and related N-acylethanolamines.


A randomized two-session crossover design for 20 healthy participants was utilised, with one condition that consisted of TM and a rest condition as control. TM increased the perceived pleasantness rating of the participants, and both TM and rest reduced the basal anxiety level as assessed by the State scale of the STAI-Y inventory. However, there were no significant effects of either time (pre- vs. post-treatment measures) as main effect or the interaction time x treatment for the plasma levels of the endocannabinoids anandamide and 2-arachidonoylglycerol or for eight other related lipids. Four lipids showed acceptable relative reliabilities, and for two of these (linoleoyl ethanolamide and palmitoleoyl ethanolamide) a significant correlation was seen between the TM-related change in levels, calculated as (post-TM value minus pre-TM value) − (post-rest value minus pre-rest value), and the corresponding TM-related change in perceived pleasantness.


It is concluded that in the participants studied here, there are no overt effects of TM upon plasma endocannabinoid levels. Possible associations of related N-acylethanolamines with the perceived pleasantness should be investigated further.


For centuries, extracts of the Cannabis sativa plant have been used for their therapeutic and recreational properties. ∆9-Tetrahydrocannabinol, the main psychoactive ingredient of cannabis produces its effects via activation of two cannabinoid (CB) receptors, CB1 and CB2. CB1 receptors are the most abundant G protein-coupled receptor expressed in the brain but are also found in cells and tissues outside the central nervous system. CB2 receptors are mainly expressed in cells of the immune tissue but have also been identified in the CNS [1, 2]. The most well-characterised endogenous CB receptor ligands (endocannabinoids, eCBs) are the arachidonic acid derivatives anandamide (AEA) [3] and 2-arachidonglycerol (2-AG) [4, 5]. AEA belongs to the N-acylethanolamine (NAE) class of lipids, which includes other putative eCB ligands such as docosatetraenoyl ethanolamide (DEA) [6] as well as the anti-inflammatory compounds palmitoylethanolamide (PEA) and stearoylethanolamide (SEA) [7, 8], the satiety factor oleoylethanolamide (OEA) [9] and other less well-characterised compounds such as linoleoyl ethanolamide (LEA), palmitoleoyl ethanolamide (POEA) and eicosapentaenoyl ethanolamide (EPEA). NAEs are catabolised by hydrolysis to the corresponding long-chain fatty acids, but AEA can also be metabolised to other compounds including N-arachidonoylglycine (NA-Gly) [10], which has biological properties of its own [11].

The endocannabinoid (eCB) system is involved in a wide variety of effects in the body, including the inhibition of central transmitter release and the regulation of appetite and pain perception [12]. A substantial amount of research has linked the eCB system to stress, fear and anxiety. One hypothesised function for the eCB system is to buffer or dampen the endocrine and behavioural effects of acute stress and negative stimuli [13, 14] and in humans increased levels of eCB have been detected in serum after exposure to acute stress [15]. Both decreases [16] and increases [17] in plasma eCB levels have been reported for individuals with post-traumatic stress disorder. Loss of eCB tone can result in mood disturbances in man: the CB1 receptor antagonist/inverse agonist rimonabant was briefly marketed for weight reduction. The medication was effective, patients lost weight and experienced food less pleasurably. However, after a short time on the market, the drug was withdrawn due an unfavourable risk of developing depressive disorders, mood alterations with depressive symptoms, and anxiety [18].

In a recent study [19], it was reported that in healthy women, the eCB system was responsive to mood changes produced by the ambient environment. These authors utilised a randomised cross-over study whereby the plasma levels of AEA, 2-AG, PEA and SEA were measured during the partaking of a meal under conditions of either pleasant or unpleasant ambience. The authors found that levels of SEA and PEA were changed during the experience, and that the lipids were correlated with measures of mood [19]. This raises the question as to whether other acute mood-changing interventions can affect these lipids in a similar manner. Human touch is involved in many rewarding behaviours and is considered to be an inherited, gene specified reward [20]. Touch massage (TM) and human touch reduce stress and anxiety and increase well-being and relaxation in humans [2124]. In the present randomised cross-over study, we have investigated whether touch massage also affects the plasma levels of eCBs and related lipids.


Influence of the interventions upon behavioural and stress-related measures

Twenty healthy individuals participated in the study. The individuals were randomly divided into two groups, one group (N = 10) experiencing a rest phase (60 min) between blood sampling on the first occasion and 60 min of TM between blood sampling on the second occasion. The other group (N = 10) was treated in the reverse order to give a maximum of 20 cases in this randomised cross-over study (for schematic, showing the sampling times, see Fig. 1). The values for all rest and TM conditions were then collated. The initial stress levels (i.e. prior to blood sampling) on a 5 point scale (5 worst) were not significantly different between the rest conditions and the TM conditions, median (range) scores of 1.5 (0–5) and 2 (0–3.5) being found for rest and TM, respectively (P = 0.61, Wilcoxon’s matched pairs signed rank test). Similarly, pre-treatment saliva cortisol levels (nmol) were not different, median (range) scores of 8.8 (4.4–17.2) and 8.7 (3.4–18.6) being found for rest and TM, respectively (P = 0.77, Wilcoxon’s matched pairs signed rank test).

Fig. 1
figure 1

Flow chart of the study showing the sampling times

The effects of TM upon perceived pleasantness, the heart rate, the total scores in the state and trait arms of the State-Trait Anxiety Inventory Form Y (STAI-Y) and the Montgomery-Åsberg Depression Rating scale (MADRS) are shown in Fig. 2. Each panel shows the individual data points plotted with the x-axis representing the pre-treatment value, and the y-axis representing the post-treatment value. The correlations between the pre- and post-treatment values were very high for all variables (Additional file 1). However, the primary aim of the plots is to provide a simple visual representation of the individual data. In this respect, the dotted 45° lines of identity represent the situation where the intervention had not changed the behavioural scores at all: a treatment effect would result in a consistent deviation away from this 45° line of identity. Significance for treatment effect was assessed using two-way ANOVA for matched samples (Table 1). There was, as expected, a significant treatment × time (pre- vs. post-treatment) interaction for the perceived pleasantness scale (Fig. 2a; Table 1). The median value did not change following rest (median 0, range −2 to +1 units for post- minus pre-treatment values), whereas it increased by a median of 0.5 units (range −1 to +3 units) following TM. A reverse pattern was seen for the heart rate (Fig. 2b): the heart rate dropped by a median of 8 beats/min (range of decrease 1–27 beats/min) following the rest period, whereas the corresponding drop for TM was smaller (3 beats/min, range of decrease −14 to 12 beats/min). Measures of anxiety and depression were also obtained for the participants, with the exception of one case pre-TM, and the behavioural scores are shown graphically in Fig. 2c–e. For the total state scores in the STAI-Y inventory, a significant main effect of time (i.e. pre vs. post-treatment) was seen, but no significant effects of either time × treatment or treatment per se. This can be seen as a downward shift relative to the dotted line in Fig. 2c. Corresponding data for the Trait arm of the STAI-Y scale and the MADRS scale are shown Fig. 2d, e. Taken together, the main effect of TM over and above rest per se in this population was upon perceived pleasantness.

Fig. 2
figure 2

Effects of rest and TM upon a perceived pleasantness; b heart rate, c STAI-state; d STAI-trait and e MADRS scores of the subjects. The graphs show the individual data points, and the dotted lines are the 45° line of identity

Table 1 Statistical evaluation of the datasets shown in Figs. 2 and 3

Influence of the interventions upon blood levels of eCBs and related lipids

Plasma levels of 2-AG, AEA, the related NAEs PEA, SEA, POEA, OEA and LEA, and the related arachidonoyl derivative NA-Gly are shown in Fig. 3; the statistical treatment of the data are given in Table 1 and the pre-: post-treatment correlation coefficients are given in Additional file 1. One sample was lost for the rest post-treatment group and one for the TM post-treatment group. Thus, eighteen cases were evaluated on all four occasions. For only one compound did all four sample groups pass the D’Agostino and Pearson test for normality (NA-Gly), but this issue was rectified by use of logged values for 2-AG, PEA, SEA, OEA, LEA, POEA and square root values for AEA. Log10-transformed NA-Gly also passed the D’Agostino and Pearson test for normality. Two-way ANOVA for repeated measures conducted on the transformed data indicated no effect of time, treatment or the interaction time x treatment other than a significant effect of treatment (P = 0.029) for PEA, due to higher values for the rest occasions than on the TM occasions (Table 1). EPEA and DEA were also measured (Fig. 3h, i). These two compounds failed the D’Agostino and Pearson test for normality both for untransformed data and upon square root transformation. Log10 transformation was not undertaken due to the presence of zero (below the limit of detection) values. Two-way ANOVA for repeated measures on the rank-transformed values indicated that there were no effects of time, treatment or time × treatment on the levels of these two lipids (Table 1).

Fig. 3
figure 3

Effects of rest and TM upon the plasma levels of a 2-AG; b AEA; c PEA; d SEA; e OEA; f LEA; g POEA; h EPEA; i DEA and j NA-Gly (N = 18). The graphs show the individual data points, and the dotted lines are the 45° line of identity

The availability of both biochemical and behavioural data allows analysis of associations to be made using non-parametric Spearman’s correlations. However, use of such a measure assumes that the relative reliability, defined by Baumgartner [25] as the “indication of the degree to which people maintain their position within a group” is acceptable. Here, the definition means the position of the lipid concentration for a given individual within the range seen for the whole sample.

We assessed the relative reliability of the lipids by comparing their pre-treatment values on the first visit (when they were blinded to the forthcoming treatment) and the second visit. The standard methods used to assess relative reliability are intraclass correlation coefficients (ICCs) and Lin’s concordance coefficient (Lin’s CCC), for which the sample sizes here are sufficient [2528]. These tests are parametric, and so the transformed data were used (Table 2). Fleiss [28] suggested that ICC values <0.4, 0.4–0.75 and >0.75 can be taken to represent poor, fair to good, and excellent reliability, respectively. Using this measure, the data suggest that the measurements for LEA, OEA, POEA and NA-Gly have sufficient relative reliability for investigation of potential associations with the behavioural variables.

Table 2 Relative reliability of measurements of endocannabinoids and six related lipids in human plasma samples

In Table 3, Spearman’s rho values for the correlation between the post- and pre-values for LEA, OEA, POEA and NA-Gly, and the corresponding changes for either the state scores or the heart rates, are presented. These were chosen since they are associated with anxiety and stress, and show the largest changes post- : pre-treatment of the variables investigated. None of the zero order correlations reached significance. Table 2 also shows the correlations for “TM-selective” (sTM) change in lipid concentrations vs. the sTM for perceived pleasantness. The sTM values were calculated as the change over time for the TM group minus the change for the rest group, thereby reducing potential influences of time, food intake etc. [19, 2933]. In this case, significant negative correlations for LEA and POEA were found (Table 3).

Table 3 Correlations between plasma lipid, state, heart rate and perceived pleasantness scores


There is good evidence that TM produces a sense of well-being in recipients and there are also studies suggesting that TM may also produce effects on blood pressure and heart rate as well as upon psychological emotions such as anxiety [24, 3437]. An advantage of the cross-over design used here is that the individuals act as their own controls, and the first author of the present study has previously used this model, with a similar sample size (22 individuals) to demonstrate changes in heart rate and heart rate variability following TM [36]. In the present study, an effect on the perceived pleasantness rating scale for TM was observed. This is consistent with a previous study from the first author, who showed that the increased level of perceived pleasantness was associated with an increased blood flow in the pregenual anterior cingulate cortex [24], a brain region that is activated by pleasurable stimuli [38]. In contrast, effects upon the state arm of the STAI-Y scale over and above that seen for rest was not seen following TM in the present study. The most likely explanation of the lack of effect on this measure is that a floor effect is occurring, given that the volunteers were not overly anxious. Certainly, the TM protocol used here produces a notable reduction in the state scores in patients scheduled for elective aortic surgery [37].

With respect to the biochemical measurements, no significant TM-induced changes, which would be apparent as an interaction time × treatment, were seen at the group level for plasma concentrations of 2-AG, AEA or related lipids following either TM or rest. To our knowledge, the only previous study investigating the response to the eCB system to a hands-on technique is that of McPartland et al. [39], who investigated osteopathic manipulative treatment in non-naïve individuals. They did not find any significant change at the 5 % level in post-treatment serum levels of AEA, 2-AG or OEA compared to pre-treatment levels following either active or sham treatment. The association of the eCB system with stress and anxiety is complex, and a synthesis of animal experimental data has suggested that the activity of this system needs to be kept within a narrow window [40]. In humans, blockade of CB1 receptors by rimonabant is associated with an unacceptable incidence of depressive disorders and mood alterations with depressive symptoms and anxiety [18]. Rimonabant given acutely does not affect basal anxiety levels, but increases anxiety under aversive conditions (simulated public speaking) [41]. Conversely, in animal studies, potentiation of the eCB system by compounds inhibiting eCB catabolism produce potentially beneficial effects in models of anxiety using conditions of high, but not low, aversiveness [4244]. These data may suggest that the eCB system is only engaged once a threshold of anxiety is reached. The type of stress may also be important. Thus, the stress produced by the Trier Social Stress Task (a series of short stressful events including a mock job interview) is sufficient to produce increases in blood levels of AEA, although this may reflect a general increase in NAE levels, since PEA and OEA levels were increased, whereas 2-AG and the related monoacylglycerol 2-oleoylglycerol levels were not changed [15]. Levels of these NAEs were not changed during unstressed conditions [15]. In healthy males participating for the first time in repeat parabolic flights, individuals responding with motion sickness showed an increased stress score, and this was accompanied by a reduction in AEA compared to individuals who did not feel motion sickness after the 10th parabola [45]. These examples are of fairly strong stressors, whereas the individuals in our studywere not unduly stressed or anxious. A lack of effect upon blood eCB levels is thus consistent with the “threshold of anxiety” hypothesis described above.

With respect to the related NAEs, the situation may be different: Schrieks et al. [19] reported that in their randomised cross-over study there were state-dependent changes in PEA and SEA levels in their study investigating the effects of ambience upon plasma levels of these lipids. These authors followed the levels of AEA, 2-AG, PEA, SEA, OEA and docohexanoylethanolamine (DHEA) in plasma at baseline, and 30 and 120 min following the partaking of a meal with or without alcohol and in either a pleasant or an unpleasant atmosphere, created by environmental cues (lighting, cleanliness, choice of music etc.). The authors found a significant interaction effect of time × ambience, so that at the 30 min time point, both PEA and SEA levels (expressed as change from baseline) were lower in the pleasant ambience arms of the study than in the unpleasant ambience arms (both with or without alcohol for SEA, with alcohol only for PEA) whereas the changes in 2-AG and AEA levels were not. It would certainly be worth designing future studies to investigate the effects of TM upon plasma eCB and NAE levels in individuals with higher stress/anxiety levels.

The present study has also investigated whether changes in the concentrations of LEA, OEA, POEA and NA-Gly were associated with changes in the state scores, the heart rates and the “TM-selective” change in perceived pleasantness. At first sight, the results suggest that this may be so, a result in line with the study of Schrieks et al. [19], who found significant correlations with changes between baseline and 120 min of their lipid measurements and mood scores in their questionnaire for the pleasant ambiance without alcohol arm of study. Mood scores in that study were assessed using the Profile of Mood States questionnaire with additional questions from the Brunel mood scale. All in all, seven measures (anger, depression, fatigue, tension, vigour, happiness and calmness) were quantitated and significant positive correlations were found between ∆AEA as well as ∆DHEA and ∆fatigue, and between ∆SEA and ∆calmness, whilst negative correlations were found between ∆2-AG, ∆PEA and ∆happiness, and for ∆OEA and ∆vigour [19]. However, it is wise to remember that the significant correlations from both the present study and that of Schrieks et al. [19] result from multiple testing. In the case of [19] the total number of correlation coefficients that can be obtained from the pleasant ambience without alcohol arm of the study is presumably 42 (6 lipids, seven behavioural scores), and in our case there were 12 zero-order correlation coefficients. There is a trade-off between running the risk of reporting Type I errors (since some of the correlation coefficients would be expected to be significant on a purely random basis) and of reporting Type II errors when conservative corrections, such as the Bonferroni correction, are used. It can be argued that in exploratory studies of this type, the Bonferroni correction is too draconian (for [19] the presumed 42 comparisons for the arm in question would, upon implementation of the Bonferroni correction, reduce the minimum level of significance to 0.0012, which would leave but a single correlation as significant; the same would be true for the less conservative false discovery rate method of Benjamini and Hochberg [46]). Thus, a reasonable conclusion to make is that the significant correlations found in that study and here are consistent with the hypothesis that the circulating levels of endocannabinoids and related NAEs are responsive to environmental cues, but confirmation of these exploratory findings is necessary.


It is concluded that in the participants studied here, there are no overt effects of TM upon plasma endocannabinoid levels. Possible associations of related N-acylethanolamines with the perceived pleasantness should be investigated further.


Ethics statement

This study was approved by the Regional Ethical Review Board in Umeå (Dnr 2014-133-32 M) and the procedures were conducted in accordance with the Declaration of Helsinki. All subjects were informed of the purpose and risks of the experiment before the study and gave their written consent. The study was carried out from September 2014 to February 2015.


Twenty healthy participants [10 men and 10 women, mean age = 25.2 years; standard deviation (SD) = 4.3 and mean body mass index = 23.8; SD = 3.1 at the start of the study] were recruited via posters. Nineteen of the participants were students and one participant worked at the University. Participants were excluded if they were smokers, took medications, or had a history of drug abuse.


The study had a randomized two-session crossover design with one condition that consisted of TM and a rest condition as control. A sequential number was assigned to an opaque envelope, which contained the assigned treatment group for the participant (half of the participants started with the rest condition and the other half started with the TM condition). The subjects participated twice with a minimum of 2 days and maximum of 56 days apart (median value 14 days; the data was positively skewed). Participants were instructed to have the same breakfast on both occasions and not to take any caffeine in the morning when the study occurred.


The study took place at the same settings and at the same time each day (see Fig. 1 for a flow chart). The participants were prepared with a venous catheter (PVC) at 8.30 a.m. and rested in supine position for 10 min before the blood samples (“pre-measures”) were collected into BD Vacutainer® acrylic based polymer gel tubes containing spray dried lithium heparin. At the same time blood pressure and heart rate were measured. The participants completed questionnaires (see below) before the rest/TM condition, starting at 9.15 a.m. During the “rest” session, the participants rested for 60 min. During the TM session, the participants received TM on legs, back, head, arms and hands for 60 min. The body parts not being massaged were covered with a blanket. The same person that conducted the TM was also present in the room for the rest session.

At 10.20 a.m., blood samples were collected, and blood pressure and heart rate were measured, after which the participants again filled in the questionnaires. All blood samples were centrifuged (15 min at RT) within 10 min of collection, and the plasma was then kept on dry ice until storage at −80 °C. Thus, the pre-treatment samples were placed on dry ice at ~9:10 a.m., the post-treatment samples at ~10:40 a.m. and then delivered to the analysis unit at 11–11:30 a.m. where they were rapidly frozen at −80 °C. Under these conditions, degradation of eCBs has been reported to be minimal, although extended times at room temperature are not recommended [47].


TM is a calm gentle massage with slow (0–5 cm/min) stroking movement at a pressure of about 2.5 N [24, 37]. TM is similar to tactile massage, gentle or soft massage and includes effleurage. There are several different massage techniques, some involving the muscles, others involving the skin. In Sweden, the term TM is used to cover different forms of light massage involving the skin. These forms of massage are characterized by gentle touches of the skin involving light pressure effleurage and long, calm stroking movements intended to increase the well-being [48], pleasure [49] and relaxation [35] of the recipients.

Subjective measures

To assess subjective effects upon perceived pleasantness and behaviours associated with anxiety, stress and depression, several different instruments were used. Perceived pleasantness was self assessed by the participants by rating how they experienced the moving touch stimulation on a visual analogue scale ranging from −5 (very unpleasant), 0 neutral, and +5 (very pleasant). This scale has previously been used for evaluation of pleasant touch stimulation [24, 37].

The Swedish version of the STAI-Y instrument was used to assess self-reported anxiety levels pre and post intervention/control. The STAY-Y was designed to measure State-Anxiety as a temporary condition of anxiety and Trait-Anxiety as a more underlying condition. The instrument has good construct validity, internal consistency, and test–retest validity [50, 51]. The 20 items are rated on a 4-point scale and after reversing positive items, the total score ranges from 20 to 80, with higher scores reflecting higher levels of anxiety. The STAI-Y instrument is sensitive enough to evaluate anxiety in small groups.

The participants subjective stress levels were rated on a visual analogue scale ranging from 0 = no perceived stress to 5 = very high perceived stress. MADRS was used for subjective experience of depression [52]. The MADRS consists of 9 items rated on a 7-point scale (0–6), which are summated to give a total score.

Heart rate was measured using an automatic blood pressure monitor, M6 Comfort, Omron Healthcare Europe, Netherlands.

Saliva cortisol

Salivary samples were collected by chewing a cylindrical cotton swab (Salivette, Sarstedt, Numbrecht, Germany). The samples were frozen at −20 °C until assayed. High performance liquid chromatography (HPLC) (Shimadzu Nexera) coupled to liquid chromatography tandem mass spectrometry (LC–MS/MS) (Sciex Qtrap 5500), was used for the analyses. The analyses were conducted in the Clinical Chemical Laboratory at the University Hospital, Umeå, Sweden.

Chemical and reagents for the lipid analyses

The following native and deuterated standards were purchased from Cayman Chemicals (Ann Arbor, MI, USA): 2-AG, AEA, PEA, SEA, OEA, LEA, POEA, EPEA, DEA, NA-Gly, 2-AG-d8, AEA-d8, PEA-d4, SEA-d3, OEA-d4, 12-[[(cyclohexylamino)carbonyl]amino]-dodecanoic acid (CUDA). Acetonitrile and methanol were from Merck (Darmstadt, Germany). Isopropanol was from VWR PROLABO (Fontenay-sous-Bois, France). Acetic acid was purchased from Aldrich Chemical Company, Inc. (Milwaukee, WI, USA). Butylhydroxytoluene (BHT) was from Cayman Chemical (Ann Arbor, MI, USA) and ethylenediaminetetraacetic acid (EDTA) from Fluka Analytical, Sigma-Aldrich (Buchs, Switzerland). Glycerol was from Fischer Scientific (Loughborough, UK). All solvents and chemicals were of HPLC grade or higher. Water was purified by a Milli-Q Gradient system (Millipore, Milford, MA, USA).

Preparation of native lipid standards and deuterated internal standard curves

Analytical quantification standards were used as ready-made standard stock solutions or as solutions prepared from solid substances and stored at −80 °C. 2-AG was prepared and stored in acetonitrile and the other standards were prepared and stored in ethanol. Final stock solution concentrations of native standards were: 250 µg/mL (2-AG, LEA and POEA), 125 µg/mL (AEA, PEA, OEA, DEA, EPEA, NA-Gly) and 83 µg/mL (SEA). Stock solutions of deuterated internal standards were prepared to a final concentration of 40 µg/mL (2-AG-d8) and 10 µg/mL (AEA-d4, PEA-d4, SEA-d3, OEA-d4).

For each native compound, a suitable internal standard was selected based on structural similarities and added to samples before extraction. The internal standards were: 2-AG-d8, AEA-d4, SEA-d3, OEA-d4 (used for POEA, EPEA and NA-Gly as well as OEA) and PEA-d4 (used for LEA and DEA as well as PEA). Recovery rates of each internal standard were calculated by normalisation against the recovery standard CUDA added in the last step before injection.

Sample preparation for lipid analysis

Sample preparation was undertaken according to our previously validated method [53]. Briefly, plasma samples (500 µL) were thawed on ice and spiked with 20 µL internal standard solution (800 ng/mL 2-AG-d8, and 20 ng/mL AEA-d4, OEA-d4, SEA-d3 and PEA-d4) and 10 μL antioxidant solution [0.2 mg/mL BHT/EDTA in methanol/water (1:1)] and then applied to pre-washed and conditioned solid phase extraction SPE cartridges (Waters Oasis HLB cartridges; 60 mg of sorbent, 30 µm particle size). After loading the sample containing internal standard and antioxidant solution, the cartridges were washed, dried under high vacuum for about 1 min, and eluted with 3 mL acetonitrile, followed by 2 mL of methanol and 1 mL of ethyl acetate into polypropylene tubes containing 6 μL of a glycerol solution (30 % in methanol). Eluates were concentrated with a MiniVac system (Farmingdale, NY, USA), reconstituted in 100 μL of methanol and vortexed (if necessary centrifuged to remove any residuals). Solutions were then transferred to LC vials with low-volume inserts, 10 μL of CUDA (50 ng/mL) was added, and LC–MS/MS analysis was performed immediately.

Instrumentation used for chromatographic separation and detection of lipids

Lipid analysis was performed using an Agilent UPLC system (Infinity 1290) coupled with an electrospray ionization source (ESI) to an Agilent 6490 Triple Quadrupole system equipped with the iFunnel Technology (Agilent Technologies, Santa Clara, CA, USA) operating in positive mode [54]. Analyte separation was performed using a Waters BEH C18 column (2.1 mm × 150 mm, 1.7 µm particle size) with 10 µL injection volume. The mobile phase consisted of (A) 0.1 % acetic acid in MilliQ water and (B) acetonitrile : isopropanol (90:10). The following gradient was employed: 0.0–2.0 min 30–45 % B, 2.0–2.5 min 45–79 % B, 2.5–11.5 min 79 % B, 11.5–12 min 79–90 % B, 12–14 min 90 % B, 14–14.5 min 90–79 % B, 14.5–15.5 min 79 % B, 15.6–19 min 30 % B.

Precursor ions, [M+H]+ product ions, multiple reaction monitoring (MRM) transitions and optimal collision energies were optimized for each analyte. ESI conditions were: capillary and nozzle voltage at 4000 and 1500 V, drying gas temperature 230 °C with a gas flow of 15 L/min, sheath gas temperature 400 °C with a gas flow of 11 L/min. The nebulizer gas flow was 35 psi, and iFunnel high and low pressure RF at 150 and 60 V. The dynamic MRM option was used for all compounds with optimized transitions and collision energies. The MassHunter Workstation software was used to control acquisition and manually integrate all peaks. The intra- and interday precision and accuracy measures of the lipids using this QqQ system are reported in detail elsewhere [54].


Wilcoxon’s signed ranks tests, paired signed ranks tests, non-parametric bivariate correlations, D’Agostino and Pearson normality tests were undertaken using v. 6 for the Macintosh (GraphPad Software Inc., San Diego, CA, USA). Intraclass correlation coefficients and two-way repeated- measures ANOVA were obtained using SPSS software v.22 for the Macintosh (IBM). Lin’s concordance correlation coefficients were calculated online using the facility provided by the National Institute of Water and Atmospheric Research (; URL checked 17 August 2015). The confidence limits for the differences between non-overlapping correlations (Additional file 1) were calculated using a Microsoft Spreadsheet provided by Dr. Thom Baguley at, URL checked 17 August 2015).

Availability of supporting data

The data sets supporting the results of this article are included within the article.











12-[[(cyclohexylamino)carbonyl]amino]-dodecanoic acid


docosatetraenoyl ethanolamide




eicosapentaenoyl ethanolamide


electrospray ionization source


high performance liquid chromatography


intraclass correlation coefficient


liquid chromatography tandem mass spectrometry


linoleoyl ethanolamide

Lin’s CCC:

Lin’s concordance coefficient


multiple reaction monitoring












square root transformation


State-trait anxiety inventory form Y




touch massage


  1. Mackie K. Distribution of cannabinoid receptors in the central and peripheral nervous system. Handb Exp Pharmacol. 2005;168:299–325.

    Article  CAS  PubMed  Google Scholar 

  2. Atwood BK, Mackie K. CB2: a cannabinoid receptor with an identity crisis. Br J Pharmacol. 2010;160:467–79.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  3. Devane W, Hanus L, Breuer A, Pertwee R, Stevenson L, et al. Isolation and structure of a brain constituent that binds to the cannabinoid receptor. Science. 1992;258:1946–9.

    Article  CAS  PubMed  Google Scholar 

  4. Sugiura T, Kondo S, Sugukawa A, Nakane S, Shinoda A, et al. 2-Arachidonylglycerol: a possible endogenous cannabinoid receptor ligand in brain. Biochem Biophys Res Commun. 1995;215:89–97.

    Article  CAS  PubMed  Google Scholar 

  5. Mechoulam R, Ben-Shabat S, Hanus L, Ligumsky M, Kaminski N, et al. Identification of an endogenous 2-monoglyceride, present in canine gut, that binds to cannabinoid receptors. Biochem Pharmacol. 1995;50:83–90.

    Article  CAS  PubMed  Google Scholar 

  6. Hanuš L, Gopher A, Almog S, Mechoulam R. Two new unsaturated fatty acid ethanolamides in brain that bind to the cannabinoid receptor. J Med Chem. 1993;36:3032–4.

    Article  PubMed  Google Scholar 

  7. Lambert DM, Vandevoorde S, Jonsson K-O, Fowler CJ. The palmitoylethanolamide family: a new class of anti-inflammatory agents? Curr Med Chem. 2002;9:663–74.

    Article  CAS  PubMed  Google Scholar 

  8. Dalle Carbonare M, Del Giudice E, Stecca A, Colavito D, Fabris M, et al. A saturated N-acylethanolamine other than N-palmitoyl ethanolamine with anti-inflammatory properties: a neglected story. J Neuroendocrinol. 2008;20(Suppl 1):26–34.

    Article  CAS  PubMed  Google Scholar 

  9. Rodríguez de Fonseca F, Navarro M, Gómez R, Escuredo L, Nava F, et al. An anorexic lipid mediator regulated by feeding. Nature. 2001;414:209–12.

    Article  PubMed  Google Scholar 

  10. Bradshaw HB, Rimmerman N, Hu SS-J, Benton VM, Stuart JM, et al. The endocannabinoid anandamide is a precursor for the signaling lipid N-arachidonoyl glycine by two distinct pathways. BMC Biochem. 2009;10:14.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Huang S, Bisogno T, Petros T, Chang S, Zavitsanos P, et al. Identification of a new class of molecules, the arachidonoyl amino acids, and characterization of one member that inhibits pain. J Biol Chem. 2001;276:42639–44.

    Article  CAS  PubMed  Google Scholar 

  12. Pacher P, Bátkai S, Kunos G. The endocannabinoid system as an emerging target of pharmacotherapy. Pharmacol Rev. 2006;58:389–462.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Hillard CJ. Stress regulates endocannabinoid-CB1 receptor signaling. Semin Immunol. 2014;26:380–8.

    Article  CAS  PubMed  Google Scholar 

  14. McLaughlin RJ, Hill MN, Gorzalka BB. A critical role for prefrontocortical endocannabinoid signaling in the regulation of stress and emotional behavior. Neurosci Biobehav Rev. 2014;42:116–31.

    Article  CAS  PubMed  Google Scholar 

  15. Dlugos A, Childs E, Stuhr KL, Hillard CJ, de Wit H. Acute stress increases circulating anandamide and other N-acylethanolamines in healthy humans. Neuropsychopharmacology. 2012;37:2416–27.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  16. Hill MN, Bierer LM, Makotkine I, Golier JA, Galea S, et al. Reductions in circulating endocannabinoid levels in individuals with post-traumatic stress disorder following exposure to the world trade center attacks. Psychoneuroendocrinology. 2013;38:2952–61.

    Article  CAS  PubMed  Google Scholar 

  17. Hauer D, Schelling G, Gola H, Campolongo P, Morath J, et al. Plasma concentrations of endocannabinoids and related primary fatty acid amides in patients with post- traumatic stress disorder. PLoS One. 2013;8:e62741.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  18. Van Gaal L, Pi-Sunyer X, Després J-P, McCarthy C, Scheen A. Efficacy and safety of rimonabant for improvement of multiple cardiometabolic risk factors in overweight/obese patients: pooled 1-year data from the Rimonabant in Obesity (RIO) program. Diabetes Care. 2008;31(Suppl 2):S229–40.

    Article  PubMed  Google Scholar 

  19. Schrieks IC, Ripken D, Stafleu A, Witkamp RF, Hendriks HFJ. Effects of mood inductions by meal ambiance and moderate alcohol consumption on endocannabinoids and N-acylethanolamines in humans: a randomized crossover trial. PLoS One. 2015;10:e0126401.

    Article  Google Scholar 

  20. Moyer CA, Rounds J, Hannum JW. A meta-analysis of massage therapy research. Psychol Bull. 2004;130:3–18.

    Article  PubMed  Google Scholar 

  21. Ågren A, Berg M. Tactile massage and severe nausea and vomiting during pregnancy—women’s experiences. Scand J Caring Sci. 2006;20:169–76.

    Article  PubMed  Google Scholar 

  22. Coan JA, Schaefer HS, Davidson RJ. Lending a hand: social regulation of the neural response to threat. Psychol Sci. 2006;17:1032–9.

    Article  PubMed  Google Scholar 

  23. Harris M, Richards KC. The physiological and psychological effects of slow-stroke back massage and hand massage on relaxation in older people. J Clin Nurs. 2010;19:917–26.

    Article  PubMed  Google Scholar 

  24. Lindgren L, Westling G, Brulin C, Lehtipalo S, Andersson M, et al. Pleasant human touch is represented in pregenual anterior cingulate cortex. Neuroimage. 2012;59:3427–32.

    Article  PubMed  Google Scholar 

  25. Baumgartner TA. Norm-referenced measurement: reliability. In: Safarit MJ, Wood TM, editors. Measurement concepts in physical education and exercise science. Champaign; 1989. pp 45–72.

  26. Lin LI-K. A concordance correlation coefficient to evaluate reproducibility. Biometrics. 1989;45:255–68.

    Article  CAS  PubMed  Google Scholar 

  27. Shrout PE, Fleiss JL. Intraclass correlations: uses in assessing rater reliability. Psychol Bull. 1979;86:420–8.

    Article  CAS  PubMed  Google Scholar 

  28. Fleiss JL. Reliability of measurement. In: Fleiss JL, editor. Design and analysis of clinical experiments. New York: Wiley; 1986. p. 1–32.

    Google Scholar 

  29. Hanlon EC, Tasali E, Leproult R, Stuhr KL, Doncheck E, et al. Circadian rhythm of circulating levels of the endocannabinoid 2-arachidonoylglycerol. J Clin Endocrinol Metab. 2015;100:220–6.

    Article  CAS  PubMed  Google Scholar 

  30. Jones PJH, Lin L, Gillingham LG, Yang H, Omar JM. Modulation of plasma N-acylethanolamine levels and physiological parameters by dietary fatty acid composition in humans. J Lipid Res. 2014;55:2655–64.

    Article  CAS  PubMed  Google Scholar 

  31. Monteleone P, Piscitelli F, Scognamiglio P, Monteleone AM, Canestrelli B, et al. Hedonic eating is associated with increased peripheral levels of ghrelin and the endocannabinoid 2-arachidonoyl-glycerol in healthy humans: a pilot study. J Clin Endocrinol Metab. 2012;97:E917–24.

    Article  CAS  PubMed  Google Scholar 

  32. Joosten MM, Balvers MGJ, Verhoeckx KCM, Hendriks HFJ, Witkamp RF. Plasma anandamide and other N-acylethanolamines are correlated with their corresponding free fatty acid levels under both fasting and non-fasting conditions in women. Nutr Metab (Lond). 2010;7:49.

    Article  PubMed Central  PubMed  Google Scholar 

  33. Hansen HS. Effect of diet on tissue levels of palmitoylethanolamide. CNS Neurol Disord Drug Targets. 2013;12:17–25.

    Article  CAS  PubMed  Google Scholar 

  34. Garner B, Phillips LJ, Schmidt HM, Markulev C, O’Connor J, et al. Pilot study evaluating the effect of massage therapy on stress, anxiety and aggression in a young adult psychiatric inpatient unit. Aust N Z J Psychiatry. 2008;42:414–22.

    Article  PubMed  Google Scholar 

  35. Billhult A, Määttä S. Light pressure massage for patients with severe anxiety. Complement Ther Clin Pract. 2009;15:96–101.

    Article  PubMed  Google Scholar 

  36. Lindgren L, Rundgren S, Winsö O, Lehtipalo S, Wiklund U, et al. Physiological responses to touch massage in healthy volunteers. Autonom Neurosci Nas Clin. 2010;158:105–10.

    Article  CAS  Google Scholar 

  37. Lindgren L, Lehtipalo S, Winsö O, Karlsson M, Wiklund U, et al. Touch massage: a pilot study of a complex intervention. Nurs Crit Care. 2013;18:269–77.

    Article  PubMed  Google Scholar 

  38. Rolls ET, O’Doherty J, Kringelbach ML, Francis S, Bowtell R, et al. Representations of pleasant and painful touch in the human orbitofrontal and cingulate cortices. Cereb Cortex. 2003;13:308–17.

    Article  CAS  PubMed  Google Scholar 

  39. McPartland JM, Giuffrida A, King J, Skinner E, Scotter J, et al. Cannabimimetic effects of osteopathic manipulative treatment. J Am Osteopath Assoc. 2005;105:283–91.

    PubMed  Google Scholar 

  40. Moreira FA, Lutz B. The endocannabinoid system: emotion, learning and addiction. Addict Biol. 2008;13:196–212.

    Article  CAS  PubMed  Google Scholar 

  41. Bergamaschi MM, Queiroz RHC, Chagas MHN, Linares IMP, Arrais KC, et al. Rimonabant effects on anxiety induced by simulated public speaking in healthy humans: a preliminary report. Hum Psychopharmacol. 2014;29:94–9.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  42. Naidu P, Varvel S, Ahn K, Cravatt B, Martin B, et al. Evaluation of fatty acid amide hydrolase inhibition in murine models of emotionality. Psychopharmacology. 2007;192:61–70.

    Article  CAS  PubMed  Google Scholar 

  43. Haller J, Barna I, Barsvari B, Gyimesi Pelczer K, Yasar S, et al. Interactions between environmental aversiveness and the anxiolytic effects of enhanced cannabinoid signaling by FAAH inhibition in rats. Psychopharmacology. 2009;204:607–16.

    Article  CAS  PubMed  Google Scholar 

  44. Sciolino NR, Zhou W, Hohmann AG. Enhancement of endocannabinoid signaling with JZL184, an inhibitor of the 2-arachidonoylglycerol hydrolyzing enzyme monoacylglycerol lipase, produces anxiolytic effects under conditions of high environmental aversiveness in rats. Pharmacol Res. 2011;64:226–34.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  45. Choukèr A, Kaufmann I, Kreth S, Hauer D, Feuerecker M, et al. Motion sickness, stress and the endocannabinoid system. PLoS One. 2010;5:e10752.

    Article  PubMed Central  PubMed  Google Scholar 

  46. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57:289–300.

    Google Scholar 

  47. Jian W, Edom R, Weng N, Zannikos P, Zhang Z, et al. Validation and application of an LC-MS/MS method for quantitation of three fatty acid ethanolamides as biomarkers for fatty acid hydrolase inhibition in human plasma. J Chromatogr B Analyt Technol Biomed Life Sci. 2010;878:1687–99.

    Article  CAS  PubMed  Google Scholar 

  48. Cronfalk BS, Strang P, Ternestedt B-M, Friedrichsen M. The existential experiences of receiving soft tissue massage in palliative home care—an intervention. Support Care Cancer. 2009;17:1203–11.

    Article  PubMed  Google Scholar 

  49. Bergsten U, Petersson IF, Arvidsson B. Perception of tactile massage as a complement to other forms of pain relief in rheumatic disease. Musculoskeletal Care. 2005;3:157–67.

    Article  PubMed  Google Scholar 

  50. Spielberger CD. Anxiety: current trends in theory and research, vol. 2. New York: Academic Press; 1972.

    Google Scholar 

  51. Spielberger CD, Gorssuch RL, Lushene PR, Vagg PR, Jacobs GA. Manual for the state-trait anxiety inventory. Palo Alto: Consulting Psychologists Press; 1983.

    Google Scholar 

  52. Montgomery SA, Åsberg M. A new depression scale designed to be sensitive to change. Br J Psychiat. 1979;134:382–9.

    Article  CAS  Google Scholar 

  53. Gouveia-Figueira S, Nording ML. Development and validation of a sensitive UPLC-ESI-MS/MS method for the simultaneous quantification of 15 endocannabinoids and related compounds in milk and other biofluids. Anal Chem. 2014;86:1186–95.

    Article  CAS  PubMed  Google Scholar 

  54. Gouveia-Figueira S, Nording ML. Validation of a tandem mass spectrometry method using combined extraction of 42 oxylipins and 15 endocannabinoid-related compounds including prostamides from biological matrices. Prostaglandins Other Lipid Med. 2015. doi:10.1016/j.prostaglandins.2015.06.003.

    Google Scholar 

  55. Conover WJ, Iman RL. Rank transformations as a bridge between parametric and nonparametric statistics. Am Statistician. 1981;35(3):124–9.

    Google Scholar 

  56. Lehmann R. General derivation of partial and multiple rank correlation coefficients. Biom J. 1977;19(4):229–36.

    Article  Google Scholar 

  57. Nieuwenhuis S, Forstmann BU, Wagenmakers E-J. Erroneous analyses of interactions in neuroscience: a problem of significance. Nat Neurosci. 2011;14:1105–7.

    Article  CAS  PubMed  Google Scholar 

  58. Zou GY. Toward using confidence intervals to compare correlations. Psychol Meth. 2007;12:399–413.

    Article  Google Scholar 

Download references

Authors’ contributions

LL was involved in the conception of the study, undertook the TM and blood sampling and contributed to the writing of the paper. SG-F and MN were responsible for the analyses of the endocannabinoids and related lipids contributed to the writing of the paper. CJF was involved in the conception of the study, the analyses and the writing of the paper. All authors read and approved the final manuscript.


The authors thank the Swedish Research Council (Grant no. 12158, medicine, to CJF); the Swedish Research Council Formas (to MN); Avtal om Läkarutbildning och Forskning (Agreement on medical training and research, to LL), Västerbotten County Council (LL); Umeå University and the Research Funds of the Medical Faculty (CJF) for financial support. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Compliance with ethical guidelines

Competing interests The authors declare that they have no competing interests.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Christopher J. Fowler.

Additional files


Additional file 1. Correlation coefficients for the pre- vs. post-treatment data summarised in Figs. 2 and 3.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lindgren, L., Gouveia-Figueira, S., Nording, M.L. et al. Endocannabinoids and related lipids in blood plasma following touch massage: a randomised, crossover study. BMC Res Notes 8, 504 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: