Next Article in Journal
Measurement of Fumonisins in Maize Using a Portable Mass Spectrometer
Next Article in Special Issue
Industrial-Scale Cleaning Solutions for the Reduction of Fusarium Toxins in Maize
Previous Article in Journal
Equatorial Spitting Cobra (Naja sumatrana) from Malaysia (Negeri Sembilan and Penang), Southern Thailand, and Sumatra: Comparative Venom Proteomics, Immunoreactivity and Cross-Neutralization by Antivenom
Previous Article in Special Issue
The Role of Nitrogen Fertilization on the Occurrence of Regulated, Modified and Emerging Mycotoxins and Fungal Metabolites in Maize Kernels
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multi-Mycotoxin Long-Term Monitoring Survey on North-Italian Maize over an 11-Year Period (2011–2021): The Co-Occurrence of Regulated, Masked and Emerging Mycotoxins and Fungal Metabolites

1
Council for Agricultural Research and Economics (CREA), Research Centre for Cereal and Industrial Crops, Via Stezzano, 24, 24126 Bergamo, Italy
2
Department of Agricultural, Forest and Food Sciences (DISAFA), Università degli Studi di Torino, Largo Braccini 2, 10095 Grugliasco, Italy
3
Council for Agricultural Research and Economics (CREA), Research Centre for Engineering and Agro-Food Processing, Via Milano 43, 24047 Treviglio, Italy
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Toxins 2022, 14(8), 520; https://doi.org/10.3390/toxins14080520
Submission received: 29 June 2022 / Revised: 26 July 2022 / Accepted: 27 July 2022 / Published: 29 July 2022

Abstract

:
Maize is considered one of the most susceptible crops to mycotoxin-producing fungi throughout the world, mainly belonging to the Fusarium spp. and Aspergillus spp. Maize is mainly used as animal feeds in Italy, as well as for human consumption, being essential for all the protected designation of origin (DOP) products. Our study investigated the occurrence of regulated mycotoxins in 3769 maize grain samples collected from 88 storage centers by the National Monitoring Network over an 11-year period (2011–2021). Moreover, an in-depth survey over a 4-year period, characterized by extremely different meteorological conditions, was conducted to investigate the co-occurrence of regulated, masked, and emerging mycotoxins. The survey confirmed that Fusarium spp. was the most frequent fungi and fumonisins were the main mycotoxins that were constantly detected in the different years and areas. Moreover, the areas characterized by high fumonisin levels were also the most prone to contamination by emerging mycotoxins produced by the same Fusarium species of the Liseola section. On the other hand, as a result of climatic changes, maize grains have also been affected by the increased frequency of aflatoxin accumulation. Deoxynivalenol, zearalenone, and other emerging mycotoxins produced by the same Fusarium species as the Discolor section occurred more abundantly in some areas in Northern Italy and in years characterized by predisposing meteorological conditions.
Key Contribution: The mycotoxin contaminations of both regulated and emerging mycotoxins showed a pronounced year-by-year variation due to both the geographical characteristics of the maize growing area and the meteorological conditions during sensitive periods of maize grain development, especially during the ripening period. The definition of risk areas/zones is essential to understand the impact of mycotoxin contaminations on the food and feed chains, in consideration of a changing climate change scenario.

1. Introduction

Maize (Zea mays) is the most important grain and forage crop in Italy. It is mainly used as an energy ingredient in livestock feeds, but its utilization in many food products and bioindustrial derivatives is increasing to a great extent [1]. Maize grains are subject to infection from several fungal diseases of the Aspergillus, Fusarium, and Penicillium genera. The colonization and growth of these fungi can lead to an altered grain quality, and this mainly occurs through mycotoxin contamination [2].
Because of the frequent occurrence and negative impact of mycotoxins on human and animal health, several countries have established mycotoxin regulations for food and feeds. In the European Union (EU), maximum concentrations of aflatoxin B1 (AF B1), deoxynivalenol (DON), zearalenone (ZEA), ochratoxin A (OTA), and the sum of fumonisin B1 (FB1) and B2 (FB2) are in place for foodstuffs, while guidance levels have been introduced for T2-HT2 [3,4]. As of yet, maximum limits for feedstuffs are only in force for AF B1, while guidance levels have been set for DON, ZEN, OTA and the sum of FB1 and FB2 [5,6]. Currently, the EU regulations do not take into account the occurrence of masked and emerging mycotoxins or secondary fungal metabolites and the presence of multiple mycotoxins. Although mycotoxin regulations globally focus on a few mycotoxins individually, maize grains are frequently contaminated simultaneously by several different mycotoxins, including conjugated and emerging mycotoxins.
Maize, due to its sensitivity to infection by toxigenic molds and to the ecological conditions during the growth and ripening of the grain, is the crop that is the most sensitive to high mycotoxin contamination in Southern Europe [7]. This is particularly evident in the Po Valley, where 95% of all the Italian grain production is concentrated. The maize grown in the plains of Northern Italy is cultivated in an environment with high air humidity and high temperatures, which are extremely favorable for the development of the main toxigenic molds. In fact, high contaminations of mycotoxins frequently occur in commercial lots under these environmental conditions [8]. For this reason, over the last two decades, thanks to the increasing attention paid to safety aspects by cereal supply chains, it has been found that mycotoxin contamination is actually the first cause of the decrease in the market value of maize grains [9]. Furthermore, a frequent and severe contamination by mycotoxins has been indicated as one of the main causes of a significant contraction of the maize cultivation area in Italy; in the past 2 decades, 50% of the acreage has been lost and converted for the cultivation of crops less prone to contamination, such as soybean or small cereals.
In the past decades, several surveys have been published around the world, including in Northern Italy, to highlight mycotoxin contamination profiles [7,10,11,12,13,14]. However, only a few surveys have investigated the mycotoxin profiles of regulated mycotoxins, masked, emerging mycotoxins, and secondary fungal metabolites [15,16]. Moreover, no other study has focused on maize production in Northern Italy by analyzing 3769 grain samples from 88 storage centers over an 11-year period.
Thus, the aims of this survey were twofold. The first aim was to draw up mycotoxin profiles of the maize grain grown in Northern Italy to highlight the different occurrences between years and production areas in order to point out the environmental conditions that aggravate the contaminations of regulated mycotoxins. The second aim was to investigate the role of emerging mycotoxins and their co-occurrence with the regulated ones.

2. Results

2.1. National Maize Monitoring Network

Our study investigated the occurrence of regulated mycotoxins in maize grain samples collected by the Italian National Monitoring Network, supported by the Italian Ministry of Agriculture, over an 11-year period (2011–2021). This network, which is coordinated by the CREA Research Centre for Cereal and Industrial Crops, is composed of 88 maize storage centers distributed throughout five macro-geographical areas (Figure 1). The Po Valley is the main maize growing area in Italy.

2.2. Metereological Trends

We analyzed the different meteorological trends observed over the 2011–2021 period by measuring the daily mean precipitation and temperature during the ripening period of maize kernels (Figure 2A).
The rainfall range observed over the years, where the average of 10 meteorological stations was considered, was between 152.08 mm in 2011 and 350.33 mm in 2014. Differences were also observed regarding temperatures. The minimum, 20.75 °C, was observed during 2014, while the maximum, 23.76 °C, was noted during 2018. The 2015 growing season also showed high temperatures (23.46 °C), which were counterbalanced by low rainfall (155.62 mm) with respect to 2018 (209.83 mm). As can be seen in Figure 2, the 2014 growing period underwent both more abundant rainfall (350.33 mm) and lower temperatures (20.75 °C) than the other years. The 2017 and 2021 growing seasons also showed low temperature levels of 22.49 °C and 22.53 °C, respectively.
The mean meteorological trend was also observed considering the five geographical areas (Figure 2B). The precipitation and temperature data from the two meteorological stations chosen to represent each geographical area were averaged over the 2011–2021 years.
A different distribution of rainfall was observed within the different areas. The W, SP, and A areas appear to be the areas with the lowest levels of precipitation, with 142.75 mm (W), 126.09 mm (SP), and 170.85 mm (A), respectively. However, these geographical areas also showed substantial differences in temperature. The W area was in fact cooler than the SP and A ones. The mean temperature for W was 21.86 °C, while it was 22.69 °C and 23.09 °C for the SP and A areas, respectively. C and E were instead the geographical areas in which the highest mean precipitation levels were observed: 260.31 mm for C and 305.51 mm for E.

2.3. Regulated Mycotoxins

The occurrence of FBs and AF B1 in maize lots is summarized in Table 1. FBs were the most common mycotoxins, being present in at least 97% of the samples analyzed from 2011 to 2021. The statistical analysis showed that the years with the highest and lowest contaminations were 2019 and 2011. During these years, the contamination level ranged between 25 and 47,300 µg/kg. By analyzing the spread in different geographical areas, E appeared to be most contaminated while W was the cleanest over the years.
Throughout the years in which this study was carried out, the presence of AF B1 was lower than that of FBs, with an average level of between 0.6 in 2014 and 12.5 µg/kg in 2012.
In 2014, the cleanest year, AF B1 was only detected in 7% of the samples, while in 2019, the most infected year, 59% of them showed an AF B1 level above 1 µg/kg. The statistical analysis showed that the areas of greatest incidence were SP and A while the W area was the least contaminated.
The effects of the year and area on the occurrence of DON and ZEA in the maize lots were also detected (Table 2).
The years with the higher levels of contamination were 2013 and 2014 while the year 2021 was found to be the least contaminated.
In 2014, the maximum concentration reached 36,583 µg/kg, and 100% of the analyzed samples were positive.
When analyzing the geographical areas, W and A showed the highest and lowest DON concentration values.
Similarly, 2014 was the most contaminated year for ZEA, with an average value of 364 µg/kg, a maximum value of 1395 µg/kg, and 91% of positive samples; 2017 was the cleanest year. When considering the spread in different areas, W and E were found to be the most contaminated. The other areas had low levels of contamination, ranging from 43 µg/kg in C to 21 µg/kg in A.
The year with the lowest incidence and contamination of each mycotoxin, which we referred to as the best year, was compared with the worst year, namely the one with the highest incidence and contamination. The diffusion of regulated mycotoxins in different maize lots in Northern Italy, as obtained from the drying and storage centres, was also determined (Figure 3).
The year 2019 resulted to be the worst one for the incidence of FBs. All the areas had percentages of samples with FB concentrations above 4000 µg/kg, which is the maximum level for maize for human food use [3], and they ranged from a minimum of 67% in the SP area to 89% in the W ones. In 2011, the best year, the percentages of samples above the maximum level set in the EU were lower, that is, over a range of 1% in W to 18% in E (Figure 3A).
In the best years, such as 2014, the incidence of AF B1 was insignificant: from 97% in the W area to 100% in the C, SP and A areas, and the samples had a lower AF B1 concentration than 2 µg/kg. In 2012, in which AF B1 was widespread on maize, 50% of the samples from area A had a larger AF B1 concentration than 20 µg/kg, thereby far exceeding the maximum level of 8 µg/kg set by the European Commission (EC) for AF B1 for maize intended for human consumption [4] (Figure 3B).
DON and ZEA showed low incidences in the maize samples analyzed over the years. In 2021, 100% of the samples, in all the areas, had a DON content of less than 750 µg/kg; in 2017, the detected ZEA concentration in all the samples was below 350 µg/kg, that is, the maximum level set by the EC for maize for food purposes [3]. The year 2014 was the worst year for both mycotoxins in all the growing areas. DON had a greater concentration than 1750 µg/kg, the maximum level set by the EC for maize for human consumption [3], and ranged between 22% and 79% in the samples from areas A and W, respectively (Figure 3C). In the same areas, ZEA showed a greater concentration than 350 µg/kg, with values of 6% and 59%, respectively (Figure 3D).
Figure 4 shows the mean contamination level of FBs in maize lots from drying and storage centers located in different areas during the 2011–2021 period. According to the statistical analysis, the E area on average showed the highest level of FB accumulation in all the areas. Instead, the area with the lowest level of FBs was generally the W one. The years with the greatest accumulation of FBs were 2012, 2013, and 2019, while 2011 was the least-contaminated year.
The mean contamination levels of AF B1, regarding the same maize lots reported for FBs, during the 2012–2021 period, are shown in Figure 5. According to the statistical analysis, the A and SP areas on average showed higher AF B1 accumulation levels than all the other areas. The area with lowest level of AF B1 was generally the W area. The year 2012 showed the greatest accumulation of AF B1 in the different maize growing areas, except in the W area, due to the meteorological conditions that occurred, which favored AF B1 production by Aspergillus spp. In fact, except for a few years and in some geographical areas, the level of AF B1 was on average lower that that recorded in 2012.
The mean contamination level of DON is shown in Figure 6. According to the statistical analysis, 2013 and 2014 were the worst years, with the greatest accumulation of DON of all the years. In both years, the W and A areas were the most and the least contaminated, respectively. The best growing season, with the lowest levels of DON contamination, was 2021.
Figure 7 shows the mean contamination level of ZEA in the investigated maize lots during the 2014–2021 period. According to the statistical analysis, 2014 was the most contaminated year and W area showed higher levels of ZEA accumulation than all the other geographical areas. The A area was on average the cleanest one. The best year for ZEA accumulation was 2017.
Figure 8 shows the spread in different geographic areas for each mycotoxin, in the year of the highest incidence, considering the different classes of incidence established by Regulations (EC) n. 1126/2007 and n. 165/2010 [3,4].
In the worst year for FBs (2019, Figure 8A), almost the entire geographical area showed levels above 4000 µg/kg, except for a small area in E, where the levels were lower (1000–4000 µg/kg).
During 2012, the worst year for aflatoxin accumulation (Figure 8B), the lowest levels of AF B1 were found in W and some areas in North E and C (0–5 µg/kg). The other areas showed distributions that fell into the highest classes. In the worst year for DON incidence (2014, Figure 8C), the SP (near Bologna) and A areas showed concentration levels included in the lowest classes, 0–1750 µg/kg. The W, C, and E areas and part of SP showed DON levels above the limit of 1750 µg/kg. A similar distribution was observed for ZEA (Figure 8D) with a lower level of incidence than that of DON.
Figure 9 shows the mean concentration of each mycotoxin detected during the investigated years. FBs (Figure 9A) were found to be endemically present in different geographical areas, with average distributed levels of between 1000 and >4000 µg/kg. As far as the AF B1 distribution is concerned, the average levels were generally distributed over the 0–10 µg/kg classes. Some small areas, where the average accumulation was greater (10–20 µg/kg), are highlighted in the figure. The distributions of DON and ZEA over the years generally showed an average distribution over all the considered incidence classes.
A Principal Component Analysis (PCA) was carried out to investigate the relationships between the main detected regulated mycotoxins, as it allows the maize samples to be clustered as similar mycotoxin contaminations for the different years and areas. The results of the PCA performed with the regulated mycotoxins detected over eleven growing seasons and in five maize areas in Northern Italy were used to reduce the number of variables to only 2, the 2 first principal components (PCs) producing a loading plot in which the two PCs, respectively, explained 39% (Principal Component 1, PC1) and 21% (Principal Component 2, PC2) of the total variations of the maize lot samples (Figure 10). The first component differentiated the years (Figure 10A) and the areas (Figure 10B) according to their mycotoxin contents, as produced by F. graminearum and F. culmorum, DON and ZEA, (loadings >|0.84|), and was also positively related to the rainfall during the ripening period and negatively related to the mean temperature during the ripening period. On the other hand, the second component differentiated the years (Figure 10A) and the areas (Figure 10B) according to their mycotoxin contents, as produced by F. verticillioides and F. proliferatum, FBs, and by Aspergillus spp. AF B1 (loadings > |0.71|). Overall, the PCA confirmed that the samples characterized by high contaminations of DON and ZEA were mainly clustered in the W area in 2014, while the SP and A areas as well as the 2012 (data not shown because not all the mycotoxins were analyzed for 2012), 2015, 2016, and 2017 years were more prone to AF B1 contamination. Furthermore, although FBs were constantly present in the different areas over the years, the E area and the 2019 year were the most contaminated ones.

2.4. Masked Mycotoxins

An in-depth study was carried out during the 2012–2015 period to evaluate the contamination levels of masked and modified mycotoxins, and emerging fungal metabolites in the same maize lots. As far as the masked mycotoxins are concerned, the mainly detected forms were those associated with DON, such as deoxynivalenol-3-glucoside (DON-3-G), and the acetylated forms, 3-acetyldeoxynivalenol (3-ADON) and 15-acetyldeoxynivalenol (15-ADON) (Table 3).
Overall, the highest contaminations were detected for DON-3-G, which on average reached about the 27% of the DON contamination level detected in the same period during the 2012–2015 period, and this was followed by 15-ADON, which showed 12% of the DON content, while 3-ADON represented only 2% of the DON contamination. The 2014 year was the most contaminated for both DON-3-G and 15-ADON and was significantly different from the other years. On the other hand, the 2013 year, which was the second-most-contaminated year, after 2013, for DON-3-G and 15-ADON, was also the most contaminated for 3-ADON.
As far as the areas are concerned, area W was significantly more prone to the contamination of all the masked forms of DON than the other areas, as previously reported for the “native” DON form (Table 3).

2.5. Emerging Mycotoxins and Other Fungal Metabolites

The emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. of the Liseola section, such as beauvericin (BEA), bikaverin (BIK), fusaric acid (FA), fusaproliferin (FUS), and moniliformin (MON), were detected for all the years and in all the areas. The effects of year and area on the occurrence of these emerging mycotoxins and fungal metabolites in maize lots from drying and storage centers located in different areas in Northern Italy during the 2012–2015 period are reported in Table 4.
The 2013 year was the one with the highest level of contamination for all the detected fungal metabolites produced by Fusarium spp. of the Liseola section, with the exception of MON, for which 2012 was significantly more contaminated than the other years. As far as the areas are concerned, no significant differences among the different areas were recorded for any of the fungal metabolites reported in Table 4, with the exception of BIK and FA, for which area A recorded a significantly higher contamination than the other areas.
Moreover, the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. of the Discolor section, such as aurofusarin (AUR), butenolide (BUT), culmorin (CULM), and nivalenol (NIV), were also detected for all the years and areas considered, as previously reported (Table 5).
The 2014 year was the one with the highest level of contamination for all the detected fungal metabolites produced by Fusarium spp. of the Discolor section, with the exception of CULM, for which 2013 was the most contaminated year. As has already been described for the regulated mycotoxins, DON and ZEA, produced by Fusarium spp. of the Discolor section and for the DON-associated masked forms, area W was also one of the most contaminated areas for these emerging fungal metabolites, with the exception of NIV, for which no significant differences were recorded between the areas (Table 5).
A PCA was also carried out to investigate the relationships between the main detected emerging mycotoxins, fungal metabolites, and regulated mycotoxins, as it allowed the maize samples to be clustered for similar mycotoxin contaminations according to the different years and areas. The results of the PCA performed with the regulated and emerging mycotoxins and other fungal metabolites detected in four growing seasons and in five maize areas in Northern Italy were used to produce a loading plot in which the two first principal components again explained 39% and 21% of the total variations of the maize lot samples (Figure 11). The first component differentiated the years (Figure 11A) and the areas (Figure 11B) according to their mycotoxin contents produced by Fusarium spp. of the Discolor section and was also positively related to the rainfall during the ripening period and negatively related to the mean temperature during the ripening period, thus confirming what had previously been reported for DON and ZEA. Indeed, the contaminations of DON-3-G, 3-ADON, 15-ADON, CULM, AUR, and BUT, which are mainly produced by F. graminearum and F. culmorum (Discolor section), were found to be associated with DON and ZEA (loadings > |0.77|).
The second component instead differentiated the years (Figure 11A) and the areas (Figure 11B) according to the mycotoxin contents produced by Fusarium spp. of the Liseola section. Indeed, the contaminations of MON, FUS, and BEA, which are mainly produced by F. subglutinans, F. proliferatum, F. verticillioides (Liseola section), and F. avenaceum (section Roseum), were found to be associated with FBs, were positively correlated with the ECB insecticide index (loadings > |0.5|), and were expressed with an arbitrarily assigned number equal to 1 if the maize lot had not been treated with an insecticide, or equal to 0.5 if the lot had instead been treated with a pyrethroid.

3. Discussion

This study evaluated, for the first time, mycotoxin profiles obtained from the contamination of regulated, masked, and emerging mycotoxins and fungal metabolites in maize grown in Northern Italy through an extensive survey conducted on maize grain lots from 88 storage centers over the 2011–2021 period, year by year.
When considering the whole area cultivated with maize, the mycotoxin concentrations showed a pronounced year-to-year variation that can be explained by the meteorological conditions during sensitive periods of maize grain development, especially during the ripening period. Several control points, such as the selection and development of resistant cultivars through conventional breeding, Good Agricultural Practices (GAP), chemical and biological control during cultivation, and proper management during storage, may in fact reduce fungal infestation and growth, as well as mycotoxin production in cereal grains. Furthermore, fungal infestation and mycotoxin production depend to a great extent on environmental factors, such as temperature, rainfall, and moisture, at any given stage [17].
Almost all the collected samples were contaminated by at least one regulated mycotoxin. A large number of samples (>85%) were co-contaminated with ≥2 mycotoxins. Single species of fungi may actually produce more than one mycotoxin concomitantly, and different fungal species may proliferate and co-exist in the same plant, both circumstances that can lead to the co-occurrence of two or more mycotoxins [18]. Thus, mycotoxin co-contamination is common and is in part governed by the climate; thus, regional trends have emerged throughout the world [7].
The most frequently observed mycotoxin mixtures were combinations of regulated mycotoxins AF B1, FBs, as well as DON and ZEA. These results agree with those previously reported by several authors [7,15,17,19,20,21], who noted higher incidences in Southern Europe than in the rest of Europe [7]. It should be highlighted that our data, pertaining to all the detected regulated mycotoxins (AF B1, FBs, DON and ZEA), showed much higher incidences and ranges of concentrations than those reported about Europe in several surveys in the literature, as summarized by Lee and Ryu [17], especially for AF B1, FBs, and DON. Indeed, the incidences reported by Lee and Ryu [17] for AF B1, FBs, and DON were 7%, 49%, and 48%, respectively, while the incidence values recorded over the years in our survey were on average equal to 41%, 98%, and 85%, respectively. Moreover, the ranges and maximum levels at which AF B1 (<LOD–145 µg/kg), FBs (< LOD–47,300 µg/kg), and DON (<LOD–36,731 µg/kg) were detected were clearly different from those reported by Lee and Ryu [17] (AF B1 = 1.6 µg/kg, FBs = 25–4438 µg/kg and DON = 0.04–2942 µg/kg). Our data confirmed that maize grown in Northern Italy is more prone to contamination than that detected by surveys conducted in other areas.
In addition to the regulated mycotoxins, twelve masked, modified, and emerging mycotoxins and fungal metabolites were also detected in the same maize lots during the in-depth survey conducted over the 2012–2015 period. The co-occurrence of regulated and emerging mycotoxins and fungal metabolites in maize has already been reported throughout Europe for maize for feed and food purposes by several authors [22,23,24,25,26], and our results therefore confirmed the risks related to the presence of emerging mycotoxins and fungal metabolites in food and feeds.
As far as the relationship between regulated and emerging mycotoxins is concerned, the FB concentrations in the maize lots were correlated with all the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. section Liseola and with ECB damage, thus confirming what has extensively been reported in literature for FBs in North Italian maize surveys [27,28], but which has only been reported for emerging mycotoxins by Blandino et al. in 2015 [29]. Indeed, damage by insects or hail may boost infections because it can create wounds that may readily be colonized by fungi.
On the other hand, the DON and ZEA concentrations in the maize were correlated with all the detected masked forms of DON and all the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. section Discolor. A positive correlation with rainfall and a negative correlation with the average temperatures measured during the ripening period were also observed. Our results are in agreement with those of Gaikpa and Miedaner, who reported that high levels of infections occurred when the weather is cool and wet [30]. Other authors have also reported that Fusarium fungi on ears and grains increase under favorable weather conditions, such as higher temperatures and higher rainfall levels in summer, and this results in a longer growing season [31]. Blandino et al. [32] reported that the highest DON and other Fusarium spp. Discolor section mycotoxin contaminations of maize grain occurred in growing seasons characterized by a high level of precipitation and lower temperatures in the period from maize flowering to harvest.
AF B1 was instead detected at the highest levels of contamination in the years and areas characterized by less rainfall and higher temperatures during summer. Indeed, Aspergillus ear rot (AER) has frequently been reported in hot and dry environments, such as Southern USA, Africa, and tropical Asia [33], but has also recently been observed in Southern Europe, especially in Northern Italy, Serbia, Slovenia, Croatia, Romania, and Hungary [34,35].
The evaluation of the co-occurrence of the main regulated mycotoxins in maize (AF B1, FBs, DON, and ZEA) has also made it possible to define risk areas within the Po Valley. These areas, which are particularly prone to heavy contamination of different mycotoxins, are characterized by different average meteorological conditions, concerning both rainfall and temperatures. Furthermore, these differences could be aggravated by the seasonal meteorological trend. As previously reported, the meteorological conditions in an area are one of the most important key factors for mycotoxin contamination, together with insect damage, to which maize is known to be prone [27,28,29,30,31,32]. Since a correlation has been highlighted between regulated mycotoxins and emerging and masked mycotoxins, which mainly depends on their affinity to the same fungi that produce them, the risk zones/areas described in the present study could be extended to include emerging mycotoxins and fungal metabolites related to the regulated ones.
The definition of risk areas/zones is essential to understand the repercussions of mycotoxin contamination on the food and feed chains, especially in view of a climate change scenario. Indeed, mycotoxins are present at every step of the agri-food chain due to their ability to move through different levels of the food pyramid (carry-over). However, their negative effects are more severe in agricultural crops and dairy farms due to direct losses in harvest, production profitability, animal health, and safety of their products [2,36].
The results of our survey could also be used to identify Good Agricultural Practices (GAP) and structural investments to help rule out mycotoxin contamination.
In conclusion, this survey not only showed the distribution of mycotoxin contamination in a relevant maize region of the EU but also provided information on the role of the environment on mycotoxins in order to enhance the maize quality of commercial lots, considering the required specifications of the food and feed chains.

4. Materials and Methods

4.1. Meteorological Data

Climatic parameters, such as the mean rainfall (mm) and mean temperature (°C), were collected during the ripening period (July–September) of maize kernels, the period that appears to be the most closely correlated with the accumulation of mycotoxins. Data were taken from different ARPAE (Agenzia Regionale per la Protezione Ambientale) meteorological stations located in the main Italian maize growing regions. Northern Italy has been geographically divided into five main growing areas: West (W), Center (C), East (E), Adriatic (A), and South Po (SP) to cover the regions with high maize production. Ten meteorological stations were contacted to obtain information: Trino (VC) and Cussanio (CN)—West, Alessandria (AL) and Gonzaga (MN)—South Po River, Bergamo (BG) and Soncino (CR)—Center, Sant’Apollinare (RO) and Ospital Monacale (FE)—Adriatic and Treviso (TV), and Udine (UD)—East. The daily temperatures and precipitation rates were measured at the 10 different stations. Each storage center was connected to the nearest meteorological station, which was always located within a radius of 30 km.

4.2. Samples

A total of 3769 representative maize grain samples were collected, after dry processing, over an 11-year period (2011–2021) from 88 storage centers. These storage centers are distributed in five geographical areas (West, Center, South Po, Adriatic, and East), which represent the regions in Northern Italy with the greatest cultivation of maize. A total of 5–10 samples were collected from each storage center each year. A dynamic grain sampling strategy [37] was adopted on the moving product to obtain a representative sample [38]. Accordingly, we took a 15–20 kg pre-sample from all the production components (whole and broken grains and small parts). By applying the standard procedures, sorter samples of 1–1.5 kg each were obtained from these pre-samples for laboratory analyses. These samples were placed in sealed bags and stored in a cool room for 1 day until milling, according to European recommendation [39]. All the samples were subsequently milled to a 1 mm sieve size with a ZM 200 Retsch Ultra-Centrifugal mill equipped with a DR 100 vibratory feeder (Retsch GmbH, Haan, Germany) and stored at 4 °C until chemical analyses.

4.3. Regulated Mycotoxins Analyses

The regulated mycotoxin concentration levels in the maize samples collected during the 2011–2021 period were determined by means of an Enzyme-Linked Immunoassorbent Assay (ELISA). Ridascreen® R-Biopharm kit tests (R3401 RIDASCREEN Fumonisin, R1211 RIDASCREEN Aflatoxin B1 30/15, R5906 RIDASCREEN DON, R1401 RIDASCREEN Zearalenon) were performed using Chemwell Automatic Awareness Engineer (inc.). The mycotoxin extraction and tests were performed according to the manufacturers’ instructions. The Limits of Detection (LODs) are declared by the manufacturer: LOD of FBs = 25 µg/kg, LOD of AF B1 = 1 µg/kg, LOD of DON = 18.5 µg/kg, and LOD of ZEA = 1.75 µg/kg.

4.4. Multi-Mycotoxin LC-MS/MS Analysis

The maize lot samples from the drying and storage centers, collected during the 2012–2015 period, were analyzed by means of a multi-mycotoxin LC-MS/MS analysis. The extraction phase was described in detail by Sulyok et al. in 2006 [40]. The chromatographic and mass spectrometric parameters of the investigated analytes were described by Malachova et al. in 2014 [41]. Briefly, representative 5 g sub-samples of the milled material were extracted using 20 mL of a mixture of acetonitrile/water/acetic acid 79/20/1 (v/v/v). After extraction, the samples were centrifuged, diluted 1 + 1 by a mixture of acetonitrile/water/acetic acid 20/79/1 (v/v/v), and injected.
Detection and quantification were performed with a QTrap 5500 LC–MS/MS System (Applied Biosystems, Foster City, CA, USA), equipped with a TurboIonSpray electrospray ionization (ESI) source and a 1290 Series HPLC System (Agilent, Waldbronn, Germany). Chromatographic separation was performed, at 25 °C, on a Gemini® C18-column, 150 × 4.6 mm i.d., 5 μm particle size, equipped with a C18 security guard cartridge, 4 × 3 mm i.d. (all from Phenomenex, Torrance, CA, USA). Quantification was performed with an external calibration, and the results were corrected for any apparent recoveries determined in the maize. The accuracy of the method was verified by participating in proficiency testing schemes organized by BIPEA (Gennevilliers, France), with 160 out of the 168 results submitted for maize and maize-based feed exhibiting a z-score of between −2 and 2.

4.5. Statistical Analysis and Data Elaboration

The Kolmogorov–Smirnov normality test and the Levene test were carried out to verify the normal distribution and the homogeneity of variances, respectively. Since the interaction between the area and year was significant, an analysis of variance (ANOVA) was performed separately for each year with the area as an independent factor for mycotoxin contamination. Multiple comparison tests were performed, by means of the Ryan–Einot–Gabriel–Welsh F (REGW-F) test, on the treatment means. Principal Component Analysis (PCA) was carried out to investigate the relationships between the main detected mycotoxins. All the data obtained were joined together before the PCA was performed with the purpose to reduce the dimensionality of the data set by finding a new set of variables, the two principal components (PC1 and PC2), smaller than the original set of variables, that nonetheless retained most of the sample’s information and the total variance associated. A statistical data analysis was carried out with the SPSS software package, version 24.0.
The data from the monitoring of mycotoxins from the various geo-referenced locations were processed with R software [42] with the sp, raster, rgdal, and leaflet packages, which is specifically used for geostatistical analysis and geolocalized representation.
The distribution maps were created through the Kriging interpolation process by setting the processing grid with a mesh of 50 × 50 (2500) nodes and a linear function.

Author Contributions

Conceptualization, S.L. and V.S.; formal analysis, S.L., V.S., E.R. and C.L.; investigation, S.L., V.S. and C.L.; data curation, S.L., V.S., C.L. and E.R.; writing—original draft preparation, S.L., V.S. and C.L.; writing—review and editing, A.R.; supervision, S.L. and A.R.; funding acquisition, S.L. and A.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was undertaken with the financial assistance of the Italian Ministry of Agricultural Food and Forestry Policies-MiPAAF as part of the following research programs: MICOPRINCEM, RQC MAIS (D.D. N. 88666, 3 December 2014) and RETI2020 (prot. 198541, 30 April 2021). Moreover, in 2018, CREA carried out this work in collaboration with ISMEA as part of the “Cereal Quality Territorial Observatory” project of the National Cereal Plan, funded by the Italian Ministry of Agricultural, Food and Forestry Policies (MIPAAF).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the corresponding author, [S.L.], upon reasonable request.

Acknowledgments

The authors would like to thank Stefania Mascheroni and Iride Mascheretti for their technical assistance. Thanks are due to the storage centers belonging to the Italian National Maize Monitoring Network. The authors would also like to thank Michael Sulyok, for his precious help and cooperation in the multi-mycotoxin analysis; Massimo Blandino and Francesca Vanara for their precious help and cooperation in the supervision and laboratory work, respectively and Marguerite Jones for the English editing. The meterological data were provided by ARPAE (Agenzia Regionale per la Protezione Ambientale), Arpae Regione Piemonte-Rete Agrometeorologia regionale (RAM)-Settore Fitosanitario e Servizi Tecnico-scientific—Sezione di Agrometeorologia, Arpae Regione Emilia-Romagna-Osservatorio clima, Arpae Regione Lombardia-meterologia, and Arpae Friuli Venezia Giulia-osservatorio meterologico regionale.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
15-DON15-acetyldeoxynivalenol
3-DON3-acetyldeoxynivalenol
AAdriatic
AF B1Aflatoxin B1
AURAurofusarin
BEABeauvericin
BIKBikaverin
BUTButenolide
CCenter
CULMCulmorin
DONDeoxynivalenol
DON-3-GDeoxynivalenol-3-glucoside
DOPProtected designation of origin
EEast
ECEuropean Commission
ECBEuropean Corn Borer
ELISAEnzyme linked immunosorbent assay
EUEuropean Union
FAFusaric acid
FBsFumonisins (the sum of FB1 and FB2)
FUSFusaproliferin
GAPGood Agricultural Practices
LC–MS/MSLiquid Chromatography tandem Mass Spectrometry
LODLimit of Detection
MONMoniliformin
NIVNivalenol
OTAOchratoxin A
PCAPrincipal Component Analysis
SPSouth Po
WWest
ZEAZearalenone

References

  1. Rai, S.; Kaur, A.; Chopra, C.S. Gluten-free products for celiac susceptible people. Front. Nutr. 2018, 5, 116. [Google Scholar] [CrossRef] [PubMed]
  2. Munkvold, G.P.; Arias, S.; Taschl, I.; Gruber-Dorninger, C. Mycotoxins in corn: Occurrence, impacts, and management. In Corn: Chemistry and Technology, 3rd ed.; Serna-Saldivar, S.O., Ed.; AACC International Press: Washington, DC, USA, 2019; pp. 235–287. ISBN 9780128119716. [Google Scholar] [CrossRef]
  3. European Commission. Commission Regulation (EC) 1126/2007 of the European Parliament and of the Council of 28 September 2007 amending Regulation (EC) No 1881/2006 setting maximum levels for certain contaminants in foodstuffs as regards Fusarium toxins in maize and maize products. Off. J. Eur. Union 2007, L255, 14–17. [Google Scholar]
  4. European Commission. Commission Regulation (EC) No 165/2010 of 26 February 2010 amending Regulation (EC) No 1881/2006 setting maximum levels for certain contaminants in foodstuffs as regards aflatoxins. Off. J. Eur. Union 2010, L50, 8–12. [Google Scholar]
  5. European Commission. Directive 2002/32/EC of the European Parliament and of the Council of 7 May 2002 on undesirable substances in animal feed. Off. J. Eur. Union 2002, L140, 10–22. [Google Scholar]
  6. Commission Recommendation. European Commission Recommendation N. 2006/576 on the presence of deoxynivalenol, zearalenone, ochratoxin A, T-2 and HT-2 and fumonisins in products intended for animal feeding. Off. J. Eur. Union 2006, 229, 7–9. [Google Scholar]
  7. Gruber-Dorninger, C.; Jenkins, T.; Schatzmayr, G. Global mycotoxin occurrence in feed: A ten-year survey. Toxins 2019, 11, 375. [Google Scholar] [CrossRef]
  8. Blandino, M.; Reyneri, A.; Vanara, F.; Tamietti, G.; Pietri, A. Influence of agricultural practices on Fusarium infection, fumonisin and deoxynivalenol contamination of maize kernels. World Mycotoxin J. 2009, 2, 409–418. [Google Scholar] [CrossRef]
  9. Pinotti, L.; Dell’Orto, V. Feed safety in the feed supply chain. Biotechnol. Agron. Soc. 2011, 15, 9–14. [Google Scholar]
  10. Castellá, G.; Bragulat, M.R.; Cabañes, F.J. Surveillance of fumonisins in maize-based feeds and cereals from Spain. J. Agric. Food Chem. 1999, 47, 4707–4710. [Google Scholar] [CrossRef]
  11. Alexa, E.; Pop, G.; Sumalan, R.M.; Radulov, I.; Poiana, M.; Tulcan, C. Fusarium species and fusarium mycotoxins in cereals from West Romania: Preliminary survey. Commun. Agric. Appl. Biol. Sci. 2011, 76, 661–666. [Google Scholar]
  12. Twarużek, M.; Skrzydlewski, P.; Kosicki, R.; Grajewski, J. Mycotoxins survey in feed materials and feedingstuffs in years 2015–2020. Toxicon 2021, 202, 27–39. [Google Scholar] [CrossRef]
  13. Pietri, A.; Bertuzzi, T.; Pallaroni, L.; Piva, G. Occurrence of mycotoxins and ergosterol in maize harvested over 5 years in Northern Italy. Food Addit. Contam. 2004, 21, 479–487. [Google Scholar] [CrossRef]
  14. Berardo, N.; Lanzanova, C.; Locatelli, S.; Laganà, P.; Verderio, A.; Motto, M. Levels of total fumonisins in maize samples from Italy during 2006–2008. Food Addit. Contam. B Surveill. 2011, 4, 116–124. [Google Scholar] [CrossRef]
  15. Weaver, A.C.; Weaver, D.M.; Adams, N.; Yiannikouris, A. Co-Occurrence of 35 Mycotoxins: A Seven-Year Survey of Corn Grain and Corn Silage in the United States. Toxins. 2021, 13, 516. [Google Scholar] [CrossRef]
  16. Raj, J.; Farkaš, H.; Jakovčević, Z.; Medina, A.; Magan, N.; Čepela, R.; Vasiljević, M. Comparison of multiple mycotoxins in harvested maize samples in three years (2018–2020) in four continents. Food Addit. Contam.—Chem. Anal. Control Expo. Risk Assess. 2022, 39, 599–608. [Google Scholar] [CrossRef]
  17. Lee, H.J.; Ryu, D. Worldwide Occurrence of Mycotoxins in Cereals and Cereal-Derived Food Products: Public Health Perspectives of Their Co-occurrence. J. Agr. Food Chem. 2017, 65, 7034–7051. [Google Scholar] [CrossRef]
  18. Feijó, C.J.A.; Orso, P.B.; Bordin, K.; Hara, R.V.; Luciano, F.B. Toxicological effects of fumonisin B1 in combination with other Fusarium toxins. Food Chem. Toxicol. 2018, 121, 483–494. [Google Scholar] [CrossRef]
  19. Marin, S.; Ramos, A.J.; Cano-Sancho, G.; Sanchis, V. Mycotoxins: Occurrence, toxicology, and exposure assessment. Food Chem. Toxicol. 2013, 60, 218–237. [Google Scholar] [CrossRef]
  20. Eskola, M.; Kos, G.; Elliott, C.T.; Hajšlová, J.; Mayar, S.; Krska, R. Worldwide contamination of food-crops with mycotoxins: Validity of the widely cited ‘FAO estimate’ of 25%. Crit. Rev. Food Sci. Nutr. 2020, 60, 2773–2789. [Google Scholar] [CrossRef]
  21. Khodaei, D.; Javanmardi, F.; Khaneghah, A.M. The global overview of the occurrence of mycotoxins in cereals: A three-year survey. Curr. Opin. Food Sci. 2021, 39, 36–42. [Google Scholar] [CrossRef]
  22. Streit, E.; Schatzmayr, G.; Tassis, P.; Tzika, E.; Marin, D.; Taranu, I.; Tabuc, C.; Nicolau, A.; Aprodu, I.; Puel, O.; et al. Current Situation of Mycotoxin Contamination and Co-occurrence in Animal Feed—Focus on Europe. Toxins 2012, 4, 788–809. [Google Scholar] [CrossRef]
  23. Tolosa, J.; Font, G.; Mañes, J.; Ferrer, E. Natural occurrence of emerging Fusarium mycotoxins in feed and fish from aquaculture. J. Agric. Food Chem. 2014, 62, 12462–12470. [Google Scholar] [CrossRef]
  24. Scarpino, V.; Vanara, F.; Sulyok, M.; Krska, R.; Blandino, M. Fate of regulated, masked, emerging mycotoxins and secondary fungal metabolites during different large-scale maize dry-milling processes. Food Res. Int. 2021, 140, 109861. [Google Scholar] [CrossRef]
  25. Tolosa, J.; Rodríguez-Carrasco, Y.; Ruiz, M.J.; Vila-Donat, P. Multi-mycotoxin occurrence in feed, metabolism and carry-over to animal-derived food products: A review. Food Chem. Toxicol. 2021, 158, 112661. [Google Scholar] [CrossRef]
  26. Mihalcea, A.; Amariei, S. Study on Contamination with Some Mycotoxins in Maize and Maize-Derived. Foods Appl. Sci. 2022, 12, 2579. [Google Scholar] [CrossRef]
  27. Alma, A.; Lessio, F.; Reyneri, A.; Blandino, M. Relationships between Ostrinia nubilalis (Lepidoptera: Crambidae) feeding activity, crop technique and mycotoxin contamination of corn kernel in northwestern Italy. Int. J. Pest Manag. 2005, 51, 165–173. [Google Scholar] [CrossRef]
  28. Camardo, L.M.; Bertuzzi, T.; Pietri, A.; Battilani, P. Mycotoxin occurrence in maize produced in Northern Italy over the years 2009–2011: Focus on the role of crop related factors. Phytopathol. Mediterr. 2015, 54, 212–221. [Google Scholar] [CrossRef]
  29. Blandino, M.; Scarpino, V.; Vanara, F.; Sulyok, M.; Krska, R.; Reyneri, A. The Role of the European Corn Borer (Ostrinia Nubilalis) on contamination of maize with 13 Fusarium mycotoxins. Food Addit. Contam. Part A 2015, 32, 533–543. [Google Scholar] [CrossRef] [PubMed]
  30. Gaikpa, D.S.; Miedaner, T. Genomics-assisted breeding for ear rot resistances and reduced mycotoxin contamination in maize: Methods, advances and prospects. Theor. Appl. Genet. 2019, 132, 2721–2739. [Google Scholar] [CrossRef] [PubMed]
  31. Piacentini, K.C.; Rocha, L.O.; Savi, G.D.; Carnielli-Queiroz, L.; Fontes, L.D.C.; Corrêa, B. Assessment of Toxigenic Fusarium Species and Their Mycotoxins in Brewing Barley Grains. Toxins 2019, 11, 31. [Google Scholar] [CrossRef] [PubMed]
  32. Blandino, M.; Scarpino, V.; Giordano, D.; Sulyok, M.; Krska, R.; Vanara, F.; Reyneri, A. Impact of Sowing Time, Hybrid and Environmental Conditions on the Contamination of Maize by Emerging Mycotoxins and Fungal Metabolites. Ital. J. Agron. 2017, 12, 215–224. [Google Scholar] [CrossRef]
  33. Okoth, S.; Rose, L.J.; Ouko, A.; Netshifhefhe, N.E.I.; Sila, H.; Viljoen, A. Assessing genotype-by-environment interactions in Aspergillus ear rot and pre-harvest aflatoxin accumulation in maize inbred lines. Agronomy 2017, 7, 86. [Google Scholar] [CrossRef]
  34. Giorni, P.; Magan, N.; Pietri, A.; Bertuzzi, T.; Battilani, P. Studies on Aspergillus section Flavi isolated in northern Italy from maize. Int. J. Food Microbiol. 2007, 113, 330–338. [Google Scholar] [CrossRef]
  35. Szabo, B.; Toth, B.; Toth, T.E.; Varga, M.; Kovacs, N.; Varga, J.; Kocsube, S.; Palagyi, A.; Bagi, F.; Budakov, D.; et al. A new concept to secure food safety standards against Fusarium species and Aspergillus flavus and their toxins in maize. Toxins 2018, 10, 372. [Google Scholar] [CrossRef]
  36. Neme, K.; Mohammed, A. Mycotoxin occurrence in grains and the role of postharvest management as a mitigation strategy. A review. Food Control 2017, 78, 412–425. [Google Scholar] [CrossRef]
  37. Munkvold, G.P. Mycotoxins in corn-occurrence, impact, and management. In Corn: Chemistry and Technology, 2nd ed.; White, P.J., Johnson, L.A., Eds.; AAC: St. Paul, MN, USA, 2003; pp. 811–881. [Google Scholar]
  38. CAST. Mycotoxins: Risks in Plant, Animal, and Human Systems, Task Force Report 139; Council for Agricultural Science and Technology: Ames, IA, USA, 2003. [Google Scholar]
  39. Commission Regulation. European Commission Regulation (EC) N 401/2006 of 23 February 2006. Laying down the methods of sampling and analysis for the official control of the levels of mycotoxins in foodstuffs. Off. J. Eur. Union 2006, L70, 12–34. [Google Scholar]
  40. Sulyok, M.; Berthiller, F.; Krska, R.; Schuhmacher, R. Development and validation of a liquid chromatography/tandem mass spectrometric method for the determination of 39 mycotoxins in wheat and maize. Rapid Commun. Mass Spectrom. 2006, 20, 2649–2659. [Google Scholar] [CrossRef]
  41. Malachova, A.; Sulyok, M.; Beltran, E.; Berthiller, F.; Krska, R. Optimization and validation of a quantitative liquid chromatography—Tandem mass spectrometric method covering 295 bacterial and fungal metabolites including all relevant mycotoxins in four model food matrices. J. Chromatogr. A 2014, 1362, 145–156. [Google Scholar] [CrossRef]
  42. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022; Available online: https://www.R-project.org/ (accessed on 28 June 2022).
Figure 1. Distribution of the Italian National Maize Monitoring Network in five macro-geographical areas: W, west (Torino) in red; C, center in green; SP, south of the Po River (blue); A, Adriatic (yellow); and E, east (Venezia) in purple.
Figure 1. Distribution of the Italian National Maize Monitoring Network in five macro-geographical areas: W, west (Torino) in red; C, center in green; SP, south of the Po River (blue); A, Adriatic (yellow); and E, east (Venezia) in purple.
Toxins 14 00520 g001
Figure 2. (A) Meteorological trend during the 2011–2021 period. (B) Meteorological trend in the different investigated geographical maize areas (averaged for 2011–2021 years).
Figure 2. (A) Meteorological trend during the 2011–2021 period. (B) Meteorological trend in the different investigated geographical maize areas (averaged for 2011–2021 years).
Toxins 14 00520 g002
Figure 3. Diffusion of regulated mycotoxins in maize lots from drying and storage centers located in different areas in Northern Italy. For each mycotoxin, the best year with the lowest incidence and contamination was compared with the worst year, with the highest one. Each mycotoxin was expressed as a percentage of incidence over different ranges of increasing contamination, pertaining to Regulations (EC) no. 1126/2007 and no. 165/2010 [3,4]: (A) diffusion of FBs; (B) diffusion of AF B1; (C) diffusion of DON; (D) diffusion of ZEA.
Figure 3. Diffusion of regulated mycotoxins in maize lots from drying and storage centers located in different areas in Northern Italy. For each mycotoxin, the best year with the lowest incidence and contamination was compared with the worst year, with the highest one. Each mycotoxin was expressed as a percentage of incidence over different ranges of increasing contamination, pertaining to Regulations (EC) no. 1126/2007 and no. 165/2010 [3,4]: (A) diffusion of FBs; (B) diffusion of AF B1; (C) diffusion of DON; (D) diffusion of ZEA.
Toxins 14 00520 g003
Figure 4. Contamination levels of FBs referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2011–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the FB contamination level is concerned).
Figure 4. Contamination levels of FBs referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2011–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the FB contamination level is concerned).
Toxins 14 00520 g004
Figure 5. Contamination levels of AF B1 referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2012–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the AF B1 contamination level is concerned).
Figure 5. Contamination levels of AF B1 referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2012–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the AF B1 contamination level is concerned).
Toxins 14 00520 g005
Figure 6. Contamination levels of DON referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2013–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the DON contamination level is concerned; b * = the areas A, W, SP and C were not statistically different from each other and they were assigned the letter b; *** = the areas were significantly different (p-value < 0.05) and they were assigned different letters as follows: area W→a; areas E and A→ab; areas C and SP→b).
Figure 6. Contamination levels of DON referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2013–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the DON contamination level is concerned; b * = the areas A, W, SP and C were not statistically different from each other and they were assigned the letter b; *** = the areas were significantly different (p-value < 0.05) and they were assigned different letters as follows: area W→a; areas E and A→ab; areas C and SP→b).
Toxins 14 00520 g006
Figure 7. Contamination levels of ZEA referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2014–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the ZEA contamination level is concerned).
Figure 7. Contamination levels of ZEA referring to the average toxin levels found in maize lots from drying and storage centers located in different areas in Northern Italy during the 2014–2021 period. Means followed by different letters are significantly different (p-value < 0.05), according to the REGW-F test, conducted separately for each year (ns = not significant, that is, the different areas were not statistically different from each other as far as the ZEA contamination level is concerned).
Toxins 14 00520 g007
Figure 8. Diffusion of regulated mycotoxins in maize lots from drying and storage centers located in the Po Valley. The worst year, that is, the one with the highest incidence and contamination, was expressed for each mycotoxin as a percentage of incidence over different ranges of increasing contamination, pertaining to Regulations (EC) no. 1126/2007 and no. 165/2010 [3,4]: (A) diffusion of FBs; (B) diffusion of AF B1; (C) diffusion of DON; (D) diffusion ZEA.
Figure 8. Diffusion of regulated mycotoxins in maize lots from drying and storage centers located in the Po Valley. The worst year, that is, the one with the highest incidence and contamination, was expressed for each mycotoxin as a percentage of incidence over different ranges of increasing contamination, pertaining to Regulations (EC) no. 1126/2007 and no. 165/2010 [3,4]: (A) diffusion of FBs; (B) diffusion of AF B1; (C) diffusion of DON; (D) diffusion ZEA.
Toxins 14 00520 g008
Figure 9. Diffusion of regulated mycotoxins in maize lots from drying and storage centers located in the Po Valley. The mean concentration, detected during the investigated years, was expressed for each mycotoxin as a percentage of incidence over different ranges of increasing contamination, pertaining to Regulations (EC) no. 1126/2007 and no. 165/2010 [3,4]: (A) diffusion of FBs; (B) diffusion of AF B1; (C) diffusion of DON; (D) diffusion of ZEA.
Figure 9. Diffusion of regulated mycotoxins in maize lots from drying and storage centers located in the Po Valley. The mean concentration, detected during the investigated years, was expressed for each mycotoxin as a percentage of incidence over different ranges of increasing contamination, pertaining to Regulations (EC) no. 1126/2007 and no. 165/2010 [3,4]: (A) diffusion of FBs; (B) diffusion of AF B1; (C) diffusion of DON; (D) diffusion of ZEA.
Toxins 14 00520 g009
Figure 10. Score plots of the considered maize lots (N = 3769) from drying and storage centers located in different areas in Northern Italy during the 2014–2021 period, marked according to the year (A) and area (B). The loadings (highlighted values > |0.5|) of each parameter (regulated mycotoxins and meteorological parameters) are reported with the two first principal components.
Figure 10. Score plots of the considered maize lots (N = 3769) from drying and storage centers located in different areas in Northern Italy during the 2014–2021 period, marked according to the year (A) and area (B). The loadings (highlighted values > |0.5|) of each parameter (regulated mycotoxins and meteorological parameters) are reported with the two first principal components.
Toxins 14 00520 g010
Figure 11. Score plots of the maize lots (N = 200) from drying and storage centers located in different areas in Northern Italy during the 2012–2015 period, marked according to the year (A) and area (B). The loadings (highlighted values > |0.5|) of each parameter (masked and emerging mycotoxins, other fungal metabolites, meteorological and insect parameters) are reported with the two first principal components.
Figure 11. Score plots of the maize lots (N = 200) from drying and storage centers located in different areas in Northern Italy during the 2012–2015 period, marked according to the year (A) and area (B). The loadings (highlighted values > |0.5|) of each parameter (masked and emerging mycotoxins, other fungal metabolites, meteorological and insect parameters) are reported with the two first principal components.
Toxins 14 00520 g011
Table 1. The occurrence of FBs and AF B1 in maize lots from drying and storage centers.
Table 1. The occurrence of FBs and AF B1 in maize lots from drying and storage centers.
FactorSource of VariationSamples
(No.)
FBsAF B1
Average (µg/kg)Min–Max (µg/kg)Pos * (%)Average (µg/kg)Min–Max (µg/kg)Pos * (%)
Year (Y)20114561886 g<LOD **–18,99097
20123626685 b<LOD **–47,3009712.5 a<LOD ***–11545
20133886447 b<LOD **—27,000994.4 cde<LOD ***—11922
20143555035 cd<LOD **—28,9201000.6 f<LOD ***–77
20153843815 ef<LOD **–38,260988.7 b<LOD ***–14548
20163144157 def<LOD **–23,220 997.6 b<LOD ***–14561
20172594212 de80–21,2701006.4 bc<LOD ***–9150
20183194209 de<LOD **–16,280992.7 ef<LOD ***–7836
20193078073 a320—23,6001006.2 bcd<LOD ***–6259
20203323224 f<LOD **–19,590983.1 def<LOD ***–7545
20212935330 c<LOD **–29,852974.0 cde<LOD ***–5134
p-value<0.001<0.001
Area (A)W793FBs/716AF B13507 d<LOD **–28,920991.3 c<LOD ***–4418
C990FBs/868AF B14820 b<LOD **–29,852995.7 b<LOD ***–11941
E793FBs/716AF B16577 a<LOD **–47,300985.3 b<LOD ***–14539
SP566FBs/472AF B14206 c<LOD **–24,360978.3 a<LOD ***–14552
A751FBs/673AF B14500 bc<LOD **–22,830999.2 a<LOD ***–11756
p-value-<0.001--<0.001--
Y × Ap-value-<0.001--<0.001--
Means followed by different letters, reported in superscript, are significantly different (the level of significance of the p-value is reported in the table), according to the REGW-F test. * Percentage of samples above the Limit of Detection (LOD). ** LOD of FBs = 25 µg/kg. *** LOD of AF B1 = 1 µg/kg.
Table 2. The occurrence of DON and ZEA in maize lots from drying and storage centers.
Table 2. The occurrence of DON and ZEA in maize lots from drying and storage centers.
FactorSource of VariationSamples
(No.)
DONZEA
Average (µg/kg)Min–Max (µg/kg)Pos * (%)Average (µg/kg)Min–Max (µg/kg)Pos * (%)
Year (Y)20133882049 b<LOD **–36,73194
20143554646 a50–36,583100364 a<LOD ***–139591
2015384221 d<LOD **–27407711 b<LOD ***–31531
2016314917 c<LOD **–17,2369415 b<LOD ***–25254
2017259264 d<LOD **–7623758 b<LOD ***–17525
2018319336 d<LOD **–38658210 b<LOD ***–15844
2019307426 cd<LOD **–90649323 b<LOD ***–32466
2020332166 d<LOD **–29038515 b<LOD ***–25842
202129350 d<LOD **–3756511 b<LOD ***–18536
p-value< 0.001< 0.001
Area (A)W655DON/568ZEA2232 a<LOD **–36,73197112 a<LOD ***–139570
C774DON/686ZEA895 b<LOD **–24,4078443 b<LOD ***–114245
E581DON/481ZEA1148 b<LOD **–26,0648895 a<LOD ***–118154
SP425DON/378ZEA424 c<LOD **–13,3257927 bc<LOD ***–84035
A516DON/450ZEA372 c<LOD **–13,9257221 c<LOD ***–128034
p-value<0.001<0.001
Y × Ap-value<0.001<0.001
Means followed by different letters, reported in superscript, are significantly different (the level of significance of the p-value is reported in the table), according to the REGW-F test. * Percentage of samples above the Limit of Detection (LOD). ** LOD of DON = 18.5 µg/kg. *** LOD of ZEA = 1.75 µg/kg.
Table 3. The occurrence of the masked mycotoxins DON-3-G, 3-ADON, and 15-ADON in maize lots from drying and storage centers.
Table 3. The occurrence of the masked mycotoxins DON-3-G, 3-ADON, and 15-ADON in maize lots from drying and storage centers.
FactorSource of VariationDON-3-G3-ADON15-ADON
Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)
Year (Y)2012133 c895 c10026 c100
2013595 b10087 a100290 b100
20141247 a10064 b91444 a94
201554 c602 c625 c23
p-value<0.001-<0.001-<0.001-
Area (A)W791 a10068 a80315 a89
C469 b8539 b75206 b73
E476 b7921 bc67155 bc64
SP33 c785 c5222 c70
A466 b7822 bc78124 bc87
p-value<0.001-<0.001-<0.001-
Y × Ap-value<0.001-<0.001-<0.001-
Means followed by different letters, reported in superscript, are significantly different (the level of significance of the p-value is reported in the table), according to the REGW-F test. * Percentage of samples above the LOD (DON-3-G = 0.8 µg/kg; 3-ADON = 1.2 µg/kg; 15-ADON = 12 µg/kg).
Table 4. The occurrence of the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. section Liseola (BEA, BIK, FA, FUS, and MON) in maize lots during the 2012–2015 period.
Table 4. The occurrence of the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. section Liseola (BEA, BIK, FA, FUS, and MON) in maize lots during the 2012–2015 period.
FactorSource of VariationBEABIKFAFUSMON
Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)
Year (Y)2012187 ab100294 b100356 bc100959 bc100852 a100
2013195 a100853 a1001236 a1001551 a100344 c100
2014135 bc100175 c100483 b981321 ab98505 b100
2015101 c100102 c100201 c87697 c79574 b100
p-value<0.001-<0.001-<0.001-0.003-<0.001-
Area (A)W189 a100306 b100492 bc971187 a100554 a100
C131 a100335 ab100604 ab961423 a98622 a100
E126 a100394 ab100588 abc97848 a88468 a100
SP144 a100277 b100344 c89835 a85686 a100
A144 a100470 a100825 a1001011 a87482 a100
p-value0.196-<0.001-<0.001-0.075-0.190-
Y × Ap-value0.136-<0.001-<0.001-0.040-0.011-
Means followed by different letters, reported in superscript, are significantly different (the level of significance of the p-value is reported in the table), according to the REGW-F test. * Percentage of samples above the LOD (BEA = 0.008 µg/kg; BIK = 8 µg/kg; FA = 16 µg/kg; FUS = 40 µg/kg; MON = 1.6 µg/kg).
Table 5. The occurrence of the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. section Discolor (AUR, BUT, CULM, and NIV) in maize lots from drying and storage centers.
Table 5. The occurrence of the emerging mycotoxins and fungal metabolites mainly produced by Fusarium spp. section Discolor (AUR, BUT, CULM, and NIV) in maize lots from drying and storage centers.
FactorSource of VariationAURBUTCULMNIV
Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)Average (µg/kg)Pos * (%)
Year (Y)2012161 c7741 b66109 c532 b100
20133929 b100567 a942438 a10021 b100
20149642 a100592 a100970 b10049 a89
2015180 c7589 b7558 c605 b11
p-value<0.001-<0.001-<0.001-<0.001-
Area (A)W4620 a98390 a1001595 a10015 a75
C4512 a88337 a851037 b7114 a75
E3213 ab85355 a79298 c7919 a70
SP481 c7437 b5642 c5235 a67
A2561 ab78393 a78241 c6528 a78
p-value<0.001-<0.001-<0.001-0.383-
Y × Ap-value0.002-0.005-<0.001-0.282-
Means followed by different letters, reported in superscript, are significantly different (the level of significance of the p-value is reported in the table), according to the REGW-F test. * Percentage of samples above the LOD (AUR = 2.4 µg/kg; BUT = 5.6 µg/kg; CULM = 8 µg/kg; NIV = 1.2 µg/kg).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Locatelli, S.; Scarpino, V.; Lanzanova, C.; Romano, E.; Reyneri, A. Multi-Mycotoxin Long-Term Monitoring Survey on North-Italian Maize over an 11-Year Period (2011–2021): The Co-Occurrence of Regulated, Masked and Emerging Mycotoxins and Fungal Metabolites. Toxins 2022, 14, 520. https://doi.org/10.3390/toxins14080520

AMA Style

Locatelli S, Scarpino V, Lanzanova C, Romano E, Reyneri A. Multi-Mycotoxin Long-Term Monitoring Survey on North-Italian Maize over an 11-Year Period (2011–2021): The Co-Occurrence of Regulated, Masked and Emerging Mycotoxins and Fungal Metabolites. Toxins. 2022; 14(8):520. https://doi.org/10.3390/toxins14080520

Chicago/Turabian Style

Locatelli, Sabrina, Valentina Scarpino, Chiara Lanzanova, Elio Romano, and Amedeo Reyneri. 2022. "Multi-Mycotoxin Long-Term Monitoring Survey on North-Italian Maize over an 11-Year Period (2011–2021): The Co-Occurrence of Regulated, Masked and Emerging Mycotoxins and Fungal Metabolites" Toxins 14, no. 8: 520. https://doi.org/10.3390/toxins14080520

APA Style

Locatelli, S., Scarpino, V., Lanzanova, C., Romano, E., & Reyneri, A. (2022). Multi-Mycotoxin Long-Term Monitoring Survey on North-Italian Maize over an 11-Year Period (2011–2021): The Co-Occurrence of Regulated, Masked and Emerging Mycotoxins and Fungal Metabolites. Toxins, 14(8), 520. https://doi.org/10.3390/toxins14080520

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