ORIGINAL RESEARCH article
Sec. Plant Breeding
Genomics and Pathways Involved in Maize Resistance to Fusarium Ear Rot and Kernel Contamination With Fumonisins
- 1Misión Biológica de Galicia (CSIC), Pontevedra, Spain
- 2Agrobiología Ambiental, Calidad de Suelos y Plantas (UVIGO), Unidad Asociada a la MBG (CSIC), Pontevedra, Spain
Fusarium verticillioides is a causal agent of maize ear rot and produces fumonisins, which are mycotoxins that are toxic to animals and humans. In this study, quantitative trait loci (QTLs) and bulk-segregant RNA-seq approaches were used to uncover genomic regions and pathways involved in resistance to Fusarium ear rot (FER) and to fumonisin accumulation in maize kernels. Genomic regions at bins 4.07–4.1, 6–6.01, 6.04–6.05, and 8.05–8.08 were related to FER resistance and/or reduced fumonisin levels in kernels. A comparison of transcriptomes between resistant and susceptible inbred bulks 10 days after inoculation with F. verticillioides revealed 364 differentially expressed genes (DEGs). In the resistant inbred bulks, genes involved in sink metabolic processes such as fatty acid and starch biosynthesis were downregulated, as well as those involved in phytosulfokine signaling and many other genes involved in cell division; while genes involved in secondary metabolism and compounds/processes related to resistance were upregulated, especially those related to cell wall biosynthesis/rearrangement and flavonoid biosynthesis. These trends are indicative of a growth–defense trade-off. Among the DEGs, Zm00001d053603, Zm00001d035562, Zm00001d037810, Zm00001d037921, and Zm00001d010840 were polymorphic between resistant and susceptible bulks, were located in the confidence intervals of detected QTLs, and showed large differences in transcript levels between the resistant and susceptible bulks. Thus, they were identified as candidate genes involved in resistance to FER and/or reduced fumonisin accumulation.
Fusarium verticillioides has long been known to be a causal agent of maize ear rot but has received more attention since the discovery that fumonisins, the mycotoxins it produces, can accumulate in maize kernels (Gelderblom et al., 1988). Fumonisin consumption can lead to serious disorders in humans and animals; for example, leukoencephalomalacia in horses and pulmonary edema in swine, both are accompanied by injuries to the liver and heart; hepatic necrosis and, eventually, kidney and liver cancers in rodents; impaired liver function in cattle; and reduced growth in poultry (Munkvold and Desjardins, 1997; Marasas, 2001; Voss et al., 2007). In addition, fumonisins were classified as possibly carcinogenic to humans by the International Agency for Research on Cancer (IARC, 2002). Epidemiological studies have suggested a likely relationship between human consumption of fumonisin-contaminated maize and incidence of esophageal cancer and neural tube defects in human embryos (Rheeder et al., 2002; Gelineau-van Waes et al., 2005; Missmer et al., 2006; Sun et al., 2007; Alizadeh et al., 2012).
Natural infection of maize kernels by F. verticillioides results in fumonisin accumulation in field-grown maize plants. Environmental conditions during kernel development and drying affect the level of fumonisin contamination at harvest; however, maize genetics may also affect both infection by the pathogen and fumonisin accumulation (Cao et al., 2013; Santiago et al., 2015). Ratings for Fusarium ear rot (FER) have extensively been used to estimate kernel infection by F. verticillioides. In previous studies, severity of FER has shown significant and moderate to high correlation coefficients with kernel contamination with fumonisins across different environments and maize materials (Bolduan et al., 2009; Eller et al., 2010; Löffler et al., 2010, 2011; Presello et al., 2011; Santiago et al., 2013). Therefore, most studies focused on reducing kernel contamination with fumonisins through maize resistance management have been devoted to uncovering the genetic basis of resistance to FER (Pérez-Brito et al., 2001; Zila et al., 2014; Chen et al., 2016; Ju et al., 2017). However, F. verticillioides has been defined as a maize endophyte, and fungal infection is not always symptomatic. Even symptom-less plants can accumulate fumonisin mycotoxins (Bacon and Hinton, 1996; Rodríguez Estrada et al., 2012).
As an alternative to laboratory quantification of fumonisins, which is tedious and expensive, molecular markers have been proposed as valuable tools in breeding programs for reducing kernel contamination with these mycotoxins. Moreover, selection based on phenotype alone could be inefficient because of the complexity of maize resistance to fumonisin accumulation (Eller et al., 2008). This complexity is mainly due to the different performance of genotypes across environments and the relative relevance of non-additive genetic effects (Butron et al., 2006, 2015; Robertson et al., 2006; Padua et al., 2016; Maschietto et al., 2017). Nevertheless, although dominance effects could be important in inheritance of resistance to fumonisin accumulation, the results of several studies on the genetics of resistance suggest that it is better to focus on inbred materials than on hybrids. This is because genetic variation is greater among inbred lines than among hybrids, and additive effects are more stable than dominance effects. Moreover, moderate to high genotypic correlation coefficients for resistance to fumonisin accumulation have been reported between inbred lines and test crosses (Löffler et al., 2011; Hung and Holland, 2012; Padua et al., 2016). The search for quantitative trait loci (QTLs) for resistance to fumonisin accumulation has been only performed on few mapping populations, most of which are bi-parental populations (Robertson-Hoyt et al., 2006; Maschietto et al., 2017; Morales et al., 2019; Samayoa et al., 2019; Gesteiro et al., 2021). In this sense, this study intends to gain additional insight into the genetic basis of resistance to kernel infection by exploring genetic variability in a new mapping population derived from a cross between two inbreds with contrasting values for both FER and fumonisin accumulation in kernels.
Other molecular tools besides DNA markers can be used at the genome-wide level to study the genetic basis of resistance to FER and kernel contamination with fumonisins. Gene expression analyses have been conducted to identify genes involved in different responses between resistant and susceptible maize genotypes to kernel infection by F. verticillioides, as described in several reviews (Lanubile et al., 2017; Santiago et al., 2020). These studies have shed light on genes involved in maize kernel, silk, and sheath responses to infection by F. verticillioides. However, there has been limited success in discovering genes involved in resistance, because most studies have used unrelated inbreds with contrasting resistance levels (Lanubile et al., 2010, 2012a,2014; Campos-Bermudez et al., 2013; Yuan et al., 2013; Wang et al., 2016). In studies using such materials, the chance of detecting false positive associations between genes and resistance is increased, because different genetic backgrounds of inbreds could underlie differential gene expression. Resistance/susceptibility is the result of continuous fungus–plant communication.
Most studies have focused on the early response of maize plants to infection. However, there is evidence that the biosynthesis of secondary metabolites that are further transformed into antifungal compounds is preferentially induced in the late stage of F. verticillioides infection (Wang et al., 2016). In addition, Maschietto et al. (2015) demonstrated that plant resistance responses were mounted between 7 and 14-days post-inoculation (dpi). In this study, we used a bulk-segregant RNA-seq (BSR) approach to compare different inbred pools 10 dpi, since it would allow us to simultaneously identify differentially expressed genes (DEGs) and polymorphic DNA markers related to fumonisin accumulation (Du et al., 2017; Lin et al., 2021).
Materials and Methods
Quantitative Trait Loci Mapping for Fusarium Ear Rot and Fumonisin Content
A linkage mapping population of 144 recombinant inbred lines (RILs) derived from a cross between the European flint inbred line EP42 and the American dent inbred line A637 was used for QTL mapping for FER and fumonisin content in kernels. EP42 is susceptible to FER and fumonisin contamination, while A637 is partially resistant to both. RIL genotypes at 121 simple sequence repeat (SSR) markers and the population linkage map have been previously published by Samayoa et al. (2014).
The 144 RILs were evaluated in Pontevedra (42°24′N, 8°38′W, and 20 m above sea level), Spain in 2016 and 2017 using a 12 × 12 lattice design with two replications. Trials were hand-planted, and each experimental plot consisted of one row spaced 0.8 m apart from the other row with 13 one-kernel hills spaced 0.18 m apart. Evaluations were performed under artificial inoculation with F. verticillioides. Five plants of each plot were inoculated in August 12th and 8th in 2016 and 2017, respectively. In each plot, individual ears with still fresh silks were chosen to be inoculated with 2 ml of a spore suspension (106 spores/ml) of a characterized toxigenic isolate of F. verticillioides (Cao et al., 2014). The inoculum was injected into the center of the ear using a four-needle vaccinator that perforated the husks and injured three to four kernels (Reid et al., 1996).
The inoculated ears from each plot were collected 2 months after inoculation and were individually rated for FER using a 7-point scale (1 = no visible disease symptoms, 2 = 1–3%, 3 = 4–10%, 4 = 11–25%, 5 = 26–50%, 6 = 51–75%, and 7 = 76–100% of kernels exhibiting visual symptoms of infection) devised by Reid and Zhu (2005). Then, the inoculated ears from each plot were dried at 35°C for 1 week and shelled, and a kernel sample of approximately 50 g was taken and stored at 4°C until chemical analyses for fumonisin quantification were performed. The kernel samples were ground through a 0.75-mm screen in a Pulverisette 14 rotor mill (Fritsch GmbH, Oberstein, Germany), and total fumonisin (fumonisins B1, B2, and B3) quantification was performed using a commercial ELISA kit (R-Biopharm Rhône Ltd., Glasgow, Scotland, United Kingdom, at the Food Technology Department of the University of Lleida, Spain. Test recovery rate was approximately 60%; assay variation coefficient was approximately 8%; specificities for B1, B2, and B3 were 100%, around 40%, and almost 100%, respectively, and detection limit was 0.025 mg kg–1. Extraction and preparation of samples, as well as test performance, were carried out as described in the commercial kits.
Phenotypic data were analyzed using a mixed model with the PROC MIXED procedure of the SAS statistical package (SAS Institute Inc, 2011). Years, replications, blocks within replications, and years and RILs were considered as random effects. Heritabilities () across environments were estimated for each trait on a family mean basis as described previously by Holland et al. (2003). The genetic (rg) and phenotypic (rp) correlations between traits were computed following Holland (2006). QTL analyses were performed with best linear unbiased predictor (BLUP) values of the RILs. Log scaling of BLUPs for fumonisin content was conducted to reduce the effect of large peaks in the data analysis (Xi et al., 2014).
The QTL analysis for FER and fumonisin content was conducted with the 144 RIL using the software PlabMQTL (Utz, 2012). The composite interval mapping (CIM) approach was used for QTL detection and to estimate QTL effects. Additive and epistatic (only for fumonisin content) effects were included in the model. Empirical LOD thresholds of 2.43, 3.83, and 1.93 were established to declare significant a putative QTL for days to silking, fumonisin content, and FER, respectively, according to executed permutation tests with 1,000 random reshuffles assuming an experiment-wise error of 0.3. A fivefold cross validation (CV/G) approach was employed for each trait using 1,000 CV runs to determine the frequency of QTL detection (Utz et al., 2000).
Four RILs with the lowest values (resistant) for fumonisin content (10–15 μg/g) and FER (∼2) along with four RILs with the highest values (susceptible, 55–75 μg/g for fumonisin content and ∼4 for FER) were tested in 2018. Fifteen seeds from each RIL were sown in a single row of 3.5 m, with the distance between adjacent rows being 1 m. The RILs were arranged in four pairs, each pair formed with one resistant and one susceptible RIL. RILs in each pair were planted in adjacent rows in order to minimize the contribution of plot differences to gene expression variations. The day on which each plant showed silks was recorded, and the plant was self-crossed; 15 days later, the main ear of the plant was inoculated as previously described. Ears were individually collected 10 days after inoculation; undamaged immature kernels around the inoculation point were removed and immediately stored in a Falcon tube, kept in liquid nitrogen, and stored at −80°C until RNA extraction.
Three replicates (ears) were used for each RIL. Twenty-four sample RNA extractions were performed using Maxwell® 16 LEV Plant RNA Kit with TRIzol in Maxwell® 16 Research Instrument (Promega, Madison, WI, United States). An RNA-seq analysis was conducted using bulks of resistant and susceptible RILs. Each of three resistant RNA bulks were made mixing equal amounts of one-ear RNA from each resistant RIL; and different ears from each RIL were assigned to different bulks. Three susceptible RNA bulks were made in a similar way. Therefore, six libraries were constructed, three for each bulk (resistant and susceptible). Bulk samples were barcoded and prepared for sequencing at Novogene Co., Ltd., where 150-bp paired-end (PE) reads were obtained on Illumina 1.9.
Quality control of Illumina reads was performed with Trimmomatic (Bolger et al., 2014) by removing poor quality reads and trimming poor quality bases. The percentages of retained reads after quality control ranged from 93.9 to 95.3% (Supplementary Table 1). Gene expression was quantified using Kallisto (Bray et al., 2016), whose algorithm uses a pseudo alignment approach to speed up the alignment procedure. The “pseudo alignment” approach makes it possible to quantify reads without recording actual alignments. The bootstrap resampling option to estimate count variability and the accuracy of quantification of each transcript per sample was set to 100. A DEG analysis was performed with Sleuth (Pimentel et al., 2017) by Wald Test to obtain Benjamini–Hochberg false discovery rate (FDR)-corrected p-values for each gene. DEGs were considered as those with FDR adjusted p-value < 0.05. The EnhancedVolcano R package was used to generate a volcano plot for gene expression (Blighe et al., 2021). To validate the RNA-seq results, the expression levels of seven genes were determined by qRT-PCR. cDNA was synthesized using Go Script Reverse Transcription System (Promega, Madison, WI, United States) following the manufacturer’s instructions and finally diluted 1:2 with nucleotide-free water. The primer sequences used for qPCR of target genes were designed based on gene sequences published in Maize Genomic Database using PRIMER 3 (Rozen and Skaletsky, 2000) and are listed in Supplementary Table 2. Real-time PCR was carried out using SYBR Green Reagents and BrightGREEN 2X qPCR master mix-low ROX (Applied Biological Materials, Richmond, BC, Canada), and three technical replicates were carried out to estimate variations due to the experimental procedure. The reaction conditions were as follows: 94°C for 1 min, followed by 38 cycles of 95°C for 10 s, 58°C for 10 s, and 72C° for 15 s. After amplification, the LinReg PCR V12.17 software (Ruijter et al., 2009; Tuomi et al., 2010) was used to calculate the efficiencies for each gene from the slopes of the standard curves. Three housekeeping genes were also amplified, and Gap C was chosen as the internal control to normalize expression data once results from the three internal controls were compared using the BestKeeper software to normalize the data (Pfaffl et al., 2004). The relative expressions for each gene were calculated following the mathematical model proposed by Pfaffl (2001) for relative quantification in RT-PCR.
Gene Ontology (GO) terms search for each DEG and functional enrichment were performed by combining results of the PANNZER (Koskinen et al., 2015) and ShinyGO v0.61 (Ge et al., 2020) web tools1; both programs were run with default parameters. Functions showing FDR-corrected p-values < 0.05 were considered enriched. The online search tool STRING was used for the retrieval of interacting genes2 (Version 11.0) (Szklarczyk et al., 2019) in order to find possible protein–protein interaction networks among the DEGs. High confidence interactions were only considered (interaction score >0.7).
Polymorphisms in Differentially Expressed Genes
The variant calling and Mpileup commands of the bcftools package (Li et al., 2009) were used to detect single nucleotide polymorphisms (SNPs) in the DEGs, which were filtered by quality score (>20) and a minimum depth read of 8. Only SNPs genotyped in all of the six samples were retained. Genomic regions containing DEGs with SNPs for which different bases were fixed in the resistant and susceptible bulks and/or INDELs with presence/absence in all the resistant bulk replicates and absence/presence in all the susceptible bulk replicates were considered as involved in resistance. In addition, polymorphic DEGs between the resistant and susceptible RILs located in the confidence intervals of the QTLs mapped were highlighted.
Quantitative Trait Loci for Fusarium Ear Rot and Kernel Fumonisin Content
Genotypic variability among the RILs was significant for days to silking, FER, and fumonisin content, while the genotype × year interaction was only significant for days to silking and FER (Supplementary Table 3). Histograms of RIL distributions for each trait are shown in Supplementary Figure 1, as well as the intervals where the RIL parents (A637 and EP42) were located and the heritability estimates for each trait on a mean basis. For each trait, many RILs showed phenotypic values higher or lower than those of the best and worst parents, respectively. Heritability estimates were always significantly different from zero. The mean heritability estimates ranged from 0.63 to 0.86.
The phenotypic correlation coefficient between FER and fumonisin content was moderate (r = 0.59 ± 0.04), but the genotypic correlation between both traits was high (0.85 ± 0.14). No significant genotypic correlation was found between FER and days to silking (0.1 ± 0.012), but the correlation coefficient between days to silking and fumonisin content, although low, was positive and significant (0.28 ± 0.12).
Models with two, three, and two QTL explained approximately 40, 30, and 17% of the genetic variation in days to silking, FER, and fumonisin content, respectively (Table 1). Additive models were appropriate for days to silking and FER, but significant additive and additive–additive effects were detected for QTLs for fumonisin content. QTLs for FER were detected at bins 4.07–4.1, 6–6.01, and 8.05–8.07, and those for fumonisin content were detected at bins 6.04–6.05 and 8.07–8.08. QTLs for fumonisin content at 8.07–8.08 clearly co-located with QTLs for FER at 8.05–8.07. Individually, A637 alleles at all QTLs reduced FER and fumonisin content, except at the QTL for FER at 4.08. However, the epistatic interaction between A637 alleles at the two QTLs for fumonisin content had positive and significant effects on that trait. No co-localization between QTLs for days to silking and fumonisin content was detected, although the genotypic correlation coefficient between these two traits was significantly different from zero. Although the genotype × environment interaction was significant for days to silking and FER, the QTLs detected showed similar or at least proportional effects on both years. The QTLs for days to silking at bins 6–6.02 and 8.04–8.05 and the QTLs for FER at bins 4.07–4.1 and 8.05–8.07 were identified in a moderate to high percentage of cross-validation runs. However, the QTLs at bin 6.00–6.01 for FER and QTLs for fumonisin content were less reliable, because they were identified in a low percentage of cross-validation runs.
Table 1. Location of quantitative trait loci (QTLs) mapped in the recombinant inbred line (RIL) population derived from EP42 × A637 for days to silking (DS), Fusarium ear rot (FER), and kernel fumonisin content (Fum) under inoculation with Fusarium verticillioides1.
A total of ∼371 million PE reads were generated; the average for each sample analyzed was approximately 62 million reads (Supplementary Table 1). After filtering, ∼365 million reads were retained, ranging from 48 to 66 million reads for each sample. Out of the 44,301 maize genes of the B73 RefGen_v4, 26, 864 genes were expressed in at least one RIL bulk (Supplementary Data Sheet 1).
Differential Gene Expression
Comparison of transcript counts between the resistant and susceptible RIL bulks revealed 364 DEGs (Supplementary Data Sheet 1); the Log2 of -fold change values ranged from −5.27 to 6.89 (Supplementary Data Sheet 1 and Figure 1). In the resistant RILs, 212 and 152 genes were upregulated and downregulated, respectively, compared with susceptible RILs. Most of the DEGs (87.9%) were annotated with information in the NCBI and UniProt non-redundant protein databases (Supplementary Data Sheet 2). In a GO enrichment analysis, none of the gene categories were significantly enriched with DEGs (corrected p-value < 0.05) (Supplementary Data Sheet 2).
Figure 1. Volcano plot for gene expression in resistant inbred bulks compared to susceptible inbred bulks using the EnhancedVolcano R package (https://github.com/kevinblighe/EnhancedVolcano). The dots located in the positive area stand for genes expressed higher in resistant inbreds, and dots located in the negative area stand for genes expressed higher in susceptible inbreds. As shown in graphic symbol, red stands for expression fold changes higher than 1.5-fold and adjusted p-value < 0.05; black stands for no differentially expressed genes (adjusted p-value > 0.05 or fold change smaller than 1.5).
Two genes, Zm00001d018341 (encoding a protein probably involved in splicing) and Zm00001d044975 (encoding a regulator of the flavonoid pathway), were only expressed in the resistant RILs, while Zm00001d048444 (encoding a protein related to sulfur starvation) and Zm00001d042228 (encoding a protein involved in the maintenance of undifferentiated cells), were only expressed in susceptible RILs. In addition, infected kernels of the resistant RILs were clearly different from those of the susceptible RILs in terms of transcript levels of genes involved in hormone biosynthesis and signaling pathways (Figure 2). In general, genes involved in strigolactone, brassinosteroid, and ethylene signaling were upregulated in the resistant RILs compared with the susceptible ones, while genes involved in the phytosulfokine, jasmonic acid, and salicylic acid signaling pathways were downregulated. Genes involved in auxin biosynthesis were downregulated, and those involved in gibberellin deactivation and auxin signaling repression and deactivation were upregulated in the resistant RILs. Alongside these patterns of gene expression, the resistant RILs showed changes in the transcript levels of genes implicated in the cell cycle, RNA metabolic processes, proteasome-ubiquitination system, post-translational modification of proteins, respiration, and primary and secondary metabolism. Some proteins encoded by these genes were strongly connected. Out of the proteins encoded by the 364 DEGs, 55 were connected with 77 edges (Figure 3). The interconnected proteins included those involved in breakdown of storage proteins, protein transport, proteasomal degradation of damaged or unneeded proteins, L-asparagine synthesis, cell division and differentiation, and RNA metabolism.
Figure 2. Expression Log2 fold change between inbred bulks resistant and susceptible to Fusarium ear rot and kernel contamination with fumonisins for differentially expressed genes (DEGs) related to hormone metabolism and signaling.
Figure 3. High-confidence protein–protein interaction networks (interaction score >0.7) among DEGs using the online search tool STRING for the retrieval of interacting genes (https://string-db.org/) (Version 11.0) (Szklarczyk et al., 2019).
Compared with the susceptible RILs, the resistant RILs showed downregulation of genes involved in cell proliferation and upregulation of genes involved in cell extension and differentiation, especially those related to cell wall biosynthesis and modification (Figure 4). The resistant RILs also showed upregulation of genes involved in amino acid biosynthesis, lipid catabolism, and respiratory electron transport, and downregulation of genes related to sink metabolism and glyoxylate and the tricarboxylic acid (TCA) cycle. Specifically, the resistant lines showed upregulation of genes implicated in fungal defense signaling and biosynthesis of secondary metabolites likely involved in defense but not genes involved in programmed cell death (Figure 5). The qRT-PCR results confirmed the RNA-seq results in that the transcripts of genes Zm00001d010840, Zm00001d053603, Zm00001d035562, Zm00001d037810, Zm00001d037921, and Zm00001d044975 were significantly more abundant in the resistant bulks than in the susceptible ones; and vice versa for Zm00001d048444 (Figure 6).
Figure 4. Distribution of DEGs between resistant and susceptible inbred bulks visualized with MapMan. Genes up and down-expressed in resistant vs. susceptible inbred bulks are in blue and red, respectively.
Figure 5. Expression Log2 fold change between inbred bulks resistant and susceptible to Fusarium ear rot and kernel contamination with fumonisins for differentially expressed genes (DEGs) related to detoxification, defense signaling, and biosynthesis of defensive compounds.
Figure 6. qRT-PCR expression levels of susceptible and resistant inbred bulks for seven genes (A–G). Each bar represents the normalized (referred to gap gene expression) expression level of the gene. Student’s t-tests were performed for comparing inbred bulk expressions, and differences between resistant and susceptible inbred bulks were significant at p < 0.05 for the seven genes.
Differential Allelic Fixation
One or more SNPs were detected between the resistant and susceptible bulks in 37 DEGs across chromosomes 1, 3, 5, 6, 8, and 9; INDELs with presence/absence in resistant/susceptible bulks or vice versa were detected in 112 DEGs located across all the chromosomes (Supplementary Data Sheet 3). We detected SNPs and INDEL polymorphisms in 3, 4, 13, and 4 DEGs located in the confidence intervals of QTL for FER and/or fumonisin content at bins 4.07–4.1, 6–6.01, 6.04–6.05, and 8.05–8.08, respectively (Table 2).
Table 2. Candidate genes (more probable in bold) for QTL for fumonisin content and/or FER referred to the reference genome B73 RefGen_v4.
Genomics of Fusarium Ear Rot and Kernel Contamination With Fumonisins
This mapping population showed genetic variability for FER and fumonisin content in the kernels, with moderate heritability values for both traits. The heritability estimates for FER were in the range of values recently estimated for other bi-parental mapping populations, while those for fumonisin content were higher (Morales et al., 2019). A wide genetic variability for FER and fumonisin content was expected in this mapping population, because it was derived from the cross between two inbreds with contrasting values for FER and fumonisin content, with A637 being more resistant to FER and having lower fumonisin content than EP42 (Santiago et al., 2013).
Consistent with previous studies, the phenotypic and genotypic correlation coefficients between FER and fumonisin content were moderate to high. This suggested that, if selecting by phenotype, FER would be preferable to fumonisin content as the selection trait, because the former is less laborious and expensive to evaluate (Bolduan et al., 2009; Löffler et al., 2010, 2011; Samayoa et al., 2019). In addition, neither significant phenotypic nor genotypic correlations were found between FER and days to silking, indicating that phenological differences among the inbreds would not affect differences in FER among the RILs. In contrast, fumonisin accumulation may be affected by kernel humidity or other environmental conditions, as a previous study has shown that fumonisin contamination risk began earlier in plants sown later (Cao et al., 2013). Nevertheless, as expected based on the results of previous studies, the different performance of the RILs across environments greatly contributed to phenotypic differences among the RILs for FER, providing further evidence that molecular markers would be more useful than phenotypic values in selection programs (Santiago et al., 2020).
Most of the genetic variations in FER and fumonisin content could not be explained by the detected QTLs. No single QTL explained more than 10% of the phenotypic variability for FER or fumonisin content. Transgressive segregation was common, because almost half of the individuals in the mapping population had FER values lower than that of the resistant parent, A637. Transgressive segregation could result from the complementarity of favorable alleles from different parents but may also be indicative of additive × additive interaction effects. According to these two non-exclusive hypotheses, previous studies have shown that favorable alleles for FER do not always come from the resistant parent, and that there are QTLs with significant epistatic effects for reduced fumonisin levels (Robertson-Hoyt et al., 2006; Butron et al., 2015; Maschietto et al., 2017; Morales et al., 2019).
Overall, the results indicate that resistance to FER and reduced fumonisin accumulation are highly polygenic traits, and that the genetic variation in FER and fumonisin accumulation in the inbred lines is due to many QTLs with additive and additive × additive minor effects (Santiago et al., 2020). Therefore, the value of selection by QTL-linked markers will be limited because most genetic variability is not captured by QTL models because of the highly polygenic nature of the two traits (resistance to FER and reduced fumonisin contamination). However, selection gains would occur if genomic selection for the polygenic background for the target trait is performed alongside introgression of resistance alleles. In this context, Guo et al. (2020) showed that the prediction accuracy of genomic selection for FER was improved using trait-related markers. In this study, although the proportion of phenotypic variability in FER explained by each QTL was low, the QTLs were stable across environments and were moderately reliable, because two of them occurred in more than half of the CV runs. Therefore, molecular markers flanking QTLs for FER on chromosomes 4 and 8 could be included as fixed effects in models to predict genotypic values for genomic selection, although previously, QTL mapping using higher marker density would be implemented to narrow down the confidence intervals of the QTLs. In addition, the genomic regions detected in this study as relevant for resistance to FER and/or to contamination with fumonisins (at bins 4.07–4.1, 6.00–6.01, 6.04–6.05, and 8.05–8.08) have been previously highlighted in QTL mapping studies using bi-parental populations (Robertson-Hoyt et al., 2006; Maschietto et al., 2017; Morales et al., 2019; Wu et al., 2020) and in genome-wide association studies (Chen et al., 2016; Coan et al., 2018; Butron et al., 2019; Samayoa et al., 2019; Guo et al., 2020; Wu et al., 2020). Moreover, a QTL for FER found at bin 4.08 in a previous study was confirmed using near-isogenic lines (NILs) (Wu et al., 2020). It has been suggested that the region at 4.07–4.08 is a genomic hotspot for maize resistance to multiple pathogens and to mycotoxin accumulation, because QTLs related to eight different diseases and a validated QTL for aflatoxin content co-localize in that region (Wisser et al., 2006; Mideros et al., 2014). Therefore, this region deserves special attention, because it may contain genes related to resistance to multiple diseases across diverse genetic backgrounds. At bin 8.07, Samayoa et al. (2019) found an SNP significantly associated with fumonisin accumulation and proposed Zm00001d012327 (putative ubiquitin-like-specific protease 2B gene) as the candidate gene.
Pathways Involved in Fusarium Infection
Quantitative trait loci studies not only identify markers for use in marker-assisted selection programs but also represent an intermediate step in locating genes involved in the inheritance of a trait. Because it is very difficult to find genes involved in resistance to FER and/or to fumonisin accumulation via QTL mapping, other approaches could be implemented to find genes related to these traits. Previous studies have detected DEGs between two unrelated inbred lines with contrasting resistance to FER in control plants and those inoculated with F. verticillioides. DEGs have been proposed as candidate genes for resistance; but differences in genetic backgrounds of inbred lines may have biased candidate gene identification in those studies (Lanubile et al., 2010, 2012a,b, 2014, 2017; Yuan et al., 2012, 2013; Campos-Bermudez et al., 2013; Yao et al., 2020). Our approach overcomes this problem, because we performed BSR, which combines BSA and RNA-seq. This minimizes biases due to background differences as long as each bulk contains individuals that are identical in one or more genomic regions implicated in a particular trait but arbitrary in all unlinked regions (Michelmore et al., 1991). Although BSA is more often used to map genes involved in monogenic inherited traits, it can be also a useful tool for mapping genes involved in quantitative traits such as resistance to FER and kernel fumonisin contamination (Liu et al., 2012; Du et al., 2017; Hou et al., 2018; Hao et al., 2019; Qin et al., 2020).
We detected 364 DEGs in immature kernels (25 days after pollination) between two groups of recombinant RILs with contrasting values for FER and fumonisin content after inoculation with F. verticillioides. However, no GO categories were significantly enriched among the DEGs. Similarly, Lanubile et al. (2014) found that DEGs between resistant and susceptible inbreds to FER were distributed among all functional classes.
However, DEGs in some specific pathways were, generally, upregulated or downregulated in resistant RILs compared with susceptible ones. Several genes involved in hormone biosynthesis, signaling, and homeostasis were differentially expressed between resistant and susceptible RILs, in accordance with the finding of Zhou et al. (2020) that genes involved in hormone signaling were among the target genes of F. verticillioides-responsive miRNAs. The plant hormones jasmonic acid, salicylic acid, and ethylene are directly involved in plant defense and regulate plant defense responses by cross-talk with other hormones such as abscisic acid, auxins, cytokinins, gibberellins, and strigolactones (Kazan and Lyons, 2014). In this study, ethylene signaling was upregulated in resistant vs. susceptible RILs; meanwhile salicylic and jasmonic acid signaling were slightly downregulated.
Upregulation of ethylene signaling was accompanied by increased transcript levels of genes involved in strigolactone and brassinosteroid signaling. Kapulnik et al. (2011) showed that strigolactones and ethylene affect some developmental processes through a common regulatory pathway in which ethylene is epistatic to strigolactones. The coordination of the ethylene-brassinosteroid pathway could be responsible for balancing resistance and growth (Liu et al., 2017). Therefore, ethylene signaling may regulate the functions of strigolactones and brassinosteroids in resistance to FER and fumonisin accumulation in kernels. All the three hormones can negatively affect the growth of Fusarium species: strigolactones were shown to strongly inhibit the growth of Fusarium oxysporum and Fusarium solani; brassinosteroid was shown to enhance the resistance of barley to Fusarium diseases; and ethylene signaling was shown to affect the relative abundance of metabolites contributing to resistance to Fusarium graminearum in maize (Dor et al., 2011; Ali et al., 2013; Zhou et al., 2019).
The downregulation of genes involved in jasmonic acid signaling in the resistant RILs was not unexpected. Christensen et al. (2014) provided genetic evidence for the involvement of jasmonic acid in maize defense against F. verticillioides, while Lanubile et al. (2015) showed that maize resistance depends on early activation of LOX genes, even though the activation of LOX genes was stronger in susceptible inbred lines than in resistant ones 14 dpi.
Auxin biosynthesis and signaling were downregulated in the resistant RILs compared with the susceptible ones. Again, this was not surprising, because pathogens have developed strategies via their effector repertoire to either interfere, or hijack phytohormone pathways to induce host susceptibility. Pathogen-induced auxin biosynthesis or modulated auxin signaling has been associated with increased host susceptibility. In this context, Brauer et al. (2019) found that the auxin pathway was one of the few susceptibility-associated pathways that functioned during wheat kernel infection by F. graminearum.
Among DEGs between resistant and susceptible RILs, there were also genes involved in cell cycle, respiration, and primary and secondary metabolism, suggesting that resistant RILs make a growth–defense trade-off to optimize their fitness during F. verticillioides infection (Yao et al., 2020). According to a possible growth–defense trade-off, sink metabolic processes such as fatty acid and starch biosynthesis were downregulated as well as phytosulfokine signaling, which regulates cell proliferation and growth, and many other genes involved in cell division. However, genes related to cell expansion and differentiation were upregulated, especially those involved in cell wall biosynthesis and modification that could be controlled and organized, in part, by the cytoskeleton and mediated by auxin transport (Perrot-Rechenmann, 2010; Bashline et al., 2014). Lanubile et al. (2014) also reported that the response of one resistant inbred line to F. verticillioides infection included activation of genes related to lignin, a structural component of the plant cell wall. Induced systemic resistance (ISR) in maize is activated by the beneficial fungus Trichoderma atroviride, and this process involves reinforcement of plant cell walls, which reduces the development of disease after inoculation with F. verticillioides (Agostini et al., 2019).
In this study, glycolysis, tricarboxylic acid, and glyoxylate cycles were downregulated concomitantly with reduction in sink processes such as starch and fatty acid biosynthesis in the resistant RILs. However, genes involved in electron transfer during respiratory metabolism were upregulated, as were those involved in amino acid biosynthesis. This suggested that the electron transport chain could be activated toward amino acid synthesis under conditions of high TCA cycle flux (Fernie et al., 2004; Li et al., 2020). Ciasca et al. (2020) recently showed that levels of the amino acid tyrosine were significantly higher in an inbred line resistant to F. verticillioides infection than in its susceptible counterpart, and that the opposite was true for L-tryptophan.
In terms of plant defense pathways, genes directly involved in defense mechanisms were activated in the resistant RILs, while genes involved in cell death were repressed. In previous transcriptomic studies, activation of genes involved in secondary metabolism, such as those in shikimate, lignin, flavonoid, and terpenoid pathways, has been observed in resistant vs. susceptible genotypes, or in plants colonized by beneficial fungi (Lanubile et al., 2014; Agostini et al., 2019). In this study, one of the most strongly activated secondary pathways was the flavonoid branch of the phenylpropanoid pathway. This is consistent with reports that flavonoid compounds such as flavones, flavonols, and anthocyanins rather than hydroxycinnamic acids contribute to differences in FER resistance and fumonisin accumulation among maize cultivars (Venturini et al., 2016; Bernardi et al., 2018). The antioxidant properties of flavonoids could also explain why genes involved in reactive oxygen species scavenging and detoxification were, in general, downregulated in resistant RILs, although another study detected them among upregulated genes in resistant vs. susceptible inbred lines (Maschietto et al., 2016).
Overall, the increased secondary metabolism in resistant RILs could be supported by decreased sink metabolism and increased catabolic processes to produce substrates for secondary metabolism. However, some authors have highlighted more direct roles of plant protein and pectin degradation, for example, plant proteolytic systems are crucial for defense against pathogens, and oligogalacturonides are elicitors of plant defense responses (Baek and Choi, 2008; Lionetti, 2015). In addition, some DEGs encode proteins that may be involved in degradation of fungal defenses, for example, secretory lipases that degrade the fungal arsenal of lipids and specific glucanases that target the fungal cell wall (Lozovaya et al., 1998; Debona et al., 2017; Lee and Park, 2019).
Besides giving an overview on differences in gene expression between resistant and susceptible inbred lines, RNA-seq analyses can help to uncover genomic regions involved in resistance. To decrease the probability of detecting false positive associations, only polymorphisms in DEGs resulting from fixation of one allele in resistant inbreds and fixation of the alternative allele in susceptible inbreds, along with INDEL polymorphisms (presence/absence in resistant/susceptible inbreds or vice versa), were used to map genomic regions likely associated with resistance. Resistant RILs presented resistant alleles for at least one marker located in the confidence interval of each QTL, while susceptible RIL lines carried susceptible alleles for at least one marker contained in the confidence interval of each QTL. DEGs with polymorphisms between resistant and susceptible inbreds were scattered across 43 bins of the 10 chromosomes, consistent with the highly polygenic nature of the studied traits (Zila et al., 2014; Samayoa et al., 2019).
The DEGs located within confidence intervals of detected QTLs and showing polymorphisms between resistant and susceptible RILs were identified as valuable candidate genes for those QTLs and deserve special attention. Differences in expression were especially noticeable for Zm00001d053603, Zm00001d035562, Zm00001d037810, Zm00001d037921, and Zm00001d010840 and were confirmed by qRT-PCR. Zm00001d053603 encodes an ABC transporter A family member involved in lipid transport, and Zm00001d010840 encodes triacylglycerol lipase-like 1, which is involved in degradation of oil bodies in seeds. The upregulation of these genes could be related to increased resistance to FER, because they would participate in mobilization of lipids from oil bodies to phytoalexin synthesis (Shimada et al., 2014). Similarly, Zm00001d035562, encoding a dirigent protein, 100502041, that yields optically active lignans from two molecules of coniferyl alcohol, could be implicated in resistance to FER through its involvement in lignan biosynthesis and/or other aspects of secondary cell wall biosynthesis (Corbin et al., 2018). Zm00001d037810 encodes a G-type lectin S-receptor-like serine/threonine-protein kinase, At2g19130. Lectin receptor-like kinases are cell surface receptors that play important roles in perceiving and processing signals that arrive in cells, such as fungi, because their extracellular domains can bind to fungal cell wall components (Vaid et al., 2013; Liu et al., 2018). Finally, Zm00001d037921 is orthologous to the MWD9.13 gene in Arabidopsis and encodes an NF-κB inhibitor-like protein. In humans, the NF-κB family of transcription factors regulates the expression of hundreds of genes associated with cell proliferation, differentiation, and death, and immune responses (Golan-Goldhirsh and Gopas, 2014). In Arabidopsis, homologs of NF-κB transcription factors have been proposed to participate in redox sensing (Li et al., 2006). Therefore, polymorphisms in the gene encoding NF-κB inhibitor-like protein may interfere in oxidative response to F. verticillioides and contribute to resistance to FER and fumonisin accumulation, since increased oxidative stress favors mycotoxin biosynthesis (Picot et al., 2010). According to the hypothesis of reduced oxidative stress in resistant kernels, genes related to cell death were downregulated in the resistant compared with the susceptible RILs. These candidate genes should be validated in future studies.
Functions of the DEGs between the resistant and susceptible inbred bulks suggest a growth–defense trade-off, because sink metabolic processes such as fatty acid and starch biosynthesis were downregulated, as well as phytosulfokine signaling and many other genes involved in cell division, while genes related to substances and processes of the secondary metabolism associated to resistance, such as those involved in cell wall biosynthesis and its adjustment, and flavonoid biosynthesis, were upregulated. Among the DEGs, five (Zm00001d053603, Zm00001d035562, Zm00001d037810, Zm00001d037921, and Zm00001d010840) were polymorphic between the resistant and susceptible bulks, were located within the confidence intervals of detected QTLs and showed large differences in transcript abundance between the resistant and susceptible bulks. Thus, they were identified as candidate genes for the detected QTLs.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: http://hdl.handle.net/10261/258356, https://doi.org/10.20350/digitalCSIC/14486.
AB conceived the study. AC took care of the field and greenhouse experiments, data recording, and sample collection with the assistance of AB, RS, and RM. MF led all the molecular biology analyses. AB and MF performed the statistical analyses of data. AB and NG drafted the initial manuscript. NG edited the manuscript. All authors read and approved the final version of the manuscript.
This research was funded by a coordinated project (subprojects RTI2018-096776-B-C21 and RTI2018-096776-B-C22) and also financed by MCIU/AEI/FEDER, UE. Open access publication fees were covered by the CSIC and the project IN607A 2021/03 (Xunta de Galicia).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
AC acknowledges a PTA contract financed by MCIU/AEI/FEDER, UE. Thanks to Yolanda Arnaiz for her invaluable help in the RNA extraction and qRT-PCR analyses. NG acknowledges an FPI contract (“Contrato predoctoral para la formación de doctores”) funded by project RTI2018-096776-B-C21. We thank Jennifer Smith, Ph.D., from Edanz (https://www.edanz.com/ac) for editing a draft of this manuscript. We also acknowledge the support of the CSIC Interdisciplinary Thematic Platform (PTI) Optimization of agricultural and forestry systems (PTI-AGROFOR).
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.866478/full#supplementary-material
Agostini, R. B., Postigo, A., Rius, S. P., Rech, G. E., Campos-Bermudez, V. A., and Vargas, W. A. (2019). Long-lasting primed state in maize plants: salicylic acid and steroid signaling pathways as key players in the early activation of immune responses in silks. Mol. Plant-Microb. Interact. 32, 95–106. doi: 10.1094/MPMI-07-18-0208-R
Alizadeh, A. M., Rohandel, G., Roudbarmohammadi, S., Roudbary, M., Sohanaki, H., Ghiasian, S. A., et al. (2012). Fumonisin B1 contamination of cereals and risk of esophageal cancer in a high risk area in northeastern Iran. Asian Pac. J. Cancer Prev. 13, 2625–2628. doi: 10.7314/apjcp.2012.13.6.2625
Bernardi, J., Stagnati, L., Lucini, L., Rocchetti, G., Lanubile, A., Cortellini, C., et al. (2018). Phenolic profile and susceptibility to Fusariuminfection of pigmented maize cultivars. Front. Plant Sci. 9:13. doi: 10.3389/fpls.2018.01189
Blighe, K., Rana, S., and Lewis, M. (2021). EnhancedVolcano: Publication-ready volcano plots with enhanced colouring and labeling. Available online at: https://github.com/kevinblighe/EnhancedVolcano (Accessed October 29, 2021).
Bolduan, C., Miedaner, T., Schipprack, W., Dhillon, B. S., and Melchinger, A. E. (2009). Genetic variation for resistance to ear rots and mycotoxins contamination in early European maize inbred lines. Crop Sci. 49, 2019–2028. doi: 10.2135/cropsci2008.12.0701
Brauer, E. K., Rocheleau, H., Balcerzak, M., Pan, Y. L., Fauteux, F., Liu, Z. Y., et al. (2019). Transcriptional and hormonal profiling of Fusarium graminearum-infected wheat reveals an association between auxin and susceptibility. Physiol. Mol. Plant Pathol. 107, 33–39. doi: 10.1016/j.pmpp.2019.04.006
Butron, A., Reid, L. M., Santiago, R., Cao, A., and Malvar, R. A. (2015). Inheritance of maize resistance to gibberella and fusarium ear rots and kernel contamination with deoxynivalenol and fumonisins. Plant Pathol. 64, 1053–1060. doi: 10.1111/ppa.12351
Butron, A., Santiago, R., Cao, A., Samayoa, L. F., and Malvar, R. A. (2019). QTLs for resistance to Fusarium ear rot in a multiparent advanced generation intercross (MAGIC) maize population. Plant Dis. 103, 897–904. doi: 10.1094/PDIS-09-18-1669-RE
Butron, A., Santiago, R., Mansilla, P., Pintos-Varela, C., Ordas, A., and Ana Malvar, R. (2006). Maize (Zea mays L.) genetic factors for preventing fumonisin contamination. J. Agricult. Food Chem. 54, 6113–6117. doi: 10.1021/jf0611163
Campos-Bermudez, V. A., Fauguel, C. M., Tronconi, M. A., Casati, P., Presello, D. A., and Andreo, C. S. (2013). Transcriptional and metabolic changes associated to the Infection by F. verticillioides in maize inbreds with contrasting ear rot resistance. PLos One 8:10. doi: 10.1371/journal.pone.0061580
Cao, A., Butron, A., Ramos, A. J., Marin, S., Souto, C., and Santiago, R. (2014). Assessing white maize resistance to fumonisin contamination. Eur. J. Plant Pathol. 138, 283–292. doi: 10.1007/s10658-013-0328-y
Cao, A., Santiago, R., Ramos, A. J., Marin, S., Reid, L. M., and Butron, A. (2013). Environmental factors related to fungal infection and fumonisin accumulation during the development and drying of white maize kernels. Int. J. Food Microbiol. 164, 15–22. doi: 10.1016/j.ijfoodmicro.2013.03.012
Chen, J., Shrestha, R., Ding, J. Q., Zheng, H. J., Mu, C. H., Wu, J. Y., et al. (2016). Genome-wide association study and QTL mapping reveal genomic loci associated with Fusarium ear rot resistance in Tropical maize germplasm. G3-Genes Genomes Genet. 6, 3803–3815. doi: 10.1534/g3.116.034561
Christensen, S. A., Nemchenko, A., Park, Y. S., Borrego, E., Huang, P. C., Schmelz, E. A., et al. (2014). The Novel Monocot-Specific 9-Lipoxygenase ZmLOX12 is required to mount an effective jasmonate-mediated defense against fusarium verticillioides in maize. Mol. Plant-Microbe Interact. 27, 1263–1276. doi: 10.1094/MPMI-06-13-0184-R
Ciasca, B., Lanubile, A., Marocco, A., Pascale, M., Logrieco, A. F., and Lattanzio, V. M. T. (2020). Application of an integrated and open source workflow for LC-HRMS plant metabolomics studies. Case-control study: metabolic changes of maize in response toFusarium verticillioides infection. Front. Plant Sci. 11:15. doi: 10.3389/fpls.2020.00664
Coan, M. M. D., Senhorinho, H. J. C., Pinto, R. J. B., Scapim, C. A., Tessmann, D. J., Williams, W. P., et al. (2018). Genome-wide association study of resistance to ear rot by Fusarium verticillioides in a Tropical field maize and popcorn core collection. Crop Sci. 58, 564–578. doi: 10.2135/cropsci2017.05.0322
Corbin, C., Drouet, S., Markulin, L., Auguin, D., Laine, E., Davin, L. B., et al. (2018). A genome-wide analysis of the flax (Linum usitatissimum L.) dirigent protein family: from gene identification and evolution to differential regulation. Plant Mol. Biol. 97, 73–101. doi: 10.1007/s11103-018-0725-x
Debona, D., Cruz, M. F. A., and Rodrigues, F. A. (2017). Calcium-triggered accumulation of defense-related transcripts enhances wheat resistance to leaf blast. Tropic. Plant Pathol. 42, 309–314. doi: 10.1007/s40858-017-0144-6
Dor, E., Joel, D. M., Kapulnik, Y., Koltai, H., and Hershenhorn, J. (2011). The synthetic strigolactone GR24 influences the growth pattern of phytopathogenic fungi. Planta 234, 419–427. doi: 10.1007/s00425-011-1452-6
Du, H. W., Zhu, J. X., Su, H., Huang, M., Wang, H. W., Ding, S. C., et al. (2017). Bulked segregant RNA-seq reveals differential expression and SNPs of candidate genes associated with waterlogging tolerance in maize. Front. Plant Sci. 8:13. doi: 10.3389/fpls.2017.01022
Eller, M. S., Payne, G. A., and Holland, J. B. (2010). Selection for reduced Fusarium ear rot and fumonisin content in advanced backcross maize lines and their topcross hybrids. Crop Sci. 50, 2249–2260. doi: 10.2135/cropsci2009.11.0683
Fernie, A. R., Carrari, F., and Sweetlove, L. J. (2004). Respiratory metabolism: glycolysis, the TCA cycle and mitochondrial electron transport. Curr. Opin. Plant Biol. 7, 254–261. doi: 10.1016/j.pbi.2004.03.007
Gelderblom, W. C. A., Jaskiewicz, K., Marasas, W. F. O., Thiel, P. G., Horak, R. M., Vleggaar, R., et al. (1988). Fumonisins - Novel mycotoxins with cancer-promoting activity produced by Fusarium moniliforme. Appl. Environ. Microbiol. 54, 1806–1811. doi: 10.1128/aem.54.7.1806-1811.1988
Gelineau-van Waes, J., Starr, L., Maddox, J., Aleman, F., Voss, K. A., Wilberding, J., et al. (2005). Maternal fumonisin exposure and risk for neural tube defects: mechanisms in an in vivo mouse model. Birth Defect. Res. Clin. Mol. Teratol. 73, 487–497. doi: 10.1002/bdra.20148
Gesteiro, N., Cao, A., Santiago, R., Malvar, R. A., and Butrón, A. (2021). Genomics of maize resistance to kernel contamination with fumonisins using a multiparental advanced generation InterCross maize population (MAGIC). BMC Plant Biol. 21:596. doi: 10.1186/s12870-021-03380-0
Guo, Z. F., Zou, C., Liu, X. G., Wang, S. H., Li, W. X., Jeffers, D., et al. (2020). Complex genetic system involved in fusarium ear rot resistance in maize as revealed by GWAS, bulked sample analysis, and genomic prediction. Plant Dis. 104, 1725–1735. doi: 10.1094/PDIS-07-19-1552-RE
Hao, Z. M., Geng, M. M., Hao, Y. R., Zhang, Y., Zhang, L. J., Wen, S. M., et al. (2019). Screening for differential expression of genes for resistance to Sitodiplosis mosellana in bread wheat via BSR-seq analysis. Theor. Appl. Genet. 132, 3201–3221. doi: 10.1007/s00122-019-03419-9
Holland, J. B. (2006). Estimating genotypic correlations and their standard errors using multivariate restricted maximum likelihood estimation with SAS Proc MIXED. Crop Sci. 46, 642–656. doi: 10.2135/cropsci2005.0191
Hou, X. G., Guo, Q., Wei, W. G., Guo, L. L., Guo, D. L., and Zhang, L. (2018). Screening of genes related to early and late flowering in tree peony based on bulked segregant RNA sequencing and verification by quantitative real-time PCR. Molecules 23:19. doi: 10.3390/molecules23030689
Ju, M., Zhou, Z. J., Mu, C., Zhang, X. C., Gao, J. Y., Liang, Y. K., et al. (2017). Dissecting the genetic architecture of Fusarium verticillioides seed rot resistance in maize by combining QTL mapping and genome-wide association analysis. Sci. Rep. 7:46446. doi: 10.1038/srep46446
Kapulnik, Y., Resnick, N., Mayzlish-Gati, E., Kaplan, Y., Wininger, S., Hershenhorn, J., et al. (2011). Strigolactones interact with ethylene and auxin in regulating root-hair elongation in Arabidopsis. J. Exp. Bot. 62, 2915–2924. doi: 10.1093/jxb/erq464
Koskinen, P., Toronen, P., Nokso-Koivisto, J., and Holm, L. (2015). PANNZER: high-throughput functional annotation of uncharacterized proteins in an error-prone environment. Bioinformatics 31, 1544–1552. doi: 10.1093/bioinformatics/btu851
Lanubile, A., Bernardi, J., Battilani, P., Logrieco, A., and Marocco, A. (2012a). Resistant and susceptible maize genotypes activate different transcriptional responses against Fusarium verticillioides. Physiol. Mol. Plant Pathol. 77, 52–59. doi: 10.1186/1471-2164-15-710
Lanubile, A., Bernardi, J., Marocco, A., Logrieco, A., and Paciolla, C. (2012b). Differential activation of defense genes and enzymes in maize genotypes with contrasting levels of resistance to Fusarium verticillioides. Environ. Exp. Bot. 78, 39–46.
Lanubile, A., Ferrarini, A., Maschietto, V., Delledonne, M., Marocco, A., and Bellin, D. (2014). Functional genomic analysis of constitutive and inducible defense responses to Fusarium verticillioides infection in maize genotypes with contrasting ear rot resistance. BMC Genom. 15:710.
Lanubile, A., Maschietto, V., Borrelli, V. M., Stagnati, L., Logrieco, A. F., and Marocco, A. (2017). Molecular basis of resistance to Fusarium ear rot in maize. Front. Plant Sci. 8:1774. doi: 10.3389/fpls.2017.01774
Lanubile, A., Maschietto, V., De Leonardis, S., Battilani, P., Paciolla, C., and Marocco, A. (2015). Defense responses to mycotoxin-producing fungi Fusarium proliferatum, F. subglutinans, and Aspergillus flavus in kernels of susceptible and resistant maize genotypes. Mol. Plant-Microbe Interact. 28, 546–557. doi: 10.1094/MPMI-09-14-0269-R
Lanubile, A., Pasini, L., and Marocco, A. (2010). Differential gene expression in kernels and silks of maize lines with contrasting levels of ear rot resistance after Fusarium verticillioides infection. J. Plant Physiol. 167, 1398–1406. doi: 10.1016/j.jplph.2010.05.015
Li, P. H., Mane, S. P., Sioson, A. A., Robinet, C. V., Heath, L. S., Bohnert, H. J., et al. (2006). Effects of chronic ozone exposure on gene expression in Arabidopsis thaliana ecotypes and in Thellungielia halophila. Plant Cell Environ. 29, 854–868. doi: 10.1111/j.1365-3040.2005.01465.x
Li, Y. W., Wang, X. J., Zeng, Y. L., and Liu, P. (2020). Metabolic profiling reveals local and systemic responses of kiwifruit to Pseudomonas syringae pv. actinidiae. Plant Direct 4:e00297. doi: 10.1002/pld3.297
Lin, M. M., Sun, S. H., Fang, J. B., Qi, X. J., Sun, L. M., Zhong, Y. P., et al. (2021). BSR-Seq analysis provides insights into the cold stress response of Actinidia arguta F1 populations. BMC Genom. 22:13. doi: 10.1186/s12864-021-07369-9
Lionetti, V. (2015). PECTOPLATE: the simultaneous phenotyping of pectin methylesterases, pectinases, and oligogalacturonides in plants during biotic stresses. Front. Plant Sci. 6:331. doi: 10.3389/fpls.2015.00331
Liu, H., Dong, S. Y., Gu, F. W., Liu, W., Yang, G. L., Huang, M., et al. (2017). NBS-LRR protein Pik-H4 interacts with OsBIHD1 to balance rice blast resistance and growth by coordinating ethylene-brassinosteroid pathway. Front. Plant Sci. 8:127. doi: 10.3389/fpls.2017.00127
Liu, P. L., Huang, Y., Shi, P. H., Yu, M., Xie, J. B., and Xie, L. L. (2018). Duplication and diversification of lectin receptor-like kinases (LecRLK) genes in soybean. Sci. Rep. 8:5861. doi: 10.1038/s41598-018-24266-6
Löffler, M., Kessel, B., Ouzunova, M., and Miedaner, T. (2011). Covariation between line and testcross performance for reduced mycotoxin concentrations in European maize after silk channel inoculation of two Fusarium species. Theor. Appl. Genet. 122, 925–934. doi: 10.1007/s00122-010-1499-y
Löffler, M., Miedaner, T., Kessel, B., and Ouzunova, M. (2010). Mycotoxin accumulation and corresponding ear rot rating in three maturity groups of European maize inoculated by two Fusarium species. Euphytica 174, 153–164. doi: 10.1007/s10681-009-0080-8
Maschietto, V., Colombi, C., Pirona, R., Pea, G., Strozzi, F., Marocco, A., et al. (2017). QTL mapping and candidate genes for resistance to Fusarium ear rot and fumonisin contamination in maize. BMC Plant Biol. 17:20. doi: 10.1186/s12870-017-0970-1
Maschietto, V., Lanubile, A., De Leonardis, S., Marocco, A., and Paciolla, C. (2016). Constitutive expression of pathogenesis-related proteins and antioxydant enzyme activities triggers maize resistance towards Fusarium verticillioides. J. Plant Physiol. 200, 53–61. doi: 10.1016/j.jplph.2016.06.006
Maschietto, V., Marocco, A., Malachova, A., and Lanubile, A. (2015). Resistance to Fusarium verticillioides and fumonisin accumulation in maize inbred lines involves an earlier and enhanced expression of lipoxygenase (LOX) genes. J. Plant Physiol. 188, 9–18. doi: 10.1016/j.jplph.2015.09.003
Michelmore, R. W., Paran, I., and Kesseli, R. V. (1991). Identification of markers linked to disease-resistance genes by bulked segregant analysis - a rapid method to detect markers in specific genomic regions by using segregating populations. Proc. Nat. Acad. Sci. USA 88, 9828–9832. doi: 10.1073/pnas.88.21.9828
Mideros, S. X., Warburton, M. L., Jamann, T. M., Windham, G. L., Williams, W. P., and Nelson, R. J. (2014). Quantitative trait loci influencing mycotoxin contamination of maize: analysis by linkage mapping, characterization of near-isogenic lines, and meta-analysis. Crop Sci. 54, 127–142. doi: 10.2135/cropsci2013.04.0249
Missmer, S. A., Suarez, L., Felkner, M., Wang, E., Merrill, A. H., Rothman, K. J., et al. (2006). Exposure to fumonisins and the occurrence of neural tube defects along the Texas-Mexico border. Environ. Health Persp. 114, 237–241. doi: 10.1289/ehp.8221
Morales, L., Zila, C. T., Mejia, D. E. M., Arbelaez, M. M., Balint-Kurti, P. J., Holland, J. B., et al. (2019). Diverse components of resistance to fusarium verticillioides infection and fumonisin contamination in four maize recombinant inbred families. Toxins 11:17. doi: 10.3390/toxins11020086
Padua, J. M. V., Dias, K. O. D., Pastina, M. M., de Souza, J. C., Queiroz, V. A. V., da Costa, R. V., et al. (2016). A multi-environment trials diallel analysis provides insights on the inheritance of fumonisin contamination resistance in tropical maize. Euphytica 211, 277–285. doi: 10.1007/s10681-016-1722-2
Pérez-Brito, D., Jeffers, D., González-de-León, D., Khairallah, M., Cortés-Cruz, M., Velázquez-Cardelas, G., et al. (2001). QTL Mapping of Fusarium moniliforme ear rot resistance in highland maize. Mexico. Agrociencia 35, 181–196.
Pfaffl, M. W., Tichopad, A., Prgomet, C., and Neuvians, T. P. (2004). Determination of stable housekeeping genes, differentially regulated target genes and sample integrity: BestKeeper - Excel-based tool using pair-wise correlations. Biotechnol. Lett. 26, 509–515. doi: 10.1023/b:bile.0000019559.84305.47
Picot, A., Barreau, C., Pinson-Gadais, L., Caron, D., Lannou, C., and Richard-Forget, F. (2010). Factors of the Fusarium verticillioides-maize environment modulating fumonisin production. Crit. Rev. Microbiol. 36, 221–231. doi: 10.3109/10408411003720209
Pimentel, H., Bray, N. L., Puente, S., Melsted, P., and Pachter, L. (2017). Differential analysis of RNA-seq incorporating quantification uncertainty. Nat. Methods 14, 687–690. doi: 10.1038/nmeth.4324
Presello, D. A., Pereyra, A. O., Iglesias, J., Fauguel, C. M., Sampietro, D. A., and Eyherabide, G. H. (2011). Responses to selection of S-5 inbreds for broad-based resistance to ear rots and grain mycotoxin contamination caused by Fusarium spp. in maize. Euphytica 178, 23–29. doi: 10.1007/s10681-010-0255-3
Qin, X. E., Zhang, W. L., Dong, X., Tian, S. K., Zhang, P. P., Zhao, Y. X., et al. (2020). Identification of fertility-related genes for maize CMS-S via Bulked Segregant RNA-Seq. Peerj 8:22. doi: 10.7717/peerj.10015
Robertson, L. A., Kleinschmidt, C. E., White, D. G., Payne, G. A., Maragos, C. M., and Holland, J. B. (2006). Heritabilities and correlations of fusarium ear rot resistance and fumonisin contamination resistance in two maize populations. Crop Sci. 46, 353–361.
Robertson-Hoyt, L. A., Jines, M. P., Balint-Kurti, P. J., Kleinschmidt, C. E., White, D. G., Payne, G. A., et al. (2006). QTL mapping for fusarium ear rot and fumonisin contamination resistance in two maize populations. Crop Sci. 46, 1734–1743. doi: 10.2135/cropsci2005.12-0450
Rodríguez Estrada, A. E., Jonkers, W., Kistler, H. C., and May, G. (2012). Interactions between Fusarium verticillioides, Ustilago maydis, and Zea mays: an endophyte, a pathogen, and their shared plant host. Fungal Genet. Biol. 49, 578–587. doi: 10.1016/j.fgb.2012.05.001
Ruijter, J. M., Ramakers, C., Hoogaars, W. M. H., Karlen, Y., Bakker, O., van den Hoff, M. J. B., et al. (2009). Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 37:e45. doi: 10.1093/nar/gkp045
Samayoa, L. F., Cao, A., Santiago, R., Malvar, R. A., and Butron, A. (2019). Genome-wide association analysis for fumonisin content in maize kernels. BMC Plant Biol. 19:166. doi: 10.1186/s12870-019-1759-1
Santiago, R., Cao, A., and Butron, A. (2015). Genetic factors involved in fumonisin accumulation in maize kernels and their implications in maize agronomic management and breeding. Toxins 7, 3267–3296. doi: 10.3390/toxins7083267
Santiago, R., Cao, A., Malvar, R. A., Reid, L. M., and Butron, A. (2013). Assessment of corn resistance to fumonisin accumulation in a broad collection of inbred lines. Field Crops Res. 149, 193–202. doi: 10.1016/j.fcr.2013.04.011
Shimada, T. L., Takano, Y., Shimada, T., Fujiwara, M., Fukao, Y., Mori, M., et al. (2014). Leaf oil body functions as a subcellular factory for the production of a phytoalexin in arabidopsis. Plant Physiol. 164, 105–118. doi: 10.1104/pp.113.230185
Sun, G., Wang, S., Hu, X., Su, J., Huang, T., Yu, J., et al. (2007). Fumonisin B1 contamination of home-grown corn in high-risk areas for esophageal and liver cancer in China. Food Addit. Contam. 24, 181–185. doi: 10.1080/02652030601013471
Szklarczyk, D., Gable, A. L., Lyon, D., Junge, A., Wyder, S., Huerta-Cepas, J., et al. (2019). STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 47, D607–D613. doi: 10.1093/nar/gky1131
Tuomi, J. M., Voorbraak, F., Jones, D. L., and Ruijter, J. M. (2010). Bias in the C-q value observed with hydrolysis probe based quantitative PCR can be corrected with the estimated PCR efficiency value. Methods 50, 313–322. doi: 10.1016/j.ymeth.2010.02.003
Utz, H. F. (2012). PlabMQTL - Software for meta-QTL analysis with composite interval mapping. Version 0.9. Institute of Plant Breeding, Seed Science, and Population Genetics. Stuttgart: University of Hohenheim.
Utz, H. F., Melchinger, A. E., and Schön, C. C. (2000). Bias and sampling error of the estimated proportion of genotypic variance explained by quantitative trait loci determined from experimental data in maize using cross validation and validation with independent samples. Genetics 154, 1839–1849. doi: 10.1093/genetics/154.4.1839
Venturini, G., Babazadeh, L., Casati, P., Pilu, R., Salomoni, D., and Toffolatti, S. L. (2016). Assessing pigmented pericarp of maize kernels as possible source of resistance to fusarium ear rot, Fusarium spp. infection and fumonisin accumulation. Int. J. Food Microbiol. 227, 56–62. doi: 10.1016/j.ijfoodmicro.2016.03.022
Wang, Y. P., Zhou, Z. J., Gao, J. Y., Wu, Y. B., Xia, Z. L., Zhang, H. Y., et al. (2016). The mechanisms of maize resistance to Fusarium verticillioides by comprehensive analysis of RNA-seq data. Front. Plant Sci. 7:1654. doi: 10.3389/fpls.2016.01654
Wisser, R. J., Balint-Kurti, P. J., and Nelson, R. J. (2006). The genetic architecture of disease resistance in maize: a synthesis of published studies. Phytopathology 96, 120–129. doi: 10.1094/PHYTO-96-0120
Wu, Y. B., Zhou, Z. J., Dong, C. P., Chen, J. F., Ding, J. Q., Zhang, X. C., et al. (2020). Linkage mapping and genome-wide association study reveals conservative QTL and candidate genes for Fusarium rot resistance in maize. BMC Genomics 21:11. doi: 10.1186/s12864-020-6733-7
Xi, B., Gu, H., Baniasadi, H., and Raftery, D. (2014). Statistical analysis and modeling of mass spectrometry-based metabolomics data. Methods Mol. Biol. 1198, 333–353. doi: 10.1007/978-1-4939-1258-2_22
Yao, L. S., Li, Y. M., Ma, C. Y., Tong, L. X., Du, F. L., and Xu, M. L. (2020). Combined genome-wide association study and transcriptome analysis reveal candidate genes for resistance toFusariumear rot in maize. J. Integ. Plant Biol. 62, 1535–1551. doi: 10.1111/jipb.12911
Yuan, G. S., Zhang, Z. M., Xiang, K., Shen, Y. O., Du, J., Lin, H. J., et al. (2013). Different gene expressions of resistant and susceptible maize inbreds in response to Fusarium verticillioides infection. Plant Mol. Biol. Rep. 31, 925–935. doi: 10.1007/s11105-013-0567-2
Yuan, G. S., Zhang, Z. M., Xiang, K., Zhao, M. J., Shen, Y. O., and Pan, G. T. (2012). Large-scale identification of differentially expressed genes in maize inbreds susceptible and resistant to Fusarium ear rot. Plant Omics 5, 471–475.
Zhou, S. Q., Zhang, Y. K., Kremling, K. A., Ding, Y. Z., Bennett, J. S., Bae, J. S., et al. (2019). Ethylene signaling regulates natural variation in the abundance of antifungal acetylated diferuloylsucroses and Fusarium graminearum resistance in maize seedling roots. New Phytol. 221, 2096–2111. doi: 10.1111/nph.15520
Zhou, Z. J., Cao, Y., Li, T., Wang, X. H., Chen, J. F., He, H., et al. (2020). MicroRNAs are involved in maize immunity against fusarium verticillioides ear rot. Genomics Proteom. Bioinform. 18, 241–255. doi: 10.1016/j.gpb.2019.11.006
Zila, C. T., Ogut, F., Romay, M. C., Gardner, C. A., Buckler, E. S., and Holland, J. B. (2014). Genome-wide association study of Fusarium ear rot disease in the U.S.A. maize inbred line collection. BMC Plant Biol. 14:372. doi: 10.1186/s12870-014-0372-6
Keywords: maize, Fusarium verticillioides, fumonisins, resistance, QTL, differentially expressed genes (DEG)
Citation: Cao A, de la Fuente M, Gesteiro N, Santiago R, Malvar RA and Butrón A (2022) Genomics and Pathways Involved in Maize Resistance to Fusarium Ear Rot and Kernel Contamination With Fumonisins. Front. Plant Sci. 13:866478. doi: 10.3389/fpls.2022.866478
Received: 31 January 2022; Accepted: 25 March 2022;
Published: 02 May 2022.
Edited by:Thomas Miedaner, University of Hohenheim, Germany
Reviewed by:Alessandra Lanubile, Catholic University of the Sacred Heart, Italy
Daniel Presello, Estación Experimental Pergamino, Instituto Nacional de Tecnología Agropecuaria (INTA), Argentina
Mingliang Xu, China Agricultural University, China
Copyright © 2022 Cao, de la Fuente, Gesteiro, Santiago, Malvar and Butrón. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ana Butrón, email@example.com
†These authors share first authorship