Plasma Levels of Endocannabinoids and Their Analogues Are Related to Specific Fecal Bacterial Genera in Young Adults: Role in Gut Barrier Integrity
PROFITH (PROmoting FITness and Health through Physical Activity) Research Group, Sport and Health University Research Institute (iMUDS), University of Granada, 18071 Granada, Spain; lortizalvarez7@ugr.es (L.O.-A.); elihuixu@gmail.com (H.X.); fj.osuna.prieto@gmail.com (F.J.O.-P.); franciscoacosta838@gmail.com (F.M.A.); b.martinez-tellez@lumc.nl (B.M.-T.)
Department of Biochemistry and Molecular Biology II, School of Pharmacy, University of Granada, 18071 Granada, Spain; jrplaza@ugr.es (J.P.-D.); agil@ugr.es (A.G.)
Leiden Academic Centre for Drug Research, Division of Systems Biomedicine and Pharmacology, Leiden University, 2300 Leiden, The Netherlands; x.di@lacdr.leidenuniv.nl
Division of BioAnalytical Chemistry, Institute of Molecular and Life Sciences (AIMMS), Vrije Universiteit Amsterdam, 1081 Amsterdam, The Netherlands; i.kohler@vu.nl
Center for Analytical Sciences Amsterdam, 1098 Amsterdam, The Netherlands
Center for Biomedical Research, Department of Analytical Chemistry, Institute of Nutrition and Food Technology, University of Granada, 18071 Granada, Spain
Research and Development of Functional Food Center (CIDAF), Health Sciences Technology Park, 18071 Granada, Spain
Turku PET Centre, University of Turku, 20014 Turku, Finland
Turku PET Centre, Turku University Hospital, 20521 Turku, Finland
InFLAMES Research Flagship Centre, University of Turku, 20014 Turku, Finland
Department of Gastroenterology, Hepatology and Infectious Diseases, Otto von Guericke University Magdeburg, 39106 Magdeburg, Germany; ramiro.vilchez@med.ovgu.de (R.V.-V.); alexander.link@med.ovgu.de (A.L.)
Children’s Hospital of Eastern Ontario Research Institute, Ottawa, ON K1H 8L1, Canada
Department of Molecular Physiology, Leiden Institute of Chemistry, Leiden University, 2300 Leiden, The Netherlands; m.van.der.stelt@chem.leidenuniv.nl
Leiden Academic Centre for Drug Research (LACDR), Department of Systems Biomedicine and Pharmacology, Leiden University, 2300 Leiden, The Netherlands; hankemeier@lacdr.leidenuniv.nl
Department of Cell Biology, Physiology and Immunology, Maimónides Biomedical Research Institute of Córdoba (IMIBIC), Reina Sofia University Hospital, University of Córdoba, 14004 Córdoba, Spain; mer.cp@hotmail.com
Unidad de Gestión Clínica Endocrinología y Nutrición, Instituto de Investigación Biomédica de Málaga-IBIMA, Hospital Universitario Virgen de la Victoria, Universidad de Málaga, 29016 Malaga, Spain; fjtinahones@uma.es
Centro de Investigación Biomédica En Red (CIBER), Fisiopatología de la Obesidad y Nutrición (CIBEROBN), Instituto de Salud Carlos III (ISCIII), 28029 Malaga, Spain
Biomedical Research Center, Institute of Nutrition and Food Technology “José Mataix”, University of Granada, Parque Tecnológico Ciencias de la Salud, 18071 Granada, Spain
Instituto de Investigación Biosanitaria, 18014 Granada, Spain
Einthoven Laboratory for Experimental Vascular Medicine, Department of Medicine, Division of Endocrinology, Leiden University Medical Center, 2300 Leiden, The Netherlands; p.c.n.rensen@lumc.nl
Department of Physical and Sports Education, School of Sports Science, University of Granada, 18071 Granada, Spain
CERNEP Research Center, Department of Education, Faculty of Education Sciences and SPORT Research Group (CTS-1024), University of Almería, 04120 Almeria, Spain
Abstract
Objective: To investigate the association of plasma levels of endocannabinoids with fecal microbiota. Methods: Plasma levels of endocannabinoids, anandamide (AEA) and 2-arachidonoylglycerol (2-AG), as well as their eleven analogues, and arachidonic acid (AA), were measured using liquid chromatography-tandem mass spectrometry in 92 young adults. DNA extracted from stool samples was analyzed using 16S rRNA gene sequencing. Lipopolysaccharide levels were measured in plasma samples. Results: Plasma levels of endocannabinoids and their analogues were not related to beta or alpha diversity indexes. Plasma levels of AEA and related N-acylethanolamines correlated positively with the relative abundance of Faecalibacterium genus (all rho ≥ 0.26, p ≤ 0.012) and Akkermansia genus (all rho ≥ 0.22, p ≤ 0.036), and negatively with the relative abundance of Bilophila genus (all rho ≤ −0.23, p ≤ 0.031). Moreover, plasma levels of 2-AG and other acylglycerols correlated positively with the relative abundance of Parasutterella (all rho ≥ 0.24, p ≤ 0.020) and Odoribacter genera (all rho ≥ 0.27, p ≤ 0.011), and negatively with the relative abundance of Prevotella genus (all rho ≤ −0.24, p ≤ 0.023). In participants with high lipopolysaccharide values, the plasma levels of AEA and related N-acylethanolamines, as well as AA and 2-AG, were negatively correlated with plasma levels of lipopolysaccharide (all rho ≤ −0.24, p ≤ 0.020). Conclusion: Plasma levels of endocannabinoids and their analogues are correlated to specific fecal bacterial genera involved in maintaining gut barrier integrity in young adults. This suggests that plasma levels of endocannabinoids and their analogues may play a role in the gut barrier integrity in young adults.
Untitled section
Keywords: endocannabinoid system, gut barrier, gastrointestinal microbiome, inflammation, short-chain fatty acids
Article notes
Untitled section
Received 2022 Apr 9; Accepted 2022 May 18; Collection date 2022 May.
1. Introduction
The endocannabinoid system is mainly composed of cannabinoid receptors 1 and 2 (CB1 and CB2), several endogenous lipids called endocannabinoids, and metabolic enzymes involved in their synthesis and degradation [1]. Anandamide (AEA), a member of the N-acylethanolamines class, and 2-arachidonoylglycerol (2-AG), as part of the acylglycerol class, are the two primary endocannabinoids [2]. Structural analogues of these two endocannabinoids, as well as arachidonic acid (AA), which is the downstream metabolite of the primary endocannabinoids AEA and 2-AG, also belong to the endocannabinoid system. Endocannabinoids and their analogues are produced “on demand” in the brain, liver, adipose tissue, skeletal muscle, and pancreas [3]. Among other actions, they can exert either pro- or anti-inflammatory actions and modulate the immune response [4]. Interestingly, the endocannabinoid system seems to be involved in gut physiology [5].
The gastrointestinal tract is colonized by microbial communities called gut microbiota [6], which play a relevant role in regulating the innate and adaptive immune system, gut motility, gut barrier homeostasis, nutrient absorption and fat distribution [7]. Dysbiosis, defined as an imbalance in microbial communities [8], can be induced by multiple factors, such as an unhealthy diet [9] or the use of antibiotics and other drugs [10,11], and is related to obesity and non-communicable chronic diseases [12]. Dysbiosis induces gut barrier dysfunction, which leads to the translocation from the gut into the bloodstream of Gram-negative bacterial components, such as lipopolysaccharide [13]. This fact leads to metabolic endotoxemia and systemic inflammation, characterized by an increase in gut permeability [14]. Therefore, changes in the gut microbiota composition result in gut barrier dysfunction [15].
Recent evidence has shown that the endocannabinoid system can modulate the gut microbiota composition in mice models [5,16]. Conditional adipocyte-specific deficiency in mice of N-acyl phosphatidylethanolamine phospholipase D, an enzyme involved in the synthesis of N-acylethanolamines, increased body weight and fat mass, inducing gut microbiota dysbiosis [17]. Plasma levels of certain N-acylethanolamines were positively correlated with the relative abundance of fecal bacteria in middle adults [18]. Several studies have shown that the endocannabinoid system plays a role in regulating gastric secretion [5], motility [5], and gut permeability [19] by preserving or damaging the gut barrier integrity [20], by yet unknown mechanisms.
To date, there is no scientific evidence in humans on whether plasma levels of endocannabinoids and their analogues are associated with fecal microbiota genera involved in the maintenance of gut barrier integrity. This study aimed to investigate the association of plasma levels of endocannabinoids and their analogues with fecal microbiota diversity and composition in young adults.
2. Material and Methods
2.1. Study Design and Participants
This is a cross-sectional study that was conducted within the framework of the ACTIBATE study [21], an exercise-based randomized controlled trial (ClinicalTrials.gov ID: NCT02365129). The present study included baseline data of 92 healthy young adults (27 men and 65 women, age: 18–25 years old). All assessments were performed in Granada (Spain) between October and November 2016. The included participants (i) were physically inactive and had a sedentary lifestyle (less of 20 min moderate to vigorous physical activity on less of 3 days per week), (ii) had stable body weight during the last 3 months (<3 kg change), (iii) were not smokers, (iv) did not take any medication (including antibiotics) in the last 3 months, (v) did not present any acute neither chronic illness, and (vi) were not pregnant.
The participants signed an informed consent that together with study protocol were performed in accordance with the Declaration of Helsinki, as revised in 2013. Moreover, both were approved by the Human Research Ethics Committee of the University of Granada (n°924) and the Centro de Granada, CEI-Granada.
2.2. Anthropometry and Body Composition Assessment
We measured weight and height with a SECA scale and a stadiometer, respectively (model 799, Electronic Column Scale, Hamburg, Germany). Dual Energy X-ray Absorptiometry (DEXA, HOLOGIC, Discovery Wi, Marlborough, MA, USA) was used for determining lean mass, fat mass and visceral adipose tissue mass. Body mass index (BMI), lean mass index and fat mass index were computed as weight, lean body mass and fat body mass (kg), respectively, divided by height squared (m2). The fat mass percentage was determined as the body fat mass divided by the total body mass and multiplied by 100. We measured the waist circumference (cm) in the minimum perimeter, at the end of a normal expiration, with the arms of the participants were relaxed on both sides of the body. However, if the minimum perimeter could not be detected (i.e., participants with overweight/obesity), it was measured just above the umbilicus, in a horizontal plane. Waist circumference was determined as the average of two measures with a plastic tape measure.
2.3. Determination of Plasma Levels of Endocannabinoids and Their Analogues
Endocannabinoids (i.e., AEA, 2-AG) and their analogues including docosahexaenoylethanolamide (DHEA), dihomo-gamma-linolenoylethanolamide (DGLEA), linoleylethanolamine (LEA), alpha-linolenoylethanolamide (alpha LEA), palmitoylethanolamine (PEA), pentadecanoylethanolamide (PDEA), palmitoleoylethanolamide (POEA), oleoylethanolamine (OEA), and stearoylethanolamine (SEA), which belong to the N-acylethanolamines class, and 2-linoleoylglycerol (2-LG) and oleoylglycerol (2-OG), which belong to the acylglycerol class [22,23], as well as AA were analyzed in plasma samples using liquid chromatography-tandem mass spectrometry (LC-MS/MS). The method followed the FDA bioanalytical method validation guidelines [24].
2.3.1. Sample Preparation
The sample preparation procedure was carried out on ice to prevent analyte degradation. Endocannabinoids and their analogues were extracted using liquid–liquid extraction. Prior to the extraction, 5 µL of an antioxidant solution composed of 0.4 mg/mL of butylated hydroxytoluene (BHT) and 10 µL of an internal standard solution containing the isotopically labelled analogues were added to 150 µL of thawed plasma samples. Then, 150 µL of a buffer solution composed of 0.2 M citric acid and 0.1 M disodium hydrogen phosphate at pH 4.5 were added to the samples. Subsequently, 1 mL extraction solution composed of butanol (ButOH) and methyl-tertbutylether (MTBE) in a ratio 50:50 (v/v) was added prior to agitation for 5 min using a bullet blender and centrifugation at 16,000× g for 10 min at 4 °C. The organic supernatant (900 µL) was collected and evaporated to dryness using a SpeedVac instrument at room temperature. The dry residues were then reconstituted in 50 µL of an ice-cold solution of methanol (MetOH) and acetonitrile (ACN) 70:30 (v/v), prior to agitation (5 min) and centrifugation at 16,000× g for 10 min at 4 °C. Finally, 40 µL of the supernatant was transferred into a glass vial for further liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis.
2.3.2. Liquid Chromatography-Tandem Mass Spectrometry Analysis
Relative quantitation of endocannabinoids and their analogues were performed using a SCIEX QTRAP® LC-ESI-MS/MS System (SCIEX, Framingham, MA, USA). The separation was performed using a BEH C18 column (50 mm × 2.1 mm, 1.7 μm) from Waters Technologies (Milford, MA, USA), maintained at 40 °C. The mobile phase was composed of 0.1% acetic acid in water (A), ACN/0.1% acetic acid in MeOH (90:10, v/v) (B), and 0.1% acetic acid in isopropanol (C). The gradient was the following: starting condition 20% B and 1% C; increase in B from 20% to 26% between 0.75 min and 0.95 min; increase in B from 26% to 34% between 0.95 min and 6 min; increase in B from 34% to 40% between 6 min and 8 min; increase in B from 40% to 54% between 8 min and 10 min; increase in B from 54% to 56% between 10 min and 12 min; increase in C from 1% to 3% between 11 min and 12 min; increase in B from 56% to 78% between 12 min and 13 min; increase in C from 3% to 6% between 12 min and 13 min; increase in B from 78% to 85% between 13 min and 14 min; increase in C from 6% to 15% between 13 min and 14 min; conditions kept for 0.5 min prior to returning to initial conditions at 14.8 min and re-equilibration for 1.2 min. The flow rate was 0.7 mL/min. The injection volume was 10 μL, preceded by an injection of 20 μL of mobile phase A as stacking solution to improve the peak shape and increase sensitivity. MS acquisition was carried out in positive mode, with the following electrospray ionization (ESI) parameters: source temperature, 600 °C; Gas 1 (nebulizer gas), 50 L/min; Gas 2 (heater gas), 50 L/min; curtain gas, 30 L/min; collision gas, medium; ion spray voltage, 5500 V. Selected reaction monitoring was used for tandem MS experiments.
2.3.3. Data Pre-Processing
SCIEX OS-MQ Software was used for peak detection and integration. The ratio of the analyte peak area to the peak area of the corresponding isotopically labelled internal standard, referred to as peak area ratio, was used for further data analysis. Quality control samples were prepared using plasma samples from healthy subjects. Quality control samples were prepared simultaneously with the study samples and regularly injected throughout the sequences. The results obtained for quality control samples were used to evaluate the quality of the data, including blank effect, retention time shifts, and peak area ratios, as well as peak area ratios, and correct for between batch variations using the in-house developed mzQuality workflow (available at http://www.mzQuality.nl, (accessed on 1 February 2020)) [25]. We calculated relative standard deviations for each internal standard and analyte present in the quality control samples. Metabolites showing relative standard deviations higher than 30% on peak area ratios in quality control samples were excluded, whereas metabolites with relative standard deviations between 15% and 30% were treated with caution (Table 1). The acylglycerols are biologically present under two isomeric forms, namely, 1-AG and 2-AG, 1-LG and 2-LG, as well as 1-OG and 2-OG, respectively. For each isomer pair, a baseline separation was obtained between the two isomers with the developed LC-MS/MS method. However, the isomeric conversion due to acyl transformation of the 2-isomer into 1-isomer after sampling is a known mechanism, notably reported for 2-AG [26]. Since this conversion could not be experimentally controlled, we summed the peak areas of both isomers 1-AG and 2-AG before calculating the peak area ratio, and labeled the isomer pair “2-AG”. The same strategy was applied for the isomer pair 1-LG and 2-LG, as well as 1-OG and 2-OG, labeled “2-LG” and “2-OG”, respectively.
| N | Mean | ± | SD | |
|---|---|---|---|---|
| Sex (women, %) | 92 | 71% | ||
| Age (years) | 92 | 22 | ± | 2 |
| Anthropometry and body composition assessment | ||||
| Body mass index (kg/m2) | 90 | 24.87 | ± | 4.72 |
| Lean mass (kg) | 83 | 41.13 | ± | 8.98 |
| Fat mass (kg) | 83 | 24.91 | ± | 9.09 |
| Lean mass index (kg/m2) | 83 | 14.43 | ± | 2.24 |
| Fat mass index (kg/m2) | 83 | 8.84 | ± | 3.14 |
| Fat mass percentage (%) | 83 | 35.94 | ± | 7.88 |
| Visceral adipose tissue (g) | 83 | 321.05 | ± | 177.19 |
| Waist circumference (cm) | 90 | 80.60 | ± | 13.97 |
| Plasma levels of endocannabinoids and their analogues (peak area ratio) | ||||
| AEA and related N-acylethanolamines | ||||
| AEA | 88 | 0.14 | ± | 0.06 |
| AA | 88 | 65.18 | ± | 21.75 |
| DHEA * | 88 | 0.08 | ± | 0.16 |
| DGLEA * | 87 | 0.02 | ± | 0.01 |
| LEA | 88 | 0.22 | ± | 0.09 |
| alpha LEA | 88 | 0.008 | ± | 0.004 |
| PEA | 87 | 1.78 | ± | 0.28 |
| PDEA * | 88 | 0.03 | ± | 0.01 |
| POEA | 88 | 0.28 | ± | 0.22 |
| OEA | 88 | 0.70 | ± | 0.21 |
| SEA | 88 | 1.33 | ± | 0.23 |
| 2-AG and related acylglycerols | ||||
| 2-AG * | 88 | 0.02 | ± | 0.09 |
| 2-LG * | 88 | 0.22 | ± | 1.41 |
| 2-OG * | 88 | 0.01 | ± | 0.06 |
| Fecal microbiota parameters | ||||
| Alpha diversity indexes | ||||
| Richness Chao | 92 | 430.68 | ± | 150.99 |
| Shannon diversity | 92 | 4.21 | ± | 0.37 |
| Inverse Simpson diversity | 92 | 34.20 | ± | 14.16 |
| Evenness Camargo | 92 | 0.24 | ± | 0.05 |
| Composition (phylum) | ||||
| Actinobacteria (%) | 92 | 1.59 | ± | 1.55 |
| Bacteroidetes (%) | 92 | 39.80 | ± | 9.06 |
| Firmicutes (%) | 92 | 48.65 | ± | 10.10 |
| Proteobacteria (%) | 92 | 6.81 | ± | 5.29 |
| Verrucomicrobia (%) | 92 | 2.09 | ± | 4.12 |
| Plasma levels of lipopolysaccharide (EU/mL) | 85 | 0.98 | ± | 1.10 |
2.4. Fecal Microbiota Analyses
2.4.1. Stool Collection and DNA Extraction
The participants collected a fecal sample (~50–60 g) in plastic sterile containers, which they transported in a portable cooler to the research center. Fecal samples were stored at −80 °C until the extraction of DNA. We used a QIAamp DNA Stool Mini Kit (QIAGEN, Barcelona, Spain) to extract DNA, following the manufacturer’s instructions, and incubating the samples at 95 °C to ensure lysis of both Gram-positive and Gram-negative bacteria. The quantification of DNA was performed using a NanoDrop ND1000 spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA). The ratio of absorbance at A260/280 nm and A260/230 nm was used for measuring DNA purity.
2.4.2. Sequencing Analysis
DNA extracted was amplified by PCR by primer pairs, 16S Amplicon PCR Forward Primer: 5′CCTACGGGNGGCWGCAG, and 16S Amplicon PCR Reverse Primer: 5′GACTACHVGGGTATCTAATCC targeting the V3 and V4 hypervariable regions of the bacterial 16S rRNA gene [27]. We performed all PCRs using 25 µL reaction volumes, incorporating 12.5 µL 2X KAPA HiFi Hotstart ready mix (KAPA Biosystems, Woburn, MA, USA), 5 µL of each primer (1 µM) and 2.5 µL of extracted DNA (10 ng). The cycle steps were as follows: (i) denaturation at 95 °C for 3 min, (ii) 8 cycles of denaturation at 95 °C for 30 s, (iii) annealing at 55 °C for 30 s, (iv) elongation at 72 °C for 30 s, (v) final extension at 72 °C for 5 min. AMPure XP beads (Beckman Coulter, Indianapolis, IN, USA) were used to purify the 16S V3 and V4 amplicons from free primers and primer dimers. The index PCR was performed using the Nextera XT Index Kit (Illumina, San Diego, CA, USA), which attaches dual indices and Illumina sequencing adapters, on a thermal cycler using the requirements previously mentioned. Before quantification, AMPure XP beads (Beckman Coulter, Indianapolis, IN, USA) were used for purifying the pooled PCR products. The resulting amplicons were sequenced at MiSeq (Illumina, San Diego, CA, USA), using a paired-end (2 × 300 nt) Illumina MiSeq sequencing system (Illumina, San Diego, CA, USA).
2.4.3. Bioinformatics Analysis
FastQ files were analyzed with the “dada2” [28] package in R software [29]. We obtained 11,659,014 paired-ends with an average of 126,728 ± 33,395 reads per sample. All samples surpassed a cut-off of 10,000 reads. Samples were resampled to an equal sequencing depth of 30,982 reads using the “phyloseq” [30] package in R [29], retrieving 11,158 phylotypes.
The taxonomic affiliation of phylotypes was assigned using the “classifier” function from the Ribosomal Database Project (RDP), based on the naive Bayesian classification [31] with a pseudo-bootstrap threshold of 80%. A total of 209 genera that belong to 16 different phyla were obtained. In order to determinate further annotation of the phylotypes to species assignments, the “seqmatch” [32] function from RDP was performed to define the discriminatory power of each sequence read; we executed annotation according to criteria previously published [33]. Microbial communities were analyzed at different taxonomic levels (phylum to species), calculating the relative abundances, which were expressed as percentages. We performed the analyses using the genera with an average relative abundance of more than 0.5%. Only phylotypes found in at least 50% of the participants were submitted to species taxonomy level assignation.
Beta and alpha diversities were estimated based on the identified microbial communities. Beta diversity shows the differences in microbial community composition between individuals, i.e., the degree to which samples differ from one another [34], whereas the alpha diversity considers the number of different phylotypes and relative abundances within a given individual [35]. The beta diversity was quantitatively measured using permutational multivariate analysis of variance (PERMANOVA) based on the Bray–Curtis [36] similarity. The alpha diversity was assessed using richness Chao, Shannon, inverse Simpson and evenness Camargo indexes with the “microbiome” [37] package using R software [29]. The Richness Chao index estimates the diversity according to the number of different phylotypes in the community [38]; the Shannon diversity increases as both the richness and the evenness of the community increase [39]; the inverse of Simpson diversity is derived from the classical Simpson diversity and indicates the richness in a community with uniform evenness [40]; the evenness Camargo index determines the equitability of phylotypes frequencies in a community [41].
2.5. Determination of Plasma Levels of Lipopolysaccharides
Plasma levels of lipopolysaccharide were determined by a chromogenic limulus amebocyte lysate (LAL) assay according to the manufacturer guidelines (LAL Chromogenic Endpoint Assay; HycultBiotech, Uden, The Netherlands) as previously reported [42]. With the purpose of inactivating endotoxin neutralizing agents, we diluted the plasma samples in pyrogen-free water and heated them to 70 °C for 10 min. All samples were measured in duplicate, and results were accepted when the intra-assay coefficient of variation was <10%. Internal recovery controls were included in the assessment.
2.6. Statistical Analysis
Data are presented as means ± standard deviations unless otherwise stated. All variables were tested for normality using the D’Agostino and Pearson omnibus test. Since most of the variables did not display a normal distribution, non-parametric tests were used for all analyses. It has been unequivocally demonstrated that human plasma levels of endocannabinoids and their analogues are related to adiposity markers such as BMI and visceral adipose tissue mass [43,44]. Based on that, we performed and presented all the analyses adjusting for BMI and the data from both sexes pooled together, since we did not detect any sex interaction (all p > 0.05). To calculate beta diversity, plasma levels of endocannabinoids and their analogues were divided into tertiles (low, intermediate and high groups) and compared using a two-way PERMANOVA with 9999 permutations for significance testing with the Paleontological Statistics (Past3) software. Partial Spearman correlations were used to investigate the correlation of plasma levels of endocannabinoids and their analogues with fecal microbiota alpha diversity and composition, using the “psych” [45] and “corrplot” [46] packages in R software. Next, we divided the participants into quartiles of plasma levels of lipopolysaccharide, and we performed partial Spearman correlations, as described above, between the plasma levels of endocannabinoids and their analogues with plasma levels of lipopolysaccharide in Q1, Q4, and the whole cohort. Kruskal–Wallis tests were performed to investigate whether there were significant differences in body composition, plasma levels of endocannabinoids and their analogues and gut microbiota composition outcomes between the different quartiles of plasma levels of lipopolysaccharide with GraphPad Prism. The level of significance was assumed at p < 0.05. R software (V.3.6.0; http://www.R-project.org, (accessed on 1 June 2019)) and GraphPad Prism version 8.0.0 for Windows (GraphPad Software, San Diego, CA, USA) were used for the statistical analysis and graphical plots.
3. Results
3.1. Characteristics of the Participants
Table 1 shows the descriptive characteristics of the participants. AEA and seven of its analogues (i.e., AA, LEA, alpha LEA, PEA, POEA, OEA and SEA) showed relative standard deviations values of ≤15%, whereas 2-AG and five of its analogues (i.e., DHEA, DGLEA, PDEA, 2-LG and 2-OG) showed relative standard deviations values between 15% and 30%.
3.2. The Plasma Levels of Endocannabinoids and Their Analogues Are Not Related to Fecal Microbiota Beta and Alpha Diversities
The beta diversity was similar across the tertiles of the plasma endocannabinoids and their analogues at both the phylum and genus taxonomic levels (all p ≥ 0.312; Table 2). Overall, the plasma levels of endocannabinoids and their analogues were not related to the alpha diversity indexes (Figure 1). However, we observed that plasma levels of DHEA (related N-acylethanolamines) were positively correlated with the richness Chao index (rho = 0.23, p = 0.032; Figure 1), whereas plasma levels of 2-LG and 2-OG (related acylglycerols) were negatively correlated with the evenness Camargo index (rho = −0.30, p = 0.004, and rho = −0.31, p = 0.003, respectively; Figure 1). We repeated all the analyses adjusting for sex, BMI and energy intake in separate models, and the results remained the same (data not shown).
| Tertiles of the Plasma Levels of Endocannabinoids and Their Analogues (Peak Area Ratio) | Phylum | Genus | ||
|---|---|---|---|---|
| Pseudo-F | p-Value | Pseudo-F | p-Value | |
| AEA and related N-acylethanolamines | ||||
| AEA | −6.626 | 0.383 | −6.789 | 0.336 |
| AA | −5.918 | 0.348 | −6.275 | 0.397 |
| DHEA | −6.581 | 0.548 | −6.620 | 0.684 |
| DGLEA | −6.339 | 0.657 | −6.146 | 0.579 |
| LEA | −6.469 | 0.405 | −6.670 | 0.418 |
| alpha LEA | −6.211 | 0.361 | −6.854 | 0.746 |
| PEA | −6.032 | 0.591 | −5.996 | 0.723 |
| PDEA | −5.818 | 0.312 | −6.370 | 0.485 |
| POEA | −6.495 | 0.423 | −6.516 | 0.336 |
| OEA | −6.410 | 0.372 | −6.645 | 0.369 |
| SEA | −6.460 | 0.373 | −6.797 | 0.442 |
| 2-AG and related acylglycerols | ||||
| 2-AG | −6.281 | 0.640 | −5.808 | 0.473 |
| 2-LG | −6.449 | 0.570 | −6.326 | 0.595 |
| 2-OG | −7.041 | 0.578 | −6.927 | 0.629 |
3.3. The Plasma Levels of Endocannabinoids and Their Analogues Levels Are Related to the Relative Abundance of Specific Fecal Microbiota Genera
The plasma levels of AEA and related N-acylethanolamines, and AA were positively correlated with the relative abundance of Verrucomicrobia and Firmicutes phyla, and negatively with the relative abundance of Proteobacteria phylum. Similarly, the plasma levels of 2-AG and related acylglycerols were positively correlated with the relative abundance of Proteobacteria phylum. Specifically, we observed that the plasma levels of LEA, alpha LEA, PDEA, POEA and OEA levels were positively correlated with the relative abundance of Akkermansia genus (Verrucomicrobia phylum; all rho ≥ 0.22, p ≤ 0.036; Figure 2A). Additionally, the plasma levels of PEA and SEA were positively correlated with the relative abundance of Lachnospiraceae incertae sedis (Firmicutes phylum; all rho ≥ 0.26, p ≤ 0.010; Figure 2A) and Faecalibacterium genera (Firmicutes phylum; all rho ≥ 0.26, p ≤ 0.012; Figure 2A). On the opposite, the plasma levels of AEA and related N-acylethanolamines (i.e., DGLEA, LEA, PDEA and OEA), as well as AA, were negatively correlated with the relative abundance of Bilophila genus (Proteobacteria phylum; all rho ≤ −0.23, p ≤ 0.031; Figure 2A). Similarly, the plasma levels of alpha LEA, PEA and SEA were negatively correlated with the relative abundance of Succinivibrio genus (Proteobacteria phylum; both rho ≤ −0.23, p ≤ 0.027; Figure 2A). Further, the plasma levels of 2-AG and the related acylglycerols 2-LG and 2-OG were positively correlated with the relative abundance of Parasutterella (Proteobacteria phylum; all rho ≥ 0.24, p ≤ 0.020; Figure 2B) and Odoribacter genera (Bacteroidetes phylum; all rho ≥ 0.27, p ≤ 0.011; Figure 2B). In contrast, the plasma levels of 2-AG and 2-LG only were negatively correlated with the relative abundance of Prevotella genus (Bacteroidetes phylum; rho ≤ −0.24, p ≤ 0.023; Figure 2B). The results were similar when visceral adipose tissue mass was included as a confounder instead of BMI or physical activity levels (data not shown).
3.4. The Plasma Levels of Endocannabinoids and Their Analogues Are Related to the Relative Abundance of Specific Bacterial Species
We investigated whether the relative abundance of specific bacterial species could drive the significant correlations observed between the plasma levels of endocannabinoids and their analogues with the relative abundance of different genera (Figure 3). We found that the plasma levels of AEA, AA, OEA and SEA were positively correlated with the relative abundance of Bacteroides dorei (Bacteroidetes phylum; all rho ≥ 0.206, p < 0.05; Figure 3A), whereas the plasma levels of alpha LEA and POEA were positively associated with the relative abundance of Parasutterella excrementihominis (Proteobacteria phylum; all rho ≥ 0.236, p ≤ 0.024; Figure 3A). Additionally, the plasma levels of SEA were positively associated with the relative abundance of Faecalibacterium prausnitzii (Firmicutes phylum; rho = 0.245, p = 0.019; Figure 3A). Opposite to this, the plasma levels of AEA and DHEA, as well as AA, were negatively correlated with the relative abundance of Bacteroides vulgatus (Bacteroidetes phylum; all rho ≤ −0.211, p ≤ 0.045; Figure 3A). Moreover, the plasma levels of AA and OEA, as well as 2-AG and the related acylglycerol 2-LG were negatively correlated with the relative abundance of Roseburia inulinivorans (Firmicutes phylum; all rho ≤ −0.219, p ≤ 0.037; Figure 3A,B). All associations remained significant when visceral adipose tissue mass was included as a confounder instead of BMI or physical activity levels (data not shown).
3.5. The Plasma Levels of Endocannabinoids and Their Analogues Are Negatively Correlated to Plasma Levels of Lipopolysaccharides but Only in Those Participants with High Plasma Levels of Lipopolysaccharides
We quantified plasma levels of lipopolysaccharide since it is an indirect marker of gut permeability [47]. We observed that plasma levels of endocannabinoids and their analogues were not correlated with the plasma levels of lipopolysaccharide (p ≥ 0.385; Tables S1 and S2). We then divided the cohort into quartiles of plasma levels of lipopolysaccharide: Q1 (0.14–0.35 EU/mL), Q2 (0.36–0.50 EU/mL), Q3 (0.50–1.22 EU/mL) and Q4 (1.35–5.45 EU/mL). We followed this approach because we wanted to focus on the comparison of extreme values, as well as to explore non-lineal associations. Overall, body composition, plasma levels of endocannabinoids and their analogues, and fecal microbiota composition were similar across the quartiles (Table S3). We found that plasma levels of AEA and related N-acylethanolamines (except for DHEA), as well as AA and 2-AG were positively correlated with plasma levels of lipopolysaccharide (all rho ≥ 0.22, p ≤ 0.035; Table 3; Table 4) only in those participants in Q1. On the opposite, plasma levels of AEA and related N-acylethanolamines (i.e., DGLEA, LEA, PEA, PDEA, POEA and OEA), as well as AA and 2-AG were negatively correlated with plasma levels of lipopolysaccharide (all rho ≤ −0.24, p ≤ 0.020; Table 3; Table 4) in participants in Q4. Overall, we found no correlation between the relative abundance of different genera in the feces with the plasma levels of lipopolysaccharide (Tables S4 and S5).
| AEA and Related N-acylethanolamines (Peak Area Ratio) | ||||||||||||||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
| AEA | AA | DHEA | DGLEA | LEA | Alpha LEA | PEA | PDEA | POEA | OEA | SEA | ||||||||||||
| rho | p | rho | p | rho | p | rho | p | rho | p | rho | p | rho | p | rho | p | rho | p | rho | p | rho | p | |
| Q1: Low LPS (N = 19) | 0.262 | 0.012 | 0.228 | 0.030 | 0.123 | 0.246 | 0.221 | 0.035 | 0.404 | <0.001 | 0.428 | <0.001 | 0.451 | <0.001 | 0.234 | 0.026 | 0.467 | <0.001 | 0.440 | <0.001 | 0.428 | <0.001 |
| Q4: High LPS (N = 21) | −0.464 | <0.001 | −0.389 | <0.001 | −0.001 | 0.992 | −0.472 | <0.001 | −0.300 | 0.004 | 0.107 | 0.315 | −0.304 | 0.003 | −0.281 | 0.007 | −0.504 | <0.001 | −0.408 | <0.001 | −0.156 | 0.139 |
| 2-AG and Related Acylglycerols (Peak Area Ratio) | ||||||
|---|---|---|---|---|---|---|
| 2-AG | 2-LG | 2-OG | ||||
| rho | p | rho | p | rho | p | |
| Q1: Low LPS (N = 19) | 0.374 | <0.001 | 0.111 | 0.295 | 0.171 | 0.106 |
| Q4: High LPS (N = 21) | −0.243 | 0.020 | 0.312 | 0.003 | 0.053 | 0.615 |
4. Discussion
The present study shows that the plasma levels of N-acylethanolamines were are positively correlated with the relative abundance of Akkermansia and Faecalibacterium genera and negatively correlated with the relative abundance of Bilophila genus. The plasma levels of acylglycerols were positively correlated with the relative abundance of Parasutterella and Odoribacter genera and negatively correlated with the relative abundance of Prevotella genus. Moreover, we observed a negative and significant correlation between the plasma levels of endocannabinoids and their analogues with the plasma levels of lipopolysaccharide, but only in participants within Q4. These findings suggest that the plasma levels of endocannabinoids and their analogues may be related in the gut barrier integrity through the gut microbiota in young adults [16]. Future studies are needed to ascertain whether endocannabinoids and their analogues could be considered markers of gut barrier integrity in humans.
4.1. Role of Endocannabinoids and Their Analogues in Gut Microbiota Diversity
In agreement with previous studies [48], our findings show that the beta diversity was similar across the different plasma endocannabinoids and their analogues tertiles. This could be partially explained by the homogeneity in age and cardiometabolic status of our participants (Table 1). Curiously, the plasma levels of DHEA were positive, although weakly, correlated with the richness Chao, an alpha diversity estimator, according to the number of species present in the fecal sample [38]. However, the plasma levels of 2-LG and 2-OG were negatively related to the evenness Camargo, which is a proxy of the equality in the numbers of species within the bacterial community [41]. This might indicate that an increase in the plasma levels of 2-LG and 2-OG may be related to the evenness in the bacterial community of young adults, which warrants further investigation.
4.3. Limitations and Strengths
The present study suffers the limitation of its cross-sectional design, which precludes the establishment of cause-effect relationships. Moreover, sample size was not calculated a priori, and thus, this study follows a hypothesis generation approach. The study population only included young adults, which does not enable extrapolation of results to older, younger, or unhealthy populations. Moreover, the associations observed were relatively weak. The strengths of this study are that we sequenced using the latest technology (Illumina platform) and annotated with RDP at the level of species taxon. Furthermore, this study is one of the few studies to date that have studied the relationship between endocannabinoids and their analogues with fecal microbiota diversity and composition in young adults.
5. Conclusions
The plasma levels of LEA, alpha LEA, PDEA, POEA, OEA, PEA, SEA, 2-AG, 2-LG, and 2-OG were positively correlated with the relative abundance of Akkermansia, Faecalibacterium, Parasutterella and/or Odoribacter genera, which may contribute to the preservation of the gut barrier integrity. Moreover, the plasma levels AEA and related N-acylethanolamines (i.e., DGLEA, LEA, PEA, PDEA, POEA and OEA), as well as AA and 2-AG, were negatively correlated with the plasma levels of lipopolysaccharide, but only in participants with high levels of lipopolysaccharide. These findings suggest that the plasma levels of endocannabinoids and their analogues are related to specific fecal bacteria genera and may play a role in the gut barrier integrity in young adults.
Acknowledgments
This study is part of a PhD thesis conducted within the Biomedicine Doctoral Studies Program of the University of Granada, Spain.
Abbreviations
| 2-AG | 2-arachidonoylglycerol; |
| 2-LG | 2-linoleoylglycerol; |
| 2-OG | 2-oleoylglycerol; |
| AA | arachidonic acid; |
| AEA | anandamide; |
| alpha LEA | alpha-Linolenoyl ethanolamide; |
| DGLEA | dihomo-gamma-linolenoyl ethanolamide; |
| DHEA | docosahexaenoyl ethanolamide; |
| LEA | linoleoyl ethanolamide; |
| NAEs | N-acylethanolamines; |
| OEA | oleoyl ethanolamine; |
| PDEA | pentadecanoyl ethanolamide; |
| PEA | palmitoyl ethanolamide; |
| POEA | palmitoleoyl ethanolamide; |
| SEA | stearoyl ethanolamide; |
| UC | unclassified. |
Supplementary Materials
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nu14102143/s1, Table S1: Partial Spearman correlation of plasma levels of AEA and related N-acylethanolamines with plasma levels of lipopolysaccharide adjusted for body mass index. Table S2: Partial Spearman correlation of plasma levels of 2-AG and related acylglycerols with plasma levels of lipopolysaccharide adjusted for body mass index. Table S3: Characteristics of the individuals according to quartiles of plasma levels of lipopolysaccharide. Table S4: Spearman correlation of relative abundance of fecal microbiota composition at genera taxonomic level with plasma levels of lipopolysaccharide. Table S5: Spearman correlation of relative abundance of fecal microbiota composition at genera taxonomic level with plasma levels of lipopolysaccharide in participants who showed very low and high values (EU/mL).
Institutional Review Board Statement
The study was conducted in accordance with the Declaration of Helsinki, as revised in 2013, and approved by the Human Research Ethics Committee of the University of Granada (n°924) and the Centro de Granada, CEI-Granada. ClinicalTrials.gov (accessed on 1 February 2015) no. NCT02365129 (registered 18 February 2015).
Informed Consent Statement
Informed consent was obtained from all subjects involved in the study.
Data Availability Statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.
Conflicts of Interest
The authors have no conflict of interest.
Funding Statement
The study was supported by the Spanish Ministry of Economy and Competitiveness via Fondo de Investigación Sanitaria del Instituto de Salud Carlos III (PI13/01393) and PTA 12264-I, Retos de la Sociedad (DEP2016-79512-R), and European Regional Development Funds (ERDF), by the Spanish Ministry of Education (FPU16/05159, FPU16/02828 and FPU17/01523), the Fundación Iberoamericana de Nutrición (FINUT), the Redes Temáticas De Investigación Cooperativa RETIC (Red SAMID RD16/0022), the Chinese Scholarship Council (CSC, No. 201707060012 to XD), InFLAMES Flagship Programme of the Academy of Finland (decision number: 337530), Fundación Alfonso Martin Escudero and Maria Zambrano fellowship by the Ministerio de Universidades y la Unión Europea—NextGenerationEU (RR_C_2021_04). AstraZeneca HealthCare Foundation, the University of Granada Plan Propio de Investigación 2016-Excellence actions: Unit of Excellence on Exercise and Health (UCEES), and by the Junta de Andalucía, Consejería de Economía, Conocimiento, Empresas y Universidad (ref. P18-RT-4455). AL and RVV are supported by the funds of European Commission through the “European funds for regional development” (EFRE) as well as by the regional Ministry of Economy, Science and Digitalization of Saxony-Anhalt as part of the “Autonomy in old Age“ (AiA) research group for “LiLife” Project (Project ID: ZS/2018/11/95324). M.C-P was recipient of a postdoctoral grant (DOC_00448) from the Consejeria de Economía, Industria, Conocimiento y Universidades (PAIDI 2020, Junta de Andalucía), Spain, cofunded by the Fondo Europeo de Desarrollo Regional (FEDER). We would like to thank the team of Data Integration Center of University Medicine Magdeburg for local data-analysis solutions; they are supported by MIRACUM and funded by the German Federal Ministry of Education and Research (BMBF) within the “Medical Informatics Funding Scheme” (FKZ 01ZZ1801H).
Footnotes
Footnote Group
References
Untitled section
References
- 1.Muccioli G.G. Endocannabinoid Biosynthesis and Inactivation, from Simple to Complex. Drug Discov. Today. 2010;15:474–483. doi: 10.1016/j.drudis.2010.03.007.
- 2.Devane W.A., Hanuš L., Breuer A., Pertwee R.G., Stevenson L.A., Griffin G., Gibson D., Mandelbaum A., Etinger A., Mechoulam R. Isolation and Structure of a Brain Constituent That Binds to the Cannabinoid Receptor. Science. 1992;258:1946–1949. doi: 10.1126/science.1470919.
- 3.Cristino L., Becker T., Di Marzo V. Endocannabinoids and Energy Homeostasis: An Update. BioFactors. 2014;40:389–397. doi: 10.1002/biof.1168.
- 4.Chiurchiù V., Leuti A., Maccarrone M. Bioactive Lipids and Chronic Inflammation: Managing the Fire Within. Front. Immunol. 2018;9:1. doi: 10.3389/fimmu.2018.00038.
- 5.Meletis C. The Important Role of the Endocannabinoid System and the Endocannabinoidome in Gut Health. Altern. Ther. Health Med. 2019;25:24–28.
- 6.Rosenbaum M., Knight R., Leibel R.L. The Gut Microbiota in Human Energy Homeostasis and Obesity. Trends Endocrinol. Metab. 2015;26:493–501. doi: 10.1016/j.tem.2015.07.002.
- 7.Everard A., Cani P.D. Diabetes, Obesity and Gut Microbiota. Best Pract. Res. Clin. Gastroenterol. 2013;27:73–83. doi: 10.1016/j.bpg.2013.03.007.
- 8.Iebba V., Totino V., Gagliardi A., Santangelo F., Cacciotti F., Trancassini M., Mancini C., Cicerone C., Corazziari E., Pantanella F., et al. Eubiosis and Dysbiosis: The Two Sides of the Microbiota. New Microbiol. 2016;39:1–12.
- 9.David L.A., Maurice C.F., Carmody R.N., Gootenberg D.B., Button J.E., Wolfe B.E., Ling A.V., Devlin A.S., Varma Y., Fischbach M.A., et al. Diet Rapidly and Reproducibly Alters the Human Gut Microbiome. Nature. 2014;505:559–563. doi: 10.1038/nature12820.
- 10.Maier L., Pruteanu M., Kuhn M., Zeller G., Telzerow A., Anderson E.E., Brochado A.R., Fernandez K.C., Dose H., Mori H., et al. Extensive Impact of Non-Antibiotic Drugs on Human Gut Bacteria. Nature. 2018;555:623–628. doi: 10.1038/nature25979.
- 11.Forslund S.K., Chakaroun R., Zimmermann-Kogadeeva M., Markó L., Aron-Wisnewsky J., Nielsen T., Moitinho-Silva L.B., Schmidt T.S., Falony G., Vieira-Silva S., et al. Combinatorial, Additive and Dose-Dependent Drug–Microbiome Associations. Nataliya Sokolovska. 2021;600:500. doi: 10.1038/s41586-021-04177-9.
- 12.Gérard P. Gut Microbiota and Obesity. Cell. Mol. Life Sci. 2016;73:147–162. doi: 10.1007/s00018-015-2061-5.
- 13.Cani P.D., Bibiloni R., Knauf C., Waget A., Neyrinck A.M., Delzenne N.M., Burcelin R. Changes in Gut Microbiota Control Metabolic Endotoxemia-Induced Inflammation in High-Fat Diet–Induced Obesity and Diabetes in Mice. Diabetes. 2008;57:1470. doi: 10.2337/db07-1403.
- 14.Cani P.D., Possemiers S., Van de Wiele T., Guiot Y., Everard A., Rottier O., Geurts L., Naslain D., Neyrinck A., Lambert D.M., et al. Changes in Gut Microbiota Control Inflammation in Obese Mice through a Mechanism Involving GLP-2-Driven Improvement of Gut Permeability. Gut. 2009;58:1091–1103. doi: 10.1136/gut.2008.165886.
- 15.Cani P.D., Everard A., Duparc T. Gut Microbiota, Enteroendocrine Functions and Metabolism. Curr. Opin. Pharmacol. 2013;13:935–940. doi: 10.1016/j.coph.2013.09.008.
- 16.Cani P.D., Plovier H., Van Hul M., Geurts L., Delzenne N.M., Druart C., Everard A. Endocannabinoids—At the Crossroads between the Gut Microbiota and Host Metabolism. Nat. Rev. Endocrinol. 2016;12:133–143. doi: 10.1038/nrendo.2015.211.
- 17.Geurts L., Everard A., Van Hul M., Essaghir A., Duparc T., Matamoros S., Plovier H., Castel J., Denis R.G.P., Bergiers M., et al. Adipose Tissue NAPE-PLD Controls Fat Mass Development by Altering the Browning Process and Gut Microbiota. Nat. Commun. 2015;6:6495. doi: 10.1038/ncomms7495.
- 18.Tagliamonte S., Laiola M., Ferracane R., Vitale M., Gallo M.A., Meslier V., Pons N., Ercolini D., Vitaglione P. Mediterranean Diet Consumption Affects the Endocannabinoid System in Overweight and Obese Subjects: Possible Links with Gut Microbiome, Insulin Resistance and Inflammation. Eur. J. Nutr. 2021;60:3703–3716. doi: 10.1007/s00394-021-02538-8.
- 19.Oláh A., Szekanecz Z., Bíró T. Targeting Cannabinoid Signaling in the Immune System: “High”-Ly Exciting Questions, Possibilities, and Challenges. Front. Immunol. 2017;8:1. doi: 10.3389/fimmu.2017.01487.
- 20.Muccioli G.G., Naslain D., Bäckhed F., Reigstad C.S., Lambert D.M., Delzenne N.M., Cani P.D. The Endocannabinoid System Links Gut Microbiota to Adipogenesis. Mol. Syst. Biol. 2010;6:392. doi: 10.1038/msb.2010.46.
- 21.Sanchez-Delgado G., Martinez-Tellez B., Olza J., Aguilera C.M., Labayen I., Ortega F.B., Chillon P., Fernandez-Reguera C., Alcantara J.M.A., Martinez-Avila W.D., et al. Activating Brown Adipose Tissue through Exercise (ACTIBATE) in Young Adults: Rationale, Design and Methodology. Contemp. Clin. Trials. 2015;45:416–425. doi: 10.1016/j.cct.2015.11.004.
- 22.Ben-Shabat S., Fride E., Sheskin T., Tamiri T., Rhee M.-H., Vogel Z., Bisogno T., De Petrocellis L., Di Marzo V., Mechoulam R. An Entourage Effect: Inactive Endogenous Fatty Acid Glycerol Esters Enhance 2-Arachidonoyl-Glycerol Cannabinoid Activity. Eur. J. Pharmacol. 1998;353:23–31. doi: 10.1016/S0014-2999(98)00392-6.
- 23.Lambert D.M., Di Marzo V. The Palmitoylethanolamide and Oleamide Enigmas: Are These Two Fatty Acid Amides Cannabimimetic? Curr. Med. Chem. 1999;6:757–773. doi: 10.2174/0929867306666220401153732.
- 24.FDA . Bioanalytical Method Validation Guidance. Vol. 1043. FDA; Silver Spring, MD, USA: 2018. p. 25.
- 25.Van der Kloet F.M., Bobeldijk I., Verheij E.R., Jellema R.H. Analytical Error Reduction Using Single Point Calibration for Accurate and Precise Metabolomic Phenotyping. J. Proteome Res. 2009;8:5132–5141. doi: 10.1021/pr900499r.
- 26.Martin J.B. The Equilibrium between Symmetrical and Unsymmetrical Monoglycerides and Determination of Total Monoglycerides. J. Am. Chem. Soc. 1953;75:5483–5486. doi: 10.1021/ja01118a005.
- 27.Pr Herlemann D., Labrenz M., Jü Rgens K., Bertilsson S., Waniek J.J., Andersson A.F. Transitions in Bacterial Communities along the 2000 Km Salinity Gradient of the Baltic Sea. ISME J. 2011;5:1571–1579. doi: 10.1038/ismej.2011.41.
- 28.Callahan B.J., Mcmurdie P.J., Rosen M.J., Han A.W., Johnson A.J.A., Holmes S.P. DADA2: High Resolution Sample Inference from Illumina Amplicon Data. Nat. Methods. 2016;13:581–583. doi: 10.1038/nmeth.3869.
- 29.R Core Team . A Language and Environment for Statistical Computing. R Core Team; Vienna, Austria: 2019. [(accessed on 1 March 2021)]. Available online: https://www.R-project.org/
- 30.McMurdie P.J., Holmes S. Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PLoS ONE. 2013;8:e61217. doi: 10.1371/journal.pone.0061217.
- 31.Wang Q., Garrity G.M., Tiedje J.M., Cole J.R. Naïve Bayesian Classifier for Rapid Assignment of RRNA Sequences into the New Bacterial Taxonomy. Appl. Environ. Microbiol. 2007;73:5261–5267. doi: 10.1128/AEM.00062-07.
- 32.Cole J.R., Wang Q., Fish J.A., Chai B., McGarrell D.M., Sun Y., Brown C.T., Porras-Alfaro A., Kuske C.R., Tiedje J.M. Ribosomal Database Project: Data and Tools for High Throughput RRNA Analysis. Nucleic Acids Res. 2014;42:D633–D642. doi: 10.1093/nar/gkt1244.
- 33.Schulz C., Schütte K., Koch N., Vilchez-Vargas R., Wos-Oxley M.L., Oxley A.P.A., Vital M., Malfertheiner P., Pieper D.H. The Active Bacterial Assemblages of the Upper GI Tract in Individuals with and without Helicobacter Infection. Gut. 2018;67:216–225. doi: 10.1136/gutjnl-2016-312904.
- 34.Lozupone C.A., Knight R. Species Divergence and the Measurement of Microbial Diversity. FEMS Microbiol. Rev. 2008;32:557–578. doi: 10.1111/j.1574-6976.2008.00111.x.
- 35.Finotello F., Mastrorilli E., Di Camillo B. Measuring the Diversity of the Human Microbiota with Targeted Next-Generation Sequencing. Brief. Bioinform. 2018;19:679–692. doi: 10.1093/bib/bbw119.
- 36.Bray J.R., Curtis J.T. An Ordination of the Upland Forest Communities of Southern Wisconsin. Ecol. Monogr. 1957;27:325–349. doi: 10.2307/1942268.
- 37.Leo L., Shetty S. Microbiome R Package. 2019. [(accessed on 1 March 2021)]. Available online: http://microbiome.github.io.
- 38.Gotelli N.J., Colwell R.K. Estimating Species Richness. Biol. Divers. Front. Meas. Assess. 2011;12:39–54.
- 39.Magurran A.E. Measuring Biological Diversity. Blackwell; Oxford, UK: 2004.
- 40.Simpson E.H. Measurement of Diversity. Nature. 1949;163:688. doi: 10.1038/163688a0.
- 41.Camargo J.A. New Diversity Index for Assessing Structural Alterations in Aquatic Communities. Bull. Environ. Contam. Toxicol. 1992;48:428–434. doi: 10.1007/BF00195643.
- 42.Clemente-Postigo M., Oliva-Olivera W., Coin-Aragüez L., Ramos-Molina B., Giraldez-Perez R.M., Lhamyani S., Alcaide-Torres J., Perez-Martinez P., El Bekay R., Cardona F., et al. Metabolic Endotoxemia Promotes Adipose Dysfunction and Inflammation in Human Obesity. Am. J. Physiol. Metab. 2019;316:E319–E332. doi: 10.1152/ajpendo.00277.2018.
- 43.Castonguay-Paradis S., Lacroix S., Rochefort G., Parent L., Perron J., Martin C., Lamarche B., Raymond F., Flamand N., Di Marzo V., et al. Dietary Fatty Acid Intake and Gut Microbiota Determine Circulating Endocannabinoidome Signaling beyond the Effect of Body Fat. Sci. Rep. 2020;10:15975. doi: 10.1038/s41598-020-72861-3.
- 44.Fanelli F., Mezzullo M., Repaci A., Belluomo I., Ibarra Gasparini D., Di Dalmazi G., Mastroroberto M., Vicennati V., Gambineri A., Morselli-Labate A.M., et al. Profiling Plasma N-Acylethanolamine Levels and Their Ratios as a Biomarker of Obesity and Dysmetabolism. Mol. Metab. 2018;14:82–94. doi: 10.1016/j.molmet.2018.06.002.
- 45.Revelle W.R. Psych: Procedures for Personality and Psychological Research. R Package Version 1.9.12. 2017. [(accessed on 1 March 2021)]. Available online: http://personality-project.org/r.
- 46.Wei T., Simko V., Levy M., Xie Y., Jin Y., Zemla J. R Package “Corrplot”: Visualization of a Correlation Matrix. 2017. [(accessed on 1 March 2021)]. Available online: https://github.com/taiyun/corrplot.
- 47.Massier L., Blüher M., Kovacs P., Chakaroun R.M. Impaired Intestinal Barrier and Tissue Bacteria: Pathomechanisms for Metabolic Diseases. Front. Endocrinol. 2021;12:1. doi: 10.3389/fendo.2021.616506.
- 48.Yatsunenko T., Rey F.E., Manary M.J., Trehan I., Dominguez-Bello M.G., Contreras M., Magris M., Hidalgo G., Baldassano R.N., Anokhin A.P., et al. Human Gut Microbiome Viewed across Age and Geography. Nature. 2012;486:222–227. doi: 10.1038/nature11053.
- 49.Cani P.D., de Vos W.M. Next-Generation Beneficial Microbes: The Case of Akkermansia Muciniphila. Front. Microbiol. 2017;8:1765. doi: 10.3389/fmicb.2017.01765.
- 50.Payahoo L., Khajebishak Y., Alivand M.R., Soleimanzade H., Alipour S., Barzegari A., Ostadrahimi A. Investigation the Effect of Oleoylethanolamide Supplementation on the Abundance of Akkermansia Muciniphila Bacterium and the Dietary Intakes in People with Obesity: A Randomized Clinical Trial. Appetite. 2019;141:104301. doi: 10.1016/j.appet.2019.05.032.
- 51.Everard A., Belzer C., Geurts L., Ouwerkerk J.P., Druart C., Bindels L.B., Guiot Y., Derrien M., Muccioli G.G., Delzenne N.M., et al. Cross-Talk between Akkermansia Muciniphila and Intestinal Epithelium Controls Diet-Induced Obesity. Proc. Natl. Acad. Sci. USA. 2013;110:9066–9071. doi: 10.1073/pnas.1219451110.
- 52.Ferreira-Halder C.V., Faria A.V.d.S., Andrade S.S. Action and Function of Faecalibacterium Prausnitzii in Health and Disease. Best Pract. Res. Clin. Gastroenterol. 2017;31:643–648. doi: 10.1016/j.bpg.2017.09.011.
- 53.Burger-van Paassen N., Vincent A., Puiman P.J., van der Sluis M., Bouma J., Boehm G., van Goudoever J.B., van Seuningen I., Renes I.B. The Regulation of Intestinal Mucin MUC2 Expression by Short-Chain Fatty Acids: Implications for Epithelial Protection. Biochem. J. 2009;420:211–219. doi: 10.1042/BJ20082222.
- 54.Tolhurst G., Heffron H., Lam Y.S., Parker H.E., Habib A.M., Diakogiannaki E., Cameron J., Grosse J., Reimann F., Gribble F.M. Short-Chain Fatty Acids Stimulate Glucagon-Like Peptide-1 Secretion via the G-Protein–Coupled Receptor FFAR2. Diabetes. 2012;61:364–371. doi: 10.2337/db11-1019.
- 55.Lacroix S., Pechereau F., Leblanc N., Boubertakh B., Houde A., Martin C., Flamand N., Silvestri C., Raymond F., Di Marzo V., et al. Rapid and Concomitant Gut Microbiota and Endocannabinoidome Response to Diet-Induced Obesity in Mice. mSystems. 2019;4:1–14. doi: 10.1128/mSystems.00407-19.
- 56.Ju T., Kong J.Y., Stothard P., Willing B.P. Defining the Role of Parasutterella, a Previously Uncharacterized Member of the Core Gut Microbiota. ISME J. 2019;13:1520–1534. doi: 10.1038/s41396-019-0364-5.
- 57.Lee J.S., Wang R.X., Alexeev E.E., Lanis J.M., Battista K.D., Glover L.E., Colgan S.P. Hypoxanthine Is a Checkpoint Stress Metabolite in Colonic Epithelial Energy Modulation and Barrier Function. J. Biol. Chem. 2018;293:6039–6051. doi: 10.1074/jbc.RA117.000269.
- 58.Etxeberria U., Hijona E., Aguirre L., Milagro F.I., Bujanda L., Rimando A.M., Martínez J.A., Portillo M.P. Pterostilbene-Induced Changes in Gut Microbiota Composition in Relation to Obesity. Mol. Nutr. Food Res. 2017;61:1500906. doi: 10.1002/mnfr.201500906.
- 59.Vijay A., Kouraki A., Gohir S., Turnbull J., Kelly A., Chapman V., Barrett D.A., Bulsiewicz W.J., Valdes A.M. The Anti-Inflammatory Effect of Bacterial Short Chain Fatty Acids Is Partially Mediated by Endocannabinoids. Gut Microbes. 2021;13:1997559. doi: 10.1080/19490976.2021.1997559.
- 60.Lassenius M.I., Pietiläinen K.H., Kaartinen K., Pussinen P.J., Syrjänen J., Forsblom C., Pörsti I., Rissanen A., Kaprio J., Mustonen J., et al. Bacterial Endotoxin Activity in Human Serum Is Associated with Dyslipidemia, Insulin Resistance, Obesity, and Chronic Inflammation. Diabetes Care. 2011;34:1809–1815. doi: 10.2337/dc10-2197.
- 61.Bariani M.V., Domínguez Rubio A.P., Cella M., Burdet J., Franchi A.M., Aisemberg J. Role of the Endocannabinoid System in the Mechanisms Involved in the LPS-Induced Preterm Labor. Reproduction. 2015;150:463–472. doi: 10.1530/REP-15-0211.
Associated Data
Supplementary Materials
Data Availability Statement
The data that support the findings of this study are available from the corresponding author upon reasonable request.