Next Article in Journal
Biological Control and Mitigation of Aflatoxin Contamination in Commodities
Next Article in Special Issue
A 65-Day Fumonisin B Exposure at High Dietary Levels Has Negligible Effects on the Testicular and Spermatological Parameters of Adult Rabbit Bucks
Previous Article in Journal
Ebulin l Is Internalized in Cells by Both Clathrin-Dependent and -Independent Mechanisms and Does Not Require Clathrin or Dynamin for Intoxication
Previous Article in Special Issue
Effects of Single and Repeated Oral Doses of Ochratoxin A on the Lipid Peroxidation and Antioxidant Defense Systems in Mouse Kidneys
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Exposure Assessment of Multiple Mycotoxins and Cumulative Health Risk Assessment: A Biomonitoring-Based Study in the Yangtze River Delta, China

1
College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
2
Institute for Agro-Food Standards and Testing Technology, Shanghai Key Laboratory of Protected Horticultural Technology, Laboratory of Quality and Safety Risk Assessment for Agro-Products (Shanghai), Ministry of Agriculture, Shanghai Academy of Agricultural Sciences, Shanghai 201403, China
*
Authors to whom correspondence should be addressed.
Submission received: 6 January 2021 / Revised: 26 January 2021 / Accepted: 27 January 2021 / Published: 1 February 2021
(This article belongs to the Special Issue Effects of Feedborne Mycotoxins on Animal Health)

Abstract

:
The extensive exposure to multiple mycotoxins has been demonstrated in many countries; however, realistic assessments of the risks related to cumulative exposure are limited. This biomonitoring study was conducted to investigate exposure to 23 mycotoxins/metabolites and their determinants in 227 adults (aged 20–88 years) in the Yangtze River Delta, China. Eight mycotoxins were detected in 110 urine samples, and multiple mycotoxins co-occurred in 51/227 (22.47%) of urine samples, with deoxynivalenol (DON), fumonisin B1 (FB1), and zearalenone (ZEN) being the most frequently occurring. For single mycotoxin risk assessment, FB1, ZEN, aflatoxin B1 (AFB1), and ochratoxin A (OTA) all showed potential adverse effects. However, for the 12 samples containing DON and ZEN, in which none had a hazard risk, the combination of both mycotoxins in two samples was considered to pose potential endocrine disrupting risks to humans by hazard index (HI) method. The combined margin of exposure (MOET) for AFB1 and FB1 could constitute a potential health concern, and AFB1 was the main contributor. Our approach provides a blueprint for evaluating the cumulative risks related to different types of mycotoxins and opens a new horizon for the accurate interpretation of epidemiological health outcomes related to multi-mycotoxin exposure.
Key Contribution: This biomonitoring-based study was performed to assess the cumulative health risks of multiple mycotoxins in the Yangtze River Delta, China.

1. Introduction

Mycotoxins are toxic secondary metabolites produced by various toxic fungal species. More than 400 mycotoxins were identified, including aflatoxins (AFs), fumonisins (FUMs), ochratoxin A (OTA), zearalenone (ZEN), and trichothecenes (TCTs), which exhibited high toxicities and impacted human and animal health (e.g., mutagenic, teratogenic, and even carcinogenic activities) [1,2]. Therefore, the maximum acceptable levels [3] and the health-based guidance values (i.e., tolerable daily intakes (TDI)) have been set by the European Commission and the European Food Safety Authority (EFSA) [4,5,6,7,8].
Humans are exposed to multiple mycotoxins via the environment and daily consumer products. These mycotoxins are quickly metabolized as metabolites/prototypes and are excreted in urine [9]. Recently, urinary metabolites/prototypes have been successfully used as indicators of mycotoxins exposure in many human biomonitoring studies (HBM), which have been performed in Africa, Europe, Asia, and America [10,11,12,13]. In China, most previously conducted investigations were mainly directed toward a single class of mycotoxins, and only one study conducted in 2019 investigated multiple mycotoxins exposure in 260 rural residents (aged 18–66 years) in Nanjing [12]. More data covering larger representative areas are necessary to reveal the real health risks associated with mycotoxins in China.
Given the fact that consumers are exposed to multiple mycotoxins through various routes [14], the cumulative risks related to multiple mycotoxins have attracted more and more attention. Previous studies have proposed the idea that, even if two or more toxins are present in concentrations below those considered safe, exposure to the mixture may result in a cumulative toxic effect [15,16]. However, these bio-monitoring paradigms only considered the adverse effects of one mycotoxin and did not adequately address the overall mycotoxin exposure or assess the realistic risks of combined exposure to multiple mycotoxins. Recently, several models including hazard index (HI), combined margin of exposure (MOET), and relative potency factors/toxic equivalency factors (RPF/TEF) have been established to assess the cumulative risk of co-occurring contaminants based on their chemical composition and knowledge of the toxicities [17,18,19]. Of these, the HI model is defined as the sum of the hazard quotient (HQ) for each individual compound in the mixture with no genotoxic or carcinogenic properties, while the margin of exposure (MOE) is proposed for the risk assessment of genotoxic and carcinogenic substances, and MOET is usually used for the cumulative risk assessment [20].
In the present study, we carried out a study to determine the levels of 23 mycotoxins/metabolites in first morning urine samples in a cohort of 227 adults living in the Jiangsu-Zhejiang-Shanghai area. The relationships between mycotoxins exposure and the dietary habits of residents as well as other elements, e.g., age, sex, and body mass index (BMI), were investigated to better understand the causal factors of the health risks to Yangtze River Delta populations. In order to evaluate the cumulative mycotoxins exposure, HQ, HI, MOE, and MOET models based on the targeted organs and the toxicity effect of typical mycotoxins were proposed to accurately interpret epidemiological health outcomes so as to set up effective prevention and control measures.

2. Results

2.1. Population Characteristics

A total of 294 volunteers, aged 7–88 years, were enrolled in the study. Based on the exclusion criteria, 38 recruited participants were excluded (age < 19 years, incomplete questionnaire, health problems). Overall, 256 healthy subjects returned urine samples and completed questionnaires between March and September 2019. After further validation, 29 urine samples were not eligible due to Cr concentrations <30 mg dL−1 or >300 mg dL−1. Finally, urine samples from 227 adults in the Yangtze River Delta, China were selected and used for the biomonitoring study. A flow chart describing the selection procedures is illustrated in Figure 1. The distributions of participants by sex, age, and BMI in the three sampling regions are presented in Table S1. Of 227 participants recruited for this study, 50.2% were male (n = 114) and 49.8% were female (n = 113). The mean age of all participants was 45.4 ± 17.5 years ranging from 20 to 88 years. The average BMI of the population was 23.1 ± 3.2 kg m−2.

2.2. Prevalence of Mycotoxins in Urine

In total, 227 urine samples were analyzed for the presence of 23 mycotoxins/metabolites, and eight analytes were detected where FB1, ZEN, and DON and its metabolites were the most frequently measured in the urine samples. The prevalence and the levels of these mycotoxin biomarkers are illustrated in Table 1 (uncorrected and corrected by Cr). Only data above the limit of detection (LOD) were included. The quantification of DON-15-GlcA was based on the standard concentration curve of DON-3-GlcA as described in a previous study [21].
DON and its metabolites including DON-3-GlcA and DON-15-GlcA were the most prevalent (prevalence of 12.33%, 10.13%, and 30.40%, respectively), and their concentrations ranged from 0.50 to 35.00 ng mg−1. Their prevalence and concentrations were lower than those reported in populations in China (Henan, Sichuan), UK, Portugal, and Belgium [21,22,23,24]. The low levels of exposure to DON in this study may be due to the minimal occurrence and concentration of DON in cereal and oil products in the Yangtze Delta region [25].
Urinary AFM1 is well established as a biomarker of recent ingestion of AFB1 [26] and, in this study, the prevalence of AFM1 (2.20%) was lower than that of other mycotoxins, confirming critical and efficient control of AFB1 in the Yangtze River Delta. The mean urinary AFM1 concentration in our study was higher than the levels reported in Brazilian subjects (0.02 ng mg−1 Cr) [27] but similar to those reported in adults in Chinese studies (0.278 ng mg−1 Cr, 0.250 ng mg−1) [12] as well as in Nigeria (0.4 ng mg−1) [28] and Malaysia (0.64 ng mg−1) [29]. Considering that AFM1 is a specific biomarker of AFB1 exposure, is a class I carcinogen [30], and is harmful even at lower concentrations, the concentration of aflatoxins should be as low as reasonably achievable [31]. These results confirmed that it is necessary to continuously pay attention to the prevention and the control of AFB1 pollution in the Jiangsu-Zhejiang-Shanghai area.
With regard to OTA, its prevalence (3.52%) in the current study was markedly higher than that found in Nanjing, where OTA was detected in 1.2% of urine samples [12]. Differences in exposure may be attributed to different sample collection times and variable environmental conditions, e.g., high temperate and moist climate [27]. In European countries, such as Germany, Czech Republic, and Belgium, OTA was detected in urine samples at a prevalence ranging from 35% to 100% [11,21,32]. This higher prevalence could be explained by the different dietary habits in European countries, such as coffee and red wine consumption, favorable matrices for fungal growth and OTA production [33,34], which are different to those in China.
FB1, classified as a group 2B human carcinogen by International Agency for Research on Cancer (IARC) [30], was detected in 12.33% of urine samples with a mean value of 1.12 ng mL−1, which was higher than that detected in Portugal (3%, 0.24 ng mL−1) and Haiti (4%, 0.44 ng mL−1) [22,35] but significantly lower than those reported in children in Cameroon (11%, 2.96 ng mL−1) and Nigeria (13.3%, 4.6 ng mL−1) [26,28]. The higher prevalence and concentration in Africa and Asia could be explained by the different dietary habits as the populations in these regions consume a high percentage of their calories as corn meal, and this often coincides with the growing areas where mycotoxin contamination may be high [36]. In addition, children exhibited higher consumption relative to body weight when compared to adults, which is a cause for concern in Cameroon and Nigeria.
The results also showed that 11.89% of urine samples were contaminated by ZEN. ZEN-14-GlcA, the metabolite of ZEN, was detected in two urine samples in this study. The presence of ZEN and its metabolite was significantly lower than that found in previous biomonitoring studies in Bangladesh (ZEN 100%, 0.028 ng mL−1; α-ZEL 100%, 0.198 ng mL−1; β-ZEL 18%, 0.013 ng mL−1) [37], Germany (ZEN 100%, 0.10 ng mL−1; α-ZEL 100%, 0.16 ng mL−1; β-ZEL 18%, 0.05 ng mL−1) [38], and Nigeria (ZEN 81.7%, 0.75 ng mL−1; α-ZEL 4.2%, 1.27 ng mL−1; and β-ZEL 5.8%, 0.88 ng mL−1) [39]. These differences in prevalence rates, average concentration, and concentration range may be the result of the analytical methods employed as well as environmental and food storage conditions.
In addition, one or more analytes were detected in 110/227 (48.46%) urine samples (Table S2), verifying frequent exposures to mycotoxins. Detectable levels of two co-occurrent mycotoxins were found in 30 samples (13.22%). The most common combinations were DON+DON-15-GlcA, DON-3-GlcA+DON-15-GlcA, and DON-15-GlcA+ZEN, which were found in 12 (5.29%), 7 (3.08%), and 4 (1.76%) samples, respectively. Among the 13 samples (5.73%) with three mycotoxins, ZEN and DON and their metabolites were the most frequently detected. Four different mycotoxins were simultaneously found in eight urine samples, demonstrating that humans are commonly exposed to multiple mycotoxins.

2.3. Differences by Residential Area, BMI and Demographic Variables

The exposure levels to some mycotoxins differed due to residential area, sex, age, and BMI. These findings might be attributed, at least in part, to dietary characteristics, body fitness, and levels of mycotoxin contamination in food. The concentrations of mycotoxins in Jiangsu, Zhejiang, and Shanghai were slightly different (Table 1). OTA was less likely to be detected in Shanghai than in Jiangsu and Zhejiang. Jiangsu showed higher levels of urinary FB1, DON, and ZEN than the other two regions. The differences in the mean concentrations of DON and its metabolites in the three regions were significant (p < 0.05) when the urinary concentrations were corrected and uncorrected for Cr. This might be because Jiangsu province is an important wheat growing region located in China’s north-south transition zone and spans two major ecological regions suitable for wheat cultivation [40], which are easily infected by Fusarium species and are vulnerable to contamination by mycotoxins. In comparison, Zhejiang with fertile soils and abundant water resources is one of the world’s most productive rice growing regions [41], while in Shanghai, there is only a small amount of cultivated land, which is also mainly used for rice cultivation, certainly resulting in a lower level of contamination by Fusarium mycotoxins. There are no obvious differences in the mean concentrations of other mycotoxins in the three regions (p > 0.05). In addition, the higher prevalence and concentration of mycotoxins in the current study compared to those reported by Fan et al. [12] may be due to the different environmental conditions and climates between March–September 2019 and January–March 2017.
With regard to age, generally, DON and ZEN were more likely to be detected in young adults (20–45 years) than in the other age groups, and their prevalence decreased with increasing age (Figure 2a). The two main metabolites of DON, DON-3-GlcA, and DON-15-GlcA, were less likely to be detected in older persons (65–88 years) and middle aged adults (46–64 years), respectively, while the metabolite of ZEN, ZEN-14-GlcA was detected in trace amounts in all samples and was not found in middle-aged adults. In general, mycotoxins were more frequently detected in young persons than in middle-aged adults or older persons with significant differences observed (p < 0.05), because younger adults have higher food consumption compared to older adults [42]. However, the differences in the concentrations of mycotoxins were not significant between the different age groups (Table S3). Different from other mycotoxins, the lowest prevalence and concentration of FB1 in urine were observed in 20–45 years group due to the rapid elimination of this toxin resulting in a low bioavailability [43].
The effect of sex was necessarily considered [44], and slight differences in mycotoxin exposure were observed between men and women (Figure 2b and Table S3). Compared to males, females had a higher prevalence of AFM1 and FB1 but a lower prevalence ZEN. The differences were not significant expect FB1 (p > 0.05). The prevalence of OTA, DON, DON-15-GlcA, and ZEN was more likely to be detected in the underweight group (BMI < 18.5 kg m−2) than in the normal weight group (BMI 18.5–23.9 kg m−2) and the overweight group (BMI > 24 kg m−2). FB1 was more likely to be detected in the overweight group than in the other two groups (Figure 2c), which might be due to higher food consumption. However, the underweight group (BMI < 18.5 kg m−2) in our study preferred to drink beverages such as coffee and tea, which are sensitive to OTA contamination. In addition, due to the short half-lives of mycotoxins (DON and ZEN and its metabolites), they were less likely to be found in overweight adults, which could be caused by the relatively higher basal metabolic rate [45].

2.4. Correlation between Mycotoxins and Food Consumption

To identify the possible sources and the determinants of mycotoxins exposure, the correlations between urinary concentrations and food consumption were analyzed based on the information in the food questionnaires. As shown in Table 2, urinary AFM1 concentrations showed a high correlation (rs = 0.866) with the consumption of nuts and seeds, as these foods are frequently contaminated by aflatoxins [46]. Urinary OTA concentrations also showed a high positive correlation (rs = 0.778, p < 0.05) with the consumption of beverages such as coffee and tea, which are easily and frequently contaminated by OTA [33,34]. Only a few moderate correlations were observed between AFM1 levels and the consumption of milk and dairy products (rs = 0.289), wheat (rs = 0.224), and beverages such as coffee and tea (rs = −0.224), which was consistent with previous studies [47]. Milk and dairy consumption was negatively correlated with the OTA levels (rs = −0.275) but positively correlated with FB1 levels (rs = 0.362). There were no statistically significant correlations (p < 0.05) between mycotoxins exposure and food consumption except OTA and beverages such as coffee and tea, which was consistent with several previous studies [48,49] and might be caused by the heterogeneous distributions of mycotoxins in foodstuffs, various routes of mycotoxin exposure, and the inaccurate dietary intake in questionnaires [50].

2.5. Risk Assessment of Single Mycotoxins

The dietary intake of AFB1, FB1, OTA, DON, and ZEN was estimated using the concentrations of biomarkers measured in urine samples (Table 3). For AFB1, a mean probable daily intake (PDI) of 0.62 µg/kg·bw/day and a maximum PDI of 1.03 µg/kg·bw/day were determined when the excretion rate of 1.5% was considered [51], which were higher than those in Brazil (mean intake = 0.001 µg/kg·bw/day) and in Haiti (mean intake = 0.03 µg/kg·bw/day) [27,35] but similar to Nigeria (mean intake = 0.67 µg/kg·bw/day) [28], which posed a potential public concern for residents in the Jiangsu-Zhejiang-Shanghai area. For OTA, the mean PDI value was 0.14 µg/kg·bw/day, which was higher than the report on a Portuguese population [22]. However, Franco et al. [27] reported a greater prevalence of OTA and a lower mean PDI value (0.031 µg/kg·bw/day) due to a high excretion rate (50%). The MOE values of AFB1 and OTA were lower than 200 in a high risk-range, leading to high concern for these mycotoxins.
In adults with positive detection of mycotoxins, 18.93% showed potential health risks with HQ values > 1. FB1 exhibited the highest hazard effect, and 12.33% of participants were exposed to this mycotoxin with HQ values > 1 and emerged with higher adverse effects than in Italy, Portugal, and Brazil [22,27,35]. This could be associated with the high occurrence of FB1 (57%) in maize-based products in China [52]. However, in this study, none of the participants showed a hazard risk in relation to DON, and the mean PDI was 0.08 µg/kg·bw/day, which was considered to be as safe as that in the overall population studied. In addition, higher HQ levels (> 1) for ZEN (0.88% of participants) were less common in adults in the current study than in the Portuguese (24%) [22] and the German (1.67%) [38] populations, which might be due to the various dietary habits as well as toxicokinetics and toxicodynamics of the populations from different countries.

2.6. Cumulative Risk Assessment of Co-Occurring Mycotoxins

As previously referenced, humans are simultaneously exposed to multiple mycotoxins, and our results on the occurrence of various mycotoxins in urine samples confirmed this. However, until now, the risk assessment of mycotoxins was normally performed on each single substance, leading to an underestimation of risk. To accurately evaluate the cumulative risks of multiple mycotoxins, HI and MOET models were proposed in this study, which are primary and frequently used regulatory approaches for the risk assessments of mixtures of toxicologically similar chemicals [20]. Considering the number of mycotoxins present in urine samples and all possible mixtures, the co-occurring mycotoxins were grouped with respect to their similar toxicological effects on a target organ (or system). To accommodate the consistency required for the risk assessment of mixtures, the reference values for mixed mycotoxins should be in the orientation to the same toxicity endpoints, which might be different to the values used for the individual risk assessment [53]. The toxicological data and the key studies/critical effects were obtained from EFSA [4,5,6,7,8] and are shown in Table S4. AFB1 and FB1 displayed qualitatively similar toxicological profiles with hepatotoxicity being the common toxic effect, which could lead to hepatocellular carcinomas or an increased incidence of megalocytic hepatocytes at low dose. With regard to the endocrine activities of DON and ZEN, ZEN had the strongest estrogenic potency, and DON showed estrogenic activity in vitro [54]. DON and ZEN could promote germ cell degeneration, sperm retention, and abnormal nuclear morphology or disturbance of the estrous cycle, ovulation, conception, and implantation in vivo. Cumulative risk assessments were subsequently performed for hepatotoxicity (AFB1 and FB1) and endocrine activity effects (DON and ZEN), respectively.
Two urine samples co-contaminated with AFB1 and FB1 (Group 1) and twelve samples containing both DON and ZEN (Group 2) were selected for cumulative health risk assessment based on the MOET and the HI models, respectively (Table 4). Due to the risk decreasing as the MOE increased [55], in two urine samples with co-occurrence of AFB1 and FB1, AFB1, which is considered the most potent aflatoxin, revealed a MOE below 10,000 in a high risk-range, and FB1 affecting the same target organ would further strengthen the harm. A number of articles also have reported the combined effects of AFB1 and FB1 in experimental animals [56,57]. Although, in the present experiment, there are few samples which were co-contaminated by AFB1 and FB1, their co-existence and combined toxicity still need to be closely monitored due to the fact that they are two of the toxicologically most relevant mycotoxins occurring in maize [56]. Due to the co-occurrence of DON and ZEN in different foodstuffs [25,58], these two mycotoxins were also the most frequently co-occurring mycotoxins in urine samples. When only one mycotoxin was considered, all 12 samples were found to be safe (HQ < 1), but when the combination of both were considered, two samples showed potential endocrine activity risks to humans. Due to the high prevalence and the low TDI value, ZEN was obviously the larger contributor to the total risk compared to DON. The above observations clearly indicated that risk assessments of single mycotoxins were strikingly insufficient and underestimated numerous adverse effects; therefore, the cumulative risks of different mycotoxins should be an important research core, which would attract more attention and further studies.
The main limitations of this study also needed to be underlined. First, dietary intake, weight, socio-demographic, and lifestyle information were self-reported; therefore, misreporting could not be excluded and may be related to other factors, including socioeconomic status, occupation, perceived weight status, and level of physical activity. Second, due to the short half-lives of some mycotoxins (such as DON, DON-3-GlcA, and α-ZEL) [22,59] included in this study, one spot urine is unlikely to monitor every exposure to mycotoxins, which could also lead to an underestimation of mycotoxins exposure. Thirdly, the pending elucidation of the mode of actions for the genotoxicity/carcinogenicity of OTA and FB1, an MOE of 10,000 to the neoplastic reference point, is warranted for the risk characterization of OTA and calculating the MOET with AFB1 and FB1. The assessment is more likely to overestimate than to underestimate the risk [8]. Further investigations are necessary to confirm the findings of this study.

3. Conclusions

This biomonitoring-based study was performed to assess the cumulative health risks of multiple mycotoxins in the Yangtze River Delta, China. The ubiquitous exposure of a broad segment of the investigated population to mycotoxins was observed, and 22.47% of adults were found to be co-exposed to two or more contaminants, mainly ZEN, DON, and their metabolites. The heterogeneity of mycotoxins exposure was identified, and several socio-demographic lifestyle characteristics, such as place of residence, age, sex, and BMI, were associated with mycotoxin body burden. Considering the high prevalence of co-contamination with multiple mycotoxins, MOE or HQ and HI or MOET models were put forward for risk assessment of exposure to single and multiple mycotoxins. Potential hazard effects were shown to be related to single mycotoxins, and furthermore, the integration of two mycotoxins clearly contributed to a greater health concern, which endowed the risks related to the mycotoxins within the TDI limits. For the 12 samples containing DON and ZEN, in which none had a hazard risk, the combination of both mycotoxins in two samples was considered to pose potential endocrine disrupting risks to humans. This is the first time an efficient strategy has been proposed for the quantitative assessment of cumulative health risks of multiple mycotoxins exposure using urinary biomarkers. This strategy successfully avoided underestimation of the health risks when exposed to multiple mycotoxins. Thus, an overhaul of legislation and control programs is required, as not only the risks of individual mycotoxins should be considered but also their synergic effects.

4. Materials and Methods

4.1. Reagents and Chemicals

Methanol and acetonitrile, both HPLC grade, were obtained from Merck (Darmstadt, Germany). Formic acid and ammonium acetate were obtained from Sigma-Aldrich (St. Louis, MO, USA). Water was purified using a Milli-Q Plus apparatus (18.2 Ω, Millipore, Billerica, MA, USA) prior to use. The standards of aflatoxin B1 (AFB1), aflatoxin B2 (AFB2), aflatoxin G1 (AFG1), aflatoxin G2 (AFG2), aflatoxin M1 (AFM1), aflatoxin M2 (AFM2), OTA, ochratoxin α (OTα), fumonisin B1 (FB1), T-2 toxin (T-2), HT-2 toxin (HT-2), deoxynivalenol (DON), 3-acetyldeoxynivalenol (3-ADON), 15-acetyldeoxynivalenol (15-ADON), fusarenon X (FUS-X), ZEN, α-zearalenol (α-ZEL), β-zearalenol (β-ZEL), α-zearalanol (α-ZAL), and β-zearalanol (β-ZAL) were purchased from Romer Labs (Union, MO, USA). The mycotoxin conjugates deoxynivalenol-3-glucuronide (DON-3-GlcA) and zearalenone-14-glucuronide (ZEN-14-GlcA) were prepared in our laboratory as described in a previous study [12]. The quantification of DON-15-GlcA was based on the standard concentration curve of DON-3-GlcA as described in a previous study [21].

4.2. Study Population

Morning urine samples were obtained from 227 adults residing in the Yangtze River Delta, China between March 2019 and September 2019. Attempts were made to achieve a representative distribution for sex (male and female), age (20–45, 46–64, 65–88-years-old) [60], and region (Jiangsu, Zhejiang, and Shanghai). The participants were asked to complete a detailed questionnaire, which included demographics, health, and lifestyle factors. A food questionnaire on dietary intake the previous day was administered, in which typical foods such as wheat, maize, rice, vegetables and fruit, meat, nuts and seeds, milk and dairy products, and drink (including beverages such as coffee and tea) consumed by Jiangsu, Zhejiang, and Shanghai residents were clearly listed. The results for daily food intake in the participants are shown in Supplementary Table S5.

4.3. Urine Sample Preparation

The morning urine samples (~20 mL) were collected from the participants using 50 mL sterile centrifuge tubes (Corning, NY, USA). All urine samples were transported to the laboratory in chilled conditions after collection and stored in a freezer (−20 °C) within 6 h.

4.4. Biomarker Analysis

The extraction and the purification of 23 mycotoxins/metabolites in urine samples were performed as described in our previous study [12]. Briefly, each urine sample (200 μL) was measured accurately into a 1.5 mL centrifuge tube, and 1.0 mL of acetonitrile containing 1% formic acid was added. The sample was vortexed for 2 min and centrifuged at 12,000× g for 10 min at 4 °C. The acquired supernatant was dried by nitrogen stream, and then the residues were re-dissolved in 200 μL of acetonitrile/water containing 5 mmol L−1 ammonium acetate (20/80, v/v) and passed through a 0.22 μm syringe filters before ultra-high-performance-liquid chromatography tandem mass spectrometry (UHPLC-MS/MS) analysis. The instrumental conditions are detailed in Supplementary S1, and observed mass transitions and target-specific parameters are listed in Table S6.
To account for the variable dilutions of urine, urinary mycotoxin concentration was also reported as the adjusted concentration on the basis of creatinine (ng mg−1 creatinine (Cr)), the concentration of which was measured using an enzymatic reaction on a Roche Hitachi 912 Chemistry Analyzer (Roche Hitachi, Basel, Switzerland). According to the proposal of the World Health Organization (WHO) [61], samples too diluted (Cr. concentration < 30 mg dL−1) or too concentrated (Cr. concentration > 300 mg dL−1) were excluded.

4.5. Estimated Dietary Exposure of Mycotoxins

Probable daily intake (PDI) of mycotoxin was calculated according to Equation (1) [62]:
PDI   ( µ g / kg · bw / day )   =   C × V × 100 W × E
where C is urinary biomarker concentration (μg L−1) (for DON, C(DON) = C (free DON) + (C (DON-3-GlcA) + C (DON-15-GlcA))*0.63), V is the mean volume of daily urine production (1.5 L day−1), W is body weight (kg), and E is mean urinary excretion rate (in percent) (1.5 for AFM1 [50], 2.5 for OTA [63], 0.5 for FB1 [64], 9.4 for ZEN [14], and 72 for DON [65]).

4.6. Risk Characterization

4.6.1. Risk Assessment of Single Mycotoxins

Two different approaches were performed for the risk assessment of individual mycotoxins based on their carcinogenicity. The MOE was applied for aflatoxins and OTA and was calculated as a ratio of benchmark dose lower confidence limit (BMDL10) and their exposure (as Equation (2)). Considering the carcinogenicity of OTA and AFB1, no level of exposure is considered safe, and the MOE was just calculated to help prioritize risk management actions. BMDL10 value were 0.4 µg/kg·bw/day for aflatoxins and 14.5 µg/kg·bw/day for OTA with neoplastic effects [7,8]. The magnitude of the MOE gives an indication of the risk level, and the Scientific Committee of EFSA and WHO have concluded that a MOE of 10,000 or more is of low concern for public health [20].
MOE = BMDL10/PDI
For the remaining mycotoxins, HQ, used for assessing the adverse effects related to single mycotoxins, is the ratio of the PDI to the reference dose for relevant health endpoints [66], which could be the TDI (Equation (3)) or the value of the no observed adverse effect level (NOAEL) (Equation (4)). The uncertainty factor (UF) was set as 100 (10 × 10), which was interpreted as reflecting extrapolation from experimental animals to human (factor 10 for inter-species variability between man and animal) and extrapolation from an average human NOAEL to a sensitive human NOAEL (factor 10 for human or intra-species variability) [67]. When the HQ was <1, the exposure was considered to be within safe limits [20].
HQTDI = PDI/TDI
HQNOVEL = PDI/NOAEL × UF

4.6.2. Cumulative Risk Assessment of Co-Occurring Mycotoxins

Considering the co-occurrence of different mycotoxins, two groups were established with regard to toxicological similarity: Group 1, AFB1 and FB1 based on hepatotoxicity; Group 2, DON and ZEN based on endocrine activity (testicular toxicity and strong estrogenic potency). The MOET was utilized for Group 1 and the HI method for Group 2. The MOET was calculated as the reciprocal of the sum of the reciprocals of the individual MOEs (Equation (5)), and it was considered to be high concern for public health if MOET < 100. HI was the sum of the respective HQs for the individual mixture components of the same toxicity (Equation (6)). Significant adverse effects were speculated if HI > 1. The values of BMDL10, NOAELs, and TDIs used in the current study are presented in Table S4.
MOET = 1/(1/MOEAFB1 + 1/MOEFB1)
HI = HQTDI,ZEN + HQNOVEL,DON

4.7. Statistical Analysis

All statistical analyses were performed using IBM SPSS Statistics software version 25 (SPSS Inc., Chicago, IL, USA) and Origin version 8.0 (Microcal, Inc., Northampton, MA, USA). For data analysis, urinary mycotoxin concentrations below the limit of detection (LOD) were replaced by zero, and the mycotoxins concentration values below the limit of quantification (LOQ) but above the LOD were replaced by 1/2 LOQ. As the data were not normally distributed, Fisher’s exact test and rank sum test were used to investigate the distributions of urinary mycotoxins by region, sex, age, and BMI. Spearman’s rank correlation coefficient (rs, ≤ 0.2 = poor; 0.2 ≤ 0.5 = fair; 0.5 ≤ 0.7 = moderate; 0.7 ≤ 0.9 = very strong) [22] was used for analysis of the relationships between two subgroups. When the p-value was <0.05, the differences were considered statistically significant.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/13/2/103/s1, S1: Detailed UHPLC-MS/MS method, Table S1: The demographic characteristics of participants in the three sampling areas, Table S2: Co-occurrence of investigated mycotoxins and their metabolites in urine samples from Jiangsu, Zhejiang, and Shanghai., Table S3: P-value of mycotoxin concentration distribution by age, sex, and BMI using the rank sum test (n = 227),Table S4: Toxicities of FB1, OTA, DON, and ZEN from EFSA data, Table S5: Food consumption in the total population stratified by location, Table S6: LC-MS/MS parameters for the detection of targeted mycotoxins.

Author Contributions

Conceptualization, Z.H. and Y.W.; methodology, Q.H., K.J., Z.T., D.N. and Z.Z.; validation, Q.H., K.J., Z.T. and K.F.; formal analysis, Q.H., K.J., Z.T. and K.F.; investigation, Q.H., K.J., K.F., J.M., D.N. and Z.Z.; resources, J.M., Z.H. and Y.W.; writing—original draft preparation, Q.H.; writing—review and editing, Z.H. and Y.W.; visualization, Q.H., J.M., D.N. and Z.Z.; supervision, Z.H. and Y.W.; project administration, Z.H. and Y.W.; funding acquisition, Z.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Shanghai Agriculture Applied Technology Development Program, China (Grant No. 2019-02-08-00-12-F01148).

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and European Commission recommendations.

Informed Consent Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Goud, K.Y.; Kailasa, S.K.; Kumar, V.; Tsang, Y.F.; Lee, S.E.; Gobi, K.V.; Kim, K. Progress on nanostructured electrochemical sensors and their recognition elements for detection of mycotoxins: A review. Biosens. Bioelectron. 2018, 121, 205–222. [Google Scholar] [CrossRef]
  2. Yang, Y.; Li, G.; Wu, D.; Liu, J.; Li, X.; Luo, P.; Hu, N.; Wang, H.; Wu, Y. Recent advances on toxicity and determination methods of mycotoxins in foodstuffs. Trends Food Sci. Technol. 2020, 96, 233–252. [Google Scholar] [CrossRef]
  3. Commission Regulation 1881/2006. Commission Regulation (EC) No 1881/2006 of 19 December 2006 setting maximum levels for certain contaminants in foodstuffs. Off. J. Eur. Commun. 2006, L364, 5–24. Available online: https://data.europa.eu/eli/reg/2006/1881/oj (accessed on 20 May 2019).
  4. European Food Safety Authority (EFSA). Scientific opinion on the risks for public health related to the presence of zearalenone in food. EFSA J. 2011, 9, 2197. [Google Scholar] [CrossRef]
  5. European Food Safety Authority (EFSA). Risks to human and animal health related to the presence of deoxynivalenol and its acetylated and modified forms in food and feed. EFSA J. 2017, 15, 4718. [Google Scholar] [CrossRef]
  6. European Food Safety Authority (EFSA). Appropriateness to set a group health-based guidance value for fumonisins and their modified forms. EFSA J. 2018, 16, 5172. [Google Scholar] [CrossRef]
  7. European Food Safety Authority (EFSA). Risk assessment of aflatoxins in food. EFSA J. 2020, 18, 6040. [Google Scholar] [CrossRef]
  8. European Food Safety Authority (EFSA). Risk assessment of ochratoxin A in food. EFSA J. 2020, 18, 6113. [Google Scholar] [CrossRef]
  9. Slobodchikova, I.; Sivakumar, R.; Rahman, M.S.; Vuckovic, D. Characterization of Phase I and glucuronide Phase II metabolites of 17 mycotoxins using liquid chromatography-high-resolution mass spectrometry. Toxins 2019, 11, 433. [Google Scholar] [CrossRef] [Green Version]
  10. Chen, C.; Mitchell, N.J.; Gratz, J.; Houpt, E.R.; Gong, Y.; Egner, P.A.; Groopman, J.D.; Riley, R.T.; Showker, J.L.; Svensen, E.; et al. Exposure to aflatoxin and fumonisin in children at risk for growth impairment in rural Tanzania. Environ. Int. 2018, 115, 29–37. [Google Scholar] [CrossRef]
  11. 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]
  12. Fan, K.; Xu, J.; Jiang, K.; Liu, X.; Meng, J.; Di Mavungu, J.D.; Guo, W.; Zhang, Z.; Jing, J.; Li, H.; et al. Determination of multiple mycotoxins in paired plasma and urine samples to assess human exposure in Nanjing, China. Environ. Pollut. 2019, 248, 865–873. [Google Scholar] [CrossRef]
  13. Jager, A.V.; Tonin, F.G.; Baptista, G.Z.; Souto, P.C.; Oliveira, C.A. Assessment of aflatoxin exposure using serum and urinary biomarkers in Sao Paulo, Brazil: A pilot study. Int. J. Hyg. Environ. Health 2016, 219, 294–300. [Google Scholar] [CrossRef]
  14. Warth, B.; Sulyok, M.; Berthiller, F.; Schuhmacher, R.; Krska, R. New insights into the human metabolism of the Fusarium mycotoxins deoxynivalenol and zearalenone. Toxicol. Lett. 2013, 220, 88–94. [Google Scholar] [CrossRef] [Green Version]
  15. Alassane-Kpembi, I.; Schatzmayr, G.; Taranu, I.; Marin, D.; Puel, O.; Oswald, I.P. Mycotoxins co-contamination: Methodological aspects and biological relevance of combined toxicity studies. Crit. Rev. Food Sci. Nutr. 2016, 57, 3489–3507. [Google Scholar] [CrossRef]
  16. Wielogorska, E.; Mooney, M.; Eskola, M.; Ezekiel, C.N.; Stranska, M.; Krska, R.; Elliott, C. Occurrence and human health impacts of mycotoxins in Somalia. J. Agric. Food Chem. 2019, 67, 2052–2060. [Google Scholar] [CrossRef]
  17. De Gavelle, E.; De Lauzon-Guillain, B.; Charles, M.; Chevrier, C.; Hulin, M.; Sirot, V.; Merlo, M.; Nougadère, A. Chronic dietary exposure to pesticide residues and associated risk in the French ELFE cohort of pregnant women. Environ. Int. 2016, 92–93, 533–542. [Google Scholar] [CrossRef]
  18. Garí, M.; Koch, H.M.; Pälmke, C.; Jankowska, A.; Wesołowska, E.; Hanke, W.; Nowak, D.; Bose-O’Reilly, S.; Polańska, K. Determinants of phthalate exposure and risk assessment in children from Poland. Environ. Int. 2019, 127, 742–753. [Google Scholar] [CrossRef]
  19. Fallahzadeh, R.A.; Khosravi, R.; Dehdashti, B.; Ghahramani, E.; Omidi, F.; Adli, A.; Miri, M. Spatial distribution variation and probabilistic risk assessment of exposure to chromium in ground water supplies; a case study in the east of Iran. Food Chem. Toxicol. 2018, 115, 260–266. [Google Scholar] [CrossRef]
  20. European Food Safety Authority (EFSA). International frameworks dealing with human risk assessment of combined exposure to multiple chemicals. EFSA J. 2013, 11, 1–69. [Google Scholar] [CrossRef]
  21. Heyndrickx, E.; Sioen, I.; Huybrechts, B.; Callebaut, A.; De Henauw, S.; De Saeger, S. Human biomonitoring of multiple mycotoxins in the Belgian population: Results of the BIOMYCO study. Environ. Int. 2015, 84, 82–89. [Google Scholar] [CrossRef]
  22. Martins, C.; Vidal, A.; De Boevre, M.; De Saeger, S.; Nunes, C.; Torres, D.; Goios, A.; Lopes, C.; Assunção, R.; Alvito, P. Exposure assessment of Portuguese population to multiple mycotoxins: The human biomonitoring approach. Int. J. Hyg. Environ. Health 2019, 222, 913–925. [Google Scholar] [CrossRef]
  23. Turner, P.C.; Rothwell, J.A.; White, K.L.; Gong, Y.; Cade, J.E.; Wild, C.P. Urinary deoxynivalenol is correlated with cereal intake in individuals from the United Kingdom. Environ. Health Perspect. 2008, 116, 21–25. [Google Scholar] [CrossRef]
  24. Wang, X.; Liang, J.; Cao, P.; Zhou, S.; Wu, A.; Gao, P.; Xu, H.; Liu, Z.; Gong, Y. Biomonitoring study of deoxynivalenol exposure in Chinese inhabitants. Int. J. Environ. Res. Public Health 2019, 16, 2169. [Google Scholar] [CrossRef] [Green Version]
  25. Li, R.; Wang, X.; Zhou, T.; Yang, D.; Wang, Q.; Zhou, Y. Occurrence of four mycotoxins in cereal and oil products in Yangtze Delta region of China and their food safety risks. Food Control 2014, 35, 117–122. [Google Scholar] [CrossRef]
  26. Njumbe Ediage, E.; Diana Di Mavungu, J.; Song, S.; Sioen, I.; De Saeger, S. Multimycotoxin analysis in urines to assess infant exposure: A case study in Cameroon. Environ. Int. 2013, 57–58, 50–59. [Google Scholar] [CrossRef] [Green Version]
  27. Franco, L.T.; Petta, T.; Rottinghaus, G.E.; Bordin, K.; Gomes, G.A.; Alvito, P.; Assunção, R.; Oliveira, C.A.F. Assessment of mycotoxin exposure and risk characterization using occurrence data in foods and urinary biomarkers in Brazil. Food Chem. Toxicol. 2019, 128, 21–34. [Google Scholar] [CrossRef]
  28. Ezekiel, C.N.; Warth, B.; Ogara, I.M.; Abia, W.A.; Ezekiel, V.C.; Atehnkeng, J.; Sulyok, M.; Turner, P.C.; Tayo, G.O.; Krska, R.; et al. Mycotoxin exposure in rural residents in northern Nigeria: A pilot study using multi-urinary biomarkers. Environ. Int. 2014, 66, 138–145. [Google Scholar] [CrossRef]
  29. Siti Husna, S.R.J.M. Association between urinary aflatoxin (AFM1) and dietary intake among adults in Hulu Langat District, Selangor, Malaysia. Nutrients 2018, 10, 460. [Google Scholar] [CrossRef] [Green Version]
  30. International Agency for Research on Cancer (IARC). Traditional herbal medicines, some mycotoxins, napthalene, and styrene. IARC Monogr. Eval. Carcinog. Risk Hum. 2002, 82, 1–556. [Google Scholar]
  31. Škrbić, B.; Živančev, J.; Antić, I.; Godula, M. Levels of aflatoxin M1 in different types of milk collected in Serbia: Assessment of human and animal exposure. Food Control 2014, 40, 113–119. [Google Scholar] [CrossRef]
  32. Ali, N.; Muñoz, K.; Degen, G.H. Ochratoxin A and its metabolites in urines of German adults—An assessment of variables in biomarker analysis. Toxicol. Lett. 2017, 275, 19–26. [Google Scholar] [CrossRef] [PubMed]
  33. Bellver Soto, J.; Fernández-Franzón, M.; Ruiz, M.; Juan-García, A. Presence of ochratoxin A (OTA) mycotoxin in alcoholic drinks from Southern European countries: Wine and beer. J. Agric. Food Chem. 2014, 62, 7643–7651. [Google Scholar] [CrossRef]
  34. Casal, S.; Vieira, T.; Cruz, R.; Cunha, S.C. Ochratoxin A in commercial soluble coffee and coffee substitutes. Food Res. Int. 2014, 61, 56–60. [Google Scholar] [CrossRef]
  35. Gerding, J.; Ali, N.; Schwartzbord, J. A comparative study of the human urinary mycotoxin excretion patterns in Bangladesh, Germany, and Haiti using a rapid and sensitive LC-MS/MS approach. Mycotoxin Res. 2015, 31, 127–136. [Google Scholar] [CrossRef] [PubMed]
  36. Soriano, J.M.; Dragacci, S. Occurrence of fumonisins in foods. Food Res. Int. 2004, 37, 985–1000. [Google Scholar] [CrossRef]
  37. Ali, N.; Degen, G.H. Biomonitoring of zearalenone and its main metabolites in urines of Bangladeshi adults. Food Chem. Toxicol. 2019, 130, 276–283. [Google Scholar] [CrossRef]
  38. Ali, N.; Degen, G.H. Urinary biomarkers of exposure to the mycoestrogen zearalenone and its modified forms in German adults. Arch. Toxicol. 2018, 92, 2691–2700. [Google Scholar] [CrossRef]
  39. Šarkanj, B.; Ezekiel, C.N.; Turner, P.C.; Abia, W.A.; Rychlik, M.; Krska, R.; Sulyok, M.; Warth, B. Ultra-sensitive, stable isotope assisted quantification of multiple urinary mycotoxin exposure biomarkers. Anal. Chim. Acta 2018, 1019, 84–92. [Google Scholar] [CrossRef]
  40. Ji, F.; Xu, J.; Liu, X.; Yin, X.; Shi, J. Natural occurrence of deoxynivalenol and zearalenone in wheat from Jiangsu province, China. Food Chem. 2014, 157, 393–397. [Google Scholar] [CrossRef]
  41. Van den Berg, M.M.; Hengsdijk, H.; Wolf, J.; Van Ittersum, M.K.; Guanghuo, W.; Roetter, R.P. The impact of increasing farm size and mechanization on rural income and rice production in Zhejiang province, China. Agr. Syst. 2007, 94, 841–850. [Google Scholar] [CrossRef]
  42. Giezenaar, C.; Chapman, I.; Luscombe-Marsh, N.; Feinle-Bisset, C.; Horowitz, M.; Soenen, S. Ageing is associated with decreases in appetite and energy intake: A meta-analysis in healthy adults. Nutrients 2016, 8, 28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Silva, L.J.; Pena, A.; Lino, C.M.; Fernandez, M.F.; Manes, J. Fumonisins determination in urine by LC-MS-MS. Anal. Bioanal. Chem. 2010, 396, 809–816. [Google Scholar] [CrossRef] [PubMed]
  44. Schiebinger, L. Scientific research must take gender into account. Nature 2014, 507, 9. [Google Scholar] [CrossRef]
  45. Wahrlich, V.; Teixeira, T.M.; Anjos, L.A.D. Validity of a population-specific BMR predictive equation for adults from an urban tropical setting. Clin. Nutr. 2018, 37, 208–213. [Google Scholar] [CrossRef]
  46. Asghar, M.A.; Ahmed, A.; Zahir, E.; Asghar, M.A.; Iqbal, J.; Walker, G. Incidence of aflatoxins contamination in dry fruits and edible nuts collected from Pakistan. Food Control 2017, 78, 169–175. [Google Scholar] [CrossRef]
  47. Montagna, M.; Napoli, C.; De Giglio, O.; Iatta, R.; Barbuti, G. Occurrence of Aflatoxin M1 in Dairy Products in Southern Italy. Int. J. Mol. Sci. 2008, 9, 2614–2621. [Google Scholar] [CrossRef] [Green Version]
  48. Di Giuseppe, R.; Bertuzzi, T.; Rossi, F.; Rastelli, S.; Mulazzi, A.; Capraro, J.; de Curtis, A.; Iacoviello, L.; Pietri, A. Plasma ochratoxin A levels, food consumption, and risk biomarkers of a representative sample of men and women from the Molise region in Italy. Eur. J. Nutr. 2012, 51, 851–860. [Google Scholar] [CrossRef]
  49. Warensjö, L.E.; Montano, M.A.; Schmidt, J.; Cramer, B.; Humpf, H.; Moraeus, L.; Olsen, M. Mycotoxins in blood and urine of Swedish adolescents—Possible associations to food intake and other background characteristics. Mycotoxin Res. 2020, 36, 193–206. [Google Scholar] [CrossRef] [Green Version]
  50. Fan, K.; Cheng, X.; Guo, W.; Liu, X.; Zhang, Z.; Yao, Q.; Nie, D.; Yao, B.; Han, Z. Ochratoxin A in human blood plasma samples from apparently healthy volunteers in Nanjing, China. Mycotoxin Res. 2020. [Google Scholar] [CrossRef]
  51. Zhu, J.Q.; Zhang, L.S.; Hu, X.; Xiao, Y.; Chen, J.S.; Xu, Y.C.; Fremy, J.; Chu, F.S. Correlation of dietary aflatoxin B1 levels with excretion of aflatoxin M1 in human urine. Cancer Res. 1987, 47, 1848. [Google Scholar] [PubMed]
  52. Yang, X.; Yang, X.; Gao, J.; Liu, Q.; Yang, D. Co-occurrence of mycotoxins in maize and maize-derived food in China and estimation of dietary intake. Food Addit. Contam. B 2019, 12, 124–134. [Google Scholar] [CrossRef] [PubMed]
  53. Apel, P.; Kortenkamp, A.; Koch, H.M.; Vogel, N.; Rüther, M.; Kasper-Sonnenberg, M.; Conrad, A.; Brüning, T.; Kolossa-Gehring, M. Time course of phthalate cumulative risks to male developmental health over a 27-year period: Biomonitoring samples of the German Environmental Specimen Bank. Environ. Int. 2020, 137, 105467. [Google Scholar] [CrossRef] [PubMed]
  54. Demaegdt, H.; Daminet, B.; Evrard, A.; Scippo, M.; Muller, M.; Pussemier, L.; Callebaut, A.; Vandermeiren, K. Endocrine activity of mycotoxins and mycotoxin mixtures. Food Chem. Toxicol. 2016, 96, 107–116. [Google Scholar] [CrossRef] [PubMed]
  55. Wang, J.; Zhao, X.; Wang, Y.; Shi, Z. Tetrabromobisphenol A, hexabromocyclododecane isomers and polybrominated diphenyl ethers in foodstuffs from Beijing, China: Contamination levels, dietary exposure and risk assessment. Sci. Total Environ. 2019, 666, 812–820. [Google Scholar] [CrossRef]
  56. Hove, M.; Van Poucke, C.; Njumbe-Ediage, E.; Nyanga, L.K.; De Saeger, S. Review on the natural co-occurrence of AFB1 and FB1 in maize and the combined toxicity of AFB1 and FB1. Food Control 2016, 59, 675–682. [Google Scholar] [CrossRef]
  57. Theumer, M.G.; Cánepa, M.C.; López, A.G.; Mary, V.S.; Dambolena, J.S.; Rubinstein, H.R. Subchronic mycotoxicoses in Wistar rats: Assessment of the in vivo and in vitro genotoxicity induced by fumonisins and aflatoxin B1, and oxidative stress biomarkers status. Toxicology 2010, 268, 104–110. [Google Scholar] [CrossRef]
  58. Da, S.L.; Madureira, F.; de Azevedo, V.E.; Faria, A.F.; Augusti, R. Development and validation of a multianalyte method for quantification of mycotoxins and pesticides in rice using a simple dilute and shoot procedure and UHPLC-MS/MS. Food Chem. 2019, 270, 420–427. [Google Scholar] [CrossRef]
  59. Vidal, A.; Claeys, L.; Mengelers, M.; Vanhoorne, V.; Vervaet, C.; Huybrechts, B.; De Saeger, S.; De Boevre, M. Humans significantly metabolize and excrete the mycotoxin deoxynivalenol and its modified form deoxynivalenol-3-glucoside within 24 hours. Sci. Rep. 2018, 8, 8. [Google Scholar] [CrossRef] [Green Version]
  60. Wang, H.; Yang, J.; Yu, X.; Zhao, G.; Zhao, Q.; Wang, N.; Jiang, Y.; Jiang, F.; He, G.; Chen, Y.; et al. Exposure of adults to antibiotics in a Shanghai suburban area and health risk assessment: A biomonitoring-based study. Environ. Sci. Technol. 2018, 52, 13942–13950. [Google Scholar] [CrossRef]
  61. World Health Organization (WHO). Biological Monitoring of Chemical Exposure in the Workplace; World Health Organization: Geneva, Switzerland, 1996. [Google Scholar]
  62. Gerding, J.; Cramer, B.; Humpf, H.U. Determination of mycotoxin exposure in Germany using an LC-MS/MS multibiomarker approach. Mol. Nutr. Food Res. 2014, 58, 2358–2368. [Google Scholar] [CrossRef] [PubMed]
  63. Degen, G.H. Are we ready to estimate daily ochratoxin A intake based on urinary concentrations? Environ. Int. 2016, 97, 254–255. [Google Scholar] [CrossRef] [PubMed]
  64. Riley, R.T.; Torres, O.; Showker, J.L.; Zitomer, N.C.; Matute, J.; Voss, K.A.; Gelineau-van Waes, J.; Maddox, J.R.; Gregory, S.G.; Ashley-Koch, A.E. The kinetics of urinary fumonisin B1 excretion in humans consuming maize-based diets. Mol. Nutr. Food Res. 2012, 56, 1445–1455. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Turner, P.C.; White, K.L.M.; Burley, V.J.; Hopton, R.P.; Rajendram, A.; Fisher, J.; Cade, J.E.; Wild, C.P. A comparison of deoxynivalenol intake and urinary deoxynivalenol in UK adults. Biomarkers 2010, 15, 553–562. [Google Scholar] [CrossRef] [PubMed]
  66. Teuschler, L.K.; Hertzberg, R.C. Current and future risk assessment guidelines, policy and methods development for chemical mixtures. Toxicology 1995, 105, 137–144. [Google Scholar] [CrossRef]
  67. European Food Safety Authority (EFSA). Guidance on Default assumptions used by the EFSA Scientific Panels and Committee, and EFSA Units in the absence of actual measured data. EFSA J. 2011, 10, 2579. [Google Scholar]
Figure 1. Flow chart describing the selection procedure of the studied samples.
Figure 1. Flow chart describing the selection procedure of the studied samples.
Toxins 13 00103 g001
Figure 2. Comparison of the prevalence of eight mycotoxins in urine (a) among different age groups: 20–45 years (n = 116), 46–64 years (n = 71), 65–88 years (n = 40); (b) between males (n = 114) and females (n = 113) in the three regions and (c) in the different BMI groups: <18.5 kg m−2 (n = 13), 18.5–23.9 kg m−2 (n = 128), >24 kg m−2 (n = 86). (AFM1: aflatoxin M1; OTA: ochratoxin A; FB1: fumonisin B1; DON: deoxynivalenol; DON-3-GlcA: deoxynivalenol-3-glucuronide; DON-15-GlcA: deoxynivalenol-15-glucuronide; ZEN: zearalenol; ZEN-14-GlcA: zearalenone-14-glucuronide; ALL: mycotoxin contamination in any case; *: 0.01 < p-value < 0.05; ** p-value < 0.01 by the chi-square test).
Figure 2. Comparison of the prevalence of eight mycotoxins in urine (a) among different age groups: 20–45 years (n = 116), 46–64 years (n = 71), 65–88 years (n = 40); (b) between males (n = 114) and females (n = 113) in the three regions and (c) in the different BMI groups: <18.5 kg m−2 (n = 13), 18.5–23.9 kg m−2 (n = 128), >24 kg m−2 (n = 86). (AFM1: aflatoxin M1; OTA: ochratoxin A; FB1: fumonisin B1; DON: deoxynivalenol; DON-3-GlcA: deoxynivalenol-3-glucuronide; DON-15-GlcA: deoxynivalenol-15-glucuronide; ZEN: zearalenol; ZEN-14-GlcA: zearalenone-14-glucuronide; ALL: mycotoxin contamination in any case; *: 0.01 < p-value < 0.05; ** p-value < 0.01 by the chi-square test).
Toxins 13 00103 g002
Table 1. Occurrence of mycotoxins/metabolites in human urine samples in the Yangtze River Delta, China.
Table 1. Occurrence of mycotoxins/metabolites in human urine samples in the Yangtze River Delta, China.
MycotoxinsAll (n = 227)Jiangsu (n = 76)Zhejiang (n = 76)Shanghai (n = 75)p-Value 1
AFM1
n (prevalence, %)5 (2.20)1 (1.32)1 (1.32)3 (4.00)
Mean ± SD 2 (ng mL−1)0.35 ± 0.170.100.430.41 ± 0.120.344
Range (ng mL−1)0.10–0.550.100.430.31–0.55
Mean ± SD (ng mg−1 Cr) 30.23 ± 0.070.120.270.26 ± 0.040.344
Range (ng mg−1 Cr)0.12–0.290.120.270.21–0.29
OTA
n (prevalence, %)8 (3.52)2 (2.63)1 (1.32)5 (6.67)
Mean ± SD (ng mL−1)0.14 ± 0.060.16 ± 0.080.160.13 ± 0.060.717
Range (ng mL−1)0.05–0.220.11–0.220.160.05–0.21
Mean ± SD (ng mg−1 Cr)0.21 ± 0.120.14 ± 0.110.280.22 ± 0.120.699
Range (ng mg−1 Cr)0.05–0.370.06–0.220.280.05–0.37
FB1
n (prevalence, %)28 (12.33)15 (19.74)10 (13.16)3 (4.00)
Mean ± SD (ng mL−1)1.12 ± 0.401.18 ± 0.361.01 ± 0.351.17 ± 0.680.416
Range (ng mL−1)0.51–1.960.51–1.760.51–1.680.75–1.96
Mean ± SD (ng mg−1 Cr)1.12 ± 0.740.88 ± 0.441.33 ± 0.981.56 ± 0.870.326
Range (ng mg−1 Cr)0.30–3.340.30–1.820.37–3.340.57–2.20
DON
n (prevalence, %)28 (12.33)10 (13.16)6 (7.89)12 (16.00)
Mean ± SD (ng mL−1)1.79 ± 1.863.29 ± 2.391.63 ± 0.750.64 ± 0.340.000 *
Range (ng mL−1)0.50–8.610.50–8.611.09–3.000.50–1.56
Mean ± SD (ng mg−1 Cr)1.50 ± 2.193.17 ± 3.050.95 ± 0.340.38 ± 0.150.000 *
Range (ng mg−1 Cr)0.17–9.880.52–9.880.56–1.460.17–0.68
DON-3-GlcA
n (prevalence, %)23 (10.13)10 (13.16)7 (9.21)6 (8.00)0.014 *
Mean ± SD (ng mL−1)1.15 ± 1.001.79 ± 1.230.60 ± 0.260.73 ± 0.35
Range (ng mL−1)0.50–4.200.50–4.200.50–1.180.50–1.23
Mean ± SD (ng mg−1 Cr)1.10 ± 1.431.89 ± 1.900.45 ± 0.370.56 ± 0.310.004*
Range (ng mg−1 Cr)0.20–6.990.53–6.990.20–1.220.32–1.18
DON-15-GlcA
n (prevalence, %)69 (30.40)22 (28.95)19 (25.00)28 (37.33)
Mean ± SD (ng mL−1)2.55 ± 4.685.09 ± 7.621.52 ± 1.471.26 ± 1.000.016*
Range (ng mL−1)0.50–35.000.50–35.000.50–6.500.50–3.39
Mean ± SD (ng mg−1 Cr)2.16 ± 3.504.43 ± 5.451.19 ± 1.111.04 ± 0.820.033*
Range (ng mg−1 Cr)0.17–20.770.30–20.770.22–3.790.17–3.53
ZEN
n (prevalence, %)27 (11.89)19 (25.00)4 (5.26)4 (5.33)
Mean ± SD (ng mL−1)0.97 ± 3.491.22 ± 4.160.40 ± 0.440.38 ± 0.560.868
Range (ng mL−1)0.10–18.350.10–18.350.10–1.040.10–1.23
Mean ± SD (ng mg−1 Cr)0.65 ± 2.070.81 ± 2.460.22 ± 0.200.31 ± 0.270.771
Range (ng mg−1 Cr)0.03–10.930.04–10.930.03–0.510.04–0.68
ZEN-14-GlcA
n (prevalence, %)2 (0.88)-1 (1.32)1 (1.33)
Mean ± SD (ng mL−1)0.69 ± 0.62-1.130.250.317
Range (ng mL−1)0.25–1.13-1.130.25
Mean ± SD (ng mg−1 Cr)0.37 ± 0.18-0.500.240.317
Range (ng mg−1 Cr)0.24–0.50-0.500.24
1p-values obtained using the Kruskal–Wallis test, which compared mycotoxin concentrations in the three regions, * p-value < 0.05. 2 Only positive samples (≥limit of detection (LOD)) were considered for calculation of mean values. (LOD: aflatoxin M1 (AFM1), 0.05 ng mL−1; ochratoxin A (OTA), 0.05 ng mL−1; fumonisin B1 (FB1), 0.1 ng mL−1; deoxynivalenol (DON), 0.5 ng mL−1; DON-3-GlcA, 0.5 ng mL−1; zearalenone (ZEN), 0.1 ng mL−1; ZEN-14-GlcA, 0.2 ng mL−1). 3 Creatinine-adjusted concentration (ng mg−1 Cr).
Table 2. Correlations (rs and p-values) between urinary mycotoxin concentrations and food consumption.
Table 2. Correlations (rs and p-values) between urinary mycotoxin concentrations and food consumption.
FoodsAFM1OTAFB1DON 1ZEN
rsprsprsprsprsp
Wheat0.2240.7180.3090.457−0.3550.064−0.1250.2920.0950.638
Maize--−0.2470.5550.1920.3280.0570.631−0.1190.554
Rice−0.0530.933−0.1260.766−0.1110.573−0.0740.5360.2900.142
Vegetables and fruit−0.1540.805−0.2470.555−0.1910.331−0.1060.374−0.0890.660
Meat0.0001.000−0.1790.6720.1910.3320.0670.5740.1510.451
Nuts and seeds0.8660.0580.1090.797−0.1980.312−0.1780.1320.1270.529
Milk and dairy produce0.2890.638−0.2750.5100.3620.058−0.0670.576−0.2060.304
Beverages, coffee and tea−0.2240.7180.7780.023 *20.0030.9900.1860.1160.1530.447
1 DON: the sum concentration of DON, DON-3-GlcA, and DON-15-GlcA. *2 p < 0.05.
Table 3. Risk assessments of individual mycotoxins by the margin of exposure (MOE) and the hazard quotient (HQ) in adults from the Yangtze River Delta, China (n = 227).
Table 3. Risk assessments of individual mycotoxins by the margin of exposure (MOE) and the hazard quotient (HQ) in adults from the Yangtze River Delta, China (n = 227).
BiomarkerModelMean PDI 1N (%) 2MedianP95Max
AFB1AFM1MOE0.622.200.672.142.50
OTAOTAMOE0.143.52110.31217.05270.67
FB1FB1HQ5.5112.335.248.5810.5
DONDON+DON-3-GlcA+DON-15-GlcAHQ0.0800.030.280.85
ZENZENHQ0.200.880.131.1813.78
1 Only positive samples (≥ LOD) were considered for the calculation of dietary intake. 2 N (%), percentage of samples posing health risks (percent in all, %).
Table 4. Cumulative risk assessment of co-occurring mycotoxins in adults using combined margin of exposure (MOET) and hazard index (HI) models, respectively.
Table 4. Cumulative risk assessment of co-occurring mycotoxins in adults using combined margin of exposure (MOET) and hazard index (HI) models, respectively.
GroupMycotoxinToxicological EffectMOE/HQ
MeanMaxN(n) 1
1AFB1hepatotoxicity0.610.672 (2)
FB1hepatotoxicity19.4620.98-
MOET 2hepatotoxicity0.590.652 (2)
2DONmale reproductive0.020.040 (12)
ZENestrogenic0.370.990 (12)
HI 3endocrine effect0.391.022 (12)
1 N(n): Number of HQ/HI values > 1 or MOE < 10,000 (total number of co-contaminated samples). 2 MOET = 1/(1/MOEAFB1 + 1/MOEFB1). 3 HIendocrine effect = HQNOAEL, DON + HQTDI, ZEN.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Huang, Q.; Jiang, K.; Tang, Z.; Fan, K.; Meng, J.; Nie, D.; Zhao, Z.; Wu, Y.; Han, Z. Exposure Assessment of Multiple Mycotoxins and Cumulative Health Risk Assessment: A Biomonitoring-Based Study in the Yangtze River Delta, China. Toxins 2021, 13, 103. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020103

AMA Style

Huang Q, Jiang K, Tang Z, Fan K, Meng J, Nie D, Zhao Z, Wu Y, Han Z. Exposure Assessment of Multiple Mycotoxins and Cumulative Health Risk Assessment: A Biomonitoring-Based Study in the Yangtze River Delta, China. Toxins. 2021; 13(2):103. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020103

Chicago/Turabian Style

Huang, Qingwen, Keqiu Jiang, Zhanmin Tang, Kai Fan, Jiajia Meng, Dongxia Nie, Zhihui Zhao, Yongjiang Wu, and Zheng Han. 2021. "Exposure Assessment of Multiple Mycotoxins and Cumulative Health Risk Assessment: A Biomonitoring-Based Study in the Yangtze River Delta, China" Toxins 13, no. 2: 103. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13020103

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