Next Article in Journal
Staphylococcal Enterotoxin A Induces Intestinal Barrier Dysfunction and Activates NLRP3 Inflammasome via NF-κB/MAPK Signaling Pathways in Mice
Next Article in Special Issue
Development of an Extraction Method of Aflatoxins and Ochratoxin A from Oral, Gastric and Intestinal Phases of Digested Bread by In Vitro Model
Previous Article in Journal
Trial Assay for Safe First-Aid Protocol for the Stinging Sea Anemone Anemonia viridis (Cnidaria: Anthozoa) and a Severe Toxic Reaction
Previous Article in Special Issue
Survey of Freshly Harvested Oat Grains from Southern Brazil Reveals High Incidence of Type B Trichothecenes and Associated Fusarium Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxins
Volume 14
Issue 1
10.3390/toxins14010028
toxins-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mycotoxin Interactions along the Gastrointestinal Tract: In Vitro Semi-Dynamic Digestion and Static Colonic Fermentation of a Contaminated Meal

by
Maria Madalena Costa Sobral
1,*,
Tiago Gonçalves
1,
Zita E. Martins
1,
Christine Bäuerl
2,
Erika Cortés-Macías
2,
Maria Carmen Collado
2 and
Isabel M. P. L. V. O. Ferreira
1,*
1
LAQV/REQUIMTE—Laboratório de Bromatologia e Hidrologia, Faculdade de Farmácia da Universidade do Porto, 4050-313 Porto, Portugal
2
IATA-CSIC—Institute of Agrochemistry and Food Science-National Research Council, Agustín Escardino 7, 46980 Valencia, Spain
*
Authors to whom correspondence should be addressed.
Toxins 2022, 14(1), 28; https://doi.org/10.3390/toxins14010028
Submission received: 28 November 2021 / Revised: 16 December 2021 / Accepted: 23 December 2021 / Published: 1 January 2022
(This article belongs to the Special Issue Mycotoxin Contamination and Food Safety)
Browse Figures
Versions Notes

Abstract

:
Aflatoxin B1 (AFB1) and ochratoxin A (OTA) naturally co-occur in several foods, but no studies have followed the fate of mycotoxins’ interactions along the gastrointestinal tract using in vitro digestion models. This study used a novel semi-dynamic model that mimics gradual acidification and gastric emptying, coupled with a static colonic fermentation phase, in order to monitor mycotoxins’ bioaccessibility by the oral route. AFB1 and OTA bioaccessibility patterns differed in single or co-exposed scenarios. When co-exposed (MIX meal), AFB1 bioaccessibility at the intestinal level increased by ~16%, while OTA bioaccessibility decreased by ~20%. Additionally, a significant increase was observed in both intestinal cell viability and NO production. With regard to mycotoxin–probiotic interactions, the MIX meal showed a null effect on Lactobacillus and Bifidobacterium strain growth, while isolated AFB1 reduced bacterial growth parameters. These results were confirmed at phylum and family levels using a gut microbiota approach. After colonic fermentation, the fecal supernatant did not trigger the NF-kB activation pathway, indicating reduced toxicity of mycotoxins. In conclusion, if single exposed, AFB1 will have a significant impact on intestinal viability and probiotic growth, while OTA will mostly trigger NO production; in a co-exposure situation, both intestinal viability and inflammation will be affected, but the impact on probiotic growth will be neglected.
Keywords:
aflatoxin B1; ochratoxin A; gut microbiota; inflammation; in vitro models
Key Contribution: This study configures the first approach to determining mycotoxins’ fate and their interactions along the gastrointestinal tract by simulating human ingestion of a contaminated meal. AFB1 and OTA’s time of residence in the stomach is short, being more rapidly bioaccessible in the intestine, where their toxicological effects are well acknowledged. This study shows that mycotoxin release in single- or co-exposure scenarios differently impacts their intestinal health outcomes.

1. Introduction

Mycotoxins are toxic metabolites usually found in several food products as a result of fungi spoilage. Aflatoxin B1 (AFB1) and ochratoxin A (OTA) are mycotoxins with recognized human toxicity, being classified as carcinogenic (group 1) and possibly carcinogenic (group 2B) to humans, respectively, by the International Agency for Research on Cancer (IARC) [1]. AFB1 and OTA are prevalent mycotoxins found in a wide range of food products, including cereals and dairy products [2,3,4]. In addition, as these two mycotoxins co-occur naturally in several foods [5,6,7] and as meals combine food products with distinct origins, multiple-mycotoxin contamination is a real situation. Indeed, the European Food Consumption Validation (EFCOVAL) project recently revealed that 600 individuals from five European countries were exposed to mixtures of mycotoxins or their metabolized forms by assessing their content in urine and serum [8]. AFB1 and OTA were among the detected mycotoxins at concentrations ranging from 0.001 to 10 ng/L and 0.1 to 100 ng/L, respectively, in urine and serum samples.
Understanding the fate of contaminants along the gastrointestinal tract (GIT) and interactions with nutrients is a challenge. During the past decade, in vitro digestion methodologies have facilitated the simulation of in vivo conditions. The international INFOGEST network published two distinct in vitro methodologies. The first, in 2014, is a static approach that mimics the physiological conditions of human digestion. This recommends the use of a constant ratio between food and enzyme/electrolytes and a constant pH, time, and temperature [9]. Although useful, it lacks other relevant physiological conditions, such as the gradual acidification of the stomach and gastric emptying [10]. Thus, in 2020, a standardized semi-dynamic digestion model was published [10], considering also the energy value of food/meal for the simulation of the gastric phase, since it impacts the time of residence in the stomach.
No studies have been found concerning the release of mycotoxins from the stomach to the duodenum nor the impact of the gradual acidification on their matrix release. Thus, this study introduces the first approach to using a semi-dynamic model to evaluate AFB1 and OTA release during GI digestion, as well as mycotoxin–food interactions. Moreover, AFB1 and OTA bioaccessibility is also influenced by their water solubility. AFB1 and OTA have a lipophilic character with a coefficient of partitioning of 1.23 and 4.74, respectively, which suggests that part of these mycotoxins will not become bioaccessible during GI digestion, reaching the colon. In addition, some food constituents (i.e., fibers, proteins) have shown to impact mycotoxins’ bioaccessibility [11]. Concerning the bioaccessibility of mycotoxins along the GIT, the few studies published over the past years focus on a specific food matrix (e.g., meat, cereals, fish) [5,12,13], with only two studies evaluating mycotoxins’ bioaccessibility using a meal approach [14,15]. However, none of these reports explores the interactions of mycotoxins after ingestion, and no information has been found regarding AFB1 and OTA interactions during the digestion process of a snack meal.
Beyond digesting foods and allowing nutrient uptake, the GIT also plays an important role in the immune system, which can be compromised by mycotoxins’ presence, as these have shown to damage the intestinal epithelium and contribute to oxidative stress [16]. AFB1 and OTA are known to increase the permeability of the epithelial barrier and decrease epithelial cells’ viability [16]; these mycotoxins have also been shown to modulate the host’s immune system by suppressing immune functions: (i) AFB1 alone or together with other AFs has been shown to reduce reactive oxygen species (ROS) and nitric oxide (NO) production; (ii) in the same way, OTA reduced the gene expression of nuclear factor kB (NF-kB) and several cytokines in the duodenum after piglets’ 30-day exposure to sub-chronic OTA levels (0.05 mg/kg feed) [17].
The gut microbiota plays a crucial role in human health by maintaining immune and metabolic homeostasis; however, some external factors, including exposure to xenobiotics, may imbalance the gut bacteria composition, leading to unhealthy outcomes [18,19]. In fact, it is believed that continuous exposure to toxins benefits the predominance of harmful bacteria over beneficial ones, which may contribute to the development of several disorders, such as obesity, autoimmune diseases, and infections [20]. Exposure of the gut microbiota to mycotoxins may trigger two different scenarios: (i) modulation of the gut microbiota profile due to mycotoxins’ antimicrobial properties or toxic effects toward bacteria or (ii) changing of mycotoxins’ toxicity (increasing or decreasing) as result of degradation/metabolization by some fecal bacteria [21]. Several animal studies have shown that AFB1 and OTA are able to change the microbiota profile, having a significant impact on taxonomic abundance at family and genus levels, as reviewed by Guerre [21]. Specifically, exposure of mice to OTA caused marked changes in gut microbial diversity, decreasing the relative abundance of Firmicutes at the phylum level [22]. In contrast, Lactobacillus casei Shirota was suggested as a detoxifying player when treating rats [23], given their ability to reduce AFB1 toxicity. In contrast, Wang et al. [24] suggested that AFB1 could alter the gut microbiota in a dose-dependent manner, decreasing phylogenic diversity.
The literature shows studies focused on gastric and duodenal bioaccessibility of mycotoxins using in vitro static methods and on numerous food matrices [5,12,15,25,26,27], but the fate of mycotoxins by simulating GI digestion of a contaminated snack meal has not been assessed. Moreover, mycotoxin–microbiota interaction studies found in the literature are animal-based studies and do not relate to the human gut microbiota [21]. An in vitro semi-dynamic digestion methodology followed by static in vitro colonic fermentation was used to assess (i) the bioaccessibility of isolated AFB1 and OTA or as a mixture during GI digestion of a meal, monitoring their release from the stomach to the duodenum through gastric emptying simulations; (ii) the impact of non-bioaccessible mycotoxins on important probiotics; and (iii) the impact of both bioaccessible and non-bioaccessible mycotoxins on intestinal epithelium viability and inflammation.

2. Results

2.1. Mycotoxins’ Analytical Method Performance

The matrix had a significant effect (>20% or <−20%) on the analytes’ quantification, either enhancing or decreasing their analytical response. The percentage of matrix effect measured according to SANTE guidelines [28] (see Section 4.3) was 71% for AFB1 and 56% for OTA. Therefore, the linearity of the chromatographic method was determined by preparing matrix-matched calibration curves using five points at amounts ranging from 2 to 40 ng in the meal and 1 to 20 ng in digested matrices (Table 1). A good linear response was obtained in all matrices (R2 > 0.993).
The recoveries of mycotoxins were within recommended values (80–120%) at all tested concentrations, with relative standard deviations (% RSD) lower than 25% (Table 2), except for OTA at the highest spiked level (55 ng/mL).
The validated method was then applied to confirm the absence of mycotoxins in the meal, before artificial contamination, and to follow the mycotoxins’ fate and interactions along the gastrointestinal tract.

2.2. Mycotoxins Bioaccessibility during Gastrointestinal Digestion

Isolated AFB1 and OTA were released along GI digestion, exhibiting a maximum cumulative isolated bioaccessibility of 51.6 and 72.4%, respectively, after intestinal digestion (Figure 1a,b). Moreover, the earlier gastric emptyings (GEs; E1 to E3) had the greatest contribution to AFB1 and OTA release into the duodenum, meaning that the time of residence of these mycotoxins in the stomach is short (~30 min), becoming quickly available for absorption or for exerting their toxic effects in the intestine. Figure 1c shows the mycotoxins’ distribution concerning the final content that became bioaccessible in the intestinal phase (BIO fraction) and the content that remained non-bioaccessible (NBIO fraction) in relation to its initial content in 10 g of meal.
Regarding interactions, the release of AFB1 or OTA was significantly affected by the presence of the other mycotoxin. The digestion of a meal contaminated with both mycotoxins promoted a release of AFB1 over the gastric phase significantly (p < 0.01) lower (41.5%) than its isolated bioaccessibility (45.5%). However, in the intestinal phase, AFB1 exhibited significantly (p < 0.01) higher final intestinal bioaccessibility in the meal contaminated with both mycotoxins (67.3%) when compared to its isolated bioaccessibility (51.6%). In contrast, OTA’s bioaccessibility was lower in the meal contaminated with both mycotoxins: in the gastric phase, the bioaccessibility of isolated OTA was 33.6% vs. 27.2% in the meal that contained both mycotoxins. The intestinal bioaccessibility of OTA in all gastric emptyings was significantly lower in the meal contaminated with both mycotoxins, exhibiting a maximum bioaccessibility of 52.8% when its isolated bioaccessibility was 72.4%. Thus, meals co-contaminated with AFB1 and OTA will likely increase AFB1 bioaccessibility by ~16% and decrease OTA bioaccessibility by ~20% at the intestinal level.

2.3. Impact of Bioaccessible Mycotoxins on Intestinal Cells’ Viability and Inflammation

There is an impact of BIO fractions after 3 h exposure of 7-day-differentiated Caco-2 monolayers on cell viability and the formation of nitric oxide and reactive oxygen species (Figure 2).
When compared to the control meal, the digested samples contaminated with AFB1 or both mycotoxins (MIX) significantly increased cell death (Figure 2a) by releasing lactate dehydrogenase from the cytosol into the culture medium, while at the studied doses, bioaccessible OTA did not present a significant impact on cells’ viability. These results were corroborated by strong Pearson correlations (r > 0.6) between LDH measurements in Caco-2 cells and AFB1 contents in the meal before digestion and in the bioaccessible fraction (Table 3).
Nitric oxide formation only differed from the control meal (p < 0.05) in digested meals with OTA or MIX (Figure 2b). These results could be supported by a moderate Pearson correlation (r = 0.535) between NO and OTA content in the meal; however, this correlation was not observed between NO and bioaccessible OTA. When summing up AFB1 and OTA bioaccessible contents, moderate correlations (r > 0.525 and r > 0.446) were observed between AFB1 + OTA bioaccessible content and NO and LDH, respectively (Table 3). The bioaccessible mycotoxins did not impact reactive oxygen species formation (Figure 2c).

2.4. Impact of Non-Bioaccessible Mycotoxins on Probiotic Lactobacillus casei and Bifidobacterium lactis

The impact of different amounts of non-bioaccessible (NBIO) fraction (0.5%, 1%, and 2%) on two important probiotic strains (Lactobacillus casei BL23 and Bifidobacterium lactis NCC 2818) was evaluated to understand how non-bioaccessible mycotoxins affect bacterial growth and to select the best NBIO percentage for colonic fermentation assays. Data are presented in Figure 3 and Table 4.
The NBIO fraction with AFB1, OTA, or MIX at 0.5% did not exert any significant changes in relation to the control, except for the maximal optical density (MOD) of the Bifidobacterium strain in the NBIO fraction containing AFB1 or OTA. The increase of NBIO fractions to 1% or 2% significantly impacted bacterial growth parameters: at 1%, both bacterial growth and lag time of the Lactobacillus strain significantly reduced in the presence of AFB1 and MIX, while OTA only affected (increased) the growth rate; with regard to the Bifidobacterium strain, only the growth rate and MOD were affected by AFB1 and OTA, while MIX did not differ from the control.
The exposure of both bacterial strains to 2% NBIO containing AFB1 or OTA significantly impacted all growth parameters in relation to the control, except for OTA in the case of the growth rate of Bifidobacterium. The mycotoxin type differently impacted bacterial growth: (i) AFB1 significantly reduced the measured parameters, (ii) OTA significantly increased them, while (iii) MIX had a null effect on bacterial growth.
Pearson correlations between the mycotoxin content in NBIO fractions and bacterial growth parameters showed the same opposite behavior between mycotoxins (Table 5):
(i)
Strong negative Pearson correlations between the AFB1 content in the meal and all growth parameters of both L. casei and B. lactis, suggesting a negative impact on bacterial growth (i.e., reduction in the MOD, growth rate, and lag time). Even stronger negative correlations were observed between the non-bioaccessible AFB1 content and probiotic growth parameters, indicating an inhibition of growth in both strains by AFB1.
(ii)
Strong positive correlations between the OTA content in the meal and L. casei growth parameters and the B. lactis MOD. Moderate correlations were observed between the non-bioacessible OTA content and L. casei growth, confirming an enhancement of the growth rate.
(iii)
No correlation was found between bacterial growth and mycotoxin content when summing up their contents (AFB1 + OTA), suggesting a null effect in the co-exposure situation.
Given that the NBIO fraction at 2% showed the greatest impact on bacterial growth, this percentage was selected for colonic fermentation assays.

2.5. Impact of Non-Bioaccessible Mycotoxins on the Gut Microbiota, Intestinal Cells’ Viability, and Inflammation

After colonic fermentation, an increase in the total bacteria content was observed (Figure 4a), but the presence of AFB1, OTA, or MIX did not influence the total bacteria content. However, mycotoxins were able to change the gut microbiota composition at the phylum and family levels. At the phylum level (Figure 4b), MIX meals followed by control meals were the samples that had the highest percentage of Firmicutes and Bacteroidetes, while AFB1 and OTA meals revealed lower abundances. At the family level (Figure 4c), the MIX meal had the highest percentage of Bacteroidaceae, Ruminococcaceae, and Lachnospiraceae, followed by control, AFB1, and OTA meals, respectively. Differently, Lactobacillaceae exhibited higher percentages in the MIX meals, followed by OTA meals, with AFB1 and control meals having a lesser contribution. With regard to Bifidobacteriaceae, the control meal had the highest contribution, followed by AFB1 and MIX meals, with OTA meals having the lower abundances.
Regarding the impact of NBIO fractions on intestinal cells’ viability and inflammation, the NBIO samples before colonic fermentation and the fecal supernatant after colonic fermentation did not significantly trigger an NF-kB inflammatory response on human colon HT29 reporter cells (Figure 5a,b).

3. Discussion

Mycotoxin bioaccessibility values represent the maximum amount that can be absorbed, which can be interpreted as the worst-case scenario of humans’ exposure. These are relevant data regarding mycotoxins’ risk assessment determination; thus, reliable models of the digestive process are used to obtain more accurate and precise information. Herein, a novel semi-dynamic digestion model, followed by static colonic fermentation, was used to understand the mycotoxins’ behavior along the GIT [10]. To the best of our knowledge, this is the first approach combining supervision of mycotoxin release from the stomach to the duodenum, together with the evaluation of the impact of non-bioaccessible mycotoxins on important gut probiotics. Moreover, this study also aimed at assessing the mycotoxins’ fate alone and in a co-exposure situation along the GIT and gain more insights into mycotoxins’ impact on intestinal epithelium viability and inflammation in intestinal and colonic phases.
The final intestinal bioaccessibility of AFB1 or OTA was affected by the presence of the other mycotoxin, increasing AFB1 release and reducing OTA’s. These interactions have already been reported between OTA and patulin during digestion, increasing their bioaccessibility compared to their single-exposure situation [5]. The promotion of AFB1 bioaccessibility represents a risk, since AFB1 belongs to the most dangerous mycotoxins in terms of carcinogenic effects [1]. In addition, the short time of residence in the stomach (~30 min) observed for both mycotoxins indicates fast availability in the duodenum to exert their toxic effects on the intestinal epithelium or on other target organs after absorption (e.g., liver and kidney) [1]. The co-ocurrence of AFB1 and OTA is a real situation with regard to several foods products, including dairy products and cereals [2,7], but limited data are found on their combined bioaccessibility: (i) Kabak et al. [27] reported 90% and 30% bioaccessibility for AFB1 and OTA, respectively, in several food matrices spiked or naturally contaminated, and (ii) Versantvoort et al. [15] reported bioaccessibility above 100% for both mycotoxins after digesting peanut slurry. These results differ from those obtained in this study, probably due to the food matrix itself that has a relevant influence on compounds’ release during digestion and the use of a different in vitro model, as Kabak et al. and Versantvoort et al. used a non-standardized in vitro static digestion model [15,27]. The use of standardized methodologies is of paramount importance to allow comparison of experiments among different research groups, as several parameters may affect the results. For example, the gradual acidification of the stomach affects protein structures and interactions with matrix constituents and/or other compounds [10], which may influence mycotoxin bioaccessibility. Indeed, AFB1 and OTA have a high capacity to establish covalent bonds with proteins [29]. The gradual acidification observed in this study agrees with previous results reported during digestion of semi-solid dairy products [10].
The intestinal epithelium represents the first barrier met by xenobiotics, such as dietary mycotoxins, known to cause several GI disorders [30,31]. Herein, AFB1 and OTA had different effects on intestinal viability and inflammation when compared to the control meal: (i) the digested meal containing AFB1 significantly reduced intestinal cells’ viability, but no effect was observed on inflammation; (ii) the digested meal with OTA did not impact cells’ viability but had a pro-inflammatory effect, increasing NO production; (iii) the meal containing both mycotoxins (MIX) significantly impacted both intestinal cells’ viability and NO production; and (iv) none of the meal samples containing mycotoxins triggered ROS formation. These data are related with each mycotoxin content in the meal or in the bioaccessible fraction, as confirmed by the strong or moderate Pearson correlations obtained for AFB1, OTA, and AFB1 + OTA contents. Thus, the increased LDH values observed in digested AFB1 and MIX meals are explained by the AFB1 content and its recognized intestinal toxicity [32,33]. Therefore, the higher LDH formation in the digested MIX meal is explained by the higher AFB1 content as a result of increased bioaccessibility in the MIX meal. The disruption of the intestinal barrier permeability of Caco-2 cells (reduced cell viability and tight-junction protein expression) was observed after exposure to increasing concentrations of AFB1 (3 to 100 µM) and OTA (1 to 100 µM) [32]. In the case of NO production, bioaccessible OTA and MIX triggered similar pro-inflammatory responses, which could suggest an interaction between AFB1 and OTA, since OTA had lower intestinal bioaccessibility (52.8%) in MIX when compared to the isolated bioaccessibility (72.4%).
Altogether, in a meal with single-mycotoxin contamination (AFB1 or OTA), only 50% of AFB1 will be bioaccessible and the remaining content will follow colonic fermentation, whereas for OTA, more than 70% will become bioaccessible. The type of mycotoxins present in the meal will likely cause different intestinal health impacts. However, in a co-exposure situation, there will be a promotion of AFB1 bioaccessibility up to 70% and OTA reduction of about 50%, tending to cause both loss of intestinal viability and increase in NO production. This means that the amount of AFB1 or OTA that will undergo colonic fermentation is dependent on the presence of other mycotoxins.
Herein, a considerable amount (32.7 to 48.4% of AFB1 and 27.6 to 47.2% of OTA) remained in the non-bioaccessible fraction, representing a different threat to the gut microbiota. AFB1 had a greater impact than OTA on both Lactobacillus casei and Bifidobacterium lactis growth curve parameters, while MIX exhibited a null effect on these probiotics. These effects were strongly correlated with their non-bioaccessible contents: negatively for AFB1, positively for OTA, and not correlated when summing up AFB1 with OTA. As AFB1 and OTA were able to modulate the growth of two important probiotics when using isolated strains, it was necessary to confirm their effects on those probiotics as well as others using a gut microbiota approach. After 24 h colonic fermentation, although AFB1, OTA, and MIX did not impact the total bacterial content, an impact on bacterial composition was observed, confirming their impact on Lactobacillaceae and Bifidobacteriaceae families as well as on other important families belonging to Firmicutes (Ruminococcaceae, and Lachnospiraceae) and Bacteroidetes (Bacteroidaceae) phyla. These data agree with previous reports showing OTA’s ability to shift the microbiota balance, increasing Lactobacillaceae and decreasing Bacteriodaceae. At the genus level, OTA decreased the population of Bacteroides, Dorea, Escherichia, Oribacterium, Ruminococcus, and Syntrophococcus, while it increased the number of Lactobacillus members [34]. With regard to AFB1, although less information was found regarding its impact on the gut microbiota, this mycotoxin has been reported as being able to decrease phylogenetical diversity but increase the evenness of the community composition. Lactic acid bacteria can also be reduced in the presence of AFB1, although there are no changes at the phylum level, as reviewed by Liew and Mohd-Redzwan [34]. When using animal models, differences in gut microbiota abundances have been reported after exposure to AFB1: (i) an increase in Bacteroidetes and a decrease in Firmicutes were also verified by Gorsu [35] after 30-day exposure of piglets to an AFB1 (0.32 mg/kg)-contaminated diet, and (ii) an increase in Lachnospiraceae family abundance was also verified in the mice gut microbiota after exposure to AFB1 [36]. In addition, probiotics from both Lactobacillus and Bifidobacteria strains have been suggested as suitable binders of AFB1 and OTA, reducing their bioaccessibility [27]. Besides mycotoxin-binding capacity, the addition of Lactobacillus casei Shirota in feed contaminated with AFB1 had ameliorative effects on gut microbiota changes induced by AFB1 [23] as increased microbiota richness.
Regarding the impact of non-bioaccessible mycotoxins on intestinal cells’ viability and inflammation, the fecal supernatant resulting from colonic fermentation of AFB1, OTA, and MIX meals did not trigger the NF-kB activation pathway after exposure of HT29 reporter cells, suggesting that the immune system is not compromised by non-bioaccessible mycotoxins. NF-kB is a transcription factor that plays a major role in the regulation of genes responsible for inflammatory and immune responses [37]. Previous studies have shown that mycotoxins such as patulin and deoxynivalenol are able to suppress NF-kB-mediated activation via Toll receptor signaling [38,39], which could compromise the immune system, weakening the response to pathogens. These results suggest that the gut microbiota is able to reduce mycotoxins’ toxicity.
In summary, if a meal is contaminated with a single mycotoxin, in the case of AFB1, 50% bioaccessible AFB1 will have a strong impact on intestinal cells’ viability and the other 50% non-bioaccessible AFB1 on specific probiotic strain growth (L. casei and B. lactis). In the case of OTA, since a higher value will be bioaccessible (~70%), it will trigger NO production at the intestinal level, and the non-bioaccessible fraction will have a less impact on probiotic growth. In a co-exposure situation, higher contents of AFB1 will become bioaccessible, potentiating loss of intestinal viability and increased NO production at the intestinal level, but the impact on probiotic growth will be neglected. After colonic fermentation, the fecal supernatant did not trigger the NF-kB activation pathway, suggesting that the gut microbiota is able to reduce mycotoxins’ toxicity.

4. Materials and Methods

4.1. Sample Preparation and Spiking

Liquid yogurt and corn cookies were obtained from local supermarkets in Porto (Portugal). Snack meals were prepared by mixing 1 yogurt (160 g) with 2 corn cookies (16 g), representing 91 and 9% of the meal composition, respectively. Prior to meal preparation, the corn cookies were ground for 12 s at 8000 rpm (Grindomix GM 200; Retsch GmbH, Haan, Germany) and then the respective amount (16 g) mixed with the yogurt. The total amount of meal (160 g) was used as a real meal proportion to allow downscaling to an in vitro approach using 10 g of the sample.
The same meal preparation (160 g) was used for all experiments to guarantee the same meal nutritional composition, appearance, and consistency. Thus, 30 g of the meal was used for each condition and contaminated with the respective amounts of mycotoxins: (a) control (without mycotoxins), (b) contamination with 9 µg of AFB1, (c) contamination with 9 µg of OTA, and (d) simultaneous contamination with 6 µg of AFB1 and 6 µg of OTA. Duplicates of gastrointestinal digestion were assessed for each condition. After contamination, the real content of mycotoxins was quantified, as described in Section 4.3. Thus, each 10 g of the sample used for digestion had 3.0 ± 0.2 µg of AFB1 or 3.2 ± 0.2 µg of OTA for single exposure and a combination of 2.2 ± 0.3 µg of AFB1 and 2.0 ± 0.2 µg of OTA for the co-exposure situation. Although this level of mycotoxins is higher than the amount usually found in foods, it was selected to counterbalance the dilution factor of digestion experiments and allow the quantification of mycotoxins at all steps of digestion (all gastric and intestinal emptyings). The samples were stored at −20 °C in plastic bags until being used for digestion experiments. Mycotoxin standards OTA (1 mg, purity 97%) and AFB1 (1 mg, purity > 98%) were purchased from Sigma-Aldrich (St. Louis, MO, USA), and individual stock solutions of 5 µg/mL were prepared and used for artificial contamination purposes.

4.2. In Vitro Upper Gastrointestinal Digestion

Upper gastrointestinal digestion composed of oral, gastric, and duodenal phases was performed according to the standardized semi-dynamic protocol described by Mulet-Cabero et al. [10]. The energy value of the meals was calculated based on the macronutrient content found in the package: the yogurt and the corn snacks had 1.6% and 1.9% fat, 3.8% and 8.1% protein, and 10% and 81.6% carbohydrates, respectively. A proportion of 1 yogurt (160 g) plus 2 cookie (16 g) snacks was defined for the study. This proportion resembled a real meal of 176 g and was scaled down using the protocol’s supplementary information. The enzymes (pepsin and pancreatin extract) and the bile extract used in this study were purchased from Sigma-Aldrich. The enzymatic activities were determined prior to digestion experiments, as recommended by Mulet-Cabero, et al. [10].
Briefly, 10 g of the snack meal was mixed with simulated salivary fluid (2.00 mL simulated salivary electrolyte fluid (eSSF), 12.5 µL of 0.3 M CaCl2(H2O)2, and 0.363 mL of ultrapure water) using human saliva as an amylase source (150 U/mL of SSF) for 2 min at 37 °C. Then, the bolus (12.5 mL) was transferred into a 70 mL glass v-form vessel thermo-stated at 37 °C containing 10% of the simulated gastric fluid corresponding to in vivo fasting conditions. The remaining 90% of the simulated gastric fluid (8.75 mL of simulated gastric electrolyte fluid (eSGF), 0.925 mL of 1.5 M HCl, 6.3 µL of 0.3 M CaCl2(H2O)2, 1.57 mL of ultrapure water, and 1.125 mL of enzyme solution (4000 U pepsin/mL SGF)) was added at a constant ratio by separated devices. A printed paddle stirrer was used at 15 rpm for agitation. Five gastric emptyings (GEs) of 5 mL each were performed every 17 min. The gastric phase had a total time of digestion of 85 min. Before each gastric emptying, the vessel content was mixed to make the sampling more accurate. The pH was measured in each emptied aliquot, and the pH was raised to 7 with NaOH 1N (300 µL) to inactivate pepsin and prepare samples for the intestinal phase. The tubes were centrifuged at a low rotation (2000× g, 5 min) to obtain a 1 mL aliquot for mycotoxins analysis in the gastric bioaccessible phase. The remaining amount was kept at 4 °C for intestinal digestion the following day.
The intestinal phase was performed using a static methodology according to the INFOGEST 2.0 protocol [40]. Briefly, to 4.3 mL of neutralized chyme of each GE, 3.7 mL of simulated intestinal fluid (2.5 mL of simulated intestinal electrolyte fluid (eSIF), 8.0 µL of 0.3 M CaCl2(H2O)2, 0.492 mL of ultrapure water, 0.200 mL of pancreatin solution (200 U trypsin/mL SIF), and 0.500 mL of bile extract (20 mmol/L SIF)) was added and incubated at 37 °C for 2 h under agitation (15 rpm). After digestion, the samples were immediately immersed in ice for 10 min, followed by centrifugation at 10,000× g for 5 min at 4 °C. Both bioaccessible (BIO) and non-bioaccesible (NBIO) fractions were collected. Mycotoxin extraction from the BIO fraction was performed on the same day of digestion and the remaining digested volume stored at −20 °C until further analysis. The NBIO fraction was stored at −80 °C until used for colonic fermentation assay.

4.3. Mycotoxin Analysis

To 200 µL of digested samples, 4 µL of the internal standard (IS,10 µg/L of aflatoxin B2, purity 98%, Sigma-Aldrich) was added and vortexed. Then, the sample was mixed with an equal volume (200 µL) of extracting solvent (ACN, w/2% formic acid). Afterward, 50 mg of MgSO4 and 35 mg of NaCl were added, vortexed immediately, and centrifuged at 10,000× g for 5 min. The organic phase was collected in a 2 mL vial, evaporated to dryness under a N2 stream, and re-dissolved in 60 µL of extracting solvent. The extraction of mycotoxins in the solid samples (before digestion) was performed in the same way but with two additional steps: to 500 mg of the samples, 4 µL of IS (10 µg/L) was added and left to homogenize for 30 min in an orbital shaker. Then, 1 mL of ultrapure water was added and homogenized for another 30 min. Thereafter, the extraction followed the same procedure, as mentioned above. In the end, 20 µL of the samples were injected into an HPLC system (Jasco LC PU-1580, Tokyo, Japan) with an auto-sampler (Jasco AS-2057 Plus, Japan), a YMC-triart C18 ExRS column (150 × 3.0 mm, 3 µm, 80 Å), and a pre-column from Phenomenex (C18, 4 × 3.0 mm (i.d.), with the column kept at 40 °C. The mycotoxins were eluted using a flow rate of 0.5 mL/min of acidified water (2% formic acid, eluent A) and acidified acetonitrile (2% formic acid, eluent B) as mobile phases, with the following gradient: 30% B in the first 2 min, increasing to 95% B in 10 min, holding this condition for 5 min, and then decreasing to 30% B in 4 min, followed by 5 min of stabilization (total run time of 26 min). AFB1, AFB2, and OTA were detected by fluorescence (FLD, Jasco FP-920, Japan) using excitation/emission wavelengths of 365/455 nm and 333/465 nm for AFs and OTA, respectively. Borwin version 1.5 Data Software (JMBS, France) was used to analyze and export the data.
The method was evaluated using an in-house validation approach assessing linearity, accuracy, and precision to perform reliable quantitative results according to EU guidelines and ICH recommendations [41,42]. The matrix effect was also evaluated for both mycotoxins according to SANTE guidelines [28].

4.4. Impact of Bioaccessible Mycotoxins on Caco-2 Cells’ Viability and on Inflammation

Caco-2 cells were grown in complete medium (CM) composed of Dulbecco’s Modified Eagle Medium (DMEM), 10% fetal bovine serum (FBS), 1% glutaMAX, 1% non-essential amino acids (NEAAs), and 1% penicillin/streptomycin (all purchased from Gibco/Life Technology, Paisley, United Kingdom) at 37 °C with 5% CO2. When the cells were 80% confluent, they were trypsinized and seeded in 96-well plates (5 × 104 cells/mL) and left for differentiation over 7 days with medium renewal every 2 days. For these assays, the 5 intestinal gastric emptyings were combined to obtain the total mycotoxins available after the intestinal phase. Then, Caco-2 cells were exposed to the digested sample (12× diluted in CM) for 3 h. After the exposure time, the cells were used for ROS determination, while the cells supernatant was used for NO, as described by Pinho et al. [43]. The viability of Caco-2 cells was also determined by measuring the release of lactate dehydrogenase in the cells’ supernatant according to Kaja et al. [44]. The pro-inflammatory effect of bioaccessible mycotoxins (BIO fraction) was carried out by measuring as a marker of inflammation nitric oxide (NO) and reactive oxygen species (ROS) formation on 7-day-differentiated Caco-2 cells [43,44].

4.5. Effect of Non-Bioaccessible Mycotoxins on Probiotics Lactobacillus casei and Bifidobacterium lactis Growth

The effect of mycotoxins that remain in the NBIO fraction after in vitro digestion was checked in two beneficial bacterial strains, Lactobacillus casei BL23 and Bifidobacterium lactis ATCC 27536, to infer the impact of mycotoxins on bacterial growth and choose the ideal proportion of the NBIO fraction for colonic fermentation, given their relevance for the human gut microbiota [45]. These were grown in MRS + 0.05% L-cysteine under anaerobic and static conditions for 24 h at 37 °C in the presence or absence of 20 µL of the sample according to Alcántara et al. [46]. Changes in optical density at 595 nm were registered in a POLARStar reader (BMG Labtech, Ortenberg, Germany), and the strains’ growth data were modeled using the Gompertz equation [47] in order to mathematically describe the microbial growth and compute the specific growth rate in the exponential phase and optical density in the stationary one.

4.6. In Vitro Colonic Fermentation

The NBIO fraction resulting from the duodenal phase followed in vitro colonic fermentation to evaluate the impact of non-bioaccessible mycotoxins on the microbiota and intestinal function. First, a pool of feces from healthy human volunteers (n = 5) was gathered and the fecal microbiota was extracted with saline solution (1:10 proportion (w/v)), homogenized in an orbital shaker over 30 min, and separated from the rest of the fecal material using a Nycodenz® density gradient without affecting the original microbiota composition [48]. Afterward, the fecal microbiota was transferred to colonic fermentation tubes along with sterile basal medium (as described by [49]) and left overnight at 37 °C under anaerobic conditions for stabilization (~24 h). The following day, the NBIO fraction (2%, w/v) was added to the colonic fermentation tubes and left for 24 h under static anaerobic conditions. Duplicates of colonic fermentation were assessed for each sample. At each time point, both supernatant and pellet were collected for further analysis. The quantity and profile of bacteria before colonic fermentation were also assessed.

4.7. Total Bacteria Quantification by qPCR and 16S rRNA Amplicon Sequencing

The DNA was extracted according to fecal-sample-based adaptation to the Maxwell® RSC PureFood GMO and Authentication Kit (Catalo No-AS1600; Promega Corporation, Madison, WI, USA). Briefly, the pellets resulting from colonic fermentation were homogenized with 850 µL of CTAB buffer and transferred to 2 mL tubes together with beads for further lysis using the FastPrep-24 5G homogenizer for 2 cycles at 6.0 m/s for 60 s, with a 30 s pause (MP Biomedicals). The samples were then heated at 95 °C for 5 min (vortex at half time), followed by the addition of Proteinase K (40 µL) and RNAase (20 µL), vortexed, and heated at 70 °C for 10 min. Then, the samples were centrifuged at 13,000 rpm for 5 min. Maxwell® RSC Cartridge preparation and loading occurred, as detailed by the manufacturer. All samples were eluted in 50 μL of the provided elution buffer. The extracted DNA was quantified (ng/µL) using a Qubit 2.0 fluorometer (Life Technology) before storing at −80 °C until further use.
Quantification of total bacteria was carried out by quantitative PCR through amplification and detection of the 16S rRNA gene. In each 96-well plate (Sorenson BioScience, Inc., Murray, UT, USA), 10 µL of the reaction mixture was created by 3.5 µL of PCR water, 0.25 µL of the forward primer (10 µM, 515F (5′-GTGCCAGCMGCCGCGGTAA), 0.25 µL of the reverse primer (10 µM, 789R (5′-GCGTGGACTACCAGGGTATCT), 5 µL of the Light Cycler 480 SYBR Green I Master mix (Roche Life Science), and 1 µL of the sample, standard, or water (negative control) (all performed in duplicate) [50,51]. The Light Cycler 480 Real-Time PCR System (Roche Life Science) was used to amplify the DNA as follows: 95 °C for 10 s, 62 °C for 10 s, and 72 °C for 20 s over 40 cycles, followed by dissociation curve analysis. The standard curve used to calculate the bacterial concentration in each sample was obtained by using tenfold dilutions of specific DNA fragments.
DNA was also used for the bacterial taxonomic composition by 16S rRNA gene sequencing of the V3–V4 region (Illumina platform using a paired-end 2 × 250 bp reading system). Raw sequences were processed in-house by Novogene Bioinformatics Technology Co. Ltd. (Tianjin, China). In brief, paired-end reads were assigned to samples based on their unique barcodes and truncated by cutting off the barcodes and primer sequences. Paired-end reads were merged using FLASH (V1.2.7) [52] (see details at [53]), a fast and accurate analysis tool, which is designed to merge paired-end reads when at least some of the reads overlap the read generated from the opposite end of the same DNA fragment, and the splicing sequences were called raw tags. Quality filtering on the raw tags was performed under specific filtering conditions to obtain high-quality clean tags [54] according to the QIIME (V1.9.0) [55] quality-controlled process. The tags were compared with the reference database (SILVA database, see [56]) using the UCHIME algorithm [57] to detect chimera sequences, and then the chimera sequences were removed. Sequences with ≥97% similarity were assigned to the same operational taxonomic unit (OTU), and taxonomy was assigned with the SILVA database (see details at [56]) for microbial annotation at each taxonomic level (phylum, family, genus levels).

4.8. Impact of Non-Bioaccessible Mycotoxins on HT29 NF-kB Reporter Cells’ Viability and on Inflammation

The pro-inflammatory effect of in vitro colonic fermentation supernatants was carried out, as detailed by Alcántara et al. [46], with some modifications. Briefly, the pro-inflammatory effect of mycotoxins after colonic fermentation of the NBIO fraction was determined using a stably transfected HT29 reporter cell line to measure secreted embryonic alkaline phosphatase (SEAP) activity upon activation of the NF-kB pathway. HT29 reporter cells were grown in DMEM supplemented with 10% FBS, 1% (v/v) sodium pyruvate, 1% P/S, and zeocin (200 µg/mL) at 37 °C and 5% CO2. For this experiment, HT29 cells were seeded in 96-well plates at 6.5 × 104 cells/well and incubated at 37 °C and 5% CO2 for 24 h. The following day, the cells were exposed to fecal supernatants in a 1:10 dilution for 24 h. After incubation, the cell supernatant was collected and SEAP activity measured using p-nitrophenyl phosphate as a phosphatase substrate according to the manufacturer’s instructions (Thermo Fisher Scientific, ref. 37620). Cell viability was measured using the resazurin reduction assay. Briefly, 150 µL of resazurin solution (10 µg/mL) was added to each well and incubated at 37 °C and 5% CO2 for 5 h so that viable cells could reduce resazurin to fluorescent resorufin, whose amount is directly proportional to the number of living cells, and cell viability was measured at 540–15 nm excitation and 590–20 nm emission using a CLARIOstar (BMG Labtech) plate reader.

4.9. Statistical Analysis

Statistical analysis was performed and graphs created using GraphPad Prism 8.1.0 for Mac. The variables’ normality was assessed by the Shapiro–Wilk test. On those variables following normal distribution, Student’s t-test was applied for comparison of AFB1 or OTA bioaccessibility between isolated and mixed digestion, while the Mann–Whitney U test was applied for variables not following normal distribution. For cell culture assays and total bacteria, one-way ANOVA and Dunnett’s post hoc test (parametric) or Kruskal–Wallis and Dunn’s post hoc tests (non-parametric) were applied to compare AFB1, OTA, or MIX effects in comparison with the control meal (without mycotoxins). Pearson correlations were used to test the association between the mycotoxin content in the meal, the bioaccessible fraction, and the non-bioaccessible fraction with results obtained from cell/bacterial assays (cell viability, inflammation, or bacterial growth). Pearson correlations were performed using IBM SPSS statistics software version 17.0.

Author Contributions

Conceptualization, M.M.C.S., I.M.P.L.V.O.F. and M.C.C.; methodology, M.M.C.S., T.G., E.C.-M. and C.B.; software, M.M.C.S., I.M.P.L.V.O.F. and M.C.C.; validation, M.M.C.S., I.M.P.L.V.O.F. and M.C.C.; formal analysis, M.M.C.S., Z.E.M. and I.M.P.L.V.O.F.; investigation, M.M.C.S., T.G., E.C.-M. and C.B.; resources, I.M.P.L.V.O.F. and M.C.C.; data curation, M.M.C.S., Z.E.M., I.M.P.L.V.O.F. and M.C.C.; writing—original draft preparation, M.M.C.S.; writing—review and editing, I.M.P.L.V.O.F. and M.C.C.; supervision, M.M.C.S., I.M.P.L.V.O.F. and M.C.C.; funding acquisition, I.M.P.L.V.O.F. and M.C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundação para a Ciência e a Tecnologia (FCT) in the framework of the project DIETxPOSOME (PTDC/SAU-NUT/6061/2020).

Institutional Review Board Statement

Ethical review and approval were waived for this study by the ethics committee of the Faculty of Pharmacy University of Porto (report no. 03-04-2021) due to the use of biological samples as human saliva for in vitro digestion experiments. The use of human feces was also approved by the local ethics committee of the Atencion Primaria-Generalitat Valenciana (CEIC-APCV).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Acknowledgments

Z.E.M. received support from QREN (NORTE-01-0145-FEDER-000052). M.M.C.S. wishes to thank INFOGEST cost action for the STSM grant (STSM-2021-002).

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Ostry, V.; Malir, F.; Toman, J.; Grosse, Y. Mycotoxins as human carcinogens—The IARC Monographs classification. Mycotoxin Res. 2017, 33, 65–73. [Google Scholar] [CrossRef] [PubMed]
  2. Becker-Algeri, T.A.; Castagnaro, D.; De Bortoli, K.; De Souza, C.; Drunkler, D.A.; Badiale-Furlong, E. Mycotoxins in Bovine Milk and Dairy Products: A Review. J. Food Sci. 2016, 81, R544–R552. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Benkerroum, N. Mycotoxins in dairy products: A review. Int. Dairy J. 2016, 62, 63–75. [Google Scholar] [CrossRef]
  4. Sarmast, E.; Fallah, A.A.; Jafari, T.; Khaneghah, A.M. Occurrence and fate of mycotoxins in cereals and cereal-based products: A narrative review of systematic reviews and meta-analyses studies. Curr. Opin. Food Sci. 2021, 39, 68–75. [Google Scholar] [CrossRef]
  5. Assunção, R.; Martins, C.; Dupont, D.; Alvito, P. Patulin and ochratoxin A co-occurrence and their bioaccessibility in processed cereal-based foods: A contribution for Portuguese children risk assessment. Food Chem. Toxicol. 2016, 96, 205–214. [Google Scholar] [CrossRef]
  6. Ibáñez-Vea, M.; González-Peñas, E.; Lizarraga, E.; de Cerain, A.L. Co-occurrence of aflatoxins, ochratoxin A and zearalenone in barley from a northern region of Spain. Food Chem. 2012, 132, 35–42. [Google Scholar] [CrossRef]
  7. Martins, C.; Assunção, R.; Cunha, S.C.; Fernandes, J.O.; Jager, A.; Petta, T.; Oliveira, C.A.; Alvito, P. Assessment of multiple mycotoxins in breakfast cereals available in the Portuguese market. Food Chem. 2018, 239, 132–140. [Google Scholar] [CrossRef]
  8. De Ruyck, K.; Huybrechts, I.; Yang, S.; Arcella, D.; Claeys, L.; Abbeddou, S.; De Keyzer, W.; De Vries, J.; Ocke, M.; Ruprich, J.; et al. Mycotoxin exposure assessments in a multi-center European validation study by 24-hour dietary recall and biological fluid sampling. Environ. Int. 2020, 137, 105539. [Google Scholar] [CrossRef]
  9. Minekus, M.; Alminger, M.; Alvito, P.; Ballance, S.; Bohn, T.; Bourlieu, C.; Carrière, F.; Boutrou, R.; Corredig, M.; Dupont, D.; et al. A standardised static in vitro digestion method suitable for food—An international consensus. Food Funct. 2014, 5, 1113–1124. [Google Scholar] [CrossRef] [Green Version]
  10. Mulet-Cabero, A.-I.; Rigby, N.M.; Brodkorb, A.; Mackie, A.R. Dairy food structures influence the rates of nutrient digestion through different in vitro gastric behaviour. Food Hydrocoll. 2017, 67, 63–73. [Google Scholar] [CrossRef] [Green Version]
  11. Ferrer, M.; Manyes, L.; Mañes, J.; Meca, G. Influence of prebiotics, probiotics and protein ingredients on mycotoxin bioaccessibility. Food Funct. 2015, 6, 987–994. [Google Scholar] [CrossRef]
  12. Sobral, M.M.C.; Cunha, S.C.; Faria, M.A.; Martins, Z.; Ferreira, I.M. Influence of oven and microwave cooking with the addition of herbs on the exposure to multi-mycotoxins from chicken breast muscle. Food Chem. 2019, 276, 274–284. [Google Scholar] [CrossRef]
  13. Tolosa, J.; Font, G.; Mañes, J.; Ferrer, E. Mitigation of enniatins in edible fish tissues by thermal processes and identification of degradation products. Food Chem. Toxicol. 2017, 101, 67–74. [Google Scholar] [CrossRef]
  14. Faria, M.; Melo, A.; Ferreira, I.I. Influence of dietary patterns on contaminants bioaccessibility and intestinal transport by in vitro assays. Food Res. Int. 2020, 137, 109358. [Google Scholar] [CrossRef]
  15. Versantvoort, C.H.; Oomen, A.G.; Van de Kamp, E.; Rompelberg, C.J.; Sips, A.J. Applicability of an in vitro digestion model in assessing the bioaccessibility of mycotoxins from food. Food Chem. Toxicol. 2005, 43, 31–40. [Google Scholar] [CrossRef]
  16. Ren, Z.; Guo, C.; Yu, S.; Zhu, L.; Wang, Y.; Hu, H.; Deng, J. Progress in Mycotoxins Affecting Intestinal Mucosal Barrier Function. Int. J. Mol. Sci. 2019, 20, 2777. [Google Scholar] [CrossRef] [Green Version]
  17. Brown, R.; Priest, E.; Naglik, J.R.; Richardson, J.P. Fungal Toxins and Host Immune Responses. Front. Microbiol. 2021, 12, 643639. [Google Scholar] [CrossRef]
  18. Chen, Y.; Zhou, J.; Wang, L. Role and Mechanism of Gut Microbiota in Human Disease. Front. Cell. Infect. Microbiol. 2021, 11, 625913. [Google Scholar] [CrossRef]
  19. Claus, S.P.; Guillou, H.; Ellero-Simatos, S. The gut microbiota: A major player in the toxicity of environmental pollutants? NPJ Biofilms Microbiomes 2016, 2, 16003. [Google Scholar] [CrossRef]
  20. Sharifi-Rad, J.; Rodrigues, C.; Stojanović-Radić, Z.; Dimitrijević, M.; Aleksić, A.; Neffe-Skocińska, K.; Zielińska, D.; Kołożyn-Krajewska, D.; Salehi, B.; Prabu, S.M.; et al. Probiotics: Versatile Bioactive Components in Promoting Human Health. Medicina 2020, 56, 433. [Google Scholar] [CrossRef]
  21. Guerre, P. Mycotoxin and Gut Microbiota Interactions. Toxins 2020, 12, 769. [Google Scholar] [CrossRef] [PubMed]
  22. Izco, M.; Vettorazzi, A.; de Toro, M.; Sáenz, Y.; Alvarez-Erviti, L. Oral Sub-chronic Ochratoxin A Exposure Induces Gut Microbiota Alterations in Mice. Toxins 2021, 13, 106. [Google Scholar] [CrossRef]
  23. Liew, W.-P.-P.; Mohd-Redzwan, S.; Than, L.T.L. Gut Microbiota Profiling of Aflatoxin B1-Induced Rats Treated with Lactobacillus casei Shirota. Toxins 2019, 11, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Wang, J.; Tang, L.; Glenn, T.C.; Wang, J.S. Aflatoxin B1 Induced Compositional Changes in Gut Microbial Communities of Male F344 Rats. Toxicol. Sci. 2016, 150, 54–63. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. De Angelis, E.; Monaci, L.; Mackie, A.R.; Salt, L.J. Bioaccessibility of T-2 and HT-2 toxins in mycotoxin contaminated bread models submitted to in vitro human digestion. Innov. Food Sci. Emerg. Technol. 2013, 25, 88–96. [Google Scholar] [CrossRef]
  26. González-Arias, C.A.; Marín, S.; Sanchis, V.; Ramos, A.J. Mycotoxin bioaccessibility/absorption assessment using in vitro digestion models: A review. World Mycotoxin J. 2013, 6, 167–184. [Google Scholar] [CrossRef] [Green Version]
  27. Kabak, B.; Brandon, E.F.; Var, I.; Blokland, M.; Sips, A.J. Effects of probiotic bacteria on the bioaccessibility of aflatoxin B1 and ochratoxin A using an in vitro digestion model under fed conditions. J. Environ. Sci. Health Part B 2009, 44, 472–480. [Google Scholar] [CrossRef]
  28. SANTE. Method Validation and Analytical Quality Control Procedures for Pesticide Residues Analysis in Food & Feed; SANTE/12682/2019; European Commission: Brussels, Belgium, 2019. [Google Scholar]
  29. Raters, M.; Matissek, R. Thermal stability of aflatoxin B1 and ochratoxin A. Mycotoxin Res. 2008, 24, 130–134. [Google Scholar] [CrossRef]
  30. Broom, L. Mycotoxins and the intestine. Anim. Nutr. 2015, 1, 262–265. [Google Scholar] [CrossRef]
  31. Carvajal-Moreno, M. Do Gastroenterologists Consider Aflatoxins as Origin of Digestive System Cancers? J. Pharmacovigil. 2017, 5, 242. [Google Scholar] [CrossRef] [Green Version]
  32. Romero, A.; Ares, I.; Ramos, E.; Castellano, V.; Martínez, M.; Martínez-Larrañaga, M.-R.; Anadón, A.; Martínez, M.-A. Mycotoxins modify the barrier function of Caco-2 cells through differential gene expression of specific claudin isoforms: Protective effect of illite mineral clay. Toxicology 2016, 353, 21–33. [Google Scholar] [CrossRef]
  33. Zhang, J.; Zheng, N.; Liu, J.; Li, F.; Li, S.; Wang, J. Aflatoxin B1 and aflatoxin M1 induced cytotoxicity and DNA damage in differentiated and undifferentiated Caco-2 cells. Food Chem. Toxicol. 2015, 83, 54–60. [Google Scholar] [CrossRef]
  34. Liew, W.-P.-P.; Mohd-Redzwan, S. Mycotoxin: Its Impact on Gut Health and Microbiota. Front. Cell. Infect. Microbiol. 2018, 8, 60. [Google Scholar] [CrossRef] [Green Version]
  35. Grosu, I.A.; Pistol, G.C.; Taranu, I.; Marin, D. The impact of dietary grape seed meal on healthy and aflatoxin B1 afflicted microbiota of pigs after weaning. Toxins 2019, 11, 25. [Google Scholar] [CrossRef] [Green Version]
  36. Ishikawa, A.; Weese, J.; Bracarense, A.; Alfieri, A.; Oliveira, G.; Kawamura, O.; Hirooka, E.; Itano, E.; Costa, M. Single aflatoxin B1 exposure induces changes in gut microbiota community in C57Bl/6 mice. World Mycotoxin J. 2017, 10, 249–254. [Google Scholar] [CrossRef]
  37. Liu, T.; Zhang, L.; Joo, D.; Sun, S.C. NF-κB signaling in inflammation. Signal Transduct. Target. Ther. 2017, 2, 17023. [Google Scholar] [CrossRef] [Green Version]
  38. Sugiyama, K.I.; Muroi, M.; Kinoshita, M.; Hamada, O.; Minai, Y.; Sugita-Konishi, Y.; Kamata, Y.; Tanamoto, K.I. NF-κB activation via MyD88-dependent Toll-like receptor signaling is inhibited by trichothecene mycotoxin deoxynivalenol. J. Toxicol. Sci. 2016, 41, 273–279. [Google Scholar] [CrossRef] [Green Version]
  39. Tsai, W.T.; Lo, Y.C.; Wu, M.S.; Li, C.Y.; Kuo, Y.P.; Lai, Y.H.; Tsai, Y.; Chen, K.C.; Chuang, T.H.; Yao, C.H.; et al. Mycotoxin Patulin Suppresses Innate Immune Responses by Mitochondrial Dysfunction and p62/Sequestosome-1-dependent Mitophagy. J. Biol. Chem. 2016, 291, 19299–19311. [Google Scholar] [CrossRef] [Green Version]
  40. Brodkorb, A.; Egger, L.; Alminger, M.; Alvito, P.; Assunção, R.; Ballance, S.; Bohn, T.; Bourlieu-Lacanal, C.; Boutrou, R.; Carrière, F.; et al. INFOGEST static in vitro simulation of gastrointestinal food digestion. Nat. Protoc. 2019, 14, 991–1014. [Google Scholar] [CrossRef]
  41. 2002/657/EC: Commission Decision of 12 August 2002 implementing Council Directive 96/23/EC concerning the performance of analytical methods and the interpretation of results. Off. J. Eur. Communities 2002, 221, 8–36.
  42. ICH. ICH Harmonized Tripartite Guideline: Validation of Analytical Procedures: Text and Methodology Q2(R1). 2005; Available online: (accessed on 27 November 2021). [Google Scholar]
  43. Pinho, S.C.; Faria, M.A.; Melo, A.; Pinto, E.; Almeida, A.; Alves, R.; Cabrita, A.R.; Fonseca, A.J.; Ferreira, I.M. Effect of skimmed milk on intestinal tract: Prevention of increased reactive oxygen species and nitric oxide formation. Int. Dairy J. 2021, 118, 105046. [Google Scholar] [CrossRef]
  44. Kaja, S.; Payne, A.J.; Singh, T.; Ghuman, J.K.; Sieck, E.G.; Koulen, P. An optimized lactate dehydrogenase release assay for screening of drug candidates in neuroscience. J. Pharmacol. Toxicol. Methods 2015, 73, 1–6. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Azad, M.; Kalam, A.; Sarker, M.; Li, T.; Yin, J. Probiotic Species in the Modulation of Gut Microbiota: An Overview. BioMed Res. Int. 2018, 2018, 9478630. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Alcántara, C.; Žugčić, T.; Abdelkebir, R.; García-Pérez, J.V.; Jambrak, A.R.; Lorenzo, J.M.; Collado, M.C.; Granato, D.; Barba, F.J. Effects of Ultrasound-Assisted Extraction and Solvent on the Phenolic Profile, Bacterial Growth, and Anti-Inflammatory/Antioxidant Activities of Mediterranean Olive and Fig Leaves Extracts. Molecules 2020, 25, 1718. [Google Scholar] [CrossRef] [Green Version]
  47. Zwietering, M.H.; Jongenburger, I.; Rombouts, F.M.; Van’T Riet, K. Modeling of the Bacterial Growth Curve. Appl. Environ. Microbiol. 1990, 56, 1875–1881. [Google Scholar] [CrossRef] [Green Version]
  48. Hevia, A.; Delgado, S.; Margolles, A.; Sánchez, B. Application of density gradient for the isolation of the fecal microbial stool component and the potential use thereof. Sci. Rep. 2015, 5, 16807. [Google Scholar] [CrossRef] [Green Version]
  49. Rubio-Del-Campo, A.; Alcántara, C.; Collado, M.C.; Rodríguez-Díaz, J.; Yebra, M.J. Human milk and mucosa-associated disaccharides impact on cultured infant fecal microbiota. Sci. Rep. 2020, 10, 11845. [Google Scholar] [CrossRef]
  50. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Lozupone, C.A.; Turnbaugh, P.J.; Fierer, N.; Knight, R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 2011, 108 (Suppl. S1), 4516. [Google Scholar] [CrossRef] [Green Version]
  51. Cruaud, P.; Vigneron, A.; Lucchetti-Miganeh, C.; Ciron, P.E.; Godfroy, A.; Cambon-Bonavita, M.A. Influence of DNA Extraction Method, 16S rRNA Targeted Hypervariable Regions, and Sample Origin on Microbial Diversity Detected by 454 Pyrosequencing in Marine Chemosynthetic Ecosystems. Appl. Environ. Microbiol. 2014, 80, 4626–4639. [Google Scholar] [CrossRef] [Green Version]
  52. Magoč, T.; Salzberg, S.L. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 2011, 27, 2957–2963. [Google Scholar] [CrossRef]
  53. FLASh: Fast Length Adjustments of Short Reads. Available online: http://ccb.jhu.edu/software/FLASH/ (accessed on 27 December 2021).
  54. Bokulich, N.A.; Subramanian, S.; Faith, J.J.; Gevers, D.; Gordon, J.I.; Knight, R.; Mills, D.A.; Caporaso, J.G. Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 2013, 10, 57–59. [Google Scholar] [CrossRef]
  55. Caporaso, J.G.; Kuczynski, J.; Stombaugh, J.; Bittinger, K.; Bushman, F.D.; Costello, E.K.; Fierer, N.; Peña, A.G.; Goodrich, J.K.; Gordon, J.I.; et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 2010, 7, 335–336. [Google Scholar] [CrossRef] [Green Version]
  56. SILVA High Quality Ribossomal RNA Databases. Available online: http://www.arb-silva.de (accessed on 27 December 2021).
  57. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [Green Version]
Toxins 14 00028 g001 550
Figure 1. Bioaccessibility of AFB1 (a) and OTA (b) during semi-dynamic gastro-duodenal digestion. The graphs represent the cumulative release of AFB1 and OTA in each gastric emptying and respective duodenal digestion when digested isolated or as a mixture. (c) Mycotoxins’ distribution along intestinal emptying (IE1 to IE5) and the content that remained in the non-bioaccessible fraction (NBIO). Data are presented as the mean ± SD of four replicates (n = 4). Parametric (t-test) or non-parametric (Mann–Whitney) statistical tests were applied to compare AFB1 or OTA release, isolated or in combination, and statistical differences are indicated as * p < 0.05, ** p < 0.01, *** p < 0.001, and **** p < 0.0001. E1 to E5 mean emptying 1 to emptying 5.
Figure 1. Bioaccessibility of AFB1 (a) and OTA (b) during semi-dynamic gastro-duodenal digestion. The graphs represent the cumulative release of AFB1 and OTA in each gastric emptying and respective duodenal digestion when digested isolated or as a mixture. (c) Mycotoxins’ distribution along intestinal emptying (IE1 to IE5) and the content that remained in the non-bioaccessible fraction (NBIO). Data are presented as the mean ± SD of four replicates (n = 4). Parametric (t-test) or non-parametric (Mann–Whitney) statistical tests were applied to compare AFB1 or OTA release, isolated or in combination, and statistical differences are indicated as * p < 0.05, ** p < 0.01, *** p < 0.001, and **** p < 0.0001. E1 to E5 mean emptying 1 to emptying 5.
Toxins 14 00028 g001
Toxins 14 00028 g002 550
Figure 2. Impact of bioaccessible mycotoxins on lactate dehydrogenase release (a), nitric oxide formation (b), and reactive oxygen species formation (c) after 3 h exposure of 7-day-differentiated Caco-2 cells. Digests were 12× diluted prior to exposure. Data represent the fold-increase related to cells’ control and are presented as the mean ± SD of two independent experiments. Parametric (Dunnett’s post hoc test) or non-parametric (Kruskal–Wallis and Dunn’s post hoc) tests were applied to compare the effects of AFB1, OTA, or MIX meals with those of the control meal (without mycotoxins). Statistical differences are indicated as * p < 0.05 and **** p < 0.0001.
Figure 2. Impact of bioaccessible mycotoxins on lactate dehydrogenase release (a), nitric oxide formation (b), and reactive oxygen species formation (c) after 3 h exposure of 7-day-differentiated Caco-2 cells. Digests were 12× diluted prior to exposure. Data represent the fold-increase related to cells’ control and are presented as the mean ± SD of two independent experiments. Parametric (Dunnett’s post hoc test) or non-parametric (Kruskal–Wallis and Dunn’s post hoc) tests were applied to compare the effects of AFB1, OTA, or MIX meals with those of the control meal (without mycotoxins). Statistical differences are indicated as * p < 0.05 and **** p < 0.0001.
Toxins 14 00028 g002
Toxins 14 00028 g003 550
Figure 3. Impact of non-bioaccessible mycotoxins on Bifidobacterium lactis and Lactobacillus casei growth over 24 h incubation at anaerobic conditions. Data are presented as the mean ± SEM of four replicates (n = 4). PBS—phosphate buffer solution; NBIO—non-bioaccessible fraction.
Figure 3. Impact of non-bioaccessible mycotoxins on Bifidobacterium lactis and Lactobacillus casei growth over 24 h incubation at anaerobic conditions. Data are presented as the mean ± SEM of four replicates (n = 4). PBS—phosphate buffer solution; NBIO—non-bioaccessible fraction.
Toxins 14 00028 g003
Toxins 14 00028 g004 550
Figure 4. Impact of non-bioaccessible mycotoxins on gut microbiota growth over 24 h incubation under anaerobic conditions: (a) total bacteria, (b) at the phylum level (Firmicutes and Bacteriodetes), and (c) at the family level (Lactobacillaceae, Bifidobacteriaceae, Lachnospiraceae, Ruminococcaceae, and Bacteroidaceae). Data for total bacteria are presented as the mean ± SD of four replicates (n = 4).
Figure 4. Impact of non-bioaccessible mycotoxins on gut microbiota growth over 24 h incubation under anaerobic conditions: (a) total bacteria, (b) at the phylum level (Firmicutes and Bacteriodetes), and (c) at the family level (Lactobacillaceae, Bifidobacteriaceae, Lachnospiraceae, Ruminococcaceae, and Bacteroidaceae). Data for total bacteria are presented as the mean ± SD of four replicates (n = 4).
Toxins 14 00028 g004
Toxins 14 00028 g005 550
Figure 5. Impact of non-bioaccessible mycotoxins on secreted embryonic alkaline phosphatase (SEAP) activity (a) and cell viability (b) after 24 h exposure of colonic HT29 reporter cells monitoring NF-kB activity. Digests were 10× diluted prior to exposure. Data represent the fold-increase related to the control (only cells) and are presented as the mean ± SD of two independent experiments.
Figure 5. Impact of non-bioaccessible mycotoxins on secreted embryonic alkaline phosphatase (SEAP) activity (a) and cell viability (b) after 24 h exposure of colonic HT29 reporter cells monitoring NF-kB activity. Digests were 10× diluted prior to exposure. Data represent the fold-increase related to the control (only cells) and are presented as the mean ± SD of two independent experiments.
Toxins 14 00028 g005
Table
Table 1. Linearity range, regression equations, coefficients of determination (R2), limits of detection (LODs), and limits of quantification (LOQs) of AFB1 and OTA in the 3 analyzed matrices (meal, gastric digest, and intestinal digest).
Table 1. Linearity range, regression equations, coefficients of determination (R2), limits of detection (LODs), and limits of quantification (LOQs) of AFB1 and OTA in the 3 analyzed matrices (meal, gastric digest, and intestinal digest).
LinearityLOD (ng)LOQ (ng)
Range (ng)Regression EquationR2
AFB1Meal2–40y = 0.0428x − 0.1110.99782.27.3
Gastric digest1–20y = 0.0479x + 0.07890.9963
Intestinal digest1–20y = 0.0697x + 0.1450.9939
OTAMeal2–40y = 0.314x − 1.010.99310.762.52
Gastric digest1–20y = 0.655x + 0.1430.9947
Intestinal digest1–20y = 0.328x − 0.1060.9993
Table
Table 2. Recoveries (%), intraday precision and interday precision (% RSD), of the analytical method applied to quantify mycotoxins in the digest fluids (gastric and duodenal) after meal digestion.
Table 2. Recoveries (%), intraday precision and interday precision (% RSD), of the analytical method applied to quantify mycotoxins in the digest fluids (gastric and duodenal) after meal digestion.
Spiked Level (ng/mL)Recovery (Interday, % RSD)Intraday (% RSD)
AFB11792 (9)5
34105 (3)8
85102 (10)9
OTA1186 (25)14
2280 (22)19
5569 (15)28
Table
Table 3. Pearson correlations between AFB1, OTA, and AFB1 + OTA contents in nanograms in the meal before digestion or the bioaccessible fraction and impact on cells (viability and inflammation). Significant correlations are highlighted in bold.
Table 3. Pearson correlations between AFB1, OTA, and AFB1 + OTA contents in nanograms in the meal before digestion or the bioaccessible fraction and impact on cells (viability and inflammation). Significant correlations are highlighted in bold.
AFB1 MealOTA MealAFB1 BIOOTA BIOSum BIO AFB1 + OTANOROSLDH
AFB1 meal1
OTA meal−0.3171
AFB1 BIO0.977−0.2641
OTA BIO−0.4310.946−0.3901
Sum BIO AFB1 + OTA0.3450.7310.4040.6851
NO0.2980.5350.2400.3380.5251
ROS0.179−0.0450.040−0.177−0.1440.6921
LDH0.6400.1870.637−0.0590.4460.7930.5791
Table
Table 4. Effect of non-bioaccessible mycotoxins on the maximal optical density (MOD), growth rate, and lag time of Lactobacillus casei BL23 and Bifidobacterium lactis ATCC 27536 strains. Data are presented as the mean ± SD of four replicates (n = 4). Differences between meals with mycotoxins (AFB1, OTA, or MIX) in relation to the control meal are indicated as * p < 0.05, ** p < 0.01, *** p < 0.001, and **** p < 0.0001.
Table 4. Effect of non-bioaccessible mycotoxins on the maximal optical density (MOD), growth rate, and lag time of Lactobacillus casei BL23 and Bifidobacterium lactis ATCC 27536 strains. Data are presented as the mean ± SD of four replicates (n = 4). Differences between meals with mycotoxins (AFB1, OTA, or MIX) in relation to the control meal are indicated as * p < 0.05, ** p < 0.01, *** p < 0.001, and **** p < 0.0001.
Lactobacillus casei
NBIO AFB1OTAMIXControl
0.5%MOD3.128 ± 0.0363.114 ± 0.0343.176 ± 0.0403.132 ± 0.023
Growth rate
(μmax) (h−1)		0.4502 ± 0.0140.4574 ± 0.0140.453 ± 0.0150.4735 ± 0.010
Lag time (h)5.29 ± 0.1375.228 ± 0.1315.326 ± 0.1495.298 ± 0.088
1%MOD3.23 ± 0.0383.235 ± 0.0243.122 ± 0.052 *3.202 ± 0.032
Growth rate
(μmax) (h−1)		0.3527 ± 0.009 ****0.5211 ± 0.011 ***0.3659 ± 0.014 ****0.4555 ± 0.012
Lag time (h)4.92 ± 0.147 **5.44 ± 0.0844.863 ± 0.210 **5.39 ± 0.121
2%MOD2.364 ± 0.169 ***3.104 ± 0.026 **2.946 ± 0.0612.924 ± 0.059
Growth rate
(μmax) (h−1)		0.131 ± 0.006 ****0.3696 ± 0.007 ****0.228 ± 0.0060.24 ± 0.007
Lag time (h)1.141 ± 1.049 ***5.383 ± 0.103 **4.308 ± 0.2584.398 ± 0.256
Bifidobacterium lactis
NBIO AFB1OTAMIXControl
0.5%MOD7.078 ± 0.025 *6.788 ± 0.041 ***7.013 ± 0.0307.014 ± 0.026
Growth rate
(μmax) (h−1)		0.7646 ± 0.0060.7319 ± 0.0090.7569 ± 0.0070.764 ± 0.006
Lag time (h)5.858 ± 0.0445.808 ± 0.0745.844 ± 0.0535.88 ± 0.046
1%MOD7.078 ± 0.025 *6.788 ± 0.041 ****7.013 ± 0.0307.014 ± 0.026
Growth rate
(μmax) (h−1)		0.7646 ± 0.0060.7319 ± 0.009 **0.7569 ± 0.0070.764 ± 0.006
Lag time (h)5.858 ± 0.0445.808 ± 0.0745.844 ± 0.0535.88 ± 0.046
2%MOD5.585 ± 0.042 ****6.224 ± 0.058 ***6.095 ± 0.0875.978 ± 0.050
Growth rate
(μmax) (h−1)		0.5585 ± 0.008 ****0.6715 ± 0.0130.6742 ± 0.0210.6902 ± 0.013
Lag time (h)5.689 ± 0.092 ****6.281 ± 0.115 *6.289 ± 0.1736.515 ± 0.101
Table
Table 5. Pearson correlations between AFB1, OTA, and AFB1 + OTA contents in nanograms in the meal or the non-bioaccessible fraction and growth curve parameters of L. casei and B. lactis strains. Significant correlations are highlighted in bold.
Table 5. Pearson correlations between AFB1, OTA, and AFB1 + OTA contents in nanograms in the meal or the non-bioaccessible fraction and growth curve parameters of L. casei and B. lactis strains. Significant correlations are highlighted in bold.
AFB1 MealOTA MealAFB1 NBIOOTA NBIOSum NBIO AFB1 + OTAL. caseiB. lactis
MODGrowth RateLag TimeMODGrowth RateLag Time
AFB1 meal1
OTA meal−0.3091
AFB1 NBIO0.954−0.4241
OTA NBIO0.0440.739−0.1811
Sum NBIO AFB1 + OTA0.7720.2560.6290.651
L. caseiMOD−0.7650.641−0.8740.556−0.2361
Growth rate−0.80.786−0.8390.437−0.3040.8711
Lag time−0.7920.63−0.8890.516−0.2790.9980.8781
B. lactisMOD−0.6860.733−0.8170.67−0.1020.9690.8960.9571
Growth rate−0.7530.373−0.8830.431−0.3420.9190.6950.9120.8831
Lag time−0.7720.216−0.8680.271−0.4570.8630.6140.860.7980.9811
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sobral, M.M.C.; Gonçalves, T.; Martins, Z.E.; Bäuerl, C.; Cortés-Macías, E.; Collado, M.C.; Ferreira, I.M.P.L.V.O. Mycotoxin Interactions along the Gastrointestinal Tract: In Vitro Semi-Dynamic Digestion and Static Colonic Fermentation of a Contaminated Meal. Toxins 2022, 14, 28. https://doi.org/10.3390/toxins14010028

AMA Style

Sobral MMC, Gonçalves T, Martins ZE, Bäuerl C, Cortés-Macías E, Collado MC, Ferreira IMPLVO. Mycotoxin Interactions along the Gastrointestinal Tract: In Vitro Semi-Dynamic Digestion and Static Colonic Fermentation of a Contaminated Meal. Toxins. 2022; 14(1):28. https://doi.org/10.3390/toxins14010028

Chicago/Turabian Style

Sobral, Maria Madalena Costa, Tiago Gonçalves, Zita E. Martins, Christine Bäuerl, Erika Cortés-Macías, Maria Carmen Collado, and Isabel M. P. L. V. O. Ferreira. 2022. "Mycotoxin Interactions along the Gastrointestinal Tract: In Vitro Semi-Dynamic Digestion and Static Colonic Fermentation of a Contaminated Meal" Toxins 14, no. 1: 28. https://doi.org/10.3390/toxins14010028

APA Style

Sobral, M. M. C., Gonçalves, T., Martins, Z. E., Bäuerl, C., Cortés-Macías, E., Collado, M. C., & Ferreira, I. M. P. L. V. O. (2022). Mycotoxin Interactions along the Gastrointestinal Tract: In Vitro Semi-Dynamic Digestion and Static Colonic Fermentation of a Contaminated Meal. Toxins, 14(1), 28. https://doi.org/10.3390/toxins14010028

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop