1. Introduction
Fungi produce many biochemicals through metabolic processes and in response to stimuli [
1]. In fact, many fungi in the plant microbiome engage in “biochemical warfare” among themselves and with other microbes in the community [
2]. They produce numerous defensive compounds to compete against other microbes and defend their space and resources [
3]. Many of the compounds have been exploited by humans for antibiotic activity [
4]; others pose a threat as “mycotoxins”, metabolites with toxic effects on mammals [
5].
Fusarium verticillioides,
F. proliferatum, and
F. graminearum are major maize pathogens and produce several mycotoxins [
6].
Fusarium poae is a common weak pathogen in maize but also produces a number of toxic metabolites [
7]. The former two produce fumonisins, while the latter two produce deoxynivalenol (DON). Both toxins are regulated by the US Food and Drug Administration [
8,
9] and may travel systemically throughout the plant beyond the
Fusarium-infected tissues [
10,
11]. Although neither mycotoxin group is necessary to the
Fusarium species for infection and disease development in maize [
12], their production has been correlated with increased severity of disease [
13,
14] and genes for pathogenicity have been found to be linked to mycotoxin biosynthesis [
15,
16,
17]. These metabolites also have been observed to interact with host resistance genes [
18] and facilitate infection of the host [
12,
13]. However, production of toxins can be specific to the pathogen strain, not just species, [
7,
19,
20], a behavior frequently obfuscating the exact functions of these biochemicals [
21].
Other microbes are known to antagonize or interact with
Fusarium in the ear microbiome and root systems of maize [
22,
23,
24]. Such microbial competition can effect changes in
Fusarium mycotoxin production [
25,
26,
27,
28,
29]. In fact, Keyser et al. [
30] postulated that fumonisins could function as defense compounds when they found that fumonisin extracts reduced growth of several
Fusarium competitors from the maize ear microbiome. Bacon et al. [
26,
31] similarly found that
F. verticilliodes produced fusaric acid toxic to the
Fusarium biocontrol agent
Bacillus mojavensis when challenged with the bacterium. Conversely, Yates et al. [
32] found that
Trichoderma viride could reduce fumonisin production by
F. verticillioides. Similarly, Chatterjee et al. [
33] found that
Clonostachys rosea inhibited fumonisin synthesis by
F. verticillioides but did not degrade the toxin. Hebbar et al. [
34] identified a number of bacteria, including several
Pseudomonas species, in maize plants producing antifungal metabolites inhibitory to
Fusarium verticillioides. In other cases, effects were more complex [
32].
Fusarium species themselves interact chemically within the complex. More than one
Fusarium species may colonize a host [
35] and some, like
Fusarium verticillioides, may colonize the host asymptomatically [
36,
37]—all producing a variety of toxins. For instance,
Fusarium graminearum and
F. verticillioides are known to interact significantly both in laboratory antagonism assays and in the field in terms of both growth and toxin production [
35,
38,
39]. However, the direction of the interaction (antagonistic or synergistic) is highly dependent on environmental parameters [
21,
38,
39].
The southern region of Idaho, USA produces 70% of the world’s sweet corn (
Zea mays convar.
saccharata var.
rugosa) seed [
40]. Here in the semi-arid, continental climate and alkaline soils [
41,
42],
Fusarium verticillioides and
F. proliferatum co-occur in the
Fusarium pathogen complex in sweet corn seed production [
43,
44]. Seed becomes infected via wind-driven and rain-splashed spores colonizing ears and re-infects successive crops systemically from infected seed surfaces and soil [
45,
46]. In addition to multiple
Fusarium species, numerous other microbes colonize the seeds through the silks, some of which are antagonistic to the
Fusarium species
Most frequently, attempts to understand the interactions of other microbes with toxigenic
Fusarium species examine the interaction of a single antagonist with a single
Fusarium species [
24,
27,
28,
47]. Yet, the complexity of the ear microbiome in maize [
22,
23,
27] indicates the biochemical behaviors may be much more complex as well. In a series of field trials, we surveyed the
Fusarium pathogens in sweet corn fields in Nampa, ID, USA.
We hypothesized that different Fusarium species present during infection of the ears at silking would differentially alter mycotoxin production in response to antagonism by five fungal species with putative antagonism toward Fusarium species First, we confirmed antagonism of F. verticillioides and F. proliferatum by our putative antagonists and tested Fusarium isolates for fumonisin and DON production. To test our hypothesis, we quantified effects of fungal antagonists on sporulation of F. verticillioides and F. proliferatum and measured fumonisin production of both Fusarium when antagonized. We hypothesized that antagonists applied to ears at silking would induce changes in Fusarium mycotoxin production in situ, and that effects would occur locally in treated ears and systemically in plants grown from treated seed. We tested this hypothesis in 2015 field trials with the same suite of fungal antagonists. We inoculated emerging silks at flowering and seed at planting to compare local and systemic effects on mycotoxin content. To confirm systemic effects, we conducted a second year of field trials in 2016 and planted seeds developed from 2015 silk inoculations trials. Finally, we used quantitative polymerase chain reaction (qPCR) on mature seed from 2015 silk inoculation trials to determine the effects of fungal antagonists on Fusarium populations in the seed.
3. Discussion
When challenged with fungal antagonists, Fusarium species from sweet corn seed fields in S. Idaho, USA, responded with detectable changes in mycotoxin production even while having no effect on the overall populations of the pathogens in seed. Not only were mycotoxin concentrations changed by direct antagonistic interactions in vitro but were altered by inoculations of antagonists locally in the ear. Moreover, inoculation with the novel Penicillium sp. WPT and Morchella snyderi effected systemic changes in accumulation of those fumonisins and DON in the seed.
Fumonisins [
10] and DON [
11] are known to move systemically in plants, concentrating in tissues without colonization of those tissues by their fungal producers. That single microbes applied either to the seed at planting (
Figure 4) or to developing ears in silk (
Figure 5) could alter those concentrations in first- and second-generation seed lots (respectively) in a complex abiotic and biotic environment would indicate that far more than the
Fusarium species themselves are involved in mycotoxin accumulation in cereals. Whether these toxins are produced by
Fusarium species as defense compounds cannot be determine by our study. However, challenging
Fusarium species with antagonists did promote changes in both fumonisin and DON concentrations relative to unchallenged
Fusarium both in vitro and in situ, indicating that other microbes in the ears of sweet corn or maize may strongly influence bioaccumulation of these toxins.
The novel species
Penicillium sp. isolate WPT mostly stimulated fumonisin production by
Fusarium both in vitro and in situ. The strongest effects of inoculation with the antagonist were systemic in the 2015 field trial in which seed was treated prior to planting. Mature seed produced on plants that developed from treated seed averaged five times the fumonisin content of seed from untreated control plants. It is interesting to note that
Penicillium WPT produces griseofulvin (J. David Miller, Carlton U. ON CA, personal communication), an anti-fungal used pharmaceutically and known to suppress fungal disease in crops [
49,
50]. Griseofulvin has been shown to move systemically in plants [
51,
52], including major grass crops such as wheat, even from soil-borne microbes [
53].
In vitro, antagonism by Penicillium WPT stimulated fumonisin production by F. verticillioides but reduced production by F. proliferatum. Penicillium griseolum had almost the opposite effects: it reduced fumonisin output by F. verticillioides while having no effect on F. proliferatum fumonisins. Penicillium griseolum had no effect either locally or systemically in situ, despite producing low concentrations of griseofulvin (J. David Miller, Carlton U. ON CA, personal communication), indicating that microbial biochemical interactions are likely specific to species. Griseofulvin accumulation was not determined either in vitro or in situ, so that no conclusions may be drawn as to its effects against Fusarium species from sweet corn and consequent effects on fumonisins. However, any biochemical antagonism between these Penicillia and Fusaria species would appear to be highly species dependent.
Just as antagonistic interactions were mostly species dependent, responses were also context dependent. Even the effects of
Penicillium WPT antagonism on systemic fumonisin accumulation were weaker in the 2016 trials compared to 2015 trials. Whether the seed was inoculated at planting or the antagonist was introduced into the seed via the silks during ear development led to different systemic effects from
M. snyderi on DON accumulation in seed. Treatment of seeds at planting with
M. snyderi had no effect on DON in the subsequent crop. However, inoculation of developing ears at silking with
M. snyderi reduced DON accumulation to half that of the untreated control in second-generation seed produced from that sown in 2016. By contrast,
Fusarium responded to
M. snyderi inoculation in developing ears by increasing DON concentrations in the seed. Different species within the genus have been known to produce anti-microbial compounds [
54,
55,
56], but little is known about role of context in their efficacy. There is no record of anti-microbial activity for
M. snyderi, but it both reduced conidial development by
F. verticilioides and
F. proliferatum in vitro and triggered reduced DON concentrations in field seed by systemic mechanisms.
Clonostachys rosea was isolated from immature kernels following silk inoculations, but we did not analyze either seed or PCR products to determine the persistence of fungal inoculants in samples tested for mycotoxins and
Fusarium abundance. Although inoculations of
Fusarium antagonists to seed and ears resulted in differing mycotoxin accumulations in subsequent seed, this response does not indicate a persistence of the inoculants within the ear microbiome. Effects could have been the result of more complicated interactions than direct antagonism between inoculants and target
Fusarium species. Moreover, effects of antagonism might outlast persistence of the inoculant either form systemic biochemical accumulations (as might be possible with
Penicillium WPT) or indirectly via changes within the existing microbial community of the developing kernel. Just as bioactive secondary metabolites like griseofulvin and the
Fusarium mycotoxin zearalenone may persist in tissues lacking the physical presence of source microbes [
51,
52,
53,
57], changes in biochemical constituency and microbial communities (e.g., relative abundance) may result from initial inoculations and persist even while the inoculant.
The consequences of the complexity of the plant microbiome on seed mycotoxin concentrations cannot be underestimated.
Fusarium proliferatum and
F. verticillioides responded very differently to the same antagonists in controlled pairings. In situ, interactions with antagonists are no longer pairings but are subject to multiple species interactions. We isolated four
Fusarium species from seed, immature kernels, and vegetative tissues from the field. Our findings were consistent with findings by Mohan and Wilson [
58] with
F. verticillioides being most common [
48] followed by
F. proliferatum.
Fusarium temperatum and
F. poae were isolated, but rare. However, none of our isolates produced DON in vitro, yet DON accumulated in mature seed and concentrations varied with antagonist interactions. It is likely that other
Fusarium species are present in the system that were not successfully isolated. For example,
F. graminearum has been repeated isolated from wheat heads in S. Idaho [
59].
Fusarium species are known to antagonize each other [
35,
38,
39]. Moreover, genetic variation among strains of the same
Fusarium sp. vary in mycotoxin production [
20]. From our own data we saw that interactions among
Fusarium species (
Figure 1) changed fumonisin concentrations. Although only one or two species may produce a particular mycotoxin in situ, other
Fusarium present may influence accumulation and concentrations in seed.
In addition to multiple
Fusarium species,
Penicillium, Aspergillus, Rhizopus, Cladosporium and
Trichoderma species are known to be present in sweet corn seed [
44,
60,
61]. We obtained our isolate of
Pichia membranifaciens as well as several unidentified bacteria and fungi from seed produced in 2014 field trials. In vitro,
P. membranifaciens effectively suppressed both growth and fumonisin production by
F. proliferatum and
F. verticillioides. What effects these natural populations of antagonists might have as “third parties” to inocula can be partially seen in our 2016 field trials (
Figure 5). Although
M. snyderi alone reduced DON output to nearly half that in plants from untreated seed, concentrations in plants in the dual
M. snyderi +
P. membranifaciens inoculation treatment was similar to untreated. Thus, other microbes besides inocula play no small role in
Fusarium output of mycotoxins.
Antagonistic inocula and perhaps naturally occurring antagonists are also likely to influence relative mycotoxin concentrations in seed by selectively antagonistic behavior toward
Fusarium species. In our qPCR analysis, we determine that
C. rosea,
P. griseolum and the combination treatment
M. snyderi +
P. membranifaciens increased the ratio of
F. verticillioides to total
Fusarium (
Figure 7).
Environmental effects, such as temperature and relative humidity, are known to contribute to variation in mycotoxin production [
21,
38,
62,
63]. To follow commercial rotation schedules, field sites changed each year of our study. We saw overall DON and fumonisin concentrations decline in 2016 compared to 2015, in the case of DON by an order of magnitude. This may have been in part due to changes in site location or the result of differences in precipitation and temperature patterns between the two years. Reductions in 2016 may also have resulted from shifts in species relative abundance of microbiota, specifically the Fusarium complex, within the ears.
Interactions in situ were mostly not predicted by in vitro results. Even where strong, antagonistic suppression of both fumonisin production and conidial development were observed in laboratory experiments, antagonists were mostly not suppressive in situ. For example, all antagonists strongly reduced conidial development of
F. verticillioides in vitro; in situ, seed treated with
Penicillium WPT at silking had over 15 times more
F. verticillioides than untreated seed. Although in vitro interactions are frequently used to screen for microbial antagonism they often fail to have predictive value for field or in situ response variables [
64,
65,
66].
If we would truly understand the variables that contribute to
Fusarium mycotoxins in sweet corn seed or maize production, we must begin to study the system as a whole in situ, where effects of the biotic environment as well as the abiotic might be studied and dissected. If antagonistic inocula can have systemic effects on fumonisin and DON concentrations in sweet corn seed in the field, it is likely that we have underestimated the effects of the entire plant microbiome and microbial “biochemical warfare” [
2] on mycotoxins.
4. Materials and Methods
4.1. Isolating the Sweet Corn Fusarium Complex
Fusarium species were isolated from seed samples collected from sweet corn field trials in Nampa ID, USA during 2014 and 2015 field trials. Seed samples were collected when kernels were immature and again at maturity. Seeds were plated onto 4% potato dextrose agar (PDA, Difco Laboratories, Inc., Detroit, MI, USA) and incubated for up to four weeks at 22 °C. Fusarium isolates were aseptically transferred to pure culture on 4% PDA.
4.2. Fusarium Species Identification
Twelve isolates of Fusarium species were selected for sequence identification to represent individual morphotypes from the pure culture collection. Genomic DNA from the 12 cultured Fusarium isolates was amplified using the ITS-1F and LR3 and the EF1 and EF2 primer pairs. PCR products were checked by agarose gel electrophoresis and purified by PEG/NaCl precipitation. These templates were sequenced with BIgDye3.1 (Applied Biosystems®, Foster, CA, USA) using the same primers used for amplification and analyzed on an ABI3130XL Genetic Analyzer (Life Technologies, Inc., Thermo Fisher Scientific, Waltham, MA, USA) with an 80cm capillary array. Forward and reverse read consensus sequences were checked against GenBank and candidate Fusarium species were identified from the EF1/2 sequences at 99–100% sequence identity and cross-validated by micromorphology.
4.3. Selecting Fungal Antagonists
Fungal antagonists of
Fusarium species were selected for laboratory and field research using in vitro antagonism assays. Fungi of interest were obtained from indigenous plant microbiomes of the intermountain west, including
Pseudotsuga menziesii var.
glauca,
Abies grandis, and
Pinus ponderosa [
67]. These fungi were previously observed to antagonize
F. proliferatum isolated from
P. menziesii in culture. A well-known mycoparasitic fungus
Clonostachys rosea, isolated from roots of
P. menziesii var.
glauca was deliberately selected for its mycoparasistic activity, and a fungistatic yeast,
Pichia membranifaciens isolated from sweet corn was also selected. New in vitro assays were designed to screen for antagonistic behavior by inoculating petri plates of 4% PDA with paired antagonists. A center line was drawn across the diameter of the plate.
Fusarium species were inoculated into the center of the space on one side of the plate approximately 1 cm from the line. Directly opposite, one antagonist was inoculated onto the plate 1 cm from the line. Control plates were grown with only the
Fusarium or an antagonist. Observations were taken of the plates for five to seven days at 22 °C. Fungi that appeared to reduce growth and spread of
Fusarium species across the center line of the plate and/or appeared to reduce growth or stimulate pigmented secondary metabolites with physical contact were subjected to a quantitative test in broth antagonism cultures.
4.4. In Vitro Quantitative Tests of Antagonism
Fusarium species were co-inoculated with selected antagonists in paired antagonism broth cultures with a Fusarium culture broth containing a solution of 2% sucrose, 2% glucose, 1.5% peptone, and 1% potato dextrose adjusted to a pH of 6.2. Each 50-mL culture contained either F. verticillioides or F. proliferatum and a single antagonist. Five antagonists observed to antagonize Fusarium species in vitro were tested: Clonostachys rosea, Morchella snyderi, Penicillium sp. nov. WPT, Penicillium griseolum, and P. membranifaciens. Control cultures contained only Fusarium species
To inoculate the broth cultures, two-week-old PDA plate cultures of each fungus (antagonists or Fusarium spp.) were flooded with sterile distilled water (SDW) and scraped with a bent, glass rod to loosen fungal propagules into solution. Propagules of M. snyderi, which forms heavy hyphal filaments and sclerotia on PDA, were loosened using a scalpel blade. One hundred microliters of Fusarium propagule solution were inoculated into each 100 mL test tube containing 50 mL of the broth solution for each antagonism. For each antagonism treatment, 100 µL of antagonist propagule solution were added. Four replicate cultures were made for each treatment with four antagonist-free Fusarium cultures for a control. Fusarium proliferatum and F. verticillioides solutions were inoculated into their respective broth cultures at 107 propagules per ml. Starting cultures contained 2 × 105 propagules per ml of each pathogen. Inocula were not standardized among antagonists since effects among antagonists were not comparable due to highly different mechanisms of antagonism and the antagonists were not being compared. Starting concentrations in propagules per ml of broth were as follows: C. rosea, 7 × 104; M. snyderi, 7 × 101; P. griseolum, 3.8 × 104; Penicillium WPT, 3.4 × 105; P. membranifaciens, 1.2 × 104.
Following inoculation, the broth cultures were incubated for three days at 22 °C on a shaker at a constant 150 rpm. On day three, counts of Fusarium conidia were made in each culture using a Neubauer cell counting chamber. Counts were computed to the concentration of conidia per ml.
4.5. Mycotoxin ELISA—Direct Antagonisms
Sequenced Fusarium isolates from the sweet corn seed isolate collection were assayed for their ability to produce mycotoxins. Three putative fumonisin producers and a single putative deoxynivalenol (DON) producer were assayed. Pure cultures of the four isolates were cultured on 4% PDA for ten days at 22 °C. Ten grams of culture material (mycelium and agar) were removed from the plates and placed in jars. Fumonisins were extracted with 70% HPLC-grade methanol at a ratio of 1:5 (10 g culture material: 50 mL methanol). DON extracts were made with deionized (DI) water at a 1:5 ratio (culture material: DI water). Extracts were then analyzed using ELISA 96-well mycotoxin enzyme kits (Agriquant®, Romer Labs Inc., Newark, DE, USA) according to manufacturer’s instructions. ELISA plates were scanned using a microplate spectrophotometer (Powerwave™, Bio-Tek Instruments, Winooski, VT, USA) to obtain optical density (OD) readings. Data were transformed using a standard curve to obtain concentrations (µg/g).
Effects of antagonisms on mycotoxin production were determined with a series of in vitro assays pairing Fusarium species with antagonists. Five antagonists were paired with two fumonisin producers: C. rosea, M. snyderi, Penicillium sp. nov. WPT, P. griseolum, and P. membranifaciens. The two fumonisin-producing Fusarium species were also paired with the other two sequenced Fusarium species to determine potential mycotoxin interactions among the isolates. Fumonisin producers were inoculated onto 4% PDA plates 1 cm from the center of the plates. On the opposing side of each plate, also 1 cm from the center of the plate, one of the antagonists or co-occurring Fusarium isolates was plated. Eight plates were made for each of the pairings with co-occurring Fusarium isolates or fungal antagonists, respectively. The positive control plates consisted of the fumonisin producers without competition in the plates. Negative controls were generated for each of the five antagonists by plating each antagonist alone. Only four plates were made for each negative control. After three weeks, a 2-cm wide by 9 cm long strip of agar and fungal material running perpendicular to the dual cultures was removed from each plate. Fumonisins were extracted by submerging each strip in 70% methanol at a ratio of 1:5 (w/v, culture material:methanol). Extracts were analyzed using the ELISA technique as previously described.
4.6. Testing Effects of Antagonism in Situ
4.6.1. 2015. Field Trial—Silk Inoculations
The summer 2015 two in situ antagonism studies were sown in Nampa ID, USA. The first antagonism study examined effects of direct antagonisms via inoculation of emerging silks of a susceptible inbred line (L2) on mycotoxin accumulation in and Fusarium colonization of developing seed. The study also functioned as a seed increase to generate antagonist-carrying F1 hybrid seed to be sown the following year (2016) to study potential carry-over and systemic effects from seed transmission of the inoculants to an F2 generation. The female inbred was sown in early June 2015 in a standard F1 hybrid increase design with four female rows bracketed by two male rows. The male line was sown approximately two weeks previous to planting the female line.
In late July (23 July 2015), emerging silks were inoculated with eight different antagonism treatments: C. rosea, M. snyderi, P. griseolum, Penicillium WPT, P. membranifaciens, M. snyderi + P. membranifaciens, a sterile water control, and an untreated control that received no wet applications. Fungal propagules (spores and macerated hyphae) were suspended in sterile distilled water at concentrations of 104 propagules per ml for M. snyderi and 107 propagules per ml for all other inocula. Inocula were not standardized among antagonists since effects among antagonists were not comparable due to highly different mechanisms of antagonism. Antagonists were applied to emerging silks using hand-held atomizers calibrated to 5 mL of suspension per ear. The M. snyderi + P. membranifaciens treatment received 2.5 mL of each antagonist for a total of 5 mL of inocula. This inoculation was repeated for three nights with inoculations beginning at approximately 6:30PM MST so that each treatment received 15 mL of suspension. Temperature and relative humidity at the silking zone were monitored over a 36-h period beginning the first evening. Treatments were arranged in an 8-block randomized block design at the head of the increase field with 10 plants treated per treatment per block. A total of 640 plants were used in the experiment.
Mature ears were harvested September 17. Ears were dried at 34 °C. Subsamples for each ear were randomly selected for quantitative PCR testing and ELISA mycotoxin assays. Later, seed from individual ears was bulked by treatment for a 2016 planting.
4.6.2. 2015. Field Trial—Seed Treatment for Systemic Effects
The second trial to investigate possible systemic antagonism was planted in May 2015. Antagonist-treated seeds from the Fusarium-susceptible inbred L2 were sown in a replicated trial in Nampa ID. Treatments were blocked in replicated ranges of 12 rows. Each treatment was separated by a buffer row, the seeds of which were treated with standard pesticide treatments. Three antagonist treatments were tested: M. snyderi, Penicillium WTP, and a control. The antagonists were propagated on food-grade yellow cornmeal. The purpose of the cornmeal was to serve as both a culture substrate and an inoculum carrier. Inoculated cornmeal was crushed, mixed, and moistened with deionized (DI) water. The inocula were applied to seeds prior to sowing with the seeds and cornmeal carrier mixed thoroughly to coat each seed. Plain moistened cornmeal was used for the control treatment.
Starting September 25, ears were harvested on a per row basis with a 15-foot row length as a replicate. Ears were dried at 34 °C, and subsamples of mature seeds from each row were taken for mycotoxin ELISA assays.
4.6.3. 2016. Field Trial—Carry-Over Effects
The F1 hybrid seeds generated from the 2015 silk inoculation trial were sown for a final field trial in May 2016. Following 2015 analyses, remaining seed was bulked by treatment. Random samples of the bulked seed were taken for sowing in the 2016 field trail. Six-row seed ranges were sown with a single treatment in a randomized block design with ranges for each treatment represented in each of four blocks. There were nine treatments sown with a total of 36 ranges. The treatments included seeds from all eight silk treatments from 2015: untreated control, water control, C. rosea, M. snyderi, P. griseolum, Penicillium WPT, P. membranifaciens, and seeds from the M. snyderi + P. membranifaciens inoculation. The seeds generated from the 2015 trial were hybrid seeds. They received no further treatment in 2016 prior to sowing. The ninth treatment was a subsample of the 2016 untreated control to which was applied the standard fungicide/insecticide seed treatment prior to sowing.
Ears were harvested by the 10-foot row (4 rows per range) in early October. Ears were dried in a commercial drier at a standard 34 °C. Seed subsamples of this F2 generation were taken from each row for ELISA analyses of mycotoxin content.
4.6.4. Mycotoxin ELISA—2015 and 2016 Field Seed
Seeds from the 2015 field trials and the 2016 field trial were analyzed for the presence and content of both fumonisins and DON. For mycotoxin analyses of the 2015 seedborne antagonism trial, samples were bulked by row and three randomly selected rows were analyzed per treatment. Samples from the 2015 silk-inoculation field trial and the 2016 field trial were bulked by row and eight random rows from each antagonism treatment were selected for analysis. Three randomly selected samples were analyzed for both DON and fumonisins from the 2015 seed treatment field trial. Samples were ground with a coffee mill, and the mill was sanitized between samples. Fumonisins were extracted from ground samples with 70% HPLC grade methanol at a ratio of 1:5 (w/v, sample:methanol). DON was extracted with DI water at a ratio of 1:5 (w/v, sample:DI water). Mycotoxin concentrations of extracts were determine using ELISA 96-well mycotoxin enzyme kits (Agriquant®, Romer Labs Inc., Newark, DE, USA) according to manufacturer’s instructions. ELISA plates were scanned using a microplate spectrophotometer (Powerwave™, Bio-Tek Instruments, Winooski, VT, USA) to obtain OD readings. Data were transformed using a standard curve to obtain concentrations (µg/g).
4.7. Molecular Determination of Fusarium Antagonism
Finally, we used the quantitative polymerase chain reaction (qPCR) technique to compare quantities of Fusarium DNA relative to seed yield (plant DNA) from a subsample of mature seeds from the 2015 silk inoculation field trial. Only 2015 seed was analyzed due to resource constraints. Sixteen randomly selected ears were taken from each of the eight treatments, and four ears were then pooled to generate four replicates per antagonist for analysis. Thirty-gram samples of dry seeds (average 240 kernels) from each treatment and control plot were lyophilized O/N and then ground using a coffee/spice grinder (Kitchen Aid model BCG11108) which has a removable stainless-steel bowl to allow convenient cleaning between samples. Sub-samples of 100 mg of the ground seeds were ground further to a powder on a MixerMill 300 with a tungsten carbide bead in a 2ml Collection Tube (Qiagen N.V., Venlo, Limburg, Netherlands).
4.7.1. Genomic DNA Extraction
DNA was extracted from the powdered seed biomass using the CTAB method [
68]. The CTAB DNA preps were then subjected the PEG/NaCl precipitation. These DNAs were characterized by UV optical absorbance using a microplate spectrophotometer (Epoch
™, BioTek Instruments, Inc., Winooski, VT, USA) and by successful endpoint PCR using tEF1 and ITS fungal primers. For greater consistency, we followed the CTAB extraction with further processing using a Qiagen DNeasy Plant kit according to the manufacturer’s instructions (Qiagen N.V., Venlo, Limburg, Netherlands) to remove residual PCR inhibitors present in the seeds. DNA quality was checked by UV optical absorbance and concentrations were measured using a Qubit
® fluorimeter (Invitrogen
™, Thermo Fisher Scientific, Waltham, MA, USA) with a QuantIt
™ broad range kit (Invitrogen
™, Thermo Fisher Scientific, Waltham, MA, USA).
Genomic DNA was also extracted from the sequenced Fusarium strains isolated from plants in the field trials conducted during 2014 and 2015. Fusarium strains were cultured on 4% PDA plates. Fungal biomass was scraped off the agar surface and DNA was extracted using the Qiagen DNeasy Plant 96 kit according to the manufacturer’s instructions (Qiagen N.V., Venlo, Limburg, Netherlands). DNA quality was checked by UV optical absorbance using an Epoch microplate spectrophotometer (BioTek Instruments, Inc., Winooski, VT, USA). Concentrations were measured using a Qubit® fluorimeter and a QuantIt™ broad range assay kit (Invitrogen™, Thermofisher Scientific, Waltham, MA, USA).
4.7.2. Fusarium Genus Primer Design
Two measures of
Fusarium content were quantified in the seeds: total
Fusarium content (genus level DNA) and the most common species,
F. verticillioides, isolated from the Nampa ID trials. Q-PCR for
F. verticillioides was handled as a separate analysis. Nicolaisen et al. [
69] developed
Fusarium species-specific primer pairs from within the tEF1α gene for 11
Fusarium species which are commonly found in cereals and other crops. During our study four species of
Fusarium were repeatedly isolated from the seeds and identified. Three of these had specific qPCR primers developed by Nicolaisen et al. [
69]; the fourth had not.
To identify regions of sequence conservation in translation elongation factor alpha gene among all Fusarium species of interest we aligned all the tEF1α sequences available in GenBank for 13 Fusarium species that are commonly found in crops in temperate regions. The species included in the alignment were F. avenaceum, F. culmorum, F. equiseti, F. graminearum, F. langsethiae, F. poae, F. proliferatum, F. pseudograminearum, F. sporotrichioides, F. subglutinans, F. temperatum, F. tricinctum, and F. verticillioides. The alignment of 504 sequences was done in CLCBio Main Workbench software (Qiagen Bioinformatics, Aarhus, Denmark). We were able to identify a region where an amplicon of between 50 and 100 bases could be produced for about 90% of the aligned sequences using an individual forward primer sequence (FsGns547F ACGCCTGGGTTCTTGACA) and a mixture of 3 similar reverse primer sequences (Rev1 FsGnsallR CGGTGACATAGTAGCGAGGAGT, Rev2 FsGnspoaeR GGTGACATAGTAGCGGGGAGT, and Rev3 FsGnstmprR TGACAGTGACATAGTAGCGAGGAG). The forward primer was used at a final concentration of 0.167 µM in the qPCR reaction and the reverse primers were used at 0.133 µM each. All primers have a predicted melting temperature (Tm) of roughly 62 °C.
4.7.3. Quantitative PCR Assay
The amplification conditions and primer concentrations were tested in endpoint PCRs and analyzed on 5% high-resolution blend agarose gels (Lonza, Basel, Switzerland). Endpoint PCRs were done using MyTaq Red 2x master mix (Bioline, London, UK). Primer dimers which were formed at typical endpoint PCR primer concentrations were avoided by dilution of the primers. Only the expected products were made when
Fusarium DNA was used as template and the selected primer concentrations were used. Some higher molecular weight material was sometimes produced especially when antagonist-treated seed DNA was used as template. We largely avoided this by keeping the extension times to a minimum and using a higher annealing temperature in the first 5 cycles, in a 2-step cycling program. All qPCRs were performed using SYBR Fast 2× Mix reagent (Bioline, London, UK) and reactions were carried out on an MJ Research Chromo4 (BioRad, Des Plaines, IL, USA) instrument. Once the cycling conditions were optimized for qPCR clean melting curves were obtained and standard curves were linear with from 10 ng down to 1 pg of
Fusarium DNA from the four
Fusarium species isolated during this study. The identities of the PCR products were checked by Sanger sequencing [
70] and found to be as expected.
The thermal cycling program for the total Fusarium genus qPCR assay was as follows:
1. 95 °C | 3 min |
2. 95 °C | 5 s |
3. 68 °C | 8 s |
4. Go to 2 | 2× |
5. 95 °C | 5 s |
6. 67 °C | 8 s |
7. Go to 5 | 1× |
8. 95 °C | 5 s |
9. 66 °C | 18 s |
10. Plate read | |
11. Go to 8 | 34× |
12. Melting curve | 70 °C to 95 °C in 0.5 °C increments, 5 s per step. |
The cycling program for the F. verticillioides primers was the same but with a 2 °C lower annealing/extension temperature.
The quantity of plant DNA in each seed DNA sample was measured by qPCR using the hor1F and hor2R plant-specific primers as described in Nicolaisen et al. [
69]. Seed DNA from one test plot sample that gave the same DNA concentrations in both the Qubit
® and Epoch™ DNA assays was used to construct the standard curve for the plant DNA assay. As determined by Qubit
® measurement, 25ng of seed DNA was used in each qPCR for
Fusarium assays.
Spiking excess pure fungal DNA into qPCRs containing seed DNA showed that there was some PCR inhibition associated with the seed DNA extracts. A correction factor was calculated by spiking two samples of seed DNA from each of the 8 sets of test plot samples with an excess of F. verticillioides DNA and comparing the calculated quantity to that obtained with the fungal DNA alone. This was done for both the Fusarium genus assay and the F. verticillioides assay.
4.8. Statistical Analyses
Propagule counts for
Fusarium species in laboratory antagonism assays were non-parametric and were analyzed using the Kruskal-Wallis Wilcoxon rank score estimate in the NPAR1WAY procedure using SAS/STAT
® software (SAS
® 9.4, 2017, SAS Institute Inc., Cary, NC, USA) [
71]. Mean comparisons were made using standard error as computed from mean score standard deviations.
All absorbance data from mycotoxin ELISA assays were transformed to concentrations using a logit-log function (Romer Labs, Inc. Newark, DE, USA) to generate a standard curve. Concentration data were then analyzed by analysis of variance and Fisher’s LSD and least squares mean comparisons using the GLM and MIXED procedures in SAS.
Results of the
Fusarium genus qPCRs and
F. verticillioides qPCRs, each measured on triplicate samples, were corrected for inhibition in each set of test plot samples and expressed as pg of total
Fusarium or
F. verticillioides DNA per ng of plant DNA. qPCR data analysis was performed using the web-based PCR Miner software package [
72]. Threshold cycle values were standardized for conversion to picograms of
Fusarium DNA per nanograms of seed DNA. These data were then analyzed using the Kruskal-Wallis Wilcoxon rank score estimate in the NPAR1WAY procedure in SAS. Mean comparisons were made using standard error as computed from mean score standard deviations.