Next Article in Journal
Development of a QuEChERS-Based UHPLC-MS/MS Method for Simultaneous Determination of Six Alternaria Toxins in Grapes
Next Article in Special Issue
Analysis of Toxigenic Fusarium Species Associated with Wheat Grain from Three Regions of Russia: Volga, Ural, and West Siberia
Previous Article in Journal
Investigation of the Metabolic Profile and Toxigenic Variability of Fungal Species Occurring in Fermented Foods and Beverage from Nigeria and South Africa Using UPLC-MS/MS
Previous Article in Special Issue
Functional Analysis of FgNahG Clarifies the Contribution of Salicylic Acid to Wheat (Triticum aestivum) Resistance against Fusarium Head Blight
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Diverse Components of Resistance to Fusarium verticillioides Infection and Fumonisin Contamination in Four Maize Recombinant Inbred Families

1
School of Integrative Plant Science, Cornell University, Ithaca, NY 14853, USA
2
Department of Crop and Soil Sciences, North Carolina State University, Raleigh, NC 27695 USA
3
Department of Entomology and Plant Pathology, North Carolina State University, Raleigh, NC 27695, USA
4
Plant Science Research Unit, United States Department of Agriculture–Agricultural Research Service, Raleigh, NC 27695, USA
*
Authors to whom correspondence should be addressed.
Submission received: 29 November 2018 / Revised: 8 January 2019 / Accepted: 22 January 2019 / Published: 1 February 2019
(This article belongs to the Special Issue Recent Advances in Fusarium Research)

Abstract

:
The fungus Fusarium verticillioides can infect maize ears, causing Fusarium ear rot (FER) and contaminating the grain with fumonisins (FUM), which are harmful to humans and animals. Breeding for resistance to FER and FUM and post-harvest sorting of grain are two strategies for reducing FUM in the food system. Kernel and cob tissues have been previously associated with differential FER and FUM. Four recombinant inbred line families from the maize nested associated mapping population were grown and inoculated with F. verticillioides across four environments, and we evaluated the kernels for external and internal infection severity as well as FUM contamination. We also employed publicly available phenotypes on innate ear morphology to explore genetic relationships between ear architecture and resistance to FER and FUM. The four families revealed wide variation in external symptomatology at the phenotypic level. Kernel bulk density under inoculation was an accurate indicator of FUM levels. Genotypes with lower kernel density—under both inoculated and uninoculated conditions—and larger cobs were more susceptible to infection and FUM contamination. Quantitative trait locus (QTL) intervals could be classified as putatively resistance-specific and putatively shared for ear and resistance traits. Both types of QTL mapped in this study had substantial overlap with previously reported loci for resistance to FER and FUM. Ear morphology may be a component of resistance to F. verticillioides infection and FUM accumulation.
Key Contribution: To our knowledge, this is the first study to document the colocalization of genetic loci underlying ear morphology and resistance to F. verticillioides infection and fumonisin contamination. We also demonstrate the value of kernel bulk density as a low-cost proxy for fumonisin levels in maize grain.

Graphical Abstract

1. Introduction

In maize (Zea mays L), the fungus Fusarium verticillioides causes Fusarium ear rot (FER) and contaminates grain with fumonisins (FUM), a family of mycotoxins produced by Fusarium spp. [1]. FER and FUM contamination can reduce yields and grain marketability [2]. FUM exposure has also been linked to adverse health outcomes in humans and domesticated animals [3,4,5,6,7,8]. The combination of the role of maize as a staple crop, conducive environmental conditions for F. verticillioides infection (FVI), and limited regulation of mycotoxins in food systems has led to high FUM exposure in the developing world [4,5,6,9,10,11].
Pre- and post-harvest management can reduce FUM levels in food systems. Agronomic practices, such as crop residue removal and soil nutrient management, have been associated with reduced risk of FUM contamination in maize [12,13]. Breeding for resistance to FVI and FUM contamination in maize varieties is a more long-term pre-harvest approach [14]. Post-harvest strategies include visual and optical grain sorting [15,16] and thermochemical treatments [17,18].
Genetic resistance to FER and FUM contamination in maize is moderately to highly heritable and quantitatively controlled, making breeding for durable resistance feasible [19,20,21,22]. Both morphological and biochemical characteristics have been shown to contribute genetic resistance. Biochemical signatures associated with differential FVI and FUM levels include lipid profiles [23] and biosynthesis of oxylipins and plant hormones, such as ethylene, abscisic acid, salicylic acid, and jasmonic acid [24,25]. The phenylpropanoid pathway, which is involved in lignin biosynthesis, has also been implicated in resistance to FVI and FUM contamination, although it is unknown whether phenylpropanoids play a passive role in resistance via pericarp hardening or whether they have an active inhibitory effect on F. verticillioides [24,26,27].
Cob and kernel morphological traits have been linked to differential FVI severity and FUM contamination. FVI and FUM levels have been negatively associated with kernel traits related to density and hardness [26,28,29,30,31,32,33], and kernel bulk density (also known as test weight) under F. verticillioides inoculation has been shown to be genetically correlated with FER and FUM [33,34]. Resistance to FUM contamination could therefore hypothetically be accomplished via indirect selection for increased kernel bulk density [35]. Cob morphology, an important component of yield [36], has also been shown to be genetically correlated with FER and FUM [33]. Variation in external symptomatology of kernels infected with F. verticillioides has been associated with differential FUM accumulation [33,37,38,39], suggesting that qualitative measures of kernel symptom severity are important components of resistance.
Understanding the variation in and the role of the maize host genetics in symptomatology of kernels under F. verticillioides inoculation could potentially reveal diverse modes of pathogenesis and/or host resistance. The maize nested association mapping (NAM) population, composed of 25 recombinant inbred line (RIL) families derived from 25 diverse inbred lines crossed to one recurrent inbred parent (B73), is a powerful tool for dissecting the genetic architecture of quantitatively inherited traits [40]. Here we employed quantitative, qualitative, external, and internal indicators of FVI severity as well as publicly available data on ear architecture phenotypes [41,42] to dissect the genetic mechanisms underlying FVI-specific and ear-mediated resistance to FVI and FUM accumulation in four NAM RIL families.

2. Results

2.1. Symptomatology Varies among Families

Four NAM RIL families (B73 × CML333, B73 × CML52, B73 × CML69, and B73 × NC358) were grown in an augmented incomplete block design across four environments. We inoculated the ears with toothpicks coated in F. verticillioides spores, a method that has previously been shown to reveal cob and kernel resistance mechanisms [33]. We evaluated the inoculated ears for external and internal indicators of disease severity. External phenotypes included FER (proportion of visibly infected kernels) and a symptom typology to account for qualitative variation in FER. Internal measures were grain FUM concentration, grain toxin load per visibly infected area (FUM:FER), and kernel bulk density (BDENinoc). As demonstrated by ANOVA, all four quantitative FVI indicators (FER, FUM, FUM:FER, BDENinoc) differed significantly among families (Table 1). On average, families that had higher BDENinoc tended to have lower FER, FUM and FUM:FER (Table 1). The B73 × CML333 family was unique in that it had high FER scores but had the least severe internal disease severity (highest BDENinoc and low FUM and FUM:FER) (Table 1).
The composition of symptom types significantly differed among families, as determined by a χ2 test comparing the proportions of plots exhibiting asymptomatic, blush, starburst, purple, moldy, or multiple symptom types across three environments among the four RIL families (χ2 = 265, p < 0.0001). Across families, the majority of plots (57% to 75%) exhibited the starburst or purple symptom types (Figure 1). The B73 × CML333 family had the greatest proportion of plots exhibiting the blush type (22%) (Figure 1). The B73 × CML52 family had the largest proportion of asymptomatic (20%) and starburst (44%) plots (Figure 1). The proportion of purple plots was greatest in the B73 × CML69 (35%) and B73 × NC358 (40%) families (Figure 1). Approximately 9% of plots exhibited multiple symptoms in all families (Figure 1).

2.2. Relationships among External Symptomatology, Kernel Bulk Density, and Toxin Load Vary among Families

BDENinoc differed significantly among symptom types in all analyses (ANOVA, p < 0.0001) (Table 2). Pairwise t-tests demonstrated that greater external symptom severity was associated with more greatly reduced BDENinoc in all families (Table 2). Specifically, asymptomatic plots had the highest BDENinoc, followed by blush, starburst, purple, and moldy plots in descending order (Table 2).
FER and FUM:FER differed significantly among symptom types in all analyses, but pairwise differences between symptom types varied by family (Table 2). Similarly, symptom types differed significantly with respect to FUM in the combined-family, B73 × CML333, and B73 × CML69 analyses, but symptom severity was not clearly linked to FUM contamination (Table 2). Unexpectedly, asymptomatic plots had the highest or second-highest FUM contamination across all families (Table 2). This phenomenon may be the result of the symptom type scoring method, which was conducted with the kernels still attached to the cob and thus kernels may have had symptoms present below the visible area.

2.3. Fusarium Ear Rot Severity, Fumonisin Contamination, and Kernel Bulk Density are Correlated

We assessed correlations among indicators of FVI severity and FUM contamination at the levels of plot phenotype and genotype means (after controlling for genetic and field effects). Across all families, FUM and FER were (1) positively phenotypically and/or genotype-means correlated with each other and (2) negatively phenotypically and genotype-means correlated with BDENinoc (Table 3). These correlations are consistent with those previously reported [33,35,43,44]. FUM tended to be more strongly correlated with BDENinoc than with FER across families (Table 3), suggesting that internal infection severity (as reflected by BDENinoc) may be a better indicator of FUM contamination that external symptom severity (FER).

2.4. Innate Ear Morphology is a Component of Resistance to FVI and FUM Contamination

Previous studies have reported associations between FVI severity and FUM contamination with characteristics of the cob and kernel tissues, such as cob morphology and kernel density-related traits [28,32,33,34,45,46,47]. Given the symptomatological variation present in the four RIL families, we sought to test the extent to which innate (uninoculated) ear architecture played a role in resistance to FVI using publicly available data on cob density (CobDen), diameter (CobDiam), length (CobLen), mass (CobMass), and volume (CobVol) as well as uninoculated kernel bulk density (BDENuninoc).
Supporting the hypothesis that maize lines with greater innate kernel bulk density (BDENuninoc) in turn have greater bulk density under F. verticillioides inoculation (BDENinoc), we found that BDENuninoc and BDENinoc were positively correlated in all families (Table S1). BDENuninoc was negatively correlated with FER in the B73 × CML52 and B73 × CML69 families and positively correlated with FUM:FER in the B73 × CML333 and B73 × CML52 families (Table S1). BDENuninoc was less strongly correlated with FER and FUM:FER than was BDENinoc (Table 3 and Table 4). Unlike the correlations found between FUM and BDENinoc, FUM was not significantly (p > 0.05) correlated with BDENuninoc in any family (Table S1).
Although the statistical significance of correlations between cob size traits and BDENinoc or BDENuninoc varied among families, the directionality of these correlations was generally conserved across families (Table S1). Both BDENinoc and BDENuninoc tended to be negatively correlated with CobDiam, CobMass, and CobVol (Table S1). Relationships between cob morphology and FER, FUM, and FUM:FER varied from family to family. In general, lines with larger cobs had greater FER severity, as demonstrated by positive correlations between FER and CobDiam or CobVol across families (Table S1). FER was negatively correlated with CobDen in the B73 × CML52 family and positively correlated with CobMass in the B73 × NC358 family (Table S1). FUM:FER correlations with cob size and/or density in the B73 × CML52 and B73 × NC358 families revealed that lines with larger, less dense cobs had reduced FUM:FER (Table S1). We expected to find similar relationships between FUM and cob morphology, but instead found that FUM was only marginally positively correlated with CobMass in the B73 × CML52 family. This may have been due to the low genetic variation for FUM in this study. Similar to previous studies [33,36], cob size traits (CobDiam, CobLen, CobVol) were significantly (p < 0.05) or marginally significantly (p < 0.1) positively correlated with each other and with CobMass. CobDen was significantly positively correlated with CobMass and significantly or marginally significantly negatively correlated with cob size traits (Table S1).

2.5. Genetic and Environmental Variation on Disease Severity Differ among Families

We assessed genetic and environmental effects on indicators of FVI severity within and among families. Across all families, broad-sense heritability of genotype-means (H) was greater for BDENinoc (0.4–0.79) and FER (0.59–0.75) than for FUM (0.12–0.17), FUM:FER (0.12–0.21), or symptom type traits (0–0.44) (Table 4, Figure 2). The variance explained by genotype, environment and the interaction between genotype and environment (G*E) varied widely by trait and family (Table 4). Of the nine FVI traits, the moldy symptom type was most affected by G*E (68% to 85% variance in all families) (Table 4). In general, G*E explained less than 20% of the variance in all other traits across families (Table 4). Traits related to fumonisin contamination (FUM and FUM:FER) had lower genetic variance (1% to 3%) and higher variance explained by environment (46% to 73%) than the other traits across all families (Table 4). In contrast, BDENinoc and FER had the highest genetic variance (BDENinoc = 23% to 64%; FER = 22% to 42%), and environment explained the least amount of variance in BDENinoc (0.3% to 8%) and the moldy symptom type (1% to 3%) (Table 4). The H for FUM found in our experiment was much lower than those previously reported (0.43–0.86) [43,44]. Resistance to mycotoxigenic fungi and kernel composition are influenced by environmental conditions and management practices [13,19,36,37,48,49,50,51], which may explain the large environmental and G*E effects here.

2.6. Genetic Architecture Differs between Ear Morphological and Resistance Traits

In total, stepwise regression selected 271 quantitative trait loci (QTL) associated with the nine FVI traits and the six ear traits in the single and joint family models from the total set of 7386 markers (Table S2). Sixty-six QTL were selected in the joint family models, and 57, 41, 59, and 48 QTL were mapped in the B73 × CML333, B73 × CML52, B73 × CML69, and B73 × NC358 families, respectively (Table S2). Traits that had greater H also had a greater number of associated QTL and greater variance explained by their respective QTL models (model R2) (Figure 2). Cob morphological traits tended to have a greater number of associated QTL (28–37 QTL each) and QTL model R2 (model R2 = 0.31–0.63) than BDENuninoc and FVI traits (3–15 QTL, R2 = 0.06–0.32), with the exception of BDENinoc (20 QTL, model R2 = 0.28–0.39), FER (26 QTL, R2 = 0.18–0.45), and the blush symptom type (10 QTL, model R2 = 0.29–0.30) (Figure 2).

2.7. Most QTL for Resistance to F. verticillioides are Trait-Specific

Within families, the majority of QTL were specific to one trait (161/271 QTL) (Table S3). A total of 110 QTL overlapped with QTL for at least one other trait within families at 44 distinct regions (“colocalized QTL”) (Table S3). Twenty-one colocalized QTL were associated with multiple ear traits but not associated with any FVI trait (Table S3). Ten colocalized QTL were associated with multiple FVI traits but no ear traits (“FVI-specific colocalized QTL”), with three each from the B73 × CML69 and B73 × CML333 families, two from the B73 × NC358 family, and one each from the B73 × CML52 and joint family analyses (Table 5). Thirteen colocalized QTL were associated with both FVI and ear traits (“FVI-ear colocalized QTL”), with five from the joint analysis, three from the B73 × CML52 family, three from the B73 × CML333 family, and one each from the B73 × CML69 and B73 × NC358 families (Table 5).

2.8. Allele Effects at FVI-Specific Loci Reflect Trait Relationships

The directionality of allele effects on quantitative FVI traits (e.g., FER, FUM) that shared QTL with each other and/or with BDENuninoc generally matched the phenotypic correlation between the corresponding traits. For example, FER was negatively correlated with both BDENinoc and BDENuninoc in the B73 × CML52 and B73 × CML69 families (Table 3 and Table 5), and allele effects on FER were opposite to allele effects on BDENinoc and BDENuninoc at three colocalized QTL in the B73 × CML52, B73 × CML69, and joint family analyses (Table 5). Similarly, FER was positively correlated with FUM and negatively correlated with FUM:FER across all families (Table 3 and Table 5), and the allele effect on FER was in the same direction as the allele effect on FUM at a colocalized QTL in the B73 × CML333 family and in the opposite direction as the allele effect on FUM:FER at colocalized QTL in the B73 × CML333 and B73 × CML52 families (Table 5).
QTL colocalized for multiple symptom type traits generally had opposite allele effects on mild (e.g., blush, starburst) vs. severe (e.g., purple, moldy) symptom types (Table 5), indicating that related mechanisms may underlie the manifestation of similar symptomatologies (blush/starburst or purple/moldy). For example, allele effects on starburst were in the same direction as those on blush and opposite to allele effects on moldy and purple at colocalized QTL (Table 5). Similarly, asymptomatic lines by definition have low FER, and allele effects on FER were opposite to asymptomatic allele effects at their two colocalized QTL (Table 5).

2.9. Characteristics of Loci Underlying Specific and Ear-Mediated Resistance to F. verticillioides

We sought to characterize and compare FVI-specific and FVI-ear QTL. First, we tested whether the two classes of resistance loci differed with respect to marker effect size, controlling for trait, family, and QTL length effects. We found that trait and family were significantly (p < 0.0001) and marginally significantly (p < 0.1) associated with absolute marker effect size, respectively. As the trait scale for FER and CobVol was larger than the other FVI and ear traits, these two traits in turn had larger absolute marker effects. Compared to the other families, the B73 × NC358 family had greater marker effect sizes and the joint family models had smaller marker effect sizes. We then investigated whether FVI-specific and FVI-ear loci tended to be found in pericentromeric regions, which have lower recombination rates than non-centromeric regions of the genome [52]. The FVI-specific (10/62 QTL) and FVI-ear (4/13 QTL) groups did not significantly differ with respect to the proportion of QTL located in pericentromeric bins (χ2 = 1.5, p = 0.2). Several other studies have mapped loci for resistance to F. verticillioides, and we identified which of our FVI-specific and FVI-ear QTL colocalized with previously described FVI loci [19,49,53,54,55,56,57,58,59]. A significantly greater proportion of our FVI-specific QTL (46/62 QTL) overlapped published FVI loci than did our FVI-ear QTL (6/13) (χ2 = 4.0, p = 0.04).

3. Discussion

Here we dissect the diversity of genetic resistance to FVI and FUM accumulation in four maize RIL families. We further establish the importance of BDENinoc as a useful proxy for FUM [33,34], especially given the poor heritability of FUM contamination in maize grain. In addition, we demonstrate that some loci underlying resistance FVI and FUM severity may have trade-offs for agronomically important kernel and cob traits.
At the phenotypic, genetic, and locus levels, BDENinoc was negatively associated with external symptom severity and FUM contamination. These results strengthen the arguments that BDENinoc could be used as a predictor of FUM contamination [28,32,33,34] and that selection for increased BDENinoc could indirectly increase resistance to FUM [35], even under conditions resulting in low heritability for FUM. The relatively low heritability for FUM in our experiment may have been due to the inoculation method (insertion of spore-coated toothpick into the developing ear), which has been shown to be less effective for screening genetic sources of resistance to FUM contamination [33,60]. However, toothpick inoculation, which has been shown to be effective at screening ear-mediated resistance mechanisms [33], allowed us to study the role of the cob in resistance to FVI and FUM. We identified 23 novel QTL associated with resistance to FVI, seven of which may be mediated by kernel and cob architecture.
Although no FUM or FUM:FER QTL colocalized with ear QTL in this study, we did find that lines with larger, less dense cobs had greater FER and lower BDENinoc at the QTL and trait correlation levels. In addition, both FER and BDENinoc were phenotypically and genetically correlated with FUM here and in previous studies [33,34] and we found substantial overlap between our ear QTL and previously described loci for resistance to FUM [19,49,53,54,55,56,57,58,59], suggesting that there may indeed be a genetic link between cob morphology and FUM but that the low heritability for FUM in our experiment may have limited our ability to detect this inferred relationship.
As cob morphology and kernel density are important components of grain yield [36] and quality [61], respectively, increased susceptibility to FVI and FUM contamination may be a by-product of the breeding process. We are currently conducting follow-up experiments with near isogenic lines (NILs) to dissect possible pleiotropy between ear architecture and resistance to FVI and FUM. We also recommend that kernel bulk density under F. verticillioides infection be more thoroughly evaluated as a proxy for FUM contamination in maize grain. The wide variation in internal vs. external symptomatology, inter-trait correlations, and QTL colocalizations in this study implies that distinct modes of resistance exist within tropical germplasm, which can be leveraged for resistance breeding.

4. Materials and Methods

4.1. Field Design and Inoculation

Four RIL families (B73 × CML333, B73 × CML52, B73 × CML69, B73 × NC358) from the maize NAM population [40] and the five inbred parents (B73, CML333, CML52, CML69, NC358) were grown at the Central Crops Research Station in Clayton, NC in four year-environments from 2012 to 2015. Each NAM family is composed of 200 RILs. The four non-B73 parents originate from breeding programs in Mexico and North Carolina—regions that are prone to mycotoxin contamination [5,62]—and have previously demonstrated greater resistance to FER than B73 [20,21]. We used an augmented incomplete block design, in which each family was grown separately with 20-plot blocks. RIL plots were replicated once per year-environment and randomized within family, and the two parental lines of each family were randomized in each 20-plot block. All RILs from the B73 × CML52 and B73 × NC358 families were grown in four year-environments from 2012 to 2015. For the B73 × CML333 and B73 × CML69 families, all RILs from these two families were grown in three year-environments from 2013 to 2015 and 40 randomly selected RILs from each of these two families were grown in 2012. We measured days to silking (DTS) in each plot as described by Buckler et al. [63]. One toothpick coated with spores from local toxigenic F. verticillioides isolates was inserted into the middle of each developing primary ear approximately 10 days after silking, and the toothpick remained in the ear until harvest [21].

4.2. Disease Phenotyping

Primary ears were harvested from each plot at maturity, dried, and then visually evaluated for FER. FER was scored based on the percentage of the kernels presenting symptoms on a 1% to 100% scale with 5% increments [19]. The average FER score of all the ears in each plot was then calculated.
We used a symptom typology to account for the qualitative variation in external FER symptoms, in which each ear was assigned a symptom type: “asymptomatic,” “blush,” “starburst,” “purple,” or “moldy” (photographic examples in Morales et al. [33]). Asymptomatic ears had no or very limited visible external symptoms [33,38,39]. The blush type was characterized by pink discoloration localized on the kernel crowns [33]. Starburst ears had whitish streaks radiating from the kernel silk scar and/or pedicel [33,38,64,65]. The kernels of purple type ears were visibly degraded and had severe purplish discoloration [33]. Moldy ears had severely degraded kernels with matted fungal growth [33,37,38]. Each plot was assigned a main symptom type based on the most frequent symptom type of the ears in the plot. Plots with less than 10% average FER were considered asymptomatic. Plots with equal representation of more than one symptom type were categorized as “multiple.”
After scoring FER and symptom typology, ears were then shelled and bulked per plot. From each plot, a random 250-mL volume of kernels was weighed, and kernel bulk density (BDENinoc) was calculated as the weight of the kernels divided by 250 mL. After weighing, the kernels were returned to the bulked plot.
All of the kernels in each bulked plot were ground into a fine powder with a Waring 7010 two-speed laboratory blender (Waring Commercial, Inc., Torrington, CT, USA). A 10-g subsample from each ground bulk was put in a 25 mL centrifuge tube. To extract fumonisins, 20 mL of 90% methanol was added to each 25 mL tube, resulting in a two-fold dilution factor at this step. The tubes were then shaken with a Lab-Line Environ Orbitol Shaker (Lab-Line Instruments, Inc., Melrose Park, IL, USA) at 150 rpm for approximately five minutes. The samples settled for 15 min, after which 0.5 mL of supernatant from each sample was transferred to a 15 mL centrifuge tube. To dilute the supernatant to a final 40-fold dilution, 9.5 mL of distilled water was added to each 15 mL tube. Fumonisin contamination (FUM) was quantified with FUM-specific enzyme-linked immunosorbent assay (ELISA) kits (Helica Biosystems, Inc., Santa Ana, CA, USA). Absorbance at 450 nm of the ELISA plates was read using a BioTek µQuant™ microplate spectrophotometer (BioTek Instruments, Inc., Winooski, VT, USA) paired with Gen5™ software (BioTek Instruments, Inc., Winooski, VT, USA). Samples that had FUM levels predicted to be above the highest standard provided by the ELISA kits, 6 µg g−1 (ppm), were serially diluted until their predicted FUM levels were within the standard curve. To approximate samples that had non-detectable FUM levels (<0.1 ppm), uniform random values between 0 and 0.1 ppm were assigned to these samples [66]. The ratio of FUM to FER (FUM:FER) was calculated as FUM/(FER+1).

4.3. Mixed Models and Heritability Estimation

For mixed model analysis, we treated each of the five symptom types as absence/presence traits. Each plot was assigned as having the absence (0) or presence (1) of each symptom type. For example, if a plot exhibited the blush symptom type, it was assigned 1 for blush and 0 for asymptomatic, starburst, purple, and moldy. As the raw BDENinoc, FER, FUM, and FUM:FER data were not normally distributed, they were Box-Cox transformed using JMP®software (SAS Institute Inc., Cary, NC, USA, 1989–2007) for further analysis.
We used JMP®software (SAS Institute Inc., Cary, NC, USA, 1989–2007) to fit mixed linear models for the five symptom type absence/presence traits and the four quantitative traits (Box-Cox-transformed BDENinoc, FER, FUM, and FUM:FER). Within and among families, mixed models were fit for each of the nine FVI traits as a separate response variable with environment, block[environment], genotype, and genotype*environment as random effects and DTS as a fixed covariate. Variance components from each model were extracted and used to calculate broad-sense heritability (H) as
σ G 2 σ G 2 + σ G E 2 e h + σ ε 2 p h ,
where σ2G, σ2GE, and σ2ε are the genotype, genotype*environment, and error variances, respectively, and eh and ph are defined as
e h = n i = 1 n 1 e i ,
p h = n i = 1 n 1 p i ,
where n is the number of genotypes and ei and pi are the number of environments and plots for the ith genotype, respectively [67].
We also sought to assess associations between FVI and kernel bulk density (BDENuninoc) and cob morphology under non-inoculated conditions. Total kernel volume, cob diameter (CobDiam), cob length (CobLen), cob mass (CobMass), and ear mass had been previously measured in the NAM in five location-environments in 2006 (Aurora, NY, USA; Clayton, NC, USA; Homestead, FL, USA; Ponce, Puerto Rico; Urbana, IL, USA) [41,42], and these publicly accessible data for our four NAM families were provided by Panzea [42]. BDENuninoc was calculated as (ear mass—CobMass)/total kernel volume. Assuming cob shape to be cylindrical, cob volume (CobVol) was calculated as π*CobLen*(CobDiam/2)2. Cob density (CobDen) was calculated as CobMass/CobVol. To estimate genotype means for BDENuninoc, CobDen, CobMass, and CobVol, we fit mixed linear models within and among families for these four ear traits with genotype and environment as random effects and DTS as a fixed effect using JMP®software (SAS Institute Inc., Cary, NC, USA, 1989–2007). Broad-sense heritability (H) was calculated as
H = σ G 2 σ G 2 + σ ε 2 e ,
where e is the number of environments and σ2G and σ2ε are the genotype and error variances, respectively. Genotype means and broad-sense heritability for CobDiam and CobLen were accessed from Brown et al. [41].

4.4. Trait Correlation Analyses

Pairwise Pearson correlations among plot phenotypes and among trait best linear unbiased predictors (BLUPs) for quantitatively measured traits (BDENinoc, BDENuninoc, FER, FUM, FUM:FER) were calculated within families using JMP®software (SAS Institute Inc., Cary, NC, USA 1989–2007). Both plot-level (phenotypic) and BLUP (genetic) correlations were assessed among the FVI traits. Genetic correlations among BLUPs for BDENuninoc and FVI traits were also calculated.

4.5. Comparison of Disease Severity among Families

We tested differences in FVI severity among families using JMP®software (SAS Institute Inc., Cary, NC, USA, 1989–2007). We used analysis of variance (ANOVA) to test the effect of family on BDENinoc, FER, FUM, and FUM:FER, and pairwise two-tailed t-tests to asses differences in these four traits between families. A likelihood ratio χ2 test was used to compare the composition of symptom types among families.

4.6. QTL Mapping, Colocalization, and Allele Effects

The NAM population had been analyzed with genotyping-by-sequencing (GBS) technology [68]. A subset of 7386 GBS markers with 0.2 cM resolution had been selected by Olukolu et al. [69]. We extracted the 7386 GBS markers for the B73 × CML333, B73 × CML52, B73 × CML69, and B73 × NC358 NAM families and used them for stepwise regression.
Stepwise regression models were fit in TASSEL version 5.2.37 [70] to identify markers associated with BDENinoc, BDENuninoc, FER, FUM, FUM:FER, and the five symptom types (asymptomatic, blush, starburst, purple, moldy). For the single-family models, a marker significance threshold of 0.001 was used [71]. For the joint-family models, marker effects were nested within family and the significance threshold was set to 0.0001 [71]. QTL were defined by the confidence intervals (CIs) of markers selected by stepwise regression [70]. Within each family, QTL were considered colocalized if their CIs overlapped.
In addition to identifying colocalized QTL, we also sought to investigate allele effect relationships at specific colocalized QTL. We extracted allele effect estimates from the TASSEL stepwise regression model outputs [70] and then compared allele effects on traits that had colocalized QTL.

4.7. Characterization of Resistance QTL

We classified the QTL mapped in our study as (1) FVI-specific if the trait-specific or colocalized QTL was only associated with FVI trait(s) or (2) FVI-ear colocalized if the QTL was associated with both ear and FVI trait(s). With JMP®software, we fit a linear model including absolute marker effect as the response and trait, family, resistance locus type, and QTL confidence interval size as fixed effects. We determined in which genetic bin(s) each resistance QTL was located using bin physical position (AGPv2) information from MaizeGDB [72], and if the QTL overlapped the bin containing the centromere on the respective chromosome, the QTL was classified as pericentromeric. We then compared the proportion of pericentromeric QTL between the two classes of resistance loci. We amassed loci associated with resistance to F. verticillioides from nine publications [19,49,53,54,55,56,57,58,59] and determined which of our QTL overlapped with said previously described loci.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/11/2/86/s1, Table S1: Genotype-means correlations among kernel bulk density (BDENinoc), Fusarium ear rot (FER), fumonisin concentration (FUM), and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation, and kernel bulk density (BDENuninoc) and cob density (CobDen), diameter (CobDiame), length (CobLen), mass (CobMass), and volume (CobVol) under uninoculated conditions within four maize RIL families.; Table S2: Summary of family-specific and joint family stepwise regression models and QTL for kernel bulk density (BDENinoc), Fusarium ear rot (FER), fumonisin concentration (FUM), and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation, and kernel bulk density (BDENuninoc) and cob density (CobDen), diameter (CobDiam), length (CobLen), mass (CobMass), and volume (CobVol) under uninoculated conditions; Table S3: Physical positions (AGPv2) and associated traits of colocalized and trait-specific QTL from family-specific and joint family analyses.

Author Contributions

Conceptualization, L.M. and R.J.N.; Methodology, L.M., C.T.Z., P.J.B.-K., J.B.H. and R.J.N.; Formal Analysis, L.M.; Investigation, L.M., C.T.Z., D.E.M.M. and M.M.A.; Resources, P.J.B.-K., J.B.H. and R.J.N.; Data Curation, L.M.; Writing-Original Draft Preparation, L.M.; Writing-Review & Editing, L.M., P.J.B.-K., J.B.H. and R.J.N.; Visualization, L.M.; Supervision, R.J.N.; Funding Acquisition, P.J.B.-K., J.B.H. and R.J.N.

Funding

This research was funded by National Science Foundation Plant Genome Research Program grant number 1127076 and McKnight Foundation project number 80267.

Acknowledgments

The authors wish to thank Jason Brewer and Judith Kolkman for technical support and Nicolas Morales for data management.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Rheeder, J.P.; Marasas, W.F.O.; Vismer, H.F. Production of Fumonisin Analogs by Fusarium Species. Appl. Environ. Microbiol. 2002, 68, 2101–2105. [Google Scholar] [CrossRef] [PubMed]
  2. Mueller, D.S.; Wise, K.A.; Sisson, A.J.; Allen, T.W.; Bergstrom, G.C.; Bosley, D.B.; Bradley, C.A.; Broders, K.D.; Byamukama, E.; Chilvers, M.I.; et al. Corn Yield Loss Estimates Due to Diseases in the United States and Ontario, Canada from 2012 to 2015. Plant Health Prog. 2016, 17, 211–222. [Google Scholar] [CrossRef]
  3. Wild, C.P.; Gong, Y.Y. Mycotoxins and Human Disease: A Largely Ignored Global Health Issue. Carcinogenesis 2009, 31, 71–82. [Google Scholar] [CrossRef] [PubMed]
  4. Marasas, W.F.O.; Riley, R.T.; Hendricks, K.A.; Stevens, V.L.; Sadler, T.W.; Gelineau-van Waes, J.; Missmer, S.A.; Cabrera, J.; Torres, O.; Gelderblom, W.C.A.; et al. Fumonisins Disrupt Sphingolipid Metabolism, Folate Transport, and Neural Tube Development in Embryo Culture and in Vivo: A Potential Risk Factor for Human Neural Tube Defects among Populations Consuming Fumonisin-Contaminated Maize. J. Nutr. 2004, 134, 711–716. [Google Scholar] [CrossRef] [PubMed]
  5. Missmer, S.A.; Suarez, L.; Felkner, M.; Wang, E.; Alfred, H.; Rothman, K.J.; Hendricks, K.A. Exposure to Fumonisins and the Occurrence of Neural Tube Defects along the Texas-Mexico Border. Environ. Health Perspect. 2006, 114, 237–241. [Google Scholar] [CrossRef] [PubMed]
  6. Kimanya, M.E.; De Meulenaer, B.; Roberfroid, D.; Lachat, C. Fumonisin Exposure through Maize in Complementary Foods is Inversely Associated with Linear Growth of Infants in Tanzania. Mol. Nutr. Food Res. 2010, 54, 1659–1667. [Google Scholar] [CrossRef]
  7. Kellerman, T.; Marasas, W.; Thiel, P.; Gelderblom, W.C.A.; Cawood, M.; Coetzer, J. Leukoencephalomalacia in Two Horses Induced by Oral Dosing of Fumonisin B1. Onderstepoort J. Vet. Res. 1990, 57, 269–275. [Google Scholar]
  8. Colvin, B.M.; Harrison, L.R. Fumonisin-Induced Pulmonary Edema and Hydrothorax in Swine. Mycopathologia 1992, 117, 79–82. [Google Scholar] [CrossRef]
  9. Mutiga, S.K.; Were, V.; Hoffmann, V.; Harvey, J.W.; Milgroom, M.G.; Nelson, R.J. Extent and Drivers of Mycotoxin Contamination: Inferences from a Survey of Kenyan Maize Mills. Phytopathology 2014, 104, 1221–1231. [Google Scholar] [CrossRef]
  10. Bankole, S.A.; Mabekoje, O.O. Occurrence of Aflatoxins and Fumonisins in Preharvest Maize from South-Western Nigeria. Food Addit. Contam. 2004, 21, 251–255. [Google Scholar] [CrossRef]
  11. Fandohan, P.; Gnonlonfin, B.; Hell, K.; Marasas, W.F.O.; Wingfield, M.J. Natural Occurrence of Fusarium and Subsequent Fumonisin Contamination in Preharvest and Stored Maize in Benin, West Africa. Int. J. Food Microbiol. 2005, 99, 173–183. [Google Scholar] [CrossRef] [PubMed]
  12. Arino, A.; Herrera, M.; Juan, T.; Estopanan, G.; Carraminana, J.J.; Rota, C.; Herrera, A. Influence of Agricultural Practices on the Contamination of Maize by Fumonisin Mycotoxins. J. Food Prot. 2009, 72, 898–902. [Google Scholar] [CrossRef] [PubMed]
  13. Blandino, M.; Reyneri, A.; Vanara, F. Influence of Nitrogen Fertilization on Mycotoxin Contamination of Maize Kernels. Crop Prot. 2008, 27, 222–230. [Google Scholar] [CrossRef]
  14. Lanubile, A.; Maschietto, V.; Borrelli, V.M.; Stagnati, L.; Logrieco, A.F.; Marocco, A. Molecular Basis of Resistance to Fusarium Ear Rot in Maize. Front. Plant Sci. 2017, 8, 1–13. [Google Scholar] [CrossRef]
  15. Afolabi, C.G.; Bandyopadhyay, R.; Leslie, J.F.; Ekpo, E.J.A. Effect of Sorting on Incidence and Occurrence of Fumonisins and Fusarium verticillioides on Maize from Nigeria. J. Food Prot. 2006, 69, 2019–2023. [Google Scholar] [CrossRef]
  16. Stasiewicz, M.J.; Falade, T.D.O.; Mutuma, M.; Mutiga, S.K.; Harvey, J.J.W.; Fox, G.; Pearson, T.C.; Muthomi, J.W.; Nelson, R.J. Multi-Spectral Kernel Sorting to Reduce Aflatoxins and Fumonisins in Kenyan Maize. Food Control 2017, 78, 203–214. [Google Scholar] [CrossRef]
  17. Dombrink-Kurtzman, M.A.; Dvorak, T.J.; Barron, M.E.; Rooney, L.W. Effect of Nixtamalization (Alkaline Cooking) on Fumonisin-Contaminated Corn for Production of Masa and Tortillas. J. Agric. Food Chem. 2000, 48, 5781–5786. [Google Scholar] [CrossRef]
  18. Humpf, H.-U.; Voss, K.A. Effects of Thermal Food Processing on the Chemical Structure and Toxicity of Fumonisin Mycotoxins. Mol. Nutr. Food Res. 2004, 48, 255–269. [Google Scholar] [CrossRef]
  19. Robertson-Hoyt, L.A.; Jines, M.P.; Balint-Kurti, P.J.; Kleinschmidt, C.E.; White, D.G.; Payne, G.A.; Maragos, C.M.; Molna, T.L.; Holland, J.B. QTL Mapping for Fusarium Ear Rot and Fumonisin Contamination Resistance in Two Maize Populations. Crop Sci. 2006, 46, 1734–1743. [Google Scholar] [CrossRef]
  20. Zila, C.T.; Ogut, F.; Romay, M.C.; Gardner, C.A.; Buckler, E.S.; Holland, J.B. Genome-Wide Association Study of Fusarium Ear Rot Disease in the U.S.A. Maize Inbred Line Collection. BMC Plant Biol. 2014, 14, 372. [Google Scholar] [CrossRef]
  21. Zila, C.T.; Samayoa, L.F.; Santiago, R.; Butrón, A.; Holland, J.B. A Genome-Wide Association Study Reveals Genes Associated with Fusarium Ear Rot Resistance in a Maize Core Diversity Panel. G3 2013, 3, 2095–2104. [Google Scholar] [CrossRef] [PubMed]
  22. Poland, J.A.; Balint-Kurti, P.J.; Wisser, R.J.; Pratt, R.C.; Nelson, R.J. Shades of Gray: The World of Quantitative Disease Resistance. Trends Plant Sci. 2009, 14, 21–29. [Google Scholar] [CrossRef] [PubMed]
  23. Dall’Asta, C.; Giorni, P.; Cirlini, M.; Reverberi, M.; Gregori, R.; Ludovici, M.; Camera, E.; Fanelli, C.; Battilani, P.; Scala, V. Maize Lipids Play a Pivotal Role in the Fumonisin Accumulation. World Mycotoxin J. 2015, 8, 87–97. [Google Scholar] [CrossRef]
  24. Lanubile, A.; Ferrarini, A.; Maschietto, V.; Delledonne, M.; Marocco, A.; Bellin, D. Functional Genomic Analysis of Constitutive and Inducible Defense Responses to Fusarium verticillioides Infection in Maize Genotypes with Contrasting Ear Rot Resistance. BMC Genom. 2014, 15, 710. [Google Scholar] [CrossRef] [PubMed]
  25. Wang, Y.; Zhou, Z.; Gao, J.; Wu, Y.; Xia, Z.; Zhang, H.; Henry, R. The Mechanisms of Maize Resistance to Fusarium verticillioides by Comprehensive Analysis of RNA-Seq Data. Front. Plant Sci. 2016, 7, 1654. [Google Scholar] [CrossRef] [PubMed]
  26. Sampietro, D.A.; Fauguel, C.M.; Vattuone, M.A.; Presello, D.A.; Catalán, C.A.N. Phenylpropanoids from Maize Pericarp: Resistance Factors to Kernel Infection and Fumonisin Accumulation by Fusarium verticillioides. Eur. J. Plant Pathol. 2013, 135, 105–113. [Google Scholar] [CrossRef]
  27. Burr, S.; Fry, S. Feruloylated Arabinoxylans Are Oxidatively Cross-Linked by Extracellular Maize Peroxidase but Not by Horseradish Peroxidase. Mol. Plant 2009, 2, 883–892. [Google Scholar] [CrossRef]
  28. Blandino, M.; Reyneri, A. Effect of Maize Hybrid Maturity and Grain Hardness on Fumonisin and Zearalenone Contamination. Ital. J. Agron. 2008, 2, 107–117. [Google Scholar] [CrossRef]
  29. Sampietro, D.A.; Vattuone, M.A.; Presello, D.A.; Fauguel, C.M. The Pericarp and Its Surface Wax Layer in Maize Kernels as Resistance Factors to Fumonisin Accumulation by Fusarium verticillioides. Crop Prot. 2009, 28, 196–200. [Google Scholar] [CrossRef]
  30. Bluhm, B.H.; Woloshuk, C.P. Amylopectin Induces Fumonisin B1 Production by Fusarium verticillioides during Colonization of Maize Kernels. Mol. Plant-Microbe Interact. 2005, 18, 1333–1339. [Google Scholar] [CrossRef]
  31. Mutiga, S.K.; Morales, L.; Angwenyi, S.; Wainaina, J.; Harvey, J.; Das, B.; Nelson, R.J. Association between Agronomic Traits and Aflatoxin Accumulation in Diverse Maize Lines Grown under Two Soil Nitrogen Levels in Eastern Kenya. Field Crop. Res. 2017, 205, 124–134. [Google Scholar] [CrossRef]
  32. Shetty, P.H.; Bhat, R.V. A Physical Method for Segregation of Fumonisin-Contaminated Maize. Food Chem. 1999, 66, 371–374. [Google Scholar] [CrossRef]
  33. Morales, L.; Marino, T.P.; Wenndt, A.J.; Fouts, J.Q.; Holland, J.B.; Nelson, R.J. Dissecting Symptomatology and Fumonisin Contamination Produced by Fusarium verticillioides in Maize Ears. Phytopathology 2018, 108, 1475–1485. [Google Scholar] [CrossRef] [PubMed]
  34. Marino, T.P. Genomic Prediction and QTL Validation for Resistance to Fusarium Ear Rot and Fumonisin in Maize. Ph.D. Thesis, North Carolina State University, Raleigh, NC, USA, October 2018. [Google Scholar]
  35. Horne, D.W.; Eller, M.S.; Holland, J.B. Responses to Recurrent Index Selection for Reduced Fusarium Ear Rot and Lodging and for Increased Yield in Maize. Crop Sci. 2016, 56, 85–94. [Google Scholar] [CrossRef]
  36. Jansen, C.; Zhang, Y.; Liu, H.; Gonzalez-Portilla, P.J.; Lauter, N.; Kumar, B.; Trucillo-Silva, I.; Martin, J.P.S.; Lee, M.; Simcox, K.; et al. Genetic and Agronomic Assessment of Cob Traits in Corn under Low and Normal Nitrogen Management Conditions. Theor. Appl. Genet. 2015, 128, 1231–1242. [Google Scholar] [CrossRef] [PubMed]
  37. Parsons, M.W.; Munkvold, G.P. Effects of Planting Date and Environmental Factors on Fusarium Ear Rot Symptoms and Fumonisin B1 Accumulation in Maize Grown in Six North American Locations. Plant Pathol. 2012, 61, 1130–1142. [Google Scholar] [CrossRef]
  38. Parsons, M.W.; Munkvold, G.P. Associations of Planting Date, Drought Stress, and Insects with Fusarium Ear Rot and Fumonisin B1 Contamination in California Maize. Food Addit. Contam. 2010, 27, 591–607. [Google Scholar] [CrossRef]
  39. Desjardins, A.E.; Plattner, R.D.; Lu, M.; Claflin, L.W. Distribution of Fumonisins in Maize Ears Infected with Strains of Fusarium moniliforme That Differ in Fumonisin Production. Plant Dis. 1998, 82, 953–958. [Google Scholar] [CrossRef]
  40. Yu, J.; Holland, J.B.; McMullen, M.D.; Buckler, E.S. Genetic Design and Statistical Power of Nested Association Mapping in Maize. Genetics 2008, 178, 539–551. [Google Scholar] [CrossRef]
  41. Brown, P.J.; Upadyayula, N.; Mahone, G.S.; Tian, F.; Bradbury, P.J.; Myles, S.; Holland, J.B.; Flint-Garcia, S.; McMullen, M.D.; Buckler, E.S.; et al. Distinct Genetic Architectures for Male and Female Inflorescence Traits of Maize. PLoS Genet. 2011, 7, e1002383. [Google Scholar] [CrossRef]
  42. Zhao, W. Panzea: A Database and Resource for Molecular and Functional Diversity in the Maize Genome. Nucleic Acids Res. 2006, 34, D752–D757. [Google Scholar] [CrossRef] [PubMed]
  43. Robertson, L.A.; Kleinschmidt, C.E.; White, D.G.; Payne, G.A.; Maragos, C.M.; Holland, J.B. Heritabilities and Correlations of Fusarium Ear Rot Resistance and Fumonisin Contamination Resistance in Two Maize Populations. Crop Sci. 2006, 46, 353–361. [Google Scholar] [CrossRef]
  44. Hung, H.Y.; Holland, J.B. Diallel Analysis of Resistance to Fusarium Ear Rot and Fumonisin Contamination in Maize. Crop Sci. 2012, 52, 2173–2181. [Google Scholar] [CrossRef]
  45. Zummo, N.; Scott, G.E. Cob and Kernel Infection by Aspergillus flavus and Fusarium moniliforme in Inoculated, Field-Grown Maize Ears. Plant Dis. 1990, 74, 627–631. [Google Scholar] [CrossRef]
  46. Presello, D.A.; Botta, G.; Iglesias, J.; Eyhérabide, G.H. Effect of Disease Severity on Yield and Grain Fumonisin Concentration of Maize Hybrids Inoculated with Fusarium verticillioides. Crop Prot. 2008, 27, 572–576. [Google Scholar] [CrossRef]
  47. Presello, D.A.; Iglesias, J.; Botta, G.; Eyherabide, G.H. Severity of Fusarium Ear Rot and Concentration of Fumonisin in Grain of Argentinian Maize Hybrids. Crop Prot. 2007, 26, 852–855. [Google Scholar] [CrossRef]
  48. Robertson-Hoyt, L.A.; Betrán, J.; Payne, G.A.; White, D.G.; Isakeit, T.; Maragos, C.M.; Molnár, T.L.; Holland, J.B. Relationships Among Resistances to Fusarium and Aspergillus Ear Rots and Contamination by Fumonisin and Aflatoxin in Maize. Phytopathology 2007, 97, 311–317. [Google Scholar] [CrossRef] [PubMed]
  49. Mideros, S.X.; Warburton, M.L.; Jamann, T.M.; Windham, G.L.; Williams, W.P.; Nelson, R.J. Quantitative Trait Loci Influencing Mycotoxin Contamination of Maize: Analysis by Linkage Mapping, Characterization of Near-Isogenic Lines, and Meta-Analysis. Crop Sci. 2014, 54, 127. [Google Scholar] [CrossRef]
  50. Fox, G.; Manley, M. Hardness Methods for Testing Maize Kernels. J. Agric. Food Chem. 2009, 57, 5647–5657. [Google Scholar] [CrossRef]
  51. Kandianis, C.B.; Michenfelder, A.S.; Simmons, S.J.; Grusak, M.A.; Stapleton, A.E. Abiotic Stress Growth Conditions Induce Different Responses in Kernel Iron Concentration across Genotypically Distinct Maize Inbred Varieties. Front. Plant Sci. 2013, 4, 1–10. [Google Scholar] [CrossRef]
  52. Rodgers-Melnick, E.; Bradbury, P.J.; Elshire, R.J.; Glaubitz, J.C.; Acharya, C.B.; Mitchell, S.E.; Li, C.; Li, Y.; Buckler, E.S. Recombination in Diverse Maize Is Stable, Predictable, and Associated with Genetic Load. Proc. Natl. Acad. Sci. USA 2015, 112, 201413864. [Google Scholar] [CrossRef] [PubMed]
  53. Chen, J.; Shrestha, R.; Ding, J.; Zheng, H.; Mu, C.; Wu, J.; Mahuku, G. Genome-Wide Association Study and QTL Mapping Reveal Genomic Loci Associated with Fusarium Ear Rot Resistance in Tropical Maize Germplasm. G3 2016, 8651, g3-116. [Google Scholar] [CrossRef] [PubMed]
  54. Eller, M.S.; Payne, G.A.; Holland, J.B. Selection for Reduced Fusarium Ear Rot and Fumonisin Content in Advanced Backcross Maize Lines and Their Topcross Hybrids. Crop Sci. 2010, 50, 2249–2260. [Google Scholar] [CrossRef]
  55. Ding, J.-Q.; Wang, X.-M.; Chander, S.; Yan, J.-B.; Li, J.-S. QTL Mapping of Resistance to Fusarium Ear Rot Using a RIL Population in Maize. Mol. Breed. 2008, 22, 395–403. [Google Scholar] [CrossRef]
  56. Xiang, K.; Zhang, Z.M.; Reid, L.M.; Zhu, X.Y.; Yuan, G.S.; Pan, G.T. A Meta-Analysis of QTL Associated with Ear Rot Resistance in Maize. Maydica 2010, 55, 281–290. [Google Scholar]
  57. Maschietto, V.; Colombi, C.; Pirona, R.; Pea, G.; Strozzi, F.; Marocco, A.; Rossini, L.; Lanubile, A. QTL Mapping and Candidate Genes for Resistance to Fusarium Ear Rot and Fumonisin Contamination in Maize. BMC Plant Biol. 2017, 17, 20. [Google Scholar] [CrossRef] [PubMed]
  58. Giomi, G.M.; Kreff, E.D.; Iglesias, J.; Fauguel, C.M.; Fernandez, M.; Oviedo, M.S.; Presello, D.A. Quantitative Trait Loci for Fusarium and Gibberella Ear Rot Resistance in Argentinian Maize Germplasm. Euphytica 2016, 211, 287–294. [Google Scholar] [CrossRef]
  59. Perez-Brito, D.; Jeffers, D.; Gonzalez-de-Leon, D.; Khairallah, M.; Cortes-Cruz, M.; Velazquez-Cardelas, G.; Azpiroz-Rivero, S.; Srinivasan, G. QTL Mapping of Fusarium moniliforme Ear Rot Resistance in High Land Maize, Mexico. Agrociencia 2001, 35, 181–196. [Google Scholar]
  60. Clements, M.J.; Kleinschmidt, C.E.; Maragos, C.M.; Pataky, J.K.; White, D.G. Evaluation of Inoculation Techniques for Fusarium Ear Rot and Fumonisin Contamination of Corn. Plant Dis. 2003, 87, 147–153. [Google Scholar] [CrossRef]
  61. Gustin, J.L.; Jackson, S.; Williams, C.; Patel, A.; Armstrong, P.; Peter, G.F.; Settles, A.M. Analysis of Maize (Zea mays) Kernel Density and Volume Using Microcomputed Tomography and Single-Kernel Near-Infrared Spectroscopy. J. Agric. Food Chem. 2013, 61, 10872–10880. [Google Scholar] [CrossRef]
  62. Shelby, R.A.; White, D.G.; Bauske, E.M. Differential Fumonisin Production in Maize Hybrids. Plant Dis. 1994, 78, 582–584. [Google Scholar] [CrossRef]
  63. Buckler, E.S.; Holland, J.B.; Bradbury, P.J.; Acharya, C.B.; Brown, P.J.; Browne, C.; Ersoz, E.; Flint-Garcia, S.; Garcia, A.; Glaubitz, J.C.; et al. The Genetic Architecture of Maize Flowering Time. Science 2009, 325, 714–718. [Google Scholar] [CrossRef] [PubMed]
  64. Costa, R.S.; Môro, F.V.; Môro, J.R.; da Silva, H.P.; Panizzi, R.D.C. Relação Entre Características Morfológicas Da Cariopse e Fusariose Em Milho. Pesqui. Agropecu. Bras. 2003, 38, 27–33. [Google Scholar] [CrossRef]
  65. Duncan, K.E.; Howard, R.J. Biology of Maize Kernel Infection by Fusarium verticillioides. Mol. Plant-Microbe Interact. 2010, 23, 6–16. [Google Scholar] [CrossRef] [PubMed]
  66. Owens, B.; Lipka, A.; Magallanes-Lundback, M.; Tiede, T.; Diepenbrock, C.; Kandianis, C.B.; Kim, E.; Cepela, J.; Mateos-Hernandez, M.; Buell, C.; et al. A Foundation for Provitamin A Biofortification of Maize: Genome-Wide Association and Genomic Prediction Models of Carotenoid Levels. Genetics 2014, 198, 1699–1716. [Google Scholar] [CrossRef] [PubMed]
  67. Holland, J.B.; Nyquist, W.E.; Cervantes-Martinez, C.T. Estimating and Interpreting Heritability for Plant Breeding. Plant Breed. Rev. 2003, 22, 9–112. [Google Scholar]
  68. Elshire, R.J.; Glaubitz, J.C.; Sun, Q.; Poland, J.A.; Kawamoto, K.; Buckler, E.S.; Mitchell, S.E. A Robust, Simple Genotyping-by-Sequencing (GBS) Approach for High Diversity Species. PLoS ONE 2011, 6, e19379. [Google Scholar] [CrossRef]
  69. Olukolu, B.A.; Wang, G.F.; Vontimitta, V.; Venkata, B.P.; Marla, S.; Ji, J.; Gachomo, E.; Chu, K.; Negeri, A.; Benson, J.; et al. A Genome-Wide Association Study of the Maize Hypersensitive Defense Response Identifies Genes That Cluster in Related Pathways. PLoS Genet. 2014, 10, e1004562. [Google Scholar] [CrossRef]
  70. Bradbury, P.J.; Zhang, Z.; Kroon, D.E.; Casstevens, T.M.; Ramdoss, Y.; Buckler, E.S. TASSEL: Software for Association Mapping of Complex Traits in Diverse Samples. Bioinformatics 2007, 23, 2633–2635. [Google Scholar] [CrossRef]
  71. Peiffer, J.A.; Romay, M.C.; Gore, M.A.; Flint-Garcia, S.A.; Zhang, Z.; Millard, M.J.; Gardner, C.A.C.; McMullen, M.D.; Holland, J.B.; Bradbury, P.J.; et al. The Genetic Architecture of Maize Height. Genetics 2014, 196, 1337–1356. [Google Scholar] [CrossRef]
  72. Andorf, C.M.; Cannon, E.K.; Portwood, J.L.; Gardiner, J.M.; Harper, L.C.; Schaeffer, M.L.; Braun, B.L.; Campbell, D.A.; Vinnakota, A.G.; Sribalusu, V.V.; et al. MaizeGDB Update: New Tools, Data and Interface for the Maize Model Organism Database. Nucleic Acids Res. 2016, 44, D1195–D1201. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Proportions of plots exhibiting distinct symptom types in four RIL families across three environments.
Figure 1. Proportions of plots exhibiting distinct symptom types in four RIL families across three environments.
Toxins 11 00086 g001
Figure 2. (A) Box-plots comparing H, number of associated quantitative trait loci (QTL), and variance explained by QTL models (model R2) among BDENinoc, FER, FUM, and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation, and kernel bulk density (BDENuninoc) and cob density, diameter, length, mass, and volume under uninoculated conditions. (B) Pearson correlations among H, number of associated QTL, and QTL model R2 for the previously described F. verticillioides infection (FVI) and ear traits (p < 0.0001 for all three pairwise correlations).
Figure 2. (A) Box-plots comparing H, number of associated quantitative trait loci (QTL), and variance explained by QTL models (model R2) among BDENinoc, FER, FUM, and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation, and kernel bulk density (BDENuninoc) and cob density, diameter, length, mass, and volume under uninoculated conditions. (B) Pearson correlations among H, number of associated QTL, and QTL model R2 for the previously described F. verticillioides infection (FVI) and ear traits (p < 0.0001 for all three pairwise correlations).
Toxins 11 00086 g002
Table 1. Comparison of kernel bulk density (BDENinoc), Fusarium ear rot (FER), fumonisin concentration (FUM), and the ratio of FUM to FER (FUM:FER) under Fusarium verticillioides inoculation among four nested association mapping (NAM) recombinant inbred line (RIL) families and among the five parental inbred lines. Analyses included (1) ANOVA for BDENinoc, FER, FUM, and FUM:FER vs. family or parent, and (2) pairwise t-tests comparing BDENinoc, FER, FUM, and FUM:FER between families or parents. FER was evaluated across four environments, FUM and FUM:FER across three environments, and BDENinoc across two environments.
Table 1. Comparison of kernel bulk density (BDENinoc), Fusarium ear rot (FER), fumonisin concentration (FUM), and the ratio of FUM to FER (FUM:FER) under Fusarium verticillioides inoculation among four nested association mapping (NAM) recombinant inbred line (RIL) families and among the five parental inbred lines. Analyses included (1) ANOVA for BDENinoc, FER, FUM, and FUM:FER vs. family or parent, and (2) pairwise t-tests comparing BDENinoc, FER, FUM, and FUM:FER between families or parents. FER was evaluated across four environments, FUM and FUM:FER across three environments, and BDENinoc across two environments.
FamilyRILs (N/env)BDENinoc ***
(g mL−1)
FER***
(%)
FUM***
(ppm)
FUM:FER ***
(ppm %−1)
B73 × CML3331860.737 ± 0.003 A43.04 ± 1.27 A24.24 ± 5.32 B0.51 ± 0.10 C
B73 × CML521770.720 ± 0.004 B34.41 ± 1.05 B11.77 ± 2.39 B0.60 ± 0.08 B
B73 × CML691860.712 ± 0.003 B40.04 ± 1.23 A31.28 ± 5.05 A1.00 ± 0.13 AB
B73 × NC3581790.683 ± 0.004 C43.27 ± 1.20 A32.62 ± 4.21 A1.19 ± 0.23 A
ParentPlots (N/env)BDENinoc ***
(g mL−1)
FER***
(%)
FUM***
(ppm)
FUM:FER
(ppm %−1)
B73390.650 ± 0.005 C47.83 ± 2.26 A72.93 ± 14.41 A2.36 ± 0.64 A
CML333110.770 ± 0.008 A26.65 ± 3.37 B7.09 ± 2.82 B0.47 ± 0.18 B
CML52110.714 ± 0.020 B13.39 ± 3.17 C3.18 ± 0.78 B0.76 ± 0.26 AB
CML69110.745 ± 0.004 AB25.42 ± 5.89 B6.69 ± 2.08 B0.48 ± 0.20 AB
NC358110.766 ± 0.007 A26.18 ± 3.67 B16.14 ± 7.36 B0.63 ± 0.25 AB
N/env denotes the number of RILs evaluated in each family per environment or the number of plots for each parental line per environment. Raw means and standard errors (SE) are reported, but ANOVA and pairwise t-tests were conducted on Box-Cox transformed data. Groups not connected by the same letter within each column are significantly different (pairwise two-tailed t-tests, p < 0.05). *** p < 0.0001 for ANOVA of trait in column header vs. family.
Table 2. Comparison of BDENinoc, FER, FUM, and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation among plots exhibiting distinct symptom types. Combined and family-specific analyses included (1) ANOVA for BDENinoc, FER, FUM, and FUM:FER vs. symptom type, and (2) pairwise t-tests comparing BDENinoc, FER, FUM, and FUM:FER between symptom types.
Table 2. Comparison of BDENinoc, FER, FUM, and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation among plots exhibiting distinct symptom types. Combined and family-specific analyses included (1) ANOVA for BDENinoc, FER, FUM, and FUM:FER vs. symptom type, and (2) pairwise t-tests comparing BDENinoc, FER, FUM, and FUM:FER between symptom types.
FamilySymptom
type
BDENinoc
(g mL−1)
FER
(%)
FUM
(ppm)
FUM:FER
(ppm %−1)
CombinedAsym.0.693 ± 0.003 A19.58 ± 0.44 E10.79 ± 1.69 A0.96 ± 0.009 A
Blush0.682 ± 0.005 A56.59 ± 1.88 C15.51 ± 6.64 C0.73 ± 0.017 C
Starburst0.656 ± 0.003 B65.54 ± 0.75 B18.40 ± 1.90 B0.75 ± 0.008 C
Purple0.641 ± 0.003 C50.51 ± 0.82 D38.39 ± 4.75 BC0.78 ± 0.009 B
Moldy0.576 ± 0.013 D76.14 ± 2.64 A97.66 ± 38.54 A0.79 ± 0.038 BC
ANOVA p<0.0001 ***<0.0001 ***0.0002 **<0.0001 ***
B73 ×Asym.0.713 ± 0.007 A19.73 ± 1.36 C4.58 ± 0.92 AB0.91 ± 0.02 A
CML333Blush0.691 ± 0.005 B58.68 ± 2.20 B14.35 ± 5.01 AB0.74 ± 0.02 B
Starburst0.685 ± 0.006 B70.49 ± 1.49 A13.80 ± 3.12 A0.73 ± 0.02 B
Purple0.658 ± 0.006 C54.04 ± 1.79 B31.19 ± 10.31 B0.71 ± 0.02 B
Moldy0.601 ± 0.021 D77.95 ± 3.77 A85.22 ± 53.51 AB0.73 ± 0.07 B
ANOVA p<0.0001 ***<0.0001 ***0.0427 *0.0001 **
B73 ×Asym.0.688 ± 0.005 A18.09 ± 0.68 D6.17 ± 0.70 A0.94 ± 0.01 A
CML52Blush0.606 ± 0.020 B46.78 ± 7.61 BC8.75 ± 8.17 AB0.70 ± 0.07 BC
Starburst0.660 ± 0.004 B59.81 ± 1.32 B10.46 ± 2.16 AB0.74 ± 0.01 C
Purple0.655 ± 0.008 B41.05 ± 1.57 C21.21 ± 7.81 B0.78 ± 0.02 B
Moldy0.503 ± 0.039 C84.19 ± 5.38 A32.09 ± 24.56 AB0.66 ± 0.08 BC
ANOVA p<0.0001 ***<0.0001 ***0.113<0.0001 ***
B73 ×Asym.0.693 ± 0.004 A21.47 ± 0.79 C14.70 ± 4.11 AB0.99 ± 0.01 A
CML69Blush0.659 ± 0.016 AB46.36 ± 6.61 B39.96 ± 25.15 BC0.78 ± 0.07 BC
Starburst0.657 ± 0.005 B67.74 ± 1.40 A21.59 ± 4.09 C0.75 ± 0.02 C
Purple0.644 ± 0.006 B54.22 ± 1.63 B42.20 ± 8.61 C0.78 ± 0.02 BC
Moldy0.588 ± 0.021 C72.49 ± 5.74 A179.86 ± 116.91 A0.88 ± 0.07 AB
ANOVA p<0.0001 ***<0.0001 ***0.015 *<0.0001 ***
B73 ×Asym.0.692 ± 0.006 A20.88 ± 0.85 D21.91 ± 7.24 A0.97 ± 0.02 A
NC358Blush0.654 ± 0.016 AB56.06 ± 4.60 BC4.93 ± 1.84 B0.67 ± 0.05 D
Starburst0.617 ± 0.006 B66.45 ± 1.68 A29.82 ± 5.56 A0.77 ± 0.02 CD
Purple0.623 ± 0.006 B51.19 ± 1.40 C53.31 ± 10.22 A0.82 ± 0.02 B
Moldy0.562 ± 0.029 C68.74 ± 7.55 AB60.22 ± 31.84 A0.89 ± 0.08 ABC
ANOVA p<0.0001 ***<0.0001 ***0.0502<0.0001 ***
ANOVA and pairwise t-tests were conducted on Box-Cox transformed data. Box-Cox-transformed means and standard errors (±) are reported. Groups not connected by the same letter within each column/family are significantly different (pairwise two-tailed t-tests, p < 0.05). ANOVA significance is denoted as * 0.05 > p ≥ 0.01; ** 0.01 > p ≥ 0.0001; *** p < 0.0001.
Table 3. Phenotypic and genotype-mean correlations among BDENinoc, FER, FUM, and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation within four NAM RIL families. Phenotypic correlations (rP) among plot-level phenotypes are in the upper diagonals. Correlations among genotype-means (rG) after controlling for genetic and field effects are in the lower diagonals.
Table 3. Phenotypic and genotype-mean correlations among BDENinoc, FER, FUM, and the ratio of FUM to FER (FUM:FER) under F. verticillioides inoculation within four NAM RIL families. Phenotypic correlations (rP) among plot-level phenotypes are in the upper diagonals. Correlations among genotype-means (rG) after controlling for genetic and field effects are in the lower diagonals.
B73 × CML333 Family
rP/rGBDENinocFERFUMFUM:FER
BDENinoc −0.46 ***−0.25 ***−0.15 *
FER−0.53 *** 0.17 **−0.18 **
FUM−0.31 **0.28 *** 0.94 ***
FUM:FER0.02−0.33***0.80 ***
B73 × CML52 family
rP/rGBDENinocFERFUMFUM:FER
BDENinoc −0.47 ***−0.12 *0.12 *
FER−0.48 *** 0.08 ms−0.36 ***
FUM−0.18 *0.20 * 0.90 ***
FUM:FER0.15 *−0.46 ***0.75 ***
B73 × CML69 family
rP/rGBDENinocFERFUMFUM:FER
BDENinoc −0.55 ***−0.30 ***0.02
FER−0.53 *** −0.01−0.33 ***
FUM−0.27 **0.21** 0.95 ***
FUM:FER0.06−0.39 ***0.81 ***
B73 × NC358 family
rP/rGBDENinocFERFUMFUM:FER
BDENinoc −0.59 ***−0.19 **0.11 ms
FER−0.58 *** 0.21 ***−0.14 **
FUM−0.25 **0.12 0.94 ***
FUM:FER0.15 ms−0.46 ***0.79 ***
Pearson correlation coefficients are reported in each cell, and significance is denoted as ms 0.1 > p ≥ 0.05 (marginally significant); * 0.05 > p ≥ 0.01; ** 0.01 > p ≥0.0001; *** p < 0.0001.
Table 4. Genotypic and field effect variance proportions, significance of days to silking (DTS), and broad-sense heritability of genotype-means (H) for BDENinoc, FER, FUM, the ratio of FUM to FER (FUM:FER), and the absence/presence of five symptom types (asymptomatic, blush, starburst, purple, moldy) under F. verticillioides inoculation within four NAM RIL families. Mixed linear models were fit with genotype, environment (env), genotype*environment (G*E), block[environment] (B[E]) as random effects and DTS as a fixed effect.
Table 4. Genotypic and field effect variance proportions, significance of days to silking (DTS), and broad-sense heritability of genotype-means (H) for BDENinoc, FER, FUM, the ratio of FUM to FER (FUM:FER), and the absence/presence of five symptom types (asymptomatic, blush, starburst, purple, moldy) under F. verticillioides inoculation within four NAM RIL families. Mixed linear models were fit with genotype, environment (env), genotype*environment (G*E), block[environment] (B[E]) as random effects and DTS as a fixed effect.
FamilyTraitRandom Effects
(Variance Proportions)
Fixed Effect
(p-Value)
H
GenotypeEnv.G*EB[E]ErrorDTS
B73 ×
CML333
BDENinoc0.640.010.0400.310.10.79
FER0.420.010.1200.450.90.75
FUM0.020.610.1600.210.30.12
FUM:FER0.030.580.1000.290.60.17
Asymptomatic0.050.070.170.020.690.20.15
Blush0.180.130.300.010.380.004 **0.44
Starburst0.120.180.170.010.530.90.33
Purple0.040.250.120.030.560.50.16
Moldy0.170.030.8000.0010.30.40
B73 ×
CML52
BDENinoc0.330.00300.040.620.002 **0.52
FER0.300.200.060.030.410.03 *0.72
FUM0.030.460.010.040.460.10.17
FUM:FER0.020.4800.010.480.30.12
Asymptomatic0.060.270.170.030.480.80.21
Blush00.01000.990.04 *0
Starburst0.110.280.180.020.420.60.35
Purple0.030.300.200.010.460.40.13
Moldy0.310.010.6801 × 10−50.80.58
B73 ×
CML69
BDENinoc0.460.030.320.020.180.04 *0.65
FER0.220.140.180.030.430.80.59
FUM0.020.700.130.010.140.10.15
FUM:FER0.020.700.070.020.190.30.21
Asymptomatic0.090.130.1600.620.20.25
Blush0.060.040.0100.880.050.18
Starburst0.080.240.150.0030.530.005 **0.27
Purple0.070.260.210.010.460.03*0.24
Moldy0.170.030.8008 × 10−50.009 **0.39
B73 ×
NC358
BDENinoc0.230.080.560.020.120.070.40
FER0.270.220.060.040.400.003 **0.70
FUM0.010.730.0400.220.70.14
FUM:FER0.030.650.0600.270.70.19
Asymptomatic0.050.0400.050.860.080.15
Blush0.110.020.8700.0010.90.28
Starburst0.060.200.030.010.700.003 **0.21
Purple0.070.2500.010.670.03 *0.24
Moldy0.130.010.8500.0020.40.32
Fixed effect p-value significance is denoted as: * 0.05 > p ≥ 0.01; ** 0.01 > p ≥ 0.0001.
Table 5. Physical positions (AGPv2) and marker effects of colocalized QTL. Marker effects were nested within family in the joint family stepwise regression model. As such, each colocalized QTL from the joint family analysis has four marker effects from the four non-B73 alleles (CML333, CML52, CML69, NC358) vs. the B73 allele. Colocalized QTL from the single-family analyses have one marker effect each (non-B73 parent allele vs. B73 allele).
Table 5. Physical positions (AGPv2) and marker effects of colocalized QTL. Marker effects were nested within family in the joint family stepwise regression model. As such, each colocalized QTL from the joint family analysis has four marker effects from the four non-B73 alleles (CML333, CML52, CML69, NC358) vs. the B73 allele. Colocalized QTL from the single-family analyses have one marker effect each (non-B73 parent allele vs. B73 allele).
AnalysisChr.Position
(bp)
Associated
Traits
Marker Effect on Trait
(Non-B73 vs. B73 Allele)
Joint117,256,105–BDENinoc+0.003CML333; +0.013CML52; +0.019CML69; −0.001NC358
family 25,068,060Blush+0.054CML333; +0.002CML52; −0.001CML69; +0.007NC358
CobLen−4.874CML333; −3.899CML52; −5.097CML69; −13.611NC358
220,081,914–Blush+0.073CML333; +0.001CML52; +0.012CML69; +0.005NC358
50,634,516CobDiam+0.531CML333; +0.380CML52; +0.510CML69; +0.907NC358
CobLen+3.033CML333; +6.651CML52; +3.954CML69; +9.236NC358
CobMass+1.197CML333; +0.655CML52; +1.605CML69; +1.843NC358
CobVol+5.456CML333; +3.441CML52; +4.325CML69; +8.825NC358
Moldy−0.026CML333; −0.013CML52; −0.021CML69; −0.006NC358
Starburst−0.001CML333; +0.023CML52; +0.044CML69; +0.047NC358
2193,400,572–BDENinoc+0.025CML333; +0.003CML52; +0.010CML69; +0.002NC358
196,813,641CobMass−1.545CML333; +0.418CML52; −1.314CML69; +0.220NC358
5205,258,663–BDENinoc+0.009CML333; −0.016CML52; −0.002CML69; −0.016NC358
208,297,027Blush−0.088CML333; −0.001CML52; −0.006CML69; −0.017NC358
CobLen−0.868CML333; −2.017CML52; −7.233CML69; −2.325NC358
FER−0.606CML333; +5.843CML52; −1.058CML69; +8.576NC358
6130,721,405–Asym+0.002CML333; +0.027CML52; +0.015CML69; +0.007NC358
133,795,910FER−4.192CML333; −6.906CML52; −6.253CML69; −4.148NC358
1079,389,587–CobDen+0.005CML333; +0.012CML52; +0.0003CML69; +0.007NC358
116,043,385Starburst−0.031CML333; −0.037CML52; −0.029CML69; −0.048NC358
B73 ×119,128,807–Blush+0.052CML333
CML333 24,311,471CobLen−5.937CML333
147,813,030–FER+6.660CML333
57,347,172FUM+0.162CML333
221,767,606–BDENinoc+0.020CML333
35,311,519Blush+0.065CML333
74,823,956–CobLen+6.465CML333
6,178,479CobVol+15.265CML333
Starburst+0.046CML333
714,152,199–FER+4.985CML333
23,748,237FUM:FER−0.009CML333
8130,983,626–BDENinoc+0.015CML333
151,773,181CobDen+0.010CML333
CobDiam−1.587CML333
Starburst+0.040CML333
B73 ×583,007,601–BDENinoc+0.018CML52
CML52 140,575,190BDENuninoc+0.009CML52
6129,994,297–Asym+0.027CML52
140,630,780CobDiam−0.744CML52
6154,530,689–FER−7.612CML52
156,104,230FUM:FER+0.010CML52
8132,405,651–BDENuninoc+0.010CML52
146,365,593CobVol−6.924CML52
FER−5.718CML52
B73 ×249,846,838–CobMass+2.016CML69
CML69 87,089,680Moldy−0.023CML69
Starburst+0.048CML69
417,330,826–BDENinoc−0.018CML69
21,710,411FER+6.347CML69
4167,066,431–Purple+0.034CML69
170,804,182Starburst−0.050CML69
7138,072,861–Asym+0.020CML69
142,429,440FER−7.435CML69
B73 ×2144,159,847–Purple−0.031NC358
NC358 160,579,525Starburst+0.046NC358
36,427,177–BDENinoc+0.021NC358
8,563,589Starburst−0.054NC358
1016,505,881–Asym+0.018NC358
77,678,052CobDen+0.009NC358

Share and Cite

MDPI and ACS Style

Morales, L.; Zila, C.T.; Moreta Mejía, D.E.; Montoya Arbelaez, M.; Balint-Kurti, P.J.; Holland, J.B.; Nelson, R.J. Diverse Components of Resistance to Fusarium verticillioides Infection and Fumonisin Contamination in Four Maize Recombinant Inbred Families. Toxins 2019, 11, 86. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11020086

AMA Style

Morales L, Zila CT, Moreta Mejía DE, Montoya Arbelaez M, Balint-Kurti PJ, Holland JB, Nelson RJ. Diverse Components of Resistance to Fusarium verticillioides Infection and Fumonisin Contamination in Four Maize Recombinant Inbred Families. Toxins. 2019; 11(2):86. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11020086

Chicago/Turabian Style

Morales, Laura, Charles T. Zila, Danilo E. Moreta Mejía, Melissa Montoya Arbelaez, Peter J. Balint-Kurti, James B. Holland, and Rebecca J. Nelson. 2019. "Diverse Components of Resistance to Fusarium verticillioides Infection and Fumonisin Contamination in Four Maize Recombinant Inbred Families" Toxins 11, no. 2: 86. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins11020086

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