Next Article in Journal
Potency Testing of Venoms and Antivenoms in Embryonated Eggs: An Ethical Alternative to Animal Testing
Next Article in Special Issue
In Vitro Mechanism Assessment of Zearalenone Removal by Plant-Derived Lactobacillus plantarum BCC 47723
Previous Article in Journal
Venom Use in Eulipotyphlans: An Evolutionary and Ecological Approach
Previous Article in Special Issue
Mycotoxin Occurrence, Toxicity, and Detoxifying Agents in Pig Production with an Emphasis on Deoxynivalenol
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Co-Occurrence of Regulated and Emerging Mycotoxins in Corn Silage: Relationships with Fermentation Quality and Bacterial Communities

1
Department of Animal Science, Food and Nutrition (DIANA), Faculty of Agricultural, Food and Environmental Sciences, Università Cattolica del Sacro Cuore, 29122 Piacenza, Italy
2
Department of Agriculture Science, University of Sassari, 07100 Sassari, Italy
3
BIOMIN Research Center, Technopark 1, 3430 Tulln, Austria
*
Author to whom correspondence should be addressed.
Submission received: 26 February 2021 / Revised: 17 March 2021 / Accepted: 19 March 2021 / Published: 23 March 2021

Abstract

:
Sixty-four corn silages were characterized for chemicals, bacterial community, and concentrations of several fungal metabolites. Silages were grouped in five clusters, based on detected mycotoxins, and they were characterized for being contaminated by (1) low levels of Aspergillus- and Penicillium-mycotoxins; (2) low levels of fumonisins and other Fusarium-mycotoxins; (3) high levels of Aspergillus-mycotoxins; (4) high levels of non-regulated Fusarium-mycotoxins; (5) high levels of fumonisins and their metabolites. Altersetin was detected in clusters 1, 3, and 5. Rugulusovin or brevianamide F were detected in several samples, with the highest concentration in cluster 3. Emodin was detected in more than 50.0% of samples of clusters 1, 3 and 5, respectively. Kojic acid occurred mainly in clusters 1 and 2 at very low concentrations. Regarding Fusarium mycotoxins, high occurrences were observed for FB3, FB4, FA1, whereas the average concentrations of FB6 and FA2 were lower than 12.4 µg/kg dry matter. Emerging Fusarium-produced mycotoxins, such as siccanol, moniliformin, equisetin, epiequisetin and bikaverin were detected in the majority of analyzed corn silages. Pestalotin, oxaline, phenopirrozin and questiomycin A were detected at high incidences. Concluding, this work highlighted that corn silages could be contaminated by a high number of regulated and emerging mycotoxins.
Keywords:
forage; mold; dairy cow; milk
Key Contribution: The current survey presents a deep characterization of mycotoxins detectable in corn silages, both for regulated and emerging mycotoxins. Abundance of strictly anaerobic bacteria might be used as biomarkers for individualizing silages characterized by low-content mycotoxins.

1. Introduction

Mycotoxins are a group of secondary metabolites produced by fungal organisms mostly belonging to the genera Aspergillus, Alternaria, Fusarium, and Penicillium strains [1,2,3]. These genera can produce different mycotoxins for which specific regulations have been established in many countries worldwide to protect consumers and livestock from their harmful effects. Approximately 100 countries have specific regulations or detailed guidelines for mycotoxins in food and feed, which are related to aflatoxins (AF) and aflatoxin M1, trichothecenes types A and B (deoxynivalenol or DON, diacetoxyscirpenol, T-2 toxin, and HT-2 toxin), fumonisins, zearalenone (ZEA), ergot alkaloids, ochratoxin A (OTA), patulin, phomopsins, and sterigmatocystin [4]. Different factors play a role in the setting of limits for a mycotoxin in a specific food category, such as availability of toxicological data and occurrence data, socio-economic issues, detailed knowledge about possibilities for adopting adequate sampling and analysis procedures. This information is often present for cereals and other concentrates as well as some of their by-products. Consequently, the regulations or recommendations are mainly issued for these feeds, whereas no legislation is available for any kinds of forages [5,6].
Some recent evidence suggests the greatest ingestion of some regulated mycotoxins by dairy cows could be mainly related to forage contamination [7,8,9], even if this aspect remains insufficiently investigated. In particular, surveys in which researchers investigated the presence of mycotoxins in home-grown forage are very limited when compared to those analyzing the problem of mycotoxin contaminations in cereals [10]. A survey by Driehuis et al. [11] revealed that the contribution of ensiled forage to mycotoxin ingestion in dairy cows on 24 farms in the Netherlands was approximately three times more than that of other feed ingredients. Furthermore, several authors [12,13,14] reviewed different studies from which an estimation of the dietary contribution of whole-plant corn silage to total mycotoxin intake in dairy cows was carried out mainly referring to regulated mycotoxins. However, many other fungal secondary metabolites besides regulated mycotoxins could be detected in forages, even if knowledge on the occurrence of these emerging and not yet regulated mycotoxins in forages appears currently very limited [15,16,17].
Even if most fungi can be eliminated during ensiling, other species such as Alternaria spp., Aspergillus fumigatus, Penicillium roqueforti, P. paneum, F. oxysporum, and Monascus ruber might be able to tolerate high levels of organic acids and carbon dioxide in addition to low availability of oxygen [6,18]. Furthermore, the presence of oxygen in some parts of silage, mainly during the first storage steps or oxygen penetration during feed-out and aerobic spoilage phases, could allow mold growth and mycotoxin production in specific areas of silage [19]. Although lactic acid bacteria are effective in hindering mold growth, just a minor increase in the oxygen concentration could provide the adequate growth conditions for these fungal strains. Indeed, when most of the acetic and lactic acids, as well as carbon dioxide decrease and more oxygen are present, nearly all cereal-associated filamentous fungi may grow [20,21].
As stated by Kung et al. [20], the evaluation of the pH and the quantification of main silage fermentation traits are the basis of evaluating silage fermentation pathways and fermentative quality. For these authors, a further step in the evaluation of silage fermentative quality is measuring other components that are commonly quantified in silages, such as nitrogenous compounds or mycotoxins. As presented previously, this is true for some conventional and regulated in other feeds mycotoxins, such as AF, fumonisins, ZEA, and DON. However, currently available data to connect quality of silage with other kinds of the so-called emerging mycotoxins, characterized by being neither routinely determined nor legislatively regulated, are very scarce [17,22,23,24].
Consequently, the aim of the current survey was to obtain data on occurrence and concentration of different regulated as well as several emerging mycotoxins detected in corn silage to provide information for a comprehensive risk assessment in this largely used forage. The hypothesis to use some microbial groups as biomarkers for characterizing silage mycotoxin contamination was also verified.

2. Results

2.1. Cluster Analysis Results

The variables used in cluster analysis were total counts of mycotoxins and concentrations of Aspergillus-, Fusarium-, Penicillium-, Alternaria-, and other mycotoxigenic fungi produced mycotoxins (Table 1). In particular, the Fusarium-produced mycotoxins were grouped as ZEA and its metabolites, trichothecenes type B, fumonisins, and their metabolites (i.e., included hidden fumonisins), enniatins, beauvericin, or other Fusarium mycotoxins. No trichothecenes type A were detected in any analyzed samples. Based on the cluster analysis, five different clusters were identified as a result of similarity in quality and quantity of mycotoxins detectable in each sample. An average root mean squared standard deviation of Mahalanobis distance among them of 1.0 was used to identify clusters. The 64 corn silages were, respectively, grouped in cluster 1 (n = 24, defined as silages contaminated by low levels of both Aspergillus- and Penicillium-produced mycotoxins), cluster 2 (n = 22, defined as silages contaminated by low levels of fumonisins, and other Fusarium-produced mycotoxins), cluster 3 (n = 2, defined as silages contaminated by high levels of Aspergillus-mycotoxins), cluster 4 (n = 9, defined as silages contaminated by high levels of Fusarium-produced mycotoxins) or cluster 5 (n = 7, defined as silages contaminated by high levels of fumonisins and their metabolites).

2.2. Alternaria Produced Mycotoxin in Corn Silages

Counts and concentrations of Alternaria-produced mycotoxins were reported in Table 2. Cluster 3 was characterized by the greatest (p < 0.05) number of Alternaria-produced mycotoxins, whereas the lowest number (i.e., 0.2; p < 0.05) was calculated for cluster 2. Altersetin was not quantified in silages belonging to clusters 2 and 4, whereas it was detected in clusters 1, 3, and 5 with incidences on total analyzed samples of 20.8, 50.0, or 28.5%. The average concentrations of altersetin in cluster 1, 3, and 4 were 6.8, 6.9, and 16.4 µg/kg dry matter or DM, respectively. The alternariol and alternariol monomethyl ether were occasionally detected in silages, with the highest incidences reported in cluster 1 (i.e., 21% and 33%, respectively) or cluster 3 (i.e., 100%). Tentoxin was not detected in cluster 3. In other clusters, tentoxin concentrations of 18.2, 46.2, 88.9, and 34.3 µg/kg DM or incidences of 4.2, 9.1, 33.3, or 14.3% were reported for cluster 1, 2, 4, and 5; respectively. Infectopyron was not detected in cluster 4.

2.3. Aspergillus Produced Mycotoxin in Corn Silages

Concerning Aspergillus produced mycotoxins (Table 2), the number of detected mycotoxins was lower (p < 0.05) in clusters 2 and 4 than clusters 3 and 5, the latter being contaminated by more than 4.0 Aspergillus-produced mycotoxins. The highest (p < 0.05) concentrations of rugulusovin was measured in cluster 3 (i.e., 310.7 µg/kg DM). Either rugulusovin or brevianamide F was detected in the majority of corn silages, with incidence values greater than 88.9%. Similar contamination levels were detected in other clusters. No aflatoxins or nigragillin were detected in any samples. Despite the concentrations of emodin appeared to be very low, being lower than 10 µg/kg DM, it was detected in more than 50.0% of samples for clusters 1, 3, and 5. Kojic acid was not detected in cluster 3, whereas concentrations ranging from 44.1 µg/kg DM to 104.0 µg/kg DM were reported for clusters 1 and 2, respectively. Incidences lower than 21% were calculated for Aspergillus produced mycotoxins such as 3-nitropropionic acid, averufin, bis(methylthio)gliotoxin, or asperphenamate. In particular, averufin was detected in one sample of cluster 1 and bis(methylthio)gliotoxin in one sample of cluster 5.

2.4. Fusarium Produced Mycotoxin in Corn Silages

Data of different regulated Fusarium produced mycotoxins are reported in Table 3.
Concerning fumonisins, clusters 1 and 2 were characterized by lower (p < 0.05) fumonisins counts than other clusters (i.e., 4.8 and 5.7 in clusters 1 and 2 with respect to 6.5, 6.7, and 7.7 in clusters 3, 4, and 5; respectively). The regulated FB1 and FB2 were detected in the majority of analyzed corn silages (i.e., occurrence values higher than 95.5%), with the greatest concentrations (p < 0.05) measured for cluster 5 (1382.2 and 285.0 µg/kg DM, respectively). Further, high fumonisins occurrences were also observed for non-regulated FB3, FB4, FA1, and masked forms of FA1, their incidences being higher than 54.2% in all clusters. The average concentrations of FB6 and FA2 resulted lower than 12.4 µg/kg DM. The greatest incidences for these two fumonisins were calculated in cluster 5 (i.e., 71.4 and 100.0%, respectively).
Concerning ZEA, this mycotoxin was detected in all clusters, with incidences ranging from 22.2% of cluster 4 to 100.0% of cluster 3. The ZEA concentrations were lower than 42 µg/kg DM without differences among clusters. The zearalenone-sulfone was detected in corn silages belonging to clusters 1, 3, and 5, respectively with concentrations of 91.2, 145.7 and 284.9 µg/kg DM or incidences of 8.3%, 100.0%, and 14.3%.
No trichothecenes type-A were detected in any samples, such as regulated T-2 and HT-2 toxins. Among trichothecenes type-B, DON was detected in 62.5%, 68.2%, 100.0%, 100.0% and 71.4% for cluster 1, 2, 3, 4 and 5; respectively. The greatest concentration was measured in cluster 3 (i.e., 151.7 µg/kg DM). Both nivalenol and deoxynivalenol-3-glucoside were detected at low concentrations and low incidences.
Concerning other than regulated Fusarium produced mycotoxins, siccanol, moniliformin, equisetin, epiequisetin, and bikaverin were detected in the majority of analyzed con silages, with incidences higher than 79.2%, except for siccanol in cluster 1 having an incidence of 58.3%. The greatest siccanol concentration (p < 0.05) was measured in cluster 4, with a value of 1,269.2 µg/kg DM. Concentrations lower than 20 µg/kg DM were observed for other mycotoxins mentioned above. Chrysogin and rubellin D were not detected in any samples, whereas monocerin was detected in all clusters at low contamination levels. The antibiotic Y was not detected in cluster 4, whereas the 15-hydroxyculmorin was detected in one corn silage belonging to cluster 1. Butenolide and apicidin were detected in clusters 3 and 5. All these aforementioned mycotoxins occurred in concentrations lower than about 20 µg/kg DM.
Enniatins counts were higher (p < 0.05) in cluster 3 than others (i.e., 4.5 vs. 1.7, respectively), Concentrations of all enniatins were low (i.e., < 15.2 µg/kg DM). Beauvericin was quantified in the majority of analyzed corn silages (i.e., incidences higher than 83.3%). Lower (p < 0.05) beauvericin was measured in clusters 1 and 2 with respect to others.

2.5. Penicillium Produced Mycotoxin in Corn Silages

Concerning Penicillium produced mycotoxins (Table 4), the clusters 3, 4 and 5 had higher counts than clusters 1 and 2 (i.e., 6.5, 5.4, 6.3 vs. 4.6 and 4.5; respectively).
Ochratoxin A and barceloneic acid were not detected in any samples, whereas 7-hydroxypestalonic, pestalotin, oxaline, phenopirrozin, and questiomycin A were detected at low concentrations (<50.0 µg/kg DM) in several corn silages, with incidences higher than 60%, except for pestalotin that was detected in 45.8% and 27.3% of corn silages belonging to cluster 1 and 2, respectively. The fellutanine A concentrations were higher (p < 0.05) in cluster 3 than others (i.e., 577.5 vs. 67.1 µg/kg DM). Mycophenolic acid, as well as its metabolite mycophenolic acid IV were detected only in one corn silage belonging to cluster 1 and one in cluster 4. Other Penicillium produced mycotoxins such as skyrin, asperglaucide, chlorocitreorosein, and cyclopenin were sporadically detected in some corn silages belonging to different clusters, even only in very low concentrations (i.e., <10 µg/kg DM).

2.6. Other Mycotoxigenic Fungi Strain Produced Mycotoxin in Corn Silages

No mycotoxins produced by fungal strains differing from the previously presented ones were detected in cluster 4, whereas the greatest number (p < 0.05) was observed for cluster 3, being 1.50 metabolites from different fungal strains. Ilicicolin E, ilicicolin C, ternatin, fungerin, calphostin, and N-Benzoyl-phenylalanine were not detected in any samples. Even when concentration levels were low (i.e., <10 µg/kg DM), these other mycotoxins were quantified mainly in some corn silages belonging to clusters 1 and 4. Macrosporin was detected in 12.5% of silages grouped in cluster 1.

2.7. Corn Silage Characteristics

In Table 5, the chemical, biological, and fermentative parameters characterizing each cluster were reported.
The chemicals did not differ among clusters and appeared in the normal range for corn silage. The starch tended to be greater (i.e., p = 0.067) in cluster 3 than others. Concerning fermentative parameters, the pH numerically tended (p = 0.156) to be higher in clusters 3, 4, and 5 (i.e., 3.84, 4.02, and 3.83, respectively) than clusters 1 and 2 (i.e., 3.67 and 3.76, respectively). Furthermore, the cluster 4 was characterized by the lowest (p = 0.057) acetic acid content (i.e., 2.54% DM), whereas the highest (p < 0.05) was measured for cluster 5 (i.e., 3.82% DM). The highest propionic acid concentration was measured for cluster 4, with values of 0.41% DM. Similar concentrations among clusters were reported for concentrations of butyric acid. Numeric differences (p = 0.246) were observed for the concentration of lactic acids, being higher than 4.00% DM in clusters 2 or lower than 3.00% DM in cluster 3. Consequently, a lactic to acetic ratio higher than 2 was calculated in cluster 2. The lowest 1, 2 propanediol concentration (p < 0.05) was measured in cluster 5 (i.e., 0.17% DM).

2.8. Corn Silage Bacterial Communities

Illumina (MiSeq) sequencing of 16S rRNA amplicons generated a data set ranging between 62,451 and 121,586 sequences per sample, after de-noising with dada2 plugin the number of sequences for sample varied between 12,088 and 55,155 raw sequences per sample with a mean length of 417.03 bp. After bioinformatics analysis, a total of 2,005,150 bacterial sequences were classified. Rarefaction analysis showed that the sequencing depth truly reflected the diversity of the microbial communities, as all the samples reached the sequencing plateau (data not shown). A total of 854 Amplicon Sequence Variances (ASVs) were found in corn silage samples after filtering and singletons removing.
Firmicutes, Proteobacteria, Bacteroidetes, and Actinobacteria phyla (Figure 1) dominated the corn silage samples with the Firmicutes phylum that overcame in abundance all other phyla.
Taxonomic bacterial community profiles at the family level (Figure 2) were dominated by Lactobacillaceae, with proportion ranging from 4.3% to 99.9%. Acetobacteraceae was the other dominant family with a proportion ranging from 0.3% to 95.6%. In 30 corn silage samples, the Lactobacillaceae related ASVs represented 90% or more of all ASVs present. The core microbiome (ASVs present in at least 75% of the samples) of corn silage was composed of four unique ASVs. The species were Lactobacillus acetotolerans, Acetobacter pasteurianus¸ Lentilactobacillus buchneri, and ASVs belong to the Lachnospiraceae family. Noteworthy, ASVs affiliated to Clostridium diolis that ferment glycerol to 1–3 propanediol [25] were found in 57% of the samples.
Based on the last approved lactic acid bacteria taxonomy, the most dominant species in the microbiome of corn silages in a decreased percentage order (Figure 3) were Lactobacillus acetotolerans, Lentilactobacillus buchneri (formerly Lactobacillus buchneri), Limosilactobacillus pontis (formerly Lactobacillus pontis), Limosilactobacillus frumenti, Lactobacillus helveticus, Limosilactobacillus panis, Lentilactobacillus farraginis/parafarraginis/diolivorans, Lactobacillus.amylovorus, Lentilactobacillus parafarraginis, Limosilactobacillus secaliphilus Limosilactobacillus reuteri, and ASVs belong to Limosilactobacillus spp. and Lactobacillus spp., only a minor number of ASVs were assigned to Limosilactobacillus coleohominis, Limosilactobacillus vaginalis, Levilactobacillus brevis, Lactobacillus hamsteri, Companilactobacillus paralimentarius, Lactobacillus delbrueckii subsp. bulgaricus, Lacticaseibacillus manihotivorans, Liquorilactobacillus vini.
The percentage of L. acetotolerans varied between 43% (in cluster 1) and 74% (in cluster 5), Limosilactobacillus pontis between 0% (cluster 5) and 33% (cluster 3), Lentilactobacillus buchneri between 0% (cluster 3) and 14% (cluster 4) and Lactobacillus helveticus between 0% (cluster 5) and 13% (cluster 1). The cluster 1 was dominated by heterofermentative LAB (L. pontis, L. panis, and Lentilactobacillus farraginis/parafarraginis/diolivorans). Conversely, the homofermentative L. acetotolerans dominated the samples of cluster 2. In the clusters 4 and 5 a higher percentage of L. buchneri respect to other clusters was detected, while an elevated percentage of the heterofermentative L. pontis was found in cluster 3.
The dominant ASVs belonging to Acetobacteraceae family were affiliated to Acetobacter pasteurianus and dominated the microbiome in only 6 samples (>50% of total reads in the sample). In cluster 5, four samples out of 7 had a percentage higher of 25% of Ac. pasteurianus. Besides Lactobacillaceae, and Acetobacteraceae, in 75% of corn silages, the ASVs found were affiliated to Lachnospiraceae family, with 57% of samples having an ASV related to Clostridium dioli. Cluster 3 was characterized by a great abundance of three ASVs (Figure 4) belonging to family Flavobacteriaceae, at genus Trabulsiella, and at species Alcaligenes faecalis (Figure 4, panel a-c). The Linear discriminant analysis Effect Size (LEfSe) algorithm allowed to identify 9 ASVs as the genomic features characterizing the differences between clusters 1 and 3 (Figure 5 and supplementary figure S1). Particularly Prevetolla, Clostridium, and Turicibacter genus and Peptostreptococcaceae family were enriched in samples inside of cluster 1, while the species Lactobacillus pontis characterized samples from cluster 3.

3. Discussion

3.1. Corn Silage Clusterization

To provide comprehensive mycotoxin occurrence data for risk assessment analysis due to corn silage intake, the contamination of several regulated or emerging mycotoxins was determined. Based on mycotoxin screening, the corn silages were grouped in five clusters characterizing by being contaminated by low levels of Aspergillus- and Penicillium-produced mycotoxins (cluster 1), contaminated by low levels of fumonisins and other Fusarium-produced mycotoxins (cluster 2), contaminated by high levels of Aspergillus- mycotoxins (cluster 3), contaminated by high levels of no regulated Fusarium-produced mycotoxins (cluster 4) or contaminated by fumonisins and their metabolites (cluster 5).

3.2. Alternaria and Aspergillus Produced Mycotoxins in Corn Silages

Alternaria spp. have been detected in fruits, vegetables, cereals, and derived products whereas few studies focused on silage, as reported in a recent review of its occurrence and role [27]. Storm et al. [28] reported low occurrences and concentrations of alternariol and alternariol monomethyl ether in forages sampled in Denmark, even if occurrence data were reported for other Alternaria produced secondary metabolites. Some Alternaria species were detected in hay and silage, as A. alternata, A. arborescens and A. tenuissima [18,29]. A. alternata was also considered as rare [30], most strains originally identified and attributed to A. alternata belonging to A. tenuissima, A. arborescens or other Alternaria species-groups. Alternariol, altertoxins, altenuene, tentoxin, and tenuazonic acid, and a broad range of compounds are produced by these species even with unconfirmed toxic properties. Escriva et al. [27] recently focused on altertoxins and macrosporin. However, Alternaria spp. can produce more than 70 different toxins including other secondary metabolites, such as 4Z-infectopyrone, phomapyrones, novae-zelandins, dehydrocurvularin, pyrenochaetic acid or alternarienonic acid produced by A. infectoria [30,31]. The main Alternaria produced mycotoxins detected were altersetin, alternariol, alternariol monomethyl ether and 4Z-Infectopyron, even only at low contamination level (<30 µg/kg DM) and low incidence. However, tentoxin contaminated 3 corn silages grouped in cluster 4 at an average concentration of 88.9 µg/kg DM. Streit et al. [23] reported also the occurrence of tentoxin, altenuene, and tenuazonic acid in feed and silages, thus confirming present findings. Alternaria produced toxins have been shown to have harmful effects in animals, including cytotoxicity, fetotoxicity, and teratogenicity, and they can cause also from hematological disorders to cancer due to their mutagenic nature [6,27]. However, toxicological data on Alternaria produced mycotoxins are limited.
The most important mycotoxins produced by Aspergillus flavus and A. parasiticus were AF and they were identified occasionally and at low levels in silage. No aflatoxins were detected in any corn silage samples, which is in line with previous findings [11,12,32]. Kojic and 3-nitropropionic acid are also produced by A. flavus and other Aspergillus species [33]. Their incidences were low in all samples, except in corn silages belonging to cluster 5. Their concentration could be considered low, being equal to 105 µg/kg DM for kojic acid and 6.8 µg/kg DM for 3-nitropropionic acid. In the Netherlands, Santos and Fink-Gremmels [34] detected in one of three sampled silages a concentration of 3-nitropropionic acid much higher of 1360 µg/kg.
Silages are often contaminated by Aspergillus fumigatus related toxins [18,35,36]. Among 226 potentially bioactive secondary metabolites produced by this mycotoxigenic fungus, gliotoxin is the most toxic one and it often reported to contaminate corn silages [37,38]. No gliotoxin was detected in any samples, whereas its metabolite bis(methylthio)gliotoxin was found only in one sample of cluster 5, which corresponds to findings from our previous survey [39,40]. Either rugulusovin or brevianamide F were detected in the majority of analyzed samples and could be produced by different mycotoxigenic filamentous fungi belonging to Aspergillus spp. and other strains. In particular, rugulusovin is a secondary metabolite of A. halophilicus [41] and brevianamide F could be produced by A. clavatus [42]. Two corn silages belonging to cluster 3 were highly contaminated by these mycotoxins, resulting in average concentrations of 310.7 and 246.4 µg/kg DM. As their toxic effect is mostly unknown in livestock, and their presence in corn silages is notable, future studies are needed to clarify the actual risk of their ingestion by farm animals. Other Aspergillus-produced mycotoxins detected in corn silages (i.e., emodin, averufin, and asperphenamate) contaminated only irregularly the analyzed samples and at very low concentrations.

3.3. Fusarium Produced Mycotoxins in Corn Silages

Fusarium mycotoxins are primarily produced by F. proliferatum and F. verticillioides and the FB1 and FB2 appeared to be the most predominant among fumonisins [43]. Their incidences were reported to be higher than 30% in corn silages sampled in North America [29,44,45], whereas lower incidences were reported in the Netherlands and France [12,46]. In this current survey, FB1 and FB2 were detected in almost all collected silages. High contamination levels of FB3 were also detected in corn silages [45,47]. Concerning other fumonisin metabolites, the FB4, FB6, FBA1, and FBA2, as well as masked fumonisin forms, the data regarding their incidences in corn silages are still scarce. As fumonisin contamination in corn silage is usually related to pre-harvest crop conditions [11,18], the specific environmental conditions characterizing the sampling years, as well as the crop management strategies before harvest, can have promoted fumonisin contaminations, with the exception for FB6 and FBA2 that showed low incidences data. Despite this, the level of fumonisin contamination appeared to be low and safe for lactating cows [9], except in cluster 5.
Different authors did not detect ZEA in corn silages [38,40], on the contrary, previous surveys often identified ZEA in corn silages, although at a concentration lower than 500 µg/kg [12,48,49,50]. These results seem to confirm present findings, in which the ZEA concentration was lower than 300 µg/kg DM for all clusters.
Concerning trichothecenes type B, they are primarily produced by F. culmorum and F. graminearum before harvest [43] and are usually detected in corn silages. In particular, DON is the main mycotoxin in silages and other forages within a broad range of concentration and occurrence [18,51], whereas nivalenol and fusarenon-X, as well as their acetylated and deacetylated analogues (3- acetyl-DON, 15-acetyl-DON and others) were detected at a less extent. Varying contamination levels of DON in forages were reported in North America and North Europe [12,28,49], sometimes exceeding 2000 µg/kg and with incidences higher than 80% [47]. On the contrary, the occurrence of nivalenol was much lower in corn silage [28,49]. These data appeared in line with those reported in a current survey, in which DON was detected in the majority of analyzed samples even at low concentrations (i.e., at most 151.7 µg/kg DM in cluster 3), whereas nivalenol occasionally contaminated corn silages independently by clusters.
Type A trichothecenes such as diacetoxyscirpenol (DAS), T2, HT2 toxins and their de-acetylated analogues from F. poae, F. sporotrichioides, and F. langsethiae [43] were often detected in corn silages [28,49,50] but were not detected in analyzed samples from this survey, thus showing the risk related to their ingestion is very low in specific sampling areas.
Among other Fusarium-produced mycotoxins, a high incidence of siccanol, moniliformin, equisetin, epiequisetin, and fusaric acid was observed. The corn silages belonging to cluster 4 were contaminated by these kinds of mycotoxins at high level, in particular regarding siccanol and fusaric acid. Streit et al. [23] reported that 29% of feed samples were contaminated with siccanol with median and maximum values of 9607 and 39,850 μg/kg, respectively. Similarly, fusaric acid was detected in 26% of feed samples with median and maximum values of 643 and 13,593 μg/kg. The fusaric acid (mean = 765 µg/kg) was detected in high prevalence and concentration in Israeli corn silage and it was defined as the toxin with the highest potential concern to dairy cow performance [52]. More recently, Ekwomadu et al. [53] found that siccanol occurred in 74% of 120 maize samples with a mean value of 64 μg/kg (i.e., range from 35 to 252 μg/kg). In the same survey, fusaric acid was found in 20% of 120 samples with a mean value of 85 μg/kg (i.e., range from 58 to 195 μg/kg). As stated by these authors, the co-occurrence of these emerging mycotoxins with other major mycotoxins and many other Fusarium metabolites of unknown toxicity is still a source of concern and should be studied in more detail. In particular, the ability of the rumen to denature these toxins might be impaired by the high passage rate in high-yielding dairy cows and mixture of toxins could end in a synergistic interaction. The presence of DON and fusaric acid in vitro caused depression of both Ruminococcus albus and Methanobrevibacter ruminantium microbial activity, whereas these effects were not observed in the presence of DON alone [54]. Custódio et al. [55] found levels of fusaric acid in corn silage used in beef cattle diets of 619 µg/kg with high incidence (i.e., 80%). The same authors highlighted that fusaric acid increased toxicity of other Fusarium-produced toxins through a synergistic mechanism [1] and this required further studies to be completely elucidated.
Similar to results reported in current survey, beauvericin and enniatins were detected at very low contamination levels in silages from Ireland and Denmark [28,56,57]. Ekwomadu et al. [53] also found an occurrence of beauvericin in 87% of 120 samples of maize grains (7.2 and 142 µg/kg for mean and max levels), whereas enniatins were detected only in less than 2% of samples.

3.4. Penicillium Produced Mycotoxins in Corn Silages

Predominant presence of post-harvest fungi in silages is associated to the presence of P. roqueforti and P. paneum [46,58,59,60]. These species belong to Penicillium section Roqueforti [61] and are able to tolerate typical silage conditions, such as high CO2, low pH, and oxygen. In particular, unfavorable weather or storage conditions during plant growth or at different ensiling phases could stimulate fungal proliferation and mycotoxin production in silage [46,62]. Nielsen et al. [63] listed the main mycotoxins produced from Penicillium strains, which are also more intensely studied and highly detected in forages. In particular, P. roqueforti was isolated from 89% of visibly moldy and from 85% of visibly un-moldy silages [64]. In the Netherlands, Driehuis et al. [12] reported an incidence of 50% for mycotoxins produced by P. roqueforti in corn silages. Among mycotoxins produced from Penicillium in silage, mycophenolic acid and roquefortines are most frequently detected [7,64,65]. On the contrary, Storm et al. [28,60] reported less than 3% of collected samples resulted being contaminated by these P. roqueforti produced mycotoxins, and these findings resulted in line with those reported in the current survey. In particular, mycophenolic acid and mycophenolic acid IV were detected in some samples of clusters 1 and 5. The concentrations appeared to be low and far to those levels of contamination associated to animal disorders [6]. Other exometabolites of Penicillium spp., such as the P. paneum biomarker marcfortine A, festuclavine, agroclavine, andrastin A, and citreoisocumarin, were also found in silages [7,28]. Some toxins related to secondary metabolism of Penicillium strains were detected in the analyzed corn silages, being secalonic acid, 7-hydroxypestalonic, pestalotin, oxaline, fellutanine A, phenopyrrozin, or questiomycin A. The concentration of these Penicillium produced mycotoxins appeared negligible, except for fellutanine A, a bio-active diketopiperazine alkaloid isolated from the cultures of Penicillium fellutanum that belongs to a class of naturally occurring 2,5-diketopiperazines. The fellutanines A−C, are non-annulated analogues of cyclo(L-Trp-L-Trp), but unlike their diannulated analogue fellutanine D are not cytotoxic. Consequently, further research is required to verify the real effects of these Penicillium related mycotoxins on animal health and on production performance.

3.5. Other Fungi Strain Produced Metabolites in Corn Silages

Illicolins are fungal metabolites produced by different strains of Fusarium spp. and other pathogenic fungi of different cereal crops, such as Ascochyta and Cylindrocladium spp. Illicolins led to cytotoxic effects on HeLa cells, but exerted only limited antibacterial activity. Subsequent studies showed that they were also toxic to various yeasts, thus showing an antibiotic effect [66]. Illicolin A and B were sporadically detected in corn silages. Ascochlorin is an isoprenoid antibiotic that is produced by the phytopathogenic fungus Ascochyta viciae [67]. Macrosporin is a bioactive metabolite produced by Stemphylium lycopersici causing leaf necrosis in plants [68]. These compounds were detected in corn silages as secondary metabolites of different fungal strains, but their effects on animal health are still unknown.

3.6. Corn Silages Chemical Characterization

The chemical composition of corn silages belonging to different clusters did not differ and the average values, as well as ranges of measurements could be considered in line with our previous survey [39,69] and with those reported in other nutritional tables. The fermentative traits slightly differed among clusters showing different fermentation quality of corn silages. In particular, average values of pH higher than 4.0 were measured for cluster 4 and this limit appears critical to identify poorly made corn silage contaminated by pathogenic bacteria [70]. For this cluster, the numeric (p > 0.158) high presence of butyric acid could be related to the metabolic activity of both saccharolytic and proteolytic Clostridium spp. [71,72]. The other clusters characterized by having a high level of mycotoxin contamination were clusters 3 and 5. Even if the pH resulted below the 4.0 threshold value, it appeared to be higher than that measured for clusters 1 and 2. Probably due to higher DM contents characterizing the two silages belonging to cluster 3, the concentrations of lactic acid were lower than for other clusters, thus reducing the lactic to acetic ratio to values below 1.0, usually defined as a threshold value for characterizing abnormal fermentation in silages [20]. For cluster 5, the high pH and the low lactic to acetic ratio as well as ammonia-nitrogen concentrations could be associated with a poor fermentation quality of these corn silages. In particular, Huhtanen et al. [73] reported the negative correlation between propionic acid and ammonia-nitrogen with digestible organic matter in lactating dairy cows. These altered traits could be probably associated with poor fermentation that produce these end-products during fermentation. Furthermore, the presence of F. verticillioides was associated with production of aldehydes and other volatile compounds [56].
The corn silages belonging to clusters 1 and 2 were less contaminated by different types of mycotoxins and their fermentation profile, characterized by low pH, low clostridia associated end-products, and high concentrations of both lactic and acetic acids, could be considered adequate [20]. These interpretations of fermentative profile appeared in line with the higher use of bacterial inoculants registered for corn silages belonging to these clusters, as well as the use of commercial products containing L. buchneri most common in cluster 1 than cluster 2. The level of ethanol measured in these corn silages did not differ with other poor preserved samples, thus confirming previous findings reported for corn silages and for which bacterial inoculants did not reduce ethanol [74,75,76,77].

3.7. Corn Silage Bacterial Communities

3.7.1. Lactobacillaceae Community in Corn Silages

Several authors showed that the dominant bacteria at the end of the fermentation phase of different types of silages belong to the family of Lactobacillaceae [78,79,80,81,82]. In general, in properly ensiled forage, facultative anaerobic bacteria gradually decrease silage pH, promoting the growth of acid-tolerant LAB, which typically dominate the eubacterial microbiome of terminal silage [72,78]. Corn silage samples analyzed in this work showed a low microbial biodiversity, and several samples were dominated by unique ASVs, and no significant differences were noticed between the five clusters in regards to the Lactobacillaceae community. This was more pronounced at the feed-out phase, as the aerobic exposure causes a decrease in microbial biodiversity [78,83]. Duniere et al. [78] observed an increase in Lactobacillaceae Operational Taxonomic Unit (out) in barley silage and hypothesized that some of the species were in a viable, but nonculturable state, with these populations increasing during aerobic exposure. Taxonomic identification of key ASVs sequences was performed by comparing the sequences to the National Center for Biotechnology Information (NCBI) database. After a blast search of dominant ASV (935,209 reads), the highest scoring hits against the sequence of the ASV (427 bp, % identity 100%) affiliated to Lactobacillus acetotolerans. This species is a facultative heterofermentative lactic acid bacteria (some strains utilize ribose), that presents a high resistance to acetic acid [26]. Xu et al. [82], sequencing the full-length 16S rRNA gene of whole crop corn silage bacteria using PacBio SMRT technology, found Lb. acetotolerans as the dominant species after 90 days of fermentation. Interestingly, the same authors also found a positive correlation between Lb. acetotolerans with 4-hydroxybutyrate, azelaic acid and 3-phenyllactic acid contents. All these compounds showed antimicrobial activity. In particular, the latter has a strong antifungal activity [84]. However, further studies will be needed to fully elucidate Lb. acetotolerans related to its antifungal activity in silages.
The latest approved taxonomy was used for the classification at species level of lactic acid bacteria [26]. The most dominant species in the microbiome of corn silage were L. acetotolerans, L. helveticus and L. buchneri. Indeed, the latter is one of the most regularly used LAB in silage inoculants due to its capacity to synthesize acetic acid even under low pH [83]. Acetic acid as a strong antifungal activity increases the aerobic stability of silage. Noteworthy, according to Drouin et al. [81], we found the Lb. helveticus as part of the common lactic bacteria community. As components of lactic bacteria community, we also found Lim. secaliphilus, Lim. panis, and Lim. pontis, which were isolated from type I and II sourdoughs, and the last species also from intestinal microbiota of swine, silage, dairy products, mezcal fermentation, and wet wheat distillers’ grains [85,86]. The habitat of Lim. secaliphilus remains unknown [87]. These species, which belong to the group of Lim. reuteri, might play a role in the metabolism of 1–2-propanediol (Sriramulu et al., 2008). Lim. panis can degrade lactic acid and acetic acid. However, Li and Nishino [88] found that Lim. panis was present in the low lactic acid samples. Other minor species found in corn silage were Len. farraginis and Len. parafarraginis, belonging to Len. buchneri group. The last species was already used as a silage inoculant and was found to be able to increase acetic acid content and improving the aerobic stability of silages in a wide range of temperature [89,90,91].

3.7.2. Acetobacteriaceae Community in Corn Silages

Concerning other bacterial communities, the presence of Acetobacteraceae is considered undesirable because it can initiate aerobic deterioration and degrade lactic acid and acetic acid to produce carbon dioxide and water [92], even if contrasting results were presented in the literature. A. pasteurianus does not produce acetic acid under anoxic conditions [93] and only in the presence of residual oxygen and ethanol, it can convert ethanol to acetic acid [94,95]. A. pasteurianus, together with species such as Klebsiella variicola, Enterobacter hormaechei, and Bacillus gibsonii, were found in deteriorated Italian ryegrass, guinea grass, and whole-crop maize silages (Li and Nishino, 2013). Li et al. (2011) found that A. pasteurianus was uniformly distributed in the three maize silage stored in bunker silos. The same authors asserted that none of the farmers acknowledged problems with these silages. Li et al. [96] suggested that acetic acid bacteria could be used to improve silage aerobic stability. Previous authors [95,96] reported that A. pasteurianus treatment did not affect the aerobic stability of corn silage. Given the ability to fixing nitrogen of this bacteria [93], it would be worth evaluating their role during the feed out on the silage microbiome.

3.7.3. Other Bacterial Family Community in Corn Silages

The differential family ASVs abundance analysis among the five clusters showed that in cluster 3, the cluster with the higher content of Aspergillus-produced mycotoxins, had a higher abundance of three ASVs belonging to the Flavobacteriaceae family (Figure 3), to the genus Trabulsiella (Figure 4) and to the species Alcaligenes faecalis (Figure 5) respectively. In particular, the family of Flavobacteriaceae is the largest family in the phylum Bacteroidetes [97] in which most are aerobic bacteria. Overall, in our samples, the ASVs abundance of this family is below 1%. Flavobacterium, a genus belonging to Flavobacteriaceae family, was found in silage because of its ability to produce cellulase [98] and release soluble carbohydrates, substrates available for other microorganisms. Alcaligenes faecalis is a gram-negative species that degrades urea and releases ammonia, thus increasing pH. This might decrease the aerobic stability of the silage and favor the fungi development. Dutkiewicz et al. [99] found that Alcaligenes faecalis is the most common species in bulk silage. Intriguing, Zhang et al. [100] reported that this bacterium has the ability to degrade ochratoxin A through the activity of N-acyl-L-amino acid amidohydrolase enzyme. Alcaligenes faecalis N1-4 isolated from tea rhizosphere soil is able to produce abundant antifungal volatiles, such as dimethyl disulfide (DMDS) and methyl isovalerate, and greatly inhibited the growth of A. flavus [101]. Therefore, the higher concentration of mycotoxin in cluster 3 samples might have stimulated the growth of these bacteria.
Trabulsiella, a genus of the family Enterobacteriaceae [102], has been found in different habitats. Trabulsiella strains have the genetic machinery for heterotrophic, aerobic metabolism and, like most close relatives, all genes for glycolysis, the tricarboxylic acid (TCA) cycle, the pentose phosphate pathway, and oxidative phosphorylation. Like their close relatives, Trabulsiella strains have the machinery for dissimilatory nitrate reduction (nitrate to nitrite to ammonia). Sapountzis et al. [103] found that T. odontotermitis, a facultative symbiont of termite’s gut, might play a role in carbohydrate metabolism and aflatoxin degradation. Suman et al. [104] found the strain of IIPTG13 Trabulsiella also might play a role in the lignin degradation. Recently, Olvera-García et al. [105] described the genome of five strains of Trabulsiella and found that this species has a pool of genes related to cellulose degradation. Maybe at the feed-out an aerobic microbial community capable of degrading cellulose settles down after the levels of soluble carbohydrates decrease.
The LEfSe analyzes allowed to detect several ASVs as biomarkers of samples of cluster 1, in particular, several strictly anaerobic species were enriched in samples of cluster 1 with respect to clusters 3, 4, and 5. The decreasing number of strictly anaerobe species in samples silage of others clusters might be due to the higher aerobic exposition of these samples, which leads to their decrease or disappearance. Therefore, the strictly anaerobic species could be used as biomarkers relate to the aerobic exposure time of silage.

3.8. Priorities for Future Investigations

Summarizing the contribution knowledge advance from this survey, a description of the presence and concentrations of emerging mycotoxins in silages has been presented. The limits of this work are mainly related to its descriptive structure not specifically linked to a focused hypothesis on the growth and occurrence of emerging mycotoxins. On the other hand, considering the little information reported in literature about these compounds, the correlation of their presence and incidence with silage characteristics and microbial communities, can be considered as a relevant contribution to the field of silage evaluation and mycotoxin risk analysis for dairy cows. More information is strongly requested to better describe the possibility to prevent the occurrence of these compounds and also to define their biological role, since the knowledge is very limited [17].
From an applicative point of view related to ruminant farms and nutritional effects, priorities of future investigation should be oriented toward two main areas: understanding of the field and silage management conditions favoring the development of emerging mycotoxins, their relations with biomass composition and environmental conditions, both at an agro-technical level and during the ensiling processes in the farm bunk. In this sense, as provided for regulated compounds [6], it is essential to deepen knowledge on the physicochemical environment of ensiling conditions favoring fungal development and mycotoxins production and moreover, on practices that could help to prevent silage contamination. Specific focus should be on the effects of pH and volatile fatty acid (VFA) concentrations in the modulation of fungal growth. As reported for regulated compounds [19], further investigations are also needed on the specific areas of the bunk that could be more affected by emerging mycotoxins.
More research is needed for a better understanding of the effects on the animal performance (intake and feed behavior, production, reproduction, health) of different ingestion doses of emerging mycotoxins, or on interactions and synergistic effects of emerging and regulated compounds that could be detrimental for animal production and human health [1].

4. Conclusions

The manuscript presents comprehensive mycotoxin occurrence data useful for risk assessment analysis due to the ingestion of contaminated corn silage. Based on mycotoxins detected in silage, different groups of contamination were defined and they were silage contaminated by low levels of Aspergillus-, Penicillium-produced mycotoxins; by low levels of fumonisins and other Fusarium-produced mycotoxins; by high levels of Aspergillus- mycotoxins; by no regulated Fusarium-produced mycotoxins or contaminated by fumonisins. More than chemicals, some fermentative traits, and different bacterial communities could be related to a different quality of these silages, thus supporting the idea that the safety of such type of forage and the concentration of these regulated and emerging mycotoxins can be modulated by different dynamics of forage colonization of specific bacterial groups. In particular, the dominant bacteria detected in all silages belong to the family of Lactobacillaceae, with the most abundant one being Lactobacillus acetolerans, followed by Acetobacteraceae. Furthermore, we observed that different not dominant bacterial communities are related to the quantity and quality of regulated and emerging mycotoxins contaminated silages. In particular, the production of organic acids with antifungal properties could improve the aerobic stability of the silages representing a positive aspect. Otherwise, some microbial groups can use organic acids as substrates for their metabolisms, favoring the aerobic deterioration of silages and mold proliferation. As a result of this survey, the use of some microbial groups as biomarkers for characterizing silage quality and mycotoxin contamination, such as strictly anaerobic bacteria, was supported even if further investigations are needed. Other important aspects are the role of good practices in silage cultivation, harvesting, and conservation to reduce the risk of fungal growth and (emerging) mycotoxin production. Finally, information is still lacking regarding the role of emerging mycotoxins, alone or in combination with regulated ones, on animal performance and their health status, which represents a novel and unexplored area of research.

5. Materials and Methods

5.1. Sample Collection, Preparation and Analysis

Sixty-four dairy farms located in the Po Valley (Italy) and Sardinia were randomly selected and visited in the 2017–2019 harvest seasons to collect corn silage samples. Corn silages were sampled at least 10 weeks after ensiling from horizontal bunker silos in agreement with the sampling procedure already described by Gallo et al. [40]. Samples of about 1.5–2 kg on a wet weight basis were sampled from at least 4 points of each silage feed-out face and analyzed as describe below.
All collected samples were split into 2 sub-samples of about 750–1000 g/sample of fresh matter. A sub-sample was treated as previously described by Gallo et al. [69] for chemical, biological, and mycotoxin analysis or by Gallo et al. [40] for characterizing fermentative traits. Briefly, samples were dried at 60 °C in a ventilated oven for 48 h until constant weight, then milled through a 1-mm screen using a laboratory mill (Thomas-Wiley, Arthur H. Thomas Co., Philadelphia, PA, USA) and stored until analysis.
All corn silages were analyzed for the presence and concentrations of fungal metabolites by LC–MS/MS at the Department of Agrobiotechnology according to Sulyok et al. [106]. As recently described by these authors, the analytical method has been extended to cover more than 500 metabolites. Briefly, 5 g of sample was weighed and extracted with 20 mL acetonitrile/water/acetic acid (79:20:1, v/v/v) for 90 min on a rotary shaker (GFL, Burgwedel, Germany). Extracts were diluted in extraction solvent (ratio 1:1) and directly injected into the LC–MS/MS instrument. All chromatographic characteristics were previously described [106]. Mycotoxins were quantified by external calibration (1/x weighted) using a multi-component standard prepared from authentic standards.
The DM was determined by gravimetric loss of free water by heating at 105 °C for 3 h (Association of Official Analytical Chemists or AOAC 1995, method 945.15); ash was determined as gravimetric residue after incineration at 550 °C for 2 h (AOAC 1995, method 942.05) and ether extract (EE) was obtained following the method 920.29 of AOAC (1995). The crude protein (CP, N × 6.25) was determined using the Kjeldahl method (AOAC 1995, method 984.13). The soluble fraction of CP (expressed on a DM basis) was determined according to Licitra et al. [107]. The neutral detergent fiber (NDF), acid detergent fiber (ADF), and lignin (ADL) were determined using the AnkomII Fiber Analyzer (Ankom Technology Corporation, Fairport, NY, USA) according to the method described by Van Soest et al. [108]. The NDF analyzes utilized a neutral detergent solution containing sodium sulfite and a heat-stable amylase (activity of 17.400 Liquefon units/ml, Ankom Technology). The NDF, ADF and ADL contents were corrected for the residual ash content. Starch was measured by polarimetry (Polax 2L, Atago®, Tokyo, Japan). The disappearance of NDF after 24 h of rumen incubation (24 h NDFD) was measured in situ by incubating nylon bags in the rumen of two cannulated dry cows for 24 h [109]. Cows were daily fed about 10 kg DM of a total mixed ration consisting of alfalfa hay, ryegrass hay, corn silage and concentrate (i.e., 250, 350, 300, and 100 g/kg DM, respectively) in two portions at 8:00 a.m. and 6:00 p.m. The diet contained 120 g CP/kg DM and 550 g NDF/kg DM.
Immediately after collection, the second sub-sample was stored as wet at 4 °C and analyzed within 24 h for acetic acid, propionic acid, butyric acid, lactic acid, NH3-N, and pH. In particular, about 50 g of wet samples were extracted using a Stomacher blender (Seward Ltd., West Sussex, UK) for 3 min in distilled water at a water-to-sample (fresh weight) ratio of 3:1. The water mixture was then filtered through a gauze and an aliquot (10 mL approximately) of the resulting solution was centrifuged at 4500× g for 15 min. The liquid supernatant was micro-filtered on a 0.45 µm syringe filter and an aliquot of micro-filtered supernatant (1.5 mL) was transferred in a vial and added with 10 µL of internal standard solution (pivalic acid at 1.5% in distilled water) for VFA gas-chromatographic analysis [76]. The analysis was carried out using a 2025 GC with Flame Ionization Detector (FID) system (Shimadzu S.r.l., Milano, Italy) equipped with auto-sampler (AOC-20i Shimadzu S.r.l., Milano, Italy) and capillary column DB-FFAP (30 m 0.250 mm; 0.25 µm; Agilent Technologies S.p.A., Milano, Italy). The analysis conditions were as follows: 200 and 230 °C injector and detector temperatures, respectively; 1.5 mL/min constant flow of hydrogen as carrier gas; 1 µL injection volume; and 30:1 split ratio. The temperature program was 40 °C 5 min, after 60 °C at 10 °C/min 5 min, after 160 °C at 5 °C/min, after 200 °C at 3 °C min. The VFA in silage samples were identified with the aid of external standards and calculated through peak areas corrected by factors of response instrumental and using pivalic acid as an internal standard. The micro-filtered extract used for volatile organic compounds (VOC) analysis was diluted (1 + 9) with 8.0 mM H2SO4 in a vial and injected into HPLC for lactic acid determination. The HPLC system consisting of PU-2080 pump, an AS-2055 sampling system, and a UV-2070 detector. An organic acid analysis column (Amine® HPX-87H Ion exchange column, 300 mm × 7.8 mm i.d., Bio-Rad, Hercules, CA, USA) was used at 35 °C with a mobile phase of 8.0 mM H2SO4 at 0.6 ml/min. The UV detector was set at 210 nm. The quantification of lactic acid was done using an external standard.
The NH3-N was determined from about 20 g of fresh samples in a slurry mix composed of 150 ml of distilled water and magnesium oxide (10 g/sample). The NH3-N content, expressed on a total nitrogen (TN) basis, was determined after steam distillation of this solution by the Kjeldahl method as reported above. Finally, the pH value was measured on the previously filtered water mixture.

5.2. Bacterial Genomic DNA Extraction and Bioinformatics

About 50 g of wet samples were homogenized as described above, using a Stomacher blender (Seward Ltd., West Sussex, UK) for 3 min in distilled water at a water-to-sample (fresh weight) ratio of 3:1. After homogenization, the sample was centrifuged for 15 min at 4600× g. The supernatant was discarded and the pellets were resuspended in 3 mL of double distilled sterile water. Aliquots of 1.5 mL of this suspension were centrifuged at 21,000× g for 3 min and the obtained pellets were used to extract the bacterial genomic DNA by the DNeasy® PowerSoil® Kit (QIAGEN, Santa Clarita, CA, USA), according to the manufacturer’s instructions. The extracted DNAs were stored at −80 °C until further analysis. The variable region V3-V4 of the 16S rRNA gene was amplified from DNA samples by PCR using the primer pair 341F (CCTACGGGNGGCWGCAG) and 805R (GACTACHVGGGTATCTAATCC) and sequenced by the 400 Illumina MiSeq platform (Paired-end, 2 × 301 bp, Macrogen, Seoul, Korea).
The Paired-End Illumina MiSeq reads (demultiplexed FASTQ files) were analyzed with the Quantitative Insights into Microbial Ecology tool (QIIME 2 version 2019.10). Sequence quality filtering and Operational Taxonomic Unit (OTU) picking were carried out using the Dada2 plugin wrapped in QIIME2 [110]. Taxonomic classification of 16S Amplicon Sequence Variant (ASV) obtained by denoising analysis was carried out using a pre-trained Naive Bayes classifier trained on the Greengenes 13_8_99% OTUs. The phylogenetic tree for diversity analysis was constructed using the q2-fragment-insertion plugin. To account for uneven sequencing depth between samples the data were normalized (the rarefaction depth was set at 12088 sequences/sample) using the diversity QIIME alpha-rarefaction plugin. Diversity analysis was carried out using the QIIME diversity core-metrics-phylogenetic plugin. To investigate the microbial diversity differences between clusters the QIIME diversity alpha-group-significance and beta-group-significance plugin were used. Data were visualized using the QIIME dedicated plugin. Linear discriminant analysis Effect Size (LEfSe) algorithm [111] was applied to identify differentially ASVs for characterizing the differences among the five clusters calculated as reported below in the Statistical Analysis section.

5.3. Statistical Analyzes

To categorize the maize silage samples into their quantity and quality of mycotoxin contents, we used a hierarchical cluster analysis using main variables related to mycotoxin contamination (i.e., total count of mycotoxins and concentrations of Aspergillus-, Fusarium-, Penicillium-, Alternaria-, and other mycotoxigenic fungi-produced mycotoxins) by the unweighted pair group mean with the arithmetic averages (UPGMA) method by the CLUSTER procedure of SAS (2003). Then, chemical, biological, fermentative measurements, mycotoxins contaminations as well as ASVs at the phylum, order, family, genus, and species levels characterizing corn silages were analyzed as a completely randomized design by using the Generalized Linear Model (GLM) procedure of SAS (2003). The fixed effect of the model was the cluster (n = 1 to 5). Mean post-hoc comparisons were performed by using Least-Squares Means (LSMEANS) option of SAS (2003), except for microbiological data for which a Benjamini-Hochberg False Discovery Rate (FDR) was used as a multiple test correction method. To understand the dynamic for individual ASVs, the sequences were manually compared to the National Center for Biotechnology Information (NCBI) database in order to annotate the phylogenic identification. Significance was declared at a p < 0.05 and, if biologically retained important, tendency at a 0.05 < p < 0.10.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/13/3/232/s1, Figure S1: relative abundance of Amplicon Sequence Variance (ASVs) as the genomic features characterizing the differences between clusters for specific bacteria.

Author Contributions

Conceptualization, A.G. and F.F.; methodology, A.G., F.G., A.S.A., S.Z., B.N., J.F., F.F.; software, A.G., F.G., S.Z., F.F.; formal analysis, A.G.; data curation, A.G., F.G., F.F.; writing—original draft preparation, A.G., A.S.A., F.F.; writing—review and editing, F.G., S.Z., B.N., J.F.; supervision, A.G., F.F.; project administration, A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

The authors are grateful to Fondazione Romeo and Enrica Invernizzi for the partial contribution to this study and to BIOMIN Holding GmbH for the assistance in mycotoxin analysis. This survey was part of Safety of Silage research project.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yiannikouris, A.; Jouany, J.-P. Mycotoxins in feeds and their fate in animals: A review. Anim. Res. 2002, 51, 81–99. [Google Scholar] [CrossRef] [Green Version]
  2. Kabak, B.; Dobson, A.D.W.; Var, I. Strategies to prevent mycotoxin contamination of food and animal feed: A review. Crit. Rev. Food Sci. Nutr. 2006, 46, 593–619. [Google Scholar] [CrossRef]
  3. Rasmussen, R.R.; Storm, I.M.L.D.; Rasmussen, P.H.; Smedsgaard, J.; Nielsen, K.F. Multi-mycotoxin analysis of maize silage by LC-MS/MS. Anal. Bioanal. Chem. 2010, 397, 765–776. [Google Scholar] [CrossRef] [Green Version]
  4. Van Egmond, H.P.; Schothorst, R.C.; Jonker, M.A. Regulations relating to mycotoxins in food Perspectives in a global and European context. Anal. Bioanal. Chem. 2007, 389, 147–157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Fink-Gremmels, J. Mycotoxins in cattle feeds and carry-over to dairy milk: A review. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess. 2008, 25, 172–180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. 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]
  7. O’Brien, M.; Nielsen, K.F.; O’Kiely, P.; Forristal, P.D.; Fuller, H.T.; Frisvad, J.C. Mycotoxins and other secondary metabolites produced in vitro by Penicillium paneum Frisvad and Penicillium roqueforti Thom isolated from baled grass silage in Ireland. J. Agric. Food Chem. 2006, 54, 9268–9276. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Driehuis, F. Silage and the safety and quality of dairy foods: A review. Agric. Food Sci. 2013, 22, 16–34. [Google Scholar] [CrossRef]
  9. Gallo, A.; Minuti, A.; Bani, P.; Bertuzzi, T.; Cappelli, F.P.; Doupovec, B.; Faas, J.; Schatzmayr, D.; Trevisi, E. A mycotoxin-deactivating feed additive counteracts the adverse effects of regular levels of Fusarium mycotoxins in dairy cows. J. Dairy Sci. 2020, 12, 11314–11331. [Google Scholar] [CrossRef] [PubMed]
  10. Gallo, A.; Rocchetti, G.; Piccioli Cappelli, F.; Pavone, S.; Mulazzi, A.; van Kuijk, S.; Han, Y.; Trevisi, E. Effect of a Commercial Bentonite Clay (Smectite Clay) on Dairy Cows Fed Aflatoxin-Contaminated Feed. Dairy 2020, 1, 135. [Google Scholar] [CrossRef]
  11. Driehuis, F.; Spanjer, M.C.; Scholten, J.M.; Te Giffel, M.C. Occurrence of mycotoxins in maize, grass and wheat silage for dairy cattle in the Netherlands. Food Addit. Contam. Part B 2008, 1, 41–50. [Google Scholar] [CrossRef]
  12. Driehuis, F.; Spanjer, M.C.; Scholten, J.M.; te Giffel, M.C. Occurrence of mycotoxins in feedstuffs of dairy cows and estimation of total dietary intakes. J. Dairy Sci. 2008, 91, 4261–4671. [Google Scholar] [CrossRef] [PubMed]
  13. 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]
  14. Ogunade, I.M.; Drouin, P.; Adesogan, A.T. Silage review: Mycotoxins in silage: Occurrence, effects, prevention, and mitigation. J. Dairy Sci. 2018, 101, 4034–4059. [Google Scholar] [CrossRef] [PubMed]
  15. Wambacq, E.; Vanhoutte, I.; Audenaert, K.; De Gelder, L.; Haesaert, G. Occurrence, prevention and remediation of toxigenic fungi and mycotoxins in silage: A review. J. Sci. Food Agric. 2016, 96, 2284–2302. [Google Scholar] [CrossRef] [PubMed]
  16. Panasiuk, L.; Jedziniak, P.; Pietruszka, K.; Piatkowska, M.; Bocian, L. Frequency and levels of regulated and emerging mycotoxins in silage in Poland. Mycotoxin Res. 2019, 35, 17–25. [Google Scholar] [CrossRef] [Green Version]
  17. Reisinger, N.; Schürer-Waldheim, S.; Mayer, E.; Debevere, S.; Antonissen, G.; Sulyok, M.; Nagl, V. Mycotoxin occurrence in maize silage-a neglected risk for bovine gut health? Toxins 2019, 11, 577. [Google Scholar] [CrossRef] [Green Version]
  18. Storm, I.M.L.D.; Sørensen, J.L.; Rasmussen, R.R.; Nielsen, K.F.; Thrane, U. Mycotoxins in silage. Stewart Postharvest Rev. 2008, 4, 1–12. [Google Scholar] [CrossRef]
  19. Weiss, K.; Kroschewski, B.; Auerbach, H. Effects of air exposure, temperature and additives on fermentation characteristics, yeast count, aerobic stability and volatile organic compounds in corn silage. J. Dairy Sci. 2016, 99, 8053–8069. [Google Scholar] [CrossRef] [Green Version]
  20. Kung, L.; Shaver, R.D.; Grant, R.J.; Schmidt, R.J. Silage review: Interpretation of chemical, microbial, and organoleptic components of silages. J. Dairy Sci. 2018, 101, 4020–4033. [Google Scholar] [CrossRef]
  21. Muck, R.E.; Nadeau, E.M.G.; Mcallister, T.A.; Santos, M.C.; Kung, L., Jr. Silage review: Recent advances and future uses of silage additives 1. J. Dairy Sci. 2018, 101, 3980–4000. [Google Scholar] [CrossRef]
  22. Streit, E.; Schatzmayr, G.; Tassis, P.; Tzika, E.; Marin, D.; Taranu, I.; Tabuc, C.; Nicolau, A.; Aprodu, I.; Puel, O.; et al. Current situation of mycotoxin contamination and co-occurrence in animal feed focus on Europe. Toxins 2012, 4, 788–809. [Google Scholar] [CrossRef] [Green Version]
  23. Streit, E.; Schwab, C.; Sulyok, M.; Naehrer, K.; Krska, R.; Schatzmayr, G. Multi-mycotoxin screening reveals the occurrence of 139 different secondary metabolites in feed and feed ingredients. Toxins 2013, 5, 504–523. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Streit, E.; Naehrer, K.; Rodrigues, I.; Schatzmayr, G. Mycotoxin occurrence in feed and feed raw materials worldwide: Long-term analysis with special focus on Europe and Asia. J. Sci. Food Agric. 2013, 93, 2892–2899. [Google Scholar] [CrossRef]
  25. Biebl, H.; Spröer, C. Taxonomy of the glycerol fermenting clostridia and description of Clostridium diolis sp. nov. Syst. Appl. Microbiol. 2002, 25, 491–497. [Google Scholar] [CrossRef]
  26. Zheng, J.; Wittouck, S.; Salvetti, E.; Franz, C.M.A.P.; Harris, H.M.B.; Mattarelli, P.; O’Toole, P.W.; Pot, B.; Vandamme, P.; Walter, J.; et al. A taxonomic note on the genus Lactobacillus: Description of 23 novel genera, emended description of the genus Lactobacillus Beijerinck 1901, and union of Lactobacillaceae and Leuconostocaceae. Int. J. Syst. Evol. Microbiol. 2020, 70, 2782–2858. [Google Scholar] [CrossRef]
  27. Escrivá, L.; Oueslati, S.; Font, G.; Manyes, L. Alternaria Mycotoxins in Food and Feed: An Overview. J. Food Quality 2017, 2017, 1569748. [Google Scholar] [CrossRef] [Green Version]
  28. Storm, I.M.L.D.; 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] [Green Version]
  29. Mansfield, M.A.; Kuldau, G.A. Microbiological and molecular determination of mycobiota in fresh and ensiled maize silage. Mycologia 2007, 99, 269–278. [Google Scholar] [CrossRef]
  30. Andersen, B.; Nielsen, K.F.; Fernández Pinto, V.; Patriarca, A. Characterization of Alternaria strains from Argentinean blueberry, tomato, walnut and wheat. Int. J. Food Microbiol. 2015, 196, 1–10. [Google Scholar] [CrossRef] [PubMed]
  31. Christensen, K.B.; Van Klink, J.W.; Weavers, R.T.; Larsen, T.O.; Andersen, B.; Phipps, R.K. Novel chemotaxonomic markers of the Alternaria infectoria species-group. J. Agric. Food Chem. 2005, 53, 9431–9435. [Google Scholar] [CrossRef]
  32. Sassahara, M.; Pontes Netto, D.; Yanaka, E.K. Aflatoxin occurrence in foodstuff supplied to dairy cattle and aflatoxin M1 in raw milk in the North of Paraná state. Food Chem. Toxicol. 2005, 43, 981–984. [Google Scholar] [CrossRef]
  33. Chang, P.-K.; Ehrlich, K.C.; Fujii, I. Cyclopiazonic acid biosynthesis of Aspergillus flavus and Aspergillus oryzae. Toxins 2009, 1, 74–99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Santos, R.R.; Fink-Gremmels, J. Mycotoxin syndrome in dairy cattle: Characterisation and intervention results. World Mycotoxin J. 2014, 7, 357–366. [Google Scholar] [CrossRef]
  35. Frisvad, J.C.; Rank, C.; Nielsen, K.F.; Larsen, T.O. Metabolomics of Aspergillus fumigatus. Med. Mycol. 2009, 47 (Suppl. S1), S53–S71. [Google Scholar] [CrossRef]
  36. Alonso, V.; Díaz Vergara, L.; Aminahuel, C.; Pereyra, C.; Pena, G.; Torres, A.; Dalcero, A.; Cavaglieri, L. Physiological behaviour of gliotoxigenic Aspergillus fumigatus sensu stricto isolated from maize silage under simulated environmental conditions. Food Addit. Contam. Part A 2015, 32, 236–244. [Google Scholar] [CrossRef]
  37. Pereyra, C.M.; Alonso, V.A.; Rosa, C.A.R.; Chiacchiera, S.M.; Dalcero, A.M.; Cavaglieri, L.R. Gliotoxin natural incidence and toxigenicity of Aspergillus fumigatus isolated from corn silage and ready dairy cattle feed. World Mycotoxin J. 2008, 1, 457–462. [Google Scholar] [CrossRef]
  38. Richard, E.; Heutte, N.; Bouchart, V.; Garon, D. Evaluation of fungal contamination and mycotoxin production in maize silage. Anim. Feed Sci. Technol. 2009, 148, 309–320. [Google Scholar] [CrossRef]
  39. Gallo, A.; Bertuzzi, T.; Giuberti, G.; Moschini, M.; Bruschi, S.; Cerioli, C.; Masoero, F. New assessment based on the use of principal factor analysis to investigate corn silage quality from nutritional traits, fermentation end products and mycotoxins. J. Sci. Food Agric. 2016, 96, 437–448. [Google Scholar] [CrossRef]
  40. Gallo, A.; Giuberti, G.; Bruschi, S.; Fortunati, P.; Masoero, F. Use of principal factor analysis to generate a corn silage fermentative quality index to rank well- or poorly preserved forages. J. Sci. Food Agric. 2016, 96, 1686–1696. [Google Scholar] [CrossRef] [PubMed]
  41. Sklenář, F.; Jurjević, Ž.; Zalar, P.; Frisvad, J.C.; Visagie, C.M.; Kolařík, M.; Houbraken, J.; Chen, A.J.; Yilmaz, N.; Seifert, K.A.; et al. Phylogeny of xerophilic aspergilli (subgenus Aspergillus) and taxonomic revision of section Restricti. Stud. Mycol. 2017, 88, 161–236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Zutz, C.; Gacek, A.; Sulyok, M.; Wagner, M.; Strauss, J.; Rychli, K. Small chemical chromatin effectors alter secondary metabolite production in Aspergillus clavatus. Toxins 2013, 5, 1723–1741. [Google Scholar] [CrossRef] [Green Version]
  43. Chelkowski, J. Fusarium: Mycotoxins, Taxonomy, Pathogenicity; Elsevier Science Publishers B. V.: Amsterdam, The Netherlands, 2014; ISBN 1483297853. [Google Scholar]
  44. Yu, W.; Yu, F.-Y.; Undersander, D.J.; Chu, F.S. Immunoassays of Selected Mycotoxins in Hay, Silage and Mixed Feed. Food Agric. Immunol. 1999, 11, 307–319. [Google Scholar] [CrossRef]
  45. Kim, S.C.; Adesogan, A.T. Influence of ensiling temperature, simulated rainfall, and delayed sealing on fermentation characteristics and aerobic stability of corn silage. J. Dairy Sci. 2006, 89, 3122–3132. [Google Scholar] [CrossRef]
  46. Richard, J.L. Some major mycotoxins and their mycotoxicoses—An overview. Int. J. Food Microbiol. 2007, 119, 3–10. [Google Scholar] [CrossRef]
  47. Mansfield, M.A. Fungi and Mycotoxins in Fresh and Ensiled Maize and the Effects of Agronomic Practices, Wheatear Conditions and Silage Characteristics. Ph.D. Thesis, Pennsylvania State University, State College, PA, USA, 2005. [Google Scholar]
  48. Cavallarin, L.; Borreani, G.; Tabacco, E.; Lúscher, A.; Jeangros, B.; Kessler, W.; Huguenin, O.; Lobsiger, M.; Millar, N.; Suter, D. Mycotoxin occurrence in farm maize silages in northern Italy. In Land Use Systems in Grassland Dominated Regions, Proceedings of the 20th General Meeting of the European Grassland Federation, Luzern, Switzerland, 21–24 June 2004; vdf Hochschulverlag AG an der ETH Zurich: Zürich, Switzerland, 2004; pp. 1023–1025. [Google Scholar]
  49. Schollenberger, M.; Müller, H.M.; Rüfle, M.; Suchy, S.; Plank, S.; Drochner, W. Natural occurrence of 16 Fusarium toxins in grains and feedstuffs of plant origin from Germany. Mycopathologia 2006, 161, 43–52. [Google Scholar] [CrossRef] [PubMed]
  50. Eckard, S.; Wettstein, F.E.; Forrer, H.R.; Vogelgsang, S. Incidence of Fusarium species and mycotoxins in silage maize. Toxins 2011, 3, 949–967. [Google Scholar] [CrossRef] [Green Version]
  51. European Food Safety Authority. Deoxynivalenol in food and feed: Occurrence and exposure. EFSA J. 2016, 11. [Google Scholar] [CrossRef]
  52. Shimshoni, J.A.; Cuneah, O.; Sulyok, M.; Krska, R.; Galon, N.; Sharir, B.; Shlosberg, A. Mycotoxins in corn and wheat silage in Israel. Food Addit. Contam. Part A Chem. Anal. Control Expo. Risk Assess. 2013, 30, 1614–1625. [Google Scholar] [CrossRef]
  53. Ekwomadu, T.I.; Dada, T.A.; Nleya, N.; Gopane, R.; Sulyok, M.; Mwanza, M. Variation of fusarium free, masked, and emerging mycotoxin metabolites in maize from Agriculture Regions of South Africa. Toxins 2020, 12, 149. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. May, H.D.; Wu, Q.; Blake, C.K. Effects of the Fusarium spp. mycotoxins fusaric acid and deoxynivalenol on the growth of Ruminococcus albus and Methanobrevibacter ruminantium. Can. J. Microbiol. 2000, 46, 692–699. [Google Scholar] [CrossRef]
  55. Custódio, L.; Prados, L.F.; Yiannikouris, A.; Holder, V.; Pettigrew, J.; Kuritza, L.; de Resende, F.D.; Siqueira, G.R. Mycotoxin contamination of diets for beef cattle finishing in feedlot. Rev. Bras. Zootec. 2019, 48, 1–12. [Google Scholar] [CrossRef]
  56. Glenn, A.E. Mycotoxigenic Fusarium species in animal feed. Anim. Feed Sci. Technol. 2007, 137, 213–240. [Google Scholar] [CrossRef]
  57. McElhinney, C.; Danaher, M.; Elliott, C.T.; O’Kiely, P. Mycotoxins in farm silages—A 2-year Irish national survey. Grass Forage Sci. 2015, 71, 339–352. [Google Scholar] [CrossRef]
  58. Nout, M.J.R.; Bouwmeester, H.M.; Haaksma, J.; Van Dijk, H. Fungal growth in silages of sugarbeet press pulp and maize. J. Agric. Sci. 1993, 121, 323–326. [Google Scholar] [CrossRef]
  59. Pereyra, M.L.G.; 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 selected mycotoxins from pre- and postfermented corn silage. J. Appl. Microbiol. 2008, 104, 1034–1041. [Google Scholar] [CrossRef]
  60. Storm, I.M.L.D.; 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]
  61. Samson, R.A.; Seifert, K.A.; Kuijpers, A.F.A.; Houbraken, J.A.M.P.; Frisvad, J.C. Phylogenetic analysis of Penicillium subgenus Penicillium using partial β-tubulin sequences. Stud. Mycol. 2004, 49, 175–200. [Google Scholar]
  62. Boudra, H.; Morgavi, D.P. Mycotoxin risk evaluation in feeds contaminated by Aspergillus fumigatus. Anim. Feed Sci. Technol. 2005, 120, 113–123. [Google Scholar] [CrossRef]
  63. Nielsen, K.F.; Sumarah, M.W.; Frisvad, J.C.; Miller, J.D. Production of metabolites from the Penicillium roqueforti complex. J. Agric. Food Chem. 2006, 54, 3756–3763. [Google Scholar] [CrossRef] [PubMed]
  64. Auerbach, H.; Maas, R.F.M.; Op Den Camp, H.J.M.; Pol, A.; Fink Gremmels, J. Biodegradation of aflatoxin B1 by bovine rumen microorganisms in vitro and its effects on rumen fermentation. In Proceedings of the Mycotox 98. Mycotoxins in Food Chain: Processing and Toxicological Aspects. Mycotox 98. Les Mycotoxines dans la Chaine Alimentaire: Aspects Toxicologiques et Technologiques, Toulouse, France, 2–4 July 1998. [Google Scholar]
  65. Mansfield, M.A.; Jones, A.D.; Kuldau, G.A. Contamination of fresh and ensiled maize by multiple penicillium mycotoxins. Phytopathology 2008, 98, 330–336. [Google Scholar] [CrossRef] [Green Version]
  66. Gutierrez-Cirlos, E.B.; Merbitz-Zahradnik, T.; Trumpower, B.L. Inhibition of the Yeast Cytochrome bc1 Complex by Ilicicolin H, a Novel Inhibitor that Acts at the Qn Site of the bc1 Complex. J. Biol. Chem. 2004, 279, 8708–8714. [Google Scholar] [CrossRef] [Green Version]
  67. Berry, E.A.; Huang, L.S.; Lee, D.W.; Daldal, F.; Nagai, K.; Minagawa, N. Ascochlorin is a novel, specific inhibitor of the mitochondrial cytochrome bc1 complex. Biochim. Biophys. Acta Bioenerg. 2010, 1797, 360–370. [Google Scholar] [CrossRef] [Green Version]
  68. Trigos, Á.; Mendoza, G.; Espinoza, C.; Salinas, A.; Fernández, J.J.; Norte, M. The role of macrosporin in necrotic spots. Phytochem. Lett. 2011, 4, 122–125. [Google Scholar] [CrossRef]
  69. Gallo, A.; Moschini, M.; Cerioli, C.; Masoero, F. Use of principal component analysis to classify forages and predict their calculated energy content. Animal 2013, 7, 930–939. [Google Scholar] [CrossRef] [Green Version]
  70. Queiroz, O.C.M.; Ogunade, I.M.; Weinberg, Z.; Adesogan, A.T. Silage review: Foodborne pathogens in silage and their mitigation by silage additives. J. Dairy Sci. 2018, 101, 4132–4142. [Google Scholar] [CrossRef]
  71. Buxton, D.R.; Muck, R.E.; Harrison, J.H.; Pahlow, G.; Muck, R.E.; Driehuis, F.; Elferink, S.J.O.; Spoelstra, S.F. Microbiology of Ensiling. Silage Sci. Technol. 2003, 42. [Google Scholar] [CrossRef]
  72. Dunière, L.; Sindou, J.; Chaucheyras-Durand, F.; Chevallier, I.; Thévenot-Sergentet, D. Silage processing and strategies to prevent persistence of undesirable microorganisms. Anim. Feed Sci. Technol. 2013, 182, 1–15. [Google Scholar] [CrossRef]
  73. Huhtanen, P.; Rinne, M.; Nousiainen, J. Evaluation of the factors affecting silage intake of dairy cows: A revision of the relative silage dry-matter intake index. Animal 2007, 1, 758–770. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Hafner, S.D.; Howard, C.; Muck, R.E.; Franco, R.B.; Montes, F.; Green, P.G.; Mitloehner, F.; Trabue, S.L.; Rotz, C.A. Emission of volatile organic compounds from silage: Compounds, sources, and implications. Atmos. Environ. 2013, 77, 827–839. [Google Scholar] [CrossRef]
  75. Bonifacio, H.F.; Rotz, C.A.; Hafner, S.D.; Montes, F.; Cohen, M.; Mitloehner, F.M. A process-based emission model of volatile organic compounds from silage sources on farms. Atmos. Environ. 2017, 152, 85–97. [Google Scholar] [CrossRef] [Green Version]
  76. Gallo, A.; Bernardes, T.F.; Copani, G.; Fortunati, P.; Giuberti, G.; Bruschi, S.; Bryan, K.A.; Nielsen, N.G.; Witt, K.L.; Masoero, F. Effect of inoculation with Lactobacillus buchneri LB1819 and Lactococcus lactis O224 on fermentation and mycotoxin production in maize silage compacted at different densities. Anim. Feed Sci. Technol. 2018, 246, 36–45. [Google Scholar] [CrossRef]
  77. Saylor, B.A.; Fernandes, T.; Sultana, H.; Gallo, A.; Ferraretto, L.F. Influence of microbial inoculation and length of storage on fermentation profile, N fractions, and ruminal in situ starch disappearance of whole-plant corn silage. Anim. Feed Sci. Technol. 2020, 267, 114557. [Google Scholar] [CrossRef]
  78. Duniere, L.; Xu, S.; Long, J.; Elekwachi, C.; Wang, Y.; Turkington, K.; Forster, R.; McAllister, T.A. Bacterial and fungal core microbiomes associated with small grain silages during ensiling and aerobic spoilage. BMC Microbiol. 2017, 17, 50. [Google Scholar] [CrossRef] [Green Version]
  79. Gharechahi, J.; Kharazian, Z.A.; Sarikhan, S.; Jouzani, G.S.; Aghdasi, M.; Hosseini Salekdeh, G. The dynamics of the bacterial communities developed in maize silage. Microb. Biotechnol. 2017, 10, 1663–1676. [Google Scholar] [CrossRef] [PubMed]
  80. Ni, K.; Minh, T.T.; Tu, T.T.M.; Tsuruta, T.; Pang, H.; Nishino, N. Comparative microbiota assessment of wilted Italian ryegrass, whole crop corn, and wilted alfalfa silage using denaturing gradient gel electrophoresis and next-generation sequencing. Appl. Microbiol. Biotechnol. 2017, 101, 1385–1394. [Google Scholar] [CrossRef] [PubMed]
  81. Drouin, P.; Tremblay, J.; Chaucheyras-Durand, F. Dynamic Succession of Microbiota during Ensiling of Whole Plant Corn Following Inoculation with Lactobacillus buchneri and Lactobacillus hilgardii Alone or in Combination. Microorganisms 2019, 7, 595. [Google Scholar] [CrossRef] [Green Version]
  82. Xu, D.; Ding, W.; Ke, W.; Li, F.; Zhang, P.; Guo, X. Modulation of metabolome and bacterial community in whole crop corn silage by inoculating homofermentative lactobacillus plantarumand heterofermentative lactobacillus buchneri. Front. Microbiol. 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  83. Ogunade, I.M.; Jiang, Y.; Pech Cervantes, A.A.; Kim, D.H.; Oliveira, A.S.; Vyas, D.; Weinberg, Z.G.; Jeong, K.C.; Adesogan, A.T. Bacterial diversity and composition of alfalfa silage as analyzed by Illumina MiSeq sequencing: Effects of Escherichia coli O157:H7 and silage additives. J. Dairy Sci. 2018, 101, 2048–2059. [Google Scholar] [CrossRef]
  84. Debonne, E.; Vermeulen, A.; Bouboutiefski, N.; Ruyssen, T.; Van Bockstaele, F.; Eeckhout, M.; Devlieghere, F. Modelling and validation of the antifungal activity of DL-3-phenyllactic acid and acetic acid on bread spoilage moulds. Food Microbiol. 2020, 88, 103407. [Google Scholar] [CrossRef] [PubMed]
  85. Vogel, R.F.; Bocker, G.; Stolz, P.; Ehrmann, M.; Fanta, D.; Ludwig, W.; Pot, B.; Kersters, K.; Schleifer, K.H.; Hammes, W.P. Identification of lactobacilli from sourdough and description of Lactobacillus pontis sp. nov. Int. J. Syst. Bacteriol. 1994, 44, 223–229. [Google Scholar] [CrossRef] [Green Version]
  86. Leser, T.D.; Amenuvor, J.Z.; Jensen, T.K.; Lindecrona, R.H.; Boye, M.; Moøller, K. Culture-independent analysis of gut bacteria: The pig gastrointestinal tract microbiota revisited. Appl. Environ. Microbiol. 2002, 68, 673–690. [Google Scholar] [CrossRef] [Green Version]
  87. Guan, H.; Ke, W.; Yan, Y.; Shuai, Y.; Li, X.; Ran, Q.; Yang, Z.; Wang, X.; Cai, Y.; Zhang, X. Screening of natural lactic acid bacteria with potential effect on silage fermentation, aerobic stability and aflatoxin B1 in hot and humid area. J. Appl. Microbiol. 2020, 128, 1301–1311. [Google Scholar] [CrossRef] [PubMed]
  88. Li, Y.; Nishino, N. Monitoring the bacterial community of maize silage stored in a bunker silo inoculated with Enterococcus faecium, Lactobacillus plantarum and Lactobacillus buchneri. J. Appl. Microbiol. 2011, 110, 1561–1570. [Google Scholar] [CrossRef] [PubMed]
  89. Endo, A.; Okada, S. Lactobacillus farraginis sp. nov. and Lactobacillus parafarraginis sp. nov., heterofermentative lactobacilli isolated from a compost of distilled shochu residue. Int. J. Syst. Evol. Microbiol. 2007, 57, 708–712. [Google Scholar] [CrossRef]
  90. Liu, Q.H.; Yang, F.Y.; Zhang, J.G.; Shao, T. Characteristics of Lactobacillus parafarraginis ZH1 and its role in improving the aerobic stability of silages. J. Appl. Microbiol. 2014, 117, 405–416. [Google Scholar] [CrossRef]
  91. Xu, Z.; He, H.; Zhang, S.; Kong, J. Effects of inoculants Lactobacillus brevis and Lactobacillus parafarraginis on the fermentation characteristics and microbial communities of corn stover silage. Sci. Rep. 2017, 7. [Google Scholar] [CrossRef] [Green Version]
  92. Spoelstra, S.F.; Courtin, M.G.; Van Beers, J.A.C. Acetic acid bacteria can initiate aerobic deterioration of whole crop maize silage. J. Agric. Sci. 1988, 111, 127–132. [Google Scholar] [CrossRef]
  93. Nanda, K.; Taniguchi, M.; Ujike, S.; Ishihara, N.; Mori, H.; Ono, H.; Murooka, Y. Characterization of acetic acid bacteria in traditional acetic acid fermentation of rice vinegar (komesu) and unpolished rice vinegar (kurosu) produced in Japan. Appl. Environ. Microbiol. 2001, 67, 986–990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Yakushi, T.; Matsushita, K. Alcohol dehydrogenase of acetic acid bacteria: Structure, mode of action, and applications in biotechnology. Appl. Microbiol. Biotechnol. 2010, 86, 1257–1265. [Google Scholar] [CrossRef]
  95. Queiroz, O.C.M.; Arriola, K.G.; Daniel, J.L.P.; Adesogan, A.T. Effects of 8 chemical and bacterial additives on the quality of corn silage. J. Dairy Sci. 2013, 96, 5836–5843. [Google Scholar] [CrossRef] [PubMed]
  96. Li, Y.; Nishino, N. Effects of ensiling fermentation and aerobic deterioration on the bacterial community in italian ryegrass, guinea grass, and whole-crop maize silages stored at high moisture content. Asian-Australas. J. Anim. Sci. 2013, 26, 1304–1312. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  97. McBride, M.J. The family flavobacteriaceae. In The Prokaryotes: Other Major Lineages of Bacteria and the Archaea; Springer: Berlin/Heidelberg, Germany, 2014; Volume 9783642389542, pp. 643–676. ISBN 9783642301230. [Google Scholar]
  98. Navarrete-Bolaños, J.L.; Jiménez-Islas, H.; Botello-Alvarez, E.; Rico-Martínez, R. Mixed culture optimization for marigold flower ensilage via experimental design and response surface methodology. J. Agric. Food Chem. 2003, 51, 2206–2211. [Google Scholar] [CrossRef] [PubMed]
  99. Dutkiewicz, J.; Olenchock, S.A.; Sorenson, W.G.; Gerencser, V.F.; May, J.J.; Pratt, D.S.; Robinson, V.A. Levels of bacteria, fungi, and endotoxin in bulk and aerosolized corn silage. Appl. Environ. Microbiol. 1989, 55, 1093–1099. [Google Scholar] [CrossRef] [Green Version]
  100. Zhang, H.H.; Wang, Y.; Zhao, C.; Wang, J.; Zhang, X.L. Biodegradation of ochratoxin A by Alcaligenes faecalis isolated from soil. J. Appl. Microbiol. 2017, 123, 661–668. [Google Scholar] [CrossRef]
  101. Gong, A.D.; Wu, N.N.; Kong, X.W.; Zhang, Y.M.; Hu, M.J.; Gong, S.J.; Dong, F.Y.; Wang, J.H.; Zhao, Z.Y.; Liao, Y.C. Inhibitory effect of volatiles emitted from alcaligenes faecalis N1-4 on aspergillus flavusand aflatoxins in storage. Front. Microbiol. 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  102. Mcwhorter, A.C.; Haddock, R.L.; Nocon, F.A.; Steigerwalt, A.G.; Brenner, D.J.; Aleksit, S.; Bockemuhl, J.; Farmer, J.J. Trabulsiella guamensis, a new genus and species of the family Enterobacteriaceae that resembles Salmonella subgroups 4 and 5. J. Clin. Microbiol. 1991, 29, 1480–1485. [Google Scholar] [CrossRef] [Green Version]
  103. Sapountzis, P.; Gruntjes, T.; Otani, S.; Estevez, J.; Da Costa, R.R.; Iii, G.P.; Perna, N.T.; Poulsen, M. The Enterobacterium Trabulsiella odontotermitis Presents Novel Adaptations Related to Its Association with Fungus-Growing Termites. Appl. Environ. Microbiol. 2015, 81, 6577–6588. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Suman, S.K.; Dhawaria, M.; Tripathi, D.; Raturi, V.; Adhikari, D.K.; Kanaujia, P.K. Investigation of lignin biodegradation by Trabulsiella sp. isolated from termite gut. Int. Biodeterior. Biodegrad. 2016, 112, 12–17. [Google Scholar] [CrossRef]
  105. Olvera-García, M.; Fontes-Perez, H.; Chávez-Martínez, A.; Barrera, O.R.; Rodríguez-Almeida, F.A.; Sanchez-Flores, A.; Corral-Luna, A. Draft genome sequences for five strains of Trabulsiella odontotermitis, isolated from Heterotermes sp. termite gut. Genome Announc. 2015, 3, e01289-15. [Google Scholar] [CrossRef] [Green Version]
  106. Sulyok, M.; Stadler, D.; Steiner, D.; Krska, R. Validation of an LC-MS/MS-based dilute-and-shoot approach for the quantification of >500 mycotoxins and other secondary metabolites in food crops: Challenges and solutions. Anal. Bioanal. Chem. 2020, 412, 2607–2620. [Google Scholar] [CrossRef] [Green Version]
  107. Licitra, G.; Hernandez, T.M.; Van Soest, P.J. Standardization of procedures for nitrogen fractionation of ruminant feeds. Anim. Feed Sci. Technol. 1996, 57, 347–358. [Google Scholar] [CrossRef]
  108. Van Soest, P.J.; Robertson, J.B.; Lewis, B.A. Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 1991, 74, 3583–3597. [Google Scholar] [CrossRef]
  109. Gallo, A.; Giuberti, G.; Bruschi, S.; Fortunati, P.; Masoero, F. Technical note: Relationship between in situ NDF degradability and enzymatic NDF hydrolysis in forages, nonforage fibrous feeds, and crop residues. J. Anim. Sci. 2017, 95, 4172–4180. [Google Scholar] [CrossRef] [PubMed]
  110. Bokulich, N.A.; Kaehler, B.D.; Rideout, J.R.; Dillon, M.; Bolyen, E.; Knight, R.; Huttley, G.A.; Gregory Caporaso, J. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome 2018, 6, 1–17. [Google Scholar] [CrossRef] [PubMed]
  111. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Relative abundance at phylum level of bacterial communities of the corn silage samples grouped in the five different clusters (based on the content and type of mycotoxins).
Figure 1. Relative abundance at phylum level of bacterial communities of the corn silage samples grouped in the five different clusters (based on the content and type of mycotoxins).
Toxins 13 00232 g001
Figure 2. Relative abundance of the top 15 families of bacterial communities of the corn silage samples grouped in the five different clusters (based on the content and type of mycotoxins).
Figure 2. Relative abundance of the top 15 families of bacterial communities of the corn silage samples grouped in the five different clusters (based on the content and type of mycotoxins).
Toxins 13 00232 g002
Figure 3. Relative abundance of the lactic acid bacteria (LAB) species of the corn silage samples grouped in the five different clusters (based to the content and type of mycotoxins). The species of LAB were named according to the new classification proposed by Zheng et al. [26].
Figure 3. Relative abundance of the lactic acid bacteria (LAB) species of the corn silage samples grouped in the five different clusters (based to the content and type of mycotoxins). The species of LAB were named according to the new classification proposed by Zheng et al. [26].
Toxins 13 00232 g003
Figure 4. ANOVA analysis with Benjamini–Hochberg False Discovery Rate (FDR) used as multiple test correction method revealed the absolute abundance. The abundances of three Amplicon Sequence Variance (ASVs) (belong to the following taxa: Trabulsiella, Flavobact, and Flavobacteriaceae) were significantly different (p < 0.05).
Figure 4. ANOVA analysis with Benjamini–Hochberg False Discovery Rate (FDR) used as multiple test correction method revealed the absolute abundance. The abundances of three Amplicon Sequence Variance (ASVs) (belong to the following taxa: Trabulsiella, Flavobact, and Flavobacteriaceae) were significantly different (p < 0.05).
Toxins 13 00232 g004
Figure 5. Linear discriminant analysis (LDA) effect size (LEfSe) analysis revealed significant Amplicon Sequence Variances (ASVs) differences in silage microbiota between cluster 1 (red bar) and cluster 3 (green bar). LDA scores (log10) > 2 and p < 0.05 are listed.
Figure 5. Linear discriminant analysis (LDA) effect size (LEfSe) analysis revealed significant Amplicon Sequence Variances (ASVs) differences in silage microbiota between cluster 1 (red bar) and cluster 3 (green bar). LDA scores (log10) > 2 and p < 0.05 are listed.
Toxins 13 00232 g005
Table 1. Counts (n) and sums (µg/kg dry matter or DM) of mycotoxins of corn silages belonging to different clusters.
Table 1. Counts (n) and sums (µg/kg dry matter or DM) of mycotoxins of corn silages belonging to different clusters.
ItemsCluster 1Cluster 2Cluster 3Cluster 4Cluster 5√MSEp Value
n = 24n = 22n = 2n = 9n = 7
Counts of mycotoxins24.723.542.525.432.75.93<0.05
Aspergillus toxins3.12.64.02.24.10.99<0.05
Alternaria toxins1.00.22.50.31.11.07<0.05
Zearalenoneand its metabolites0.40.22.00.20.60.55<0.05
Trichothecenes type B0.80.71.51.00.90.560.256
Fumonisins and their metabolites4.85.86.56.77.71.46<0.05
Enniatins0.80.33.50.21.01.18<0.05
Beauvericin0.81.01.01.01.00.240.133
Other Fusarium toxins6.56.911.57.28.91.62<0.05
Penicillium toxins4.64.56.55.46.31.10<0.05
Other fungi toxins0.60.11.50.00.90.970.103
Unspecified fungi toxins0.80.13.00.00.70.79<0.05
Sums of mycotoxins
Aspergillus toxins147.084.5565.270.3186.7104.04<0.05
Alternaria toxins5.84.418.729.618.732.670.308
Zearalenoneand its metabolites8.84.0152.80.541.446.27<0.05
Trichothecenes type B28.815.4192.633.557.641.67<0.05
FB and their metabolites215.4339.1475.3473.51944.9289.56<0.05
Enniatins0.60.33.10.55.74.460.075
Beauvericin4.18.530.819.727.113.15<0.05
Other Fusarium toxins229.9755.3619.71564.8675.1172.65<0.05
Penicillium toxins154.691.6708.287.3142.2107.34<0.05
Other fungi toxins1.10.14.30.04.02.850.013
Unspecified fungi toxins17.81.8102.00.026.023.51<0.05
√MSE: root mean square error. When not detectable, the limit of detection of specific mycotoxins was used to compute statistical analysis.
Table 2. Concentration 1 (µg/kg DM) and incidence (%, within parenthesis) of Alternaria- and Aspergillus-produced mycotoxin in corn silages belonging to different clusters.
Table 2. Concentration 1 (µg/kg DM) and incidence (%, within parenthesis) of Alternaria- and Aspergillus-produced mycotoxin in corn silages belonging to different clusters.
ItemsCluster 1Cluster 2Cluster 3Cluster 4Cluster 5√MSEp Value
n = 24n = 22n = 2n = 9n = 7
Alternaria mycotoxin count, n10.22.50.31.11.070.007
Altersetin6.8 (20.8)nd (-)6.9 (50.0)nd (-)16.4 (28.6)7.40.358
Alternariol1.9 (20.8)nd (-)1.2 (100)nd (-)nd (-)1.410.659
Alternariol, methyl-ether1.3 (33.3)0.6 (4.5)0.9 (100)nd (-)2.9 (14.3)0.840.267
Tentoxin18.2 (4.2)46.2 (9.1)nd (-)88.9 (33.3)34.3 (14.3)111.550.929
4Z-Infectopyron15.6 (16.7)3.9 (4.5)27.5 (50.0)nd15.3 (57.1)12.040.617
Aspergillus mycotoxin count, n3.12.642.24.10.99<0.05
Rugulusovin65 (100)31.1 (100)310.7 (100)24.6 (100)34.7 (100)53.56<0.05
Emodin4.9 (50)1.6 (31.8)7.2 (100)1 (11.1)3.7 (57.1)2.790.066
Brevianamide F74.3 (95.8)30.6 (95.5)246.4 (100)40.5 (88.9)100.1 (100)65.65<0.05
Kojic acid44.1 (16.7)104 (22.7)nd (-)83.5 (11.1)105.0 (42.9)86.310.734
3-Nitropropionic acid5.2 (12.5)1.3 (4.5)nd (-)2 (11.1)6.8 (57.1)6.250.824
Averufin0.3 (4.2)nd (-)nd (-)nd (-)nd (-)--
Bis(methylthio)gliotoxinnd (-)nd (-)nd (-)nd (-)3 (14.3)--
Asperphenamatend (-)0.3 (4.5)nd (-)nd (-)1.7 (14.3)--
nd: mycotoxin not detectable in any samples. √MSE: root mean square error. When not detectable, the limit of detection of specific mycotoxins was used to compute statistical analysis. 1 Concentration represented the average values of only detectable mycotoxins.
Table 3. Concentration 1 (µg/kg DM) and incidence (%, within parenthesis) of regulated and emerging Fusarium-produced mycotoxin in corn silages belonging to different clusters.
Table 3. Concentration 1 (µg/kg DM) and incidence (%, within parenthesis) of regulated and emerging Fusarium-produced mycotoxin in corn silages belonging to different clusters.
ItemsCluster 1Cluster 2Cluster 3Cluster 4Cluster 5√MSEp Value
n = 24n = 22n = 2n = 9n = 7
Fumonisins mycotoxins count, n4.85.76.56.77.71.46<0.05
Fumonisin B1 2155 (95.8)228.7 (100)181.8 (100)321.9 (100)1382.2 (100)216.24<0.05
Fumonisin B2 236.6 (95.8)67.7 (95.5)139.1 (100)96.1 (100)285 (100)74.58<0.05
Fumonisin B316.3 (75)19.1 (100)69 (100)28.2 (100)149.6 (100)39.05<0.05
Fumonisin B410.1 (70.8)17 (95.5)28.9 (100)14.9 (100)68.6 (100)21.14<0.05
Fumonisin B62.9 (8.3)4.2 (13.6)nd (-)2.5 (22.2)4.2 (71.4)1.730.607
Fumonisin A14 (54.2)4.4 (77.3)9.9 (100)4.5 (100)18.1 (100)6.78<0.05
Fumonisin A22.5 (16.7)4.3 (40.9)5.6 (50)2.9 (66.7)12.4 (100)3.45<0.05
Masked Fumonisin A116.3 (58.3)8.8 (54.5)43.8 (100)7 (77.8)26 (100)13.66<0.05
Zearalenone metabolites count, n0.40.220.20.60.550.001
Zearalenone 24.1 (29.2)22 (18.2)7.1 (100)2.5 (22.2)1.7 (42.9)19.040.579
Zearalenone-sulfone91.2 (8.3)nd (-)145.7 (100)nd (-)284.9 (14.3)132.850.584
Trichothecenes type-B count, n0.80.71.510.90.560.256
Deoxynivalenol 238.4 (62.5)22.6 (68.2)151.7 (100)33.5 (100)73.8 (71.4)41.690.002
Nivalenol55 (8.3)nd (-)81.9 (50)nd (-)34.4 (14.3)27.730.635
Deoxynivalenol-3-glucoside6.7 (4.2)ndndndnd--
other Fusarium mycotoxins count, n6.56.911.57.28.91.62<0.05
Siccanol30 (58.3)464.4 (100)210 (100)1269.2 (100)255.8 (100)208.26<0.05
Monocerin3.9 (16.7)5.4 (40.9)19.9 (100)5.2 (33.3)3.3 (57.1)6.7260.088
Moniliformin10.9 (95.8)9 (95.5)39.9 (100)16.8 (100)39 (100)16.53<0.05
Equisetin8.7 (91.7)17.7 (90.9)141.1 (100)5.4 (100)8.9 (100)29.09<0.05
Epiequisetin8.4 (79.2)15.1 (90.9)123.4 (100)5.1 (100)9.4 (100)24.48<0.05
Culmorin25.7 (42.7)14.3 (50)47.5 (100)14.6 (55.6)44.2 (71.4)19.46<0.05
15-Hydroxyculmorin4.7 (4.2)nd (-)nd (-)nd (-)nd (-)--
Butenolide12.9 (33.3)nd (-)17.8 (50)nd (-)17.9 (28.6)1.92<0.05
Bikaverin10 (91.7)18.1 (90.9)19.4 (100)15.2 (100)49.4 (100)12.75<0.05
Apicidinnd (-)nd (-)4.2 (50)nd (-)1.3 (14.3)--
Antibiotic Y11 (12.5)1.3 (4.5)15.1 (50)nd (-)18.5 (28.6)15.70.833
7-Hydroxykaurenolide2 (8.3)2.4 (4.5)7.8 (100)nd (-)3.5 (14.3)4.760.681
Aurofusarin8.5 (16.7)7 (22.7)5.5 (100)6.3 (33.3)6 (71.4)5.680.963
Fusaric acid159.3 (100)225.0 (100)195.5 (100)241 (100)263.7 (100)155.940.428
Enniatins-Beauvericin count, n1.71.34.51.221.21<0.05
Enniatin A3.2 (8.3)0.4 (4.5)nd (-)nd (-)nd (-)4.030.676
Enniatin A10.3 (12.5)0.5 (4.5)0.3 (100)nd (-)0.6 (14.3)0.090.176
Enniatin B0.6 (33.3)0.9 (9.1)0.8 (100)2.6 (11.1)3.7 (28.6)1.440.15
Enniatin B11 (29.2)2.2 (9.1)1.5 (100)1.8 (11.1)15.2 (28.6)6.690.202
Enniatin B2nd (-)nd (-)0.9 (50)nd (-)0.6 (28.6)0.330.562
Beauvericin4.9 (83.3)8.5 (100)30.8 (100)19.7 (100)27.2 (100)13.57p < 0.05
nd: mycotoxin not detectable in any samples. √MSE: root mean square error. When not detectable, the limit of detection of specific mycotoxins was used to compute statistical analysis. 1 Concentration represented the average values of only detectable mycotoxins. 2 Regulated mycotoxins by EU legislations.
Table 4. Concentration 1 (µg/kg of DM) and incidence (%, within parenthesis) of Penicillium- and different fungal strains-produced mycotoxins in corn silages belonging to different clusters.
Table 4. Concentration 1 (µg/kg of DM) and incidence (%, within parenthesis) of Penicillium- and different fungal strains-produced mycotoxins in corn silages belonging to different clusters.
ItemsCluster 1Cluster 2Cluster 3Cluster 4Cluster 5√MSEp Value
n = 24n = 22n = 2n = 9n = 7
Penicillium mycotoxins count, n4.64.56.55.46.31.10.001
Skyrinnd (-)nd (-)nd (-)nd (-)6.2 (14.3)--
Asperglaucidend (-)nd (-)nd (-)1.2 (11.1)nd (-)--
Secalonic acid4.1 (29.2)5.3 (50)nd (-)7.3 (33.3)17.9 (57.1)10.780.214
7-Hydroxypestalonic1.2 (4.2)1.4 (22.7)nd (-)1.7 (44.4)2.6 (28.6)0.420.045
Chlorocitreoroseinnd (-)nd (-)1.3 (50)nd (-)nd (-)--
Pestalotin3.4 (45.8)3.9 (27.3)26.5 (100)3 (100)5.2 (85.7)4.29-
Oxaline5.7 (62.5)8 (63.6)3.8 (100)4.8 (55.6)9.9 (85.7)5.230.314
Flavoglaucin10.4 (8.3)2.5 (9.1)nd (-)nd (-)nd (-)9.460.488
Cyclopenin1 (25)nd (-)1 (100)nd (-)nd (-)0.580.972
Fellutanine A109.1 (100)54 (100)577.6 (100)49.9 (100)55.5 (100)79.6-
Mycophenolic acid5.7 (4.2)nd (-)nd (-)nd (-)43.5 (28.6)21.520.388
Mycophenolic acid IV0.5 (4.2)nd (-)nd (-)nd (-)2.3 (28.6)2.090.6
Phenopyrrozin23.3 (100)18.1 (95.5)47.1 (100)10 (100)15.6 (100)24.670.339
Questiomycin A18.3 (79.2)12.7 (86.4)51.7 (100)18.4 (100)33.2 (100)17.670.015
Metabolites from different fungal strains count, n0.60.11.5-0.90.970.103
Ilicicolin A0.8 (4.2)0.3 (4.5)nd (-)nd (-)nd (-)--
Ilicicolin B1.1 (4.2)nd (-)nd (-)nd (-)1.7 (14.3)--
Citreorosein3.4 (12.5)nd (-)3.9 (50)nd (-)4.6 (28.6)1.740.775
Ascochlorin1.3 (20.8)1.4 (4.5)3.9 (50)nd (-)7.7 (14.3)0.880.011
Bassianolide5.2 (4.2)nd (-)nd (-)nd (-)6.3 (14.3)--
Ascofuranone0.5 (12.5)0.4 (4.5)0.8 (50)nd (-)3.2 (14.3)0.240.028
Macrosporin1.5 (12.5)nd (-)nd (-)nd (-)nd (-)--
Unspecific Metabolites count, n0.80.13nd0.70.79<0.05
Iso-rhodoptilometrin0.8 (33.3)1.1 (4.5)1.8 (100)nd (-)0.8 (28.6)0.630.294
Norlichexanthone1 (8.3)nd (-)2 (100)nd (-)5.4 (14.3)0.640.059
nd: mycotoxin not detectable in any samples. √MSE: root mean square error. When not detectable, the limit of detection of specific mycotoxins was used to compute statistical analysis. 1 Concentration represented the average values of only detectable mycotoxins.
Table 5. Chemical, biological and fermentative traits, fermentative quality scores, and silage making procedures characterizing corn silages belonging to different clusters.
Table 5. Chemical, biological and fermentative traits, fermentative quality scores, and silage making procedures characterizing corn silages belonging to different clusters.
ItemsCluster 1Cluster 2Cluster 3Cluster 4Cluster 5√MSEp Value
n = 24n = 22n = 2n = 9n = 7
Chemical parameters (% DM)
DM (% fresh matter)34.734.337.933.531.34.560.342
Ash5.75.75.76.05.80.530.553
CP8.18.27.78.48.90.870.230
EE2.93.03.02.83.00.210.351
NDF37.737.135.539.036.92.400.180
ADF25.224.723.325.824.71.860.389
ADL2.82.93.23.13.00.320.281
NDIP1.11.10.81.11.00.250.507
ADIP0.80.80.60.70.80.160.441
24 h NDFD (% NDF)52.150.947.651.950.73.190.303
Starch31.231.536.128.730.23.440.067
Fermentative parameters (% DM)
pH (dmnl)3.673.763.844.023.830.3310.156
Acetic acid3.422.873.282.543.821.300.223
Propionic acid0.140.130.040.170.410.19<0.05
Butyric acid0.0080.0030.0010.0050.0060.0110.750
Lactic acid3.294.021.913.063.781.6120.246
Lactic to Acetic1.22.20.61.51.11.700.273
Ethanol0.60.50.60.60.30.480.779
Lactic to (Acetic/Ethanol)0.911.780.501.300.991.570.376
1,2 propanediol0.600.490.730.500.170.289<0.05
N-NH3 (% TN)11.3811.717.787.8513.587.7370.569
dmnl = dimensionless; CP = crude protein; EE = ether extract; NDF = neutral detergent fiber; ADF = acid detergent fiber; ADL = lignin; NDIP = neutral detergent insoluble protein; ADIP = acid detergent insoluble protein; 24 h NDFD = NDF disappearance evaluated in situ after 24 h of rumen incubation; N-NH3 = ammonia nitrogen; TN = total nitrogen. √MSE: root mean square error.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Gallo, A.; Ghilardelli, F.; Atzori, A.S.; Zara, S.; Novak, B.; Faas, J.; Fancello, F. Co-Occurrence of Regulated and Emerging Mycotoxins in Corn Silage: Relationships with Fermentation Quality and Bacterial Communities. Toxins 2021, 13, 232. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13030232

AMA Style

Gallo A, Ghilardelli F, Atzori AS, Zara S, Novak B, Faas J, Fancello F. Co-Occurrence of Regulated and Emerging Mycotoxins in Corn Silage: Relationships with Fermentation Quality and Bacterial Communities. Toxins. 2021; 13(3):232. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13030232

Chicago/Turabian Style

Gallo, Antonio, Francesca Ghilardelli, Alberto Stanislao Atzori, Severino Zara, Barbara Novak, Johannes Faas, and Francesco Fancello. 2021. "Co-Occurrence of Regulated and Emerging Mycotoxins in Corn Silage: Relationships with Fermentation Quality and Bacterial Communities" Toxins 13, no. 3: 232. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13030232

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