Next Article in Journal
Protein-Bound Uremic Toxins Lowering Effect of Sevelamer in Pre-Dialysis Chronic Kidney Disease Patients with Hyperphosphatemia: A Randomized Controlled Trial
Next Article in Special Issue
Role of Sesamia nonagrioides and Ostrinia nubilalis as Vectors of Fusarium spp. and Contribution of Corn Borer-Resistant Bt Maize to Mycotoxin Reduction
Previous Article in Journal
Emergence of New ST301 Shiga Toxin-Producing Escherichia coli Clones Harboring Extra-Intestinal Virulence Traits in Europe
Previous Article in Special Issue
Environmental Conditions Affecting Ochratoxin A during Solar Drying of Grapes: The Case of Tunnel and Open Air-Drying
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distribution, Genetic Diversity and Biocontrol of Aflatoxigenic Aspergillus flavus in Serbian Maize Fields

1
Faculty of Technology, University of Novi Sad, Bulevar Cara Lazara 1, 21000 Novi Sad, Serbia
2
Faculty of Agriculture, University of Novi Sad, Trg Dositeja Obradovića 8, 21000 Novi Sad, Serbia
*
Authors to whom correspondence should be addressed.
Submission received: 28 June 2021 / Revised: 4 August 2021 / Accepted: 4 August 2021 / Published: 27 September 2021
(This article belongs to the Special Issue Occurrence and Integrated Management of Mycotoxins)

Abstract

:
Maize is one of the leading export products in the Republic of Serbia. As a country where economic development depends on agriculture, maize production plays a critical role as a crop of strategic importance. Potential aflatoxin contamination of maize poses a risk to food and feed safety and tremendous economic losses. No aflatoxin contamination of maize samples harvested in 2019 and 2020 in different localities in the Republic of Serbia was detected by the Enzyme-Linked Immunosorbent Assay (ELISA) test and High-Performance Liquid Chromatography (HPLC) method. On the other hand, the Cluster Amplification Patterns (CAP) analyses of the isolated Aspergillus flavus strains from 2019 maize samples confirmed the presence of key biosynthesis genes responsible for aflatoxin production. Artificial inoculation and subsequent HPLC analysis of the inoculated maize samples confirmed the high capacity of the A. flavus strains for aflatoxin production, pointing to a high risk of contamination under favorable conditions. Prevention of aflatoxin contamination is primarily based on A. flavus control, where biocontrol agents play a significant role as sustainable disease management tools. In this study, antagonistic activity screening of the novel strains belonging to the Bacillus genus indicated superior suppression of A. flavus strains by two Bacillus strains isolated from the rhizosphere of Phaseolus vulgaris.

1. Introduction

The share of agricultural production in the gross domestic product (GDP) in the Republic of Serbia accounts for approximately 10%. The country’s economic development heavily depends on the agricultural sector due to its importance in the food industry, where crops are used as raw materials, and the contribution of agricultural products in international trade. In 2020, maize was ranked second on the list of export products in the Republic of Serbia. It supports the fact that maize is the crop of strategic importance for the country [1]. Speaking globally, there are three leading agricultural crops: wheat, rice, and maize [2]. The maize taking the dominant position in the world agricultural system is explained by the possibility of being used in not only food but other industry branches as a multipurpose crop. High nutritional values (carbohydrates content 70–75%) make maize suitable raw material for food and feed production and define this crop as a critical factor for world nutrition and livelihood security [3].
Numerous fungal plant diseases and their effect on the qualitative and quantitative aspects of food production led to significant economic losses measured in billions of US dollars worldwide [4]. On the other hand, a limited number of fungal pathogens can cause severe problems affecting food safety and profitability of the plant production comparable to the negative impact caused by the A. flavus species [5,6]. The estimation of the European Commission (EC) confirms the influence of mycotoxin contamination resulting in annual global crop losses of 5 to 10% [7]. A. flavus is an opportunistic fungal pathogen of crops, predominantly maize, peanuts, and cotton, characterized by the high potential for aflatoxin production [5]. Aflatoxins are potent, highly toxic secondary metabolites that can compromise food and feed security and cause severe health issues [8]. Even exposure to low concentrations of aflatoxin increases the risk of immune suppression, malabsorption of nutrients, infertility, and reduction in life expectancy [9,10]. Aflatoxin class considers four major types, including B1, B2, G1, and G2, and from the food safety point of view, the most relevant is aflatoxin B1 (AFB1) [11]. According to the International Agency for Research on Cancer (IARC) AFB1 is characterized as a carcinogen (Group 1a) [11,12].
The aflatoxigenic potential of A. flavus varies from atoxigenic to highly toxigenic strains [12]. The incidence of a toxigenic A. flavus species has been shown to be associated with geographic origin and substrate characteristics [13]. Previous studies on the distribution of A. flavus in maize fields reported a high incidence of atoxigenic strains (above 70%) [14]. The significance of A. flavus as a plant pathogen lies in aflatoxin contamination, and less in resulting yield losses as a consequence of plant infection. Additionally, contamination levels and manifested infection symptoms are commonly disproportional. It means that even barely noticeable infection signs could be followed by high-level aflatoxin contamination [11]. Regarding the distribution in the ecological niches, A. flavus species are the most common habitants of the tropical environment, with a relatively high temperature range of 28 °C to 37 °C and high relative humidity of about 95% [9]. Under the common climate conditions, the occurrence of the A. flavus on the territory of the Republic of Serbia is not typical [11]. However, global climate changes triggered the more frequent occurrence of the A. flavus species in the regions generally characterized by the low risk of contamination [15]. In the past decade, several aflatoxin outbreaks worldwide raised awareness about the importance of defining strategies effective in preventing the development of aflatoxigenic strains [16]. As a consequence of the extremely drought 2012 producing season, the Republic of Serbia faced an aflatoxin outbreak in maize fields resulting in significant economic losses [11].
The scientific community recognized the application of biocontrol agents as a promising sustainable solution to address the emerging issue [17,18]. Bacteria belonging to the Bacillus genus are conceded as a plant beneficial species in agricultural practice due to high antagonistic activity against phytopathogenic organisms, the ability to promote plant growth, and improve the soil quality [19]. The primary role as a biological tool for plant disease management is defined by the outstanding ability of the Bacillus species to out-compete the target pathogens by the synergistic activity of antimicrobial compounds, competition for nutrient sources and space, and triggering induction of plant’s defense response [20]. Besides antibiotics, Bacillus strains are well known for the production of many other metabolites of interest, including biosurfactants and enzymes [21]. Additional criteria making Bacillus species ideal candidates for the application in biocontrol are spore-forming ability, rapid replication, and resistance to adverse environmental conditions [22].
Numerous studies have proven the high potential of Bacillus species in the suppression of fungal pathogens. Besides promising scientific results and high potential, the full commercialization of microbial biopesticides is still in the preliminary phase, with a limited number of products available on the market. The existing disbalance between the incontestable potential and the current market scenario requires constant research and efforts to isolate novel strains with high antagonistic activity. Estimating their antimicrobial potential and suitability to be used as a central point for designing viable bioprocess solutions is a necessary precondition for boosting the commercialization and broader use of microbial biopesticides [23]. Development of the biotechnological solution for the production of biocontrol agents starts with selecting the strain expressing the highest potential in suppression of the target pathogen. The beginning of the screening procedure typically includes a significant number of potentials producing strains. The rhizosphere soil is a well-known source rich in antagonistic strains. Bacterial strains originating from the rhizosphere naturally coexist in the dynamic environment and constantly interfere with numerous microbial community members [24]. The great advantage of Bacillus strains as biocontrol agents is that they are core soil inhabitants, well adapted to the environmental conditions where they should be lately applied in the form of biopesticide [21].
In the Republic of Serbia, an agricultural practice still relies on the usage of chemical pesticides. The biopesticides market in the country is still underdeveloped and accounts for only 1.3% of the overall market of plant protection products. Data on the biopesticides import in the period from 2015 to 2020 indicates that the largest amount of imported biopesticides refers to bioinsecticides, followed by biofungicides, bioacaricides, and biobactericides [25,26]. Currently, there is a limited number of microbial biopesticides registered and available on the Serbian market. To establish the basis for eco-friendly future in plant protection, co-operation between the scientific community, government, and private sector is needed [23].
Since identification and quantification of aflatoxins in food sources are significant steps in food safety management [8], the principal aim of this study was to evaluate the aflatoxigenic potential of A. flavus species isolated from the maize originated from the different localities of the Republic of Serbia. Aflatoxin contamination of the maize samples harvested in 2019 and 2020 was evaluated by the HPLC and the ELISA method. The genetic potential of isolated A. flavus species in terms of aflatoxin production capacity was examined by CAP analysis. The second aim was to find the effective biological response to the phytopathogenic fungi development by isolating and evaluating the potential of Bacillus spp. strains to be used as biocontrol agents, against toxigenic A. flavus strains. The novel strains were isolated from the rhizosphere soil of different vegetable plants, originated from localities in the Autonomous Province of Vojvodina, Republic of Serbia. The antagonistic effect of Bacillus spp. against toxigenic and atoxigenic A. flavus strains was determined by the well diffusion assay.

2. Results

2.1. Determination of Aflatoxin B1 by the ELISA Test and Total Aflatoxins and Aflatoxin B1 by the HPLC Method

The High-Performance Liquid Chromatography (HPLC) method was applied to determine the aflatoxin B1 (AFB1) and total aflatoxins content in collected samples of maize originated from 10 selected localities harvested in 2019 (Štitar, Valjevo, Pančevo, Sabanta, Subotica, Nadalj, Loznica, Bečej, Sombor, Rogojevac) (Table 1) and 2020 (Rumenka, Oparić, Kuzmin, Lepojević, Martinci, Krušedol, Valjevo, Beška, Bečej, Sombor) (Table 2). Additionally, the aflatoxin B1 content in the maize samples was determined by the Enzyme-Linked Immunosorbent Assay (ELISA) test (Table 3 and Table 4). The obtained results of the ELISA method pointed out that only one sample (LO, 2019) from the Loznica locality has tested positive for the presence of aflatoxin B1 (0.0046 mg/kg). The other samples (95% of the total number) were not contaminated by aflatoxin B1. The results of the HPLC method confirmed the ELISA testing outcome since only the sample originated from Loznica locality (LO) showed a positive result of aflatoxin contamination (0.002 mg/kg). The aflatoxins (including aflatoxin B1) were detected in 1 out of 20 samples (5%).

2.2. Macro and Micromorphological Characterization of A. flavus

Identification of Aspergillus spp. single-spore strains isolated from the collected maize samples from ten localities (Štitar, Valjevo, Pančevo, Sabanta, Subotica, Nadalj, Loznica, Bečej, Sombor, Rogojevac) in 2019 to the species level considered macro and micromorphological characterization of the strains, after isolation using the selective medium for A. flavus [27]. The results of macroscopic observations of the Aspergillus strains are presented in Figure A1 (Appendix C). Initially, the mycelia of A. flavus strains were white. After three days of incubation, the colour changes were observed when the sporulation started and progressed radially over the colonies. White soft velvety colonies turned into the yellow-green compact powdery mass with a whitish margin by the end of five days of incubation. The colonies were flat at the borders and raised in the middle. The 5-days old colony diameter ranged from 3.5 cm to 4 cm. The strains also produced exudates that were brown or colourless, while the reverse of the colonies was pale.
The results of microscopic observations of the Aspergillus strains are presented in Figure A2 (Appendix D). Micromorphology of the isolated strains indicated the presence of colourless, smooth, or finely roughened thick-walled conidiophores. The conidiophores were unbranched and non-septated. The conidia shape was radial to elliptic, while vesicules were globose to sub-globose. Phialides were loosely packed, radiating in all directions from metulae. Based on the presented morphological characteristics and previous isolation of potential aflatoxin producers using the selective medium, all 10 isolated strains were identified as members of the A. flavus species.

2.3. CAP Analysis of the Genetic Profiles of A. flavus Strains

Ten monosporial strains of A. flavus (further designated as SS—single spore) were selected to determine the aflatoxigenic potential by molecular characterization based on Cluster Amplification Patterns (CAP) analysis. The applied molecular technique considers screening deletions in the aflatoxin biosynthesis gene cluster [11]. A total number of 32 CAP markers spaced approximately every 5 kb along 157 kb of the subtelomere region were amplified in four multiplex PCRs [28]. Figure 1 represents the results of the multiplex PCR analysis for 10 monosporial strains of A. flavus isolated from the maize sampled at different localities in the Republic of Serbia in 2019. The obtained results pointed out that nine out of 10 strains show a genetic potential for aflatoxin synthesis. In contrast, only one strain originated from the locality Rogojevac (RO2BSS) possesses significant deletions in the target region, implying atoxigenic character. Making a comparison between nine aflatoxigenic strains, it could be noticed there is a difference in the genetic profile of genes responsible for aflatoxin synthesis among tested strains. The obtained genetic profile of strain from Loznica corresponds to the registered contamination of the maize samples from which it was previously isolated.

2.4. Assessment of Aflatoxigenic Potential of A. flavus Strains by Artificial Inoculation

Artificial inoculation of the maize seed samples was performed to assess the aflatoxigenic potential of the A. flavus strains characterized as potential aflatoxin producers based on the CAP genetic profiles analysis done in the previous investigation step. The artificial inoculation of the maize seed samples aimed to prove the ability of the selected A. flavus strains to produce aflatoxins under simulated favorable environmental conditions, i.e., to confirm expression of the genes responsible for aflatoxin biosynthesis and activation of a corresponding metabolic pathway in the presence of suitable environmental induction factors. After seven days of incubation, the HPLC method was employed to determine aflatoxin presence and content in the infected samples. The results of the HPLC analysis (Table 5) revealed seven samples that tested positive for the presence of aflatoxins: VA1BSS, LO1ASS, SO1ASS, SA2BSS, SU1ASS, PA2DSS, NA2BSS. On the other hand, when it comes to A. flavus strains BČ1CSS, ŠT2BSS, and RO2BSS the expression of genes responsible for aflatoxin biosynthesis did not occur, and no aflatoxin contamination was registered in maize inoculated with these strains. The strain RO2BSS was previously characterized as atoxigenic due to detected gene deletions and results obtained after the artificial inoculation confirmed the lack of potential for the aflatoxin synthesis.

2.5. Potential Bacillus spp. Antagonistic Strains Isolation

In the present study, 76 potential producing Bacillus spp. presented in Table A1 (Appendix A), were isolated from the rhizosphere soil of different vegetable plants, sampled from localities in the Autonomous Province of Vojvodina, Republic of Serbia. The identification was based on conventional techniques according to Bergey’s manual of determinative bacteriology [29].

2.6. Screening of the Bacillus spp. Antagonistic Activity against Aflatoxigenic A. flavus Strain SA2BSS

The preliminary screening included evaluation of the antagonistic effect of 76 Bacillus spp. strains against one aflatoxigenic A. flavus (SA2BSS) strain, which previously showed the potential to produce the largest amount of aflatoxin B1 among the ten tested strains. The cultivation broth samples of 76 Bacillus strains, obtained after four days of cultivation, were tested in triplicates using the well diffusion method.
The One-way ANOVA results (Table 6) pointed out the significant effect of the producing strain on the obtained inhibition zone diameters, confirming the variations of the antagonistic activity against the A. flavus phytopathogen among the tested Bacillus spp. strains (p ≤ 0.05).
Mean values and standard deviations of the inhibition zone diameters obtained by testing cultivation broth samples of 76 producing strains against aflatoxigenic A. flavus SA2BSS isolate are presented in Table A1 (Appendix B). Duncan’s multiple range test was used to define homogenous groups of producing strains at the same level of statistical significance. The highest inhibitory effect was expressed by strains Mah 1a and Kro 4a, which were at the same level of statistical significance, followed by 23 more strains that showed inhibitory effect against the tested phytopathogen. The remaining 51 strains did not show any antagonistic activity.
Ten strains with the highest inhibitory activity registered against aflatoxigenic A. flavus SA2BSS were selected to investigate broader spectrum antimicrobial activity against a larger number of A. flavus strains isolated from corn samples in 2019 from 10 localities in the Republic of Serbia to select an appropriate antagonist for suppression of aflatoxin producers.

2.7. Selection of Bacillus Antagonist for Suppression of Aflatoxigenic A. flavus Strains

The following screening step included 10 Bacillus spp. with the highest antagonistic potential selected after the preliminary testing of inhibitory activity against an aflatoxigenic strain A. flavus SA2BSS. The strains with the highest suppressive effect were tested against strains of A. flavus isolated from maize samples harvested at 18 different localities in the Republic of Serbia in 2019). This screening step included both toxigenic and atoxigenic A. flavus strains.
The testing of the cultivation broth samples of 10 Bacillus producing strains was performed in the same manner as in the previous investigation step, followed by a similar statistical analysis of the obtained experimental data. The One-way ANOVA results, given in Table 7, again confirmed the significant effect of the producing strain on the obtained inhibition zone diameters, with p-value less than 0.05.
Mean values and standard deviations of the inhibition zone diameters obtained by testing cultivation broth samples of 10 selected Bacillus spp. against A. flavus strains obtained from 18 different localities in the Republic of Serbia during 2019 are presented in Table 8, grouped using Duncan’s multiple range test in homogenous groups of the same statistical significance. Six out of ten strains expressed antimicrobial activity against all tested A. flavus strains. The most intensive suppressive activity was exhibited by the strains designated as Mah 1a (Figure 2) and Mah 1b, which are classified in the group of the same level of statistical significance.

3. Discussion

The sampling of maize harvested in 2019 and 2020 was performed at 18 and 10 different localities in the Republic of Serbia, respectively. The samples were investigated for the presence of aflatoxin B1 using the ELISA test (Table 3 and Table 4), and the content of aflatoxin B1 and total aflatoxins content were determined using the HPLC method (Table 1 and Table 2). The results indicated only one sample with positive result originated from the territory of Loznica, whose total aflatoxin content was below the limit defined by the legislative in the Republic of Serbia (0.002 mg/kg), and the sample was considered as safe from the aspect of food safety. Therefore, based on the presented results, it could be concluded that these two growing seasons resulted in the production of aflatoxin-free or health-safe maize in the selected localities in the Republic of Serbia. In terms of weather conditions, both years were characterized as warm seasons with average annual precipitation rate, and heavy rains during May and June in 2019, and June of 2020 [30,31]. Described weather conditions are defined as convenient for undisturbed maize production and timely harvest [32]. Previous studies indicated the weather conditions influence the incidence and level of aflatoxin contamination of maize grown in the Republic of Serbia. For instance, in the period from 2009–2011 occurrence of aflatoxins in maize samples was not detected [33]. A significantly different scenario happened only a year later. Weather conditions changes, including hot and dry spring and summer in 2012, and drought period that lasted from June to September, resulted in heavy infections of maize by A. flavus and, consequently, significant aflatoxin contamination [34]. The contamination level of the maize samples was in the range from 1.01 to even 86.1 µg/kg [33]. Similar weather conditions but with an absence of prolonged drought period occurred in growing season 2013. The occurrence of aflatoxins in maize this season indicated a lower contamination frequency of aflatoxins (24.7%) in comparison to 2012 (72.2%) [35]. In contrast with the weather conditions in 2012 and 2013, in 2014 an extreme amount of precipitation was recorded. The increased moisture created unfavorable conditions for A. flavus infections and resulted in absence of aflatoxins in maize samples. A year later, aflatoxin contamination was again recorded (36.5%), but in 2016 high precipitation rate limited growth of the aflatoxigenic fungi (5% of contaminated maize samples) [35]. Vegetation season in 2017 was warmer and dryer above average weather conditions, and results of the analyses for the aflatoxin presence in maize samples collected in the Autonomous Province of Vojvodina, northern agricultural part of the Republic of Serbia, indicated a contamination level of 67% [36].
Afterwards, isolation of potential aflatoxin producers was performed using the obtained maize samples from 2019 from 18 locations. The selective medium was used during isolation to target A. flavus strains, which show the greatest potential to produce aflatoxins [27]. Macromorphological (Figure A1Appendix C) and micromorphological (Figure A2Appendix D) characterization was applied to confirm the belonging of the isolated strains to the species A. flavus. Despite originating from different localities, most strains have shown similarities in morphological traits, which correspond to the morphological characteristics specific for the A. flavus species [37,38].
The aflatoxigenic potential of the isolated A. flavus strains was confirmed by the CAP analyses, previously successfully applied to address the genetic potential for aflatoxins production [15]. The results of CAP analyses have suggested the high distribution of the strains with the genetic potential to produce aflatoxins on the territory of the Republic of Serbia. The results pointed out that even 90% of the strains had a genetic basis for the aflatoxins synthesis, while only one strain isolated from maize sample from the locality Rogojevac (RO2B) had significant deletions in the aflatoxin biosynthesis gene cluster. The inability of this strain to produce aflatoxins is lately confirmed by the artificial inoculation of the maize seeds and the HPLC analyses to determine the content of the produced aflatoxins during artificial inoculation. Methods for monitoring indels within gene clusters required for the biosynthesis of aflatoxins and cyclopiazonic acid (CPA) are used for detecting intraspecies variability of A. flavus, but also for the selection of isolates with atoxigenic properties as potential biocontrol agents [28,39,40,41]. Based on deletions and insertions of nucleotides in the sequence of an aflatoxigenic gene, a pattern that implicates stability of toxigenic properties is created. Therefore, cluster amplification pattern (CAP) analysis provides information about the stability of atoxigenic isolates [28], but also for potential and stability in the synthesis of aflatoxins and CPA. The absence of deletions in both aflatoxin and CPA clusters may be a criterion for the selection of toxigenic isolates given that many authors stated additive or even synergistic effects of aflatoxins and CPA [42,43,44]. The aim of this research was to determine the isolates’ capacity for aflatoxin biosynthesis, to select the most stable and potent isolate in aflatoxin production, and to test the efficacy of biocontrol agent based on Bacillus spp. against primarily toxigenic isolates. Therefore, genes for aflatoxins and CPA were observed in the first place. Sugar cluster, adjacent to aflatoxin and CPA clusters was also monitored for deletions, however, according to available literature sources, the sugar gene cluster has no direct influence on aflatoxin biosynthesis or the expression on genes in the aflatoxin cluster. Nevertheless, there are data about the possible indirect relationship between these two clusters [45]. Aflatoxin formation relies upon carbon source in a way that simple sugars (glucose, sucrose, fructose, and maltose) support aflatoxin synthesis, while peptone, sorbose, or lactose does not. Additionally, close proximity between the two gene clusters indicates a linkage between them in the metabolism of carbohydrates leading to the induction of aflatoxin biosynthesis. Further, the nadA gene in the aflatoxin biosynthetic pathway was considered to be a part of the sugar cluster, however, gene profiling studies using microarray proved that this gene belongs to the aflatoxin gene cluster and has a role in AFG1/AFG2 formation [45,46].
Variations in profiles obtained by the CAP analysis indicate different toxigenic profiles that are in relation to the stability in the biosynthesis of aflatoxin in artificially inoculated samples. Isolates that expressed the highest potential for aflatoxin production (VA1BSS, SA2BSS, PA2DSS) and those with lower detected aflatoxin levels (LO1ASS, SO1ASS, SU1ASS, NA2BSS) have similar profiles. Also, a similar CAP profile in the aflatoxin cluster (AC01-AC13) have isolates with no aflatoxin detected in the maize sample (BČ1CSS, ŠT2BSS). This could be explained by the fact that the synthesis of aflatoxins depends on various factors that can modulate the expression of genes responsible for coding enzymes that control the biosynthesis pathway [10]. These may also include environmental factors which may activate different cell signaling pathways that can affect the expression of the genes involved in toxin production. The inability of aflatoxin production is a result of deletions which are common for the genes involved in the early stages of aflatoxin biosynthesis. In contrast, genes responsible for the later stages are usually characterized by the presence of single nucleotide polymorphisms (SNPs) [47]. Generated genetic profiles of the aflatoxin biosynthesis gene cluster indicated intraspecies variability between the aflatoxigenic strains, which could be classified into four groups. Genetic diversity among A. flavus strains isolated from different localities could be explained by the difference in cropping practice employed in a certain field [48]. Additional reasons are gene flow as a result of human activities as well as different competition strategies of A. flavus strains depending on the environmental conditions [48].
Climate conditions are critical factors for the growth and development of A.flavus and subsequent aflatoxin biosynthesis [15]. The 2019 and 2020 seasons in the Republic of Serbia were similar and unfavorable for the development of the A. flavus in terms of weather conditions. The results of this study confirmed the lack of maize contamination, except for one sample (from the locality Loznica from 2019), but with the aflatoxin content below the permissible limit. However, proven genetic potential and confirmed gene expression resulting in the high amount of produced AFB1 and total aflatoxins after artificial inoculation of maize seeds have demonstrated the remarkable capacity of the A. flavus strains present in the fields of the Republic of Serbia to produce aflatoxins. Favorable climate conditions, drought, and heat stress increase the probability of A. flavus development and pose a risk for a high level of aflatoxin contamination in maize as an entry point of a food chain. This kind of scenario combined with unpredictable consequences of climate changes implies the necessity of strict control of A. flavus distribution in the fields. Considering the maize as a crop of strategic importance in the Republic of Serbia, the consequences of the potential damage due to the outbreak of aflatoxin contamination would dramatically influence the country’s economy and food and feed safety, as previously seen in season 2012 [49]. If climate conditions changes initiate more regular aflatoxin contamination in the United States of America, as the largest maize producer, it was estimated that losses to the maize industry could reach from $52.1 million to $1.86 billion annually [50]. Therefore, this study was also focused on the investigation of sustainable biocontrol methods for suppression of A. flavus, as a means of preventive action for aflatoxin contamination emergence. On the other hand, the production of aflatoxins brings the more severe consequence of A. flavus presence in the maize fields, but what also should be taken into account are economic losses due to maize fungal infection. Fungal diseases of cereals can cause a yield reduction in the range from 15–20%, but even more in some extreme cases (60%) [51].
Screening of the bacteria belonging to the Bacillus genus as a promising biocontrol agent revealed intensive suppressive activity exhibited by the Mah 1a and Mah 1b strains against toxigenic and atoxigenic A. flavus strains, isolated from the maize grown in the Republic of Serbia. The preliminary screening included 76 Bacillus strains as potential antagonists against one aflatoxigenic A. flavus (SA2BSS) strain, with the highest potential of aflatoxin B1 production. On the other hand, the main screening experiment included all isolated A. flavus strains to evaluate if there is a difference in the activity of Bacillus strains on A. flavus population. Both antagonistic strains characterized by the highest inhibitory activity (Mah1a and Mah1b) were isolated from the rhizosphere soil of the Phaseolus vulgaris. Additionally, antagonistic and phytopathogenic strains originate from the same region, which contributes to the efficiency and adaptation capability of selected biocontrol agents to the environmental conditions of the potential application site [21]. The rhizosphere is a great source of beneficial bacterial strains, and almost 95% of the soil Gram-positive bacilli are a member of the Bacillus genus [52]. The Bacillus strains isolated from the rhizosphere soil of the Phaseolus vulgaris were previously defined as strains of agricultural interest due to their plant-beneficial and pathogen-suppressing activities [22,53,54]. Production of extracellular enzymes is of great importance since it contributes to biocontrol activity and adaptation to the environmental conditions, giving those strains a more competitive advantage over the other microbial inhabitants of the particular ecosystem [52]. The Bacillus strains are marked out as core members of the microbiome in Phaseolus vulgaris rhizosphere [55]. Moreover, Bacillus spp. Isolated from the Phaseolus vulgaris rhizosphere stood out by their superiority among other strains thanks to their plant growth promotion characteristics and potential for antimicrobial metabolites production [56]. Hence, the isolated Bacillus strains Mah 1a and Mah 1b successfully inhibited the growth and development of A. flavus in vitro. These strains are currently being further investigated as potential biocontrol agents for the suppression of fungal maize diseases and aflatoxin contamination. Previous studies also indicated the use of Bacillus spp. For the suppression of fungal pathogens [57,58,59,60], including Aspergillus species [18,61,62,63]. Further research from the aspect of biocontrol product development should first include identifying the selected antagonists and the precise determination of the mechanism of antifungal activity. Depending on the previously defined mechanism of action, further steps in bioprocess development should be determined to achieve the maximization of microbial biomass or metabolites production. Development of the bioprocess solution should include optimization of the medium composition, bioprocess parameters, and downstream procedure for the production of the microbial biocontrol agent. The initial investigation steps should be performed at the laboratory scale with the perspective of scaling up the developed production technology to pilot or industrial level [21,64,65]. All these phases should be followed by in planta testing under field conditions to obtain a highly efficient biocontrol product for suppression of A. flavus and an eco-friendly tool for preventive action against aflatoxin contamination outbreaks.

4. Materials and Methods

4.1. Isolation of Fungal Strains

All A. flavus strains examined in this study were isolated from maize sampled during the 2019 harvest season. Collected maize samples from 18 localities (Pančevo, Užice, Loznica, Subotica, Valjevo, Sirig, Novi Sad, Bečej, Sombor, Maglić, Karavukovo, Nadalj, Kulpin, Sivac, Sabanta, Štitar, Lepojević, Rogojevac) were ground, placed in sterile paper bags and stored at 4 °C until use. Ground maize aliquots of 5 g each were suspended in 25 mL of sterile water and spread on the selective isolation medium—Clean up (CU) medium, supplemented with 5 mg/mL of Bengal Rose and 1 mg/mL of Dichloran, and amended with the antibiotics (10 mg/mL of Chloramphenicol and 10mg/mL Streptomycin) [66]. Isolations were performed in three replicates per sample and CU plates were incubated at 31 °C for three days. Three-day-old plates were examined, and the total number of A. flavus colonies were recorded. Plates with ten or fewer colonies were selected, and pick-ups were conducted by lightly touching one conidiophore of the discrete A. flavus colony and by single point inoculation of the center of 5–2 medium (5% V-8 juice, 2% agar, pH 5.2, 1000 mL H2O) [27]. Plates were incubated for 7–10 days at 31 °C. Aspergillus section Flavi strains were identified to the species level based on macroscopic and microscopic morphological characteristics [27]. Pure cultures were stored in sterile water (six colonized agar plugs with a diameter of 3mm added to 1.5 mL of sterile water) and deposited at the Microbial culture collection of Laboratory for detection of pathogens, pests, and weeds of Faculty of Agriculture, University of Novi Sad.

4.2. Selection of A. flavus Strains and Single Spore Isolation

For further analysis and accurate identification, phenotypic purity of A. flavus strains was achieved through single-spore culturing. Single spore isolations of selected strains were conducted by seral dilutions of spore suspension of strains cultivated at 31 °C for 7–10 days on 5-2 agar medium—100 µL of the most diluted suspensions (10−5 to 10−8 dilutions) were transferred to 1% Malt agar (HiMedia Laboratories, Mumbai, India) medium and incubated at 31 °C for 24–48 h. Plates with ten or fewer colonies were selected and single colonies were transferred to the center of 5-2 medium at 31 °C for 7–10 days. Single spore strains were stored in sterile water and deposited at Microbial culture collection of Laboratory for detection of pathogens, pests, and weeds of Faculty of Agriculture, University of Novi Sad, Novi Sad, Republic of Serbia.

4.3. DNA Extraction from the A. flavus Strains

Ten A. flavus single spore strains characterized based on the colony and spore morphology were selected for further study. To prepare DNA, A. flavus strains were center-point inoculated onto 5-2 agar medium and incubated 8–10 days at 31 °C. The total genomic DNA was extracted using Cetyltrimethyl Ammonium Bromide (CTAB) method [67]. Purified DNA was used as a template for the PCR (polymerase chain reaction) amplification for Cluster Amplification Patterns (CAP) analysis.

4.4. Monitoring Deletions in the Aflatoxin Biosynthesis Gene Cluster of Selected A. flavus Strains

Cluster amplification pattern markers were amplified in four multiplex PCRs [28]. PCR was performed in 10 μL volumes using 6 ng of genomic DNA, 0.08 µmol of each primer and SuperHot MasterMix 2x (Bioron, Römerberg, Germany) on Surecycler 8800 Thermal Cycler (Agilent Technologies, Santa Clara, CA, USA) under the following conditions: 94 °C for 1 min, followed by 30 cycles at 94 °C for 30 s, 62 °C for 90 s, 72 °C for 90 s and the final extension step of 72 °C for 10 min. Products were visualized on 1.4% agarose gel in 1 × sodium boric acid buffer [68].

4.5. Quantitative Determination of Aflatoxin B1 and Total Aflatoxins

4.5.1. ELISA Analysis of Aflatoxin B1 Presence in Maize Samples

Determination of Aflatoxin B1 content in maize samples collected at different locations in the Republic of Serbia (Table 3 and Table 4) in 2019 and 2020 was performed using AgraQuant® Aflatoxin B1 ELISA Test Kit (Romer Labs GmbH, Tulln an der Donau, Austria). Each ground maize sample (20 g) was extracted with 100 mL of 70% methanol and vigorously shook using Benchtop Shaking Incubator 222DS (Labnet International Inc., Edison, NJ, USA) at 200 rpm for 5 min. The sample was allowed to settle, and the top layer of the extract was filtered through a Whatman #1 filter paper (Whatman, UK). Afterwards, 100 µL of the collected filtrate was diluted using the assay buffer. To perform the analysis, 100 µL of each sample or standard was mixed with 200 µL of the conjugate in individual dilution wells, and then 100 µL from each dilution well was transferred to a respective antibody-coated microwell. After 15 min of incubation at room temperature, the plate was washed with distilled water and 100 µL of substrate solution was added to each well, allowed to incubate for 5 min, and then stop solution (100 µL for each well) was added. The absorbance of each well was read at 450 nm (reference wavelength 630 nm) within 10 min after the addition of stop solution using the Microplate Photometer HiPo MPP-96 (BioSan, Rīga, Latvia).

4.5.2. HPLC Quantitative Analysis of Aflatoxins’ Content in Maize Samples

Chromatographic determination of aflatoxins was carried out on a 1260 series HPLC system (Agilent Technologies, Santa Clara, CA, USA) with a DAD (Diode-Array Detector) and FLD (Fluorescence Detector) detectors (Agilent Technologies, USA) and a Hypersil ODS (150 × 4.6 mm i.d., particle size 5 μm) column (Agilent Technologies, USA). Exactly 12.5 g of samples were extracted using 50 mL of acetonitrile (St. Louis, MO, United States) and water mixture (84:16, v/v). The extracts were then cleaned up on MycosepTM 224 column (Romer Labs. Inc., Union, MO, USA). Thereafter, 3 mL of cleaned-up extract was evaporated just to dryness at 60 °C under gentle steam of nitrogen. The residue was dissolved in 300 µL of the mobile phase. HPLC conditions were determined according to Oliveira et al. (2009) [69]. All analyses were conducted in duplicate.

4.6. Potential for Aflatoxin Production in Maize

Artificial inoculation of healthy, uncontaminated (previously analysed by HPLC and ELISA) and undamaged, sterile maize seeds was carried out according to the method described by Probst and Cotty [70] with a few modifications, using the A. flavus strains isolated from maize samples collected during 2019 (ŠT2BSS, VA1BS, PA2DSS, SA2BSS, SU1ASS, NA2BSS, LO1ASS, BČ1CSS, SO1ASS, RO2BSS). Maize seeds (100 g) were placed in 1 L glass jars with perforated lids. Grain moisture was adjusted to 25% moisture content by adding sterile water. Grains were periodically shaken to achieve uniform moisture distribution. Thereafter, grains were autoclaved for 20 min at 121 °C and 2.1 bar. For the inoculum preparation, conidia from 6-day old cultures of A. flavus strains were suspended in sterile distilled water, adjusted to a concentration of 105 spore/mL, and added to each glass jar. The inoculated seeds were incubated for 5–7 days at 30 °C and analysed for aflatoxin contamination (total aflatoxin, aflatoxin B1) by HPLC (as previously described in Section 4.5.2.). Aflatoxin formation is directly affected by temperature. Optimal aflatoxin production is observed at temperatures near 30 °C (28 °C to 35 °C), hence this temperature was selected for the assay [45]. Uninoculated seeds of maize harvested in 2020 were used as control samples and were previously characterized as non-contaminated by aflatoxins using the HPLC method and ELISA test (Section 2.1).

4.7. Isolation of Antagonists

A total number of 76 potential producing strains (antagonists against A. flavus) was isolated from the rhizosphere soil, sampled from different localities in the Autonomous Province of Vojvodina, Republic of Serbia. The selective isolation of sporogenic strains was performed by resuspending 1 g of rhizosphere soil samples in 9 mL of saline and incubating at 28 °C for 15 min at 150 rpm. After homogenisation, thermal treatment at 100 °C for 7 min was performed. Dilution series (10−1, 10−2, 10−3) were prepared, and 500 µL of each dilution was inoculated on nutrient agar (Himedia Laboratories, Mumbai, India) plates and incubated at 28 °C for 48 h. The next step included obtaining pure cultures of morphologically different strains, which were selected according to their morphological and biochemical traits [29] indicating belonging to the Bacillus genus. The selected colonies were picked by a sterile loop and inoculated on nutrient agar plates, followed by incubation at 28 °C for 48 h. The isolated pure cultures were stored on nutrient agar slant at 4 °C at the culture collection of the Laboratory for Biochemical Engineering, Faculty of Technology Novi Sad, University of Novi Sad.

4.8. Inoculum Preparation and Cultivation of Bacillus spp. Antagonists

The inoculum preparation started by incubation of the potential producing Bacillus spp. strains on nutrient agar for 48 h at 28 °C. The second step included transferring the loopful biomass of antagonists into the liquid media (nutrient broth—Himedia, Laboratories, Mumbai, India) and incubating at 28 °C for 48 h on a rotary shaker with an agitation rate of 170 rpm. Inoculation of cultivation media was performed by adding 10% (v/v) of the prepared inocula (5 mL) to the Erlenmeyer flasks containing 50 mL of nutrient broth (Himedia Laboratories, Mumbai, India). The cultivation was carried out on a rotary shaker for 96 h, with a temperature of 28 °C, with an agitation rate of 170 rpm.

4.9. Antimicrobial Activity Assay

The biomass of phytopathogenic A. flavus strains was suspended in sterile saline to achieve a spores’ concentration of 105 CFU/mL. Sabourad maltose agar media (Himedia Laboratories, Mumbai, India) were melted and tempered (50 ± 1 °C) and, before pouring into the Petri dishes, inoculated using 1 mL of the prepared suspensions. The well diffusion method in triplicate tests was employed to evaluate the antagonistic effect of the cultivation broth samples (100 µL) obtained after 96 h of cultivation of the selected Bacillus sp. antagonists against the tested phytopathogens. The incubation was performed at 30 °C for 96 h, and followed by the inhibition zone diameters measurement.

4.10. Statistical Analysis

The analysis of the obtained experimental results included calculating the average values and standard deviations of the measured inhibition zone diameters using Microsoft® Excel 2010 software (Microsoft Corporation, DC, USA). Statistical data analyses were performed using Statistica 13.5 software (Tibco Software Inc., Carslbad, CA, USA), and the employed methods were ANOVA and post hoc testing using Duncan’s multiple range test. All statistical analyses were performed at the significance level of 0.05.

Author Contributions

Conceptualization, M.G., D.B., J.G. and V.V.; methodology, J.G., M.G., D.B., D.M., V.V. and I.P.; software, J.G., I.P. and T.N.; validation, V.V., M.L., T.N. and I.P.; formal analyses, M.L., V.V. and I.P.; investigation, V.V., M.L. and T.N.; resources, J.G., M.G., D.B. and D.M.; data curation, J.G., V.V., M.G., D.B., M.L. and D.M.; writing—original draft preparation, V.V. and M.L.; writing—reviewing and editing, I.P., M.G. and D.B.; visualization, M.L. and V.V.; supervision, J.G., M.G., D.B. and D.M.; project administration, J.G.; funding acquisition, J.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by The Science Fund of the Republic of Serbia, PROMIS, #6064541, BioSolAfla.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Table A1. Potential antagonists isolated from the rhizosphere soil of different vegetable plants originated from localities in the Autonomous Province of Vojvodina, Republic of Serbia.
Table A1. Potential antagonists isolated from the rhizosphere soil of different vegetable plants originated from localities in the Autonomous Province of Vojvodina, Republic of Serbia.
Bacillus sp. Potential AntagonistRhizosphere Soil Samples Used for IsolationLocality of Strain Origin
Cve 2bBeta vulgaris subsp. vulgarisSelenča
Luk 1bAllium cepaMaglić
Luk 2Allium cepaMaglić
Luk 3Allium cepaMaglić
Luk 4aAllium cepaMaglić
Luk 4bAllium cepaMaglić
Kra 3Cucumis sativusBač
Pas 1bPhaseolus vulgarisObrovac
Pas 2aPhaseolus vulgarisObrovac
Pas 2bPhaseolus vulgarisObrovac
Pas 3aPhaseolus vulgarisObrovac
Kra 2cCucumis sativusBač
Pas 4aPastinaca sativaObrovac
Pas 4bPastinaca sativaObrovac
Gra 1Pisum sativumBačko Novo Selo
Gra 2aPisum sativumBačko Novo Selo
Gra 2bPisum sativumBačko Novo Selo
Gra 3aPisum sativumBačko Novo Selo
Gra 3bPisum sativumBačko Novo Selo
Gra 4aPisum sativumBačko Novo Selo
Gra 4bPisum sativumBačko Novo Selo
Kra 2bCucumis sativusBač
Kra 2aCucumis sativusBač
Mah 2aPhaseolus vulgarisTovariševo
Mah 2bPhaseolus vulgarisTovariševo
Mah 3aPhaseolus vulgarisTovariševo
Mah 3bPhaseolus vulgarisTovariševo
Paš 3Pastinaca sativaGložan
Mah 4aPhaseolus vulgarisTovariševo
Mah 4bPhaseolus vulgarisTovariševo
Paš 1aPastinaca sativaGložan
Šar 2aDaucus carota subsp. sativusBegeč
Pap 4bCapsicum annuumOdžaci
Pap 4aCapsicum annuumOdžaci
Pap 2bCapsicum annuumOdžaci
Pap 1bCapsicum annuumOdžaci
Cve 1Beta vulgaris subsp. vulgarisSelenča
Cve 2aBeta vulgaris subsp. vulgarisSelenča
Kra 4Cucumis sativusBač
Cve 3Beta vulgaris subsp. vulgarisSelenča
Cve 4Beta vulgaris subsp. vulgarisSelenča
Kro 1aSolanum tuberosumGajdobra
Kro 1bSolanum tuberosumGajdobra
Kro 2Solanum tuberosumGajdobra
Pap 1aCapsicum annuumOdžaci
Par 4Solanum lycopersicumDeronje
Kup 3bBrassica oleracea var. capitataFutog
Kro 4bSolanum tuberosumGajdobra
Kup 1Brassica oleracea var. capitataFutog
Kup 2Brassica oleracea var. capitataFutog
Kup 3aBrassica oleracea var. capitataFutog
Par 1Solanum lycopersicumDeronje
Kra 1aCucumis sativusBač
Šar 2bDaucus carota subsp. sativusBegeč
Kra 1bCucumis sativusBač
Kro 3aSolanum tuberosumGajdobra
Par 2Solanum lycopersicumDeronje
Paš 2Pastinaca sativaGložan
Paš 4Pastinaca sativaGložan
Mah 3cPhaseolus vulgarisTovariševo
Luk 1aAllium cepaMaglić
Pas 1aPhaseolus vulgarisObrovac
Kro 3bSolanum tuberosumGajdobra
Pas 3bPhaseolus vulgarisObrovac
Kup 4Brassica oleracea var. capitataFutog
Par 3Solanum lycopersicumDeronje
Šar 4Daucus carota subsp. sativusBegeč
Paš 1bPastinaca sativaGložan
Pap 3Capsicum annuumOdžaci
Šar 3bDaucus carota subsp. sativusBegeč
Šar 3aDaucus carota subsp. sativusBegeč
Pap 2aCapsicum annuumOdžaci
Mah 1bPhaseolus vulgarisTovariševo
Šar 1Daucus carota subsp. sativusBegeč
Mah 1aPhaseolus vulgarisTovariševo
Kro 4aSolanum tuberosumGajdobra

Appendix B

Table A2. Duncan’s multiple range test results—mean values and standard deviations of inhibition zone diameters obtained using cultivation broth samples of 76 Bacillus strains against A. flavus SA2BSS.
Table A2. Duncan’s multiple range test results—mean values and standard deviations of inhibition zone diameters obtained using cultivation broth samples of 76 Bacillus strains against A. flavus SA2BSS.
Bacillus sp. AntagonistInhibition Zone Dimeter (mm)
Cve 2b0.00 ± 0.00 a
Luk 1b0.00 ± 0.00 a
Luk 20.00 ± 0.00 a
Luk 30.00 ± 0.00 a
Luk 4a0.00 ± 0.00 a
Luk 4b0.00 ± 0.00 a
Kra 30.00 ± 0.00 a
Pas 1b0.00 ± 0.00 a
Pas 2a0.00 ± 0.00 a
Pas 2b0.00 ± 0.00 a
Pas 3a0.00 ± 0.00 a
Kra 2c0.00 ± 0.00 a
Pas 4a0.00 ± 0.00 a
Pas 4b0.00 ± 0.00 a
Gra 10.00 ± 0.00 a
Gra 2a0.00 ± 0.00 a
Gra 2b0.00 ± 0.00 a
Gra 3a0.00 ± 0.00 a
Gra 3b0.00 ± 0.00 a
Gra 4a0.00 ± 0.00 a
Gra 4b0.00 ± 0.00 a
Kra 2b0.00 ± 0.00 a
Kra 2a0.00 ± 0.00 a
Mah 2a0.00 ± 0.00 a
Mah 2b0.00 ± 0.00 a
Mah 3a0.00 ± 0.00 a
Mah 3b0.00 ± 0.00 a
Paš 30.00 ± 0.00 a
Mah 4a0.00 ± 0.00 a
Mah 4b0.00 ± 0.00 a
Paš 1a0.00 ± 0.00 a
Šar 2a0.00 ± 0.00 a
Pap 4b0.00 ± 0.00 a
Pap 4a0.00 ± 0.00 a
Pap 2b0.00 ± 0.00 a
Pap 1b0.00 ± 0.00 a
Cve 10.00 ± 0.00 a
Cve 2a0.00 ± 0.00 a
Kra 40.00 ± 0.00 a
Cve 30.00 ± 0.00 a
Cve 40.00 ± 0.00 a
Kro 1a0.00 ± 0.00 a
Kro 1b0.00 ± 0.00 a
Kro 20.00 ± 0.00 a
Pap 1a0.00 ± 0.00 a
Par 40.00 ± 0.00 a
Kup 3b0.00 ± 0.00 a
Kro 4b0.00 ± 0.00 a
Kup 10.00 ± 0.00 a
Kup 20.00 ± 0.00 a
Kup 3a0.00 ± 0.00 a
Par 117.00 ± 0.00 b
Kra 1a17.33 ± 1.52 bc
Šar 2b18.00 ± 1.00 bcd
Kra 1b18.33 ± 0.57 bcde
Kro 3a19.00 ± 0.00 bcde
Par 219.00 ± 1.73 bcde
Paš 219.67 ± 0.57 bcdef
Paš 419.67 ± 0.57 bcdef
Mah 3c20.00 ± 0.00 bcdefg
Luk 1a20.00 ± 1.00 bcdefg
Pas 1a20.17 ± 0.76 cdefg
Kro 3b20.33 ± 0.57 cdefg
Pas 3b20.67 ± 0.76 defg
Kup 420.83 ± 1.15 defg
Par 321.33 ± 2.30 efgh
Šar 422.67 ± 0.57 fghi
Paš 1b23.00 ± 1.00 ghi
Pap 324.00 ± 1.00 hi
Šar 3b25.33 ± 4.50 i
Šar 3a25.33 ± 7.81 i
Pap 2a29.67 ± 2.51 j
Mah 1b30.00 ± 5.00 j
Šar 135.00 ± 4.50 k
Mah 1a40.00 ± 0.00 l
Kro 4a42.33 ± 6.80 l
Superscript letters (a–l) represent different levels of statistical significance. Values marked with the same letter are at the same level of significance.

Appendix C

Figure A1. Macromorhological characteristics of A. flavus strains. (a)—ŠT2BSS, (b)—VA1BSS, (c)—PA2DSS, (d)—SA2BSS, (e)—SU1ASS, (f)—NA2BSS, (g)—LO1ASS, (h)—BČ1CSS, (i)—SO1ASS, (j)—RO2BSS.
Figure A1. Macromorhological characteristics of A. flavus strains. (a)—ŠT2BSS, (b)—VA1BSS, (c)—PA2DSS, (d)—SA2BSS, (e)—SU1ASS, (f)—NA2BSS, (g)—LO1ASS, (h)—BČ1CSS, (i)—SO1ASS, (j)—RO2BSS.
Toxins 13 00687 g0a1aToxins 13 00687 g0a1b

Appendix D

Figure A2. Micromorhological characteristics of A. flavus strains. (a)—ŠT2BSS, (b)—VA1BSS, (c)—PA2DSS, (d)—SA2BSS, (e)—SU1ASS, (f)—NA2BSS, (g)—LO1ASS, (h)—BČ1CSS, (i)—SO1ASS, (j)—RO2BSS.
Figure A2. Micromorhological characteristics of A. flavus strains. (a)—ŠT2BSS, (b)—VA1BSS, (c)—PA2DSS, (d)—SA2BSS, (e)—SU1ASS, (f)—NA2BSS, (g)—LO1ASS, (h)—BČ1CSS, (i)—SO1ASS, (j)—RO2BSS.
Toxins 13 00687 g0a2aToxins 13 00687 g0a2bToxins 13 00687 g0a2c

References and Note

  1. Novaković, T. Analiza bruto dodate vrednosti poljoprivrede u Republici Srbiji (The analysis of gross value added of agricultural production in the Republic of Serbia). Ekon. Ideje Praksa 2019, 32, 39–56. [Google Scholar]
  2. Sethuraman, G.; Zain, N.A.M.; Yusoff, S.; Ng, Y.M.; Baisakh, N.; Cheng, A. Revamping ecosystem services through agroecology—The case of cereals. Agriculture 2021, 11, 204. [Google Scholar] [CrossRef]
  3. Shiferaw, B.; Prasanna, B.M.; Hellin, J.; Bänziger, M. Crops that feed the world 6. Past successes and future challenges to the role played by maize in global food security. Food Secur. 2011, 3, 307–327. [Google Scholar] [CrossRef] [Green Version]
  4. Shuping, D.S.S.; Eloff, J.N. The use of plants to protect plants and food against fungal pathogens: A review. Afr. J. Tradit. Complement Altern. Med. 2017, 14, 120–127. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Klich, M.A. Aspergillus flavus: The major producer of aflatoxin. Mol. Plant Pathol. 2007, 8, 713–722. [Google Scholar] [CrossRef] [PubMed]
  6. Massomo, S.M.S. Aspergillus flavus and aflatoxin contamination in the maize value chain and what needs to be done in Tanzania. Sci. Afr. 2020, 10, e00606. [Google Scholar] [CrossRef]
  7. Focker, M.; van der Fels-Klerx, H.J.; Magan, N.; Edwards, S.G.; Grahovac, M.; Bagi, F.; Budakov, D.; Suman, M.; Schatzmayr, G.; Krska, R.; et al. The impact of management practices to prevent and control mycotoxins in the european food supply chain: MyToolBox project results. World Mycotoxin J. 2021, 14, 139–154. [Google Scholar] [CrossRef]
  8. Mahato, D.K.; Lee, K.E.; Kamle, M.; Devi, S.; Dewangan, K.N.; Kumar, P.; Kang, S.G. Aflatoxins in food and feed: An overview on prevalence, detection and control strategies. Front. Microbiol. 2019, 10, 2266. [Google Scholar] [CrossRef]
  9. Okayo, R.O.; Andika, D.O.; Dida, M.M.; K’otuto, G.O.; Gichimu, B.M. Morphological and molecular characterization of toxigenic Aspergillus flavus from groundnut kernels in Kenya. Int. J. Microbiol. 2020, 2020, 8854718. [Google Scholar] [CrossRef] [PubMed]
  10. Caceres, I.; Al Khoury, A.; Khoury, R.; Lorber, S.; Oswald, I.P.; El Khoury, A.; Atoui, A.; Puel, O.; Bailly, J.D. Aflatoxin biosynthesis and genetic regulation: A review. Toxins 2020, 12, 150. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  11. Savić, Z.; Dudaš, T.; Loc, M.; Grahovac, M.; Budakov, D.; Jajić, I.; Krstović, S.; Barošević, T.; Krska, R.; Sulyok, M.; et al. Biological control of aflatoxin in maize grown in Serbia. Toxins 2020, 12, 162. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Khan, R.; Ghazali, F.M.; Mahyudin, N.A.; Samsudin, N.I.P. Biocontrol of aflatoxins using non-aflatoxigenic Aspergillus flavus: A literature review. J. Fungi 2021, 7, 381. [Google Scholar] [CrossRef]
  13. Rodrigues, P.; Venâncio, A.; Kozakiewicz, Z.; Lima, N. A Polyphasic approach to the identification of aflatoxigenic and non-aflatoxigenic strains of Aspergillus section Flavi isolated from Portuguese almonds. Int. J. Food Microbiol. 2009, 129, 187–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Razzaghi-Abyaneh, M.; Shams-Ghahfarokhi, M.; Allameh, A.; Kazeroon-Shiri, A.; Ranjbar-Bahadori, S.; Mirzahoseini, H.; Rezaee, M.B.A. Survey on distribution of Aspergillus section Flavi in corn field soils in Iran: Population patterns based on aflatoxins, cyclopiazonic acid and sclerotia production. Mycopathologia 2006, 161, 183–192. [Google Scholar] [CrossRef]
  15. Medina, A.; Rodriguez, A.; Magan, N. Effect of climate change on Aspergillus flavus and aflatoxin b1 production. Front. Microbiol. 2014, 5, 348. [Google Scholar] [CrossRef]
  16. Battilani, P.; Toscano, P.; Van Der Fels-Klerx, H.J.; Moretti, A.; Camardo Leggieri, M.; Brera, C.; Rortais, A.; Goumperis, T.; Robinson, T. Aflatoxin B1 Contamination in Maize in Europe Increases Due to Climate Change. Sci. Rep. 2016, 6, 24328. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Udomkun, P.; Wiredu, A.N.; Nagle, M.; Müller, J.; Vanlauwe, B.; Bandyopadhyay, R. Innovative Technologies to Manage Aflatoxins in Foods and Feeds and the Profitability of Application—A Review. Food Control 2017, 76, 127–138. [Google Scholar] [CrossRef] [PubMed]
  18. Lagogianni, C.S.; Tsitsigiannis, D.I. Effective Biopesticides and Biostimulants to Reduce Aflatoxins in Maize Fields. Front. Microbiol. 2019, 10, 2645. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Radhakrishnan, R.; Hashem, A.; Abd Allah, E.F. Bacillus: A Biological Tool for Crop Improvement through Bio-Molecular Changes in Adverse Environments. Front. Physiol. 2017, 8, 667. [Google Scholar] [CrossRef]
  20. Cawoy, H.; Bettiol, W.; Fickers, P.; Ongena, M. Bacillus -Based Biological Control of Plant Diseases. In Pesticides in the Modern World—Pesticides Use and Management; Stoytcheva, M., Ed.; InTech: Rijeka, Croatia, 2011; pp. 274–302. [Google Scholar]
  21. Pajčin, I.; Vlajkov, V.; Frohme, M.; Grebinyk, S.; Grahovac, M.; Mojićević, M.; Grahovac, J. Pepper bacterial spot control by Bacillus velezensis: Bioprocess solution. Microorganisms 2020, 8, 1463. [Google Scholar] [CrossRef]
  22. Shafi, J.; Tian, H.; Ji, M. Bacillus species as versatile weapons for plant pathogens: A Review. Biotechnol. Biotechnol. Equip. 2017, 31, 446–459. [Google Scholar] [CrossRef] [Green Version]
  23. Damalas, C.A.; Koutroubas, S.D. Current status and recent developments in biopesticide use. Agriculture 2018, 8, 13. [Google Scholar] [CrossRef] [Green Version]
  24. Mota, M.S.; Gomes, C.B.; Souza Júnior, I.T.; Moura, A.B. Bacterial selection for biological control of plant disease: Criterion determination and validation. Braz. J. Microbiol. 2017, 48, 62–70. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Plant Protection Directorate, Ministry of Agriculture, Forestry and Water Management of the Republic of Serbia, Belgrade, Republic of Serbia (2020). Data on imports of plant protetction products and active substances for the period 2015–2020. obtained on request from the author.
  26. Ministry of Agriculture, Water Management and Forestry. List of Products for Plant Protection in Organic Agricaltural Production. Available online: https://novi.uzb.minpolj.gov.rs/wp-content/uploads/2021/05/Lista_sredstava_za_zas_bilja_za_org_proizvodnju_na_dan_14 (accessed on 15 June 2021).
  27. Cotty, P.J. Virulence and cultural characteristics of two Aspergillus flavus strains pathogenic on cotton. Phytopathology 1989, 79, 808–814. [Google Scholar] [CrossRef] [Green Version]
  28. Callicott, K.A.; Cotty, P.J. Method for monitoring deletions in the aflatoxin biosynthesis gene cluster of Aspergillus flavus with multiplex PCR. Lett. Appl. Microbiol. 2015, 60, 60–65. [Google Scholar] [CrossRef] [PubMed]
  29. Vos, P.; Garrity, G.; Jones, D.; Krieg, N.R.; Ludwig, W.; Rainey, F.A.; Schleifer, K.-H.; Whitman, W. Bergey’s Manual of Systematic Bacteriology: The Firmicutes. In Bergey’s Manual of Systematic Bacteriology, 2nd ed.; Springer: New York, NY, USA, 2009; Volume 3, ISBN 978-0-387-95041-9. [Google Scholar]
  30. RHMZ: Annual Bulletin for Serbia 2019. Available online: http://www.hidmet.gov.rs/data/klimatologija/latin/2019.pdf (accessed on 25 June 2021).
  31. RHMZ: Annual Bulletin for Serbia 2020. Available online: http://www.hidmet.gov.rs/data/klimatologija/latin/2020.pdf (accessed on 25 June 2021).
  32. RHMZ: Annual Agrometeorological Analysis Archives, 2019–2020. Hydrometeorological Service of Serbia. Belgrade (Serbia). Available online: http://www.hidmet.gov.rs/data/agro/godina.pdf (accessed on 25 June 2021).
  33. Kos, J.; Mastilović, J.; Hajnal, E.J.; Šarić, B. Natural occurrence of aflatoxins in maize harvested in Serbia during 2009–2012. Food Control 2013, 34, 31–34. [Google Scholar] [CrossRef]
  34. Papić Brankov, T.; Jovanović, M.; Grujić, B. Aflatoxin standards and maize trade. Ekon. Poljopr. 2013, 60, 595–607. [Google Scholar]
  35. Kos, J.; Janić Hajnal, E.; Šarić, B.; Jovanov, P.; Mandić, A.; Đuragić, O.; Kokić, B. Aflatoxins in maize harvested in the Republic of Serbia over the period 2012–2016. Food Addit. Contam. Part B Surveill. 2018, 11, 246–255. [Google Scholar] [CrossRef] [PubMed]
  36. Glamočić, D.; Jajić, I.; Polovinski-Horvatović, M.; Krstović, S.; Guljaš, D. Aflatoxin in maize silage collected from AP Vojvodina, Serbia. Sci. Pap. Anim. Sci. Biotechnol. 2018, 51, 28–32. [Google Scholar]
  37. Mangal, M.; Bansal, S.; Sharma, M. Macro and micromorphological characterization of different Aspergillus isolates. Legum. Res. 2014, 37, 372–378. [Google Scholar] [CrossRef]
  38. Khan, R.; Ghazali, F.M.; Mahyudin, N.A.; Samsudin, N.I.P. Morphological characterization and determination of aflatoxigenic and non-aflatoxigenic Aspergillus flavus isolated from sweet corn kernels and soil in Malaysia. Agriculture 2020, 10, 450. [Google Scholar] [CrossRef]
  39. Cotty, P.J. Influence of field application of an atoxigenic strain of Aspergillus flavus on the populations of A. flavus infecting cotton bolls and on the aflatoxin content of cottonseed. Phytopathology 1994, 84, 1270–1277. [Google Scholar] [CrossRef]
  40. Grubisha, L.C.; Cotty, P.J. Genetic analysis of the Aspergillus flavus vegetative compatibility group to which a biological control agent that limits aflatoxin contamination in US crops belongs. Appl. Environ. Microbiol. 2015, 81, 5889–5899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Dorner, J.W. Efficacy of a biopesticide for control of aflatoxins in corn. J. Food Prot. 2010, 37, 495–499. [Google Scholar] [CrossRef] [PubMed]
  42. Smith, E.E.; Kubena, L.F.; Braithwaite, C.E.; Harvey, R.B.; Phillips, T.D.; Reine, A.H. Toxicological evaluation of aflatoxin and cyclopiazonic acid in broiler chickens. Poult. Sci. 1992, 71, 1136–1144. [Google Scholar] [CrossRef] [PubMed]
  43. Kumar, R.; Balachandran, C. Histopathological changes in broiler chickens fed aflatoxin and cyclopiazonic acid. Vet. Arhiv. 2009, 79, 31–40. [Google Scholar]
  44. Ostry, V.; Toman, J.; Grosse, Y.; Malir, F. Cyclopiazonic acid: 50th anniversary of its discovery. World Mycotoxin J. 2018, 11, 135–148. [Google Scholar] [CrossRef] [Green Version]
  45. Yu, J. Current understanding on aflatoxin biosynthesis and future perspective in reducing aflatoxin contamination. Toxins 2012, 4, 1024–1057. [Google Scholar] [CrossRef] [Green Version]
  46. Price, M.S.; Yu, J.; Nierman, W.C.; Kim, H.S.; Pritchard, B.; Jacobus, C.A.; Bhatnagar, D.; Cleveland, T.E.; Payne, G.A. The aflatoxin pathway regulator AflR induces gene transcription inside and outside of the aflatoxin biosynthetic cluster. FEMS Microbiol. Lett. 2006, 255, 275–279. [Google Scholar] [CrossRef] [Green Version]
  47. Adhikari, B.N.; Bandyopadhyay, R.; Cotty, P.J. Degeneration of aflatoxin gene clusters in Aspergillus flavus from Africa and North America. AMB Express 2016, 6, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Acur, A.; Arias, R.S.; Odongo, S.; Tuhaise, S.; Ssekandi, J.; Adriko, J.; Muhanguzi, D.; Buah, S.; Kiggundu, A. Genetic Diversity of Aflatoxin-Producing Aspergillus Flavus Isolated from Selected Groundnut Growing Agro-Ecological Zones of Uganda. BMC Microbiol. 2020, 20, 252. [Google Scholar] [CrossRef]
  49. Kos, J.; Skrinjar, M.; Mandic, A.; Misan, A.; Bursic, V.; Saric, B.; Janic-Hajnal, E. Presence of Aflatoxins in Cereals from Serbia. Food Feed Res. 2014, 41, 31–38. [Google Scholar] [CrossRef] [Green Version]
  50. Mitchell, N.J.; Bowers, E.; Hurburgh, C.; Wu, F. Potential Economic Losses to the US Corn Industry from Aflatoxin Contamination. Food Addit. Contam. Part A Chem. Anal. Control Expo Risk Assess 2016, 33, 540–550. [Google Scholar] [CrossRef] [PubMed]
  51. Różewicz, M.; Wyzińska, M.; Grabiński, J. The most important fungal diseases of cereals—Problems and possible solutions. Agronomy 2021, 11, 714. [Google Scholar] [CrossRef]
  52. Nabti, E. Growth stimulation of Phaseolus vulgaris L. plantules by strain Bacillus amylolique faciens hla producer of beneficial agricultural enzymes. JOJ Hortic. Arboric. 2018, 2, 19–25. [Google Scholar] [CrossRef]
  53. Pérez-García, A.; Romero, D.; de Vicente, A. Plant protection and growth stimulation by microorganisms: Biotechnological applications of bacilli in agriculture. Curr. Opin. Biotechnol. 2011, 22, 187–193. [Google Scholar] [CrossRef]
  54. Santoyo, G.; del Orozco-Mosqueda, M.C.; Govindappa, M. Mechanisms of biocontrol and plant growth-promoting activity in soil bacterial species of Bacillus and Pseudomonas: A Review. Biocontrol Sci. Technol. 2012, 22, 855–872. [Google Scholar] [CrossRef]
  55. Stopnisek, N.; Shade, A. persistent microbiome members in the common bean rhizosphere: An integrated analysis of space, time, and plant genotype. ISME J. 2021. [Google Scholar] [CrossRef] [PubMed]
  56. Sendi, Y.; Pfeiffer, T.; Koch, E.; Mhadhbi, H.; Mrabet, M. Potential of common bean (Phaseolus vulgaris L.) root microbiome in the biocontrol of root rot disease and traits of performance. J. Plant Dis. Prot. 2020, 127, 453–462. [Google Scholar] [CrossRef]
  57. Aliye, N.; Fininsa, C.; Hiskias, Y. Evaluation of rhizosphere bacterial antagonists for their potential to bioprotect potato (Solanum tuberosum) against bacterial wilt (Ralstonia solanacearum). Biol. Control 2008, 47, 282–288. [Google Scholar] [CrossRef]
  58. Kumar, P.; Dubey, R.C.; Maheshwari, D.K. Bacillus strains isolated from rhizosphere showed plant growth promoting and antagonistic activity against phytopathogens. Microbiol. Res. 2012, 167, 493–499. [Google Scholar] [CrossRef]
  59. Chen, K.; Tian, Z.; He, H.; Long, C.; Jiang, F. Bacillus species as potential biocontrol agents against citrus diseases. Biol. Control 2020, 151, 104419. [Google Scholar] [CrossRef]
  60. Rong, S.; Xu, H.; Li, L.; Chen, R.; Gao, X.; Xu, Z. Antifungal activity of endophytic Bacillus safensis b21 and its potential application as a biopesticide to control rice blast. Pestic. Biochem. Physiol. 2020, 162, 69–77. [Google Scholar] [CrossRef] [PubMed]
  61. Bharose, A.; Gajera, H. Antifungal activity and metabolites study of bacillus strain against aflatoxin producing aspergillus. J. Appl. Microbiol. Biochem. 2018, 2, 8. [Google Scholar] [CrossRef]
  62. Rajkumar, K.; Naik, M.K.; Amaresh, Y.S.; Chennappa, G. Bioefficacy of Bacillus subtilis against Aspergillus flavus, the cause of aflatoxin contamination in chilli. Int. J. Chem. Sci. 2018, 6, 2050–2053. [Google Scholar]
  63. Siahmoshteh, F.; Hamidi-Esfahani, Z.; Spadaro, D.; Shams-Ghahfarokhi, M.; Razzaghi-Abyaneh, M. Unraveling the mode of antifungal action of bacillus subtilis and Bacillus amyloliquefaciens as potential biocontrol agents against aflatoxigenic Aspergillus parasiticus. Food Control 2018, 89, 300–307. [Google Scholar] [CrossRef]
  64. Grahovac, J.; Pajčin, I.; Vlajkov, V.; Rončević, Z.; Dodić, J.; Cvetković, D.; Jokić, A. Xanthomonas campestris Biocontrol agent: Selection, medium formulation and bioprocess kinetic analysis. Chem. Ind. Chem. Eng. Q. 2020, 27, 131–142. [Google Scholar] [CrossRef]
  65. Montesinos, E. Development, registration and commercialization of microbial pesticides for plant protection. Int. Microbiol. 2003, 6, 245–252. [Google Scholar] [CrossRef] [PubMed]
  66. Cotty, P.J. Comparison of four media for the isolation of Aspergillus flavus group fungi. Mycopathologia 1994, 125, 157–162. [Google Scholar] [CrossRef]
  67. Padmalatha, K.; Prasad, M.N.V. Optimization of DNA isolation and PCR protocol for RAPD analysis of selected medicinal and aromatic plants of conservation concern from peninsular India. Afr. J. Biotechnol. 2006, 5, 230–234. [Google Scholar] [CrossRef]
  68. Brody, J.R.; Kern, S.E. Sodium boric acid: A Tris-free, cooler conductive medium for DNA electrophoresis. Biotechniques 2004, 36, 214–216. [Google Scholar] [CrossRef] [PubMed]
  69. Oliveira, C.A.F.; Gonçalves, N.B.; Rosim, R.E.; Fernandes, A.M. Determination of aflatoxins in peanut products in the northeast region of São Paulo, Brazil. Int. J. Mol. Sci. 2009, 10, 174–183. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Probst, C.; Cotty, P.J.; Harris, S. Relationships between in vivo and in vitro aflatoxin production: Reliable prediction of fungal ability to contaminate maize with aflatoxins. Fungal Biol. 2012, 116, 503–510. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Multiplex PCR amplicons—CAP analysis of the genetic profiles of A. flavus isolates. S—GeneRuler 1 kb Plus DNA ladder (Thermo Fischer), 1—VA1BSS, 2—LO1ASS, 3—RO2BSS, 4—BČ1CSS, 5—SO1ASS, 6—SA2BSS, 7—SU1ASS, 8—PA2DSS, 9—ŠT2BSS, 10—NA2BSS. Primers used in multiplex PCR: SC01, IC01, AC01, AC02, AC03, AC04, AC05, AC06, AC07, AC08, AC09, AC10, AC11, AC12, AC13, IC02, Iac, CC01, CC02, CC03, CC04, ST01, ST02, ST03, ST04, ST05, ST06, ST07, ST08, ST09, ST10, ST11, ST12 [28].
Figure 1. Multiplex PCR amplicons—CAP analysis of the genetic profiles of A. flavus isolates. S—GeneRuler 1 kb Plus DNA ladder (Thermo Fischer), 1—VA1BSS, 2—LO1ASS, 3—RO2BSS, 4—BČ1CSS, 5—SO1ASS, 6—SA2BSS, 7—SU1ASS, 8—PA2DSS, 9—ŠT2BSS, 10—NA2BSS. Primers used in multiplex PCR: SC01, IC01, AC01, AC02, AC03, AC04, AC05, AC06, AC07, AC08, AC09, AC10, AC11, AC12, AC13, IC02, Iac, CC01, CC02, CC03, CC04, ST01, ST02, ST03, ST04, ST05, ST06, ST07, ST08, ST09, ST10, ST11, ST12 [28].
Toxins 13 00687 g001
Figure 2. Inhibiton zones obtained using cultivation broth sample of Bacillus sp. Mah 1 strain against aflatoxigenic A. flavus SA2BSS.
Figure 2. Inhibiton zones obtained using cultivation broth sample of Bacillus sp. Mah 1 strain against aflatoxigenic A. flavus SA2BSS.
Toxins 13 00687 g002
Table 1. Aflatoxins content determination in maize samples collected from 10 different localities in the Republic of Serbia in 2019 by the HPLC method.
Table 1. Aflatoxins content determination in maize samples collected from 10 different localities in the Republic of Serbia in 2019 by the HPLC method.
LocalityStrain MarkAF B1 (mg/kg)Total (mg/kg)
LoznicaLO0.0020.002
SomborSO<0.001<0.001
SuboticaSU<0.001<0.001
PančevoPA<0.001<0.001
Bečej<0.001<0.001
SabantaSI<0.001<0.001
NadaljNA<0.001<0.001
ValjevoVA<0.001<0.001
RogojevacRO<0.001<0.001
ŠtitarŠT<0.001<0.001
Table 2. Aflatoxins content determination in maize samples collected from 10 different localities in the Republic of Serbia in 2020 by the HPLC method.
Table 2. Aflatoxins content determination in maize samples collected from 10 different localities in the Republic of Serbia in 2020 by the HPLC method.
LocalityStrain MarkAF B1 (mg/kg)Total (mg/kg)
RumenkaRU<0.001<0.001
OparićOP<0.001<0.001
KuzminKU<0.001<0.001
LepojevićLE<0.001<0.001
MartinciMC<0.001<0.001
Krušedol<0.001<0.001
ValjevoVA<0.001<0.001
Beška<0.001<0.001
Bečej<0.001<0.001
SomborSO<0.001<0.001
Table 3. Aflatoxin B1 presence determination in maize samples harvested in 2019 collected from 10 different localities in the Republic of Serbia by the ELISA test.
Table 3. Aflatoxin B1 presence determination in maize samples harvested in 2019 collected from 10 different localities in the Republic of Serbia by the ELISA test.
LocalityStrain MarkAF B1 (mg/kg)
ŠtitarŠT0.0000
ValjevoVA0.0000
PančevoPA0.0000
SabantaSI0.0000
SuboticaSU0.0000
VršacNA0.0000
LoznicaLO0.0046
Bečej0.0000
SomborSO0.0000
RogojevacRO0.0000
Table 4. Aflatoxin B1 presence determination in maize samples harvested in 2020 collected from 10 different localities in the Republic of Serbia by the ELISA test.
Table 4. Aflatoxin B1 presence determination in maize samples harvested in 2020 collected from 10 different localities in the Republic of Serbia by the ELISA test.
LocalityStrain MarkAF B1 (mg/kg)
RumenkaRU0.0000
OparićOP0.0000
KuzminKU0.0000
LepojevićLE0.0000
MartinciMC0.0000
Krušedol0.0000
ValjevoVA0.0000
Beška0.0000
Bečej0.0000
SomborSO0.0000
Table 5. Aflatoxins content determination in artificially inoculated maize seed samples by the HPLC method as assessment of aflatoxigenic potential of 10 A. flavus strains.
Table 5. Aflatoxins content determination in artificially inoculated maize seed samples by the HPLC method as assessment of aflatoxigenic potential of 10 A. flavus strains.
Locality of Strain OriginStrain Used for InoculationAF B1 (mg/kg)Total AF (mg/kg)
ŠtitarŠT2BSS<0.001<0.001
ValjevoVA1BSS989.42217.6
PančevoPA2DSS1281.31891.0
SabantaSA2BSS1354.42147.0
SuboticaSU1ASS445.7838.8
NadaljNA2BSS102.7321.8
LoznicaLO1ASS347.9962.4
BečejBČ1CSS<0.001<0.001
SomborSO1ASS330.4564.6
RogojevacRO2BSS<0.001<0.001
Uninoculated control<0.001<0.001
Table 6. One-way ANOVA of inhibition zone diameters for cultivation broth samples of Bacillus spp. antagonists used for suppression of aflatoxigenic A. flavus SA2BSS.
Table 6. One-way ANOVA of inhibition zone diameters for cultivation broth samples of Bacillus spp. antagonists used for suppression of aflatoxigenic A. flavus SA2BSS.
EffectSSDFMSF-valuep-Value
Intercept13,678.75113,678.755142.2190.00
Antagonist31,264.4175416.86156.7090.00
Error404.331522.66
SS—sum of squares, MS—mean squares, DF—degree of freedom.
Table 7. One-way ANOVA of inhibition zone diameters for cultivation broth samples of 10 selected Bacillus spp. antagonists used for suppression of toxigenic and atoxigenic A. flavus strains.
Table 7. One-way ANOVA of inhibition zone diameters for cultivation broth samples of 10 selected Bacillus spp. antagonists used for suppression of toxigenic and atoxigenic A. flavus strains.
EffectSSDFMSF-Valuep-Value
Intercept174,074.11174,074.11422.8280.00
Antagonist200,611.7922,290.2182.1930.00
Error123,567.21010122.3
SS—sum of squares, MS—mean squares, DF—degree of freedom.
Table 8. Duncan’s multiple range test results—mean values and standard deviations of inhibition zone diameters obtained using cultivation broth samples of 10 selected Bacillus spp. against toxigenic and atoxigenic A. flavus strains.
Table 8. Duncan’s multiple range test results—mean values and standard deviations of inhibition zone diameters obtained using cultivation broth samples of 10 selected Bacillus spp. against toxigenic and atoxigenic A. flavus strains.
AntagonistInhibition Zone Diameter (mm)
Šar 3b0.00 ± 0.00 a
Šar 10.00 ± 0.00 a
Pap 2a0.00 ± 0.00 a
Pap 30.00 ± 0.00 a
Paš 1b7.97 ± 12.75 b
Par 311.27 ± 16.74 c
Šar 3a14.74 ± 16.50 d
Kro 4a22.18 ± 18.49 e
Mah 1b36.96 ± 9.81 f
Mah 1a37.52 ± 8.82 f
Superscript letters (a–f) represent different levels of statistical significance. Values marked with the same letter are at the same level of significance.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vlajkov, V.; Grahovac, M.; Budakov, D.; Loc, M.; Pajčin, I.; Milić, D.; Novaković, T.; Grahovac, J. Distribution, Genetic Diversity and Biocontrol of Aflatoxigenic Aspergillus flavus in Serbian Maize Fields. Toxins 2021, 13, 687. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13100687

AMA Style

Vlajkov V, Grahovac M, Budakov D, Loc M, Pajčin I, Milić D, Novaković T, Grahovac J. Distribution, Genetic Diversity and Biocontrol of Aflatoxigenic Aspergillus flavus in Serbian Maize Fields. Toxins. 2021; 13(10):687. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13100687

Chicago/Turabian Style

Vlajkov, Vanja, Mila Grahovac, Dragana Budakov, Marta Loc, Ivana Pajčin, Dragan Milić, Tihomir Novaković, and Jovana Grahovac. 2021. "Distribution, Genetic Diversity and Biocontrol of Aflatoxigenic Aspergillus flavus in Serbian Maize Fields" Toxins 13, no. 10: 687. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins13100687

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