Skip to main content

Linkage mapping and genome-wide association study reveals conservative QTL and candidate genes for Fusarium rot resistance in maize

Abstract

Background

Fusarium ear rot (FER) caused by Fusarium verticillioides is a major disease of maize that reduces grain yield and quality globally. However, there have been few reports of major loci for FER were verified and cloned.

Result

To gain a comprehensive understanding of the genetic basis of natural variation in FER resistance, a recombinant inbred lines (RIL) population and one panel of inbred lines were used to map quantitative trait loci (QTL) for resistance. As a result, a total of 10 QTL were identified by linkage mapping under four environments, which were located on six chromosomes and explained 1.0–7.1% of the phenotypic variation. Epistatic mapping detected four pairs of QTL that showed significant epistasis effects, explaining 2.1–3.0% of the phenotypic variation. Additionally, 18 single nucleotide polymorphisms (SNPs) were identified across the whole genome by genome-wide association study (GWAS) under five environments. Compared linkage and association mapping revealed five common intervals located on chromosomes 3, 4, and 5 associated with FER resistance, four of which were verified in different near-isogenic lines (NILs) populations. GWAS identified three candidate genes in these consistent intervals, which belonged to the Glutaredoxin protein family, actin-depolymerizing factors (ADFs), and AMP-binding proteins. In addition, two verified FER QTL regions were found consistent with Fusarium cob rot (FCR) and Fusarium seed rot (FSR).

Conclusions

These results revealed that multi pathways were involved in FER resistance, which was a complex trait that was controlled by multiple genes with minor effects, and provided important QTL and genes, which could be used in molecular breeding for resistance.

Background

Fusarium ear rot (FER) is one of the most important food and feed safety challenges in global maize production [1]. FER not only reduces the yield and quality of harvested maize but also is fatal to humans and animals, which consume the contaminated grain containing mycotoxins from some of the Fusarium spp. [2]. More than 10 Fusarium spp. can cause ear rot, among them, Fusarium verticillioides [synonym F. moniliforme Sheldon] and F. graminearum are the two most important species which can cause FER and Gibberella ear rot, respectively [3,4,5].

Fusarium verticillioides is an important maize pathogen in the world, which can lead to serious economic losses [6], particularly in China [7,8,9], the United States [10] and Southern Europe [11, 12]. Fusarium verticillioides can survive in plant residue, healthy seeds and soil, and initiate the infection of maize from seedborne or airborne inoculum, causing seedling disease, Fusarium stalk rot and FER [10, 13]. FER usually occurs on physically injured kernels, random kernels, or groups of kernels, and consists of a light pink or white mold [10]. Infected maize kernels contain toxic fumonisin that is carcinogenic in humans and livestock and even causes porcine pulmonary edema, equine leukoencephalomalacia and rat hepatocarcinoma [14, 15]. Chemical and agronomic measures are not very effective in controlling FER [1]. The best strategy to control FER and to reduce the incidence of fumonisin contamination is breeding and promoting maize varieties with genetic resistance in order [16]. Moreover, in a RIL population from NC300 × B104, Robertson et al. [17] found the genomic and phenotypic correlations between FER and fumonisin accumulation is 0.87 and 0.64, respectively, indicating that it was possible to select lines with reduced FER and fumonisin contamination at the same time [17]. These strategies require us to understand the genetics of resistance clearly, and identify the alleles that can significantly reduce the hazard from F. verticillioides [16].

Resistance to FER is complex because it is characterized by a quantitative inheritance in which additive, dominant, epistatic, and genotype by environment interaction effects are important [18,19,20,21]. Based on biparental populations, Mapping studies have shown that resistance to FER is controlled by many genes with relatively small effects that vary in environments and populations [4, 22]. Although different maize inbred lines and hybrids own different genetic variation for resistance to FER, there is no evidence of maize materials with complete resistance to either FER or fumonisin contamination in maize [23,24,25]. It is very important to identify novel resistance genes against F. verticillioides in order to find a lasting solution to FER problems in maize production. Several studies have identified QTL associated with resistance to F. verticillioides and subsequent reduced fumonisin accumulation using cross-populations, such as F2, F2:3, RILs. Zhang et al. detected six and four QTL in a F2 population of 230 individuals in two environments, respectively, and two QTL were identified consistently in both environments [26]. Using a F2:3 population, Pérez-Brito et al. [21] detected 13 QTL for kernel resistance to FER, which displayed significant QTL × environment interactions, and Chen et al. [19] discovered a QTL for FER resistance affecting approximately 18% of the phenotypic variation and accounting for up to 35% of the phenotypic effect in near isogenic lines when in homozygosity. In two additional studies based on RIL populations, Ding et al. [20] detected two QTL on chromosome 3 (bin3.04), which were consistently identified across all environments, and found significant epistatic effects among QTL and interactions effects between mapped loci and environments, and Li et al. [27] detected a resistance QTL with 10.2% of the phenotypic variation, but no epistatic effects were detected. In addition, complexity of FER could be associated with grain moisture content (GM) and European corn borer (Ostrinia nubilalis) [28, 29].

Recently, to uncover genomic regions associated with reduced FER and fumonisin B1 (FB1) mycotoxin contamination and identify molecular markers to perform marker-assisted selection, Maschietto et al. [30] used an F2:3 population of 188 progenies developed by crossing CO441 (resistant) and CO354 (susceptible) genotypes and evaluated FER severity and FB1 contamination content, and detected 15 QTL for FER and 17 QTL for FB1 contamination. Eight QTL located on chromosomes 1, 2, 3, 6, 7, and 9 were in common between FER and FB1, making the selection of genotypes possible with resistance to FER and low fumonisin contamination [30]. Certainly, there are many other studies on resistance to FER based on linkage mapping. This approach is widely used because linkage mapping generates lower false positive results which make up the defect of few alleles in offspring populations [31,32,33]. However, no genes have been isolated by map-based cloning to date, and few stable QTL have been verified for molecular breeding.

GWAS has shown enormous potential for detecting QTL with high resolution in diverse germplasm [34]. A large number of recombinational events and tens of thousands of SNPs increase the accuracy and shorten the confidence interval of QTL mapping. Now many quantitative traits have been successfully studied by GWAS in maize [35]. In 2016, Chen and his colleagues presented 45 SNPs that were significantly related to FER resistance, each of which had relatively small additive effects on FER resistance and explained 1–4% of phenotypic variation [36]. In addition, works on GWAS for FER have been performed by many other research teams, such as [37,38,39,40], and so on. Compared to traditional linkage analysis, association mapping offers higher mapping resolution and eliminates the time and cost of developing synthetic mapping populations, which make up the defect of false positive [41, 42]. So combining GWAS and linkage mapping could play a great role in identifying casual loci [43, 44].

In this study, we reused linkage mapping to identify genomic regions associated with FER resistance in a biparental RIL population that was evaluated across four environments. Then, GWAS was performed based on the data collected from five environments to detect alleles associated with resistance to FER. Next, we validated four common genomic regions in NIL populations and analyzed the candidate genes in these regions. Last, we discussed the probable mechanism of resistance to FER and stable QTL for molecular breeding.

Results

Phenotypic analysis

First of all, we determined the best time of inoculation for ear rot. For determining the proper inoculation time, we evaluated the phenotype of six inoculation periods of the resistant materials, BT-1 and CML295, and susceptible N6 from the 5th to the 35th day after silking (Fig. S1). The resistant materials BT-1 and CML295 showed susceptibility from the 5th to 10th day after silking, but were stable and resistant from the 15th to 35th day. The FER resistance of susceptible N6 became more and more resistant from the 5th to 35th day after silking. However, the most significant difference in resistance between N6 and BT-1 or CML295 was from the inoculation on 15th day after silking; thus, it was effortless to evaluate the materials inoculated at this time.

Descriptive statistics for FER resistance in the RIL and GWAS populations are presented in Table 1. There was a visible difference in resistance between parent lines BT-1 and N6, which had combined means 1.10 and 6.11, respectively (Fig. S2). The wide variations were also observed in each environment in the RIL and GWAS population, which ranged from 1 to 7. The frequency of phenotypic value of the GWAS population for resistance followed an approximately normal distribution, but a skewed distribution in the RIL (Fig. S3). The genotypic variance (σ2g) and genotype-by-environment variance (σ2ge) of resistance were significant in both populations. Heritability for resistance was 0.81 in the RIL population, 0.79 in the GWAS population. The high heritability indicated that much of the phenotypic variance was genetically controlled in the populations and suitable for QTL mapping.

Table 1 Descriptive statistics of FER resistance for the RIL and GWAS populations in different environments

QTL mapping analysis

A total of 10 QTL were identified for FER resistance (Table 2, Fig. S4), which were located on Chr. 1 (bin 1.02/03), Chr. 2 (bin 2.00/01), Chr. 3 (bin 3.01/02, 3.06/07), Chr. 4 (bin 4.05, 4.05/06, 4.08), Chr. 5 (bin 5.00, 5.03/04), and Chr. 10 (bin 10.6/07). The increasing resistance effect of eight QTL originated from the resistant parent BT-1, whereas only two QTL from the susceptible parent N6. Among these QTL, three QTL were located on chromosome 4 (bin 4.05/08) and the one WQ6 (bin 4.05/06) between markers mmc0371 and A007339 had the largest resistance effect for Fusarium ear rot, which could explain more than 9.3% of the phenotypic variation. Then the QTL, on bin 3.06/07 had the second largest resistance effect explaining about 4.5%. These 10 QTL showed both additive effects (A) and additive by environment effects (AbyE), but additive effects explained 25.1% of the phenotypic variation, whereas interaction effects explained only 5.5%.

Table 2 Quantitative trait loci (QTL) for FER resistance identified in the RIL population using the ICIM-ADD method under MET

To determine the epistatic effect, epistatic QTL mapping was performed. A total of three pairs of QTL interactions were detected by the ICIM-EPI method at an LOD value of 7, which explained 3.2, 2.4, and 2.4% of the phenotypic variation (Table S1, Fig. S5). The epistatic effect between QTL with flanking markers umc2256 and bnlg1144 and QTL with umc1791 and IDP4548 had the largest effect, and explained 3.2%. Although each QTL had the negative additive effect, the interaction effect showed a positive effect, which revolved the complexity of the resistance to FER.

GWAS for FER

Single marker-based GWAS was performed using a mixed linear model (MLM) incorporating both the population structure (first three PCs) and K into the model. A total of 18 SNPs were significantly associated with FER resistance with p ≤ 1.0 × 10− 4 (Table 3, Fig. 1). These SNPs explaining 5.6 -10.2% of phenotypic variation was distributed on five chromosomes, and the number of SNPs per chromosome ranged from 1 on chromosome 3 to 6 on chromosomes 4. The most significant SNP was located on chromosome 7 (S7_153,838,246) with the lowest P value (3.38 × 10− 6) and it explained 10.2% of the phenotypic variation. The second SNP with the lowest P value was located on chromosome 4 and explained 6.8% of the phenotypic variation. Detailed information of 18 SNPs significantly associated with FER resistance is provided in Table 3. The QQ Plot showed that the observed P value was in agreement with the expected P value, whereas the observed P value was lower than the expected P value at a threshold greater than three (Fig. S6). As a result, FER may not be explained by a major gene. Some loci with lower significance may not have been detected, but this should not have affected the identification of loci significantly associated with FER resistance. Based on the physical position of the significant SNPs in the B73 version 2 reference genome, these significant SNPs were associated with 11 candidate genes, some of which were directly involved in resistance according to gene annotation, GRMZM2G150179, for example.

Table 3 The significant single nucleotide polymorphisms (SNPs) and their candidate genes associated with FER resistance identified in this study
Fig. 1
figure 1

Manhattan plots of GWAS for the F. verticillioides ear rot resistance in maize

Gene Ontology (GO) annotation was carried out on 11 candidate genes identified by GWAS. The process of growth, stress response, and cell formation was significantly enriched. These processes feel into four main categories, including seven associated candidate genes. The first was the cytoskeleton process, including cytoskeleton and cellular component organization, and involved candidate genes GRMZM2G449160 and GRMZM2G463471. The second was the process of protein binding, which involved the most genes, including GRMZM2G107686, GRMZM2G086072, GRMZM2G463471, GRMZM2G134980, and GRMZM5G818643, which indicated the significance in FER in posttranslational regulation. The third category was the process of regulation of cellular processes, and contained GRMZM2G107686, GRMZM2G086072, and GRMZM2G449160. The last category was the process of stress response, involving GRMZM2G059381 and GRMZM2G134980.

Conjoint analysis for FER resistance

Ten QTL identified through linkage mapping and 18 significant single SNPs detected by GWAS were integrated to analyze the resistance, and four consistent loci were found (Table 4), located on bin3.01/3.02 (WQ4), bin4.05/4.06 (WQ5, WQ6), bin4.08 (WQ7), and bin5.00 (WQ8). These SNPs were further studied in the following experiment. From the conjoint analysis, identification of consistent loci suggested that there were resistance loci for FER with stable effects at different genetic backgrounds and environmental conditions.

Table 4 The consistent loci from linkage mapping and GWAS

QTL verification

To fine map the QTL (WQ5, WQ6, and WQ7) on chromosome 4, a NIL population with the genetic background of susceptible parent N6 was developed using a backcross and marker assistance selection with flanking markers. The percentage of infected kernels (PIK) was brought into the phenotype evaluation. The lines N-44 and N-54 with positive homozygous alleles (WQ5, WQ6, and WQ7) from the resistant parent BT-1 showed a lower PIK compared with N-55 and N6, and N-55 with only WQ5 and WQ6 was more resistant than parent N6, but more susceptible than lines N-44 and N-54, regardless of Zhengzhou and Xuchang. This indicated that WQ7 could decrease approximately 8 and 6 PIK. WQ5 and WQ6 together improved approximately 7 and 8% in resistance compared with N6 in Zhengzhou and Xuchang, respectively (Table 5). The analysis of variance also showed the same result, which indicated that these three QTL could increase resistance to FER. The detailed genotypes and phenotypes of the three NILs can been found in the supplementary materials (Table S2, Fig. 2).

Table 5 The genotype and phenotype of NILs in two environments
Fig. 2
figure 2

Phenotypic variation in FER severity at harvest among the NILs in artificially inoculated ears with F. verticillioides. N-44 is represented by the two ears on the left (a), N-54 (b), N-55 (c) in the middle, and N6 on the right (d)

A segregation population was constructed for WQ3 by a backcross between the NIL, CP-1 with the target the fragment linked with umc2101 and umc2256, and recurrent parent N6. Finally, WQ3 was verified by a family with a total of 58 plants in 2017 (Table S3).

Furthermore, the GWAS also showed a total of four significant SNPs with p < 1 × 10− 4 in these verified QTL and represented three candidate genes: GRMZM2G449160 for WQ3, GRMZM2G463471 for WQ5 and WQ6, and GRMZM2G059381 for WQ7 (Table 4).

Discussion

QTL analysis and GWAS for FER resistance

QTL analysis is a well-established and widely-used tool for dissecting the genetic basis of complex traits in plants [45]. Many QTL associated with important agronomic traits have been mapped but only a few causal genes were cloned in cereals [46, 47]. Similarly, to date, many QTL have been mapped, but no causal genes have been cloned underlying QTL for FER resistance controlled by many minor-effect QTL that play a great role in maize [48]. These indicate that the positional cloning of minor-effect QTL is still difficult because of their low heritability. Compared to traditional linkage-based analyses, GWAS offers higher mapping effects containing mapping resolution and a greater number of loci, because of many polymorphic SNPs, and eliminates the time and cost associated with developing synthetic mapping populations [41, 42]. However, GWAS easily generates false positive results because of the population structure. Thus, combining GWAS and linkage mapping could exploit the complementary strengths of both approaches to identify casual loci or genes [43, 44].

To decrease the loss from FER and explore the genetic mechanism, we begin to study resistance to FER more than a decade years ago. Today, we have formed a series of relatively perfect inoculating systems and phenotypic identification methods [49], and have achieved some degree of success [19, 27, 36, 50,51,52]. In this study, 10 QTL and 18 SNPs (P < 1 × 10− 4) were detected on the whole genome. Among them, four significant SNPs were located in four QTL, which represented three candidate genes, GRMZM2G449160, GRMZM2G463471, and GRMZM2G059381. GRMZM2G449160 is a member of glutaredoxins (GRXs), which belongs to the antioxidants involved in cellular stress responses. Proteomic analysis found that homologous OsGRX20 increased by 2.7-fold after infection by bacterial blight in rice [53]. GRMZM2G059381 belongs to the AMP-binding protein and the homologous OsBIABP1 is involved in the regulation of the defense response through salicylic acid (SA) and/or jasmonic acid (JA) / ethylene (ET) signaling pathways [54]. GRMZM2G463471 is a member of the actin-depolymerizing factors (ADFs), whose primarily function is to regulate the severing and depolymerization of actin filaments. However, in recent years, the activity of ADFs proteins has been linked to a variety of cellular processes, including those associated with responses to stress. Zhang et al. [55] found a member of ADFs, e.g., TaADF4, from wheat, was required for resistance to the stripe rust pathogen Puccinia striiformis f. sp. Tritici. These results indicate that the three candidate genes in this study may be associated with FER resistance in maize, which will be focused on in the following study.

Phenotypic evaluation for FER resistance

An accurate phenotype is the key to the study of FER. The acquisition of the phenotype was influenced by the inoculation method, date, and the inoculation dose. At present, there are three common inoculation methods used for the study of FER resistance, namely the silk channel inoculation method [56, 57], silk sprayed with inoculum method [4], and the sponge and nail punch method [58]. Among them, the last method is widely used because of easy control of the inoculation dose.

In the long-term study of FER, we explored and optimized the nail punch method [49]. The key to this approach is the operation timing of inoculation. This method is most suitable for inoculation in the milk ripening period, the 15th day after silking, because earlier or later inoculation can not accurately reflect the resistance of the materials. The most significant difference in resistance between susceptible inbred line N6 and resistant inbred line BT-1 or CML295 was from the inoculation on 15th day after silking; thus, it was effortless to evaluate the materials inoculated at this time. To assess the resistance of polymorphic GWAS population, it was divided into two parts according to the date after silking and planted at two different times to ensure the same time of inoculation.

Stable QTL for Fusarium resistance in different tissues and studies

For more than 10 years, our group studied Fusarium resistance in different maize tissues [19, 27, 36, 50,51,52]. We confirmed that the resistance loci and mechanism of different tissues were different. In the GWAS population, we found some lines showed different resistance between different tissues, for example some lines had high FER resistance with weak Fusarium cob rot resistance (FCR) and Fusarium seed rot resistance (FSR). Therefore, we compared the QTL identified for Fusarium resistance in ear, cob (FCR), and seed (FSR) (Fig. 3a). These studies used the same resistance parent line and similar susceptible lines, but had different results [50, 52]. Some QTL were independent, for example, the QTL located on bin 3.01/02, bin 5.00, and bin 10.06/07 were specific for FER. However, there were still some QTL that appeared conservative; for example, QTL in bin 3.06/3.07 and 4.05/4.06 were identified in FER, FCR, and FSR. These conservative QTL were not identified for Fusarium seeding rot, which may cause by different resistant mechanism or different research populations [59]. The candidate genes identified by GWAS were different between these three traits, as well [50, 52, 60]. It is worth mentioning that many candidate genes for FCR were involved in cell and tissue development, whereas none of this kind of gene was identified in FER. These results suggested that the resistance against Fusarium in maize is a complex trait, which requires continued research to resolve this important production problem.

Fig. 3
figure 3

QTL and SNPs for Fusarium resistance in different tissues and studies. a) Comparison of QTL for Fusarium rot resistance detected by different tissues. b) Comparison of QTL and SNPs for Fusarium ear rot resistance detected by previous reports. Colored lines represent different QTL or SNPs in different studies

In our previous study [27], four QTL located in bin3.01/4.06/5.01/6.04 were detected using composite interval mapping method across two environments (2007ZZ and 2008ZZ). Among these QTL, two QTL in bin3.01/5.01 were also detected in this study, but the QTL in bin4.06 in the previous study was not detected; rather bin4.05/4.06 and bin4.08 in this study were verified by NILs. Moreover, no significant loci were detected on chromosome 6 in the RIL population and GWAS population in this analysis.

Previous studies [19, 30, 36, 37, 61, 62] together with the present study were compared, which showed that loci for FER resistance are widely distributed, indicating FER resistance is a complex trait controlled by many minor QTL (Fig. 3b). The five common loci detected in RIL and GWAS populations were also identified by other researchers showing that these QTL were stable. These results also indicated that maybe it was effective to increase resistance according to selecting lines with a few stable QTL. Furthermore, the QTL verified by NILs on bin 3.06/07 and 4.05/06 will be the focus on our following study plan, because they were identified by different research and Fusarium resistance traits.

Conclusions

In this study, a RIL populations and one GWAS population were used to identify and map QTL for FER resistance. A total of 10 QTL for FER resistance were detected by QTL mapping and 18 SNPs were identified by GWAS at P < 1 × 10− 4. Four QTL of the five common loci in the RIL and GWAS were verified in NILs and three candidate genes may be associated with FER resistance. These results confirmed that FER resistance was strongly controlled by multiple genes with low effect and the QTL and candidate genes identified in this study could help to better understand the genetic basis and explore the mechanism of FER resistance. At the same time, it was feasible to select maize lines with higher Fusarium resistance because of the two stable QTL in different tissues and studies.

Methods

Germplasm and experimental design

A biparental population composed of 250 recombinant inbred lines (RILs) was constructed by a cross between inbred lines BT-1 and N6 (Fig. S2) by the single-seed descent method. The line BT-1 was reformed by tropical Asia material with highest resistance to FER, which was screened out of 90 inbred lines inoculated with Fusarium verticillioides [63], whereas the susceptible N6 was a Tangsipingtou line. The RIL population together with the two parents were grown in Zhengzhou (34°52′N 113°37′E) in 2007 and 2008 (2007ZZ, 2008ZZ), Wenxian (34°57′N 113°2′E) in 2015 (2015WX), and Xuchang (34°3′N 113°41′E) in 2016 (2016XC). These three areas have a temperate continental climate, with adequate rainfall and enough light-temperature resource in the growing period of maize. NILs were developed by crossing a recipient line N6 with the donor line BT-1 through five cycles of advanced backcrosses. These NILs were evaluated in 2016 and 2017.

The GWAS population consisting of 265 maize inbred lines from the heterotic populations Tangsipingtou, Reid, Lancaster, P group, and some tropical lines from the International Maize and Wheat Improve Center (CIMMYT), was evaluated in Zhengzhou in 2014, 2015, and 2016 (2014ZZ, 2015ZZ, and 2016ZZ) and Wenxian in 2015 (2015WX) and Xuchang in 2016 (2016XC).

All the populations above were laid out in a randomized complete block design with two replications in each environment. Sixteen plants were planted in 4 m row plots with 0.67 m row spacing. Importantly, pesticide was artificially sprayed at the ten-leaf stage to control corn borers. Some other field management was performed according to the standard agronomic practices in each location.

Inoculation and FER evaluation

A single-spore isolate of F. verticillioides was reproduced artificially on sterile mature maize kernels, incubated for 7 d at 28 °C. Then, the spores were harvested, and the concentration was estimated using a hemocytometer and adjusted to 1 × 106 spores ml− 1 in sterile distilled water with 0.2 μl/ml Tween 80 surfactant. The top ear of each plant in each row was inoculated for 2 ml spore suspension on the fifteenth day after silking using the sponge and the nail punch method [58] along with the portable inoculating tool assembled by ourselves.

Resistance to FER was assessed by disease severity according to Reid et al. [64] using a rating scale from 1 to 7 [64]. Ten ears of each line in a row were selected for phenotypic evaluation and the final data for each replication was the average value. The average rating scale from the two replications represented the final value in each environment.

Phenotypic data analysis

The analysis of variance (ANOVA) of phenotype data was carried out using the multifunctional IciMapping version 4.2 [65]. Best linear unbiased predictions (BLUPs) of the combined resistance to FER for each line in Linkage analysis and GWAS were calculated by using a mixed linear model (lmer) in the R version 3.6.0 with the R stats package [66], in which used replication, environments (year-location combinations), and entries were considered random effects. The BLUP value of each line was used for GWAS analysis.

Linkage mapping

In our previous study, we constructed a linkage map of the RIL population (BT × N6), which contained 207 polymorphic SSR markers and had a length of 1820.8 cM with an average 11.7 cM distance [27]. In this study, we reconstructed a linkage map with 222 polymorphic markers using the IciMapping version 4.2, which had a total length of 1890.7 cM and an average genetic distance of 8.5 cM between markers.

To understand the QTL by environment interaction effects, the mapping strategy of MET (mapping of additive and digenic epistasis genes from multi-environmental trials) was performed in the IciMapping software. Two methods were used in QTL mapping, i.e., (i) ICIM-ADD: Inclusive Composite Interval Mapping of additive QTL (ii) ICIM-EPI: Inclusive Composite Interval Mapping of digenic epistatic QTL. The threshold value of LOD was 2.5 for ICIM-ADD and 7 for ICIM-EPI.

GWAS and candidate gene annotation

Through the genotyping-by-sequencing (GBS) method conducted in Cornell University, a total of 955,650 SNPs were identified in the GWAS population [67, 68], and the SNP flanking sequence and position information were available on the “panzea” website (http://cbsusrv04.tc.cornell.edu/users/panzea/download.aspx?filegroupid=4). The filtered parameters of SNPs, linkage disequilibrium (LD) between each pair of SNPs, the principal component analysis (PCA), Kinship matrix and population structure analysis were performed according to our previous study [69].

Candidate gene information was obtained from the MaizeGDB (http://www.maizegdb.org/) genome browser based on the physical position of significant SNPs in B73 RefGen_v2.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

Abbreviations

FER:

Fusarium ear rot

RIL:

Recombinant inbred lines

QTL:

Quantitative trait loci

SNPs:

Single nucleotide polymorphisms

GWAS:

Genome-wide association study

NILs:

Near-isogenic lines

ADFs:

Actin-depolymerizing factors

FCR:

Fusarium cob rot

FSR:

Fusarium seed rot

GO:

Gene Ontology

PIK:

Percentage of infected kernels

References

  1. Munkvold GP, Desjardins AE. Fumonisins in maize: can we reduce their occurrence. Plant Dis. 1997;81:556–65.

    Article  CAS  PubMed  Google Scholar 

  2. Missmer S, Suarez L, Felkner M. Exposure to fumonisins and the occurrence of neural tube defects along the Texas-Mexico border. Environ Health Perspect. 2006;114:237–41.

    Article  PubMed  Google Scholar 

  3. Kebebe AZ, Reid LM, Zhu X, Wu J, Woldemariam T, Voloaca C, Xiang K. Relationship between kernel drydown rate and resistance to Gibberella ear rot in maize. Euphytica. 2014;201:79–88.

    Article  CAS  Google Scholar 

  4. Mesterházy Á, Lemmens M, Reid LM. Breeding for resistance to ear rots caused by Fusarium spp. in maize – a review. Plant Breed. 2012;131:1–19.

    Article  Google Scholar 

  5. Seifert KA, Aoki T, Baayen RP, Brayford D, Burgess LW, Chulze S, Gams W, Geiser D, Gruyter JD, Leslie JF, Logrieco A, Marasas WFO, Nirenberg HI, Donnell KO, Rheeder J, Samuels GJ, Summerell B, Thrane U, Waalwijk C. The name Fusarium moniliforme should no longer be used. Mycol Res. 2003;107:643–4.

    Article  Google Scholar 

  6. Ridenour JB, Bluhm BH. The novel fungal-specific gene FUG1 has a role in pathogenicity and fumonisin biosynthesis in Fusarium verticillioides. Mol Plant Pathol. 2017;18(4):513–28.

    Article  CAS  PubMed  Google Scholar 

  7. Duan C, Qin Z, Yang Z, et al. Identification of Pathogenic Fusarium spp. Causing Maize Ear Rot and Potential Mycotoxin Production in China. Toxins (Basel). 2016;8(6):186.

    Article  CAS  Google Scholar 

  8. Fu M, Li R, Guo C, Pang M, Liu Y, Dong J. Natural incidence of Fusarium species and fumonisins B1 and B2 associated with maize kernels from nine provinces in China in 2012. Food Addit Contam Part A Chem Anal Control Expo Risk Assess. 2015;32(4):503–11.

    Article  CAS  PubMed  Google Scholar 

  9. Zhang XF, Zou CJ, Cui LN, Xiao LI, Yang XR, Luo HH. Identification of pathogen causing maize ear rot and inoculation technique in Southwest China. Southwest China Journal of Agricultural Sciences. 2012;25:2078–82.

    Google Scholar 

  10. Munkvold GP. Epidemiology of Fusarium diseases and their mycotoxins in maize ears. Eur J Plant Pathol. 2003;109:705–13.

    Article  CAS  Google Scholar 

  11. Logrieco A, Bottalico A, Mulè G, Moretti A, Perrone G, Mulé G, Moretti A, Perrone G. Epidemiology of toxigenic fungi and their associated mycotoxins for some Mediterranean crops. Eur J Plant Pathol. 2003;109:645–67.

    Article  CAS  Google Scholar 

  12. Santiago R, Cao A, Butrón A. Genetic factors involved in fumonisin accumulation in maize kernels and their implications in maize agronomic management and breeding. Toxins. 2015;7:3267–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Irma MR, De YMJ, Silva-Rojas HV, Gabino GDLS, Guzmán-De-Peña DA. Biodiversity of Fusarium species in Mexico associated with ear rot in maize, and their identification using a phylogenetic approach. Mycopathologia. 2007;163(1):31–9.

    Article  Google Scholar 

  14. Covarelli L, Stifano S, Beccari G, Raggi L, Lattanzio VM, Albertini E. Characterization of Fusarium verticillioides strains isolated from maize in Italy: Fumonisin production, pathogenicity and genetic variability. Food Microbiol. 2012;31:17–24.

    Article  CAS  PubMed  Google Scholar 

  15. Gelderblom WC, Jaskiewicz K, Marasas WF, Thiel PG, Horak RM, Vleggaar R, Kriek NP. Fumonisins-novel mycotoxins with cancer-promoting activity produced by Fusarium moniliforme. Appl Environ Microbiol. 1998;54:1806–11.

    Article  Google Scholar 

  16. Mukanga M, Derera J, Tongoona P. Gene action and reciprocal effects for ear rot resistance in crosses derived from five tropical maize populations. Euphytica. 2010;174:293–301.

    Article  CAS  Google Scholar 

  17. Robertson LA, Jines MP, Balint-Kurti P, Kleinschmidt CE, White DG, Payne G, Maragos CM, Molnár TL, Holland JB. QTL mapping for Fusarium ear rot and Fumonisin contamination resistance in two maize populations. Crop Sci. 2006;46:1734–44.

    Article  CAS  Google Scholar 

  18. Boling MB, Grogan CO. Gent action affecting host resistance to Fusarium ear rot of maize. Crop Sci. 1965;4:305–7.

    Google Scholar 

  19. Chen J, Ding J, Li H, Li Z, Sun X, Li J, Wang R, Dai X, Dong H, Song W, Chen W, Xia Z, Wu J. Detection and verification of quantitative trait loci for resistance to Fusarium ear rot in maize. Mol Breed. 2012;30:1649–56.

    Article  CAS  Google Scholar 

  20. Ding JQ, Wang XM, Subhash C, Yan JB, Li JS. QTL mapping of resistance to Fusarium ear rot using a ril population in maize. Mol Breed. 2008;22:395–403.

    Article  Google Scholar 

  21. Pérez-Brito D, Jeffers D, González-de-León D, Khairallah M, CortésCruz M, Velázquez-Cardelas M, Azpíroz-Rivero S, Srinivasan G. QTL mapping of Fusarium moniliforme ear rot resistance in highland maize, Mexico. Agrociencia. 2001;35:181–96.

    Google Scholar 

  22. Lanubile A, Maschietto V, Borrelli VM, Stagnati L, Logrieco AF, Marocco A. Molecular basis of resistance to Fusarium ear rot in maize. Front Plant Sci. 2017;8:1774.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Clements MJ, Maragos CM, Pataky JK, White DG. Sources of resistance to fumonisin accumulation in grain and Fusarium ear and kernel rot of corn. Phytopathology. 2004;94:251–60.

    Article  CAS  PubMed  Google Scholar 

  24. Clements M, Kleinschmidt C. Evaluation of inoculation techniques for Fusarium ear rot and fumonisin contamination of corn. Plant Dis. 2003;87:147–53.

    Article  CAS  PubMed  Google Scholar 

  25. Lanubile A, Pasini L, Lo Pinto M, Battilani P, Prandini A, Marocco A. Evaluation of broad spectrum sources of resistance to Fusarium verticillioides and advanced maize breeding lines. World Mycotoxin J. 2011;1:43–51.

    Article  CAS  Google Scholar 

  26. Zhang F, Wan XQ, Pan GT. Molecular mapping of QTL for resistance to maize ear rot caused by Fusarium moniliforme. Acta Agron Sin. 2007:491–6.

  27. Li ZM, Ding JQ, Wang RX, Chen JF, Sun XD, Chen W, Song WB, Dong HF, Dai XD, Xia ZL, Wu JY. A new QTL for resistance to Fusarium ear rot in maize. J Appl Genet. 2011;52:403–6.

    Article  PubMed  Google Scholar 

  28. Blandino M, Scarpino V, Vanara F, Sulyok M, Krska R, Reyneri A. Role of the European corn borer (Ostrinia nubilalis) on contamination of maize with 13 Fusarium mycotoxins. Food Addit Contam Part A Chem Anal Control Expo Risk Assess. 2015;32(4):533–43.

    Article  CAS  PubMed  Google Scholar 

  29. Xiang K, Reid LM, Zhang ZM, Zhu XY, Pan GT. Characterization of correlation between grain moisture and ear rot resistance in maize by QTL meta-analysis. Euphytica. 2012;183(2):185–95.

    Article  CAS  Google Scholar 

  30. 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.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  31. Ding J, Zhang L, Chen J, Li X, Li Y, Cheng H, Huang R, Zhou B, Li Z, Wang J, Wu J. Genomic dissection of leaf angle in maize (Zea mays L.) using a four-way cross mapping population. PLoS One. 2015;10:e0141619.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  32. Jiang C, Zeng ZB. Multiple trait analysis of genetic mapping for quantitative trait loci. Genetics. 1995;140(3):1111–27.

    CAS  PubMed  PubMed Central  Google Scholar 

  33. Yang Q, Zhang D, Xu M. A sequential quantitative trait locus fine-mapping strategy using recombinant-derived progeny. J Integr Plant Biol. 2012;54:228–37.

    Article  PubMed  Google Scholar 

  34. Buntjer JB, Sørensen AP, Peleman JD. Haplotype diversity: the link between statistical and biological association. Trends Plant Sci. 2005;10:466–71.

    Article  CAS  PubMed  Google Scholar 

  35. Xiao Y, Liu H, Wu L, Warburton M, Yan J. Genome-wide association studies in maize: praise and stargaze. Mol Plant. 2017;10:359.

    Article  CAS  PubMed  Google Scholar 

  36. 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. Genes Genom Genet. 2016;6:3803–15.

    CAS  Google Scholar 

  37. Butron A, Santiago R, Cao A, Samayoa LF, Malvar RA. QTL for resistance to Fusarium ear rot in a multiparent advanced generation intercross (MAGIC) maize population. Plant Dis. 2019;103(5):897–904.

    Article  CAS  PubMed  Google Scholar 

  38. Coan MMD, Senhorinho HJC, Pinto RJB, Scapim CA, Tessmann DJ, Paul WW, Warburton ML. Genome-wide association study of resistance to ear rot by, in a tropical field maize and popcorn core collection. Crop Sci. 2018;58(2):564.

    Article  CAS  Google Scholar 

  39. Zila CT, Ogut F, Romay MC, Gardner CA, Buckler ES, Holland JB. Genome-wide association study of Fusarium ear rot disease in the U.S.A. maize inbred line collection. BMC Plant Biol. 2014;14(1):372.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Zila CT, Samayoa LF, Santiago R, Butrón A, Holland JB. A genome-wide association study reveals genes associated with Fusarium ear rot resistance in a maize core diversity panel. G3 (Bethesda). 2013;3(11):2095–104.

    Article  CAS  Google Scholar 

  41. Flint-Garcia SA, Thuillet AC, Yu J, Pressoir G, Romero SM, Mitchell SE, Doebley J, Kresovich S, Goodman MM, Buckler ES. Maize association population: a high-resolution platform for quantitative trait locus dissection. Plant J. 2005;44:1054–64.

    Article  CAS  PubMed  Google Scholar 

  42. Yu J, Buckler EES, Mamo BE, Barber BL, Steffenson BJ. Genetic association mapping and genome organization of maize. Plant Biotechnology. 2006;17:155–60.

    CAS  Google Scholar 

  43. Motte H, Vercauteren A, Depuydt S, Landschoot S, Geelen D, Werbrouck S, Goormachtig S, Vuylsteke M, Vereecke D. Combining linkage and association mapping identifies RECEPTOR-LIKE PROTEIN KINASE1 as an essential Arabidopsis shoot regeneration gene. Proc Natl Acad Sci U S A. 2014;111:8305–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  44. Pedergnana V, Syx L, Cobat A, Guergnon J, Brice P, Fermé C, Carde P, Hermine O, Le-Pendeven C, Amiel C, Taoufik Y, Alcais A, Theodorou I, Besson C, Abel L. Combined linkage and association studies show that HLA class II variants control levels of antibodies against Epstein-Barr virus antigens. PLoS One. 2014;9:e102501.

    Article  PubMed  PubMed Central  Google Scholar 

  45. Cooper M, van Eeuwijk FA, Hammer GL, Podlich DW, Messina C. Modeling QTL for complex traits: detection and context for plant breeding. Curr Opin Plant Biol. 2009;12:231–40.

    Article  CAS  PubMed  Google Scholar 

  46. Holland JB. Genetic architecture of complex traits in plants. Curr Opin Plant Biol. 2007;10(2):156–61.

    Article  CAS  PubMed  Google Scholar 

  47. Mackay TFC, Stone EA, Ayroles JF. The genetics of quantitative traits: challenges and prospects. Nat Rev Genet. 2009;10:565–77.

    Article  CAS  PubMed  Google Scholar 

  48. Brauner PC, Melchinger AE, Schrag TA, Utz HF, Schipprack W, Kessel B, Ouzunova M, Miedaner T. Low validation rate of quantitative trait loci for gibberella ear rot resistance in European maize. Theor Appl Genet. 2017;130:175–86.

    Article  CAS  PubMed  Google Scholar 

  49. Dong C, Wu Y, Gao J, Zhou Z, Mu C, Ma P, Chen J, Wu J. Field inoculation and classification of maize ear rot caused by Fusarium verticillioides. Bio-protocol. 2018;8:e3099.

    CAS  Google Scholar 

  50. Ju M, Zhou Z, Mu C, Zhang X, Gao J, Liang Y, Chen J, Wu Y, Li X, Wang S, Wen J, Yang L, Wu J. Dissecting the genetic architecture of Fusarium verticillioides seed rot resistance in maize by combining QTL mapping and genome-wide association analysis. Sci Rep. 2017;7:46446.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Wang Y, Zhou Z, Gao J, Wu Y, Xia Z, Zhang H, Wu J. The mechanisms of maize resistance to Fusarium verticillioides by comprehensive analysis of RNA-seq data. Front Plant Sci. 2016;7:1654.

    PubMed  PubMed Central  Google Scholar 

  52. Mu C, Gao J, Zhou Z, et al. Genetic analysis of cob resistance to F. verticillioides: another step towards the protection of maize from ear rot. Theor Appl Genet. 2019;132(4):1049–59.

    Article  CAS  PubMed  Google Scholar 

  53. Ning X, Sun Y, Wang C, Zhang W, Sun M, Hu H, Liu J, Yang L. A rice cpyc-type glutaredoxin OsGRX20 in protection against bacterial blight, methyl viologen and salt stresses. Front Plant Sci. 2018;9:111.

    Article  PubMed  PubMed Central  Google Scholar 

  54. Zhang XC, Yu X, Zhang HJ, Song FM. Molecular characterization of a defense-related AMP-binding protein gene, OsBIABP1, from rice. J Zhejiang Univ Sci B. 2009;10(10):731–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Zhang B, Yuan H, Wang J, Yan H, Shimono M, Day B, Ma Q. TaADF4, an actin-depolymerizing factor from wheat, is required for resistance to the stripe rust pathogen Puccinia striiformis f. sp. tritici. Plant J. 2017;89:1210.

    Article  CAS  PubMed  Google Scholar 

  56. Clements MJ, Kleinschmidt CE, Maragos CM, Pataky JK, White DG. Evaluation of inoculation techniques for Fusarium ear rot and fumonisin contamination of corn. Plant Dis. 2003;58 Suppl :C173.

  57. Reid LM, Spaner D, Mather DE, Bolton AT, Hamilton RI. Resistance of maize hybrids and inbreds following silk inoculation with three isolates of Fusarium graminearum. Plant Dis. 1993;77:1248–51.

    Article  Google Scholar 

  58. Drepper WJ, Renfro BL. Comparison of methods for inoculation of ears and stalks of maize with Fusarium moniliforme. Plant Dis. 1990;74:952–6.

    Article  Google Scholar 

  59. Septiani P, Lanubile A, Stagnati L, Busconi M, Nelissen H, Pè ME, Dell’Acqua M, Marocco A. Unravelling the genetic basis of Fusarium seedling rot resistance in the MAGIC maize population: novel targets for breeding. Sci Report. 2019;9:5665.

    Article  CAS  Google Scholar 

  60. Stagnati L, Lanubile A, Samayoa LF, Bragalanti M, Giorni P, Busconi M, Holland JB, Marocco A. A Genome Wide Association Study Reveals Markers and Genes Associated with Resistance to Fusarium verticillioides Infection of Seedlings in a Maize Diversity Panel. G3 (Bethesda). 2019;9(2):571–9.

    CAS  Google Scholar 

  61. Kebede AZ, Woldemariam T, Reid LM, Harris LJ. Quantitative trait loci mapping for Gibberella ear rot resistance and associated agronomic traits using genotyping-by-sequencing in maize. Theor Appl Genet. 2015;129(1):17–29.

    Article  PubMed  CAS  Google Scholar 

  62. Xiang K, Zhang ZM, Reid LM, Zhu XY, Yuan GS, Pan GT. A meta-analysis of QTL associated with ear rot resistance in maize. Maydica. 2010;55(3):281–90.

    Google Scholar 

  63. Chen W, Wu JY, Yuan HX. Identification of resistance on maize germplasm to maize ear rot. Maize Science. 2002.

  64. Reid LM, Mather DE, Bolton AT, Hamilton RI. Evidence for a gene for silk resistance to Fusarium graminearum Schw. Ear rot of maize. J Hered. 1994;85:118–21.

    Article  Google Scholar 

  65. Meng L, Li H, Zhang L, Wang J. QTL IciMapping: integrated software for genetic linkage map construction and quantitative trait locus mapping in bi-parental populations. The Crop Journal. 2015;3:269–83.

    Article  Google Scholar 

  66. R Core Team. A Language and Environment for Statistical Computing. R Foundation for Statistical Computing, Vienna, China, 2019. https://www.R-project.org.

  67. Elshire RJ, Glaubitz JC, Sun Q, Poland JA, Kawamoto K, Buckler ES, Mitchell SE. A robust, simple genotyping-by-sequencing (GBS) approach for high diversity species. PLoS One. 2011;6:e19379.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  68. Glaubitz JC, Casstevens TM, Lu F, Harriman J, Elshire RJ, Sun Q, Buckler ES. TASSEL-GBS: a high capacity genotyping by sequencing analysis pipeline. PLoS One. 2014;9:e90346.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  69. Li X, Zhou Z, Ding J, Wu Y, Zhou B, Wang R, Ma J, Wang S, Zhang X, Xia Z, Chen J, Wu J. Combined linkage and association mapping reveals qtl and candidate genes for plant and ear height in maize. Front Plant Sci. 2016;7:833.

    PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

This work was financially supported by the International (Regional) Joint Research Project of the National Natural Science Foundation of China (Grant No. 31761143009) and the Co-construction State Key Laboratory of Wheat and Crop Science of China.

Funding

This study was supported by the International (Regional) Joint Research Project of the National Natural Science Foundation of China (Grant No. 31761143009) and the Co-construction State Key Laboratory of Wheat and Crop Science of China.

Author information

Authors and Affiliations

Authors

Contributions

JW, YW, and ZZ designed the experiments and wrote the manuscript. YW, CD, JD, HL, YH, RW, XS, JL, XD, WS, and WC performed the QTL identified experiments. YW, CD, ZZ, CM, and XL performed the GWAS experiments and the NIL verify experiments. YW, ZZ, and JC analyzed the data. XZ and JC contributed GWAS population and analysis tools. ZZ analyzed candidate genes and prepared figures. All these authors have read and approved this manuscript.

Corresponding author

Correspondence to Jianyu Wu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that there are no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wu, Y., Zhou, Z., Dong, C. et al. Linkage mapping and genome-wide association study reveals conservative QTL and candidate genes for Fusarium rot resistance in maize. BMC Genomics 21, 357 (2020). https://doi.org/10.1186/s12864-020-6733-7

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-020-6733-7

Keywords