Next Article in Journal
AbobotulinumtoxinA (Dysport®), OnabotulinumtoxinA (Botox®), and IncobotulinumtoxinA (Xeomin®) Neurotoxin Content and Potential Implications for Duration of Response in Patients
Next Article in Special Issue
Safety of Corn and Corn-Based Products Intended for Human Consumption Concerning Fumonisins from a Brazilian Processing Plant
Previous Article in Journal
Biosynthetic Oligoclonal Antivenom (BOA) for Snakebite and Next-Generation Treatments for Snakebite Victims
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Post-Harvest Contamination with Mycotoxins in the Context of the Geographic and Agroclimatic Conditions in Romania

1
National Research & Development Institute for Food Bioresources—IBA Bucharest, 5 Baneasa Ancuta Street, 2nd district, Bucharest 020323, Romania
2
National Meteorological Administration (METEO—Romania), 97 Bucuresti‒Ploiesti Street, 1st district, Bucharest 013686, Romania
*
Author to whom correspondence should be addressed.
Submission received: 22 October 2018 / Revised: 7 December 2018 / Accepted: 8 December 2018 / Published: 13 December 2018
(This article belongs to the Special Issue Fungal Growth and Mycotoxins: Challenges for developing countries)

Abstract

:
This study aimed to assess post-harvest contamination with mycotoxins in the context of the geographic and agroclimatic conditions in Romania in 2012–2015, a period that was characterized by extreme meteorological events and the effects of climate change. The samples were randomly sampled from five agricultural regions of Romania and analyzed for mycotoxins by enzyme-linked immunosorbent assay. An SPSS analysis was done to explore correlations between mycotoxins (deoxynivalenol—DON, aflatoxins—AF, ochratoxin A—OTA, zearalenone—ZEA), product types (raw cereal, processed cereal, cereal-based food), geographic coordinates (latitude, longitude, agricultural region), and agroclimatic factors (air temperature, precipitation, soil moisture reserve, aridity index, soil type). In the southeast part of the Southern Plain and Dobrogea (Baragan Plain, located at 44–45° N, 26–27° E), contamination with AF and OTA was higher in raw and processed cereals (maize in silo, silo corn germs) in the dry years (2012 and 2013), and contamination with DON was high in processed cereals (wheat flour type 450) in the rainy year (2014). DON and OTA contamination were significantly correlated with cumulative precipitation in all years, while AF and ZEA contamination were non-significantly correlated with climatic factors and aridity indices. The distribution of mycotoxins by product type and the non-robust correlations between post-harvest mycotoxins and agrometeorological factors could be explained by the use of quality management systems that control cereal at warehouse receptions, performant processing technologies, and the quality of storage spaces of agri-food companies. The Baragan Plain is Romania’s most sensitive area to the predicted climate change in southeast Europe, which may be associated with its increased cereal contamination with AF and OTA.
Key Contribution: This approach correlates mycotoxin contamination with geographic and current/future agroclimatic factors in Romania and southeast Europe; moreover, it can be used by risk managers in agri-food companies and integrated into the official control systems by competent authorities.

1. Introduction

Cereal contamination occurs in the field under the influence of agroclimatic conditions that are favorable to fungal attack and mycotoxin production, and may continue post-harvest in the storage–processing–commercialization chain, thus resulting in financial losses and food safety and security issues [1,2,3,4]. Awareness of the effects of current and future agroclimatic conditions as well as the implementation of integrated mycotoxin management systems are issues of significant concern at both the national and global levels as a result of grain exports by the big traders.
Romania is located in southeast Europe between latitudes 43° N and 48° N and longitudes 20° E and 29° E, and has a humid temperate continental climate with hot summers in parts of the Southern country [5,6]. Due to its agroclimatic conditions, soil types, and agricultural technologies, Romania is a significant cereal producer at the global level. In the period of 1961–2016, it belonged to the first group for maize (8,800,672.82 t) and wheat (5,724,718.45 t), the second group for soybeans (180,864.95 t), the third group for triticale (104,945 t), and the fourth group for sorghum (15,539.54 t) [7]. From a capability point of view, arable soils in Romania are classified as very fertile (chernozem), fertile (phaeozem), or moderately fertile (luvisol) [8]. Since chernozem and phaeozem have high agroclimatic aridity indices [6,9], these soils require irrigation—especially in dry years—to sustain the fertility, productivity, and quality of crops. The agricultural area of Romania is divided into regions with different agroclimatic conditions due to relief forms and significant differences in the quantity and quality of cereal crops: Zone 1—Southern Plain, Oltenia Plain, and Western Plain (moderate temperate continental climate, chernozem); zone 2—Moldavian Plain and Dobrogea Plain (arid temperate continental climate, chernozem); and zone 3—Transylvania and the Southern Hilly Area (wet temperate continental climate, luvisol).
Failure to respect agricultural technologies and the existence of extreme meteorological events (high temperature, drought, heavy rain) adversely affects crop quantity and quality, favoring fungal attack and mycotoxin production. Until 2012, most Romanian studies reported the incidence and range of mycotoxins in cereals, food, and feed in separate regions, without an assessment of these factors in the geographic and agroclimatic contexts [10,11,12,13,14,15,16]. In-field contamination with deoxynivalenol was evaluated in wheat, rye, and triticale crops from 2012 to 2016, and in the wheat chain in 2015, taking into account the geographic positions and pedoclimatic conditions of the agricultural regions [17,18].
In Romania, the climate scenario is forecasted to involve a rise in temperature by 3 °C–5 °C in the north and 4 °C–5 °C in the south by 2100, with rainfall reduction and drought especially in the south, followed by a decline in crop yields of up to 20% [19,20,21,22]. Considering the current and future climate conditions, the socio-economic factors, and the political context as a European Union member country, it was predicted that Romania will present a medium climate risk index and an average vulnerability to climate change among the countries of southeast Europe [23,24].
The aim of this article was to assess post-harvest contamination with mycotoxins in the context of the geographic and agroclimatic conditions of Romania and climate change in southeast Europe. The research objectives were (a) to quantify mycotoxins in post-harvest samples in the period of 2012–2015, taking into account extreme meteorological events; (b) to statistically analyze post-harvest contamination with mycotoxins by year, product type, geographic coordinates, and agroclimatic data; and (c) to determine possible trends in mycotoxin contamination in Romania and southeast Europe in the context of climate change. This represents the first time that post-harvest contamination with mycotoxins has been investigated in Romania and southeast Europe in this context, and the approach used in this study can be applied in any country or even at the continental level along the agri-food chain. This approach was chosen because it correlates mycotoxin contamination with local agroclimatic factors, and because it can be used by risk managers in agri-food companies and integrated into the control systems of competent authorities.

2. Results

2.1. Agrometeorological Factors

In Romania, the growth of cereal crops depends on the value of agrometeorological factors during periods of maximum water consumption—May–June for winter wheat and June–August for maize. During the agricultural period of September 2011 to August 2015, the air temperature, precipitation, and soil moisture reserve varied between agricultural regions characterized by different relief forms. Tables S1–S6 provide the results of the ANOVA and Kruskal–Wallis nonparametric test of the physical context (agroclimatic factors: Air temperature, precipitation, water reserve in the soil) in which mycotoxin contamination occurred, also indicating grouping by years and locations.

2.1.1. Air Temperature

The average annual air temperature showed very significant differences (p < 0.001) by region and year (Table S1). The highest air temperatures were recorded in the warm and dry Southern regions (Oltenia Plain, Southern Plain and Dobrogea), followed by the humid and temperate regions (Southern Hilly Area, Western Plain) and colder regions (Moldavia, Transylvania) (Table S2). The year 2012 was the warmest and was significantly different from the other years. The next warmest year was 2015, while the differences in average air temperature between 2013 and 2014 were non-significant (Table S3).

2.1.2. Precipitation

The average cumulative precipitation in March–August showed non-significant differences between regions, but very significant differences between years (p < 0.001) (Tables S4 and S5). The highest precipitation was recorded in Southern Romania (Southern Hilly Area, Southern Plain, and Dobrogea) in 2014, and the lowest precipitation was also recorded in Southern Romania (Southern Plain and Dobrogea) in 2012 and 2015. The year, 2014, was the rainiest year, followed by 2013, and finally 2012 and 2015 between which the differences were non-significant (Table S5).

2.1.3. Soil Moisture Reserve

The average soil moisture reserve was not significantly difference between regions or years according to the ANOVA test.
However, after reanalyzing the agrometeorological factors with the Kruskal–Wallis non-parametric test for independent samples, significant differences were found for the distribution of temperature and cumulative precipitation between regions and the distribution of the soil moisture reserve between regions and years (Figure 1, Table S6).
The years, 2012 and 2015, were part of a series of extreme dry years that occurred in the 21st century in Romania (2002, 2003, 2007, 2012, and 2015), in which there were more than 60 days with maximum temperatures ≥32 °C, insufficient rainfall (≤200 mm), and a soil water deficit, especially in the southeast part of the country (Baragan Plain) during the period of June to August [25].

2.2. Post-Harvest Contamination with Mycotoxins

2.2.1. Deoxynivalenol (DON)

Post-harvest contamination with DON had an incidence rate of 42.9% (39/91) and varied between <18.5 and 1269.94 µg/kg (mean 82.39 ± 158.42 µg/kg, median <18.5 µg/kg); the incidence of samples above maximum level was 1.1% (1/91). The highest DON level was detected in a sample of wheat flour type 450 (processed cereals, in the year 2015) obtained from in-field contaminated wheat in the rainy year of 2014 after 11 months of storage (Figure 2A and Figure 3A). The incidence of DON-positive samples was high in the cold and wet areas of Transylvania (36.4%) and in the Southern Hilly Area (51.4%), as well as in the warm and arid Southern Plain and Dobrogea due to the abundant rainfall in 2014 (43.3%) (Figure 3A). The highest DON incidence was recorded in the year, 2014, which was very rainy in May–July (55.6%), followed by 2013 and 2012, which were rainy in May–June and very dry in July–August (41.2%, 38.5%) (Figure 3A). The processed cereal (wheat flour type 450) with DON above maximum level was sampled in May 2015 from the Southern Plain and Dobrogea (Chirnogi, Calarasi county) (Figure 4) and was possibly obtained from wheat contaminated with 1025 µg/kg DON in the rainy year of 2014 [26]. The contamination was caused by the abundant rainfall from May to July 2014 in the south; however, the region is not usually favorable to Fusarium attack and DON production due to its agroclimatic conditions (arid temperate continental climate, chernozem soil with high aridity indices). As such, only in-field cereals subjected to extreme weather conditions were affected [15,17,19,26].

2.2.2. Total aflatoxins (AF)

Post-harvest contamination with AF had an incidence rate of 45.4% (44/97), which varied between <1.75 and 82.94 µg/kg (mean 3.85 ± 14.80 µg/kg, median <1.75 µg/kg); the incidence of samples above the maximum level was 4.1% (4/97). The highest AF levels were detected in raw cereal (2012: Maize grains in silos—max. 82.94 µg/kg; 2013: Maize grains—max. 79.99 µg/kg) and processed cereals (2012: Maize extract—max. 4.53 µg/kg; 2013: Wheat flour type 1350—max. 8.34 µg/kg) (Figure 2B and Figure 3B). The detection of AF in maize grains in silos from October 2012 and August 2013 showed that contamination had occurred in-field in the dry years of 2012 and 2013. The incidence of AF was higher in the dry areas of Moldavia (12.5% positive samples; 12.5% samples above maximum level; max. 4.53 µg/kg) and in the Southern Plain and Dobrogea (8.4%; max. 82.94 µg/kg), while its incidence was sporadic in the cold and humid areas of Transylvania and in the Southern Hilly Area due to the local conditions in each year (Figure 2A,B). The highest incidence of AF above the maximum level was recorded in 2012 (7.7%; max. 82.94 µg/kg) and 2013 (7.3%; max. 79.99 µg/kg), which had low rainfall and high temperatures in the critical period of June–August. There was no sample above maximum level in 2014 and 2015 (Figure 2B and Figure 3B). The raw and processed cereals (maize grain in silos, silo-corn germs, maize extract) with AF above maximum levels were sampled from the Southern Plain and Dobrogea (Ianca, Braila county; Tandarei, Ialomita county) and Moldavia (Iasi, Iasi county) (Figure 4), where the agroclimatic conditions are favorable to Aspergillus sp. and AF production (arid temperate continental climate, chernozem soil with high aridity indices).

2.2.3. Ochratoxin A (OTA)

Post-harvest contamination with OTA had an incidence rate of 6.8% (4/59), which varied between <2.50 and 6.72 µg/kg (mean 2.70 ± 0.43 µg/kg, median <2.50 µg/kg); the incidence of samples above the maximum level was 6.8% (4/59). The highest OTA levels were detected in processed cereals (2012: Maize grain in silos—max. 3.39 µg/kg; 2014: Three samples of silo-corn germs—max. 5.64 µg/kg, 6.24 µg/kg, and 6.72 µg/kg) (Figure 3C). OTA only exceeded the maximum level in samples from the warm and arid Southern Plain and Dobrogea (10.0% positive samples; 10% samples above maximum level; max. 6.72 µg/kg) (Figure 2C). The highest incidence of OTA above maximum level was recorded in the dry year of 2012 (8.3%; max. 3.39 µg/kg) and the rainy year of 2014 (18.8% positive samples; 18.8% samples above maximum level; max. 6.72 µg/kg) (Figure 2C). However, the silo-corn germ samples from 2014 were obtained from corn contaminated in-field in the dry year of 2013 and stored for four months. The processed cereals (maize grains in silos; silo-corn germs) with OTA above the maximum level were sampled from the Southern Plain and Dobrogea (Faurei, Braila county; Tandarei, Ialomita county) (Figure 4) where, as previously mentioned, the agroclimatic conditions are favorable to Aspergillus sp. and OTA production.

2.2.4. Zearalenone (ZEA)

Post-harvest contamination with ZEA had an incidence rate of 7.1% (6/84), which varied between <1.75 and 7.05 µg/kg (mean 1.92 ± 0.22 µg/kg, median NA); the incidence of samples above the maximum level was 0% (0/84). The highest levels were detected in processed cereals in the year, 2014 (silo-corn germs—max. 4.65 µg/kg and 7.05 µg/kg), with both samples being from maize 2013 (Figure 2D). ZEA contamination was low in the Southern Hilly Area (5.6%; max. 2.94 µg/kg), which has a humid temperate continental climate, and the Southern Plain and Dobrogea (13.8%; max. 7.05 µg/kg), which have arid temperate continental climates (Figure 2A,B). No samples were positive for ZEA in Transylvania, Moldavia, or the Oltenia Plain regions. The highest incidence of ZEA was detected in the dry year of 2013 (6.5%; max. 3.38 µg/kg) and the rainy year of 2014 (23.5%; max. 7.05 µg/kg). No samples with ZEA incidence were found in the dry years of 2012 and 2015 (Figure 2D). The processed cereals (silo-corn germs, wheat flour type 1350) that were positive for ZEA were sampled from the Southern Plain and Dobrogea (Tandarei, Ialomita county; co-occurrence of ZEA and OTA) (Figure 4).

2.3. Statistical Analysis between Post-Harvest Mycotoxins and the Geographic and Agroclimatic Conditions in Romania

2.3.1. Comparison of the Mycotoxin Distribution by Product Type, Agricultural Region, and Year

The distribution of AF in raw cereals and DON, OTA, and ZEA in processed cereals was greater in the Southern Plain and Dobrogea under the influence of the agroclimatic conditions in 2012–2015 (Figure 2 and Figure 3). The mycotoxin distribution by product type revealed a very significant difference in DON (significance level = 0.000, N = 91) (Figure 2A) and a non-significant difference in distribution for AF (significance level = 0.616, N = 98) in processed cereals, even though the AF contamination was very high in post-harvest raw cereal (maize) from the Southern Plain and Dobrogea in the dry years of 2012 and 2013 (Figure 3B).

2.3.2. Dependence of Mycotoxin on the Agricultural Region—Correlations between Mycotoxins and Geographic Position

Because the statistical analysis of the differences between contamination levels by region revealed non-significant differences, the contamination with mycotoxins was considered to be local. Therefore, the analysis took geographic position into account, considering the locations of mills and storage units of the products (cereals, processed cereals, cereal-based foods). The graphical method was used to determine the distribution of each mycotoxin (panels on lines and columns, variables = Northern latitude, Eastern longitude, agricultural region, and year), the Baragan Plain was delimited as the area of interest between 44 and 45° N and 26 and 27° E (Figure 4). From this area, the linear correlations between mycotoxins, climatic factors, and aridity indices were determined through filtration by year as well as by the total study period in order to obtain a more substantial number of measurements.

2.3.3. Correlation between Mycotoxins and Meteorological Factors

DON was inversely correlated with the cumulative precipitation in all years (rxy = −0.296 *), and it was more robustly and directly correlated with cumulative precipitation in 2012 (rxy = 0.521 *) (Table 1).
AF and ZEA were not significantly correlated with meteorological factors or aridity indices (de Martonne aridity index, Iar-dM; climatic water deficit, CWD) either over the entire study period or in filtered years (Table 2).
When the study period was taken as a whole, OTA was distinctly and significantly directly correlated with cumulative precipitation (rxy = 0.416 **) and significantly inversely correlated with the mean air temperature (rxy = −0.317 *). In 2014, OTA was directly correlated with the cumulative rainfall by year (rxy = 0.904 *, N = 6) in the Southern Plain and Dobrogea for the Tandarei locality (Table 2).

2.3.4. Dependence of Mycotoxins on Soil Type

The dependence of mycotoxins on the soil type was analyzed by the non-parametric Jonckheere–Terpstra test for independent samples, which revealed the dependence between OTA and chernozem soils (significance level = 0.003, N = 59).
All of the analyzed situations showed non-significant correlations between mycotoxins and the water reserve in soil, and the aridity indices (Iar-dM, CWD) were predicted (Table 1 and Table 2).

3. Discussion

3.1. Mycotoxins in the Context of the Geographic and Agroclimatic Conditions of Romania

In Romania, post-harvest contaminations with AF, OTA, DON, and ZEA were detected in the Southern Plain and Dobrogea (Baragan Plain), under the extreme meteorological events in the period of 2012–2015. The Baragan Plain is located at latitude 44–45° N and longitude 26–27° E. It has calcic chernozem soil and a semi-arid temperate continental climate with sub-Mediterranean influences (multiannual average air temperature of 10–>11 °C, and a moderately dry and dry pluviometric regime <400 mm). This region is considered the second-most arid area after the Black Sea coast and is predicted to be the most vulnerable area to climate change by 2100 [6,25,27]. Due to its geographic position and semi-arid agroclimatic conditions, the Baragan Plain has optimal conditions for cereal contamination with Aspergillus sp. and AF [28,29]. However, it is also an important agricultural area of Romania, providing the highest corn production and grain storage capacity (755 authorised silos and warehouses, with a total capacity of 3,949,506 t), mainly through foreign traders.
In June–August of 2012–2015, the maximum air temperature exceeded 30 °C in the Southern Plain and Dobrogea (Baragan Plain), with the highest temperatures occurring in August (44 °C in 2012, 40 °C in 2013, 38 °C in 2014, and 41 °C in 2015). Hot summers in Southern parts place Romania in the group of Carpathian countries (Croatia, Hungary, Slovakia, Czech Republic, Poland, Ukraine, Romania, and Serbia) that are subject to heat wave events [30]. Among the cereals analyzed for the period of 2012–2015, maize was the most affected in terms of both production (−28% in the dry year of 2012, relative to the period of 1961–2016) and contamination with AF, followed by feed contamination with AF and cow’s milk contamination with AFM1. Despite the drought in 2012, Romania remained a significant cereal producer at the global level (first group for maize, second group for triticale and soybeans, third group for sorghum) [7].
In-field cereal contamination with DON is frequent in the Western and Northern regions where F. graminearum and F. culmorum are predominant, while it is sporadic and low in the Southern and Eastern areas [10,11,12,13,14,15,17,26,31]. Although there were reports of high DON contamination in grains during harvest in 2012, 2013, and 2014 [17,26], the present post-harvest study found low contamination, except in wheat flour type 450 obtained from wheat in 2014. The fact that 2014 was very rainy caused cereal control to be tightened at its reception at warehouses to prevent mycotoxin from entering the food chain. Post-harvest contamination with DON showed a very significant correlation with the processed cereals, and a significant correlation with the cumulative precipitation in all years of the study period. A very significant correlation of DON with meteorological factors (air temperature, precipitation, relative humidity) and regions was reported only for in-field contaminated cereals [32,33]. The use of modern cereal processing technologies by large companies in Romania can reduce contamination with DON up to 73–75%, even in wet sub-Carpathian regions where wheat mills are located [18]. Between 2012 and 2015, DON contamination was not detected in any cereal-based food sample.
In-field and post-harvest cereal contamination with AF is more frequent in the Southeastern part of the Southern Plain and Dobrogea (Baragan Plain), which is characterized by an arid temperate continental climate and chernozem soils that require irrigation to maintain high fertility. The agroclimatic conditions of the region vary annually and favor contamination with both Aspergillus sp. (A. flavus, A. fumigatus) and Fusarium sp. (F. graminearum, F. culmorum), with maize being the most contaminated grain [15]. The extreme weather events (high temperatures, drought) in Romania in 2012 and 2013 favored in-field maize contamination with AF and its transfer to cow’s milk as AFM1, which were noted in the Rapid Alert System for Food and Feed (RASFF) [34,35,36]. Post-harvest contamination with AF showed a non-significant correlation with meteorological factors and aridity indices, and a non-significant distribution in processed cereals, even if it was very high in post-harvest raw cereals (maize) from the Southern Plain and Dobrogea in the dry years of 2012 and 2013.
In-field and post-harvest contamination with OTA were reported only in the Southeastern part of the Southern Plain and Dobrogea (Baragan Plain) in co-existence with AF contamination in corn germs under the influence of extreme drought in 2012 and 2013. Post-harvest contamination with OTA showed an inversely significant correlation with the mean annual temperature, a significant correlation with the cumulative precipitation in all years, and, finally, a dependence on chernozem soils with high aridity indices. In Romania, cereal contamination with OTA has been reported in both the Southeastern region, which features a semi-arid temperate continental climate, as well as in the Western region, with has a semi-wet temperate continental climate [11,12,15].
Post-harvest contamination with ZEA was very low in the period of 2012–2015, but it was detected both in regions with a wet temperate continental climate (Transylvania, Southern Hilly Area) and regions with an arid temperate continental climate (Southern Plain and Dobrogea, Moldavia). Post-harvest contamination with ZEA was non-significantly correlated with agrometeorological factors (air temperature, precipitation, soil water reserve, aridity indices) when the analyzed years were taken as a whole. A significant correlation of ZEA with meteorological factors has been reported for in-field contaminated cereals [11,12,13,14,31,32].
The distribution of mycotoxins by product type was found to be influenced by the agricultural region and meteorological factors in Romania, in 2012–2015. The highest distribution of mycotoxins was recorded in the Southern Plain and Dobrogea, which both experienced drought in 2012 and 2013 and abundant rainfall in 2014. The existence of certified quality control systems in companies has limited mycotoxin contamination at the initial stages of the cereal chain (raw cereal, processed cereal), so that the cereal-based food produced is safe for consumption [18]. The non-significant correlation between post-harvest contamination with mycotoxins and agroclimatic factors can be explained by the control of cereal batches at the warehouse reception, the modernization of old or construction of modern storage spaces, and the implementation of high-performance processing technologies and quality management systems in agri-food companies. Romania has been a member of the European Union since January 2007, and must comply with the regulations for cereal storage and contamination of agri-food products with mycotoxins. We believe that these criteria have a greater contribution to the post-harvest non-significant correlation between mycotoxins and climatic factors, as compared with the ecosystem inside silos and the complex interactions between abiotic and biotic factors [37,38].

3.2. Mycotoxins in the Context of Current and Future Climate Conditions in Southeastern Europe

In 2012–2015, extreme weather events were recorded in Southeastern Europe (Romania, Serbia, Croatia, Macedonia, Albania) with drought in the years of 2012, 2013, and 2015, and abundant rainfall and floods in the year, 2014 [39,40]. Although the air temperature and precipitation showed significant differences between the years, this period was included in the warmest 15 years at the global level, compared to the reference period of 1961–1990 [40,41].
In 2012 and 2013, cereals contaminated with AF were redirected to animal feed, leading to the contamination of cow’s milk with AFM1 and RASFF notification of feed corn contaminated with AF. These incidences also led to numerous scientific publications [34,35,36,42,43,44,45,46,47]. In 2014, cereal contamination with DON was reported in the Northern part of Southeastern Europe, in Romania and Serbia [26,47]. In the period of 2012–2015, RASFF published a total of 1478 notifications on mycotoxins in food and feed, most frequently for AF (86.1%) in the dry years of 2012, 2013, and 2015, for OTA (11.2%) in the dry year of 2013 and the rainy year of 2014, and for both DON (1.4%) and fumonisin B (FB) (1.3%) in the rainy year of 2014 [34,48,49]. Cereal contamination with OTA in both dry and rainy years occurred because the fungus grows under different agroclimatic conditions; Aspergillus sp. is predominant in warm and temperate regions, and Penicillium sp. is predominant in colder areas [50,51].
In Southeastern Europe (Southern Romania, Northern Bulgaria, Northern Serbia, Eastern Croatia), the most fertile agricultural areas are plains with chernozem soils, on which wheat and maize crops are mainly grown. Because of their structure, chernozems are heavily affected by extreme weather events and climate change. In this European subregion, heat and drought events have intensified and led to a high-risk climate index (51–100), apart from Romania, which has a medium climate risk (21–50) and the highest aridization trend in Baragan Plain [6,23,24,30,52].
Knowing the current and future agroclimatic conditions is very important for the prevention and control of mycotoxins throughout the agri-food chain (in-field and post-harvest) because the climate change scenario in Southern Europe has involved an increase in the frequency, intensity, and duration of extreme weather events (drought, floods), as well as temperature increases of 1.8 °C–2.1 °C and a decrease in rainfall by 6% [24]. These climatic changes will influence fungal growth and mycotoxin production, depending on the ability of fungi to adapt [29,53,54,55,56]. Climate scenarios for the Balkans have shown that maize will be most affected by AF contamination in the +2 °C scenario; this situation is also valid for the Southern part of Romania [57,58].
In Europe, cereal contamination with mycotoxins varies depending on the geographic position and agroclimatic conditions. For example, AF, OTA, and FB prevail in the Southern region, which has a hot/wet Mediterranean climate, while DON and ZEA prevail in the Northern region, which has a colder climate [11,12,15,50,51,59,60,61,62,63,64,65].

4. Conclusions

This article provides important information on contamination with mycotoxins in the agri-food chain (post-harvest raw cereal, processed cereal, and cereal-based food) in the context of geographic and agroclimatic conditions of Romania in the period of 2012–2015, which involved extreme meteorological events. This included information regarding the current and future climate conditions of Southeastern Europe.
Post-harvest contamination with DON, AF, OTA, and ZEA mycotoxins was of a local character, with the Baragan Plain (southeast part of the Southern Plain and Dobrogea region) being the area of interest due to its agroclimatic conditions—a semi-arid temperate continental climate and chernozem soil with high aridity indices—which were strongly influenced by the annual meteorological events.
Among the cereals, maize was most affected by the extreme weather conditions, with the lowest productivity and high contamination with AF in the dry year of 2012, and high contamination with OTA in the dry year of 2013. DON and OTA were significantly correlated with the cumulative precipitation in all years, while AF and ZEA were not significantly correlated with climatic factors or aridity indices. The distribution of mycotoxins by product type and the non-robust correlations between post-harvest mycotoxins and agrometeorological factors could be explained by the use of quality management systems that control cereal production at warehouse receptions, performant processing technologies, and the high quality of storage spaces of the agri-food companies. Redirection of contaminated cereal batches to feed has led to RASFF notifications of contaminated feed and cow’s milk exports from Romania and other Southeastern European countries, which have caused financial and image crises in the economic and public health sectors.
Because this European subregion is subject to climate change (temperature rise, rainfall decrease), agri-food producers and competent authorities need to be aware of its impact on cereals and food quality and safety in order to develop an integrated management system to reduce the effects on human and animal health, as well as on the economy. This is the first time that post-harvest contamination (warehouse–processing–food) has been investigated in Romania or Southeastern Europe in general, and this approach can be applied in any country or even at the continental level along the agri-food chain.

5. Materials and Methods

5.1. Post-Harvest Samples and Product Types

Between January 2012 and July 2015, samples were randomly collected (N = 152; 1 kg/sample) from private cereal traders and food manufacturers with work points in five of the six agricultural regions of Romania (the Southern Plain and Dobrogea; the Oltenia Plain; the Southern Hilly Area; Moldavia; Transylvania; none in the Western Plain), which have certified quality systems according to the European legislation in force. Cereal-based food manufacturers also own mills from which they distribute processed cereals to their work points throughout the country. Systems of quality and distribution within and between companies form rigorous control processes and reduce the risk of mycotoxin contamination along the cereal chain.
To assess contamination with mycotoxins, the samples were grouped into three product types: Raw cereals (35 samples: Wheat, maize in silos, triticale, soybean, sorghum), processed cereals (69 samples: Wheat flour and bran; maize extract and starch, corn germs), and cereal-based foods (48 samples: Bread, cakes, pasta).

5.2. Mycotoxin Analysis

The samples were milled with a Retsch ZM 200 ultracentrifugal mill (Retsch, Haan, Germany) and analyzed for deoxynivalenol (DON), total aflatoxins (AF), ochratoxin A (OTA), and zearalenone (ZEA) by the enzyme-linked immunosorbent assay (ELISA) coupled with a Sunrise microplate reader (absorbance 450 nm; Tecan, Salzburg, Austria). Mycotoxin analysis was performed with Ridascreen® DON, Ridascreen® Aflatoxin Total, Ridascreen® Ochratoxin A 30/15, and Ridascreen® Zearalelon test kits (R-Biopharm, Darmstadt, Germany). The ELISA laboratory used in the current study is accredited for the analysis of mycotoxins in cereal and food samples according to ISO 17025:2005 (by the Romanian Accreditation Association (RENAR)) and participated in the FAPAS Proficiency Testing Programme (FERA, Sand Hutton, NY, UK) with satisfactory results, demonstrating laboratory performance and results comparable to laboratories worldwide that use different methods.
The 152 post-harvest samples were analyzed for the presence of mycotoxins (total of 331 tests: DON 91; AF 97; OTA 59; ZEA 84) depending on the product type and according to the Commission Regulations (EC) no. 1881/2006 and no. 1126/2007 [66,67], which set the maximum levels of contaminants as follows:
DON: Cereal—1250 µg/kg; processed cereal—750 µg/kg; food (bread, cakes)—500 µg/kg;
AF: Cereal—4 µg/kg; processed cereal—4 µg/kg; silage maize—10 µg/kg;
OTA: Cereal—5 µg/kg; processed cereal—3 µg/kg; and
ZEA: Cereal—100 µg/kg; processed cereal—75 µg/kg; food (bread, cakes)—50 µg/kg.

5.3. Geographic Coordinates

The Northern latitude and Eastern longitude of each county were determined using Google Earth [68], and counties were grouped into agricultural regions, based on their agroclimatic factors.

5.4. Agroclimatic Data

Agrometeorological factors (air temperature, precipitation, soil moisture reserve) were recorded by the official agrometeorological network with Meteorological Automatic Weather Station (MAWS) and Ceres-Wheat and Decision Support System for Agrotechnology Transfer (DSSAT) v.3.5. software between 1 September 2011 and 31 August 2015. The dominant soil types (chernozem, phaeozem, luvisol) in each county were set at a scale of 1:1,500,000 based on the Soil Atlas of Europe [8]. The aridity indices (de Martonne aridity index, Iar-dM; climatic water deficit, CWD) of each county were set based on data published for the period of 1900–2000 [6] to determine the correlation with mycotoxins on a long-term basis.

5.5. Data Collection

All data were collected in an Excel file and transferred into an SPSS database with the following variables: Years (2012–2015, with extreme meteorological events), mycotoxins (DON, AF, OTA, ZEA), post-harvest product types (raw cereal, processed cereal, cereal-based food), geographic coordinates (Northern latitude, Eastern longitude, agricultural region), and agroclimatic data (agrometeorological factors—air temperature, precipitation, and soil moisture reserve; soil types—chernozem, phaeozem, and luvisol; aridity indices—Iar-dM and CWD).

5.6. Statistical Analysis

The influences of the geographic position and agroclimatic conditions on post-harvest contamination with mycotoxins in Romania were determined by statistical analysis using SPSS v.23 software (IBM, Armonk, NY, USA) (Statistical Package for the Social Sciences software with ANOVA, general linear model, Pearson correlation, statistical tests for means comparison, non-parametric Kruskal–Wallis, and Jonckheere–Terpstra tests, graphical methods). The probability was considered to be statistically significant at p ≤ 0.05. The comparison of agrometeorological factors by region and year was made for the March–August period because there were differences in this period that influenced the growth of cereal crops; for example, the spring season began very early in the south, starting in February–March. The statistical analysis between agrometeorological factors and mycotoxins was performed for the critical periods specific to wheat and maize (May–June for DON; July–August for AF, OTA, and ZEA).

5.7. Replicability and Results of this Approach

The replicability and results regarding post-harvest contamination with mycotoxins are supported by data on the in-field contamination with mycotoxins (DON, AF), as well as the physico-chemical quality and productivity of annual crops (wheat, rye, triticale) under the influence of geographic position and agroclimatic conditions in Romania, between 2012 and 2016. All data for the in-field cereals (≈3600 samples) were statistically analyzed by SPSS v.23 software and evaluated in the context of the climate conditions for Romania and climate change in Southeastern Europe by the year, 2100. It is anticipated that these data will contribute to the understanding of mycotoxin occurrence in the geo-agroclimatic context of this European subregion (manuscripts are being drafted).

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/10/12/533/s1. Table S1. Differences in the average air temperature by region and year in Romania in 2012–2015 (ANOVA two-factors without replications). Table S2. Class divisions of the average air temperature by region in Romania in 2012–2015 (t-test multiple). Table S3. Class division of the average air temperature by year in Romania, 2012–2015 (t-test multiple). Table S4. Differences in the average cumulative precipitation by region and year in Romania in 2012–2015 (ANOVA two-factors without replications). Table S5. Class divisions of the average cumulative precipitation by year in Romania in 2012–2015 (t-test multiple). Table S6. Comparison of the average values of agrometeorological factors between agricultural regions and years in Romania in 2012–2015 (Kruskal–Wallis non-parametric test for independent samples, when the null hypothesis was rejected).

Author Contributions

Conceptualization, V.G. (as projects’ coordinator), I.A. and E.M.; Methodology and Formal analysis, V.G., I.A. and E.M.; Investigation, A.A.D., I.S. and M.E.C. (mycotoxin analysis) and O.A.O. (agrometeorological data); Data collection, V.G. (Excell file: mycotoxins, agrometeorological data, geographic coordinates, soil types, aridity indices), E.M. (World file: agrometeorological data) and I.A. (SPSS database); Software, I.A. (SPSS analysis); Validation, V.G., I.A. and E.M.; Projects administration and Supervision, V.G. and E.M.; Supervision, Projects administration within the institutional Nucleu and ADER 2020 Programmes, Support of publishing, E.I. and N.B.; Writing—Original Draft, Writing—review & editing, V.G. and I.A.

Funding

This research was funded through Nucleu Programme (Grants no. PN 12 48 01 05, PN 16 46 04 01 and PN 18 02 01 03) by the Romanian Ministry of Research and Innovation, and the ADER 2020 Programme (Grant no. ADER 8.1.1) by the Ministry of Agriculture and Rural Development.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Paterson, R.R.M.; Lima, N. How will climate change affect mycotoxins in food? Food Res. Int. 2010, 43, 1902–1914. [Google Scholar] [CrossRef] [Green Version]
  2. Paterson, R.R.M.; Lima, N.; Taniwaki, M.H. Coffee, mycotoxins and climate change. Food Res. Int. 2014, 61, 1–15. [Google Scholar] [CrossRef]
  3. Wild, C.P.; Gong, Y.Y. Mycotoxins and human disease: a largely ignored global health issue. Carginogenesis 2010, 31, 1–82. [Google Scholar] [CrossRef] [PubMed]
  4. Tirado, M.C.; Clarke, R.; Jaykus, L.A.; McQuatters-Gollop, A.; Frank, J.M. Climate change and food safety: A review. Food Res. Int. 2010, 43, 1745–1765. [Google Scholar] [CrossRef]
  5. Chen, D.; Chen, H.W. Using the Köppen classification to quantify climate variation and change: An example for 1901-2010. Environ. Dev. 2013, 6, 69–79. [Google Scholar] [CrossRef]
  6. Paltineanu, C.; Mihailescu, I.F.; Seceleanu, I.; Dragota, C.; Vasenciuc, F. Using aridity indices to describe some climate and soil features in eastern Europe: a Romanian case study. Theor. Appl. Clim. 2007, 90, 263–274. [Google Scholar] [CrossRef]
  7. Food and Agriculture Organization (FAO); FAOSTAT. Compare Data. Available online: http://www.fao.org/faostat/en/#compare (accessed on 18 July 2018).
  8. European Commission (EC); Soil Atlas of Europe, European Soil Bureau Network. Office for Official Publications of the European Communities, L-2995 Luxembourg. 2005. Available online: https://esdac.jrc.ec.europa.eu/Projects/Soil_Atlas/Download/Atlas.pdf (accessed on 10 January 2017).
  9. Croitoru, A.E.; Piticar, A.; Imbroane, A.M.; Burada, D.C. Spatiotemporal distribution of aridity indices based on temperature and precipitation in the extra-Carpathian regions of Romania. Theor. Appl. Clim. 2012, 112, 597–607. [Google Scholar] [CrossRef]
  10. Misca, C.; Daminescu, L.; Jianu, C.; David, I.; Misca, L.; Marginean, O.; Radoi, B.; Rinovetz, A.; Bujanca, G.; Velciov, A.; et al. The incidence of the strains of Fusarium sp. and of Zearalenone in cereals analyzed from the south west of Romania. Ann. West Univ. Timisoara Ser. Biol. 2014, XVII, 137–144. [Google Scholar]
  11. Alexa, E.; Dehelean, C.A.; Poiana, M.A.; Radulov, I.; Cimpean, A.M.; Bordean, D.M.; Tulcan, C.; Pop, G. The occurrence of mycotoxins in wheat from western Romania and histopathological impact as the effect of feed intake. Chem. Cent. J. 2013, 7, 2–11. [Google Scholar] [CrossRef] [PubMed]
  12. Alexa, E.; Pop, G.; Sumalan, R.; Radulov, I.; Poiana, M.; Tulcan, C. Fusarium species and Fusarium mycotoxins in cereals from West Romania: A preliminary survey. Commun. Agric. Appl. Biol. Sci. 2011, 76, 661–666. Available online: https://0-www-ncbi-nlm-nih-gov.brum.beds.ac.uk/pubmed/22702185 (accessed on 20 January 2017). [PubMed]
  13. Curtui, V.; Usleber, E.; Dietrich, R.; Lepschy, J.; Märtlbauer, E. A survey on the occurrence of mycotoxins in wheat and maize from western Romania. Mycopathologia 1998, 143, 97–103. [Google Scholar] [CrossRef] [PubMed]
  14. Banu, I.; Aprodu, I.; Nicolau, A. The occurrence of Fusarium Mycotoxins (Deoxynivalenol and Zearalenone) in Wheat and High Fibre Wheat Bread in Eastern Romania. J. Environ. Prot. Ecol. 2011, 2, 519–525. [Google Scholar]
  15. Tabuc, C.; Taranu, I.; Calin, L. Survey of mould and mycotoxin contamination of cereals in south-eastern Romania in 2008–2010. Arch. Zootech. 2011, 14, 25–38. [Google Scholar]
  16. Stroia, C.; Tabuc, C.; Neacsu, A. Incidence of Fusarium species and its mycotoxins in cereals from western Romania. Res. J. Agric. Sci. 2010, 42, 302–309. [Google Scholar]
  17. Gagiu, V.; Mateescu, E.; Smeu, I.; Dobre, A.A.; Cucu, M.E.; Oprea, O.A.; Iorga, E.; Belc, N. A survey of the cereal contamination with deoxynivalenol in Romania, for 2011-2013 period. Romanian Biotechnol. Lett. 2017, 22, 12240–12249. Available online: https://www.rombio.eu/vol22nr1/----11_V.Gagiu_Final__17.02.2017_.pdf (accessed on 1 March 2017).
  18. Gagiu, V.; Doja, L.; Mateescu, E.; Smeu, I.; Cucu, M.E.; Oprea, O.; Iorga, E.; Belc, N. Contamination with deoxynivalenol in the milling-bakery industry under the influence of climatic conditions from Romania. J. Hyg. Eng. Des. 2016, 16, 38–44. Available online: http://www.jhed.mk/filemanager/JHED%20Vol.%2016/02.%20FQ&S/05.%20Full%20paper%20-%20Valeria%20Gagiu.pdf (accessed on 3 January 2017).
  19. Sandu, I.; Mateescu, E. Schimbari Climatice Actuale Si Previzibile in Romania. Schimbari Climatice—Provocare Majora Pentru Cercetarea Agricola [Current and predictable climate change in Romania. Climate change—A Major Challenge for Agricultural Research]; Academia Romana—Sectia de Stiinte Agricole si Silvice: Editura Academiei Romane, Romania, 2014; pp. 17–36. ISBN 978-973-27-2421-7. [Google Scholar]
  20. Busuioc, A.; Caian, M.; Bojariu, R.; Boroneant, C.; Baciu, M.; Dumitrescu, A. Scenarii de Schimbare a Regimului Climatic in Romania Pe Perioada 2001-2030 [Scenarios for Changing the Climate Regime in Romania for the Period 2001-2030]; Administratia Nationala de Meteorologie: Bucharest, Romania, 2012; pp. 1–25. Available online: http://mmediu.ro/new/wp-content/uploads/2014/02/2012-04-23_schimbari_climatice_schimbareregimclimatic2001_2030.pdf (accessed on 8 January 2018).
  21. IPCC. Climate Change 2007: Synthesis Report. Contribution of Working Groups I, II and III to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Pachauri, R.K., Reisinger, A., Eds.; IPCC: Geneva, Switzerland, 2007; 104p, Available online: http://www.ipcc.ch/pdf/assessment-report/ar4/syr/ar4_syr_full_report.pdf (accessed on 8 January 2018).
  22. World Bank. World Development Report 2010: Development and Climate Change. p. 439. Available online: https://openknowledge.worldbank.org/handle/10986/4387 (accessed on 8 January 2018).
  23. Climate Risk Index (CRI). Global Climate Risk Index 1997–2016. Available online: https://germanwatch.org/de/14638 (accessed on 8 January 2018).
  24. Fay, M.; Block, R.I.; Ebinger, J. Adapting to Climate Change in Eastern Europe and Central Asia. Europe and Central Asia Reports; World Bank: Washington, DC, USA, 2010; Available online: http://hdl.handle.net/10986/2407 (accessed on 8 January 2018).
  25. Oprea, O.A.; Mateescu, E.; Barbu, A.; Tudor, R. Extreme dry years in the 21st century at the level of the agricultural areas of Muntenia, Romania. Sciendo 2018, 1, 101–109. [Google Scholar] [CrossRef]
  26. Gagiu, V.; Cucu, M.; Dobre, A.; Mateescu, E.; Oprea, O.; Belc, N. Contaminarea cu micotoxina deoxinivalenol a recoltei de grâne din România, în anul agricol 2013–2014. In Contamination of Grain Crops with Deoxynivalenol in Romania, in Agricultural Year 2013-2014; Editura Inovativ Media: București, Romania, 2014; ISBN 978-973-0-18158-6. [Google Scholar]
  27. Vrînceanu, P.G.A. Research on the Vulnerability of Soils to the Degradation Processes—Desertification in the Bărăgan Plain. Ph.D. Thesis, USAMV Bucharest, Faculty of Land Improvements and Environment Engineering, București, Romania, 2011; 256p. Available online: https://www.usamv.ro/images/Programe_de_studii/Doctorat/Teze_de_doctorat/Arhiva_2011/Vrinceanu_Rez.pdf (accessed on 3 August 2018).
  28. Cotty, P.J.; Jaime-Garcia, R. Influence of climate on aflatoxin producing fungi and aflatoxin contamination. Int. J. Food Microbiol. 2007, 119, 109–115. [Google Scholar] [CrossRef]
  29. van der Fels-Klerx, H.J.; Liu, C.; Battilani, P. Modelling climate change impacts on mycotoxin contamination. Special Issue: Mycotoxins in a changing world. World Mycotoxin J. 2016, 9, 717–726. [Google Scholar] [CrossRef]
  30. Spinoni, J.; Lakatos, M.; Szentimrey, T.; Bihari, Z.; Szalai, S.; Vogt, J.; Antofie, T. Heat and cold waves trends in the Carpathian Region from 1961 to 2010. Int. J. Climatol. 2015, 35, 4197–4209. [Google Scholar] [CrossRef] [Green Version]
  31. Tabuc, C.; Marin, D.; Guerre, P.; Susan, T.; Bailly, J.D. Molds and Mycotoxin Content of Cereals in Southeastern Romania. J. Food Prot. 2009, 72, 662–665. Available online: http://jfoodprotection.org/doi/pdf/10.4315/0362-028X-72.3.662?code=FOPR-site (accessed on 18 April 2017). [CrossRef] [PubMed]
  32. Stanciu, O.; Juan, C.; Miere, D.; Loghin, F.; Maňes, J. Occurrence and co-occurrence of Fusarium mycotoxins in wheat grains and wheat flour from Romania. Food Control 2017, 73, 147–155. [Google Scholar] [CrossRef]
  33. Váňová, M.; Klem, K.; Matušinský, P.; Trnka, M. Prediction Model for Deoxynivalenol in Wheat Grain Based on Weather Conditions. Plant Prot. Sci. 2009, 45. Available online: https://www.agriculturejournals.cz/publicFiles/13965.pdf (accessed on 18 April 2018). [CrossRef]
  34. European Commission (EC). The Rapid Alert System for Food and Feed 2013 Annual Report. Available online: https://ec.europa.eu/food/sites/food/files/safety/docs/rasff_annual_report_2013.pdf (accessed on 3 August 2018).
  35. van Egmond, H.P. Foreword. World Mycotoxin J. 2018, 11, 1–4. [Google Scholar] [CrossRef]
  36. De Rijk, T.C.; van Egmond, H.P.; van der Fels-Klerx, H.J.; Herbes, R.; de Nijs, M.N.; Samson, R.A.; Slate, A.B.; van der Spiegel, M. A study of the 2013 Western European issue of aflatoxin contamination of maize from the Balkan area. World Mycotoxin J. 2015, 8, 641–651. [Google Scholar] [CrossRef]
  37. Dall’Asta, C.; Battilani, P. Fumonisins and their modified forms, a matter of concern in future scenario? Special Issue: Mycotoxins in a changing world. World Mycotoxin J. 2016, 9, 727–739. [Google Scholar] [CrossRef]
  38. Chakraborty, S.; Newton, A.C. Climate change, plant diseases and food security: an overview. Plant Pathol. 2011, 60, 2–14. [Google Scholar] [CrossRef] [Green Version]
  39. Spinoni, J.; Naumann, G.; Vogt, J.; Barbosa, P. European drought climatologies and trends based on a multi-indicator approach. Glob. Planet. Chang. 2015, 127, 50–57. [Google Scholar] [CrossRef]
  40. Sandu, I. Schimbari Climatice Observate Si Viitoare [Climate Changes Observed and Future]; Oral Presentation at World Meteorological Day with the theme. Climate Knowledge to Prevent Climate Change; ASAS: Bucharest, Romania, 2015; Available online: http://www.asas.ro/wcmqs/sectii/stiinta-solului/documente/ZMM%20ASAS%20FINAL.pdf (accessed on 4 March 2018).
  41. Mateescu, E. The Romanian Agrometeorological Services and Products—Current Status and Challenges in the Context of Climate Change; Workshop, Agrometeorologists for farmers in hotter, drier, wetter future 9–10 November 2016, Ljubljana, Slovenia; METEO: București, Romania, 2016; Available online: http://www.dmcsee.org/uploads/file/405_6_elenamateescu.pdf (accessed on 5 August 2018).
  42. Pleadin, J.; Vulić, A.; Perši, N.; Škrivanko, M.; Capek, B.; Cvetnić, Ž. Annual and regional variations of aflatoxin B1 levels seen in grains and feed coming from Croatian dairy farms over a 5-year period. Food Control 2015, 47, 221–225. [Google Scholar] [CrossRef]
  43. Bilandžić, N.; Božic, Ð.; Ðokic, M.; Sedak, M.; Solomun Kolanović, B.; Varenina, I.; Sanin Tanković, S.; Cvetnić, Z. Seasonal effect on aflatoxin M1 contamination in raw and UHT milk from Croatia. Food Control 2014, 40, 260–264. [Google Scholar] [CrossRef]
  44. Škrbic, 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]
  45. Kos, J.; Mastilović, J.; Hajnal, E.J.; Šarić, B. Natural occurrence of aflatoxins in maize harvested in Serbia during 2009–2012. Food Control 2013, 34, 31–34. [Google Scholar] [CrossRef]
  46. Dimitrieska-Stojković, E.; Stojanovska-Dimzoska, B.; Ilievska, G.; Uzunov, R.; Stojković, G.; Hajrulai-Musliu, Z.; Jankuloski, D. Assessment of aflatoxin contamination in raw milk and feed in Macedonia during 2013. Food Control 2016, 59, 201–206. [Google Scholar] [CrossRef]
  47. Krnjaja, V.; Tomić, Z.; Stanković, S.; Petrović, T.; Bijelić, Z.; Mandić, V.; Obradović, A. Fusarium infection and deoxynivalenol contamination in winter wheat. Biotechnol. Anim. Husb. 2015, 31, 123–131. [Google Scholar] [CrossRef]
  48. European Commission (EC). The Rapid Alert System for Food and Feed 2012 Annual Report. Available online: https://ec.europa.eu/food/sites/food/files/safety/docs/rasff_annual_report_2012_en.pdf (accessed on 3 August 2018).
  49. European Commission (EC). The Rapid Alert System for Food and Feed 2015 Annual Report. Available online: https://ec.europa.eu/food/sites/food/files/safety/docs/rasff_annual_report_2015.pdf (accessed on 3 August 2018).
  50. Gagiu, V. Triticale crop and contamination with mycotoxins under the influence of climate change—Global study. J. Hyg. Eng. Des. 2018, 23, 30–45. Available online: http://jhed.mk/filemanager/JHED%20Vol.%2023/02.%20FQ&S/04.%20Full%20paper%20-%20Valeria%20Gagiu.pdf (accessed on 3 January 2017).
  51. Schatzmayr, G.; Streit, E. Global occurrence of mycotoxins in the food and feed chain: Facts and figures. World Mycotoxin J. 2013, 6, 213–222. [Google Scholar] [CrossRef]
  52. Bojariu, R.; Velea, L.; Cica, R.D.; Dobrinescu, A.E.; Bîrsan, M.; Dumitrescu, A. Mechanisms of drought persistence in Romania in the climate change perspective. Ann. Univ. Craiova Biol. Hortic. Food Prod. Process. Technol. Environ. Eng. 2012, XVII, 523–526. Available online: http://www.euro4m.eu/Publications/090_BojariuR-craiova02.pdf (accessed on 19 March 2018).
  53. Medina, Á.; Rodríguez, A.; Sultan, Y.; Magan, N. Climate change factors and Aspergillus flavus: effects on gene expression, growth and aflatoxin production. World Mycotoxin J. 2015, 8, 171–179. [Google Scholar] [CrossRef]
  54. Medina, Á.; Gilbert, M.K.; Mack, B.M.; Obrian, G.R.; Rodríguez, A.; Bhatnagar, D.; Payne, G.; Magan, N. Interactions between water activity and temperature on the Aspergillus flavus transcriptome and aflatoxin B1 production. Int. J. Food Microbiol. 2017, 256, 36–44. [Google Scholar] [CrossRef]
  55. Akbar, A.; Medina, A.; Magan, N. Impact of interacting climate change factors on growth and ochratoxin A production by Aspergillus section Circumdati and Nigri species on coffee. World Mycotoxin J. 2016, 9, 863–874. [Google Scholar] [CrossRef]
  56. Wu, F.; Bhatnagar, D.; Bui-Klimke, T.; Carbone, I.; Hellmich, R.L.; Munkvold, G.P.; Paul, P.; Payne, G.; Tackle, E.S. Climate change impacts on mycotoxin risks in US maize. World Mycotoxin J. 2011, 4, 79–93. [Google Scholar] [CrossRef] [Green Version]
  57. Battilani, P.; Rossi, V.; Giorni, P.; Pietri, A.; Gualla, A.; van der Fels-Klerx, H.J.; Booij, C.J.H.; Moretti, A.; Logrieco, A.; Miglietta, F.; et al. Modelling and mapping the emergence of aflatoxins in cereals in the EU due to climate change. Sci. Rep. Submitt. EFSA 2012, 9, 223E. [Google Scholar] [CrossRef]
  58. Battilani, P.; Toscano, P.; van der Fels-Klerx, H.J.; Moretti, A.; Camardo Leggieri, M.; Brera, C.; Rortais, A.; Goumperis, T.; Robinson, T. Aflatoxin B1 contamination in maize in Europe increases due to climate change. Sci. Rep. 2016, 6, 24328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Demirel, R.; Yilmaz, S.N. Mycotoxigenic molds and mycotoxins in flours consumed in Turkey. J. Sci. Food Agric. 2014, 94, 1577–1584. [Google Scholar] [CrossRef] [PubMed]
  60. Aydin, A.; Aksu, H.; Gunsen, U. Mycotoxin levels and incidence of mold in Turkish rice. Environ. Monit. Assess 2011, 178, 271–280. [Google Scholar] [CrossRef] [PubMed]
  61. Kabak, B. Ochratoxin A in cereal-derived products in Turkey: Occurrence and exposure assessment. Food Chem. Toxicol. 2008, 47, 348–352. [Google Scholar] [CrossRef] [PubMed]
  62. Jurjevic, Z.; Solfrizzo, M.; Cvjetkovic, B.; De Girolamo, A.; Visconti, A. The occurrence of Beauvericin in Corn from Croatia. Food Technol. Biotechnol. 2002, 40, 91–94. Available online: http://www.ftb.com.hr/images/pdfarticles/2002/April-June/40-91.pdf (accessed on 18 April 2017).
  63. Kirinčič, S.; Škrjanc, B.; Kos, N.; Kozolc, B.; Pirnat, N.; Tavčar-Kalcher, G. Mycotoxins in cereals and cereal products in Slovenia—Official control of foods in the years 2008–2012. Food Control 2015, 50, 157–165. [Google Scholar] [CrossRef]
  64. Rafai, P.; Bata, A.; Jakab, L.; Vanyi, A. Evaluation of mycotoxin-contaminated cereals for their use in animal feeds in Hungary. Food Addit. Contam. 2000, 17, 799–808. [Google Scholar] [CrossRef]
  65. Bolechová, M.; Benešová, K.; Bélachová, S.; Čáslavský, J.; Pospíchalová, M.; Mikulíková, R. Determination of seventeen mycotoxins in barley and malt in the Czech Republic. Food Control 2015, 47, 108–113. [Google Scholar] [CrossRef]
  66. Commission Regulation (EC) No 1881/2006 of 19 December 2006 setting maximum levels for certain contaminants in foodstuffs. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32006R1881&from=EN (accessed on 1 November 2017).
  67. Commission Regulation (EC) No 1126/2007 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. Available online: https://eur-lex.europa.eu/legal-content/EN/TXT/PDF/?uri=CELEX:32007R1126&from=EN (accessed on 1 November 2017).
  68. Google Earth. Google Landsat/Copernicus. Available online: https://earth.google.com/web/ (accessed on 21 March 2017).
Figure 1. Comparisons of the distribution of agrometeorological factors by agricultural region and year: Air temperature (A), cumulative precipitation (B), soil moisture reserve (C,D) (Kruskal–Wallis non-parametric test for independent samples). °, °°, *°, ****°, *°°°°°—position of the extreme individual values.
Figure 1. Comparisons of the distribution of agrometeorological factors by agricultural region and year: Air temperature (A), cumulative precipitation (B), soil moisture reserve (C,D) (Kruskal–Wallis non-parametric test for independent samples). °, °°, *°, ****°, *°°°°°—position of the extreme individual values.
Toxins 10 00533 g001aToxins 10 00533 g001b
Figure 2. Distribution of deoxynivalenol (A), total aflatoxins (B), ochratoxin A (C), and zearalenone (D) by product type, agricultural region, and year in Romania in 2012–2015 (graphical method).
Figure 2. Distribution of deoxynivalenol (A), total aflatoxins (B), ochratoxin A (C), and zearalenone (D) by product type, agricultural region, and year in Romania in 2012–2015 (graphical method).
Toxins 10 00533 g002aToxins 10 00533 g002b
Figure 3. Distribution of deoxynivalenol (A) and total aflatoxins (B) by product type in Romania in 2012–2015 (Kruskal–Wallis non-parametric test for independent samples). * and **—position of the extreme individual values.
Figure 3. Distribution of deoxynivalenol (A) and total aflatoxins (B) by product type in Romania in 2012–2015 (Kruskal–Wallis non-parametric test for independent samples). * and **—position of the extreme individual values.
Toxins 10 00533 g003
Figure 4. The geographic origin of post-harvest samples with mycotoxins above the maximum level in Romania, in 2012–2015. DON—Calarasi county; AF—Braila, Ilfov, Valcea, and Iasi counties; OTA—Ialomita and Braila counties. Red shape: Baragan Plain (44–45° N, 26–27° E) in the Southern Plain and Dobrogea region.
Figure 4. The geographic origin of post-harvest samples with mycotoxins above the maximum level in Romania, in 2012–2015. DON—Calarasi county; AF—Braila, Ilfov, Valcea, and Iasi counties; OTA—Ialomita and Braila counties. Red shape: Baragan Plain (44–45° N, 26–27° E) in the Southern Plain and Dobrogea region.
Toxins 10 00533 g004
Table 1. Correlations between mycotoxins (deoxynivalenol (DON)) and the geographic position of the Baragan Plain, aridity indices, and agrometeorological factors by year in Romania in 2012–2015 (Pearson correlation coefficient).
Table 1. Correlations between mycotoxins (deoxynivalenol (DON)) and the geographic position of the Baragan Plain, aridity indices, and agrometeorological factors by year in Romania in 2012–2015 (Pearson correlation coefficient).
Correlations for the Baragan Plain
(44–45° N, 25–27° E)
Annual Agrometeorological Factors, 2012–2015May–June Agrometeorological Factors, 2012–2015Aridity Indices, 1900–2000
Cumulative PrecipitationAverage air TemperatureSoil moisture ReserveAverage air TemperatureCumulative PrecipitationDe Martonne (Iar-dM)Climatic Water Deficit (CWD)
DON All YearsPearson Correlation−0.296 *−0.028−0.0400.120−0.1950.003−0.029
Significance (two-tailed)0.0350.8440.7790.4040.1690.9830.837
N51515151515151
DON 2012Pearson Correlation0.521 *0.092−0.2770.3530.2790.0590.315
Significance (two-tailed)0.0220.7090.2500.1390.2470.8110.189
N19191919191919
DON 2015Pearson Correlation0.619−0.057−0.215−0.4350.801 * −0.573
Significance (two-tailed)0.1020.8930.6080.2810.017 0.138
N8888888
*—correlation is significant at the 0.05 level (two-tailed).
Table 2. Correlations between mycotoxins (aflatoxins (AF), ochratoxin A (OTA), zearalenone (ZEA)) and the geographic position of the Baragan Plain, aridity indices, and agrometeorological factors by year in Romania in 2012–2015 (Pearson correlation coefficient).
Table 2. Correlations between mycotoxins (aflatoxins (AF), ochratoxin A (OTA), zearalenone (ZEA)) and the geographic position of the Baragan Plain, aridity indices, and agrometeorological factors by year in Romania in 2012–2015 (Pearson correlation coefficient).
Correlations for the Baragan Plain
(44–45° N, 25–27° E)
Annual Agrometeorological Factors, 2012–2015July–August Agrometeorological Factors, 2012–2015Aridity Indices, 1900–2000
Cumulative PrecipitationAverage air TemperatureSoil moisture ReserveAverage air TemperatureCumulative PrecipitationDe Martonne
(Iar-dM)
Climatic Water Deficit (CWD)
AF All YearsPearson Correlation−0.0210.064−0.1420.025−0.125−0.070−0.062
Significance (two-tailed)0.8760.6420.2960.8560.3570.6080.652
N56565656565656
OTA All YearsPearson Correlation0.416 **−0.317 *−0.112−0.1560.107−0.0940.044
Significance (two-tailed)0.0060.0380.4810.3190.4950.5500.777
N43434243434343
OTA 2014Pearson Correlation0.904 *0.3290.2360.3570.796 0.239
Significance (two-tailed)0.0130.5250.6520.4870.058 0.648
N6666666
ZEA All YearsPearson Correlation0.217−0.118−0.046−0.0430.016−0.136−0.124
Significance (two-tailed)0.1420.4280.7570.7740.9150.3620.406
N47474747474747
*—correlation is significant at the 0.05 level (two-tailed); **—correlation is significant at the 0.01 level (two-tailed).

Share and Cite

MDPI and ACS Style

Gagiu, V.; Mateescu, E.; Armeanu, I.; Dobre, A.A.; Smeu, I.; Cucu, M.E.; Oprea, O.A.; Iorga, E.; Belc, N. Post-Harvest Contamination with Mycotoxins in the Context of the Geographic and Agroclimatic Conditions in Romania. Toxins 2018, 10, 533. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120533

AMA Style

Gagiu V, Mateescu E, Armeanu I, Dobre AA, Smeu I, Cucu ME, Oprea OA, Iorga E, Belc N. Post-Harvest Contamination with Mycotoxins in the Context of the Geographic and Agroclimatic Conditions in Romania. Toxins. 2018; 10(12):533. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120533

Chicago/Turabian Style

Gagiu, Valeria, Elena Mateescu, Ileana Armeanu, Alina Alexandra Dobre, Irina Smeu, Mirela Elena Cucu, Oana Alexandra Oprea, Enuta Iorga, and Nastasia Belc. 2018. "Post-Harvest Contamination with Mycotoxins in the Context of the Geographic and Agroclimatic Conditions in Romania" Toxins 10, no. 12: 533. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins10120533

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