Next Article in Journal
Brain Active Areas Associated with a Mental Arithmetic Task: An eLORETA Study
Next Article in Special Issue
A New Phenolic Acid Decarboxylase from the Brown-Rot Fungus Neolentinus lepideus Natively Decarboxylates Biosourced Sinapic Acid into Canolol, a Bioactive Phenolic Compound
Previous Article in Journal
Machine-Learning-Based Prediction of 1-Year Arrhythmia Recurrence after Ventricular Tachycardia Ablation in Patients with Structural Heart Disease
Previous Article in Special Issue
Mass Transfer Enhancement in High-Solids Anaerobic Digestion of Organic Fraction of Municipal Solid Wastes: A Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Mycotoxins in Silage on Biogas Production

1
Agricultural Research, Ltd., Zahradní 1, 664 41 Troubsko, Czech Republic
2
Department of Agricultural, Food and Environmental Engineering, Faculty of AgriSciences, Mendel University in Brno, Zemědělská 1, 613 00 Brno, Czech Republic
3
Department of Experimental Biology, Section of Microbiology, Faculty of Science, Masaryk University, Kamenice 753/5, 625 00 Brno, Czech Republic
4
Department of Agrochemistry, Soil Science, Microbiology and Plant Nutrition, Mendel University in Brno, Zemědělská 1, 613 00 Brno, Czech Republic
5
Agrovyzkum Rapotin, Ltd., Vyzkumniku 267, 788 13 Rapotin, Czech Republic
6
Department of Agrosystems and Bioclimatology, Faculty of AgriSciences, Mendel University in Brno, Zemědělská 1, 613 00 Brno, Czech Republic
*
Authors to whom correspondence should be addressed.
Bioengineering 2023, 10(12), 1387; https://doi.org/10.3390/bioengineering10121387
Submission received: 18 October 2023 / Revised: 22 November 2023 / Accepted: 26 November 2023 / Published: 2 December 2023
(This article belongs to the Special Issue From Residues to Bio-Based Products through Bioprocess Engineering)

Abstract

:
Mycotoxins can pose a threat to biogas production as they can contaminate the feedstock used in biogas production, such as agricultural crops and other organic materials. This research study evaluated the contents of deoxynivalenol (DON), zearalenone (ZEA), fumonisin (FUM), and aflatoxin (AFL) mycotoxins in maize silage prior to it being processed in a biogas plant and in digestate produced at the end of the anaerobic digestion (AD) process. In the experiment, three samples of silage were collected from one silage warehouse: Variant 1 = low contamination, Variant 2 = medium contamination, and Variant 3 = heavy contamination, which were subjected to investigation. A significantly reduced biogas production was recorded that was proportional to the increasing contamination with molds, which was primarily due to the AD of silage caused by technologically erroneous silage treatment. The AD was connected with changes in silage composition expressed by the values of VS content, sugar content, lactic acid content, acetic acid content, and the ratio of lactic acid content to acetic acid content. The production of biogas and methane decreased with the increasing contents of NDF, ADF, CF, and lignin. The only exception was Variant 2, in which the content of ADF, CF, and lignin was lower (by 8–11%) than that in Variant 1, and only the content of NDF was higher (by 9%) than that in Variant 1. A secondary factor that also correlated with changes in the composition of the substrate was the development of undesirable organisms, which further contributed to its degradation and to the production of mycotoxins. It was also demonstrated in this study that during the AD process, the tested mycotoxins were degraded, and their content was reduced by 27–100%. Only the variant with low mold contamination showed a DON concentration increase of 27.8%.

1. Introduction

Cereals (wheat, rice, maize, etc.) are the most important agricultural crops worldwide, the main reason being their relatively high nutritional value and ease of growth [1]. According to the FAO [2], the annual global growing area of maize and wheat in 2020 was ca. 202 and 219 million hectares, respectively, with maize being widely used as a food or fodder crop. The production of silage is a widely used practice to retain the nutritional value of farm animal feeds made of cereal biomass, most frequently by means of spontaneous lactic fermentation in anaerobic conditions [3].
Nevertheless, depending on growing conditions, cereals can be infested by multiple strains of toxicogenic fungi, with the majority of them producing more than one type of mycotoxin [1]. Mycotoxins are secondary metabolites of fungi with diverse chemical structures that contaminate many of the most frequently consumed food and feed products worldwide [4,5]. There are two ways for cereals to be contaminated by mycotoxins: (i) fungi growing as pathogens on plants or (ii) saprophytically on stored crops [6]. In the past, mycotoxins were recorded particularly in wheat and maize grains, which exhibited the highest concentrations of fumonisin (FB), deoxynivalenol (DON), aflatoxin (AF), and zearalenone (ZEN) [1].
The development of mycotoxins depends on a range of conditions during the pre-harvest and post-harvest handling of agricultural products that can be divided according to Adegoke and Letum [7] and Magan et al. [8]: (A) internal factors: moisture content, water activity, substrate type, plant type, and composition of nutrients; (B) external factors: climate, temperature, oxygen level; (C) factors of processing: drying, mixing, addition of preservatives, grain handling; and (D) implicit factors: interaction with insects, strains of fungi, microbiological ecosystem.
Adegoke and Letum [7] also include higher temperatures, abundant rains during harvest (biomass moisture content), stand density, mechanical harvesting, vegetation stage, drought, variety, soil structure, and temperatures below freezing as factors participating in the spread of molds and the development of mycotoxins. Warmer climatic conditions promote a more frequent occurrence of aflatoxins and fumonisins. In contrast, colder regions with high humidity support the development of ochratoxins, zearalenone, deoxynivalenol (DON), T-2 toxin, and diacetoxyscirpenol (DAS) [9]. The three most frequently occurring toxicogenic fungi in Europe are Aspergillus, Penicillium, and Fusarium [10].
In farm animals, the consumption of feeds contaminated by mycotoxins may have unfavorable consequences such as reduced feed intake, refusal of feed, poor feed conversion, reduced increment in body weight, increased occurrence of diseases (due to suppressed immunity), and impaired capacity of reproduction [11,12], which lead to economic losses [13,14]. Many batches of cereals cannot be used for food, feed, or energy because of their high contamination by mycotoxins. In recent years, increasing attention has been given to renewable sources of energy that can replace fossil fuels. The process of anaerobic digestion (AD), whose goal is the production of biogas from diverse organic waste materials and suitable energy crops, appears to be promising [15,16]. In addition to biogas, a side product of the AD process is digestate, which can have high agronomic value as fertilizer thanks to the high content of mineral nitrogen and relatively stable organic matter [17,18]. Digestate was demonstrated to have positive effects in agriculture on the growth and production of crops if applied to post-harvest residues or in combination with readily degradable organic matter [19] as compared with other kinds of organic fertilizers [20]. Processes taking place during the AD of the organic feedstock enhance the properties of the biomass source; e.g., it improves the availability of N in digestate and P in plants [21]. Thanks to these properties, fertilization based on digestate can be as effective as conventional organic supplements such as farmyard manure or compost [22].
Thus, AD may represent an interesting approach to biogas production from organic waste materials (such as silage) contaminated by mycotoxins. However, when using material contaminated by mycotoxins as a feedstock for biogas production, one has to consider the possible negative influences on the process of AD, biogas yield, and digestate quality [1,22]. These aspects were dealt with in several studies that revealed reduced contents of several mycotoxins in mesophilic and thermophilic conditions, both in batch tests and in semi-continual reactors [23,24,25]. The main problem with using silages highly contaminated with mycotoxins is their influence on microbial diversity in the fermenters of biogas plants. The process of AD itself can reduce the content of mycotoxins by more than 50% [1] because some mycotoxins are susceptible to biological degradation [26]. However, the problem during the AD of silage contaminated by mycotoxins is the decreased production of biogas [27]. The reason is simple: increased concentrations of mycotoxins impair the diversity of the microbial community [1], increase the accumulation of organic acids, and decrease the pH value in the fermenter [28]. All these effects then result in AD inhibition and, hence, the reduced production of biogas [1,27,28]. In the process of AD, individual types of mycotoxins exhibit different potentials for biodegradation [27,28]. The mycotoxins with the greatest potential for biodegradation are AF B1, FM, and DON. It was found that the process of biodegradation through AD can reduce the concentration of these mycotoxins by more than 50% [1,27,28]. It is also necessary to realize that silage AD gives rise not only to biogas but also to fermentation residues in the form of digestate, which can be contaminated by mycotoxins too [28]. Additionally, when digestate is used as fertilizer, it is necessary to investigate its influence on the quality of food and feed products not only in terms of nutritional value [29], but also in terms of soil properties [30,31], as well as the possible contamination of the environment with mycotoxins [32].
This study presents a novel combined evaluation of the content of mycotoxins in the contaminated silage, their effect on biogas and methane production, and the degradation of mycotoxins during the AD process (monitored as the residual content of mycotoxins in the digestate). This study has the following objectives: (1) to determine the degree to which the presence of mycotoxins in contaminated silage affects the production of biogas or methane and (2) to determine whether degradation of tested mycotoxins takes place in the process of silage AD. Furthermore, the hypotheses of this study are as follows: H0: Silage contamination by mycotoxins does not affect biogas and methane production during anaerobic digestion. H1: Degradation of tested mycotoxins does not occur in the process of silage AD.

2. Materials and Methods

2.1. Organization of the Experiment—Collection of Contaminated Samples

The experiment was based on testing contaminated silage under laboratory conditions with regard to biogas production and analyzing some qualitative indicators. A silage pit with visible mold infestation was selected for the experiment.
Three variants of maize silage were sampled from the silage pit (Figure 1): Variant 1 = low contamination (LC)—visibly non-infected sample; Variant 2 = medium contamination (MC); Variant 3 = heavy contamination (HC). The variants were of different quality in terms of visible mold infestation (Figure 1), with three replicates each, i.e., a total of nine samples.

2.2. Silage Characteristics

The selected qualitative indicators (N-substances, crude fiber, acid detergent fibers, etc.) were analyzed in the collected samples of contaminated silage to determine the effects of mycotoxin contamination on silage quality. N-substances were determined by the Kjeldahl method according to ISO 20483:2013 [33] using the KjeltecTM 2300 analyzer (FOSS Analytical, Hillerød, Denmark), and the protein content was calculated from the N-substance value by multiplying the N-substance value by the empirical factor 6.25. The fat content was determined gravimetrically using the water-cooled Soxhlet [34] extractor BEHR 6 (Behr Labor-Technik GmbH, Düsseldorf, Germany) by direct sample extraction with diethyl ether. The crude fiber (CF) content was determined using two-stage hydrolysis with sulfuric acid and potassium hydroxide. The ash content was determined according to ISO 6865:2000 [35]. The acid detergent fiber (ADF) content was determined using a solution of concentrated sulfuric acid and cetyltrimethylammonium bromide (CTAB). The acid detergent lignin (ADL) was determined according to ISO 13906:2008 [36]. Neutral detergent fiber (NDF) was determined using a solution of sodium lauryl sulfate and ethylenediaminetetraacetic acid according to ISO 16472:2006 [37]. The following parameters were determined during silage analysis: dry matter, VS, ash, pH, aqueous solution, lactic acid, acetic acid, propionic acid, and butyric acid. The content of dry matter (=total solids, TS/DM) and volatile solids in the samples was determined gravimetrically by drying in the electric furnace LMH 07/12 (LAC, Židlochovice, Czech Republic) at 105 °C to constant weight and by annealing the dried samples at 550 °C to constant weight in accordance with the standards ČSN EN 15934 [38] and ČSN EN 15935 [39]. All the parameters mentioned, including the description of their determination, have already been given in the work of Hunady et al. [40] and Kintl et al. [41].
The starch content of the silage was determined using the polarimetric method (Polamat S, Carl Zeiss Jena GmbH, Jena, Germany) in accordance with ISO 6493:2000 [42]. The fat content was determined gravimetrically using the water-cooled Soxhlet extractor by direct sample extraction with petroleum ether [43]. The content of sugars (reducing saccharides) was determined using the Luf–Schoorl method [44,45,46]. The pH of the H2O leachate was determined according to ČSN 46 7092-42 [47]. Lactic acid, acetic acid, propionic acid, and butyric acid were determined according to Stringer [48]—using the IONOSEP 2001 capillary isotachophoresis instrument (RECMAN—laboratory technique, Ltd., Ostrava, Czech Republic).

2.3. Digestion Batch Tests

The digestion batch tests were performed according to Kintl et al. [41], using systems for batch tests (Figure 2). The inoculum was digestate from the biogas plant in Čejč, Czech Republic, processing maize silage and slurry, with the following parameters: mesophilic temperature conditions, 38 °C; dry matter, 3.43 ± 0.06%; volatile solids, 68.54 ± 0.11%; pH, 7.0; FOS, 740 mg/L; TAC, 2200 mg/L; NH4+, 590 mg/L. The biogas composition was analyzed using a Dräeger X-am 5600 biogas analyzer (Drägerwerk AG & Co. KGaA, Lübeck, Germany), as described in the work of Kintl et al. [49] and Kintl et al. [50]. The volume of biogas generated was converted to standard temperature and pressure (273.15 K and 1 bar). The process parameters were as follows: initial organic loading rate = 4.5 gvs introduced substrate/L; retention time = 21 days; temperature during the test = 42 °C ± 0.1 °C. The silage parameters can be found in Table 1 in the Results section.

2.4. Analysis of Mycotoxins

Samples of silage and digestate were dried at 60 °C and then analyzed for the content of mycotoxins deoxynivalenol (DON), aflatoxin (AFL), zearalenone (ZEA), and fumonisin (FUM) using the enzymatic immunosorbent assay (ELISA) according to Skládanka et al. [51]. The ELISA (MyBioSource, San Diego, CA, USA) is a competitive direct enzyme test for the qualitative analysis of plant biomass focused on the content of mycotoxins. The concentrations of individual mycotoxins were expressed in micrograms per kg of silage or digestate [52].

2.5. Statistical Analysis—Data Treatment

All experimental parameters were measured in three replicates. The data were processed in the Statistica 14 program (TIBCO Software, Inc., Palo Alto, SF, USA). The procedures used to analyze the data included exploratory data analysis (EDA), one-way analysis of variance (ANOVA), and Tukey’s HSD post hoc test. The data were standardized, and the relationships between the values of the measured parameters were analyzed using correlation analysis, factor analysis (FA), principal component analysis (PCA), and cluster analysis (CA). All analyses were performed at a significance level of p < 0.05.

3. Results and Discussion

3.1. Qualitative Parameters of Silage

The quality of the individual silage variants was assessed by determining the following parameters: volatile solids content (VS), dry matter (DM), neutral detergent fiber (NDF), acid detergent fiber (ADF), crude fiber (CF), lignin, proteins, starch, sugars, lactic acid (LA), acetic acid (AA), titratable acidity (TA), and pH. The measured parameters are divided into three subsections that characterize the content or the state of selected indicators: Section 3.1.1, DM, pH, VS, proteins, lipids, sugar; Section 3.1.2, Starch, NDF, ADF, CF, lignin; Section 3.1.3, LA, AA, LA/AA, TA, propionic acid, butyric acid. In general, measurement data show that contamination of the silages with molds leads to a deterioration in their qualitative indicators. The development of the measured indicators confirms the importance of adhering to the correct ensiling process [53]: (a) appropriate harvest date (dry matter up to 35%); (b) quality of the chopped material—a length of 15–25 mm is sufficient at the optimal stage of maize maturity; (c) ensiling process—oxygen must be displaced from the harvested matter (shreddings), and then further access of oxygen into the silage must be prevented. An anaerobic environment has to be created in which no atmospheric oxygen is supplied. Only then can the respiration of cells stop and the fermentation of the lactic acid bacteria begin. A necessary factor that influences the quality of the resulting silage is also the supplementation of silage additives, most commonly lactic acid bacteria.

3.1.1. DM, pH, VS, Proteins, Lipids, and Sugar

The DM content in the silage ranged from 34.3% (MC variant) to 47.9% (HC variant). Differences between the variants were statistically significant (Table 1). According to Wilkinson [54], the ideal DM content in silage is 30–35%, with excessive DM content generally associated with susceptibility to molds. This indicates that all variants in our experiment with the exception of MC had a higher dry matter content than the value described as optimal by Wilkinson [54], creating the conditions for the development of molds (particularly in HC).
The pH value in the aqueous solution (Table 1) increased significantly with increasing optical contamination from 4.22 (LC) and reached 6.10 in the HC variant. Kung et al. [55] state that a normal pH value for silage with a DM content of 30–40% is 3.7–4.0. The pH values of our test variants were higher, and significantly higher again in the HC variant. As the pH value is a physical factor that ensures the preservation of the silage, pH values above 6 indicate the probable development of molds and are evidence of poor silage management [56].
The main cause of the drop in silage pH is the production of lactic acid (LA), a typical fermentation product during the ensiling of maize, the amount of which is significantly higher than the amount of other organic acids (acetic, propionic, and butyric acids) normally detected during ensiling [57]. This is also confirmed by the results of our study, which showed a very strong negative correlation (R = −0.94) between the pH value and the concentration of LA in the silage (Table A1).
Silage-degrading microorganisms are more effectively inhibited at lower pH values than at higher ones. In the past, silages with a pH value of more than 4.2 were considered of poor quality [58]. Currently, silages with a high DM content are produced on a large scale and are stable even at relatively high pH values. Although pH remains one of the most important indicators of good fermentation in low-DM silages, it is not a reliable indicator in high-DM silage [59]. This research also shows that DM content is very strongly correlated with pH (R = 0.95) and that the values of both parameters increase with increasing contamination of silage.
The VS content ranged from 92.05%DM (MC) to 94.38%DM (LC) and was significantly highest in the low contamination (LC) variant (Table 1). According to Heuzé et al. [60], the average VS content in maize silages with DM content of 35–40% ranges from 94.8 to 97.8%DM. Wilkinson [54] claims that an optimum VS content should be above 92%DM.
The content of proteins, lipids, and sugars showed significant differences between the variants. The highest protein content was found in the HC variant, and the highest lipid and sugar contents were found in the LC variant (Table 1). Heuzé et al. [60] state that the average content of proteins and lipids in maize silage with 35–40% dry matter is approx. 6.8%DM and 3.0%DM, respectively. Therefore, we assume that the organic substances were degraded during the ensiling process and subsequently contaminated by molds. As microorganisms, molds first utilize readily available organic substances (sugars) and only then utilize the more complex forms such as lipids [61]. Although this explains the reduced content of lipids and sugars, it does not explain the slightly increased content of proteins, which could have been caused by protein-based substances produced by molds [62]. Kung et al. [55] found that clostridial silages, for example, are often characterized by higher than normal pH and higher than normal soluble protein concentrations. The mycelial biomass of various fungi can be used as a rich source of mycoproteins [63], and the production of microbial proteins by filamentous fungi is therefore a prospective bioprocess in the food and feed industry [64]. In some fungal groups, the mycelium is able to convert the otherwise difficult-to-process lignocellulosic substrates into proteins and can therefore be used as a sustainable source [65].

3.1.2. Starch, NDF, ADF, CF, and Lignin

The experimental variants did not differ statistically in their starch content. However, the fiber (NDF, ADF, CF, Table 2) and lignin contents were significantly highest in the HC variant. Heuzé et al. [60] claim that the average starch content in maize silage with 35–40% dry matter is approx. 31.6%DM. Based on a meta-analysis, García-Chávez et al. [66] give an average value for the starch content of 23.3%DM.

3.1.3. LA, AA, LA/AA, TA, Propionic Acid, Butyric acid

The LA content ranged from 0.65%DM (HC) to 1.44%DM (LC), and the difference between LC and MC on the one hand and HC on the other was statistically significant (Table 3). Kung and Shaver [57] claim that the usual values of LA concentration in maize silages with 30–40% dry matter are between 4 and 7%DM.
Similarly, the AA content decreased with increasing contamination, although the differences between the variants were statistically significant. According to Wilkinson [54], the optimum AA content in maize silage is 2–3%DM. In our experiment, the concentrations of both LA and AA were very low in all variants (Table 3) and were below the values described as optimal.
According to Kung et al. [55], lactic acid is one of the most important substances indicating the course of fermentation in silages with a low dry matter content. During the ensiling process, the LA produced by the lactic acid bacteria usually occurs in the highest concentration and is most strongly involved in lowering the pH value during fermentation, as it is about 10 to 12 times stronger than all other major acids [55]. However, in silages with a high DM content, the LA content does not always indicate successful fermentation [59]. Fermentation during which LA is produced leads to the lowest losses of dry matter and energy during storage [57]. The LA/AA ratio is commonly used as an indicator of the quality of fermentation. During optimal silage fermentation, the ratio of these acids usually ranges from approximately 2.5 to 3.0 [55]. An LA/AA ratio below 1 is usually an indication of abnormal fermentation. Low concentrations of AA may not be sufficient to inhibit lactate-assimilating yeast [55,57].
In our experiment, the LA/AA ratio ranged from 0.94 (MC variant) to 1.99 (HC variant) and was statistically significantly highest in the HC variant compared to the other two variants (Table 3). This shows that the values of all three variants were below the interval considered to be optimal. One reason for this was the very low content of LA, which was far below the usual values in maize silage [57]. The values of all parameters thus clearly indicate an increase in aerobic degradation, which correlates with the degree of silage contamination by molds. The aerobic phase can change the chemical composition of the silage after the silo has been opened. Increased pH (>6.0) caused by the growth of yeasts responsible for silage degradation leads to active growth of toxicogenic fungi during the feeding phase, especially in poorly treated silages [56]. The presence of propionic acid and butyric acid was below the detection limit of 10 mol/L.

3.2. Content of Mycotoxins in Silage

The levels of mycotoxins (DON, ZEA, FUM, AFL) found in maize silage in this study were on average lower than the values reported in the available literature, and the degree of contamination was not able to be distinguished [67,68,69,70].

3.2.1. DON

The average DON concentration in all samples was 361.7 µg/kgDM and ranged from 180 to 670 µg/kgDM (Table 4, Figure 3). The DON content increased with increasing contamination of the samples and was statistically significantly higher in the HC variant (638.3 µg/kgDM) than in the LC and MC variants (195.7 and 251.0 µg/kgDM).
Based on the previous studies, it can be stated that DON is one of the most frequently detected mycotoxins in silage, and its concentrations can be very high [71,72]. The average values given in the literature for the DON content in maize range from 280 [72] to 3142 µg/kgDM [69]. A worldwide study of samples of feed for livestock, which lasted 3 years, showed that DON is a common threat to livestock with an occurrence of 59% and an average contamination of 1104 μg/kg [73]. In the Netherlands, a new study was carried out on the presence of mycotoxins in silages. It was found that the major sources of DON in dairy cattle feed are maize and wheat silages with average concentrations of 854 and 621 μg/kg and maximum concentrations of 3142 and 1165 μg/kg, respectively [69]. Another recent multi-year study of mycotoxins in Poland [70] revealed the presence of DON in 86% of 143 samples of maize silage with average and maximum concentrations of 223 and 7860 μg/kg, respectively.
The DON content we found in the most contaminated HC variant (638.3 µg/kgDM) was therefore rather below average compared to the average values mentioned in the above studies.

3.2.2. ZEA

The average content of ZEA was 115.3 µg/kgDM and ranged from 63 to 147 µg/kgDM (Table 4). Compared to LC and MC, the statistically significantly lowest ZEA content was found in the HC variant (Figure 3).
The values mentioned in the literature range on average from 66 [68] to 432 µg/kgDM [74]. The production of ZEA is supported by conditions with high humidity and alternating low (11–14 °C) and mild (27 °C) temperatures [75]. Whitlow and Hagler [76] reported an average ZEA concentration of 525 μg/kgDM and an occurrence of 30% in 461 samples of maize silage in the USA. In the Netherlands, Driehuis et al. [69] detected ZEA in 50% of 140 maize silage samples, with an average concentration of 146 μg/kgDM. According to Storm et al. [68], zearalenone was the most frequently detected mycotoxin in a study evaluating the exposure of Danish cattle to mycotoxins from maize silage. Furthermore, Rodrigues and Naehrer [73] reported that approximately 45% of 7049 samples of feed for livestock collected in America, Europe, and Asia contained ZEA with an average concentration of 233 μg/kgDM.

3.2.3. FUM

The FUM content also decreased with increasing contamination of the samples, but the decrease was not statistically significant (Table 5, Figure 3). The average FUM content in all silage samples was 239.4 µg/kgDM. The lowest and highest values (194 and 296 µg/kgDM) were found in the MC variant.
The lowest and highest values (194 and 296 µg/kgDM) were found in the MC variant. The lowest and highest average values for maize silage published in the literature are 11.3 µg/kgDM [70] and 3000 µg/kgDM [77], respectively. In the Netherlands, Driehuis et al. [69] report a maximum value of 34,000 µg/kgDM in maize silage. The most important factors favoring the secretion of fumonisins by Fusarium species are hot and dry periods, followed by conditions with increased humidity and damage by insects [69,77]. According to Rodrigues and Naehrer [73], fumonisins are the most common threat among mycotoxins for livestock, with an occurrence of 64% and an average concentration of 1965 μg/kgDM in 7049 samples of livestock feed collected in America, Europe, and Asia.
According to Schmidt et al. [78], approximately 48.6% of 327 maize silages sampled in Brazil were contaminated with fumonisin B1; the average concentration was 369 μg/kg. Gonzalez-Pereyra et al. [56] reported that in their study, the content of fumonisin in maize silage ranged from 340 to 2490 μg/kgDM, with higher concentrations found in samples from the upper layer and sidewalls of the silage, which are usually more susceptible to degradation due to aerobic conditions and less intensive fermentation.

3.2.4. AFL

The highest average aflatoxin content was found in the HC variant (3.7 µg/kgDM). However, the differences between the variants were not statistically significant (Table 5, Figure 4). The maximum value was 4.0 µg/kg. Limit values are specified, for example, in the European Commission Regulation No. 650/2010 with a maximum tolerated value of 5 µg AFLB1/kg.
In the literature, the average values of AFL content range from 3 µg/kgDM [78] to 33 µg/kgDM [79]. In particular, in the intercontinental research study by Rodrigues and Naehrer [73], aflatoxins were present in 33% of 7049 samples of feed for livestock at an average concentration of 63 μg/kg. And in the survey of Argentinian dairy farms that suspected problems with mycotoxins, about 30.2% of maize silage samples contained more than 20 μg/kg of aflatoxins, with a maximum concentration of 42.1 μg/kg.
AFL is one of the mycotoxins produced by Aspergillus spp. which are considered storage molds because they do not normally infect the crop before harvest. However, some species, such as Aspergillus flavus, can also infect plants in the field and produce aflatoxins during periods of high temperatures (>32 °C), high humidity (>80%), or drought stress [71,80].
The relationships among the individual mycotoxins were studied too. Table 6 shows a very strong negative dependence between DON and ZEA (R = −0.96) and ZEA and AFL (R = −0.78). This dependence indicates that higher contents of DON and AFL in silage can result in a lower content of ZEA, and to some extent also of FUM, which also has a negative (although weaker) correlation (R = −0.40) with the content of DON.
The correlation analysis further shows (Table A1) that the described relationships between the studied mycotoxins were likely to be caused by different requirements of molds for the environment [81]. In our study, a very strong positive correlation was observed between DON and pH; dry matter; and contents of NDF, ADF, CF, and lignin. On the other hand, a strong negative correlation was recorded between DON and the contents of lipids, sugar, LA, and AA. This character of correlations was to a lesser extent (moderately strong dependence) observed for AFL. By contrast, for ZEA and to a lesser extent for FUM, the correlations of the mentioned parameters were opposite.
Interactions between the mycotoxin content in the maize silage and the qualitative parameters were analyzed in detail using PCA (Appendix C). Three basic factors were identified (Table A4, Figure A2 and Figure A3) with Factor 1 and Factor 2 explaining more than 90% of the variability in the measured values (Figure A4). Factor 1 was labeled as the main factor as it explained 78.77% of the variability in the measured values. The factor correlated positively with the selected indicators of contaminated silage, but negatively with the production of biogas and its quality. It can therefore be assumed that the factor describes the level of maize silage contamination. Factor 2 also showed a negative correlation with the production of biogas and its quality, which was stronger than that of Factor 1. Considering that Factor 2 explained only 21.23% of the variability and had a demonstrably negative influence on biogas production, we assume that this factor represents a kind of mycotoxin.

3.3. Effect of Mycotoxins on the Production of Biogas and Methane

The production of biogas and methane ranged from 0.4305 (MC variant) to 0.6178 m3/kgVS (LC variant) and from 0.2248 (MC variant) to 0.3792 m3/kgVS (LC variant), respectively. The measured values are summarized in Table 7. The LC variant exhibited a statistically significantly higher biogas and methane yield than the MC and HC variants.
In particular, the methane concentration in the biogas showed no statistically significant differences between the respective variants. The final concentration of methane in the biogas ranged from 63%vol to 64%vol (Table 8).
Correlation, cluster (Table A2), factor, and principal component analyses were carried out to classify the measured parameters and mycotoxin contents according to their similarity and the degree of their influence on the production of biogas and methane. The results of these analyses and their comparison are shown in Table 9 and Table A1, Table A2, Table A3, Table A4 and Table A5, as well as in Figure 5 and Figure A1, Figure A2, Figure A3 and Figure A4.
The resulting categorization into five clusters (Table 9) is consistent with the results of the factor analysis with three factors of which Factor 1 represents the members of three clusters: Clusters 2, 3 and 4. Cluster 2 includes variables with the negative factor loadings of Factor 1, while Clusters 3 and 4 have positive factor loadings (except proteins) (Table 9). The results of the correlation analysis (Table 9 and Table A1) show that none of the four mycotoxins has a significant correlation with the production of biogas or methane. As mentioned in the above subsection, Factor 1 explained the influence of the mycotoxin concentration in the silage, and Factor 2 explained the type of mycotoxin.
Although a moderately strong negative, statistically non-significant correlation was found between biogas and methane production and DON content (R = −0.54 and R = −0.48, respectively), this dependence can be considered secondary, as much stronger significant positive correlations were found between biogas and methane production and parameters belonging to Cluster 1 (VS) and Cluster 2 (sugars, lipids, and LA). In contrast, the DON content had a strong negative correlation with the members of Cluster 2 (Table A1). This result is consistent with the conclusions of Ferrara et al. [1], according to which the high content of mycotoxins had no influence on the overall production of biogas, especially methane. The same conclusions were reached by De Gelder et al. [24].
According to Merrettig-Bruns and Sayder [22], the fermentation of moldy maize silage as the only substrate led to a lower biogas yield compared to the fermentation of maize silage without molds. The authors also report that an extended analysis of subgroups revealed a significantly lower content of easily digestible carbohydrates in silages made from moldy maize.
The comparison of the LC and HC variants shows that the biogas and methane yield was primarily reduced by the statistically significant decrease in the contents of VS, sugars, lipids, LA, and AA (members of Cluster 1 and Cluster 2) and at the same time by an associated statistically significant increase in the contents of NDF, ADF, and CF (members of Cluster 3) in the HC variant. However. The content of ADF, CF, and lignin was lower in the MC variant (by 8–11%) than in the LC variant, and the content of NDF was higher (by 9%) than that in the LC variant. The comparison of the LC and MC samples also showed only a slight difference in the content of lipids, LA, and AC and in the value of LA/AA. This process was accompanied by increased values for DM, pH, and LA/AA and a reduced value for TA.
Obviously, there were technological errors in the silage treatment, which led to changes in the substrate composition, resulting in reduced production of LA and thus insufficient preservation of the silage and its devaluation. This is also indicated by the titratable acidity (TA) values, which decreased in statistically significant proportion to the degree of contamination (Table 3).
The primary factor that reduced the production of biogas and methane is therefore the aerobic degradation of the silage in correlation with the changes in its composition expressed by contents of VS, sugars, LA, and AA and their ratios. The ratios showed a moderately strong to strong positive correlation, increasing the production of both gases. In contrast, the biogas and methane yield decreased with increasing content of NDF, ADF, CF, and lignin. The secondary factor related to the changes in substrate composition was the development of undesirable microorganisms that further contributed to substrate degradation and mycotoxin production. According to McDonald et al. [82], yeasts can oxidize lactic acid and thus increase the pH of the silage, which favors the growth of other microorganisms.
The division of mycotoxins (DON, ZEA, FUM, AFL) into four clusters (Clusters 2, 3, 4, 5) could be related to the conditions under which microorganisms producing the respective toxins can develop successfully (Figure 5, Table A1).
  • Cluster 2—ZEA
ZEA has a strong positive correlation with the contents of lipids, sugar, LA, and AA and a negative correlation with the values of pH; DM; and the content of NDF, ADF, CF, and lignin.
2.
Cluster 3—DON
DON has a strong positive correlation with the values of pH; DM; and the content of NDF, ADF, CF, and lignin and a negative correlation with the contents of lipids, sugar, LA, and AA.
3.
Cluster 4—AFL
With respect to the mutual proximity of Cluster 3 (DON) and Cluster 4 (AFL) (Figure 5), the AFL content has a moderately strong positive correlation with the values of DM; pH; and lignin, NDF, ADF, and CF content and, similar to the DON content, a negative correlation with members of Cluster 2 (ZEA) (LA, AA, lipids).
4.
Cluster 5—FUM
This cluster contains only FUM and starch. The contents of both FUM and starch show only a weak or slight dependence on the variables in the other clusters and also on each other. These two variables are the only ones represented by Factor 3 (Table 9). It can therefore be assumed that the levels of FUM and starch are related to parameters that were not considered in this study. This fact is illustrated in Figure A3.
Silages may contain a mixture of mycotoxins resulting from both pre-harvest contamination with Fusarium spp. and Aspergillus spp. [72,83] and post-harvest contamination by molds which commonly occur in silages, such as Penicillium spp. or Aspergillus spp. [84,85]. According to Zain [86], there are many common molds that do not produce mycotoxins. The presence of molds in silage does not have to always indicate the presence of mycotoxins, and the absence of molds does not indicate the absence of mycotoxins. Conditions for the growth of molds and for the development of mycotoxins do not have to be always the same [76]. Molds of Fusarium spp. for example can grow intensively at temperatures ranging from 25 to 30 °C without producing mycotoxins, while at temperatures below the freezing point, they produce large quantities of mycotoxins with minimal growth [76,87].

3.4. Content of Mycotoxins in Digestate

The comparison of the contents of mycotoxins DON, ZEA, FUM, and AFL in the silage and in the digestate (Table 10) shows that the mycotoxins were degraded during the anaerobic fermentation process, and their content was reduced by 27–100%. The only exception is the LC variant in which the concentration of DON increased by 27.8%. This can be explained by the fact that the sample contained the lowest initial concentration of DON which was not metabolized by the microorganisms. During the process of anaerobic digestion, the mass of the initial feedstock was reduced, which could have led to the higher concentration of DON in the digestate.
Similar results were also published by Ferrara et al. [1], who reported in their study that the degradation of mycotoxins was approx. 54% for aflatoxin B1 and 60% for fumonisins. De Gelder et al. [24] report that in batch tests with mycotoxins, aflatoxin B1, ochratoxin A, deoxynivalenol, zearalenone, and T-2 toxin were degraded by more than 90%. According to Richter et al. [88], the ensiling of DON-contaminated maize groats showed a significantly reduced DON content; the DON content was also reduced during the fermentation of whole-grain maize. These facts show that the use of mycotoxin-contaminated substrates for biogas production is a good way for their processing, although some studies conducted with mycotoxins in contaminated raw materials for biogas production indicated that the toxins are degraded at different rates during the anaerobic fermentation process [23,24,28]. However, further investigation of the effect of mycotoxins on the methanogenic archaea in the fermenter and the transfer of mycotoxin residues from the digestate to arable land and thus to crops is a necessary prerequisite.
The measured levels of mycotoxins in silage and the influence of these mycotoxins on the process of biogas production and possibly on the quality of digestate confirm that appropriate attention should be paid to the ensiling process. According to Rada and Vlková [89], the ensiling process requires quality control and strict organization. Compaction of the ensiled material (either maize biomass or biomass of forage plants) is elementary, and the air must be pressed out. The addition of silage inoculants (lactic acid bacteria) is also necessary. As far as the compaction of the silage and its covering are concerned, there are not too many alternative technologies that would differ significantly. With silage inoculants, however, you have the choice between a large number of homofermentative and heterofermentative lactic acid bacteria [90,91,92]. Homofermentative bacteria produce lactic acid, while heterofermentative bacteria are also able to produce acetic acid [90]. Heterofermentative inoculants in particular open the way for innovations in the ensiling process, as these bacteria can also produce enzymes, in addition to lactic and acetic acids, which contribute to better stability of the ensiling process and to the quality of the silage. According to Rada and Vlková [89], however, it is necessary to distinguish between the type of plant biomass to be ensiled (different inoculants for maize biomass and for crops difficult to ensile such as alfalfa) and the type of preparation (length of the chopped material). Therefore, the choice of an ensiling agent is an important part of the ensiling process, influencing not only the quality of the resulting silage, but also its durability.

4. Conclusions

The presence of mycotoxins in animal feed has an unfavorable effect on the performance and health of farm animals and can also endanger human health. This study has shown that the use of anaerobic digestion (AD) to produce biogas and methane is a promising method for utilizing mycotoxin-contaminated maize silage that cannot be used as feed for livestock. None of the four mycotoxins tested (deoxynivalenol, zearalenone, fumonisin, aflatoxin) in the silage had a significant effect on the production of biogas (methane), either alone or in combination with the others. This means that the alternative hypothesis H1 was rejected.
The observed decrease in biogas yield, which was proportional to the increasing contamination by molds, was primarily caused by the aerobic degradation of the silage in correlation with the changes in its nutrient composition, as reflected in the values of VS, sugar, lactic acid, and acetic acid contents and their ratios which exhibited a statistically significant moderately strong to strong positive correlation, thus increasing the production of the two gases, while an increase in the contents of NDF, ADF, CF, and lignin reduced the production of biogas and methane with increasing contamination. The secondary factor that correlated with the changes in substrate composition was the development of undesirable microorganisms that further contributed to substrate degradation and the production of mycotoxins. In addition to biogas, a by-product of the AD process is digestate, which can be of high agronomic value due to its high mineral nitrogen content and relatively stable organic matter. This study showed that the mycotoxins tested were degraded during the AD process and that their levels were reduced by 27–100%. Only in one case, in the variant with low mold contamination, the concentration of deoxynivalenol increased by 27.8%.

Author Contributions

Conceptualization, A.K., T.V., I.H. and J.E.; methodology, A.K., T.V., I.H., J.S. and J.E.; validation, A.K., T.V., I.H., J.S., T.H., M.B., M.V. and J.E.; formal analysis, A.K., T.V., I.H., J.S., J.H., M.V. and J.E.; investigation, A.K., T.V., I.H., J.S., M.V. and J.E.; resources, A.K.; data curation, A.K. and I.H.; writing—original draft preparation, A.K., I.H., T.V. and J.E.; writing—review and editing, T.V. and J.E.; visualization, A.K., T.V., I.H. and J.E.; supervision, A.K., T.V. and J.E.; project administration, A.K.; funding acquisition, A.K. All authors have read and agreed to the published version of the manuscript.

Funding

The research was supported by the Ministry of Agriculture of the Czech Republic within the frame of the Long-Term Conception of the Development of Research Organization Agricultural Research, Ltd. Troubsko.

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 on request from the corresponding author (J.E.).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Correlation Analysis

Table A1. Correlation matrix. Red-marked correlation coefficients are significant at p < 0.05.
Table A1. Correlation matrix. Red-marked correlation coefficients are significant at p < 0.05.
VariableBIOGASMETHANEDONZEAFUMAFLpHDMVSNDFADFCFLIGNINSTARCHPROTEINSLIPIDSSUGARLAAA
BIOGAS1.000.99−0.540.390.270.00−0.55−0.280.90−0.66−0.30−0.31−0.150.360.410.630.900.680.55
METHANE0.991.00−0.480.350.210.00−0.50−0.220.93−0.61−0.23−0.25−0.130.370.470.570.870.650.50
DON−0.54−0.481.00−0.96−0.400.610.990.96−0.140.990.960.960.77−0.140.54−0.96−0.82−0.93−0.99
ZEA0.390.35−0.961.000.28−0.78−0.98−0.980.01−0.94−0.95−0.94−0.800.05−0.630.900.740.900.98
FUM0.270.21−0.400.281.000.03−0.38−0.360.08−0.39−0.33−0.32−0.470.29−0.150.370.380.170.38
AFL0.000.000.61−0.780.031.000.660.720.210.570.630.630.600.350.59−0.58−0.36−0.55−0.66
pH−0.55−0.500.99−0.98−0.380.661.000.95−0.170.990.940.940.79−0.150.51−0.95−0.84−0.94−1.00
DM−0.28−0.220.96−0.98−0.360.720.951.000.130.910.980.980.83−0.040.74−0.89−0.65−0.84−0.96
VS0.900.93−0.140.010.080.21−0.170.131.00−0.290.130.120.220.280.760.270.670.340.17
NDF−0.66−0.610.99−0.94−0.390.570.990.91−0.291.000.900.910.71−0.170.40−0.97−0.90−0.96−0.99
ADF−0.30−0.230.96−0.95−0.330.630.940.980.130.901.001.000.80−0.050.75−0.90−0.64−0.85−0.94
CF−0.31−0.250.96−0.94−0.320.630.940.980.120.911.001.000.83−0.080.74−0.89−0.64−0.87−0.94
LIGNIN−0.15−0.130.77−0.80−0.470.600.790.830.220.710.800.831.00−0.270.70−0.63−0.44−0.68−0.78
STARCH0.360.37−0.140.050.290.35−0.15−0.040.28−0.17−0.05−0.08−0.271.000.150.000.240.230.14
PROTEINS0.410.470.54−0.63−0.150.590.510.740.760.400.750.740.700.151.00−0.410.02−0.35−0.51
LIPIDS0.630.57−0.960.900.37−0.58−0.95−0.890.27−0.97−0.90−0.89−0.630.00−0.411.000.870.920.96
SUGAR0.900.87−0.820.740.38−0.36−0.84−0.650.67−0.90−0.64−0.64−0.440.240.020.871.000.870.84
LA0.680.65−0.930.900.17−0.55−0.94−0.840.34−0.96−0.85−0.87−0.680.23−0.350.920.871.000.94
AA0.550.50−0.990.980.38−0.66−1.00−0.960.17−0.99−0.94−0.94−0.780.14−0.510.960.840.941.00

Appendix B. Cluster Analysis

Table A2. Euclidean distances.
Table A2. Euclidean distances.
VariableBIOGASMETHANEDONZEAFUMAFLpHDMVSNDFADFCFLIGNINSTARCHPROTEINSLIPIDSSUGARLAAA
BIOGAS0.000.334.963.113.414.004.984.531.265.154.564.584.303.203.072.431.252.262.68
METHANE0.330.004.863.233.564.004.904.421.085.074.434.474.243.172.922.621.452.362.83
DON4.964.860.005.614.742.500.310.804.260.460.810.801.904.282.705.615.405.565.65
ZEA3.113.235.610.003.395.335.635.633.975.575.585.585.363.895.101.242.051.290.60
FUM3.413.564.743.390.003.934.704.663.854.714.624.604.853.374.303.183.163.653.14
AFL4.004.002.505.333.930.002.332.103.562.632.422.432.533.232.555.024.664.985.15
pH4.984.900.315.634.702.330.000.854.320.451.000.971.844.282.795.595.435.585.66
DM4.534.420.805.634.662.100.850.003.731.220.510.571.634.082.035.495.145.435.59
VS1.261.084.263.973.853.564.323.730.004.553.723.753.523.401.973.422.323.263.65
NDF5.155.070.465.574.712.630.451.224.550.001.251.232.164.333.095.625.525.605.64
ADF4.564.430.815.584.622.421.000.513.721.250.000.231.784.112.025.515.125.445.57
CF4.584.470.805.584.602.430.970.573.751.230.230.001.644.162.055.505.125.475.57
LIGNIN4.304.241.905.364.852.531.841.633.522.161.781.640.004.502.175.114.805.195.33
STARCH3.203.174.283.893.373.234.284.083.404.334.114.164.500.003.704.003.493.513.72
PROTEINS3.072.922.705.104.302.552.792.031.973.092.022.052.173.700.004.743.954.644.92
LIPIDS2.432.625.611.243.185.025.595.493.425.625.515.505.114.004.740.001.431.150.81
SUGAR1.251.455.402.053.164.665.435.142.325.525.125.124.803.493.951.430.001.451.59
LA2.262.365.561.293.654.985.585.433.265.605.445.475.193.514.641.151.450.000.99
AA2.682.835.650.603.145.155.665.593.655.645.575.575.333.724.920.811.590.990.00
Figure A1. Results of cluster analysis: Tree diagram for 19 variables. Single linkage. Euclidean distances. K-means clustering of variables—members of clusters. NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Figure A1. Results of cluster analysis: Tree diagram for 19 variables. Single linkage. Euclidean distances. K-means clustering of variables—members of clusters. NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Bioengineering 10 01387 g0a1
Table A3. Results of cluster analysis: K-means clustering of variables. Members of clusters and distances from respective cluster center.
Table A3. Results of cluster analysis: K-means clustering of variables. Members of clusters and distances from respective cluster center.
Cluster 1 (Biogas)Cluster 2 (ZEA)Cluster 3 (DON)Cluster 4 (AFL)Cluster 5 (FUM)
VariableDistanceVariableDistanceVariableDistanceVariableDistanceVariableDistance
1BIOGAS0.1658761ZEA0.3082241DON0.0998821AFL4.001FUM0.561778
2METHANE0.1097182LIPIDS0.2151932pH0.1518742LIGNIN4.302STARCH0.561778
3VS0.2586803SUGAR0.4044573DM0.1785733PROTEINS3.07
4LA0.2403954NDF0.245168
AA0.1518885ADF0.182380
6CF0.180448

Appendix C. Factor Analysis

Table A4. Factor loadings. Extraction: principal components. Highlighted (red-marked) factor loadings > 0.70.
Table A4. Factor loadings. Extraction: principal components. Highlighted (red-marked) factor loadings > 0.70.
VariableFactor 1Factor 2Factor 3
BIOGAS−0.58−0.810.00
METHANE−0.53−0.84−0.03
DON1.00−0.05−0.02
ZEA−0.970.19−0.10
FUM−0.39−0.060.67
AFL0.63−0.420.46
pH1.00−0.030.01
DM0.94−0.320.00
VS−0.19−0.97−0.12
NDF0.990.110.02
ADF0.93−0.31−0.02
CF0.94−0.29−0.04
LIGNIN0.78−0.36−0.32
STARCH−0.16−0.350.76
PROTEINS0.49−0.86−0.08
LIPIDS−0.96−0.08−0.13
SUGAR−0.85−0.51−0.06
LA−0.95−0.16−0.10
AA−1.000.03−0.02
Table A5. Eigenvalues. Extraction: principal components.
Table A5. Eigenvalues. Extraction: principal components.
Value NumberEigenvalue% TotalCumulative EigenvalueCumulative %
111.7168.9111.7168.91
22.6315.4714.3484.38
31.388.1115.7292.49
Figure A2. Plot of eigenvalues.
Figure A2. Plot of eigenvalues.
Bioengineering 10 01387 g0a2
Figure A3. Factor analysis—factor loadings, Factor 1 vs. Factor 2 vs. Factor 3. Extraction: principal components. VS—volatile solids; NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Figure A3. Factor analysis—factor loadings, Factor 1 vs. Factor 2 vs. Factor 3. Extraction: principal components. VS—volatile solids; NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Bioengineering 10 01387 g0a3

Appendix D. Principal Components and Classification

Figure A4. Projection of the variables on the factor plane (1 × 2) and projection of the cases (means) on the factor plane (1 × 2). Merged graphs. Total variance (Factor 1, Factor 2) values (%) of cases LC, MC, and HC are calculated from variable means. LC—low contamination; MC—medium contamination; HC—heavy contamination. VS—volatile solids; NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin. Symbol * indicates supplementary variables in biplot graph (PCA).
Figure A4. Projection of the variables on the factor plane (1 × 2) and projection of the cases (means) on the factor plane (1 × 2). Merged graphs. Total variance (Factor 1, Factor 2) values (%) of cases LC, MC, and HC are calculated from variable means. LC—low contamination; MC—medium contamination; HC—heavy contamination. VS—volatile solids; NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin. Symbol * indicates supplementary variables in biplot graph (PCA).
Bioengineering 10 01387 g0a4

References

  1. Ferrara, M.; Haidukowski, M.; D’Imperio, M.; Parente, A.; De Angelis, E.; Monaci, L.; Logrieco, A.F.; Mulè, G. New insight into microbial degradation of mycotoxins during anaerobic digestion. Waste Manag. 2021, 119, 215–225. [Google Scholar] [CrossRef] [PubMed]
  2. FAO. Agricultural Production Statistics; FAO: Rome, IT, USA, 2020. [Google Scholar]
  3. Muck, R.E.; Nadeau, E.M.G.; McAllister, T.A.; Contreras-Govea, F.E.; Santos, M.C.; Kung, L. Silage review: Recent advances and future uses of silage additives. J. Dairy Sci. 2018, 10, 3980–4000. [Google Scholar] [CrossRef] [PubMed]
  4. Binder, E.M.; Tan, L.M.; Chin, L.J.; Handl, J.; Richard, J. Worldwide occurrence of mycotoxins in commodities, feeds and feed ingredients. Anim. Feed Sci. Technol. 2007, 137, 265–282. [Google Scholar] [CrossRef]
  5. Miller, J.D.; Greenhalgh, R. Metabolites of fungal pathogens and plant resistance. Acs. Sym. Ser. 1988, 379, 117–129. [Google Scholar]
  6. Glenn, A.E. Mycotoxigenic Fusarium species in animal feed. Anim. Feed Sci. Technol. 2007, 137, 213–240. [Google Scholar] [CrossRef]
  7. Adegoke, O.G.; Letum, P. Strategies for the Prevention and Reduction of Mycotoxins in Developing Countries. In Mycotoxin and Food Safety in Developing Countries; Makun, H.A., Ed.; InTech: Rijeka, Croatia, 2013. [Google Scholar]
  8. Magan, N.; Sachis, V.; Aldred, D. Role of spoilage fungi in seed deterioration. In Fungal Biotechnology in Agricultural, Food and Environmental Applications; Aurora, D.K., Ed.; Marcell Dekker: New York, NY, USA, 2004; pp. 311–323. [Google Scholar]
  9. Balazs, E.; Schepers, J.S. The mycotoxin threat to food safety. Int. J. Food Microbiol. 2007, 119, 1–2. [Google Scholar] [CrossRef] [PubMed]
  10. Creppy, E.E. Update of survey, regulation and toxic effects of mycotoxins in Europ. Toxicol. Lett. 2002, 127, 19–28. [Google Scholar] [CrossRef]
  11. Fink-Gremmels, J.; Malekinejad, H. Clinical effects and biochemical mechanisms associated with exposure to the mycoestrogen zearalenone. Anim. Feed Sci. Technol. 2007, 137, 326–341. [Google Scholar] [CrossRef]
  12. Morgavi, D.P.; Riley, R.T. An historical overview of field disease outbreaks known or suspected to be caused by consumption of feeds contaminated with Fusarium toxins. Anim. Feed Sci. Technol. 2007, 137, 201–212. [Google Scholar] [CrossRef]
  13. Focker, M.; Van der Fels-Klerx, H.J.; Oude Lansink, A.G.J.M. Financial losses for Dutch stakeholders during the 2013 aflatoxin incident in Maize in Europe. Mycotoxin Res. 2021, 37, 193–204. [Google Scholar] [CrossRef]
  14. Huwig, A.; Freimund, S.; Kappeli, O.; Dutler, H. Mycotoxin detoxication of animal feed by different adsorbents. Toxicol. Lett. 2001, 122, 179–188. [Google Scholar] [CrossRef]
  15. Weiland, P. Biogas production: Current state and perspectives. Appl. Microbiol. Biotechnol. 2010, 85, 849–860. [Google Scholar] [PubMed]
  16. Appels, L.; Lauwers, J.; Degrève, J.; Helsen, L.; Lievens, B.; Willems, K.; Van Impe, J.; Dewil, R. Anaerobic digestion in global bioenergy production: Potential and research challenges. Renew. Sustain. Energy Rev. 2011, 15, 4295–4301. [Google Scholar] [CrossRef]
  17. Cucina, M.; Zadra, C.; Marcotullio, M.C.; Di Maria, F.; Sordi, S.; Curini, M.; Gigliotti, G. Recovery of energy and plant nutrients from a pharmaceutical organic waste derived from a fermentative biomass: Integration of anaerobic digestion and composting. J. Environ. Chem. Eng. 2017, 5, 3051–3057. [Google Scholar] [CrossRef]
  18. Tambone, F.; Scaglia, B.; D’Imporzano, G.; Schievano, A.; Orzi, V.; Salati, S.; Adani, F. Assessing amendment and fertilizing properties of digestates from anaerobic digestion through a comparative study with digested sludge and compost. Chemosphere 2010, 81, 577–583. [Google Scholar] [CrossRef] [PubMed]
  19. Liedl, B.E.; Bombardiere, J.; Williams, M.L.; Stowers, A.; Postalwait, C.; Chatfield, J.M. Solid effluent from thermophilic anaerobic digestion of poultry litter as a potential fertilizer. Hort Sci. 2004, 39, 877B–877. [Google Scholar] [CrossRef]
  20. Elbl, J.; Sláma, P.; Vaverková, M.D.; Plošek, L.; Adamcová, D.; Škarpa, P.; Kynický, J.; Havlíček, Z.; Dvořáčková, H.; Brtnický, M.; et al. Jatropa seed cake and organic waste compost: The potential improvement of soil fertility. Ecol. Chem. Eng. S 2016, 23, 131–141. [Google Scholar]
  21. Gutser, R.; Ebertseder, T.; Weber, A.; Schraml, M.; Schmidhalter, U. Short-term and residual availability of nitrogen after long-term application of organic fertilizers on arable land. J. Plant Nutr. Soil Sci. 2005, 168, 439–446. [Google Scholar] [CrossRef]
  22. Merrettig-Bruns, U.; Sayder, B. Impact of Mycotoxins and Moldy Maize Silage on the Biogas Process. Heliyon 2022. Under Review. Available online: https://ssrn.com/abstract=4084331 (accessed on 2nd October 2023). [CrossRef]
  23. Salati, S.; D’Imporzano, G.; Panseri, S.; Pasquale, E.; Adani, F. Degradation of aflatoxin B1 during anaerobic digestion and its effect on process stability. Int. Biodet. Biodeg. 2014, 94, 19–23. [Google Scholar] [CrossRef]
  24. De Gelder, L.; Audenaert, K.; Willems, B.; Schelfhout, K.; De Saeger, S.; De Boevre, M. Processing of mycotoxin contaminated waste streams through anaerobic digestion. Waste Manag. 2018, 71, 122–128. [Google Scholar] [CrossRef] [PubMed]
  25. Giorni, P.; Pietri, A.; Bertuzzi, T.; Soldano, M.; Piccinini, S.; Rossi, L.; Battilani, P. Fate of mycotoxins and related fungi in the anaerobic digestion process. Bioresour. Technol. 2018, 265, 554–557. [Google Scholar] [CrossRef] [PubMed]
  26. Vanhoutte, I.; Audenaert, K.; Gelder De, L. Biodegradation of mycotoxins: Tales from known and unexplored worlds. Front. Microbiol. 2016, 7, 561. [Google Scholar] [CrossRef]
  27. Soldano, M.; Pietri, A.; Bertuzzi, T.; Fabbri, C.; Piccinini, S.; Gallucci, F.; Aureli, G. Anaerobic digestion of mycotoxin-contaminated wheat: Effects on methane yield and contamination level. BioEnergy Res. 2021, 14, 313–321. [Google Scholar] [CrossRef]
  28. Tacconi, C.; Cucina, M.; Pezzolla, D.; Zadra, C.; Gigliotti, G. Effect of the mycotoxin aflatoxin B1 on a semi-continuous anaerobic digestion process. Waste Manag. 2018, 78, 467–473. [Google Scholar] [CrossRef] [PubMed]
  29. Lošák, T.; Válka, T.; Elbl, J.; Kintl, A.; Keutgen, A.; Keutgen, N.; Demková, L.; Árvay, J.; Varga, L.; Hnátková, H.; et al. Fertilization with Magnesium-and Sulfur-Supplemented Digestate Increases the Yield and Quality of Kohlrabi. Sustainability 2020, 12, 5733. [Google Scholar] [CrossRef]
  30. Brtnicky, M.; Kintl, A.; Holatko, J.; Hammerschmiedt, T.; Mustafa, A.; Kucerik, J.; Kucerik, J.; Vitez, T.; Prichystalova, J.; Baltazar, T.; et al. Effect of digestates derived from the fermentation of maize-legume intercropped culture and maize monoculture application on soil properties and plant biomass production. Chem. Biol. Technol. Agric. 2022, 9, 43. [Google Scholar] [CrossRef]
  31. Holatko, J.; Hammerschmiedt, T.; Kintl, A.; Danish, S.; Skarpa, P.; Latal, O.; Baltazar, T.; Fahad, S.; Akça, H.; Taban, S.; et al. Effect of carbon-enriched digestate on the microbial soil activity. PLoS ONE 2021, 16, e0274148. [Google Scholar] [CrossRef]
  32. Kolackova, I.; Smolkova, B.; Latal, O.; Skalickova, S.; Skladanka, J.; Horky, P.; Knot, P.; Hammerschmiedt, T.; Kintl, A.; Holatko, J.; et al. Does Digestate Dose Affect Fodder Security and Nutritive Value? Agriculture 2022, 12, 133. [Google Scholar] [CrossRef]
  33. ISO 20483:2013; Cereals and Pulses—Determination of the Nitrogen Content and Calculation of the Crude Protein Content—Kjeldahl Method. International Organization for Standardization (ISO): Geneva, Switzerland, 2013.
  34. Ruiz-Jiménez, J.; Priego-capote, F.; Luque de Castro, M.D. Identification and quantification of trans fatty acids in bakery products by gas chromatography–mass spectrometry after dynamic ultrasound-assisted extraction. J. Chromatogr. A 2004, 1045, 203–210. [Google Scholar] [CrossRef]
  35. ISO 6865:2000; Animal Feeding Stuffs—Determination of Crude Fibre Content—Method with Intermediate Filtration. International Organization for Standardization (ISO): Geneva, Switzerland, 2000.
  36. ISO 13906:2008; Determination of Acid Detergent Fibre (ADF) and Acid Detergent Lignin (ADL) Contents. International Organization for Standardization (ISO): Geneva, Switzerland, 2008.
  37. ISO 16472:2006; Animal Feeding Stuffs—Determination of Amylase-Treated Neutral Detergent Fibre Content. International Organization for Standardization (ISO): Geneva, Switzerland, 2006.
  38. ČSN EN 15934; Sludges, Treated Biowaste, Soils and Wastes—Calculation of Dry Matter Fraction after Determination of Dry Residue or Water Content. Czech Standardization Agency (ČSN): Praha, Czech Republic, 2013.
  39. ČSN EN 15935; Sludges, Treated Biowaste, Soils and Wastes—Determination of Loss on Ignition. Czech Standardization Agency (ČSN): Praha, Czech Republic, 2013.
  40. Hunady, I.; Ondrísková, V.; Hutyrová, H.; Kubíková, Z.; Hammerschmiedt, T.; Mezera, J. Use of Wild Plant Species: A Potential for Methane Production in Biogas Plants. Int. J. Renew. Energy Res. 2021, 11, 920–932. [Google Scholar]
  41. Kintl, A.; Huňady, I.; Holátko, J.; Vítěz, T.; Hammerschmiedt, T.; Brtnický, M.; Ondrisková, V.; Elbl, J. Using the Mixed Culture of Fodder Mallow (Malva verticillata L.) and White Sweet Clover (Melilotus albus Medik.) for Methane Production. Fermentation 2022, 8, 94. [Google Scholar] [CrossRef]
  42. ISO 6493:2000; Animal Feeding Stuffs—Determination of Starch Content—Polarimetric Method. International Organization for Standardization (ISO): Geneva, Switzerland, 2000.
  43. García-Ayuso, L.E.; Luque de Castro, M.D. A multivariate study of the performance of a microwave-assisted Soxhlet extractor for olive seeds. Anal. Chim. Acta 1999, 382, 309–316. [Google Scholar] [CrossRef]
  44. Pomeranz, Y.; Meloan, E.C. Food Analysis: Theory and Practice, 3rd ed.; Chapman and Hall: New York, NY, USA, 1994. [Google Scholar]
  45. Marrubini, G.; Papetti, A.; Genorini, E.; Ulrici, A. Determination of the sugar content in commercial plant milks by near infrared spectroscopy and Luff-Schoorl total glucose titration. Food Anal. Methods 2017, 10, 1556–1567. [Google Scholar] [CrossRef]
  46. Kintl, A.; Hammerschmiedt, T.; Vítěz, T.; Brtnický, M.; Vejražka, K.; Huňady, I.; Látal, O.; Elbl, J. Possibility of using tannins to control greenhouse gas production during digestate storage. Waste Manag. 2023, 156, 75–83. [Google Scholar] [CrossRef] [PubMed]
  47. ČSN 46 7092-42; Testing Methods for Feeding Stuffs—Part 42: Quality Testing of Silages. Czech Standardization Agency (ČSN): Praha, Czech Republic, 1998.
  48. Stringer, R. ELECTROPHORESIS|Overview. In Encyclopedia of Analytical Science; Stringer, R., Worsfold, P., Townshend, A., Poole, C., Eds.; Elsevier: Oxford, UK, 2005. [Google Scholar]
  49. Kintl, A.; Elbl, J.; Vítěz, T.; Brtnický, M.; Skládanka, J.; Hammerschmiedt, T.; Vítězová, M. Possibilities of Using White Sweetclover Grown in Mixture with Maize for Biomethane Production. Agronomy 2020, 10, 1407. [Google Scholar] [CrossRef]
  50. Kintl, A.; Vítěz, T.; Elbl, J.; Vítězová, M.; Dokulilová, T.; Nedělník, J.; Skládanka, J.; Brtnický, M. Mixed culture of corn and white lupine as an alternative to silage made from corn monoculture intended for biogas production. BioEnergy Res. 2019, 12, 694–702. [Google Scholar] [CrossRef]
  51. Skladanka, J.; Nedelnik, J.; Adam, V.; Dolezal, P.; Moravcova, H.; Dohnal, V. Forage as a Primary Source of Mycotoxins in Animal Diets. Int. J. Environ. Res. Public Health 2011, 8, 37–50. [Google Scholar] [CrossRef]
  52. Kintl, A.; Zímová, N.; Brtnický, M.; Hammerschmiedt, T.; Smutný, V.; Kincl, D.; Nerušil, P.; Huňady, I.; Elbl, J. Effect of cover crops undersown in maize on the mycotoxin content in maize biomass. Acta Fytotech. 2023, 26. in press. [Google Scholar] [CrossRef]
  53. Stefanon, B.; Procida, G. Effects of including silage in the diet on volatile compound profiles in Montasio cheese and their modification during ripening. J. Dairy Res. 2004, 71, 58–65. [Google Scholar] [CrossRef]
  54. Wilkinson, J.M. Silage; Chalcombe Publications: Southampton, UK, 2005. [Google Scholar]
  55. Kung, L.; Shaver, R.; Grant, R.; Schmidt, R. Silage review: Interpretation of chemical, microbial, and organoleptic components of silages. J. Dairy Sci. 2018, 101, 4020–4033. [Google Scholar] [CrossRef] [PubMed]
  56. González-Pereyra, M.L.; Alonso, V.A.; Sager, R.; Morlaco, M.B.; Magnoli, C.E.; Astoreca, A.L.; Rosa, C.A.R.; Chiacchiera, S.M.; Dalcero, A.M.; Cavaglieri, L.R. Fungi and select mycotoxins from pre- and postfermented corn silage. J. Appl. Microbiol. 2008, 104, 1034–1041. [Google Scholar] [CrossRef] [PubMed]
  57. Kung, L.; Shaver, R. Interpretation and use of silage fermentation analysis reports. Focus Forage 2001, 3, 20–28. [Google Scholar]
  58. Meiske, J.C.; Linn, J.G.; Goodrich, R.D. Types of laboratory silos and an evaluation of their usefulness. In Proceedings of the 2nd International Silage Research Conference; National Silo Association, Inc.: Waterloo, IA, USA, 1975; pp. 99–126. [Google Scholar]
  59. Cherney, D.; Cherney, J.H.; Cox, W.J. Fermentation Characteristics of Corn Forage Ensiled in Mini-Silos. J. Dairy Sci. 2005, 87, 42–46. [Google Scholar] [CrossRef] [PubMed]
  60. Heuzé, V.; Tran, G.; Edouard, N.; Lebas, F. Maize Silage. Animal Feed Resources Information System. Feedipedia, a Program by INRAE, CIRAD, AFZ and FAO. Last Updated on 22 June 2017. Available online: https://www.feedipedia.org/node/13883 (accessed on 12 January 2023).
  61. Massefera, D.A.; Florention, L.A.; Rabelo, C.H.S.; Härter, C.J.; Rezende, A.V.; Reis, R.A. Replacement of aruana grass by gliricidia (Gliricidia sepium) on silage quality. R. Bras. Zootec. 2015, 44, 231–239. [Google Scholar] [CrossRef]
  62. Amara, A.A.; El-Baky, N.A. Fungi as a source of edible proteins and Animal Feed. J. Fungi 2023, 9, 73. [Google Scholar] [CrossRef]
  63. Chahal, D.S. Production of protein-rich mycelial biomass of a mushroom, Pleurotus sajorcaju, on corn stover. J. Ferment. Bioeng. 1989, 68, 334–338. [Google Scholar] [CrossRef]
  64. Bakratsas, G.; Polydera, A.; Nilson, O.; Chatzikonstantinou, A.V.; Xiros, C.; Katapodis, P.; Stamatis, H. Mycoprotein production by submerged fermentation of the edible mushroom Pleurotus ostreatus in a batch stirred tank bioreactor using agro-industrial hydrolysate. Food 2023, 12, 2295. [Google Scholar] [CrossRef]
  65. Scholtmeijer, K.; van den Broek, L.A.; Fischer, A.R.H.; van Peer, A. Potential protein production from lignocellulosic materials using edible mushroom forming fungi. J. Agric. Food Chem. 2023, 71, 4450–4457. [Google Scholar] [CrossRef]
  66. García-Chávez, I.; Meraz-Romero, E.; Castelán-Ortega, O.; Zaragoza Esparza, J.; Osorio Avalos, J.; Robles Jimenez, L.E.; González-Ronquillo, M. Corn Silage, Meta-Analysis of The Quality and Yield of Different Regions in the World. Preprints 2020, 2020100094. [Google Scholar]
  67. Storm, I.M.; Kristensen, N.B.; Raun, B.M.L.; Smedsgaard, J.; Thrane, U. Dynamics in the microbiology of maize silage during whole-season storage. J. Appl. Microbiol. 2010, 109, 1017–1026. [Google Scholar] [CrossRef] [PubMed]
  68. Storm, I.M.; Rasmussen, R.R.; Rasmussen, P.H. Occurrence of Pre- and Post-Harvest Mycotoxins and Other Secondary Metabolites in Danish Maize Silage. Toxins 2014, 6, 2256–2269. [Google Scholar] [CrossRef] [PubMed]
  69. Driehuis, F.; Spanjer, M.C.; Scholten, J.; Giffel, M.C. Occurrence of mycotoxins in maize, grass and wheat silage for dairy cattle in the Netherlands. Food Addit. Contam. Part B Surveill 2008, 1, 41. [Google Scholar] [CrossRef] [PubMed]
  70. Kosicki, R.; Błajet-Kosicka, A.; Grajewski, J.; Twaruzek, M. Multiannual mycotoxin survey in feed materials and feeding stuffs. Anim. Feed Sci. Technol. 2016, 215, 165–180. [Google Scholar] [CrossRef]
  71. Storm, I.M.; Sørensen, J.L.; Rasmussen, R.R.; Nielsen, K.F. Mycotoxins in silage. Stewart Postharvest Rev. 2008, 4, 1–12. [Google Scholar]
  72. Gallo, A.; Giuberti, G.; Frisvad, J.C.; Bertuzzi, T.; Nielsen, K.F. Review on mycotoxin issues in ruminants: Occurrence in forages, effects of mycotoxin ingestion on health status and animal performance and practical strategies to counteract their negative effects. Toxins 2015, 7, 3057–3111. [Google Scholar] [CrossRef] [PubMed]
  73. Rodrigues, I.; Naehrer, K. A three-year survey on the world-wide occurrence of mycotoxins in feedstuffs and feed. Toxins 2012, 4, 663–675. [Google Scholar] [CrossRef]
  74. Schollenberger, M.; Müller, H.M.; Rüfle, M.; Suchy, S.; Plank, S.; Drochner, W. Natural occurrence of 16 Fusarium toxins ingrains and feedstuffs of plant origin from Germany. Mycopathologia 2006, 161, 43–52. [Google Scholar] [CrossRef]
  75. Vigier, B.; Reid, L.M.; Seifert, K.A.; Stewart, D.W.; Hamilton, R.I. Distribution and prediction of Fusarium species associated with maize ear rot in Ontario. Can. J. Plant Pathol. 1997, 19, 60–65. [Google Scholar] [CrossRef]
  76. Whitlow, L.W.; Hagler, W.M. Mycotoxins in dairy cattle: Occurrence, toxicity, prevention and treatment. Proc. Southwest Nutr. Conf. 2005, 124–138. [Google Scholar]
  77. Mansfield, M.A.; Wolf, E.D.; Kuladau, G.A. Relationships between weather conditions, agronomic practices, and fermentation characteristics with deoxynivalenol content in fresh and ensiled maize. Plant Dis. 2005, 89, 1151–1157. [Google Scholar] [CrossRef] [PubMed]
  78. Schmidt, P.; Novins, K.I.C.O.; Junges, D.; Almeida, R.; de Souza, C.M. Concentration of mycotoxins and chemical composition of corn silage: A farm survey using infrared thermography. J. Dairy Sci. 2015, 98, 6609–6619. [Google Scholar] [CrossRef] [PubMed]
  79. Keller, L.A.M.; González Pereyra, M.L.; Keller, K.M.; Alonso, V.A.; Oliveira, A.A.; Almeida, T.X.; Barbosa, T.S.; Nunes, L.M.T.; Cavaglieri, L.R.; Rosa, C.A.R. Fungal and mycotoxinscontamination in corn silage: Monitoring risk before and after fer-mentation. J. Stored Prod. Res. 2013, 52, 42–47. [Google Scholar] [CrossRef]
  80. Tsitsigiannis, D.I.; Dimakopoulu, M.; Antoniou, P.P.; Tjamos, E.C. Biological control strategies of mycotoxigenic fungi and associated mycotoxins in Mediterranean basin crops. Phytopathol. Mediterr. 2012, 51, 158–174. [Google Scholar]
  81. Cheli, F.; Campagnoli, A.; Dell’Orto, V. Fungal populations and mycotoxins in silages: From occurrence to analysis. Anim. Feed Sci. Technol. 2013, 183, 1–16. [Google Scholar] [CrossRef]
  82. McDonald, P.; Henderson, A.R.; Heron, S.J.E. The Bio-Chemistry of Silage, 2nd ed.; Chalcombe Publications: Marlow, UK, 1991. [Google Scholar]
  83. Baath, H.; Knabe, O.; Lepom, P. Occurrence of Fusari-um species and their mycotoxins in maize silage. Studies on the Fusarium infestation of maize silage plants. Arch. Anim. Nutr. 1990, 40, 397–405. [Google Scholar]
  84. Garon, D.; Richard, E.; Sage, L.; Bouchart, V.; Pottier, D.; Lebailly, P. Microflora and multi mycotoxin detection in corn silage: Experimental study. J. Agric. Food Chem. 2006, 54, 3479–3484. [Google Scholar] [CrossRef] [PubMed]
  85. Alonso, V.A.; Pereya, C.M.; Keller, L.A.M.; Dalcero, A.M.; Rosa, C.A.R.; Chiacchiera, S.M.; Cavaglieri, L.R. Fungi and mycotoxins in silage: An overview. J. Appl. Microbiol. 2013, 115, 637–643. [Google Scholar] [CrossRef]
  86. Zain, M.E. Impact of mycotoxins on humans and animals. J. Saudi Chem. Soc. 2011, 15, 129–144. [Google Scholar] [CrossRef]
  87. Abraham, Z.J. Fusarium Species—Their Biology and Toxicology; John Wiley & Sons: New York, NY, USA; Chichester, UK, 1986. [Google Scholar]
  88. Richter, W.I.F.; Schuster, M.; Rattenberger, E. Einfluss der fermentation von silomais auf die nachweisbarkeit von Deoxynivalenol (DON). Mycotoxin Res. 2002, 18, 16–19. [Google Scholar] [CrossRef]
  89. Rada, V.; Vlková, E. Silage Inoculants; Institute of Animal Science: Laguna, Philippines, 2010; 58p, ISBN 978-80-7403-069-7. [Google Scholar]
  90. Weinberg, Z.G.; Muck, R.E.; Weimer, P.J.; Chen, Y.; Gamburg, M. Lactic acid bacteria used in inoculants for silage as probiotics for ruminants. Appl. Biochem. Biotechnol. 2004, 118, 1–9. [Google Scholar] [CrossRef] [PubMed]
  91. Weinberg, Z.G.; Muck, R.E.; Weimer, P.J. The survival of silage inoculant lactic acid bacteria in rumen fluid. J. Appl. Microbiol. 2003, 94, 1066–1071. [Google Scholar] [CrossRef] [PubMed]
  92. Xu, Y.; Aung, M.; Sun, Z.; Zhou, Y.; Xue, T.; Cheng, X.; Cheng, Y.; Hao, L.; Zhu, W.; Degen, A. Ensiling of rice straw enhances the nutritive quality, improves average daily gain, reduces in vitro methane production and increases ruminal bacterial diversity in growing Hu lambs. Anim. Feed Sci. Technol. 2023, 295, 115513. [Google Scholar] [CrossRef]
Figure 1. Localization of experimental variants characterized by the sampling layers and the depth of silage in the silage pit, determined as follows: Variant 1 = low contamination (LC)—sampled at a depth of 50 cm; Variant 2 = medium contamination (MC); Variant 3 = heavy contamination (HC).
Figure 1. Localization of experimental variants characterized by the sampling layers and the depth of silage in the silage pit, determined as follows: Variant 1 = low contamination (LC)—sampled at a depth of 50 cm; Variant 2 = medium contamination (MC); Variant 3 = heavy contamination (HC).
Bioengineering 10 01387 g001
Figure 2. Fermentation test design (Kintl et al. [50]).
Figure 2. Fermentation test design (Kintl et al. [50]).
Bioengineering 10 01387 g002
Figure 3. Contents of DON, ZEA, and FUM mycotoxins in silage. Different letters indicate significant differences (p < 0.05) between variants (LC, MC, HC) in the specific parameter (content of individual mycotoxin).
Figure 3. Contents of DON, ZEA, and FUM mycotoxins in silage. Different letters indicate significant differences (p < 0.05) between variants (LC, MC, HC) in the specific parameter (content of individual mycotoxin).
Bioengineering 10 01387 g003
Figure 4. Contents of AFL mycotoxin in silage. Different letters indicate significant differences (p < 0.05) between variants (LC, MC, HC) in the specific parameter (content of individual mycotoxin).
Figure 4. Contents of AFL mycotoxin in silage. Different letters indicate significant differences (p < 0.05) between variants (LC, MC, HC) in the specific parameter (content of individual mycotoxin).
Bioengineering 10 01387 g004
Figure 5. Factor loadings, Factor 1 vs. Factor 2. Extraction: principal components. Members of clusters (K-means clustering). VS—volatile solids; NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid; DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Figure 5. Factor loadings, Factor 1 vs. Factor 2. Extraction: principal components. Members of clusters (K-means clustering). VS—volatile solids; NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid; DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Bioengineering 10 01387 g005
Table 1. Contents of DM, pH, VS, proteins, lipids, and sugar in silage samples and results of Tukey’s HSD test.
Table 1. Contents of DM, pH, VS, proteins, lipids, and sugar in silage samples and results of Tukey’s HSD test.
VariantDMpHVSProteinsLipidsSugar
% ± SD% ± SD% ± SD% ± SD% ± SD% ± SD
LC36.55 ± 0.48 B4.22 ± 0.02 C94.38 ± 0.02 C9.8 ± 0.01 B2.13 ± 0.09 A0.82 ± 0.03 A
MC34.31 ± 0.07 C4.51 ± 0.02 B92.05 ± 0.06 A8.42 ± 0.01 C1.97 ± 0.07 A0.48 ± 0.02 B
HC47.92 ± 0.15 A6.10 ± 0.01 A93.17 ± 0.02 B10.09 ± 0.12 A1.52 ± 0.10 B0.30 ± 0.02 C
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. SD = standard deviation; HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between respective variants in the specific parameter.
Table 2. Contents of starch, NDF, ADF, CF, and lignin in silage samples and results of Tukey’s HSD test.
Table 2. Contents of starch, NDF, ADF, CF, and lignin in silage samples and results of Tukey’s HSD test.
VariantStarchNDFADFCFLignin
% ± SD% ± SD% ± SD% ± SD% ± SD
LC28.38 ± 8.29 A40.92 ± 0.51 C24.7 ± 0.80 B20.72 ± 0.04 B3.93 ± 0.74 AB
MC27.88 ± 7.76 A44.56 ± 0.29 B22.72 ± 0.34 B19.14 ± 0.08 B3.52 ± 0.15 B
HC28.17 ± 8.05 A53.89 ± 0.61 A34.54 ± 1.79 A29.48 ± 1.82 A4.94 ± 0.03 A
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. ±SD = standard deviation; HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between respective variants in the specific parameter.
Table 3. Contents of lactic acid (LA) and acetic acid (AA), LA/AA, and titratable acidity of silage samples and results of Tukey’s HSD test.
Table 3. Contents of lactic acid (LA) and acetic acid (AA), LA/AA, and titratable acidity of silage samples and results of Tukey’s HSD test.
VariantLactic Acid (LA)Acetic Acid (AC)LA/AATitratable Acidity (TA)
% ± SD% ± SD% ± SD% ± SD
LC1.44 ± 0.05 A1.42 ± 0.02 A1.01 B2063.00 ± 24.64 A
MC1.17 ± 0.15 A1.25 ± 0.01 B0.94 B991.671 ± 29.74 B
HC0.65 ± 0.14 B0.33 ± 0.01 C1.99 A803.332 ± 25.38 C
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. ±SD = standard deviation; HSD = post hoc Tukey’s test. Different letters (A, B, and C) indicate significant differences (p < 0.05) between respective variants in the specific parameter.
Table 4. Descriptive statistics of DON and ZEA concentrations in maize silage.
Table 4. Descriptive statistics of DON and ZEA concentrations in maize silage.
VariantDON (µg/kgDM)ZEA (µg/kgDM)
Mean±SDMinMaxHSDMean±SDMinMaxHSD
LC195.715.0180.0210.0B137.71.5136.0139.0A
MC251.020.1230.0270.0B138.78.0131.0147.0A
HC638.330.1610.0670.0A69.76.563.076.0B
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. ±SD = standard deviation; HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between respective variants in the specific parameter.
Table 5. Descriptive statistics of FUM and AFL concentrations in maize silage.
Table 5. Descriptive statistics of FUM and AFL concentrations in maize silage.
VariantFUM (µg/kgDM)AFL (µg/kgDM)
Mean±SDMinMaxHSDMean±SDMinMaxHSD
LC249.710.0240.0260.0A3.00.03.03.0A
MC243.751.1194.0296.0A2.70.62.03.0A
HC225.011.5214.0237.0A3.70.63.04.0A
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. ±SD = standard deviation; HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between respective variants in the specific parameter.
Table 6. Summary of correlations between the measured variables.
Table 6. Summary of correlations between the measured variables.
VariantMean±SDDONZEAFUMAFL
DON361.67209.801.00−0.96−0.400.61
ZEA115.3334.65−0.961.000.28−0.78
FUM239.4428.88−0.400.281.000.03
AFL3.110.600.61−0.780.031.00
Legend: DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin. Values shown in red indicate a strong correlation.
Table 7. Biogas and methane yield after 21 days of fermentation in three silage variants with different degrees of contamination.
Table 7. Biogas and methane yield after 21 days of fermentation in three silage variants with different degrees of contamination.
VariantBiogas (m3/kgVS)Methane (m3/kgVS)
Mean±SDHSDMean±SDHSD
LC0.61780.0032A0.37920.0119A
MC0.43050.0113B0.22480.0015B
HC0.44640.0204B0.24790.0135B
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. ±SD = standard deviation. HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between respective variants.
Table 8. Concentration of methane in biogas after 21 days of fermentation in three silage variants with different degrees of contamination.
Table 8. Concentration of methane in biogas after 21 days of fermentation in three silage variants with different degrees of contamination.
VariantMethane Content (%vol)
Mean±SDHSD
LC630.0A
MC640.5A
HC63.50.5A
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. ±SD = standard deviation. HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between respective variants.
Table 9. Results of the correlation, cluster, factor, and principal component analyses.
Table 9. Results of the correlation, cluster, factor, and principal component analyses.
VariableCluster AnalysisCorrelation (R)Factor Analysis
Euclidean DistanceR: Biogas vs. Variables 1FactorFactor Loading 2
Members of Cluster Number 1—BIOGAS
1BIOGAS0.001.00Factor 2−0.81
2METHANE0.330.99Factor 2−0.84
3VS1.260.90Factor 2−0.97
Members of Cluster Number 2—ZEA
1ZEA3.110.39Factor 1−0.97
2LIPIDS2.430.63Factor 1−0.96
3SUGAR1.250.90Factor 1−0.85
4LA2.260.68Factor 1−0.95
5AA2.680.55Factor 1−1.00
Members of Cluster Number 3—DON
1DON4.96−0.54Factor 11.00
2pH4.98−0.55Factor 11.00
3DM4.53−0.28Factor 10.94
4NDF5.15−0.66Factor 10.99
5ADF4.56−0.30Factor 10.93
6CF4.58−0.31Factor 10.94
Members of Cluster Number 4—AFL
1AFL4.000.00Factor 10.63
2LIGNIN4.30−0.15Factor 10.78
3PROTEINS3.070.41Factor 2−0.86
Members of Cluster Number 5—FUM
1FUM3.410.27Factor 30.67
2STARCH3.200.11Factor 30.76
Legend: 1 Correlation between biogas yield and other parameters. Marked correlations are significant at p < 0.05. 2 Highlight factor (red-marked) loadings ≥ 0.70. NDF—neutral detergent fiber; ADF—acid detergent fiber; CF—crude fiber. LA—lactic acid; AA—acetic acid; DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin.
Table 10. Content of mycotoxins (DON, ZEA, FUM, AFL) in the three variants (LC, MC, HC) of silage before the AD process and in the digestate after the AD process and their difference.
Table 10. Content of mycotoxins (DON, ZEA, FUM, AFL) in the three variants (LC, MC, HC) of silage before the AD process and in the digestate after the AD process and their difference.
VariantFUM (µg/kgDM)AFL (µg/kgDM)Difference
Mean± SDMinMaxHSDMean± SDMinMaxHSDMean (µg/kgDM)Mean (%)
DON-LC195.715.0180.0210.0B250.020.0230.0270.0A54.327.8
DON-MC251.020.1230.0270.0B183.315.3170.0200.0B−67.7−27.0
DON-HC638.330.1610.0670.0A150.030.0120.0180.0B−488.3−76.6
ZEA-LC137.71.5136.0139.0A26.72.125.029.0B−111.0−80.6
ZEA-MC138.78.0131.0147.0A29.34.525.034.0AB−109.3−78.7
ZEA-HC69.76.563.076.0B36.73.533.040.0A−42.8−61.4
FUM-LC249.710.0240.0260.0A82.76.077.089.0C−167.0−66.9
FUM-MC243.751.1194.0296.0A114.73.1112.0118.0A−129.0−51.7
FUM-HC225.011.5214.0237.0A98.05.692.0103.0B−107.5−47.7
AFL-LC3.00.03.03.0A0.00.00.00.0-−3.0−100.0
AFL-MC2.70.62.03.0A0.00.00.00.0-−2.7−100.0
AFL-HC3.70.63.04.0A0.00.00.00.0-−3.7−100.0
Legend: LC—low contamination; MC—medium contamination; HC—heavy contamination. DON—deoxynivalenol; ZEA—zearalenone; FUM—fumonisin; AFL—aflatoxin. ±SD = standard deviation. HSD = post hoc Tukey’s test. Different letters indicate significant differences (p < 0.05) between the respective variants.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Kintl, A.; Vítěz, T.; Huňady, I.; Sobotková, J.; Hammerschmiedt, T.; Vítězová, M.; Brtnický, M.; Holátko, J.; Elbl, J. Effect of Mycotoxins in Silage on Biogas Production. Bioengineering 2023, 10, 1387. https://doi.org/10.3390/bioengineering10121387

AMA Style

Kintl A, Vítěz T, Huňady I, Sobotková J, Hammerschmiedt T, Vítězová M, Brtnický M, Holátko J, Elbl J. Effect of Mycotoxins in Silage on Biogas Production. Bioengineering. 2023; 10(12):1387. https://doi.org/10.3390/bioengineering10121387

Chicago/Turabian Style

Kintl, Antonín, Tomáš Vítěz, Igor Huňady, Julie Sobotková, Tereza Hammerschmiedt, Monika Vítězová, Martin Brtnický, Jiří Holátko, and Jakub Elbl. 2023. "Effect of Mycotoxins in Silage on Biogas Production" Bioengineering 10, no. 12: 1387. https://doi.org/10.3390/bioengineering10121387

APA Style

Kintl, A., Vítěz, T., Huňady, I., Sobotková, J., Hammerschmiedt, T., Vítězová, M., Brtnický, M., Holátko, J., & Elbl, J. (2023). Effect of Mycotoxins in Silage on Biogas Production. Bioengineering, 10(12), 1387. https://doi.org/10.3390/bioengineering10121387

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