- Research article
- Open Access
Bacterial communities and metabolic activity of faecal cultures from equol producer and non-producer menopausal women under treatment with soy isoflavones
BMC Microbiology volume 17, Article number: 93 (2017)
Isoflavones are polyphenols with estrogenic activity found mainly in soy and soy-derived products that need to be metabolised in the intestine by the gut bacteria to be fully active. There is little knowledge about isoflavone bioconversion and equol production in the human intestine. In this work, we developed an in vitro anaerobic culture model based on faecal slurries to assess the impact of isoflavone supplementation on the overall intestinal bacterial composition changes and associated metabolic transformations.
In the faecal anaerobic batch cultures of this study bioconversion of isoflavones into equol was possible, suggesting the presence of viable equol-producing bacterial taxa within the faeces of menopausal women with an equol producer phenotype. The application of high-throughput DNA sequencing of 16S rRNA gene amplicons revealed the composition of the faecal cultures to be modified by the addition of isoflavones, with enrichment of some bacterial gut members associated with the metabolism of phenolics and/or equol production, such as Collinsella, Faecalibacterium and members of the Clostridium clusters IV and XIVa. In addition, the concentration of short-chain fatty acids (SCFAs) detected in the isoflavone-containing faecal cultures was higher in those inoculated with faecal slurries from equol-producing women.
This study constitutes the first step in the development of a faecal culturing system with isoflavones that would further allow the selection and isolation of intestinal bacterial types able to metabolize these compounds and produce equol in vitro. Although limited by the low number of faecal cultures analysed and the inter-individual bacterial diversity, the in vitro results obtained in this work tend to indicate that soy isoflavones might provide an alternative energy source for the increase of equol-producing taxa and enhancement of SCFAs production. SCFAs and equol are both considered pivotal bacterial metabolites in the triggering of intestinal health-related beneficial effects.
Isoflavones are polyphenols with estrogenic activity found at relatively high concentration in soy and soy-derived products. Epidemiological studies suggest high intakes of isoflavones to be associated with fewer and less intense menopause symptoms, and a reduced incidence of hormone-dependent and aging-associated diseases such as osteoporosis, cardiovascular diseases and cancer . Although there is a growing body of scientific evidence showing beneficial effects in counteracting symptoms such as hot flushes and vasomotor reactions in menopausal women , the European Food Safety Authority (EFSA) has recently refuted health claims about the role of isoflavones in body functions . This may be so because the functionality of isoflavones depends on their bioavailability and their conversion into (more) active metabolites within the intestinal tract. In nature, isoflavones are mostly (>80%) conjugated with sugars as isoflavone-glycosides (daidzin, genistin, glycitin) with low availability and bioactivity . For full activity, aglycones (daidzein, genistein, glycitein) need to be released from their glycosides and, occasionally, metabolized to more potent metabolites, such as equol (produced from daidzein). Equol is the isoflavone-derived compound with the strongest estrogenic activity and antioxidant capacity, mechanisms by which isoflavones seem to provide the purported beneficial health effects . The conversion of isoflavones into equol appears to be mainly carried out by intestinal microorganisms, especially bacteria inhabiting the distal portion of the gut . Both, traditional and molecular methods have revealed marked individual diversity in intestinal bacterial communities. Such differences in the bacterial community composition may, in fact, influence the fate of isoflavone metabolic transformations [7, 8], which could ultimately contribute to differences in the physiological response to isoflavone treatment.
In that sense, only 30–50% of Western individuals produce equol , and only these might fully benefit from the positive health effects of isoflavone consumption. The metabolism of isoflavones involves several steps mediated by enzymes provided by different bacterial taxa [7, 8]. Even though scientific information is accumulating on the microorganisms producing equol, their enzymes and the metabolic routes involved, our current knowledge is still limited [6, 9, 10]. Most of the equol-producing bacteria characterized so far belong to the family Coriobacteriaceae . However, it is not yet clear whether this family is the only intestinal group acting on isoflavones and producing equol.
Previously, the metabolism of daidzein by faecal bacterial consortia has been microbiologically characterized by conventional culturing methods [11–13]. The availability of high-throughput DNA sequencing techniques opens new potentials to tracking changes in the bacterial communities during isoflavone supplementation in both in vitro and in vivo systems. A better knowledge about the identity and individual variability of the intestinal bacteria that metabolize isoflavones and convert them into equol, would provide an important step for developing strategies to increase bioavailability and concentration of active compounds, e.g. by supplying suitable equol-producing probiotic bacteria.
In this study, anaerobic batch cultures of faecal samples from equol producer and non-producer menopausal women (as determined by urine equol excretion of >1000 nM) under treatment with soy isoflavones were performed, using media with and without isoflavones, to examine which bacteria take benefit from soy isoflavones and could be involved in equol production. PCR- denaturing gradient gel electrophoresis (DGGE) and high-throughput DNA sequencing were used to determine whether faecal cultures grown in the presence of isoflavones showed any change in their bacterial community composition and/or structure. These molecular analyses, based on PCR-amplified partial 16S rRNA gene sequences, were complemented by metabolic profiling of the culture supernatants using ultra-high-performance liquid chromatography (UHPLC) and gas chromatography (GC). This approach allowed the suggestion of links between structural responses in bacterial community compositions to metabolic activities.
Stool samples from isoflavone-treated menopausal women
Menopausal women who had been receiving treatment for 6 months with 80 mg/day of an isoflavone concentrate (Fisiogen; Zambon, Bresso, Italy) were recruited. Stool samples were obtained from four women whose faecal microbiota had been characterized and whose equol producer status during isoflavone treatment had been determined in a previous study . Three of the women (WC, WG and WP) were equol producers (urine equol >1000 nM as defined by Rowland et al. ), while the fourth (WE) was a non-producer (<10 nM in urine). Faeces were collected and transported to the laboratory as previously described .
Faecal batch cultures
Ten-fold faecal dilutions were prepared by homogenizing 1 g of faeces in 9 ml of a pre-reduced phosphate buffer saline solution (PBS) under strict anaerobic conditions (80% N2, 10% CO2, 10% H2) in a Whitley DG500 Workstation anaerobic chamber. A 10% (v/v) aliquot of the resulting faecal slurry was used to inoculate the medium for colonic bacteria (MCB) described by van der Meulen et al. , modified to have a lower glucose content (2 g/l) (hereafter referred to as mMCB). This medium was prepared without and with a concentration of 160 mg/l commercial isoflavone supplement (Fisiogen) by dissolving two capsules per liter of medium. One capsule of Fisiogen contains, as stated in the package, 80 mg of isoflavone of which 55–72% is genistin/genistein and 28–45% other soy isoflavones. It also includes as excipients E341, E460, E468, E551, E472, E464, E904, E171, and E124. The isoflavone-containing medium is hereafter referred to as mMCBISO. These faecal primary cultures thus prepared were incubated at 37 °C in fermentation flasks with mild stirring (250 rpm) for 24 h under anaerobic conditions. To find out whether production of equol in vitro was maintained in successive subculturings, primary cultures in mMCBISO were used as inocula (10% v/v) for a second round of cultivation in fresh mMCBISO. These new cultures were incubated as above and referred to as secondary cultures.
Positive controls for isoflavone conversion and equol production were obtained by inoculating (at 10%) mMCBISO with the strains Slackia isoflavoniconvertens DSM 22006 or Slackia equolifaciens DSM 24851. Both these strains were grown under anoxic conditions in Gifu Anaerobic Medium (GAM, Nissui Pharmaceutical, Tokyo, Japan) supplemented with 0.5% arginine (Merck, Darmstadt, Germany). After overnight culturing, cell of the two strains were washed twice in pre-reduced PBS before inoculating the media (mMCB and mMCBISO) or the faecal cultures and incubated as above.
Detection and quantification of equol and its isoflavone precursors
Equol and its isoflavone precursors [daidzin (daidzein-7-O-glucoside) and daidzein] were measured using a UHPLC procedure based on a method for equol determination in urine . Briefly, after 24 h incubation, faecal cultures were centrifuged at 800 rpm for 10 min, duplicate samples of 3 ml of each supernatant collected, and their isoflavone content extracted using Bond Elute-C18 solid-phase cartridges (Agilent Technologies, Santa Clara, CA, USA). After drying, organic extracts were dissolved in 100 μl of HPLC grade methanol, and 1 μl of each sample injected into the UHPLC apparatus (Waters, Palo Alto, CA, USA). Samples were analysed in duplicate. Equol was quantified using a fluorescence detector (excitation 280 nm, emission 310 nm), while daidzin and daidzein were identified with a photodiode array (PDA) detector by comparison with their retention times and spectral characteristics at 260 nm. Quantification was performed against calibration curves prepared using commercial standards (all from Sigma-Aldrich, St. Louis, MO, USA).
Short chain fatty acids determination
Short chain fatty acids (SCFAs) in the primary and secondary cultures were determined by GC. Culture supernatants were centrifuged for 10 min at 800 rpm, and then for 10 min further at 13,200 rpm, before filtering through 0.45 μm PTFE filters. The filtered samples were then mixed 10:1 (v/v) with 1 mg/ml of 2-ethyl butyric acid (Sigma-Aldrich) dissolved in methanol as an internal standard. A chromatographic system composed of 6890 N GC apparatus (Agilent Technologies) connected to a flame ionization detector was used to identify and quantify the SCFAs as described elsewhere . All samples were analysed in triplicate. The Shapiro–Wilk test was initially used to check for the normal distribution of the data. Since the data were not normally distributed, we used the non-parametric Wilcoxon signed-rank test to examine the differences in the content of SCFAs between the faecal cultures with (mMCBISO) and without isoflavones (mMCB). This statistic test is used to assess differences between related groups of samples and in our case the same faecal cultures were analysed under two different conditions (mMCB versus mMCBISO). All these calculations were performed using the SPSS software v.22.0.
Microbial community analysis
Total bacterial DNA was extracted from the faecal primary and secondary cultures using the phenol-chloroform-based protocol of Zoetendal et al.  with modifications as described elsewhere . The main modification consisted in the addition of an enzymatic lysis step before mechanical disruption in a FastPrep FP120 apparatus (Qbiogene, Carlsbad, CA, USA). After precipitation and drying, DNA samples were suspended in 100 μl of sterile molecular biology grade water (Sigma-Aldrich).
PCR-denaturing gradient gel electrophoresis (DGGE) analysis
Total DNA from the faecal cultures with and without isoflavones was used as a template for amplification of the variable region V3 of the 16S rRNA gene by PCR using the universal prokaryotic primers F357-GC and R518 . DGGE electrophoresis was performed as described elsewhere . DNA from S. equolifaciens and S. isoflavoniconvertens was used to provide DGGE markers. Genomic DNA from these strains was obtained from overnight cultures in GAM + 0.5% arginine using the GenElute Bacterial Genomic DNA kit (Sigma-Aldrich) and subsequently employed in PCR amplifications. Amplicons were purified using GenElute PCR Clean-Up columns (Sigma-Aldrich) and mixed in equal amounts to obtain a DGGE marker.
High-throughput DNA sequencing and analysis of 16S rRNA gene amplicons
DNA samples obtained from the primary and secondary cultures were subjected to PCR amplification of the variable region V4-V5 of the 16S rRNA gene using the universal prokaryotic primers 515F and 909R ; these were designed to include Illumina adapters. Double indexed amplicons were generated by the protocol of Caporaso et al.  with minor modifications. In brief, DNA from the samples was extracted in duplicate and amplified in triplicate; in total, 72 amplification reactions were obtained. These were subsequently paired-end sequenced in an Illumina Miseq System (Illumina, San Diego, CA, USA) by the StarSeq Company (Mainz, Germany), and treated as independent replicates.
Bioinformatic analysis was performed using the Mothur software package (v.1.34.1) following the MiSeq Standard Operating Procedure (SOP) . Briefly, sequences longer than 380 bp in length, or shorter than 370 bp, and those containing ambiguous base pairs or homopolymers regions of >8 bp were removed. All other sequences were aligned using the SINA alignment service of the SILVA 16S rRNA sequence database. Chimera removal was performed using the UCHIME algorithm . A random subset of 35,000 sequences per sample was used to balance numbers of reads among samples. Sequences were then clustered into operational taxonomic units (OTUs) using a 0.03 dissimilarity cut-off. Sequences were taxonomically classified using the Ribosomal Database Project (RDP) database. The Bayesian classifier with an 80% confidence threshold was used in the taxonomic assignment with the genus level as the lowest taxonomic unit considered. The MOTHUR program was also used to perform weighted UniFrac analysis, which was employed to assess the similarity of the microbial communities between samples. Construction of a heatmap was performed using the R statistical software. Clustering was accomplished using the complete linkage method with Euclidean distance measure. Multivariable statistical analysis was performed by principal coordinates analysis (PCoA) and non-metric multidimensional scaling (NMDS) with two dimensions. Differences in the microbial composition of faecal cultures with and without isoflavones were sought by analysis of molecular variance (AMOVA) and analysis of similarities (ANOSIM). The identification of differentially abundant taxa was assessed using the Metastats software . Multiple hypothesis tests were adjusted using the false discovery rate (FDR) correction ; an FDR q-value threshold of 0.25 was used to identify significant differences. In an attempt to assign the differential OTUs at the species level, manual sequence comparisons were performed against the Greengenes 16S rRNA gene database.
Isoflavone conversion and equol production
The uninoculated medium with isoflavones, mMCBISO, was shown by UHPLC to contain three times more conjugated precursor daidzin that its corresponding aglycone daidzein (Table 1). As expected, equol was never detected in the absence of isoflavones in the cultures. In the medium inoculated with S. equolifaciens DSM 24851, used as a positive control for equol production, daidzin levels were maintained while daidzein was transformed and converted to equol, reaching 2 μg/ml after 24 h incubation. In the primary faecal cultures from the three equol producers (WC, WG and WP) comparable amounts of equol to that produced by the positive control were recorded (1.40–3.07 μg/ml). Congruently, daidzin levels were drastically reduced in these cultures, indicating glycosidase activity from faecal bacteria towards this isoflavone glycoside. In the primary culture from the non-producing woman (WE), daidzin was also shown to decrease while daidzein accumulated, but no equol was found (limit of quantification 0.002 μg/ml). As expected, in the secondary culture from this woman equol was not detected either. In contrast, equol production was maintained to some extent in the secondary cultures from two of the three equol producers (WG and WP).
Comparing the results obtained in the mMCB medium, total SCFA production was greater in the faecal cultures containing isoflavones (mMCBISO) (Fig. 1). The concentrations of the major SCFAs (acetic, propionic and butyric acids) were particularly enhanced in the primary faecal cultures whose inocula came from the equol-producing women (Table 2); significant increases were seen for propionic acid and isobutyric and isovaleric acids (the latter two are minority, branched chain fatty acids [BCFAs]) (p < 0.05). In contrast, in the cultures from the equol non-producing woman (WE), no differences between the corresponding mMCBISO and mMCB were seen for major or minor SCFAs in the primary cultures. Consequently, the acetic/propionic ratio for this faecal culture was five times higher than that recorded for the equol producers; indeed, the concentration of propionic and butyric acids for this culture were much lower than those recorded in the primary cultures of the equol-producing women. SCFAs production in the secondary cultures from these women was maintained to some extent, but with a major reduction in the concentration of isovaleric acid.
Bacterial community analysis by PCR-DGGE
The structural diversity of the bacterial communities in the primary faecal cultures was initially investigated by PCR-DGGE. In samples from the mMCBISO medium inoculated with S. equolifaciens or S. isoflavoniconvertens, faint bands corresponding to the expected size for these species were appreciated (Additional file 1). Clear differences were seen in the bacterial DGGE patterns of the cultures for the different women. However, no changes were observed in the profiles for the mMCB and mMCBISO faecal cultures from the same woman. This suggested to us that the presence of isoflavones did not affect under the experimental culture conditions employed the dominant bacterial populations. Therefore, the bacterial communities of the secondary cultures were not further analysed using this technique.
Bacterial community analysis by high-throughput DNA sequencing
The application of the Illumina technology allowed obtaining an average of 111,542 high quality PCR-amplicon sequences for each replicate. The number of high quality sequences among all the replicates ranged between 71,138–150,292, and estimated sample coverage was considered to be good for all libraries. The mean number of quality sequences retrieved per faecal culture, together with the numbers of OTUs, as defined by ≥97% sequence identity, and α-diversity indexes calculated for each sample are summarized in Additional file 2. To allow further comparative analyses between samples, each 16S rRNA gene amplicon library was rarefied to 35,000 sequences per replicate. The rarefaction curves for the normalized sequences showed that the secondary culture for subject WC (which did not render equol in vitro) was the least diverse (Fig. 2). Further, the numbers of OTUs obtained in all cultures (primary culture in mMCB, primary culture in mMCBISO, and secondary culture in mMCBISO) from the equol non-producer subject WE were smaller (2000–2300 OTUs) than for all other cultures (3000–3500 OTUs) (Additional file 2). Compared to the faecal cultures in mMCB, the presence of isoflavones had no significant influence on the richness index Sobs (number of observed species) (Fig. 3). However, a significant reduction of the bacterial diversity was recorded for the secondary cultures (p < 0.01).
Sequences were then taxonomically assigned using the RDP classifier. A majority of sequences in all samples belonged to the phylum Firmicutes (mean 55%), while sequences from Bacteroidetes (27%), Actinobacteria (9%) and Proteobacteria (7%) phyla were less abundant. At this taxonomic range, only 0.1% of the sequences remained unclassified/unassigned. At the genus-level, between 52 and 75 different genera were identified within each faecal culture. Among the dominant genera Bacteroides presented the highest percentage with a relative abundance between 10 and 20% of assigned reads; Faecalibacterium ranged between 5 and 10%; and sequences related to the Clostridium “sensu stricto” group also reached these percentages (5–10%), but only in the primary cultures from the equol non-producing woman (Fig. 4). Sequences assigned to the genus Collinsella also ranged within these abundances in the primary cultures from the equol producers.
Subdominant genera, varying in percentage from 1 to 5%, were represented by Parabacteroides, Ruminococcus and Dorea; the latter only in the primary cultures from equol producers. In agreement with the DGGE results, the bacterial composition profile obtained by 16S rDNA sequencing was quite similar for equivalent faecal cultures when grown with and without isoflavones (Fig. 4). Nonetheless, changes in subdominant, minor and assumingly rare genera were revealed by the high throughput DNA sequencing technique.
A heatmap depicting relative abundances across samples (including primary and secondary cultures) is shown in Fig. 5. In general, primary faecal cultures in mMCB and mMCBISO from the same woman clustered together suggesting similar microbial profiles, as denoted previously by DGGE. On the other hand, the secondary cultures, with major changes in the bacterial communities’ patterns, plotted in a separate branche.
The application of the Metastats method (FDR correction 0.25) revealed that the presence of isoflavones in the culture medium altered the abundance at several taxonomic hierarchies. For example, the phylum Bacteroidetes was less represented in the presence of isoflavones (p = 0.001), in particular the family Bacteroidaceae with a decrease of reads from the genera Bacteroides and Parabacteroides (Additional file 3). On the contrary, the representation of the family Ruminococcaceae significantly increased with isoflavones (p = 0.011). At the genus level, increases were recorded for Roseburia and Odoribacter sequences. Reads of these two genera were still significantly more abundant when excluding from analysis the faecal culture from the equol non-producing woman (WE). If omitting samples from this woman, and by comparing faecal primary cultures with secondary cultures in the presence of isoflavones, two genera belonging to the family Coriobacteriaceae (capable of synthesizing equol) were significantly increased, Slackia (p = 0.019) and Eggerthella (p = 0.030). Comparison of the microbial communities in the mMCBISO primary cultures for equol producers with that of the non-producer revealed a significant increase in the number of sequences belonging to 28 different genera (Additional file 4). The same number of genera showed an increased in relative abundance when comparing secondary cultures producing equol with those that did not (Additional file 5). Fifteen of the 28 genera (Collinsella, Faecalibacterium, Dorea, Suterella, Clostridium (groups XIVa and XIVb), Alistipes, Oscillibacter, Barnesiella, Coprococcus, Finegoldia, Butyricicoccus, Asaccharobacter, Murdochiella, Allisonella and Odoribacter) were found in both the primary and their derived secondary cultures, suggesting these groups to be more abundant in cultures yielding equol. Of note was the appearance of Collinsella (a member of the family Coriobacteriaceae not described as equol producer), which accounted for 6–15% of sequences in cultures that produced equol, while reads for this genus in isoflavone-containing cultures that did not yield equol was less than 0.04%. Furthermore, reads of the genus Asaccharobacter, a member of the same family harbouring equol-producing strains, increased significantly in all cultures producing equol.
To visualize the overall differences between the bacterial community structures, multivariate statistical analyses, i.e., NMDS and PCoA were applied. This revealed that the bacterial communities of the women’s faecal primary and secondary cultures clustered together, with the exception of that of the secondary culture for woman WC (which, as stated above, showed a reduced diversity and did not produce equol) (Fig. 6). The cluster for the microbiota of this latter sample and those for the primary/secondary cultures for the equol non-producing subject (WE), lay at an appreciable distance from the cultures derived from equol producers (stress value 0.206 in NMDS). Further, the bacterial communities of the mMCBISO cultures for the equol- producing women showed more similarity to one another than to those grown in mMCB. AMOVA revealed the clustering of the two types of samples (in mMCBISO and mMCB) to be significantly different (p = 0.036).
Finally, OTUs significantly different in abundance between primary faecal cultures grown in in mMCBISO and in mMCB were manually assigned to the species level (sequence identity higher that 97%). In cultures grown with isoflavones, OTUs belonging to Faecalibacterium prausnitzii, Eubacterium halli, Subdoligranulum variabile, Ruminococcus flavefaciens/callidus, Blautia obeum and Bacteroides xylanisolvens/ovatus were increased (Additional file 6).
In this work, the bacterial diversity and metabolic changes in cultures of faecal samples from isoflavone-treated menopausal women were investigated. For that purpose, an in vitro anaerobic batch culture system in a modified medium for colonic bacteria supplemented or not with isoflavones was developed. Production of equol was determined with a recently described UHPLC method  and verified by inoculating the isoflavone-containing medium with S. equolifaciens, an equol-producing bacterium. Similar levels of equol production to those obtained with this control strain were recorded in the primary faecal cultures for women with an equol producer phenotype (as determined previously by equol urinary excretion ) - but not in the faecal cultures from a non-equol producer. This strongly indicates differences among the faecal samples in types and/or numbers of intestinal bacteria with capacity to produce equol. The stability of isoflavone bioconversion into equol was further tested by making subcultures of the faecal primary cultures in fresh media. In the presence of isoflavones, equol production was maintained in all except one of these secondary cultures, suggesting that, under the experimental conditions, bacterial types capable of metabolizing isoflavones and producing equol maintain viability and the equol-producing capability, which opens possibilities for a rational propagation and selection of these bacteria. It should be noted that with the UHPLC methodology used in this study other daidzein-derived metabolites that could be formed in the faecal cultures, such as hydroxydadidzein (DHD) and O-desmethylangolensin (O-DMA), might have passed undetected.
In our study, the bacterial composition of the faecal primary and secondary cultures was characterized by two culture-independent methods. PCR-DGGE analysis revealed no major changes in the bacterial communities between the primary cultures with and without isoflavones. This would suggest that the dominant gut populations are a priori neither involved in the metabolism of isoflavones nor influenced by the presence of these compounds. However, the 16S rDNA sequencing data revealed significant changes in different taxa in the presence of isoflavones, particularly with respect to minority groups, which cannot be tracked by the DGGE technique. Although the diversity observed in the faecal cultures was similar in the media with and without isoflavones, a larger number of OTUs was scored in samples yielding equol. This suggests a more diverse bacterial community coming from the faeces of women with an equol producer status.
Collinsella was one of the genera that increased in abundance considerably, both in the primary and secondary cultures that produced equol. This genus belongs to Coriobacteriaceae, a family that harbours most of the currently described equol-producing bacteria . Increases in the number of Collinsella sequences in response to isoflavone treatment have recently been reported in vivo . Although to a lesser extent, the number of reads associated with the genus Asaccharobacter also increased significantly in the faecal cultures producing equol. This agrees with the fact that the single species of this genus, Asaccharobacter celatus, has been reported to be an equol producer . Other intestinal bacteria recently associated with isoflavone metabolism in humans, such as Dorea and Finegoldia [28, 30], were also significantly increased in all cultures producing equol. In agreement with the present results, isoflavones have further been described to stimulate in vivo majority microorganisms of the Clostridium clusters XIVa and IV, to the latter of which F. prausnitzii belongs [31, 14]. On the other hand, like other polyphenols, isoflavones may have antimicrobial activity, which could modulate the diversity and composition of the faecal bacterial communities . In this work, addition of isoflavones to the faecal cultures drove to a decrease of some taxa. In particular, the genera Bacteroides and Parabacteroides reduced their numbers in the presence of these compounds.
Papers reporting production of equol in slurry cultures inoculated with faeces from equol producing subjects have been published previously [12, 13, 30]. However, as far as we are aware, this is the first study that making use of a soy isoflavone extract (commercial formulation) added directly to a colonic culturing medium, describes in deep the bacterial changes and response to isoflavone enrichment by high-throughput 16S DNA sequencing. Although this, the presence in the supplement of other minor components apart from isoflavones (excipients) that might influence the faecal bacterial communities cannot be completely ruled out.
Production of SCFAs (relevant gut bacterial metabolites) was determined in the batch faecal cultures in relation with the addition of isoflavones and the production of equol. SCFAs production in the colon results from bacterial fermentation of dietary starches, fibre and sugars . The SCFAs acetic (C2), propionic (C3) and butyric (C4) acids are the main fermentation end-products in the gut, and are thought to have anti-inflammatory and anti-carcinogenic activities . In contrast, BCFAs, such as isobutyric and isovaleric acids, often associated with protein breakdown, have been less studied. The present data reveal changes in SCFA production in the presence of isoflavones. The increased quantities of the main SCFAs (acetic, propionic and butyric) and BCFAs in the faecal cultures with isoflavones indicates microbial activity leading to the production of these compounds. The soy isoflavone-glycosides can be initially deconjugated by endogenous intestinal bacteria with hydrolytic activity (probably β-glucosidases), as it was observed in this study with the reduction in the daidzin levels after cultivation. This activity releases isoflavone aglycones , but also sugars that may provide an alternative energy source to the cultures. This would agree with some studies involving prebiotic-related compounds such as beta-glucan or cellobiose, which have been shown to increase BCFAs production [35, 36]. At the genus level, Roseburia and Odoribacter sequences were both significantly increased in mMCBISO as compared to the medium without isoflavones. Reduction in members of these two genera has been suggested to enhance inflammation by decreasing SCFA availability [37, 38]. Although not statistically significant, butyrate production in this work was greater in the faecal cultures that produced equol. Butyrate-producing bacteria such as Roseburia, E. hallii (both from Clostridium cluster XIVa) and F. prausnitzii were more abundant when faecal cultures were incubated with isoflavones, particularly in those that produce equol. When significant OTUs were assigned manually at the species level other cellulolytic butyrate producers, such as R. flavefaciens and S. variabile (both of the Clostridium cluster IV) , and propionate producers from the Clostridium cluster XIVa, such as B. obeum, were also identified.
Cross-feeding mechanisms among the more diverse faecal microbiota in equol producers might further be partly responsible for the increased production of propionic and butyric acids in isoflavone-containing cultures. Acetate is the predominant SCFA produced by bacterial fermentation in the gut, and both, F. prausnitzii and E. hallii (the sequences of which became more abundant when cultures were grown in the presence of isoflavones) are well known acetate-utilizing, butyrate-producing species .
Finally, isoflavone metabolism might also be affected by the presence of SCFAs; i.e., butyric acid has been shown to enhance equol production by A. celatus , and an increase production of equol by a mixed microbial culture from human faeces in the presence of propionate and butyrate has also been reported .
Faecal cultures from women with an equol producer phenotype yielded equol in vitro in faecal fermentations with isoflavones, suggesting the presence in faeces of active bacterial types able to produce this compound. The application of high-throughput DNA sequencing of 16S rRNA gene PCR amplicons revealed the composition of the faecal bacterial communities to be modified by the presence of isoflavones, including increases in equol-producing taxa. Although limited by the low number of cultures, isoflavones seem also to promote, at least in vitro, SCFAs production. This increase might be due to the growth stimulation of specific SCFA-producing bacterial types from the Ruminococaceae (members of the Clostridium cluster IV) and Lachnospiraceae (Clostridium cluster XIVa) families.
Analysis of molecular variance
Analysis of similarities
Branched chain fatty acids
Denaturing gradient gel electrophoresis
European Food Safety Authority
False discovery rate
modified medium for colonic bacteria
- mMCBISO :
modified medium for colonic bacteria supplemented with isoflavones
Non-metric multidimensional scaling
Operational taxonomic units
Phosphate buffered saline
Principal coordinates analysis
Ribosomal database project
Short-chain fatty acids
Standard operating procedure
Ultra-high-performance liquid chromatography
Pilšáková L, Riečanský I, Jagla F. The physiological actions of isoflavone phytoestrogens. Physiol Res. 2010;59:651–64.
Bolaños R, Del Castillo A, Francia J. Soy isoflavones versus placebo in the treatment of climacteric vasomotor symptoms: systematic review and meta-analysis. Menopause. 2010;17:660–6.
EFSA (European Food Science Authority). Scientific Opinion on the substantiation of health claims related to soy isoflavones and maintenance of bone mineral density (ID 1655) and reduction of vasomotor symptoms associated with menopause (ID 1654, 1704, 2140, 3093, 3154, 3590) (further assessment) pursuant to Article 13(1) of Regulation (EC) No 1924/2006. EFSA Journal. 2012;10:2847.
Crozier A, Jaganath IB, Clifford MN. Dietary phenolics: chemistry, bioavailability and effects on health. Nat Prod Rep. 2009;26:1001–43.
Setchell KD, Clerici C. Equol: history, chemistry, and formation. J Nutr. 2010;140:1355S–62S.
Atkinson C, Frankenfeld CL, Lampe JW. Gut bacterial metabolism of the soy isoflavone daidzein: exploring the relevance to human nutrition. Exp Biol Med. 2005;230:155–70.
Clavel T, Mapesa JO. Phenolics in human nutrition: importance of the intestinal microbiome for isoflavone and lignan bioavailability. In: Ramawat KG, Merillon JM, editors. Natural products. Amsterdam: Elsevier; 2013. p. 2433–63.
Sánchez-Calvo JM, Rodríguez-Iglesias MA, Molinillo JMG, Macías FA. Soy isoflavones and their relationship with microflora: beneficial effects on human health in equol producers. Phytochem Rev. 2013;12:979–1000.
Yuan JP, Wang JH, Liu X. Metabolism of dietary soy isoflavones to equol by human intestinal microbiota. Implications for health. Mol Nutr Food Res. 2007;51:765–81.
Kemperman RA, Bolca S, Loger LC, Vaughan EE. Novel approaches for analyzing gut microbes and dietary polyphenols: challenges and opportunities. Microbiology. 2010;156:3224–31.
Rafii F, Davis C, Park M, Heinze TM, Beger RD. Variations in metabolism of the soy isoflavonoid daidzein by human intestinal microfloras from different individuals. Arch Microbiol. 2003;180:11–6.
Atkinson C, Berman S, Humbert O, Lampe JW. In vitro incubation of human feces with daidzein and antibiotics suggests interindividual differences in the bacteria responsible for equol production. J Nutr. 2004;134:596–9.
Gardana C, Canzi E, Simonetti P. The role of diet in the metabolism of daidzein by human faecal microbiota sampled from Italian volunteers. J Nutr Biochem. 2009;20:940–7.
Guadamuro L, Delgado S, Redruello B, Flórez AB, Suárez A, Martínez-Camblor P, Mayo B. Equol status and changes in faecal microbiota in menopausal women receiving long-term treatment for menopause symptoms with a soy-isoflavone concentrate. Front Microbiol. 2015;6:777.
Rowland IR, Wiseman H, Sanders TA, Adlercreutz H, Bowey EA. Interindividual variation in metabolism of soy isoflavones and lignans: influence of habitual diet on equol production by the gut microflora. Nutr Cancer. 2000;36:27–32.
van der Meulen R, Adriany T, Verbrugghe K, De Vuyst L. Kinetic analysis of bifidobacterial metabolism reveals a minor role for succinic acid in the regeneration of NAD+ through its growth-associated production. Appl Environ Microbiol. 2006;72:5204–10.
Redruello B, Guadamuro L, Cuesta I, Álvarez-Buylla JR, Mayo B, Delgado S. A novel UHPLC method for the rapid and simultaneous determination of daidzein, genistein and equol in human urine. J Chromatogr B Analyt Technol Biomed Life Sci. 2015;1005:1–8.
Salazar N, Gueimonde M, Hernández-Barranco AM, Ruas-Madiedo P, de los Reyes-Gavilán CG. Exopolysaccharides produced by intestinal Bifidobacterium strains act as fermentable substrates for human intestinal bacteria. Appl Environ Microbiol. 2008;74:4737–45.
Zoetendal EG, Heilig HG, Klaassens ES, Booijink CC, Kleerebezem M, Smidt H, et al. Isolation of DNA from bacterial samples of the human gastrointestinal tract. Nat Protoc. 2006;1:870–3.
Delgado S, Cabrera-Rubio R, Mira A, Suárez A, Mayo B. Microbiological survey of the human gastric ecosystem using culturing and pyrosequencing methods. Microb Ecol. 2013;65:763–72.
Muyzer G, de Waal EC, Uitterlinden AG. Profiling of complex microbial populations by denaturing gradient gel electrophoresis analysis of polymerase chain reaction-amplified genes encoding for 16S rRNA. Appl Environ Microbiol. 1993;59:695–700.
Li J, Rui J, Pei Z, Sun X, Zhang S, Yan Z, et al. Straw-and slurry-associated prokaryotic communities differ during co-fermentation of straw and swine manure. Appl Microbiol Biotechnol. 2014;98:4771–80.
Caporaso JG, Lauber CL, Walters WA, Berg-Lyons D, Lozupone CA, Turnbaugh PJ, et al. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc Natl Acad Sci U S A. 2011;108(Suppl 1):4516–22.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.
White JR, Nagarajan N, Pop M. Statistical methods for detecting differentially abundant features in clinical metagenomic samples. PLoS Comput Biol. 2009;5:e1000352.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Series B Stat Methodol. 1995;57:289–300.
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:e108924.
Minamida K, Ota K, Nishimukai M, Tanaka M, Abe A, Sone T, et al. Asaccharobactercelatus gen. Nov., sp. nov., isolated from rat caecum. Int J Syst Evol Microbiol. 2008;58:1238–40.
Decroos K, Vanhemmens S, Cattoir S, Boon N, Verstraete W. Isolation and characterisation of an equol-producing mixed microbial culture from a human faecal sample and its activity under gastrointestinal conditions. Arch Microbiol. 2005;183:45–55.
Clavel T, Fallani M, Lepage P, Levenez F, Mathey J, Roche V, et al. Isoflavones and functional foods alter the dominant intestinal microbiota in postmenopausal women. J Nutr. 2005;135:2786–92.
Ríos-Covián D, Ruas-Madiedo P, Margolles A, Gueimonde M, de los Reyes-Gavilán CG, Salazar N. Intestinal short chain fatty acids and their link with diet and human health. Front Microbiol. 2016;7:185.
Flint HJ, Duncan SH, Scott KP, Louis P. Links between diet, gut microbiota composition and gut metabolism. Proc Nutr Soc. 2015;74:13–22.
Franke AA, Lai JF, Halm BM. Absortion, distribution, metabolism, and excretion of isoflavonoids after soy intake. Arch Biochem Biophys. 2014;59:24–8.
Nilsson U, Johansson M, Nilsson A, Björck I, Nyman M. Dietary supplementation with beta-glucan enriched oat bran increases faecal concentration of carboxylic acids in healthy subjects. Eur J Clin Nutr. 2008;62:978–84.
van Zanten GC, Krych L, Röytiö H, Forssten S, Lahtinen SJ, Abu Al-Soud W, et al. Synbiotic Lactobacillus acidophilus NCFM and cellobiose does not affect human gut bacterial diversity but increases abundance of lactobacilli, bifidobacteria and branched-chain fatty acids: a randomized, double-blinded cross-over trial. FEMS Microbiol Ecol. 2014;90:225–36.
Duncan SH, Louis P, Flint HJ. Cultivable bacterial diversity from the human colon. Lett Appl Microbiol. 2007;44:343–50.
Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, et al. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012;13:R79.
The technical assistance of B. Redruello and A.M. Hernández-Barranco with the chromatographic analyses is appreciated.
This research was funded by projects from the Spanish Ministry of Economy and Competitiveness (MINECO) (Ref. AGL2011–24300 and AGL-2014-57,820-R) and Asturias Principality (Ref. GRUPIN14-137) . L. Guadamuro was supported by a research grant of the FPI Program from MINECO (Ref. BES-2012-062502) and a FPI Short Stays Grant (Ref. EEBB-I-15-10507).
Availability of data and materials
All relevant data are within the paper and its additional files. The raw sequencing data are deposited in the Sequence Read Archive (SRA) of the NCBI (http://www.ncbi.nlm.nih.gov) database under accession no. SRP064029.
SD and BM contributed with the conception and design of the study. BM and CCT provided material and human resources. LG and SD were involved in all steps of the faecal cultures determinations. CCT was in charge of the experimental design for the high-throughput sequencing. ABD participated in the analysis of PCR-generated 16S rRNA gene libraries and interpretation of the data. SD drafted the manuscript. BM and CCT performed a critical revision of the manuscript. All authors approved the final version of the article.
LG, BM and SD work at the Department of Microbiology and Biochemistry, Instituto de Productos Lácteos de Asturias (IPLA), Consejo Superior de Investigaciones Científicas (CSIC), Paseo Río Linares s/n, 33300 Villaviciosa, Asturias, Spain.
ABD and CCT work at the Thünen-Institut für Biodiversität, Bundesallee 50, 38116 Braunschweig, Germany, where LG performed a short stay as part of her PhD formative training.
The authors declare that they have no competing interest.
Consent for publication
Ethics approval and consent to participate
Human participants providing stool samples for this research signed written informed consent. Approval was obtained from the Research Ethics Committee of the Asturias Principality, Spain (approval number: 15/2011).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Effect of isoflavones on dominant bacterial populations determined by DGGE. PCR-DGGE profiles of the primary faecal cultures from equol-producing women grown in modified medium for colonic bacteria supplemented (mMCBISO) or not (mMCB) with isoflavones; A) WC samples, B) WG samples. M: DGGE marker [comprising the species Slackia isoflavoniconvertens (1) and Slackia equolifaciens (2)]. (PPTX 164 kb)
Summary of sequence processing. Number of quality reads, sample coverage, richness estimators and diversity indexes of 16S rDNA libraries of primary and secondary faecal cultures. (DOCX 20 kb)
Effect of isoflavones in microbial abundance. Families and genera showing significant (p value <0.05) increases (grey) and decreases in their relative abundances (% sequences) when comparing primary faecal cultures in medium with and without isoflavones. (DOCX 16 kb)
Differences in microbial genera associated with equol production in primary faecal cultures. Genera showing significant increases (p value <0.05) in their relative abundances (% sequences) in medium with isoflavones when comparing primary cultures from non-producer and equol producer women. (DOCX 16 kb)
Differences in microbial genera associated with equol production in secondary faecal cultures. Genera showing significant increases (p value <0.05) in their relative abundances (% sequences) in medium with isoflavones when comparing secondary cultures that rendered equol production with those that did not. (DOCX 17 kb)
Effect of isoflavones in microbial species abundance. Identification of OTUs (3% distance level) showing significant increases (p value <0.05) in their relative abundances (%) in primary faecal cultures in mMCBISO as compared to that in mMCB. (DOCX 17 kb)