- Research article
- Open Access
Dietary lignans, plasma enterolactone levels, and metabolic risk in men: exploring the role of the gut microbiome
BMC Microbiology volume 22, Article number: 82 (2022)
The conversion of plant lignans to bioactive enterolignans in the gastrointestinal tract is mediated through microbial processing. The goal of this study was to examine the relationships between lignan intake, plasma enterolactone concentrations, gut microbiome composition, and metabolic risk in free-living male adults.
In 303 men participating in the Men’s Lifestyle Validation Study (MLVS), lignan intake was assessed using two sets of 7-day diet records, and gut microbiome was profiled through shotgun sequencing of up to 2 pairs of fecal samples (n = 911). A score was calculated to summarize the abundance of bacteria species that were significantly associated with plasma enterolactone levels. Of the 138 filtered species, plasma enterolactone levels were significantly associated with the relative abundances of 18 species at FDR < 0.05 level. Per SD increment of lignan intake was associated with 20.7 nM (SEM: 2.3 nM) higher enterolactone concentrations among participants with a higher species score, whereas the corresponding estimate was 4.0 nM (SEM: 1.7 nM) among participants with a lower species score (P for interaction < 0.001). A total of 12 plasma metabolites were also significantly associated with these enterolactone-predicting species. Of the association between lignan intake and metabolic risk, 19.8% (95%CI: 7.3%-43.6%) was explained by the species score alone, 54.5% (95%CI: 21.8%-83.7%) by both species score and enterolactone levels, and 79.8% (95%CI: 17.7%-98.6%) by further considering the 12 plasma metabolites.
We identified multiple gut bacteria species that were enriched or depleted at higher plasma levels of enterolactone in men. These species jointly modified the associations of lignan intake with plasma enterolactone levels and explained the majority of association between lignan intake and metabolic risk along with enterolactone levels and certain plasma metabolites.
Lignans are a class of secondary plant metabolites that are part of the large polyphenol family . Lignans are particularly rich in flax and sesame seeds, although whole grains, other seeds, legumes, coffee, and wine, account for the primary lignan intake in the general population . Both animal [3, 4] and human experiments [5,6,7] have demonstrated potentially beneficial effects of lignan intake on human health, yet the biological mechanisms are not fully understood. It is well elucidated that plant lignans are converted by gut microbiota to produce more bioactive enterolignans [8,9,10,11], which are readily absorbed into the circulation and exert subsequent health effects. In this process, the biosynthetic capacity of the gut microbiota to metabolize plant lignans[12, 13] includes O-deglycosylation of plant lignans to the aglycone form of lignans and the subsequent O-demythylation, dehydroxylation, and dehydrogenation of aglycones to produce enterolactone and/or enterodiol . Several gut microbiota species involved in lignan metabolism, including Ruminococcus bromii, Bacteroides ovatus, and Eggerthella lenta, were identified through in vitro experiments [15,16,17,18,19,20].
Early culture-based studies, although critical to establish bacteria-mediated capacity, may have not considered all microbes that contribute to the production of enterolignans due to the known difficulties in culturing many of the microbes comprising the human gastrointestinal microbiome. The use of metagenomics sequencing to profile the human gut microbiome allows for the identification of microbial taxonomic and functional features that are associated with enterolignan production in free-living individuals, although such data are sparse. In addition, it is well elucidated that there are significant between-person variabilities in producing enterolignans upon the consumption of same amount of plant lignans, which might be explained by human gut microbiome [13, 21,22,23], although no studies have been conducted to evaluate whether the human gut microbiome modulates the association between lignan intake and enterolignan levels, as well as subsequent cardiometabolic risk. To fill these knowledge gaps, we examined the interrelationships between gut microbiome interrogated using shotgun sequencing, plasma enterolactone concentrations, and cardiometabolic risk factors in men participating in the Men’s Lifestyle Validation Study (MLVS). Our primary hypothesis is that there is a consortium of human gut microbes that predict the levels of enterolignans and these microbes modulate the associations between lignan intake and enterolignan levels. In a secondary analysis, we explored whether enterolactone-predicting species also are associated with certain plasma metabolites that might further explain the association between lignan intake and cardiometabolic risk.
The average intake of total lignans was 3633 ug/day, which were primarily secoisolariciresinol (2957 ug/day), followed by pinoresinol (320 ug/day), lariciresinol (301 ug/day), and matairesinol (32 ug/day). The average plasma enterolactone levels were 29.3 nM (Table 1). The basic characteristics, including intake of energy and alcohol, smoking status, use of antibiotics, and intake of probiotics, were similar across different groups of dietary lignans and plasma enterolactone levels (Table 1).
Gut microbiota taxonomic profile and plasma enterolactone levels
PERMANOVA test based on Bray–Curtis dissimilarities showed that the overall structural variation of the gut microbiome was weakly but significantly associated with plasma enterolactone levels (R2 = 0.01, P < 0.001,) (Supplementary file Figure S1). After multivariate adjustment of covariates, a total of 18 classified species were significantly associated with enterolactone levels at FDR < 0.05 level (Table 2, Supplementary file Figure S2). For example, the relative abundance of species such as Coprococcus sp.ART55/1, Faecalibacterium prausnitzii, Alistipes shahii, Butyrivibrio crossotus, and Methanobrevibacter smithii, was significantly associated with higher enterolactone levels, and inverse associations were observed for Bacteroides dorei, B. fragilis, Clostridium bolteae, C. leptum, C. symbiosum, Lachnospiraceae bacterium.1.4.56FAA, and Ruminococcus. sp. 5.1.39BFAA.
Interactions between microbiome and lignan intake on plasma enterolactone levels
To summarize the microbial composition that was significantly associated with enterolactone levels, we built a rank-based, un-weighted species score based on the relative abundances or absence/presence of all 18 species. The multivariate-adjusted mean value of enterolactone was 15.6 nM (SEE: 5.0) among participants in the lowest decile of species score and 51.8 nM (SEE: 5.4) among participants in the top decile (P for trend < 0.0001, Fig. 1A). Lignan intake significantly interacted with the species score on plasma levels of enterolignans: dietary lignan (per SD) was associated with 4.0 nM higher enterolactone levels among participants with a lower species score, however the increment was 20.7 nM per SD lignan intake among participants with a higher species score (P for interaction < 0.001, Fig. 1B). The multivariate-adjusted means of enterolactone levels were 15.3 nM, 24.7 nM, 27.9 nM, and 47.1 nM, respectively, among participants with low lignan-low species score, high lignan-low species score, low lignan-high species, and participants with high lignan-high species score (Fig. 1C).
Lignan intake, plasma enterolactone levels, and metabolic risk
Dietary lignan was significantly associated with a lower level of BMI, HbA1c, CRP, and the overall metabolic risk score, and these associations were mainly mediated by both the species score and plasma enterolactone levels (Fig. 2). For example, lignan intake (per SD) was inversely associated with the metabolic risk score (β = -0.15 SD, SEE: 0.05, Ptrend = 0.003), 19.8% (95% CI: 7.3%, 43.6%) of which was explained by the species score and 54.5% (95% CI: 21.8%, 83.7%) jointly by the species score and plasma enterolactone levels. We did not detect any significant interaction between lignan intake and the species score on metabolic risk score or individual metabolic risk factors (Fig. 2).
Microbial compositions potentials of super pathways and plasma enterolactone levels
Overall, of the 110 super pathways considered in the analyses after filtering, the DNA abundance of 55 pathways was significantly associated with plasma enterolactone levels (FDR < 0.05), including 13 enriched and 42 depleted super pathways (Supplementary file Table S1), such as enriched branched amino acid biosynthesis, 5-aminoimidazole ribonucleotide biosynthesis, L-isoleucine biosynthesis I, L-serine, glycine biosynthesis I, hexuronide and hexuronate degradation and β-D-glucuronide and D-glucuronate degradation. We further examined the functional profile at the enzyme level, enterolactone was significantly associated with the relative abundance of 229 enzymes (FDR < 0.05, Supplemental Table S2), including 133 depleted enzymes, such as endo-1,3(4)-beta-glucanase (EC 126.96.36.199), and n-acetylgalactosamine-4-sulfatase (EC 188.8.131.52), and 96 enriched enzymes, such as glycogen phosphorylase (EC 184.108.40.206), 1,2-diacylglycerol 3-glucosyltransferase (EC 220.127.116.11), and glucosylceramidase (EC 18.104.22.168).
Plasma metabolites, microbiome, and metabolic risk
Of the 190 metabolites included in the analysis, 12 metabolites, including cinnamoylglycine, hippurate, pipecolic acid, C18:1 lysophosphatidylcholine (LPC), C18:2 LPC, C18:2 lysophosphatidylethanolamine, 2-aminohippuric acid, 5-acetylamino-6-amino-3-methyluracil, hydroxycotinine, N-acetylleucine, trigonelline and glycodeoxycholate/glycoche, exhibited significant associations with the relative abundance of several enterolactone-predicting species, such as M. smithii; B. crossotus, C. symbiosum, and L. bacterium 1.4.56FAA (Fig. 3, Supplementary file Table S3). The 12 metabolites, together with plasma enterolactone and the stool microbiome species score, jointly explained 79.8% (95%CI: 17.7%, 98.6%; P = 0.002) of the association between lignan intake and the metabolic risk score.
In the present study, we aimed to examine the role of the microbiome in the inter-connections between lignan intake, plasma enterolactone concentrations, and metabolic health in free-living healthy men. We identified multiple gut microbiota species and multiple microbial functional metabolic pathways that were significantly associated with plasma enterolactone concentrations. The association between dietary lignans and plasma enterolactone concentrations was significantly stronger among participants with a microbiota taxonomic profile that predicted higher enterolactone concentrations. In addition, the favorable association of dietary lignans with metabolic risk was partially mediated by enterolactone levels and the microbiome profile jointly. Lastly, the species that predicted enterolactone concentrations also predicted levels of certain plasma metabolites, such as cinnamoylglycine, glycine, and C18:2 LPC, which further explained the association between lignan intake and metabolic risk.
To date, no single bacteria species or strains that could independently convert plant lignans to enterolignans have been identified . Results from recent studies [16, 25] suggest that the metabolism of lignans involves phylogenetically diverse bacteria and multiple species from distinct bacterial phyla to produce bioactive enterolignans. Conventional culture-based studies are currently the main source of our knowledge on gut bacteria that may mediate these reactions. For example, Clavel et al. [13,14,15,16,17] demonstrated that C. saccharogumia, C. cocleatum, C. ramosum, B. fragilis, B. ovatus, and B. distasonis were able to deglycosylate plant lignans, whereas Butyribacterium methylotrophicum, Eubacterium callanderi, E. limosum, E. callanderi, Blautia product, and Butyribacterium methylotrophicum were capable of demethylation. Other bacteria, such as R. bromii, R. lactaris, Bifidobacterium strains, E. lenta, Enteroccus faecalis, C. scindens, Lactonifactor longoviformis, Peptostreptococcus productus, C. coccoides, Ruminococcus strains, Eubacterium strains, and B. pseudocatenulatum, can mediate these reactions as well [16,17,18,19,20, 26,27,28,29,30]. However, it is unknown whether these findings can be entirely generalizable to epidemiological study settings.
In the current study among free-living male adults, we replicated certain previously culture-based studies identified species and genera, which are known to be involved in the conversion of plant lignans to enterolactone . Specifically, B. ovatus is known to deglycosylate plant lignans [15,16,17,18,19,20]. M. smithii, the dominant archaea in the human intestinal ecosystem, can enhance the bacterial digestion of dietary polysaccharides and is critical for the syntrophic hydrogen metabolism, which is important to enterolactone production [31, 32]. Meanwhile, through interrogating the metagenomic profiling data, we also identified a wide range of novel species that were significantly associated with enterolactone levels. In particular, F. prausnitzii, one of the most abundant bacterial species in human gut, is well known as an anaerobic fiber metabolizer, and for its implication as a biomarker for diagnostics and prognostics of gastrointestinal diseases, such as ulcerative colitis and Crohn’s disease [24, 33, 34]. Our findings echo the notion that the enterolignan production involves a consortium of phylogenetically-distinct bacteria and also highlight the potentials of using modern metagenomic profiling to aid in the discovery of novel taxonomic features that are associated with biomarkers or disease outcomes. Nonetheless, our findings warrant further investigations and replications in future studies, especially in different populations that include females and broader race/ethnicity mix.
It is well-documented that the production of enterolignans from lignan intake is highly individualized, and such individualized responses might be largely ascribed to the microbial composition [21, 35]. Our study highlights, for the first time, that the gut microbial species that predict plasma enterolactone levels also significantly modulate the relationship between lignan intake and plasma enterolactone levels. Specifically, we observed a much stronger association between lignan intake and enterolactone levels when the microbiota was enriched with species that predicted higher enterolactone levels and/or depleted with species that predicted lower levels. This finding has a clear implication for developing a precision nutrition approach to enhancing the production of bioactive microbiota metabolites, such as enterolignans, through modulating the microbial composition, which may subsequently lead to improved metabolic health.
Of note, the greater production of enterolactone did not translate to a stronger inverse association between lignan intake and metabolic risk. Nonetheless, in the current study, both gut microbial profiles and plasma enterolactone levels accounted for significant proportion of the inverse association between lignan intake and metabolic risk. Existing literature also suggests potential beneficial effects of higher enterolactone levels or lignan intake on cardiometabolic conditions, including coronary heart disease and type 2 diabetes [36,37,38,39,40]. The mechanisms underlying lignans and enterolignans’ health benefits may be through inhibiting lipid peroxidation [41,42,43], reducing oxygen species production [44, 45], inducing gene expression of antioxidant enzymes , and reducing vitamin E catabolism . The structural similarity of enterolactone to 17β-estradiol allows enterolignans to bind to estrogen receptor alpha (ERα) and exert weak estrogenic or anti-estrogenic effects . Lignans also increase the levels of sex-hormone-binding protein [49, 50], which leads to reduced free estradiol, improved insulin resistance [51, 52], and a lower diabetes risk. Lignans may also improve insulin resistance through inhibiting pancreatic α-amylase  and decreasing inflammation .
It is not surprising that we found a few pathways of housekeeping genes, such as amino acid synthesis and nucleotide synthesis pathways, that were associated with enterolactone levels given that multiple species were significantly associated with enterolactone levels. Interestingly, several of the enterolactone-predicting species were also associated with plasma levels of other metabolites, such as hippurate and cinnamoylglycine, both of which are microbiota-dependent glycine-conjugated metabolites [54, 55]. Our finding is consistent with a previous cross-sectional analysis that showed a significant association between levels of enterolactone and hippuric acid in plasma . In addition, another cross-sectional analysis in the TwinsUK cohort showed significant inter-relationships between diet, F. prausnitzii abundance, plasma hippurate levels, and metabolic syndrome . A more recent analysis in the same cohort showed that plasma levels of cinnamoylglycine was positively associated with the microbiome α diversity and inversely associated with glucose level and obesity . Studies reported significant associations between total or specific LPC concentrations with lower metabolic and diabetes risk [59, 60], though the role of LPCs in metabolic pathophysiology is not completely understood and the role of microbiota is not clear. Enzyme level analysis provided further supports of the enterolactone related species at function level, such as glucosyltransferases that played an important role in the lignan metabolism [61, 62], was significantly associated with enterolactone level in our study. Overall, this exploratory analysis demonstrated novel, inter-relationships between microbiome and plasma metabolites that might be relevant to enterolactone metabolism, although longitudinal studies are needed to establish temporal relationships with metabolic risk.
The current analysis leveraged repeated assessments of diet, gut microbiome, and metabolic risk markers, which reduces error due to within-person variability. The primary limitation of this study was the cross-sectional design and observational nature of the data, which limited causal interpretations of the findings. Our findings shall be replicated in future epidemiological studies and culture-based in vitro experiments. In addition, since our study participants were older male health professionals who are mostly white, whether the findings can be generalized to women, younger population, individuals of other racial/ethnical groups, or populations with different demographic characteristics warrant further investigations.
In conclusion, in a group of healthy men we identified multiple gut microbiota species that were significantly associated with plasma enterolactone levels. A microbial profile defined by the species was significantly associated with plasma enterolactone concentrations and also significantly strengthened the relationship between plant lignan intake and plasma enterolactone levels. Both the enterolactone and the microbial profile that predicted higher enterolactone levels mediated association between higher lignan intake and lower metabolic risk. Plasma metabolites that were associated with the enterolactone-predicting species further explained the association of lignan intake with metabolic risk. These findings underscore a pivotal role of microbiome in modulating and mediating inter-relationships between diet, gut microbiota-derived metabolites, and metabolic risk.
The MLVS is a sub-study conducted in 2012–2013 within the Health Professionals Follow-up Study (HPFS) cohort with a goal of examining the validity of self-reported dietary and lifestyle assessments The HPFS was established in 1986, when 51,529 male US health professionals aged 40–75 years completed a mailed questionnaire about their medical history and lifestyle at baseline, and follow-up questionnaires have been administered to assess and update lifestyle, diet, and medical history.
Briefly, HPFS participants who provided blood samples during 1994–1995, participated in the 2010 HPFS follow-up survey with a valid food frequency questionnaire (FFQ), and were free of cardiovascular disease, cancer, and major neurological diseases, were eligible to participate in the MLVS. A total of 648 eligible HPFS men participated in the MLVS and accomplished multiple examinations, including dietary assessments, collection of blood and fecal samples, and assessments of lifestyle factors within a year. In the MLVS and HPFS, demographics, anthropometry, diet, lifestyle, medical history, and other data were collected/updated during follow-up. Covariates, including total energy intake, physical activity, smoking, alcohol consumption, BMI at age 21, and other variables, were derived in the MLVS and HPFS. The current analysis was restricted to 303 MLVS participants who provided fecal samples, completed 7-day diet records (7DDRs), provided blood sample for the measurement of cardiometabolic biomarkers, and had valid enterolactone measurements (Supplemental file Figure S3).
The study protocol was approved by the Harvard T.H. Chan School of Public Health Institutional Review Board, and informed consent was obtained from all participants.
Fecal sample collection
A total of 308 MLVS participants provided up to two pairs of stool samples. These participants were on average 71 years old (range: 65–82) at fecal sample collections. In addition, two blood samples were collected 6 months apart roughly at the same time of the fecal sample collections. Fecal sample collections were self-administered through deposition in a commode specimen collection system (Fisher Scientific). Participants were requested to provide two pairs of fecal samples (6 months apart) from two consecutive bowel movements and instructed to put a small sample of stool in collection tubes with RNAlater preservative. The time interval between two consecutive fecal sample collections was 24–72 h. Participants then shipped the tubes to the cohort biorepository within 24 h via prepaid overnight FedEx. A questionnaire was administered to collect information at fecal sample collections, such as the Bristol Stool Chart, use of antibiotic medications, age, and intake of any probiotics (except yogurt) in the preceding 2 months.
In the MLVS, diet was assessed twice via 7DDRs. Briefly, participants received a food scale (Escali Corporation, Burnsville, Minnesota) and ruler, an instructional DVD, and instructions via telephone. Participants measured and reported weights (gram) of foods before and after eating so actual intake could be computed. The participants also provided recipes of all home-prepared foods, including the number of servings in each recipe and the portion of the recipe they consumed. The food intake data were then combined with the Nutrition Data System for Research 2011 database to calculate the dietary intakes of total lignans as well as 4 individual lignans, including secoisolariciresinol, matairesinol, lariciresinol, and pinoresinol.
Taxonomic and functional profiling of metagenomic samples
Detailed information of sample collection and immediate ex situ conservation of metagenomic data, laboratory handling, and paired-end (100 × 100 nucleotides (nt)) shotgun sequencing of DNA can be found in previous publications of the MLVS [63, 64]. Standard protocols were applied to extract DNA from fecal samples. The Nextera XT DNA Library Preparation Kit was subsequently used to prepare DNA samples for sequencing. DNA was then sequenced to a target depth of 1-2Gnt each. The bioBakery workflow was used to generate taxonomic and functional profiles . Briefly, quality controls included depletion of duplicate reads using KneadData (http://huttenhower.sph.harvard.edu/kneaddata), the removal of human sequences, functional profiling by HUMAnN2 , and taxonomic profiling by MetaPhlAn2 . DNA reads were assigned to UniRef90 gene families by HUMAnN2; then the characterized gene families were assigned to MetaCyc pathways , as described in detail elsewhere . The current analysis included 911 metagenomes from the 303 MLVS participants. We further filtered out taxonomic features with a relative abundance less than 10−4 in greater than 10% of all samples. We filtered all gene families with a relative abundance less than 10−5 in greater than 10% of all samples.
Assessment of enterolactone and metabolic risk factors
For measurements of enterolactone and biomarkers, repeated fasting blood samples were collected 6 months apart through venipuncture into sodium heparin tubes, shipped by overnight mail with an ice pack, and then were processed to separate the plasma upon arrival at the cohort biorepository. Electrospray ionization orbitrap liquid chromatography mass spectrometry (ESI-LCMS; model Q-Exactive, Thermo Scientific Inc., Waltham, MA) in negative mode148 was used to measure enterolactone concentrations in plasma .
Detailed measurement and assessment of metabolic risk score is described previously . In brief, plasma levels of high-density lipoprotein (HDL)-cholesterol and triglycerides were assayed using enzymatic methods. Hemoglobin A1c (HbA1c) levels were measured using the turbidimetric immunoinhibition (Roche Diagnostics). Blinded quality control samples (10%) were interspersed among participants’ samples. Intra-assay coefficients of variation (CVs) were estimated to be < 7% for all plasma markers, except enterolactone for which the CV was 17%. Between-batch variations of enterolactone measurements were corrected using a standardized batch analytic method 
Four of the five individual components of metabolic syndrome , including triglycerides, HDL-C, body mass index (BMI), and HbA1c, were used to derive a metabolic risk score in the current analysis. Blood pressure or fasting glucose assessments were not measured in the MLVS. We first assigned a score of 1–5 to quintiles of each of the 4 factors, with 1 assigned to the lowest quintile and 5 to the highest quintile for triglycerides, BMI and HbA1c. The scoring algorithm was reversed for HDL-C that 1 was assigned for the highest quintile and 5 for the lowest quintile. The metabolic risk score was the sum of the 4 components scores with a theoretical range of 5 (lowest possible score, healthy) to 20 (highest possible score, unhealthy).
Plasma metabolomics profiling
Plasma metabolome was profiled using high-throughput liquid chromatography-mass spectrometry (LC–MS) techniques at the Broad Institute of MIT and Harvard . We used hydrophilic interaction liquid chromatography with positive ionization mode mass spectrometry detection (HILIC-pos) to separate polar metabolites. Raw data were processed using TraceFinder software (Thermo Fisher Scientific) and Progenesis QI (Nonlinear Dynamics, Newcastle upon Tyne, UK). Known metabolite identities were confirmed using authentic reference standards. A total of 201 known metabolites were profiled. Metabolites with missing values in more than 25% of the participant samples were excluded from the analysis (n = 11), leaving 190 metabolites. The rest of metabolites that were below level of detection were replaced using 1/2 minimum value of detected metabolites.
The analyses involving the relative abundance of microbial species and super pathways (DNA) were based on the 911 metagenome assessments. We matched each of the microbiome measurements to biomarkers, diet, and other variables collected at the time closest to the time of fecal sample collection. For example, the 1st and 2nd metagenomes (i.e., the first pair of fecal samples) were each matched to biomarkers measured at the first blood collection, and average intake of lignans from the first set of 7DDRs. The same strategy was used to match the second pair of fecal metagenomes with corresponding diet/biomarker measurements.
We applied normalized via arc-sin square root transformation for all species/pathway data and then used MaAslin 2 (https://huttenhower.sph.harvard.edu/maaslin2) to examine the associations between the enterolactone concentrations and the relative abundance of taxonomy and pathways. Within-person correlations between the up to four data points per person was accounted with a random effect built in MaAslin 2. We used Graphical Phylogenetic Analysis (GraPhlAn) (https://huttenhower.sph.harvard.edu/graphlan) to visualize the results of taxonomies. Covariates considered in multivariate analyses included age, physical activity, BMI at age 21, alcohol consumption, smoking, Bristol Stool Chart categories, use of antibiotics in preceding year, intake of any probiotics (except yogurt) in the preceding 2 months, and total energy intake. P values below 0.05 after false discovery rate (FDR) correction following the Benjamini–Hochberg method were considered statistically significant.
To further alleviate multiple comparison concerns, we developed an un-weighted score  to summarize the abundance of species that were significantly associated with enterolactone levels in the primary analyses at FDR < 0.05. To derive the species score, for species that were detected in ≥ 455 (of the 911) samples, we categorized them as being either “high” or “low” by the median relative abundant value. For species that were detected in < 455 (of the 911) samples, the binary categories were based on the presence or absence of the species. We then assigned 1 for higher abundance or presence of species that were positively associated with enterolactone or 0 otherwise, and we reversed the score for species that were inversely associated with enterolactone levels. We then summed the scores across all species to derive an overall species score.
To examine potential interactions between lignan intake and the gut microbiome on plasma levels of enterolactone, we defined the low and high species score using the median value as the cutoff point. The interaction between dietary lignans and species score on enterolactone was tested by further including the interaction term of dietary lignans and species score in a multivariate adjusted linear mixed model that included both dietary lignan intake and species score as well as the covariates listed above. To examine the extent to which plasma enterolactone levels and the species score explained the association between dietary lignans and the metabolic risk score, we estimated the magnitude of change in the regression coefficient for dietary lignans with and without adjustment for species score and/or plasma enterolactone level. The percentage of association between dietary lignans and metabolic risk that was explained by the species score or enterolactone was computed as follows: (1 –(β1 / β2)) X 100% using the %MEDIATE SAS macro (publicly available at www.hsph.harvard.edu/faculty/spiegelman/mediate.html), where β2 was the regression coefficient of dietary lignans without adjustment of the species score and enterolactones and β1 was the coefficient with such an adjustment. Similarly, in these analyses the averaged species score from the first and second metagenomic assessments were aligned with the first dietary and metabolic risk score assessments, while the averaged species score from the third and fourth metagenomic assessments were matched to the second dietary and metabolic risk score measurements.
In the secondary analysis that explored the role of other plasma metabolites in the inter-relationship between lignan intake, microbiome, and metabolic risk, all plasma metabolites were normalized after log-transformation. This analysis was focused on the metabolites that were associated with the species that significantly predicted enterolactone levels. Partial Spearman correlation coefficients adjusted for the same covariates listed above were used to evaluate the associations between the species and the metabolites. The potential mediational role of the metabolites in the association between dietary lignans and metabolic risk was examined. The correlations were visualized using Cytoscape (https://cytoscape.org/). The significant level was set at P values below 0.05 after false discovery rate (FDR) correction, which corresponds to partial Spearman correlation coefficients above 0.2 or below -0.2.
Availability of data and materials
Sequence data have been deposited in the Sequence Read Archive under BioProject accession PRJNA354235. The datasets used and analyzed during the current study are available from the corresponding author on reasonable request with standard internal approval procedure. Access is restricted due to participant confidentiality and privacy concerns. Further information including the procedures to obtain and access data from the Health Professionals Follow-up Study is described at https://sites.sph.harvard.edu/hpfs/resources-for-hpfs-investigators/.
Barton D, Nakanishi K, Meth-Cohn O, editors. Comprehensive natural products chemistry. Amsterdam: Elsevier; 1999.
Tham DM, Gardner CD, Haskell WL. Clinical review 97: Potential health benefits of dietary phytoestrogens: a review of the clinical, epidemiological, and mechanistic evidence. J Clin Endocrinol Metab. 1998;83(7):2223–35.
Prasad K. Reduction of serum cholesterol and hypercholesterolemic atherosclerosis in rabbits by secoisolariciresinol diglucoside isolated from flaxseed. Circulation. 1999;99(10):1355–62.
Prasad K. Regression of hypercholesterolemic atherosclerosis in rabbits by secoisolariciresinol diglucoside isolated from flaxseed. Atherosclerosis. 2008;197(1):34–42.
Pan A, Yu D, Demark-Wahnefried W, Franco OH, Lin X. Meta-analysis of the effects of flaxseed interventions on blood lipids. Am J Clin Nutr. 2009;90(2):288–97.
Penalvo JL, Lopez-Romero P. Urinary enterolignan concentrations are positively associated with serum HDL cholesterol and negatively associated with serum triglycerides in U.S. adults. J Nutr. 2012;142(4):751–6.
Tuohy K, Del Rio D, editors. Diet-Microbe interactions in the gut: effects on human health and disease. Italy: Elsevier Inc.; 2015.
Bowey E, Adlercreutz H, Rowland I. Metabolism of isoflavones and lignans by the gut microflora: a study in germ-free and human flora associated rats. Food Chem Toxicol. 2003;41(5):631–6.
Eeckhaut E, Struijs K, Possemiers S, Vincken JP, Keukeleire DD, Verstraete W. Metabolism of the lignan macromolecule into enterolignans in the gastrointestinal lumen as determined in the simulator of the human intestinal microbial ecosystem. J Agric Food Chem. 2008;56(12):4806–12.
Pettersson D, Aman P, Knudsen KE, Lundin E, Zhang JX, Hallmans G, et al. Intake of rye bread ileostomists increases ileal excretion of fiber polysaccharide components and organic acids but does not increase plasma or urine lignans and isoflavonoids. J Nutr. 1996;126(6):1594–600.
Setchell KD, Lawson AM, Borriello SP, Harkness R, Gordon H, Morgan DM, et al. Lignan formation in man–microbial involvement and possible roles in relation to cancer. Lancet. 1981;2(8236):4–7.
Bach Knudsen KE, Serena A, Kjaer AK, Tetens I, Heinonen SM, Nurmi T, et al. Rye bread in the diet of pigs enhances the formation of enterolactone and increases its levels in plasma, urine and feces. J Nutr. 2003;133(5):1368–75.
Clavel T, Henderson G, Alpert CA, Philippe C, Rigottier-Gois L, Dore J, et al. Intestinal bacterial communities that produce active estrogen-like compounds enterodiol and enterolactone in humans. Appl Environ Microbiol. 2005;71(10):6077–85.
Clavel T, Dore J, Blaut M. Bioavailability of lignans in human subjects. Nutr Res Rev. 2006;19(2):187–96.
Clavel T, Borrmann D, Braune A, Dore J, Blaut M. Occurrence and activity of human intestinal bacteria involved in the conversion of dietary lignans. Anaerobe. 2006;12(3):140–7.
Clavel T, Henderson G, Engst W, Dore J, Blaut M. Phylogeny of human intestinal bacteria that activate the dietary lignan secoisolariciresinol diglucoside. FEMS Microbiol Ecol. 2006;55(3):471–8.
Clavel T, Lippman R, Gavini F, Dore J, Blaut M. Clostridium saccharogumia sp. nov. and Lactonifactor longoviformis gen. nov., sp. nov., two novel human faecal bacteria involved in the conversion of the dietary phytoestrogen secoisolariciresinol diglucoside. Syst Appl Microbiol. 2007;30(1):16–26.
Jin JS, Zhao YF, Nakamura N, Akao T, Kakiuchi N, Min BS, et al. Enantioselective dehydroxylation of enterodiol and enterolactone precursors by human intestinal bacteria. Biol Pharm Bull. 2007;30(11):2113–9.
Roncaglia L, Amaretti A, Raimondi S, Leonardi A, Rossi M. Role of bifidobacteria in the activation of the lignan secoisolariciresinol diglucoside. Appl Microbiol Biotechnol. 2011;92(1):159–68.
Wang LQ, Meselhy MR, Li Y, Qin GW, Hattori M. Human intestinal bacteria capable of transforming secoisolariciresinol diglucoside to mammalian lignans, enterodiol and enterolactone. Chem Pharm Bull (Tokyo). 2000;48(11):1606–10.
Possemiers S, Bolca S, Eeckhaut E, Depypere H, Verstraete W. Metabolism of isoflavones, lignans and prenylflavonoids by intestinal bacteria: producer phenotyping and relation with intestinal community. FEMS Microbiol Ecol. 2007;61(2):372–83.
De Filippis F, Vannini L, La Storia A, Laghi L, Piombino P, Stellato G, et al. The same microbiota and a potentially discriminant metabolome in the saliva of omnivore, ovo-lacto-vegetarian and Vegan individuals. PLoS One. 2014;9(11):e112373.
Hålldin E, Eriksen AK, Brunius C, da Silva AB, Bronze M, Hanhineva K, Aura AM, Landberg R. Factors explaining interpersonal variation in plasma enterolactone concentrations in humans. Mol Nutr Food Res. 2019;63(16):e1801159.
Hullar MA, Lancaster SM, Li F, Tseng E, Beer K, Atkinson C, et al. Enterolignan-producing phenotypes are associated with increased gut microbial diversity and altered composition in premenopausal women in the United States. Cancer Epidemiol Biomarkers Prev. 2015;24(3):546–54.
Bess EN, Bisanz JE, Yarza F, Bustion A, Rich BE, Li X, et al. Genetic basis for the cooperative bioactivation of plant lignans by Eggerthella lenta and other human gut bacteria. Nat Microbiol. 2020;5(1):56–66.
Jin JS, Hattori M. Further studies on a human intestinal bacterium Ruminococcus sp. END-1 for transformation of plant lignans to mammalian lignans. J Agric Food Chem. 2009;57(16):7537–42.
Jin JS, Hattori M. Human intestinal bacterium, strain END-2 is responsible for demethylation as well as lactonization during plant lignan metabolism. Biol Pharm Bull. 2010;33(8):1443–7.
Jin JS, Kakiuchi N, Hattori M. Enantioselective oxidation of enterodiol to enterolactone by human intestinal bacteria. Biol Pharm Bull. 2007;30(11):2204–6.
Jin JS, Zhao YF, Nakamura N, Akao T, Kakiuchi N, Hattori M. Isolation and characterization of a human intestinal bacterium, Eubacterium sp. ARC-2, capable of demethylating arctigenin, in the essential metabolic process to enterolactone. Biol Pharm Bull. 2007;30(5):904–11.
Lagkouvardos I, Klaring K, Heinzmann SS, Platz S, Scholz B, Engel KH, et al. Gut metabolites and bacterial community networks during a pilot intervention study with flaxseeds in healthy adult men. Mol Nutr Food Res. 2015;59(8):1614–28.
Lampe JW, Martini MC, Kurzer MS, Adlercreutz H, Slavin JL. Urinary lignan and isoflavonoid excretion in premenopausal women consuming flaxseed powder. Am J Clin Nutr. 1994;60(1):122–8.
Samuel BS, Hansen EE, Manchester JK, Coutinho PM, Henrissat B, Fulton R, et al. Genomic and metabolic adaptations of Methanobrevibacter smithii to the human gut. Proc Natl Acad Sci U S A. 2007;104(25):10643–8.
Lopez-Siles M, Duncan SH, Garcia-Gil LJ, Martinez-Medina M. Faecalibacterium prausnitzii: from microbiology to diagnostics and prognostics. ISME J. 2017;11(4):841–52.
Benus RF, van der Werf TS, Welling GW, Judd PA, Taylor MA, Harmsen HJ, et al. Association between Faecalibacterium prausnitzii and dietary fibre in colonic fermentation in healthy human subjects. Br J Nutr. 2010;104(5):693–700.
Kuijsten A, Arts IC, Van’tVeer P, Hollman PC. The relative bioavailability of enterolignans in humans is enhanced by milling and crushing of flaxseed. J Nutr. 2005;135(12):2812–6.
Xu C, Liu Q, Zhang Q, Gu A, Jiang ZY. Urinary enterolactone is associated with obesity and metabolic alteration in men in the US national health and nutrition examination survey 2001–10. Br J Nutr. 2015;113(4):683–90.
Sun Q, Wedick NM, Pan A, Townsend MK, Cassidy A, Franke AA, et al. Gut microbiota metabolites of dietary lignans and risk of type 2 diabetes: a prospective investigation in two cohorts of U.S. women. Diabetes Care. 2014;37(5):1287–95.
Pan A, Sun J, Chen Y, Ye X, Li H, Yu Z, et al. Effects of a flaxseed-derived lignan supplement in type 2 diabetic patients: a randomized, double-blind, cross-over trial. PLoS One. 2007;2(11):e1148.
Cornish SM, Chilibeck PD, Paus-Jennsen L, Biem HJ, Khozani T, Senanayake V, et al. A randomized controlled trial of the effects of flaxseed lignan complex on metabolic syndrome composite score and bone mineral in older adults. Appl Physiol Nutr Metab. 2009;34(2):89–98.
Liu G, Li J, Li Y, Hu Y, Franke AA, Liang L, et al. Gut microbiota–derived metabolites and risk of coronary artery disease: a prospective study among US men and women. Am J Clin Nutr. 2021;114(1):238-47.
During A, Debouche C, Raas T, Larondelle Y. Among plant lignans, pinoresinol has the strongest antiinflammatory properties in human intestinal Caco-2 cells. J Nutr. 2012;142(10):1798–805.
Kitts DD, Yuan YV, Wijewickreme AN, Thompson LU. Antioxidant activity of the flaxseed lignan secoisolariciresinol diglycoside and its mammalian lignan metabolites enterodiol and enterolactone. Mol Cell Biochem. 1999;202(1–2):91–100.
Prasad K. Hydroxyl radical-scavenging property of secoisolariciresinol diglucoside (SDG) isolated from flax-seed. Mol Cell Biochem. 1997;168(1–2):117–23.
Lee B, Kim KH, Jung HJ, Kwon HJ. Matairesinol inhibits angiogenesis via suppression of mitochondrial reactive oxygen species. Biochem Biophys Res Commun. 2012;421(1):76–80.
Prasad K. Antioxidant activity of secoisolariciresinol diglucoside-derived metabolites, secoisolariciresinol, enterodiol, and enterolactone. Int J Angiol. 2000;9(4):220–5.
Cortes C, Palin MF, Gagnon N, Benchaar C, Lacasse P, Petit HV. Mammary gene expression and activity of antioxidant enzymes and concentration of the mammalian lignan enterolactone in milk and plasma of dairy cows fed flax lignans and infused with flax oil in the abomasum. Br J Nutr. 2012;108(8):1390–8.
Hanzawa F, Nomura S, Sakuma E, Uchida T, Ikeda S. Dietary sesame seed and its lignan, sesamin, increase tocopherol and phylloquinone concentrations in male rats. J Nutr. 2013;143(7):1067–73.
Yoder SC, Lancaster SM, Hullar MAJ, Lampe JW. Chapter 7: Gut microbial metabolism of plant lignans: influence on human healthhy. In Tuohy K, Del Rio D, editors. Diet-Microbe Interactions in the Gut: Effects on Human Health and Disease. Elsevier Inc.; 2014. p. 103-117. https://doi.org/10.1016/C2012-0-01316-6. Accessed 21 Mar 2022.
Martin ME, Haourigui M, Pelissero C, Benassayag C, Nunez EA. Interactions between phytoestrogens and human sex steroid binding protein. Life Sci. 1996;58(5):429–36.
Adlercreutz H, Mousavi Y, Clark J, Hockerstedt K, Hamalainen E, Wahala K, et al. Dietary phytoestrogens and cancer: in vitro and in vivo studies. J Steroid Biochem Mol Biol. 1992;41(3–8):331–7.
Winters SJ, Gogineni J, Karegar M, Scoggins C, Wunderlich CA, Baumgartner R, et al. Sex hormone-binding globulin gene expression and insulin resistance. J Clin Endocrinol Metab. 2014;99(12):E2780–8.
Ding EL, Song Y, Manson JE, Hunter DJ, Lee CC, Rifai N, et al. Sex hormone-binding globulin and risk of type 2 diabetes in women and men. N Engl J Med. 2009;361(12):1152–63.
Hano C, Renouard S, Molinie R, Corbin C, Barakzoy E, Doussot J, et al. Flaxseed (Linum usitatissimum L.) extract as well as (+)-secoisolariciresinol diglucoside and its mammalian derivatives are potent inhibitors of alpha-amylase activity. Bioorg Med Chem Lett. 2013;23(10):3007–12.
Wikoff WR, Anfora AT, Liu J, Schultz PG, Lesley SA, Peters EC, et al. Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. Proc Natl Acad Sci U S A. 2009;106(10):3698–703.
Feng YL, Cao G, Chen DQ, Vaziri ND, Chen L, Zhang J, et al. Microbiome-metabolomics reveals gut microbiota associated with glycine-conjugated metabolites and polyamine metabolism in chronic kidney disease. Cell Mol Life Sci. 2019;76(24):4961–78.
Miles FL, Navarro SL, Schwarz Y, Gu H, Djukovic D, Randolph TW, et al. Plasma metabolite abundances are associated with urinary enterolactone excretion in healthy participants on controlled diets. Food Funct. 2017;8(9):3209–18.
Pallister T, Jackson MA, Martin TC, Zierer J, Jennings A, Mohney RP, et al. Hippurate as a metabolomic marker of gut microbiome diversity: modulation by diet and relationship to metabolic syndrome. Sci Rep. 2017;7(1):13670.
Menni C, Zhu J, Le Roy CI, Mompeo O, Young K, Rebholz CM, et al. Serum metabolites reflecting gut microbiome alpha diversity predict type 2 diabetes. Gut Microbes. 2020;11(6):1632–42.
Papandreou C, Bullo M, Zheng Y, Ruiz-Canela M, Yu E, Guasch-Ferre M, et al. Plasma trimethylamine-N-oxide and related metabolites are associated with type 2 diabetes risk in the Prevencion con Dieta Mediterranea (PREDIMED) trial. Am J Clin Nutr. 2018;108(1):163–73.
Bachlechner U, Floegel A, Steffen A, Prehn C, Adamski J, Pischon T, et al. Associations of anthropometric markers with serum metabolites using a targeted metabolomics approach: results of the EPIC-potsdam study. Nutr Diabetes. 2016;6:e215.
Ono E, Waki T, Oikawa D, et al. Glycoside-specific glycosyltransferases catalyze regio-selective sequential glucosylations for a sesame lignan, sesaminol triglucoside. Plant J. 2020;101(5):1221–33. https://doi.org/10.1111/tpj.14586.
Noguchi A, Fukui Y, Iuchi-Okada A, et al. Sequential glucosylation of a furofuran lignan, (+)-sesaminol, by Sesamum indicum UGT71A9 and UGT94D1 glucosyltransferases. Plant J. 2008;54(3):415–27. https://doi.org/10.1111/j.1365-313X.2008.03428.x.
Mehta RS, Abu-Ali GS, Drew DA, Lloyd-Price J, Subramanian A, Lochhead P, et al. Stability of the human faecal microbiome in a cohort of adult men. Nat Microbiol. 2018;3(3):347–55.
Abu-Ali GS, Mehta RS, Lloyd-Price J, Mallick H, Branck T, Ivey KL, et al. Metatranscriptome of human faecal microbial communities in a cohort of adult men. Nat Microbiol. 2018;3(3):356–66.
McIver LJ, Abu-Ali G, Franzosa EA, Schwager R, Morgan XC, Waldron L, et al. bioBakery: a meta’omic analysis environment. Bioinformatics. 2018;34(7):1235–7.
Truong DT, Franzosa EA, Tickle TL, Scholz M, Weingart G, Pasolli E, et al. MetaPhlAn2 for enhanced metagenomic taxonomic profiling. Nat Methods. 2015;12(10):902–3.
Abubucker S, Segata N, Goll J, Schubert AM, Izard J, Cantarel BL, et al. Metabolic reconstruction for metagenomic data and its application to the human microbiome. PLoS Comput Biol. 2012;8(6):e1002358.
Caspi R, Billington R, Ferrer L, Foerster H, Fulcher CA, Keseler IM, et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2016;44(D1):D471–80.
Franke AA, Custer LJ, Wilkens LR, Le Marchand LL, Nomura AM, Goodman MT, et al. Liquid chromatographic-photodiode array mass spectrometric analysis of dietary phytoestrogens from human urine and blood. J Chromatogr B Analyt Technol Biomed Life Sci. 2002;777(1–2):45–59.
Poole EM, Tworoger SS, Hankinson SE, Schernhammer ES, Pollak MN, Baer HJ. Body size in early life and adult levels of insulin-like growth factor 1 and insulin-like growth factor binding protein 3. Am J Epidemiol. 2011;174(6):642–51.
Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults. Executive summary of the third report of The National Cholesterol Education Program (NCEP) expert panel on detection, evaluation, and treatment of high blood cholesterol in adults (Adult Treatment Panel III). JAMA. 2001;285(19):2486–97.
Paynter NP, Balasubramanian R, Giulianini F, Wang DD, Tinker LF, Gopal S, et al. Metabolic predictors of incident coronary heart disease in women. Circulation. 2018;137(8):841–53.
Paynter NP, Chasman DI, Pare G, Buring JE, Cook NR, Miletich JP, et al. Association between a literature-based genetic risk score and cardiovascular events in women. JAMA. 2010;303(7):631–7.
We thank the participants and staffs of the HPFS and MLVS for their participation and work.
This study is supported by research grant R01HL035464 (ER, QS), DK120870 and R01DK126698 (QS), HL060712 and DK119268 (FBH), R01CA202704 (ATC, CH), K24DK098311 (ATC), and U54DE023798 (CH), K99DK119412 (DDW), and K99DK122128 (JL) from the National Institutes of Health (NIH). The Men's Lifestyle Validation Study was supported by U01CA152904 from the National Cancer Institute (NCI). The Health Professionals Follow-Up Study was supported by U01CA167552 from NCI and R01HL35464 from NHLBI. The funding source had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and the decision to submit the manuscript for publication.
Ethics approval and consent to participate
The study protocol was approved by the Harvard T.H. Chan School of Public Health Institutional Review Board. All participants provided written informed consent.
Consent for publication
Competing of interest
All authors have no support from any organization for the submitted work; no financial relationships with any organizations that might have an interest in the submitted work in the previous three years, no other relationships or activities that could appear to have influenced the submitted work. Y.L. was supported by grants from the California Walnut Commission and Swiss Reinsurance Company.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplemental Table S2. list of enzymes that significantly associated with enterolactone concentration.
Supplementary Figure S1. Association between enterolactone and microbial community (PERMANOVA with Bray-Curtis dissimilarities: R2=0.01044, P<0.001). Supplemental Figure S2. Taxonomic tree with highlighted species that significantly associated with enterolactone concentrations. Supplemental Figure S3. Flowchart of participant enrollment. Supplementary Table S1. Multivariable-adjusted associations between plasma enterolactone and relative abundances of genetic predisposition of super pathways*. Supplementary Table S3. Plasma metabolites significantly associated with enterolactone-predicting species.
About this article
Cite this article
Li, Y., Wang, F., Li, J. et al. Dietary lignans, plasma enterolactone levels, and metabolic risk in men: exploring the role of the gut microbiome. BMC Microbiol 22, 82 (2022). https://doi.org/10.1186/s12866-022-02495-0