Skip to main content

Impact of lignans in oilseed mix on gut microbiome composition and enterolignan production in younger healthy and premenopausal women: an in vitro pilot study



Dietary lignans belong to the group of phytoestrogens together with coumestans, stilbenes and isoflavones, and themselves do not exhibit oestrogen-like properties. Nonetheless, the gut microbiota converts them into enterolignans, which show chemical similarity to the human oestrogen molecule. One of the richest dietary sources of lignans are oilseeds, including flaxseed. The aim of this pilot study was to determine the concentration of the main dietary lignans in an oilseed mix, and explore the gut microbiota-dependent production of enterolignans for oestrogen substitution in young and premenopausal women. The oilseed mix was fermented in a pH-controlled batch culture system inoculated with women’s faecal samples. The lignan content and enterolignan production were measured by ultra‐high-performance liquid chromatography–tandem mass spectrometry (UHPLC–MS/MS), and the faecal-derived microbial communities were profiled by 16S rRNA gene-based next-generation sequencing.


In vitro batch culture fermentation of faecal samples inoculated with oilseed mix for 24 h resulted in a substantial increase in enterolactone production in younger women and an increase in enterodiol in the premenopausal group. As for the gut microbiota, different baseline profiles were observed as well as different temporal dynamics, mainly related to Clostridiaceae, and Klebsiella and Collinsella spp.


Despite the small sample size, our pilot study revealed that lignan-rich oilseeds could strongly influence the faecal microbiota of both younger and premenopausal females, leading to a different enterolignan profile being produced. Further studies in larger cohorts are needed to evaluate the long-term effects of lignan-rich diets on the gut microbiota and find out how enterolactone-producing bacterial species could be increased. Diets rich in lignans could potentially serve as a safe supplement of oestrogen analogues to meet the cellular needs of endogenous oestrogen and deliver numerous health benefits, provided that the premenopausal woman microbiota is capable of converting dietary precursors into enterolignans.


Dietary lignans are phytoestrogens found in plants as aglycones or glycosides [1]. Around 70 different plant species are rich in various lignans in their roots, rhizomes, stems, leaves, seeds, and fruits [2]. In particular, grains with bran and oilseeds, especially flaxseeds, are extremely rich sources [3, 4]. The most studied aglycone forms of dietary lignans are secoisolariciresinol (Seco), lariciresinol (Lari), pinoresinol (Pino) and matairesinol (Mat) (Fig. 1), whose range of concentration varies from 88.06 mg/100 g of fresh weight (FW) to 436.99 mg/100 g FW [3]. Secoisolariciresinol diglucoside (SDG), i.e. the glycosylated form of Seco, represents approximately 1% of the dry weight of flaxseeds [3], and its concentration ranges from 11.9 to 25.9 mg/g FW [5]. This has led to flaxseed becoming the most common source of enterolignan precursors for research.

Fig. 1
figure 1

Dietary lignans. Secoisolariciresinol (Seco); Matairesinol (Mat); Secoisolariciresinol Diglucoside (SDG); Lariciresinol (Lari); Pinoresinol (Pino)

Dietary lignans do not exhibit oestrogen-like properties but are known to be metabolised by the gut microbiota into enterolignans (or mammalian lignans), i.e. enterodiol (ED) and enterolactone (EL). Given that EL, more than ED, is able to bind the oestrogen receptor ER-α [6,7,8,9] and modulate endogenous oestrogen levels [10], it would be interesting to investigate the production of enterolignans by premenopausal-related gut microbiota layouts, following a lignan-rich diet. The magnitude of the modulating effects of enterolignans is expected to be determined by the biological level of endogenous oestrogen: enterolignans would exert antagonistic activity with normal oestrogen levels, while acting as weak oestrogens when their levels are low [11]. Additionally, it has been suggested that circulating enterolignans may stimulate the synthesis of sex hormone-binding globulin, which binds sex hormones, reducing their systemic levels and biological activity and, consequently, the risk of hormone-dependent cancer [12]. EL may also inhibit aromatase, an enzyme that converts androgens to oestrogens in fat, muscle and other tissues [13]. There is not much known yet about ER-β and G protein-coupled estrogen receptor (GPR), except that the deficiency of ER-β can inhibit expression of facilitative glucose transporter type 4 GLUT4 [14] and GPR is involved in body weight regulation, glucose and lipid homeostasis and inflammation (Sharma et al. [15]). After activating these receptors, oestrogen can promote mitotic activity, neoplastic transformation and cellular proliferation. The abundance of intracellular ER receptors confirms the significance of oestrogen and explains the fluctuations of the main physiological parameters during the natural female cycle as well as the dramatic health changes that follow the general decline in oestrogen after the menopause, such as an increase of cardiovascular conditions, weight gain, breast- and colon cancers, osteoarthritis and memory loss in post-menopausal women [16,17,18]. Suggesting that enterolignans are able to ligate on ER-α, substituting endogenous oestrogen and potentially alleviating the symptoms of lowered oestrogen, it would be interesting to investigate the potential role of a diet high in enterolignan precursors together with the activity of the gut microbiota in premenopausal women.

Numerous clinical trials and animal model studies have demonstrated the association between high concentrations of ED and EL in blood and urine and various health benefits, such as normal blood lipids [19], reduced risk of breast cancer [20] and osteoporosis [9, 21], as well as improved signs of metabolic syndrome [22]. In particular, a number of studies have revealed that the anti-cancer effect of enterolignans is higher than that of their precursors (i.e. dietary lignans) and that the decreased risk of cancer correlates with blood and urinary levels of enterolignans but not with lignan-rich diets [23,24,25], thus emphasizing the relevance of the gut microbiota as a key mediator of their effect.

According to in vitro studies from Peterson et al. [26], the faecal microbiota can metabolise 100% of Lari, producing 46% of EL and 54% of ED, while the microbial conversion of other lignans has been reported as incomplete (in 24 h 72% of SDG, 55% of Pino and 62% of Mat are metabolised into enterolignans, and only 21% of SDG will be converted into EL). Furthermore, some intestinal microbes have been shown to convert most ED to EL [27], but one-third of adults’ microbiotas could directly produce EL, skipping the ED stage [28]. As for gender, the female microbiota is capable of producing more ED and EL than the male one [29]. Even age has been reported to affect enterolignan production, with the child microbiota being less capable of converting dietary lignans into ED and EL [30].

Although the lignan transformation by the gut microbiota is recognized as essential in protecting against menopausal symptoms as well as certain hormone-dependent chronic diseases (e.g. cancer, cardiovascular disease and osteoporosis) [31,32,33], to date only a few bacterial species (often subdominant) have been identified as enabling such a transformation [33, 34]. In particular, EL production has been related to the abundance of Ruminococcus species, i.e. R. bromii and R. lactaris [35], as well as to those of LactobacillusEnterococcus [36] and Methanobrevibacter, an archaeon central to the syntrophic hydrogen metabolism that may be important for EL production [37]. On the other hand, low-EL excreter phenotypes have been associated with the pathobiont, sulphite-reducing bacterium, Bilophila, as well as microbial signatures of epithelial dysfunction and inflammation [38].

In an attempt to bridge this gap, providing further insights into bacterial taxa capable of converting lignans into EL/ED, here we investigated the lignan profile of a commercially prepared blend of oilseeds as a functional dietary supplement for women, by ultra-high-performance liquid chromatography–tandem mass spectrometry (UHPLC–MS/MS), and then used this as a substrate for inoculating stool samples from young and premenopausal women in in vitro anaerobic stirred batch culture systems, to simulate the physicochemical characteristics of the distal colon. These in vitro batch systems are technological platforms that can simulate the spatial, temporal and environmental features that microbes experience within the gut environment. Being host-free systems, these are the ideal systems in which to study microbial perturbations resulting from the addition of exogenous stimuli, as microbial changes can be measured in isolation of any concurrent effects on the host [39,40,41,42,43,44,45]. The aim of these experiments was to determine the production of oilseed mix-derived enterolignans and to explore differences in bacterial taxa that could be involved in lignan transformation. We investigated the dynamics of individual lignans (Seco, Lari, Pino and Mat) as well as enterolignans (ED and EL), gut microbial profiles and metabolic end products of microbial fermentation. The faecal-derived microbial communities were profiled at 0, 5 and 24 h by 16S rRNA gene-based next-generation sequencing, and the dietary lignan content was measured by UHPLC–MS/MS were profiled at 0.2, 5 and 24 h.


Sample hydrolysis, extraction and quantification

The methodology for the extraction of lignans from food samples as well as of enterolignans from fermentation faecal samples was adapted from the work of Nørskov and Knudsen [46] as well as the work of Milder et al. [47] with minor modifications. The combination of solvent-based extraction and alkaline hydrolysis was used to release lignans from food or faecal matrix. On the other hand, enzymatic extraction with β-glucuronidase/sulfatase was performed in order to degrade the glycosidic bonds and to release the lignan aglycones. Solid Phase Extraction method was used to prepare samples for UHPLC–MS/MS analysis.

A UHPLC–MS/MS method was developed to obtain high resolution and signal accuracy. The deuterated standard Secoisolariciresinol-d6 (Seco-d6) was used as the reference standard. A multiple reaction monitoring (MRM) method was set up in negative ion mode for each analyte and quantification was established using the most intense MRM signal transition with 8–16 point calibration curves of pure analytical standards. The obtained MRM detection parameters are given in Table 1.

Table 1 Multiple reaction monitoring (MRM) detection parameters

Determination of the lignan content in the oilseed mix

Obtained concentrations (mg/100 g of dry matter) of free and total (free plus bound) Seco, Lari, Pino and Mat in the oilseed mix are listed below in Table 2.

Table 2 Free and total (free plus bound) lignans in oilseed mix and baseline concentration in in vitro fermentation medium

The UHPLC–MS/MS analysis revealed that amongst free fractions, the quantities of Lari in the oilseed mix prevail (44.9 mg out of 90.5 mg per 100 g of mix), whereas amongst total lignans, the mean values for Seco are the highest (1467.6 mg out of 1960.7 mg per 100 g of mix).

Analysis of enterolignans in fermented samples

The results of batch culture fermentation experiments revealed differences in the concentration of individual lignans and enterolignans at different time points in the samples with added oilseed mix between young healthy (YD) and premenopausal (PD) donors (Table 3). Figures 2 and 3 report the values of free lignans and enterolignans respectively, assessed during the fermentation at 3 time points, 0.2 (T0.2), 5 (T5) and 24 h (T24).

Table 3 Concentration of lignans (Seco, Lari, Mat and Pino) and enterolignans (ED and EL) in the faecal samples from younger healthy donors and premenopausal donors at 0.2 (T0.2), 5 (T5) and 24 h (T24) of fermentation
Fig. 2
figure 2

Dynamics of individual dietary lignans (Seco, Lari, Mat and Pino) over 24 h of fermentation in faecal samples from younger healthy donors (YD) and premenopausal donors (PD) with added oilseed mix. Data are expressed as the mean ± SD

Fig. 3
figure 3

Concentration of enterodiol (ED) and enterolactone (EL) in faecal samples from younger healthy donors (YD) and premenopausal donors (PD) at 0.2 (T0.2), 5 (T5) and 24 h (T24) of fermentation with added oilseed mix

The average sum of the total individual lignans in the YD samples falls about 10 times during the first 5 h of fermentation, with a slight increase at 24 h. The concentration of Mat and Pino in this group decreased significantly, whereas Seco only slightly decreased and Lari increased. At the same time, the amounts of total enterolignans in the YD samples increased more than 10 times in 24 h, with EL expanding nearly 20 times but ED being slightly reduced.

In contrast, different dynamics were found in the fermentation samples from PDs: the total lignans increased almost 4 times in 24 h with only Pino falling about 10 times, while the levels of Seco, Mat and Lari expanded. Additionally, the amount of total enterolignans halved in 24 h, with ED levels expanding nearly 60 times, but EL dropping almost 100 times from 0.2 h to 24 h, suggesting the inability of faecal bacteria from PDs to convert ED to EL. It is also worth noting that all premenopausal participants entered the experiment with very different initial levels of enterolignans (particularly, ED), resulting in a high value of standard deviation.

Figure 2 shows the dynamics of Seco, Lari, Pino and Mat in the samples from YDs (a) and PDs (b) with added oilseed mix in 24-h fermentation experiments. While very different concentrations of lignans between YDs and PDs can be observed at T0.2, at T5 their levels are almost identical as well as at T24, except for Mat, whose level increased almost 500 times in the PD group.

Figure 3 illustrates the dynamics of ED and EL in YD and PD fermentation samples during the 24-h batch culture experiment. The ED levels are equally low in YD and PD samples at T0.2 and T5, while at T24, the value in the PD samples is almost 37 times higher than in the YD samples. As for EL, the average concentration is remarkably high in the PD samples at T0.2, falling slightly at T5 and experiencing a tremendous drop at T24. Conversely, the EL levels in the YD samples are insignificant at the beginning, slightly increase at T5 and are boosted at T24. Taken together, our results clearly demonstrate different dynamics of lignan transformation in different age groups of women: adding oilseed mix to YD samples leads to an overall drop in levels of dietary lignans and to the production of enterolignans. Contrary to this, in PDs, Mat levels increase after 24 h, with increased level of ED but decreased level of EL. These observations suggest that PDs may have a different, potentially altered intestinal microbiota configuration, with very poor conversion of dietary lignans and inability to efficiently convert ED or Mat to EL.

Impact of the oilseed mix on faecal-derived microbial communities from younger healthy and premenopausal women

The faecal microbial communities from YDs and PDs were profiled over time during the fermentation experiments, to assess whether the different dynamics of conversion of the lignans of the oilseed mix were potentially attributable to different layouts and trajectories of the gut microbiota. In parallel, for each woman, two other batch cultures were set up, inoculated with the well-known prebiotic compound inulin (Raftilose P95) as a positive control or without addition of any compound (i.e. the negative control), respectively. The 16S rRNA gene-based next-generation sequencing of all fermentation samples yielded a total of 1,809,764 high-quality reads, with an average of 33,514 ± 7194 sequences per sample, binned into 2623 amplicon sequence variants (ASVs) at 99% similarity.

No significant differences were observed in alpha diversity across the entire dataset, regardless of the origin of the faecal sample (PD vs. YD), experimental condition (oilseed mix vs. inulin vs. negative control) and time point (T0 vs. T5 vs. T24) (p value > 0.05, Wilcoxon test). However, it should be noted that the biodiversity of both PD and YD faecal-derived microbial communities tended to decrease over time in all experimental conditions, except for the negative controls where it was generally kept (Additional file 1: Fig. S1). With specific regard to the oilseed mix, this decreasing trend was already appreciable after 5 h of fermentation in PD samples while only after 24 h for YD samples.

The Principal Coordinates Analysis (PCoA) of inter-sample variation (i.e. beta diversity), based on unweighted (Fig. 4) and weighted (Additional file 2: Fig. S2) UniFrac distances, showed a significant separation between the faecal-derived microbial communities of PDs and YDs, regardless of the experimental condition and time point (p value < 1 × 10−4, permutation test with pseudo-F ratio). Within each group (PD and YD), the samples still segregated significantly according to both the experimental condition and the time point (p value < 0.001), suggesting a differential impact of supplements over time, likely related to the baseline microbial community. In particular, the PD faecal microbial ecosystem underwent an early modification of its overall structure, as evidenced by the shift at T5, and continued to evolve up to 24 h under all experimental conditions. The same was basically true for YD samples based on unweighted UniFrac distances, while a significant shift in the weighted UniFrac-based PCoA was observed only at T24, suggesting a dominant configuration likely to be more resilient to external perturbations compared to the PD one. Please, see Additional file 3: Table S1 for the results of adonis and ANOSIM statistics applied to weighted and unweighted UniFrac distance-based ordination.

Fig. 4
figure 4

Faecal-derived microbial communities of premenopausal women separate from those of younger healthy women in 24-h fermentation experiments in the presence of oilseed mix, inulin or without additions. Left, Principal Coordinates Analysis (PCoA) based on unweighted UniFrac distances, showing all fermentation samples, coloured by group of women (premenopausal, light blue vs. younger healthy, orange). A significant separation between groups was found, regardless of experimental condition (oilseed mix, inulin and negative control - “Ctr-”) and time point (T0, T5 and T24) (p value < 1 × 10−4, permutation test with pseudo-F ratio). Right, PCoA plots showing the fermentation samples for premenopausal women (top panel) and younger healthy women (bottom panel). Within each group of women, the samples separate significantly by experimental condition (oilseed mix, green; inulin, red; Ctr-, yellow) and time point (T0, circle; T5, square; T24, diamond) (p value < 0.001)

At T0, a number of significant taxonomic differences were indeed observed between PDs and YDs. In particular, compared to YDs, the baseline faecal-derived microbial ecosystem of PDs showed increased relative abundance of some Bacteroidetes members, especially Bacteroides and Rikenellaceae, as well as Ruminococcaceae genera (p value ≤ 0.05, Wilcoxon test). On the other hand, Coriobacteriaceae, especially Collinsella, and Streptococcus were much more represented in the baseline microbiota of YDs vs. PDs (p value ≤ 0.05) (Additional file 4: Fig. S3).

With specific regard to the impact of oilseed mix on the faecal-derived microbial communities of PDs vs. YDs, we observed both common and unique microbial signatures of response (Fig. 5). Among the features shared between PDs and YDs, it is worth noting that 24 h of fermentation with the oilseed mix resulted in decreased proportions of Ruminococcaceae and various members of Lachnospiraceae, and increased amounts of Enterobacteriaceae (p value ≤ 0.1). Although in the absence of statistical significance, these variations were also observed in the control vessels (i.e. with the addition of inulin or without any supplement), suggesting that they could be related to the experiment per se. Nonetheless, the Enterobacteriaceae increase in the 24-h faecal-derived microbial communities of YDs in the presence of the oilseed mix was far greater than in PDs (p value = 0.012) (24-h relative abundance: YDs with the oilseed mix, 78.5%; PDs with the oilseed mix, 27.7%; YDs with inulin, 46.4%; YDs without addition, 48.2%; PDs with inulin, 21.5%; PDs without addition, 17.5%). Going down in the taxonomic classification, we found that this increase was largely attributable to Klebsiella, whose 24-h relative abundance was 51.3 vs. 1.7% in the YD and PD group, respectively. The mean values in the YD and PD controls were 6.2% (± 8.7%, SD) and 1.7% (± 1.3%), respectively. Unfortunately, the other enterobacteria were unclassified at the genus level. At the species level, Klebsiella ASVs were found to be variously assigned to K. pneumoniae and K. aerogenes, with the latter being overall more prevalent (i.e. present in 9/27 fermentation samples of YDs, of which 5 in the presence of oilseed mix, while only in 5/27 for PDs, of which 3 in the presence of oilseed mix). Furthermore, it is interesting to point out that the family Clostridiaceae increased in PD-derived microbial communities in the presence of oilseed mix but not in the YD-related ones, as well as the genus Collinsella, whose trend was exactly opposite (i.e. decreased) in YD-related samples (p value ≤ 0.2).

Fig. 5
figure 5

Families and genera of the faecal-derived microbial communities of premenopausal (PD) and younger healthy women (YD), differing over time in 24-h fermentation experiments in the presence of oilseed mix, inulin or without additions (Ctr -). Boxplots showing the relative abundance distribution of bacterial taxa in the different study groups at T0, T5 and T24 (YD inulin, light red; YD Ctr -, light yellow; YD oilseed mix, light green; PD inulin, dark red; PD Ctr -, dark yellow; PD oilseed mix, dark green). Only significantly different taxa or trends are shown (p value ≤ 0.2, Wilcoxon test)


In this work, we determined the concentration of four main dietary lignans in an oilseed mix commercially prepared from the following ingredients: pumpkin, sunflower seeds, buckwheat and millet flakes, milled flaxseed, hemp and chia. Compared to flaxseed, the oilseed mix had 4 times less free lignans with Seco being 8.2 times less and Mat almost 6 times less than in flaxseed. Comparable volumes of Lari were found in both flaxseed and oilseed mix: 44.2 mg/100 g and 44.9 mg/100 g, respectively. Seco and Mat are the main substrate for enterolignans, with Seco being converted to EL via ED, and Mat being directly converted to EL. Differences in lignan concentration have been previously discussed [5, 48] and are thought to be related to variance in crop, climate, storage conditions and other factors. Therefore, it was expected to obtain concentrations that would diverge from previously reported data.

We then used the oilseed mix in batch culture fermentation experiments to evaluate the dynamics of enterolignan production by faecal-derived microbial communities from younger healthy (YD) vs. premenopausal (PD) female donors [49]. In vitro fermentation systems have been used for many years, however most models have yet to fully characterize the reproducibility (technical replicates), repeatability (biological replicates), and stability (time to and maintenance of enviromental conditions) of cultured communities [39,40,41,42,43,44,45].

Adding 0.5 g of oilseed mix to the stool samples of YDs resulted in a tenfold increase in EL concentration and a decrease in ED in 24 h, as compared to the negative control. Contrary to this, the same amount of supplement in the presence of PD samples resulted in a significant ED boost, whose 24-h concentration increased 110 times compared to the negative control and 80 times compared to the YD group under the same experimental conditions. However, despite such a high volume of ED, EL levels continued to decline, probably indicating the failure of bacterial communities from PDs to efficiently convert ED to EL. Consistent with this, i.e. with the inability of PD-related microbiotas to produce EL, high levels of Mat were found at the end of fermentation [50, 51].

The biotransformation of dietary lignans into EL and ED is a complex process believed to involve microbial consortia capable of performing four major catalytic reactions—O-deglycosylation, O-demethylation, dihydroxylation, and dehydrogenation—with each reaction being catalyzed by a specialized group of bacteria that often constitute a minor component of the overall microbiome [34,35,36]. Although evidence has been provided of the association between specific enterolignans and certain human gut bacteria, such as that between high urinary EL excretion and increased proportions of Ruminococcus spp. while decreased amounts of Bilophila [36, 52], the complexity and diversity of the gut microbiota are considered essential to maximise lignan conversion and thus influence human exposure to enterolignans. In this regard, here we demonstrated that YDs and PDs have a different intestinal microbiota profile and that these profiles undergo different temporal dynamics in the presence of oilseed mix, potentially contributing to the different production pattern of EL/ED, as observed in our in vitro batch culture fermentations. In particular, the early changes observed in the dominant structure of the faecal-derived microbial ecosystem of PDs seem to suggest a less resilient configuration and therefore more susceptible to disturbances, compared to YDs. From a taxonomic point of view, PDs were found to be significantly depleted in Collinsella, while being enriched in a number of Bacteroidetes members, including Bacteroides and Rikenellaceae, that typically thrive on high-fat diets and are enriched in obese gut microbiomes [53]. As regards the impact of the oilseed mix, both YD and PD faecal-derived microbial communities showed an enrichment over time in Enterobacteriaceae, specifically Klebsiella, but with a 24-h abundance much greater in the former. As previously observed [54, 55], several Klebsiella spp.—closely related to K. pneumoniae—are capable of biotransforming SDG into Seco, with the latter representing a substrate for EL production through microbial metabolism. The differential representation of Klebsiella could therefore help explain the greater levels of EL in YD- vs. PD-related samples. On the other hand, it is worth noting the differential trend of Collinsella, whose relative abundance increased in PD-related microbial communities while it decreased in YD-related ones in the presence of oilseed mix, which suggests a possible role of this taxon in EL/ED conversion. Interestingly, Collinsella has recently been described as a bacterial genus involved in the biotransformation of dietary phytoestrogens in equol, possibly contributing to the desirable health benefits of isoflavones in postmenopausal women [56, 57]. Finally, the PD faecal-derived microbial ecosystem appeared to be enriched in the Clostridiaceae family, whose members are already known to be variously involved in lignan metabolism [58, 59].


Despite the small sample size and obvious limitations of an in vitro study, our findings support that the faecal microbiota of healthy younger women is more efficient in converting dietary lignans into enterolignans, especially EL, than the microbiota of older premenopausal women. Furthermore, our study points to Klebsiella and Collinsella spp. as previously overlooked actors of the human gut microbiota with a strong potential to be directly involved in EL/ED production. To this end, in vitro gut microbiota models have proved to be useful tools to study the effects of food components, diet, and pharmaceutical molecules on gut microbiota composition. The advantages and limitations of in vivo and in vitro models as well as the developments to improve the modeling of host-microbe interactions have been extensively reviewed. Future in vivo studies in larger cohorts are needed to confirm our findings and understand how we can modulate the microbiota composition, favouring microbial species capable of producing EL more effectively. If a woman’s microbiota effectively converts dietary lignans to EL, consistent consumption of lignan-rich foods will potentially preserve a woman’s health in the long term, preventing or at least delaying the onset of degenerative conditions typically associated with menopause.


Oil seed samples

Eurocan Ltd (UK) provided five different organic raw oil seeds (flaxseed, chia, hulled sunflower, hulled pumpkin, hulled hemp), organic buckwheat and millet flakes, as well as a ready-made oilseed mix “MightyMix”. Oilseed mix was milled using Hinari Genie MB280 electric grinder. Duplicates 25 mg of the mix were stored in separate 2-ml Eppendorf tubes at − 20 °C. The present oilseed mix nutritional profile was characterised by Campden BRI (Chipping Campden) Ltd, accredited to ISO17025:2005 by UKAS by official reference methods (fat: Weibull Stoldt; sugars: HPLC; total dietary fibre: AOAC; protein: Kehjdahl) as reported in Additional file 5: Table S2.

Stool sample collection and preparation

Faecal samples were donated by three healthy younger (aged 25–30 years) and three older premenopausal female donors (aged 40–55 years). All donors confirmed to be healthy of metabolic and gastrointestinal conditions, were not taking prebiotic or probiotic supplements and did not receive antibiotic treatments in the 6 months preceding the study. Information on the donors’ health status, lifestyle habits, clinical anamnesis, and medicine use was collected with a pre-informative questionnaire. All faecal samples were collected on site, kept at − 20 °C and used within a maximum of 15 min after collection. The samples were diluted 1/10 w/v in anaerobic PBS (0.1 mol/l phosphate buffer solution, pH 7.4) and homogenized (Stomacher 400, Seward, West Sussex, UK) for 2 min at 460 paddle-beats.

In vitro batch culture fermentation experiments

The batch culture fermentation method was carried out as previously described by Gibson et al. [39]. Each vessel was inoculated with 5 ml of fresh faecal slurry (1/10 w/v) for both healthy and premenopausal subjects. A known prebiotic compound inulin (Raftilose P95, 95% oligosaccharide, β(2-1)-fructan, of which 60% w/w glucose-fructose, 40% w/w fructose, degree of polymerization, 3–10) serving as a positive control was added to a separate batch-culture vessel. A further vessel was prepared under the same conditions but without the addition of any compound (negative control, “ctr-”) whereas another vessel was used to add the oilseed mix. Batch culture fermenters were ran under anaerobic conditions for a period of 24 h during which samples (5 ml) were collected at time 0, 0.2, 5 and 24 h. Specifically, T0 samples for microbiota analysis were taken from the respective vessels under operating conditions and not from the women’s feces. Samples were stored at − 80 °C until needed for 16S rRNA gene-based next-generation sequencing analysis and ultra‐high-performance liquid chromatography (UHPLC)–tandem mass spectrometry (MS/MS) quantification.

Lignan extraction methods

The methodology for the extraction of lignans from single food samples and enterolignans from faecal samples was adapted from the work of Nørskov and Knudsen [46] as well as Milder et al. [47] and optimised in regards to the weight and the character of the samples. To each oilseed sample (25 mg), 1 ml of n-hexane, was added. The samples were vortexed and left at room temperature with gentle agitation for 20 min. Samples were centrifuged at 13,200 rpm, at 4 °C for 10 min. The supernatant was discarded, and the pellets were kept for the next steps.

Extraction of free lignans from oilseed samples and fermentation samples

Oilseed free lignan samples

Pellets were extracted with 0.5 ml of 100% methanol, vortexed and left at room temperature with gentle agitation for 10 min. After centrifugation at 13,200 rpm, at 4 °C for 10 min, the supernatant was collected into clean 2-ml Eppendorf tubes and left at nitrogen stream to evaporate for dryness.

Fermentation enterolignan samples

One milliliter of fermentation sample was extracted with 0.5 ml of 100% methanol, sonicated for 10 min, and then kept with gentle agitation for 10 min. Next, samples were centrifuged at 13,200 rpm, for 10 min at 4 °C, and the supernatant was collected into clean 2-ml tubes and evaporated for dryness.

Afterwards, dried fermentation and oilseed samples were incubated for 16 h at 37 °C with added ß-glucuronidase/sulfatase from Helix pomatia (freshly dissolved in 0.05 M sodium acetate buffer with concentration of 2 mg/ml), cooled and added with 0.5 ml of acidified water (0.4% of formic acid) to stop the hydrolysis, vortexed and then centrifuged at 13,200 rpm for 10 min at 4 °C. Samples were ready for solid phase extraction (SPE).

Extraction of total lignans from oilseed samples

Following defatting, pellets were extracted with 0.5 ml of 0.3 M NaOH in 70% methanol, vortexed, and then incubated for 1 h with gentle agitation at 60 °C. After cooling down, the pH was adjusted by adding 20 µl of glacial acetic acid (pH = 5). Next, samples were centrifuged for 10 min at 13,200 rpm at 4 °C; the supernatant was collected into 2-ml clean Eppendorf tubes and evaporated for dryness under a nitrogen stream. Enzymatic hydrolysis was performed by adding 0.6 ml of ß-glucuronidase/sulfatase from Helix pomatia (same as above) to each dried sample and setting samples for overnight incubation at 37 °C coupled with gentle agitation. Afterwards, samples were cooled down, added with 0.5 ml of acidified water (0.4% of formic acid), vortexed and then centrifuged at 13,200 rpm, at 4 °C for 10 min. Supernatants were collected into clean Eppendorf tubes ready SPE.

Solid Phase Extraction (SPE) of lignans and enterolignans

SPE was performed using 1-ml cartridges Strata® C-18 (55 µm, 70Å) from Phenomenex UK and SPE Vacuum Manifold. Waste was collected into 15-ml Falcon tubes, and final samples were collected into clean 2-ml Eppendorf tubes. After assembling manifold, cartridges were prepared as follows: with locked taps, 0.5 ml of acetonitrile was added to each cartridge and left for 10 min, then drained out. Next, with locked taps, 0.5 ml of LC–MS water was added, then drained after 10 min. After this, samples were added for slow elution through C18 material, followed by 0.5 ml of methanol, added twice to each cartridge and let to elute slowly until dry, afterwards, the vacuum was applied to dry the sorbent. Each cartridge was eluted with 0.4 ml of acetonitrile, and after draining, the vacuum was applied to facilitate full elution. Samples were then evaporated under a nitrogen stream and stored at − 20 °C. Aliquots of 0.5 ml of LC–MS water containing the internal standard (seco-d6) at the final concentration of 20 ng/ml and fermentation samples were added with 100 µl of LC–MS water. Samples were vortexed and 200 µl was dispersed into well plates ready for LC–MS analysis.

Microbial DNA extraction

Total microbial DNA was extracted from around 250 mg of in vitro fermentation samples using the QIAamp DNA Stool Mini Kit (QIAGEN, Hilden, Germany) according to the manufacturer’s instructions. DNA concentration and quality were evaluated using the NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Wilmington, DE).

16S rRNA gene-based next-generation sequencing, bioinformatics and statistics

The V3–V4 hypervariable region of the 16S rRNA gene was PCR-amplified using the primer set 341F/805R, as previously reported [60]. PCR products of about 460 bp were purified using a magnetic bead-based system (Agencourt AMPure XP; Beckman Coulter, Brea, CA) and indexed by limited-cycle PCR using Nextera technology (Illumina, San Diego, CA). Indexed libraries, further cleaned up as described above, were pooled at equimolar concentration, denatured and diluted to 6 pmol/l. Sequencing was performed on an Illumina MiSeq platform using the 2 × 250 bp paired-end protocol, according to the manufacturer’s instructions. Sequencing reads were deposited in the National Center for Biotechnology Information Sequence Read Archive (NCBI SRA; BioProject ID PRJNA592433).

The obtained paired-end reads were processed using a pipeline combining PANDAseq [61] and QIIME 2 [62]. High-quality reads were filtered and clustered into amplicon sequence variants (ASVs) at 99% similarity through an open-reference strategy performed with DADA2 [63]. Taxonomy was assigned using the vsearch classifier [64] against Greengenes database as a reference (release May 2013). Alpha diversity was measured using the number of observed ASVs and the Faith’s Phylogenetic Diversity (PD whole tree). Beta diversity was computed based on weighted and unweighted UniFrac distances and visualized on a Principal Coordinates Analysis (PCoA) plot. For the identification of Klebsiella species, the ASVs assigned to the genus Klebsiella were subjected to BLAST analysis [65]. Statistics was performed using R Studio 1.0.44 on R software version 3.3.2 [66] implemented with the packages stats and vegan [67]. The significance of data separation in the PCoA plot was tested by a permutation test with ANOSIM and pseudo-F ratio statistics using the function adonis in vegan. Bar plots were built using the R packages made4 [68] and vegan. Non-parametric tests (Kruskal–Wallis test or Wilcoxon test, paired or unpaired as needed) were achieved using the stats package. A p value ≤ 0.05 was considered statistically significant; a p value between 0.05 and 0.2 was considered a tendency.

Ultra‐high-performance liquid chromatography (UHPLC)–MS/MS

Acquity H class UPL chromatography equipment was used and separations were performed on an Aquity UPLC® HSS PFP 1.8 µm 2.1 × 100 mm C18 (Waters, UK) class column with a column protection of Acquity UPLC® HSS T3 1.8 µm Van Guard™ pre-column 3/Pk 2.1 × 5 mm column (Waters, UK) at a flow rate of 0.65 ml/min at 30 °C. The mobile phases A and B consisted of 100% LC–MS water and 100% acetonitrile, respectively. The gradient started at 95% phase A and 5% phase B, was held constant for 6 min, then phase B increased to 75% during 0.9 min, with the subsequent increase to 95% for 0.1 min, followed by 99% increase of phase A during the last 2 min. The total run for each sample was 10 min. Sample injection volume was 2 µl. The negative mode was used for ionisation. Detection was performed using Xevo TQ-micro (Waters, UK) quadrupole mass spectrometer, which facilitates the detection of low concentrated analytes. Parent and daughter ions (m/z) are described for each compound in Table 1 together with cone voltage and collision energy. The data were collected and analysed using MassLynx software.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request. Sequencing reads were deposited in the National Center for Biotechnology Information Sequence Read Archive (NCBI SRA; BioProject ID PRJNA592433).



Secoisolariciresinol diglucoside














Oestrogen receptor


Solid phase extraction


Liquid chromatography–mass spectrometry


Multiple reaction monitoring


Ultra-high-performance liquid chromatography


Younger donor


Premenopausal donor


  1. Peñalvo JL, Haajanen KM, Botting N, Adlercreutz H. Quantification of lignans in food using isotope dilution gas chromatography/mass spectrometry. J Agric Food Chem. 2005;53(24):9342–7.

    Article  PubMed  CAS  Google Scholar 

  2. Landete JM. Updated knowledge about polyphenols: functions, bioavailability, metabolism, and health. Crit Rev Food Sci Nutr. 2012;52(10):936–48.

    Article  CAS  PubMed  Google Scholar 

  3. Phenol-Explorer. Accessed on 24 Mar 2020.

  4. Buck K, Zaineddin AK, Vrieling A, Heinz J, Linseisen J, Flesch-Janys D, Chang-Claude J. Estimated enterolignans, lignan-rich foods, and fibre in relation to survival after postmenopausal breast cancer. Br J Cancer. 2011;105(8):1151.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Durazzo A, Lucarini M, Camilli E, Marconi S, Gabrielli P, Lisciani S, Gambelli L, Aguzzi A, Novellino E, Santini A. Dietary lignans: definition, description and research trends in databases development. Molecules. 2018;23(12):3251.

    Article  PubMed Central  CAS  Google Scholar 

  6. Aehle E, Müller U, Eklund PC, Willför SM, Sippl W, Dräger B. Lignans as food constituents with estrogen and antiestrogen activity. Phytochemistry. 2011;72(18):2396–405.

    Article  CAS  PubMed  Google Scholar 

  7. Mali AV, Padhye SB, Anant S, Hegde MV, Kadam SS. Anticancer and antimetastatic potential of enterolactone: clinical, preclinical and mechanistic perspectives. Eur J Pharmacol. 2019;852:107–24.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Mueller SO, Simon S, Chae K, Metzler M, Korach KS. Phytoestrogens and their human metabolites show distinct agonistic and antagonistic properties on estrogen receptor α (ERα) and ERβ in human cells. Toxicol Sci. 2004;80(1):14–25.

    Article  CAS  PubMed  Google Scholar 

  9. Habauzit D, Armengaud J, Roig B, Chopineau J. Determination of estrogen presence in water by SPR using estrogen receptor dimerization. Anal Bioanal Chem. 2008;390(3):873–83.

    Article  CAS  PubMed  Google Scholar 

  10. Adlercreutz H, Hämäläinen E, Gorbach S, Goldin B. Dietary phyto-oestrogens and the menopause in Japan. The Lancet. 1992;339(8803):1233.

    Article  CAS  Google Scholar 

  11. Kajla P, Sharma A, Sood DR. Flaxseed—a potential functional food source. J Food Sci Technol. 2015;52(4):1857–71.

    Article  CAS  PubMed  Google Scholar 

  12. Zand RSR, Jenkins DJ, Diamandis EP. Flavonoids and steroid hormone-dependent cancers. J Chromatogr B. 2002;777(1–2):219–32.

    Article  CAS  Google Scholar 

  13. Allard C, Bonnet F, Xu B, Coons L, Albarado D, Hill C, Fagherazzi G, Korach KS, Levin ER, Lefante J. Activation of hepatic estrogen receptor-α increases energy expenditure by stimulating the production of fibroblast growth factor 21 in female mice. Mol Metab. 2019;22:62–70.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Rüegg J, Cai W, Karimi M, Kiss NB, Swedenborg E, Larsson C, Ekström TJ, Pongratz I. Epigenetic regulation of glucose transporter 4 by estrogen receptor β. Mol Endocrinol. 2011;25(12):2017–28.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  15. Sharma G, Mauvais-Jarvis F, Prossnitz ER. Roles of G protein-coupled estrogen receptor GPER in metabolic regulation. J Steroid Biochem Mol Biol. 2018;176:31–7.

    Article  CAS  PubMed  Google Scholar 

  16. Dosi R, Bhatt N, Shah P, Patell R. Cardiovascular disease and menopause. J Clin Diagn Res. 2014;8(2):62.

    PubMed  PubMed Central  Google Scholar 

  17. Speth RC, D’Ambra M, Ji H, Sandberg K. A heartfelt message, estrogen replacement therapy: use it or lose it. Am J Physiol Heart Circul Physiol. 2018;315(6):H1765–78.

    Article  CAS  Google Scholar 

  18. Rosano GM, Panina G. Oestrogens and the heart. Therapie. 1999;54(3):381–5.

    CAS  PubMed  Google Scholar 

  19. Bloedon LT, Balikai S, Chittams J, Cunnane SC, Berlin JA, Rader DJ, Szapary PO. Flaxseed and cardiovascular risk factors: results from a double blind, randomized, controlled clinical trial. J Am Coll Nutr. 2008;27(1):65–74.

    Article  CAS  PubMed  Google Scholar 

  20. Kwa M, Plottel CS, Blaser MJ, Adams S. The intestinal microbiome and estrogen receptor–positive female breast cancer. JNCI J Natl Cancer Instit. 2016;108(8):djw029.

    Google Scholar 

  21. Gass ML, Khan S. Estrogen and estrogen analogs for prevention and treatment of osteoporosis. In: Osteoporosis. Edited by Anonymous Elsevier; 2013, p. 1805–25.

  22. Krebs NF, Gao D, Gralla J, Collins JS, Johnson SL. Efficacy and safety of a high protein, low carbohydrate diet for weight loss in severely obese adolescents. J Pediatr. 2010;157(2):252–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Adlercreutz H, Bannwart C, Wähälä K, Mäkelä T, Brunow G, Hase T, Arosemena PJ, Kellis JT Jr, Vickery LE. Inhibition of human aromatase by mammalian lignans and isoflavonoid phytoestrogens. J Steroid Biochem Mol Biol. 1993;44(2):147–53.

    Article  CAS  PubMed  Google Scholar 

  24. Piller R, Chang-Claude J, Linseisen J. Plasma enterolactone and genistein and the risk of premenopausal breast cancer. Eur J Cancer Prev. 2006;15(3):225–32.

    Article  CAS  PubMed  Google Scholar 

  25. Aarestrup J, Kyrø C, Knudsen KE, Weiderpass E, Christensen J, Kristensen M, Würtz AM, Johnsen NF, Overvad K, Tjønneland A. Plasma enterolactone and incidence of endometrial cancer in a case–cohort study of Danish women. Br J Nutr. 2013;109(12):2269–75.

    Article  CAS  PubMed  Google Scholar 

  26. Peterson J, Dwyer J, Adlercreutz H, Scalbert A, Jacques P, McCullough ML. Dietary lignans: physiology and potential for cardiovascular disease risk reduction. Nutr Rev. 2010;68(10):571–603.

    Article  PubMed  Google Scholar 

  27. Lampe JW. Isoflavonoid and lignan phytoestrogens as dietary biomarkers. J Nutr. 2003;133(3):956S–64S.

    Article  CAS  PubMed  Google Scholar 

  28. Gaya P, Peirotén Á, Medina M, Landete JM. Bifidobacterium adolescentis INIA P784: the first probiotic bacterium capable of producing enterodiol from lignan extracts. J Funct Foods. 2017;29:269–74.

    Article  CAS  Google Scholar 

  29. Clavel T, Henderson G, Alpert C, Philippe C, Rigottier-Gois L, Doré J, Blaut M. Intestinal bacterial communities that produce active estrogen-like compounds enterodiol and enterolactone in humans. Appl Environ Microbiol. 2005;71(10):6077–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Gaya P, Sánchez-Jiménez A, Peirotén Á, Medina M, Landete JM. Incomplete metabolism of phytoestrogens by gut microbiota from children under the age of three. Int J Food Sci Nutr. 2018;69(3):334–43.

    Article  CAS  PubMed  Google Scholar 

  31. Poluzzi E, Piccinni C, Raschi E, Rampa A, Recanatini M, De Ponti F. Phytoestrogens in postmenopause: the state of the art from a chemical, pharmacological and regulatory perspective. Curr Med Chem. 2014;21(4):417–36.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Chen MN, Lin CC, Liu CF. Efficacy of phytoestrogens for menopausal symptoms: a meta-analysis and systematic review. Climacteric. 2015;18(2):260–9.

    Article  CAS  PubMed  Google Scholar 

  33. Edel AL, Patenaude AF, Richard MN, Dibrov E, Austria JA, Aukema HM, Pierce GN, Aliani M. The effect of flaxseed dose on circulating concentrations of alpha-linolenic acid and secoisolariciresinol diglucoside derived enterolignans in young, healthy adults. Eur J Nutr. 2016;55(2):651–63.

    Article  CAS  PubMed  Google Scholar 

  34. Clavel T, Henderson G, Engst W, Doré J, Blaut M. Phylogeny of human intestinal bacteria that activate the dietary lignan secoisolariciresinol diglucoside. FEMS Microbiol Ecol. 2006;55(3):471–8.

    Article  CAS  PubMed  Google Scholar 

  35. Hullar MA, Lancaster SM, Li F, Tseng E, Beer K, Atkinson C, Wähälä K, Copeland WK, Randolph TW, Newton KM. Enterolignan-producing phenotypes are associated with increased gut microbial diversity and altered composition in premenopausal women in the United States. Cancer Epidemiol Prevent Biomark. 2015;24(3):546–54.

    Article  CAS  Google Scholar 

  36. Lagkouvardos I, Kläring K, Heinzmann SS, Platz S, Scholz B, Engel K, Schmitt-Kopplin P, Haller D, Rohn S, Skurk T. 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.

    Article  CAS  PubMed  Google Scholar 

  37. Holma R, Kekkonen RA, Hatakka K, Poussa T, Vapaatalo H, Adlercreutz H, Korpela R. Low serum enterolactone concentration is associated with low colonic LactobacillusEnterococcus counts in men but is not affected by a synbiotic mixture in a randomised, placebo-controlled, double-blind, cross-over intervention study. Br J Nutr. 2014;111(2):301–9.

    Article  CAS  PubMed  Google Scholar 

  38. Samuel BS, Hansen EE, Manchester JK, Coutinho PM, Henrissat B, Fulton R, Latreille P, Kim K, Wilson RK, Gordon JI. Genomic and metabolic adaptations of Methanobrevibacter smithii to the human gut. Proc Natl Acad Sci. 2007;104(25):10643–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Gibson GR, Fuller R. Aspects of in vitro and in vivo research approaches directed toward identifying probiotics and prebiotics for human use. J Nutr. 2000;130(2):391S–5S.

    Article  CAS  PubMed  Google Scholar 

  40. Grimaldi R, Swann JR, Vulevic J, Gibson GR, Costabile A. Fermentation properties and potential prebiotic activity of Bimuno® galacto-oligosaccharide (65% galacto-oligosaccharide content) on in vitro gut microbiota parameters. Br J Nutr. 2016;116(3):480–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Cinquin C, Le Blay G, Fliss I, Lacroix C. Immobilization of infant fecal microbiota and utilization in an in vitro colonic fermentation model. Microb Ecol. 2004;48(1):128–38.

    Article  CAS  PubMed  Google Scholar 

  42. Van den Abbeele P, Grootaert C, Marzorati M, Possemiers S, Verstraete W, Gérard P, Rabot S, Bruneau A, El Aidy S, Derrien M. Microbial community development in a dynamic gut model is reproducible, colon region specific, and selective for Bacteroidetes and Clostridium cluster IX. Appl Environ Microbiol. 2010;76(15):5237–46.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  43. Pompei A, Cordisco L, Raimondi S, Amaretti A, Pagnoni UM, Matteuzzi D, Rossi M. In vitro comparison of the prebiotic effects of two inulin-type fructans. Anaerobe. 2008;14(5):280–6.

    Article  PubMed  CAS  Google Scholar 

  44. Lesmes U, Beards EJ, Gibson GR, Tuohy KM, Shimoni E. Effects of resistant starch type III polymorphs on human colon microbiota and short chain fatty acids in human gut models. J Agric Food Chem. 2008;56(13):5415–21.

    Article  CAS  PubMed  Google Scholar 

  45. Macfarlane GT, Gibson GR, Cummings JH. Comparison of fermentation reactions in different regions of the human colon. J Appl Bacteriol. 1992;72(1):57–64.

    CAS  PubMed  Google Scholar 

  46. Nørskov NP, Knudsen KEB. Validated LC–MS/MS method for the quantification of free and bound lignans in cereal-based diets and feces. J Agric Food Chem. 2016;64(44):8343–51.

    Article  PubMed  CAS  Google Scholar 

  47. Milder IE, Arts IC, Venema DP, Lasaroms JJ, Wähälä K, Hollman PC. Optimization of a liquid chromatography–tandem mass spectrometry method for quantification of the plant lignans secoisolariciresinol, matairesinol, lariciresinol, and pinoresinol in foods. J Agric Food Chem. 2004;52(15):4643–51.

    Article  CAS  PubMed  Google Scholar 

  48. Smeds AI, Eklund PC, Willför SM. Content, composition, and stereochemical characterisation of lignans in berries and seeds. Food Chem. 2012;134(4):1991–8.

    Article  CAS  PubMed  Google Scholar 

  49. Ramezani TF, Bidhendi YR, Solaymani-Dodaran M, Tohidi M, Firouzi F, Azizi F. Improving prediction of age at menopause using multiple Anti Mullerian Hormone measurements; Tehran Lipid-Glucose Study. J Clin Endocrinol Metab. 2020;105(2):1–27.

    Google Scholar 

  50. Aura A, Martin-Lopez P, O’Leary KA, Williamson G, Oksman-Caldentey K, Poutanen K, Santos-Buelga C. In vitro metabolism of anthocyanins by human gut microflora. Eur J Nutr. 2005;44(3):133–42.

    Article  CAS  PubMed  Google Scholar 

  51. Mena P, Ludwig IA, Tomatis VB, Acharjee A, Calani L, Rosi A, Brighenti F, Ray S, Griffin JL, Bluck LJ. Inter-individual variability in the production of flavan-3-ol colonic metabolites: preliminary elucidation of urinary metabotypes. Eur J Nutr. 2019;58(4):1529–43.

    Article  CAS  PubMed  Google Scholar 

  52. Lampe JW, Kim E, Levy L, Davidson LA, Goldsby JS, Miles FL, Navarro SL, Randolph TW, Zhao N, Ivanov I. Colonic mucosal and exfoliome transcriptomic profiling and fecal microbiome response to a flaxseed lignan extract intervention in humans. Am J Clin Nutr. 2019;110(2):337–90.

    Article  Google Scholar 

  53. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, Ling AV, Devlin AS, Varma Y, Fischbach MA. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559.

    Article  CAS  PubMed  Google Scholar 

  54. Wang C, Ma X, Yang D, Guo Z, Liu G, Zhao G, Tang J, Zhang Y, Ma M, Cai S. Production of enterodiol from defatted flaxseeds through biotransformation by human intestinal bacteria. BMC Microbiol. 2010;10(1):115.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  55. Zhou Y, Zhu S, Yang D, Zhao D, Li J, Liu S. Characterization of Klebsiella sp. strain S1: a bacterial producer of secoisolariciresinol through biotransformation. Can J Microbiol. 2016;63(1):1–10.

    Article  PubMed  CAS  Google Scholar 

  56. Nakatsu CH, Armstrong A, Clavijo AP, Martin BR, Barnes S, Weaver CM. Fecal bacterial community changes associated with isoflavone metabolites in postmenopausal women after soy bar consumption. PLoS ONE. 2014;9(10):e108924.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  57. Guadamuro L, Dohrmann AB, Tebbe CC, Mayo B, Delgado S. Bacterial communities and metabolic activity of faecal cultures from equol producer and non-producer menopausal women under treatment with soy isoflavones. BMC Microbiol. 2017;17(1):93.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  58. Jin J, 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.

    Article  CAS  PubMed  Google Scholar 

  59. Yoder SC, Lancaster SM, Hullar MA, Lampe JW. Gut microbial metabolism of plant lignans: influence on human health. In: Diet-microbe interactions in the gut. Edited by Anonymous Elsevier; 2015, p. 103–17.

  60. Barone M, Turroni S, Rampelli S, Soverini M, D’Amico F, Biagi E, Brigidi P, Troiani E, Candela M. Gut microbiome response to a modern Paleolithic diet in a Western lifestyle context. PLoS One. 2019;14(8):1–14.

    Article  CAS  Google Scholar 

  61. Masella AP, Bartram AK, Truszkowski JM, Brown DG, Neufeld JD. PANDAseq: paired-end assembler for illumina sequences. BMC Bioinformatics. 2012;13(1):31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  62. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  63. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Rognes T, Flouri T, Nichols B, Quince C, Mahé F. VSEARCH: a versatile open source tool for metagenomics. PeerJ. 2016;4:e2584.

    Article  PubMed  PubMed Central  Google Scholar 

  65. Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J. Mol Biol (N Y). 1990;1990(215):403–10.

    Article  Google Scholar 

  66. The R Project for Statistical Computing. Accessed 24 Mar 2020.

  67. The Vegan Package. Accessed 24 Mar 2020.

  68. Culhane AC, Thioulouse J, Perrière G, Higgins DG. MADE4: an R package for multivariate analysis of gene expression data. Bioinformatics. 2005;21(11):2789–90.

    Article  CAS  PubMed  Google Scholar 

Download references


We thank the student Oghenerukevwe Obayiuwana for the technical assistance during the in vitro fermentation experiments and lignan extraction.


We thank to the University of Roehampton for funding this project.

Author information

Authors and Affiliations



GC and AC designed the experiments and wrote the manuscript. AK performed the in vitro fermentation experiments and UHPLC–MS/MS analysis. EK prepared the samples for DNA analysis. MB, ST and PB carried out 16S rRNA gene-based next-generation sequencing and analysis and wrote the microbiota-related parts. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Adele Costabile.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Fig. S1.

The biodiversity of the faecal-derived microbial communities from premenopausal and younger healthy women tends to decrease over time in 24-h fermentation experiments in the presence of oilseed mix and inulin. Boxplots showing the distribution of alpha diversity values, according to the number of observed amplicon sequence variants (ASVs, upper panel) and Faith’s Phylogenetic Diversity index (PD whole tree, lower panel), for the faecal microbial communities from premenopausal (PD) and younger healthy (YD) women at 0 (T0), 5 (T5) and 24 h (T24) of fermentation in the presence of oilseed mix (shades of green), inulin (as a positive control; shades of red) or without additions (negative control, “Ctr -”; shades of yellow).

Additional file 2: Fig. S2.

Weighted UniFrac-based Principal Coordinates Analysis (PCoA) of the faecal-derived microbial communities of premenopausal and younger healthy women over 24 h of fermentation in the presence of oilseed mix, inulin or without additions. Left, PCoA plot showing all fermentation samples, coloured by group of women (premenopausal, light blue vs. younger healthy, orange). A significant separation between groups was found, regardless of experimental condition (oilseed mix, inulin and negative control—“Ctr-”) and time point (T0, T5 and T24) (p value < 1 × 10−4, permutation test with pseudo-F ratio). Right, PCoA plots showing the fermentation samples for premenopausal women (top panel) and younger healthy women (bottom panel). Within each group of women, the samples separate significantly by experimental condition (oilseed mix, green; inulin, red; Ctr-, yellow) and time point (T0, circle; T5, square; T24, diamond) (p value < 0.001).

Additional file 3: Table S1.

Results of adonis and ANOSIM statistics applied to weighted and unweighted UniFrac distance-based ordination (related to Fig. 4 and Additional file 2: Fig. S2). PD: premenopausal donor; YD: younger donor.

Additional file 4: Fig. S3.

Genus-level relative abundance profiles of the faecal-derived microbial communities of premenopausal (PD) and younger healthy women (YD) at 0, 5 and 24 h of fermentation in the presence of oilseed mix, inulin or without additions. *, unclassified amplicon sequence variants reported at higher taxonomic level; **, other. For each group of women (YD and PD), the profiles are shown in the following order: samples in the presence of inulin at T0, T5 and T24 (red), samples in the negative control (Ctr -) at T0, T5 and T24 (yellow), and samples in the presence of oilseed mix at T0, T5 and T24 (green). The black arrow below the histograms indicates the temporal succession for each triplet of samples, i.e. T0, T5, T24.

Additional file 5: Table S2.

Composition and nutritional profile of the seed mix (MightyMix) used in the present study.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Corona, G., Kreimes, A., Barone, M. et al. Impact of lignans in oilseed mix on gut microbiome composition and enterolignan production in younger healthy and premenopausal women: an in vitro pilot study. Microb Cell Fact 19, 82 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: