- 1Aquatic and Crop Resources Development, National Research Council of Canada, Saskatoon, SK, Canada
- 2Department of Plant Science, University of Manitoba, Winnipeg, MB, Canada
- 3Ottawa Research and Development Centre, Agriculture and Agri-Food Canada, Ottawa, ON, Canada
- 4Aquatic and Crop Resources Development, National Research Council of Canada, Ottawa, ON, Canada
- 5Swift Current Research and Development Centre, Agriculture and Agri-Food Canada, Swift Current, SK, Canada
- 6Lethbridge Research and Development Centre, Agriculture and Agri-Food Canada, Lethbridge, AB, Canada
Fusarium head blight (FHB) is a devastating fungal disease of small-grain cereals that results in severe yield and quality losses. FHB resistance is controlled by resistance components including incidence, field severity, visual rating index, Fusarium damaged kernels (FDKs), and the accumulation of the mycotoxin deoxynivalenol (DON). Resistance conferred by each of these components is partial and must be combined to achieve resistance sufficient to protect wheat from yield losses. In this study, two biparental mapping populations were analyzed in Canadian FHB nurseries and quantitative trait loci (QTL) mapped for the traits listed above. Nine genomic loci, on 2AS, 2BS, 3BS, 4AS, 4AL, 4BS, 5AS, 5AL, and 5BL, were enriched for the majority of the QTL controlling FHB resistance. The previously validated FHB resistance QTL on 3BS and 5AS affected resistance to severity, FDK, and DON in these populations. The remaining seven genomic loci colocalize with flowering time and/or plant height QTL. The QTL on 4B was a major contributor to all field resistance traits and plant height in the field. QTL on 4AL showed contrasting effects for FHB resistance between Eastern and Western Canada, indicating a local adapted resistance to FHB. In addition, we also found that the 2AS QTL contributed a major effect for DON, and the 2BS for FDK, while the 5AL conferred mainly effect for both FDK/DON. Results presented here provide insight into the genetic architecture underlying these resistant components and insight into how FHB resistance in wheat is controlled by a complex network of interactions between genes controlling flowering time, plant height, local adaption, and FHB resistance components.
Introduction
Fusarium head blight (FHB), or scab, mainly caused by the fungus Fusarium graminearum Schwabe [telomorph: Gibberella zeae Schw. (Petch)], is the most serious fungal disease affecting bread and durum wheat production in Canada (Gilbert and Tekauz, 2000). Fusarium damaged kernels (FDKs) are typically shriveled (or shrunken) and white or pink in appearance (Gilbert and Tekauz, 2000; Bai and Shaner, 2004). The lighter weight of FDKs and the high level of FDKs within a commercial wheat crop result in severe yield and quality losses (Gilbert and Tekauz, 2000; McMullen et al., 2012). FHB is also a food and feed safety concern due to the contamination of grain by the mycotoxin deoxynivalenol (DON) (Gilbert and Tekauz, 2000). The disease has become a more serious threat to wheat production for farmers with trends toward Fusarium isolates that produce higher levels of DON and more frequent FHB outbreaks across the Western Canadian Prairies (Gilbert and Tekauz, 2000). The warm and moist weather conditions experienced in 2014 and 2016, which favor Fusarium infection, caused huge losses for Saskatchewan farmers, with >50% of seed samples infected with Fusarium1. Annual losses attributed to FHB in Canada are in the hundreds of millions of dollars (Haile et al., 2019).
Chemical fungicides and agronomic practices have proven to be only partially effective at controlling FHB (Gilbert and Tekauz, 2000; McMullen et al., 2012), with genetic resistance offering the most effective approach for limiting the economic and ecological impacts of the disease (Gilbert and Tekauz, 2000; Bai and Shaner, 2004; McMullen et al., 2012; Gilbert and Haber, 2013). The genetics of FHB resistance is complex and contains multiple components. The resistance to the initial plant infection measured by the incidence of infection in the presence of natural or augmented inoculum (e.g., spray inoculation) is referred to as type I resistance (Mesterhazy, 1995; Mesterhazy et al., 1999; Bai and Shaner, 2004). Resistance to fungal spread across the wheat head and measured by the severity of infection is referred to as type II resistance (Mesterhazy, 1995; Mesterhazy et al., 1999; Bai and Shaner, 2004). Wheat can also display different levels of resistance to kernel infection, tolerance to infection, and accumulation of DON toxin (Miller et al., 1985; Mesterhazy, 1995; Mesterhazy et al., 1999; Bai and Shaner, 2004). Previous studies indicate that the genetic architecture of these different components are either only partially shared, or independent (Mesterhazy, 1995; Steiner et al., 2004; Lv et al., 2014). Therefore, it can be hypothesized that, in general, combining different types of resistance will result in greater overall resistance to FHB. This could be explored to improve the classification of new wheat cultivars from moderate resistance to a full resistance.
Two main mechanisms can explain the FHB resistance in wheat: physiological resistance (active) and disease escape (passive) (Mesterhazy, 1995). Physiological resistance decreases infection and reproduction of F. graminearum by active plant processes with gene products that contribute to plant defense; disease escape (or passive resistance) allows plants to avoid fungal infection and/or disease progression through morphological and developmental features (Mesterhazy, 1995; Prat et al., 2017). Numerous FHB resistance loci, identified as quantitative trait loci (QTL), are associated with developmental traits, such as heading date, flowering time, anther retention, plant height, and features of the wheat head that include spike length and spike density (Gervais et al., 2003; Srinivasachary et al., 2008, 2009; Buerstmayr et al., 2009; Skinnes et al., 2010; Lu et al., 2013; Szabo-Hever et al., 2014; Buerstmayr and Buerstmayr, 2016). These developmental traits contribute to FHB resistance through disease escape. A notable example is the well-documented negative association between plant height and FHB incidence and severity (Srinivasachary et al., 2008, 2009; Lu et al., 2013; Buerstmayr and Buerstmayr, 2016; He et al., 2016).
Progress in breeding FHB-resistant cultivars has been hindered by the lack of effective resistance sources and the quantitative nature of FHB resistance. FHB resistance is a polygenic trait controlled by multiple genes that can have either major or minor effects and is significantly affected by the interaction of genotype and environment (Bai and Shaner, 1994; Gilbert and Tekauz, 2000; McMullen et al., 2012; Buerstmayr et al., 2019). For hexaploid wheat, several highly resistant sources have been identified and exploited for breeding purposes, including the Chinese cultivars Sumai 3 and Wuhan 1, the Chinese landrace Wangshuibai, the Japanese landrace NyuBai, and the Brazilian cultivar Frontana (Bai and Shaner, 1994; Buerstmayr et al., 2019). However, use of alien germplasm is problematic, notably due to the introgression of undesired genes linked to resistance, known as linkage drag. In these cases, a main challenge is to achieve desirable FHB resistance while maintaining good yield, quality, and agronomic traits. The use of native QTL identified in adapted cultivars without associated linkage drag is therefore preferred by breeders to develop cultivars with better FHB resistance (McCartney et al., 2016; Steiner et al., 2017). Due to the smaller effects of known native QTL, multiple QTL need to be combined, or pyramided, to achieve a desirable level of resistance.
Genetic studies have identified as many as 556 QTL for FHB resistance on all 21 chromosomes in wheat (Venske et al., 2019). Meta-QTL analysis has consolidated these into 19–65 clusters (Liu et al., 2009; Löffler et al., 2009; Venske et al., 2019). Three of these, all originating from the very resistant cultivar Sumai 3, have been validated as major QTL contributing to FHB resistance (Buerstmayr et al., 2009, 2019), including Fhb1 on chromosome arm 3BS (Anderson et al., 2001; Liu et al., 2006), Qfhs.ifa-5A on 5AS (Fhb5) (Buerstmayr et al., 2002; Somers et al., 2003; Xue et al., 2011; Steiner et al., 2019), and Fhb2 on 6BS (Anderson et al., 2001; Cuthbert et al., 2007). All three of these major QTL have been exploited for wheat breeding, although most attention has focused on Fhb1, which is now present in several new FHB-resistant North American and European varieties (Hao et al., 2019). Fhb1 accounts for 20–60% of phenotypic variation in breeding populations, primarily by conferring strong type II resistance (Miedaner and Korzun, 2012). It is the QTL most consistently reported in the literature and three groups have recently cloned Fhb1 gene candidates, although they propose three different mechanisms of action for two different candidate genes (Rawat et al., 2016; Li et al., 2019; Su et al., 2019). The majority of FHB studies have focused on type I and II resistance, with relatively few studies of DON and FDK. Although toxin levels are ultimately the most important measure of disease to consumers, these studies have been limited largely due to the high costs of toxin analysis. The correlation between DON, FDK, and other types of FHB resistance can vary significantly and only a small number of QTL have been reported as exclusively associated with DON, including on 2AS (Semagn et al., 2007), 3A, 7B (Szabo-Hever et al., 2014), 3B, and 3D (He et al., 2019).
In the present study, we aimed to characterize and identify the genetic architecture of different FHB-resistant components including incidence, severity, FDK, and DON from a multi-parental population created by crossing FL62R1 with the two elite Canadian wheat cultivars Stettler and Muchmore. FL62R1 is a wheat line with good FHB resistance derived from the four-way cross QG22.24 / Alsen // SS Blomidon / Alsen that targeted to combine both type I and type II resistance from QG22.24 and Alsen, respectively (Comeau et al., 2008). The line FL62R1 was selected by a “systemic genetic breeding approach” implemented in Eastern Canada (Comeau et al., 2008). It has comparable FHB resistance to the best check, Sumai 3, but possesses good yield potential and agronomic traits (Comeau et al., 2008). Characterization of the resistant QTL for different resistance components from FL62R1 will facilitate achieving desirable FHB resistance by combining these resistance components. In addition, recent findings indicate that FHB resistance QTL are often associated with morphology and development-related genes, such as plant height and flower time (Mesterhazy, 1995; Prat et al., 2017; reviewed by Steiner et al., 2017). In the present study, we further investigated this correlation in different genetic backgrounds and multiple environments. Information from this analysis may provide an approach to incorporate traits related to disease escape into the management of FHB resistance by breeding practices.
Materials and Methods
Materials
Two spring wheat (Triticum aestivum L.) double haploid (DH) populations, derived from the cross of FL62R1 (as common parental line) with Stettler (Stettler population) and Muchmore (Muchmore population), were developed using microspore culture (Eudes and Amundsen, 2005). FL62R1 is an Eastern Canadian spring wheat line derived from the systemic breeding approach, where it was selected from the four-way cross QG22.24 / Alsen // SS Blomidon / Alsen that was subjected to complex stresses over for a few generations (Comeau et al., 2008). The parent Blomidon with pedigree “Weih23.1/Kokart” was bred in Germany, and the QG22.24 parent with the pedigree “5thLACOS-167 / Ae. kotschyi, 400008 // AC Pollet /3/ 5thLACOS-167,” from eastern Canada, showed type I FHB resistance (personal communication with wheat breeders, Drs. A. Comeau and S. B. Rosa). Stettler and Muchmore are two semi-dwarf, high-yielding Canada Western Red Spring (CWRS) wheat cultivars (DePauw et al., 2009, 2011a). The Stettler population contained 182 lines and the Muchmore population consisted of 202 DH lines.
Disease Inoculation and Phenotyping
The two DH populations were evaluated for FHB resistance in disease nurseries at Ottawa, Ontario, and Carman, Manitoba in 2015 and 2016 with three biological replications and a random complete block design (RCBD) in single meter rows. At the Ottawa nursery, grain spawn inoculation was used, as described by Xue et al. (2006) and McCartney et al. (2016). Briefly, a mixture of three local F. graminearum isolates (Xue et al., 2006; McCartney et al., 2016) were used to infect corn and barley kernels, which were spread on the soil surface at a rate of 80 g m2 6 weeks after planting. To promote FHB infection, plots were irrigated for 30 min every morning and 30 min each afternoon to create favorable moisture for infection (Xue et al., 2006; McCartney et al., 2016).
A spray approach was used to inoculate plants at the Carman disease nursery. When ∼50% of plants reached anthesis, a mixture of four local F. graminearum isolates were applied at a rate of 50 ml per row and sprayed onto plants at a concentration of 50,000 macroconidia L–1. Plants were sprayed a second time 2–3 days after the first application. Thereafter, a daily mist irrigation for 10 min every hour for a 12-h period (6:00 pm–6:00 am) was utilized to promote FHB infection (McCartney et al., 2016).
At 18–21 days post-inoculation, FHB incidence (INC) was scored as an estimate of the percentage of infected heads within the plot, and severity (SEV) was measured as an estimate of infected spikelets in an infected head. Visual rating index (VRI) was estimated from INC and SEV using the formula: VRI = (INC × SEV)/100. Plants were harvested with minimum air force to prevent the loss of low weight infected kernels and a visual assessment of kernel damage in the collected seed was used to estimate kernel damage, expressed as percentage of FDK. The concentration of DON in grain was evaluated by commercial enzyme linked immunosorbent assay (ELISA) kits, Neogen Veratox® (Lansing, MI, United States) for DON 5/5 ELISA kits (Ye, 2015). In addition, measurements of days to anthesis (DA) was recorded when 50% of the main tillers in the row had begun anthesis and plant height (HT) was measured from the soil surface to the top of main tiller spikes.
Genotyping
Genomic DNA was extracted from the two DH populations using the BioSprint 96 Extraction Platform and DNA Plant Kit (Qiagen, Hilden, Germany). DNA was resuspended in purified water and quantified with a Quant-iTTM PicoGreen® dsDNA Assay Kit (Thermo Fisher Scientific Inc., Bartlesville, OK, United States). All DNA samples were diluted to 50 ng/μl for SNP array genotyping. The two DH populations were genotyped with the Illumina iSelect 90K SNP array (Wang et al., 2014) and raw data were processed and cleaned as described by Wang et al. (2014) using the diploid version of GenomeStudio (Illumina, San Diego, CA, United States). The SNP marker to differentiate Rht-B1 semi dwarf allele was run following standard KASP guidelines (Rasheed et al., 2016). Additionally, 423 SSR markers were tested, including some previously reported as associated with FHB QTLs and randomly selected markers that spanned every chromosome. Markers that showed clear polymorphism and segregation were selected and added to the SNP map. High-quality SNP and SSR markers were further refined with the R/QTL package (Broman et al., 2003) using the following criteria: heterozygous SNPs were converted to missing data; SNPs with missing genotypic data higher than 15% and/or individuals with more than 20% missing genotypic data were removed; SNPs and SSRs with the minor allele frequency less than 10% were also discarded from downstream analysis.
QTL Mapping
Genetic maps were developed with Mapdisto (Lorieux, 2012) using the markers described above, a cutoff recombination value of 0.3, and a threshold logarithm of odds (LOD) of 6.5. Physical positions of mapped markers were retrieved by blast searching SNP sequences against IWGSC RefSeq v1.0 (International Wheat Genome Sequencing Consortium, 2018; CS Ref 1.0) with best hits. QTL analysis was performed by composite interval mapping (CIM) as implemented in WinQTL Cartographer (Wang et al., 2007) with a walking speed of 1 cM and a window size of 10 cM. The significance of a QTL was claimed at the threshold of 5% significance level by the 3000 permutation test in WinQTL Cartographer (Churchill and Doerge, 1994; Basten et al., 2001; Wang et al., 2007). Due to the known minor effects of most FHB resistance-related QTL, QTL with an LOD value below the 5% significance level but higher than the empirical cutoff threshold value of LOD = 2.5 were also included in the analysis. QTL were physically projected on CS Ref 1.0 (International Wheat Genome Sequencing Consortium, 2018) using the physical position of peak markers and were visualized by the R package RIdeogram (Hao et al., 2020).
Statistical Analysis
Statistical analysis of the available phenotypic data was performed in R 3.30 (R Core Team, 2014) with the lme4 package (Bates et al., 2015). A linear mixed model analysis (Bates et al., 2015) was fitted for each phenotypic trait. For within-year data analysis, the linear model used was yib = μ + Gi + Rb + eib, where yib is the trait, Gi is the effect of the i-th genotype in the b-th block, μ is the mean value, Rb is the effect of the b-th block, and eib is the residual. For across years, the model used was yijb = μ + Gi + Ej + (GXE)ij + r(j)b + eijb, where, yijb is the phenotypic value of i-th genotype in the k-th year and b-th block, μ is the mean value, Gi is the effect of the i-th genotype, Ej is the effect of the j-th year, (GXE)ij is the effect genotype-by-year interaction of the i-th genotype in j-th year, r(j)b is the effect of b-th block in the j-th year, and eijb is the residual. The variance component of each trait was determined by the restricted maximum likelihood (REML) method within lme4 (Bates et al., 2015) by setting all effects as random. Variance components were used to estimate the repeatability within year with the equation and broad-sense heritability across year was estimated using the equation , where , , and stand for variance of genotype, genotype-by-year interactions, and error; j represents number of years, and b indicates number of blocks.
Results
Phenotypic Analysis
Considering all lines from the two populations analyzed, INC ranged from 0 to 100% at both Carman and Ottawa with a mean value of approximately 70–80% (Table 1 and Figure 1). The distributed SEV observed at Carman ranged from 5 to 95%, skewed toward a lower mean value of around 30%, while at Ottawa, SEV had a near-normal distribution, with a mean around 45%. FDK was observed from 1 to 47% with a mean of approximately 15% at both sites (Table 1 and Figure 1). DON levels ranged from 0.1 to 47.0 ppm with distribution toward a lower mean value of about 10 ppm at both sites (Table 1 and Figure 1). DA and HT have high heritability and show normal distribution (Table 1 and Supplementary Figure 1). All traits tested showed moderate to high heritability, with H2 values from 0.55 to 0.88 (Table 1).
Table 1. Phenotypic variation and heritability (H2) of Fusarium head blight (FHB) and plant traits from field experiments of the Muchmore and Stettler populations in different environments
Figure 1. Distribution of phenotypic traits from field trials in different years and locations. Fusarium damaged kernels (FDKs), incidence (INC), and severity (SEV) at (A) Carman, Manitoba (W) and (B) Ottawa, Ontario (O) trials; and (C) deoxynivalenol (DON) from all environments.
Positive correlation was observed among all FHB resistance components, while DA and HT showed mostly negative correlation with FHB traits (Figure 2). When comparing the same location and year, INC, FDK, SEV, and DON correlated well with each other (r = 0.4–0.7). Correlations were still high among these four resistant components when comparing across 2 years at the same location, with r ranging from 0.3 to 0.7. However, the correlation between FHB resistance components was lower across sites, with moderate r ranging from 0.3 to 0.5 across all environments. When comparing the individual FHB measures, the visual disease resistance components (INC, SEV, and VRI) were generally better correlated with each other than with FDK and DON, and FDK and DON were very highly correlated with each other (Figure 2). A high negative correlation was observed between HT and all the resistance components at Ottawa in both years, ranging from −0.3 to −0.6. However, a substantially lower correlation was found for these traits at Carman (from 0 to 0.3). DA displayed the lowest correlation against the resistance components, with r values from −0.2 to 0.4. The HT and DA traits showed the highest correlation across all the tested environments with r values from 0.5 to 0.7, respectively.
Figure 2. Phenotypic correlation between Fusarium head blight-resistance components and associated traits at two field locations for two populations. The areas of the circles show the absolute value of the corresponding correlation coefficients (r). INC, incidence; SEV, severity; FDK, Fusarium damaged kernels; DON, deoxynivalenol; HT, plant height; DA, day to anthesis; O, Ottawa, Ontario; W, Carman, Manitoba. 15, 2015; 16, 2016; ST, Stettler population; MM, Muchmore population.
The correlation between greenhouse type II resistance in the Stettler population (Zhang et al., 2018) and field resistance was also examined (Supplementary Figure 2). The type II FHB resistance measured in the greenhouse showed significant but lower correlations with SEV in all tested environments (r = 0.23–0.37; Supplementary Figure 2). At Carman, except DON in 2015, INC, FDK, and DON also showed significant but lower correlations with greenhouse resistance (Supplementary Figures 2A,C). However, there was no correlation between the greenhouse results and Ottawa trials except for a small but significant correlation with 2015 FDK and DON (Supplementary Figures 2B,D).
QTL Analysis
Quantitative trait loci analysis identified a large number of loci that contributed the different resistance components (Table 2 and Supplementary Table 1). Using the IWGS Chinese Spring (CS) reference 1.0 (CS Ref 1.0; International Wheat Genome Sequencing Consortium, 2018), the physical position of the QTL was anchored to approximate genomic regions. This permitted us to reduce the large number of QTL down to nine genomic loci affecting all of the FHB resistance components, located on chromosome arms 2AS, 2BS, 3BS, 4AS, 4AL, 4BS, 5AS, 5AL, and 5B (Table 2 and Figure 3). The majority of the FHB resistance QTL were contributed by the common parental line FL62R1, and generally, this line was also associated with increased HT and delayed DA. Detailed information on QTL detected is provided below, prioritized by their contribution to FHB resistance.
Table 2. Nine major QTL identified for FHB-resistant components across four environments in the Muchmore and Stettler populations.
Figure 3. QTLs of Fusarium head blight resistance, flower time, and plant height anchored on the Chinese Spring reference 1.0 (CS Ref 1.0). The marker located at the peak of the QTL identified by blast on CS Ref 1.0 was used to localize QTL. INC, incidence; SEV, severity; FDK, Fusarium damaged kernels; DON, deoxynivalenol; HT, plant height; DA, day to anthesis; O, Ottawa, Ontario; W, Carman, Manitoba.
4BS: The 4BS genomic locus, from 7.04 to 42.02 Mb on the CS Ref 1.0 assembly, is the region with the highest number of QTL identified in this study. It was detected across almost all environments for the FHB resistance traits and plant height (Table 2 and Figure 3). The phenotypic effect of the 4BS region on the different resistance components was greater than observed for other QTL identified, with the percentage of phenotypic variance explained by the QTL (PV) ranging from 5.4 to 28% (Table 2). The 4BS locus had a more stable effect in the Stettler population compared to the Muchmore population. We also observed that 4BS contributed more to FDK and DON resistance at Ottawa than at Carman, with approximately twofold higher PV in both populations in 2 years of field trials. At Ottawa in 2015, 4BS QTL were found to have the largest effect on all FHB resistance components, with PV of 35.4 and 27.4% for FDK and DON, respectively. Mapping of the functional Rht-B1 KASP marker (Rasheed et al., 2016) and anchoring to the CS Ref 1.0 assembly located the 4BS QTL to the chromosomal region carrying Rht-B1.
3BS: QTL identified on 3BS were localized to an interval between 6.68 and 11.58 Mb on the CS Ref 1.0 assembly. QTL for SEV, FDK, and DON were stably detected in both populations at Carman in 2015 with relatively large effects, with PV of 7–19%. During the 2016 trial, QTL also were found for these components except for SEV in the Stettler population. No QTL was found for INC in this genomic locus. QTL for VRI were identified in the Stettler population in 2015 and the Muchmore population in 2016. At Ottawa, the only QTL in this region was identified for DON in the 2015 Stettler population. No DA or HT QTL were found in the 3BS region.
4AL and 4AS: The QTL from 4AL (physical position from 641.4 to 692.4 Mb on the CS assembly) were mainly detected in the 2015 trials. For the Muchmore population, stable QTL were found at both Carman and Ottawa for SEV, FDK, DON, and VRI, while QTL for INC were only detected at Carman in this region. For the Stettler population, QTL were also found for SEV, FDK, DON, and VRI on 4AL at Ottawa, while this region only had a QTL for DON at Carman. In 2016 the 4AL region had a smaller impact on FHB resistance, with only a minor effect for DON QTL found in the Muchmore population at Carman. One of the largest and most significant QTL in the entire dataset was found for DA in the 4AL region. It was detected in each population, year, and location, except for Stettler at Ottawa in 2016, and had a large impact with PV ranging from 10.9 to 49.8%. The common parental line, FL62R1, contributed to longer DA and taller plants, though at Ottawa, the locus contributed to FHB resistance but increased FHB susceptibility at Carman.
In the region on 4AS (∼12.2–38.7 Mb), QTL were only detected in the Stettler population. SEV was the most stable QTL detected at Ottawa and Carman. The alleles contributed by Stettler conferred better disease resistance at the 4AS QTL.
5AL and 5AS: The 5AL genomic locus covered the region from 528.6 to 584.7 Mb on the CS Ref 1.0 assembly. QTL from this region were observed for all of the traits, though each was not consistent across environments. The largest effect 5AL QTL were for DON at Ottawa; the Muchmore population in 2016 with a PV of 16%, followed by a 9.8% PV from the Stettler population in 2015. Minor QTL for DA and HT were found on the 5AL region in both populations.
FHB resistance QTL were found on 5AS in the Muchmore population only. In 2015, 5AS QTL were found in Ottawa and Carman (contributing up to 11% PV). In 2016, a 5AS QTL for DON was detected in Carman, in addition to an FDK QTL at Ottawa. VRI and SEV 5AS QTL were found at both sites in 2016.
5BL: QTL from this region were located from 551.0 to 582.5 Mb in the CS 1.0 reference genome. The 5BL QTL were only identified in the Stettler population, where they were responsible for a large phenotypic effect on FHB resistance at Ottawa in 2015 and 2016, with PV as high as 13.3% for INC and 22.7% for SEV. In all cases, alleles from Stettler contributed to disease resistance at the 5BL locus. A major QTL with large effect for DA was detected in the 5BL interval in the Stettler population, with Stettler responsible for delayed anthesis.
2AS and 2BS: The 2AS QTL were localized from 36.6 to 58.4 Mb on the CS Ref 1.0 assembly. This region conferred resistance for almost all of the resistance components in both the Stettler and Muchmore populations in both years of testing at the Carman site. No QTL was detected on 2AS at Ottawa.
The 2BS QTL region controlled the majority of FHB resistance traits in the Muchmore population at Ottawa in 2015 and 2016.
Discussion
FHB resistance in wheat is a complex trait and controlled by several resistance components including INC, SEV, DON accumulation, FDK, and spreading (greenhouse). A better understanding of the genetic architecture underlying these different components may enable the development of wheat cultivars with better resistance to FHB. In the present study, a systematic analysis of different components was undertaken with two large bi-parental mapping populations, revealing how different FHB resistance components, plant height and flowering time, interact.
Correlation Between FHB Resistance Component Traits
There was a moderate to high heritability displayed across the FHB resistance components tested, which is in agreement with previous findings, indicating that genetic variation plays a main role in the phenotypic variation of these FHB traits (Buerstmayr et al., 2000; Lu et al., 2013; He et al., 2019). In contrast, the highly complex genotype-by-environment interactions of FHB resistance are emphasized with the low to moderate correlation observed between the different FHB resistance components measured in this study (INC, SEV, FDK, and DON) when comparing across years and sites. Visual disease symptoms were highly correlated, as were FDK and DON; however, INC and SEV were found to have lower correlation with DON and FDK (Figure 2A). Our results agree with previous findings (Paul et al., 2005; Lu et al., 2013; He et al., 2019) that visual FHB disease symptoms are poorly correlated with DON content, indicating that visual FHB disease resistance contributes only in part to DON accumulation. Given the increasing priority for breeding new varieties with greater resistance to DON, Canadian cultivar registration is evaluated with an index (ISD = 0.2 × INC + 0.2 × SEV + 0.6 × DON) that gives greater weight to DON accumulation than visual measures of disease resistance combined2. There is a clear need to understand the genetic architecture of the different resistance components beyond INC and SEV.
FHB type II resistance (resistance to fungal spread in the spike) is determined by point inoculation in the greenhouse, while SEV is routinely used as an estimate for field FHB spread/type II resistance (Bai and Shaner, 2004). There was a significant but low correlation between greenhouse type II resistance using our previously published data (Zhang et al., 2018) and SEV in the present study (Supplementary Figure 2). Lu et al. (2013) observed a lack of correlation between greenhouse and field FHB resistance in their germplasm, reasoning that this discrepancy is due to a lack of Fhb1, a gene commonly reported as conferring major type II resistance. FL62R1 is moderately susceptible to point inoculation despite carrying Fhb1 (Zhang et al., 2018), supporting the interpretation from Lu et al. (2013). However, the major 2BL QTL controlling greenhouse type II resistance identified by Zhang et al. (2018) was not detected by field evaluation, indicating that greenhouse type II resistance may occur through a somewhat different mechanism than field resistance in the populations analyzed in this study.
A higher negative correlation was observed between plant height and FHB resistance components at Ottawa than at Carman. This is likely the result of the different FHB inoculation approaches utilized at each location: grain spawn at Ottawa and spray treatment at Carmen. The spawn approach relies on wind or water splashes to spread FHB spores from the ground surface to spikes and establish infection (McMullen et al., 2012). The spray approach directly applies FHB spores onto wheat spikes, thus minimizing the effect of height on FHB infection caused by the different disease pressure at different distances between the spike and soil surface. This may explain the low negative correlation between plant height and FHB resistance components at Carman. However, taller plants are also prone to have a less favorable environment (less moisture) for disease infection and development (Yan et al., 2011), and neither spawn nor spray approach can differentiate the effect of height on FHB resistance related to the different spike microclimates. The strong negative correlation between FHB resistance and plant height and flowering time is well established (Buerstmayr et al., 2009; Buerstmayr and Buerstmayr, 2016) and discussed in greater detail below.
Genetic Architecture of FHB Resistance Components and Their Association With Flowering Time and Plant Height
A large number of loci that control different FHB resistance were identified, with genomic loci on 2AS, 2BS, 3BS, 4AS, 4AL, 4BS, 5AS, 5AL, and 5BL enriched with QTL for one or more traits. The genomic locus 4BS from FL62R1 conferred the most stable FHB resistance, with INC, SEV, FDK, and DON QTL identified at all tested environments, in both the Stettler and Muchmore populations. A major QTL for plant height was also identified in this region of 4BS that can be attributed to the presence of the semi-dwarf allele of Rht-B1. Numerous studies have previously demonstrated that semi-dwarf Rht1 alleles are associated with FHB susceptibility (Srinivasachary et al., 2008, 2009; Lu et al., 2013; Buerstmayr and Buerstmayr, 2016; Prat et al., 2017), with the wild-type alleles Rht-B1a and Rht-D1a contributing to FHB resistance. Rht-B1 may contribute to resistance directly through the control of the plant height and flower microclimate through pleiotropy, hormone signaling, and modulation of flowering, anther retention and spike traits, or even through linkage drag (McCartney et al., 2007; Yan et al., 2011; Buerstmayr and Buerstmayr, 2016; Davière and Achard, 2016; He et al., 2016; Herter et al., 2018; Xu et al., 2020). The present study found that when the effect of plant height on FHB resistance caused by different disease pressure was controlled through the use of spray inoculation at Carman, the 4BS locus still had a large effect on FHB resistance, though the effect was reduced compared to the site using the spawn approach. The effect caused by different microclimates of spikes at different heights was not investigated in this study, so height differences per se could still be responsible for the FHB resistance observed, though other mechanisms are likely also involved. Future application of CRISPR technology to edit wild-type Rht-B1 alleles will likely prove very useful to resolve the association between FHB resistance and Rht-B1 dwarf alleles.
The largest effect and most extensively characterized FHB resistance QTL is Fhb1 located on chromosome arm 3BS (Anderson et al., 2001; Buerstmayr et al., 2019; Hao et al., 2019), which was also identified in this study. We found that the Fhb1 QTL from the common parental line FL62R1 affected SEV, VRI, FDK, and DON, with a large phenotypic effect on FDK and DON at Carman in both populations. Fhb1 was likely contributed from Alsen, one of the parental lines of FL62R1 (Comeau et al., 2008). Alsen was the first North American high-quality hard red spring cultivar with improved FHB resistance (Frohberg et al., 2006). Fhb1 has been proposed as the key factor for DON detoxification (Lemmens et al., 2005), and our results support this claim. In spite of Fhb1 being important to type II FHB resistance, our previous study (Zhang et al., 2018) found that Fhb1 from FL62R1 had no effect in the Stettler population during greenhouse point inoculation tests. This may be attributed to Fhb1 having distinct effects when in different genetic backgrounds (McCartney et al., 2007; Bokore et al., 2017; Li et al., 2019; Brar et al., 2019). Most research has assessed the effect of Fhb1 on non-adapted or adapted germplasm, with few studies conducted on elite cultivars. Indeed, testing of additional Canadian elite cultivars, such as Carberry (DePauw et al., 2011b) and Penhold (Cuthbert et al., 2017), corroborates the moderate susceptibility demonstrated in greenhouse point inoculation tests (L. P. Wang, unpublished). However, it is also possible that the heavy inoculum used in our greenhouse inoculation tests overwhelmed the plants and prevented the detection of any Fhb1 effect.
In addition to genetic background, the method of inoculation, the aggressiveness of the F. graminearum race tested, the presence of other resistance genes, and environmental conditions are factors that could affect the level of resistance conferred by Fhb1. As pointed out by Hao et al. (2019), the different stages of development for inoculation and the aggressiveness of races used in Fhb1 cloning experiments (Rawat et al., 2016; Li et al., 2019; Su et al., 2019) may be partially responsible for the different findings of these studies. The present study applied the spawn method of inoculation at Ottawa and spray inoculation at Carman. The aggressiveness of the races tested at each location differed, and the Fhb1 QTL was only detected at Carman. It is possible that the stronger effect of the 4BS QTL at Ottawa, associated to plant height and thus more influenced by the inoculation method used, masked the detection of Fhb1. Liu et al. (2019) reported that Fhb1 had a very minor effect in a hard red spring wheat breeding population in disease nurseries that used the spawn approach. Although these authors did not compare the effectiveness of Fhb1 using spray inoculation, their findings are consistent with results we observed at Ottawa with spawn inoculation. Finally, Bokore et al. (2017) reported different effectiveness of Fhb1 in Sumai 3-derived North American spring wheat breeding lines when tested at two locations.
The 5AS QTL was only found in the Muchmore population and showed variable effects on different resistance components. As observed for 3BS (Fhb1), no height- or flowering time-related QTL were detected in this region. The 5AS QTL harbors the SSR marker GWM304 that is linked to Fhb5, indicating that it may be Qfhs.ifa-5A (Fhb5). Fhb5 has been reported to contribute mainly to INC and SEV when tested using the spray approach (Buerstmayr et al., 2002; Somers et al., 2003; Steiner et al., 2019). In the present study, the 5AS QTL showed minor effects for SEV, FDK, and DON in a portion of years and locations tested. The 5AS resistance allele is likely derived from Alsen through FL62R1. The Canadian CWRS cultivar Carberry has notably higher FHB resistance for its class, and also carries the Fhb5 allele derived from Alsen (Bokore et al., 2017). These findings indicate that in addition to Fhb1, Fhb5 was also deployed into some North American wheat cultivars. Because of its relatively minor effect, retaining this allele through visual section will be challenging, and marker-assisted selection could be a useful approach to introduce it to new cultivars and breeding programs.
The remaining genomic loci controlling FHB resistance (2AS, 2BS, 4AS, 4AL, 5AL, and 5BL) identified in this study were all found to colocalize with flowering time or flowering time and plant height. The large number of FHB resistance loci associated with flowering time indicates that a complex flowering regulatory network was involved in controlling disease resistance in the populations tested. The establishment of FHB infection needs a favorable environment for disease establishment and the appropriate susceptible stage of plant development, namely, the early flowering stage at the time of fungal introduction (Mesterhazy, 1995). Thus, if wheat plants do not flower during the period that offers favorable conditions for FHB infection created by the local environment, they will show reduced disease symptoms. Flowering time is largely thought of as controlling FHB infection through disease escape. However, it has recently become clear that the flowering period is also a key factor for wheat spikelet development, with many pleiotropic genes/QTL influencing the timing of flowering and development of spikes (Lewis et al., 2008; Shaw et al., 2013; Liu et al., 2019; Wolde et al., 2019; Chen et al., 2020). This shows that the influence of flowering time on FHB resistance is likely to be much more complex than previously expected.
The chromosome arms 4AL and 5BL contain two major genomic loci that showed the largest effect for DA and minor effects for HT. The 4AL QTL were found to control SEV, FDK, and DON, with the largest effect on SEV. QTL from 4A have previously been reported in durum wheat (Prat et al., 2017) and in the winter wheat cultivars Arina (Paillard et al., 2004; Lu et al., 2013; Buerstmayr and Buerstmayr, 2015) and Heyne (Zhang et al., 2012). Recently, Wu et al. (2019) identified a major QTL from the same region for type II resistance (point inoculation in the greenhouse). The previous study of FL62R1 found a minor QTL from the same region that controls resistance for point inoculation (Zhang et al., 2018). Together, results indicate that the 4AL locus plays roles for both greenhouse and field resistance. The 4AL region has been characterized intensively for structure variation with the presence of two reciprocal translocations and two inversions (Mickelson-Young et al., 1995; Miftahudin et al., 2004; Dvorak et al., 2018) between a portion of 4AL to 4DL, where a portion of 4AL corresponds to 4DL. The present study also found that this region has a paracentric inversion when comparing the genetic order of our map to the physical position of this region in the CS Ref 1.0 assembly. Wolde et al. (2019) mapped the TB-A1, which may have arisen from translocation of TEOSINTE BRANCHED1 (TB1) from the 4DS region to the 4AL region between 629.3 and 634.8 Mb, at almost the same position as we identified the 4AL FHB QTL. TB1 is known to play a role in controlling morphometric inflorescence traits in wheat, and we have previously identified QTL for spike density and flowering time that colocalize to the region (Zhang et al., 2018). In addition, Dixon et al. (2018) found that TB1 interacts with FLOWERING LOCUS T1 to regulate inflorescence architecture in bread wheat (T. aestivum L.). These findings, together with ours, indicate that TB1 may be the candidate gene within this region for FHB resistance. We also found that FHB resistance from this region is intensely affected by environments and sites. The common parental line, FL62R1, contributed to longer DA and taller plants; in Ottawa, the locus contributed to FHB resistance while it increased FHB susceptibility at Carman. This signified that the factors responsible for local adaption also control the FHB resistance mechanism. Flowering time genes are critical to adjusting the timing of wheat flowering to its climate and can have opposite effects on fitness in different environments; thus, the locus expresses a genetic trade-off between environments (Mitchell-Olds, 2013).
The 5BL QTL region has a major effect for DA and a minor effect for height in the Stettler population, with the QTL from Stettler delaying time to anthesis. The major vernalization gene VrnB1 is located in this region. This 5BL region conferred FHB resistance at Ottawa, having larger effects on incidence and SEV than FDK and DON. This QTL was also found in spring wheat and in Tunisian landrace durum wheat (Qfhs.ndsu-5BL, Ghavami et al., 2011; He et al., 2019). Currently, two major FHB resistance loci on 5AL have been reported. The first maps to the VrnA1 region, around 584.7 Mb, and has been described in both hexaploid and durum wheat (Prat et al., 2017; Sari et al., 2018; He et al., 2019). The FHB resistance 5AL QTL from FL62R1 mapped to this region, though it only had a minor effect on flowering time and plant height. The second loci correspond to the Q region identified from Triticum turgidum ssp. dicoccum (Zhang et al., 2014), Triticum macha (Buerstmayr et al., 2011), and the hexaploid PI277012 (Chu et al., 2011), which localized around 658.1 Mb. Q and VrnA1 map within 70 Mb from each other, indicating the existence of two QTL, though large phenotypic bias and relatively small population size make it possible that both QTL correspond to the same locus. Given that the QTL were identified from tetraploid to hexaploid wheat species, it is probable that the resistance gene is highly conserved from an evolutionary perspective. With this assumption, Vrn1 and the Q gene, both responsible for domestication, are likely candidate genes.
We identified FHB resistance QTL for FHB syntenic region on 2AS and 2BS colocalized with the major flower control genes Ppd-A1 and Ppd-B1. There was only a minor effect on DA, and the detection of the QTL was not stable across environments: the 2BS QTL was only found in the Muchmore population, and almost exclusively at Ottawa; conversely, the 2AS QTL was found in both populations, but almost exclusively at Carman. Previous studies have also identified QTL for FHB resistance in these regions (Giancaspro et al., 2016; He et al., 2019), and we previously reported a minor QTL for resistance following greenhouse point inoculation in this region (Zhang et al., 2018). Ppd-1 has a major effect on spikelet formation, and the photoperiod insensitive allele speeds up the rate of spikelet initiation (Ochagavia et al., 2018). The Ppd genes are potential candidate genes for FHB resistance of both the 2AS and 2BS QTL, possibly conferring both greenhouse and field resistance through controlling the speed of flower development.
Finally, intensive research have been conducted for types I and II of FHB resistance in wheat. Despite the importance of FDK for grain quality and DON for food safety, little research was performed on FDK and DON. Recently, He et al. (2019) reported two DON QTL, 3BL, which also conferred a minor effect for FHB resistance, and 3DL, which showed no effect on field FHB resistance. Within the current study, no major QTL were identified exclusively for FDK and/or DON, but QTL with larger effect for FDK and/or DON were observed. The 2AS QTL contributed a major effect for DON and the 2BS QTL showed a major effect for FDK. The 3BS QTL had a larger effect on FDK and DON than other FHB measures, and the 5AL QTL primarily affected FDK and DON. These findings show that resistance genes for FDK and/or DON are partially shared with incidence and severity. Because of this as well as the importance of FDK and DON, there is a need to invest research on FDK and DON resistance in wheat. The QTL we identified here can be used by breeding programs as promising target QTL to develop low FDK and low DON wheat cultivars.
Conclusion and Perspective for Future Genetic Research and FHB Resistance Breeding
The genetic architecture of high-level FHB resistance for different resistance components from FL62R1 were identified. It is concluded that FL62R1 confers Sumai 3 levels of resistance through the interaction of Fhb1, Fhb5, RhtB1a, and a complex flower regulatory network consisting primarily of six genomic loci controlling flower timing, including photoperiod, vernalization, and the TEOSINTE BRANCHED1 gene. Different resistance components are partially controlled by different genetic networks and more research is needed to characterize the underlying genetic architecture. To achieve better overall FHB resistance, it is necessary to pyramid different resistant components, especially for the less well understood resistance to FDK and DON accumulation.
For future research into the genetic control of FHB resistance, there is a need to clearly characterize the important resistant loci, from the magnitude of their effects to the effect on different resistant components, especially in elite cultivar genetic backgrounds. This is especially true for the major resistance gene Fhb1. Given the recent cloning of two genes said to be Fhb1 by three separate groups (Rawat et al., 2016; Li et al., 2019; Su et al., 2019), another need is to clearly characterize genes at this locus and their biological function for FHB resistance. Our findings that disease escape caused by plant height and flowering time can act as the predominant resistance mechanism and even totally mask the Fhb1 effect emphasize the need to differentiate physiological resistance from escape. One approach is to optimize genetic mapping populations through haplotype analysis of parental lines and fixing alleles of key genes controlling flowering time and plant height. Alternatively, a subpopulation that does not segregate for these alleles can be developed by stratifying these loci from a larger mapping population through haplotype analysis and field data as practiced by Ren et al. (2019). Finally, since spray inoculation can minimize the effect of height and flower time, it should be helpful to precisely identify genomic loci contributing to physiological resistance and estimate their effects in the genetic study of FHB resistance.
From a breeding perspective, due to the complexity of the network controlling flowering time, its association with FHB resistance, and local adaption properties, it is important to completely identify the haplotypes of the major flowering time gene(s) in existing breeding programs. By fine-tuning the complex network from local adapted haplotypes (or alleles), and combining Fhb1 and/or Fhb5 as a base for FHB resistance with appropriate semi-dwarf alleles, desirable FHB resistance comparable to Sumai 3 levels should be achievable. The success of FL62R1 is a good example to support this breeding practice.
Data Availability Statement
The datasets generated for this study can be found in FigShare, https://doi.org/10.6084/m9.figshare.12948263.v1.
Author Contributions
PF conceived this study and acquired fund for this study. PF, WZ, KB, and RC designed the experiment. FJ and FE developed DH populations. HR contributed to seed increasing of DH populations. KB, PG, BP, and WZ performed the experiments. AB-B, GF, and ZR contributed to field trials and disease evaluation. WZ and KB analyzed the data and interpreted the results. WZ, KB, and PF wrote the manuscript. All authors contributed to the article and approved the submitted version.
Funding
This research was funded by the Sustainable Agri-Food program of National Research Council Canada.
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.
Acknowledgments
The authors thank Ron DePauw for Stettler and Muchmore germplasm, Janet Condie for technical assistance with the 90K SNP genotyping, and Daoquan Xiang for critically reviewing the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.580833/full#supplementary-material
Footnotes
- ^ https://www.grainscanada.gc.ca/en/grain-research/export-quality/cereals/wheat/western/annual-fusarium-damage/
- ^ http://www.pgdc.ca/
References
Anderson, J. A., Stack, R., Liu, S., Waldron, B., Fjeld, A. D., Coyne, C., et al. (2001). DNA markers for Fusarium head blight resistance QTLs in two wheat populations. Theor. Appl. Genet. 102, 1164–1168. doi: 10.1007/s001220000509
Bai, G., and Shaner, G. (2004). Management and resistance in wheat and barley to Fusarium head blight. Annu. Rev. Phytopathol. 42, 135–161. doi: 10.1146/annurev.phyto.42.040803.140340
Basten, C. J., Weir, B. S., and Zeng, Z. (2001). QTL Cartographer Version 2.0. Raleigh, NC: North Carolina State University.
Bates, D., Mächler, M., Bolker, B., and Walker, S. (2015). Fitting linear mixed effects models using lme4. J. Stat. Softw. 67, 1–48.
Bokore, F. E., Knox, R. E., DePauw, R. M., Clarke, F., Cuthbert, R. D., Campbell, H. L., et al. (2017). Validation of molecular markers for use with adapted sources of Fusarium head blight resistance in wheat. Plant Dis. 101, 1292–1299. doi: 10.1094/pdis-10-16-1421-re
Brar, G. S., Brûlé-Babel, A., Ruan, Y. F., Henriquez, M. A., Pozniak, C. J., Kutcher, H. R., et al. (2019). Genetic factors affecting Fusarium head blight resistance improvement from introgression of exotic Sumai 3 alleles (including Fhb1, Fhb2, and Fhb5) in hard red spring wheat. BMC Plant. Biol. 19:179. doi: 10.1186/s12870-019-1782-2
Broman, K. W., Wu, H., Sen, S., and Churchill, G. A. (2003). R/qtl: QTL mapping in experimental crosses. Bioinformatics 19, 889–890. doi: 10.1093/bioinformatics/btg112
Buerstmayr, H., Ban, T., and Anderson, J. A. (2009). QTL mapping and marker-assisted selection for Fusarium head blight resistance in wheat: a review. Plant Breed. 128, 1–26. doi: 10.1556/crc.36.2008.suppl.b.1
Buerstmayr, H., Lemmens, M., Hartl, L., Doldi, L., Steiner, B., Stierschneider, M., et al. (2002). Molecular mapping of QTLs for Fusarium head blight resistance in spring wheat. I. Resistance to fungal spread (type II resistance). Theor. Appl. Genet. 104, 84–91. doi: 10.1007/s001220200009
Buerstmayr, H., Steiner, B., Lemmens, M., and Ruckenbauer, P. (2000). Resistance to Fusarium head blight in winter wheat: heritability and trait associations. Crop Sci. 40, 1012–1018. doi: 10.2135/cropsci2000.4041012x
Buerstmayr, M., and Buerstmayr, H. (2015). Comparative mapping of quantitative trait loci for Fusarium head blight resistance and anther retention in the winter wheat population Capo x Arina. Theor. Appl. Genet. 128, 1519–1530. doi: 10.1007/s00122-015-2527-8
Buerstmayr, M., and Buerstmayr, H. (2016). The semidwarfing alleles Rht-D1b and Rht-B1b show marked differences in their associations with anther-retention in wheat heads and with Fusarium head blight susceptibility. Phytopathology 106, 1544–1552. doi: 10.1094/phyto-05-16-0200-r
Buerstmayr, M., Lemmens, M., Steiner, B., and Buerstmayr, H. (2011). Advanced backcross QTL mapping of resistance to Fusarium head blight and plant morphological traits in a Triticum macha x T. aestivum population. Theor. Appl. Genet. 123, 293–306. doi: 10.1007/s00122-011-1584-x
Buerstmayr, M., Steiner, B., and Buerstmayr, H. (2019). Breeding for Fusarium head blight resistance in wheat – Progress and challenges. Plant Breed. 139, 429–454. doi: 10.1111/pbr.12797
Chen, Z., Cheng, X., Chai, L., Wang, Z., Du, D., Wang, Z., et al. (2020). Pleiotropic QTL influencing spikelet number and heading date in common wheat (Triticum aestivum L.). Theor. Appl. Genet. 133, 1825–1838. doi: 10.1007/s00122-020-03556-6
Chu, C., Niu, Z., Zhong, S., Chao, S., Friesen, T. L., Halley, S., et al. (2011). Identification and molecular mapping of two QTLs with major effects for resistance to Fusarium head blight in wheat. Theor. Appl. Genet. 123, 1107–1119. doi: 10.1007/s00122-011-1652-2
Churchill, G. A., and Doerge, R. W. (1994). Empirical threshold values for quantitative trait mapping. Genetics 138, 963–971.
Comeau, A., Langevin, F., Caetano, V. R., Haber, S., Savard, M. E., Voldeng, H., et al. (2008). A systemic approach for the development of FHB resistant germplasm accelerates genetic progress. Cereal Res. Commun. 36, 5–9. doi: 10.1556/crc.36.2008.suppl.b.2
Cuthbert, P. A., Somers, D. J., and Brule-Babel, A. (2007). Mapping of Fhb2 on chromosome 6BS: a gene controlling Fusarium head blight field resistance in bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 114, 429–437. doi: 10.1007/s00122-006-0439-3
Cuthbert, R. D., DePauw, R. M., Knox, R. E., Singh, A. K., McCaig, T. N., McCallum, B., et al. (2017). AAC Penhold Canada prairie spring red wheat. Can. J. Plant Sci. 98, 207–214.
Davière, J. M., and Achard, P. (2016). A pivotal role of DELLAs in regulating multiple hormone signals. Mol. Plant 9, 10–20. doi: 10.1016/j.molp.2015.09.011
DePauw, R. M., Knox, R. E., Clarke, F. R., Clarke, J. M., and McCaig, T. N. (2009). Stettler hard red spring wheat. Can. J. Plant Sci. 89, 945–951. doi: 10.4141/cjps08227
DePauw, R. M., Knox, R. E., McCaig, T. N., Clarke, F. R., and Clarke, J. M. (2011a). Muchmore hard red spring wheat. Can. J. Plant Sci. 91, 1–7. doi: 10.1094/cm-2009-0415-01-rs
DePauw, R. M., Knox, R. E., McCaig, T. N., Clarke, F. R., and Clarke, J. M. (2011b). Carberry hard red spring wheat. Can. J. Plant Sci. 91, 529–534. doi: 10.4141/cjps10187
Dixon, L. E., Greenwood, J. R., Bencivenga, S., Zhang, P., Cockram, J., Mellers, G., et al. (2018). TEOSINTE BRANCHED1 regulates inflorescence architecture and development in bread wheat (Triticum aestivum). Plant Cell 30, 563–581. doi: 10.1105/tpc.17.00961
Dvorak, J., Wang, L., Zhu, T., Jorgensen, C. M., Luo, M. C., Deal, K. R., et al. (2018). Reassessment of the evolution of wheat chromosomes 4A, 5A, and 7B. Theor. Appl. Genet. 131, 2451–2462. doi: 10.1007/s00122-018-3165-8
Eudes, F., and Amundsen, E. (2005). Isolated microspore culture of Canadian 6X triticale cultivars. Plant Cell Tissue Organ Cult. 82, 233–241. doi: 10.1007/s11240-005-0867-9
Frohberg, R. C., Stack, R. W., Olson, T., Miller, J. D., and Mergoum, M. (2006). Registration of ‘Alsen’ wheat. Crop Sci. 46, 2311–2312. doi: 10.2135/cropsci2005.12.0501
Gervais, L., Dedryver, F., Morlais, J. Y., Bodusseau, V., Negre, S., Bilous, M., et al. (2003). Mapping of quantitative trait loci for field resistance to Fusarium head blight in an European winter wheat. Theor. Appl. Genet. 106, 961–970. doi: 10.1007/s00122-002-1160-5
Ghavami, F., Elias, E., Mamidi, S., Ansari, O., Sargolzaei, M., Adhikari, T., et al. (2011). Mixed model association mapping for Fusarium head blight resistance in Tunisian-derived durum wheat populations. G3 1, 209–218. doi: 10.1534/g3.111.000489
Giancaspro, A., Giove, S. L., Zito, D., Blanco, A., and Gadaleta, A. (2016). Mapping QTLs for Fusarium head blight resistance in an interspecific wheat population. Front. Plant Sci. 7:1381. doi: 10.3389/fpls.2016.01381
Gilbert, J., and Haber, S. (2013). Overview of some recent research developments in Fusarium head blight of wheat. Can. J. Plant Pathol. 35, 149–174. doi: 10.1080/07060661.2013.772921
Gilbert, J., and Tekauz, A. (2000). Review: recent developments in research on fusarium head blight of wheat in Canada. Can. J. Plant Pathol. 22, 1–8. doi: 10.1080/07060660009501155
Haile, J. K., N’Diaye, A., Walkowiak, S., Nilsen, K. T., Clarke, J. M., Kutcher, H. R., et al. (2019). Fusarium head blight in durum wheat: recent status, breeding directions, and future research prospects. Phytopathology 109, 1664–1675. doi: 10.1094/phyto-03-19-0095-rvw
Hao, Y., Rasheed, A., Zhu, Z., Wulff, B. B. H., and He, Z. (2019). Harnessing wheat Fhb1 for Fusarium resistance. Trends Plant Sci. 25, 1–3. doi: 10.1016/j.tplants.2019.10.006
Hao, Z., Lv, D., Ge, Y., Shi, J., Weijers, D., Yu, G., et al. (2020). RIdeogram: drawing SVG graphics to visualize and map genome-wide data on the idiograms. PeerJ Comput. Sci. 6:e251. doi: 10.7717/peerj-cs.251
He, X., Dreisigacker, S., Singh, R. P., and Singh, P. K. (2019). Genetics for low correlation between Fusarium head blight disease and deoxynivalenol (DON) content in a bread wheat mapping population. Theor. Appl. Genet. 132, 2401–2411. doi: 10.1007/s00122-019-03362-9
He, X., Singh, P. K., Dreisigacker, S., Singh, S., Lillemo, M., and Duveiller, E. (2016). Dwarfing genes Rht-B1b and Rht-D1b are associated with both type I FHB susceptibility and low anther extrusion in two bread wheat populations. PLoS One 11:e0162499. doi: 10.1371/journal.pone.0162499
Herter, C. P., Ebmeyer, E., Kollers, S., Korzun, V., Leiser, W. L., Würschum, T., et al. (2018). Rht24 reduces height in the winter wheat population ‘Solitär×Bussard’ without adverse effects on Fusarium head blight infection. Theor. Appl. Genet. 131, 1263–1272. doi: 10.1007/s00122-018-3076-8
International Wheat Genome Sequencing Consortium (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science 361:eaar7191. doi: 10.1126/science.aar7191
Lemmens, M., Scholz, U., Berthiller, F., Dall’Asta, C., Koutnik, A., Schuhmacher, R., et al. (2005). The ability to detoxify the mycotoxin deoxynivalenol colocalizes with a major quantitative trait locus for fusarium head blight resistance in wheat. Mol. Plant Microbe Interact. 18, 1318–1324. doi: 10.1094/mpmi-18-1318
Lewis, S., Faricelli, M. E., Appendino, M. L., Valarik, M., and Dubcovsky, J. (2008). The chromosome region including the earliness per se locus Eps−Am 1 affects the duration of early developmental phases and spikelet number in diploid wheat. J. Exp. Bot. 59, 3595–3607. doi: 10.1093/jxb/ern209
Li, G., Zhou, J., Jia, H., Gao, Z., Fan, M., Luo, Y., et al. (2019). Mutation of a histidine-rich calcium-binding-protein gene in wheat confers resistance to Fusarium head blight. Nat. Genet. 51, 1106–1112. doi: 10.1038/s41588-019-0426-7
Liu, S., Hall, M. D., Griffey, C. A., and McKendry, A. L. (2009). Meta-analysis of QTL associated with Fusarium head blight resistance in wheat. Crop Sci. 49, 1955–1968. doi: 10.2135/cropsci2009.03.0115
Liu, S., Zhang, X., Pumphrey, M. O., Stack, R. W., Gill, B. S., and Anderson, J. A. (2006). Complex microcolinearity among wheat, rice, and barley revealed by fine mapping of the genomic region harboring a major QTL for resistance to Fusarium head blight in wheat. Funct. Integr. Genomics 6, 83–89. doi: 10.1007/s10142-005-0007-y
Liu, Y., Salsman, E., Fiedler, J. D., Hegstad, J. B., Green, A., Mergoum, M., et al. (2019). Genetic mapping and prediction analysis of FHB resistance in a hard red spring wheat breeding population. Front. Plant Sci. 10:1007. doi: 10.3389/fpls.2019.01007
Lorieux, M. (2012). MapDisto: fast and efficient computation of genetic linkage maps. Mol. Breed. 30, 1231–1235. doi: 10.1007/s11032-012-9706-y
Löffler, M., Schön, C. C., and Miedaner, T. (2009). Revealing the genetic architecture of FHB resistance in hexaploid wheat (Triticum aestivum L.) by QTL meta-analysis. Mol. Breed. 23, 473–488. doi: 10.1007/s11032-008-9250-y
Lu, Q. X., Lillemo, M., Skinnes, H., He, X. Y., Shi, J. R., Ji, F., et al. (2013). Anther extrusion and plant height are associated with Type I resistance to Fusarium head blight in bread wheat line ‘Shanghai-3/Catbird’. Theor. Appl. Genet. 126, 317–334. doi: 10.1007/s00122-012-1981-9
Lv, C., Song, Y., Gao, L., Yao, Q., Zhou, R., Xu, R., et al. (2014). Integration of QTL detection and marker assisted selection for improving resistance to Fusarium head blight and important agronomic traits in wheat. Crop J. 2, 70–78. doi: 10.1016/j.cj.2013.10.004
McCartney, C. A., Brûlé-Babel, A. L., Fedak, G., Martin, R. A., Mccallum, B. D., Gilbert, J., et al. (2016). Fusarium head blight resistance QTL in the spring wheat cross Kenyon/86ISMN 2137. Front Microbiol. 7:1542. doi: 10.3389/fmicb.2016.01542
McCartney, C. A., Somers, D. J., Fedak, G., DePauw, R. M., Thomas, J., Fox, S. L., et al. (2007). The evaluation of FHB resistance QTLs introgressed into elite Canadian spring wheat germplasm. Mol. Breeding 20, 209–221. doi: 10.1007/s11032-007-9084-z
McMullen, M., Bergstrom, G., DeWolf, E., Dill-Macky, R., Hershman, D., Shaner, G., et al. (2012). A unified effort to fight an enemy of wheat and barley: Fusarium head blight. Plant Dis. 96, 1712–1728. doi: 10.1094/pdis-03-12-0291-fe
Mesterhazy, A. (1995). Types and components of resistance to Fusarium head blight of wheat. Plant Breed. 114, 377–386. doi: 10.1111/j.1439-0523.1995.tb00816.x
Mesterhazy, A., Bartok, T., Mirocha, C. G., and Komoroczy, R. (1999). Nature of wheat resistance to Fusarium head blight and the role of deoxynivalenol for breeding. Plant Breed. 118, 97–110. doi: 10.1046/j.1439-0523.2001.00580.x
Mickelson-Young, L., Endo, T. R., and Gill, B. S. (1995). A cytogenetic ladder-map of the wheat homoeologous group-4 chromosomes. Theor. Appl. Genet. 90, 1007–1011. doi: 10.1007/bf00222914
Miedaner, T., and Korzun, V. (2012). Marker-assisted selection for disease resistance in wheat and barley breeding. Phytopathology 102, 560–566. doi: 10.1094/phyto-05-11-0157
Miftahudin, Ross, K., Ma, X.-F., Mahmoud, A. A., Layton, J., Rodriguez Milla, M. A., et al. (2004). Analysis of expressed sequence tag loci on wheat chromosome group 4. Genetics 168, 651–663.
Miller, J. D., Young, J. C., and Sampson, D. R. (1985). Deoxynivalenol and Fusarium head blight resistance in spring cereals. J. Phytopathol. 113, 359–367. doi: 10.1111/j.1439-0434.1985.tb04837.x
Mitchell-Olds, T. (2013). Selection on QTL and complex traits in complex environments. Mol. Ecol. 22, 3427–3429. doi: 10.1111/mec.12345
Ochagavia, H., Prieto, P., Savin, R., Griffiths, S., and Slafer, G. (2018). Dynamics of leaf and spikelet primordia initiation in wheat as affected by Ppd-1a alleles under field conditions. J. Exp. Bot. 69, 2621–2631. doi: 10.1093/jxb/ery104
Paillard, S., Schnurbusch, T., Tiwari, R., Messmer, M., Winzeler, M., Keller, B., et al. (2004). QTL analysis of resistance to Fusarium head blight in Swiss winter wheat (Triticum aestivum L.) Theor. Appl. Genet. 109, 323–332. doi: 10.1007/s00122-004-1628-6
Paul, P. A., Lipps, P. E., and Madden, L. V. (2005). Relationship between visual estimates of Fusarium head blight intensity and deoxynivalenol accumulation in harvested wheat grain: a meta-analysis. Phytopathology 95, 1225–1236. doi: 10.1094/phyto-95-1225
Prat, N., Guilbert, C., Prah, U., Wachter, E., Steiner, B., Langin, T., et al. (2017). QTL mapping of Fusarium head blight resistance in three related durum wheat populations. Theor. Appl. Genet. 130, 13–27. doi: 10.1007/s00122-016-2785-0
R Core Team (2014). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Rasheed, A., Wen, W., Gao, F., Zhai, S., Jin, H., Liu, J., et al. (2016). Development and validation of KASP assays for genes underpinning key economic traits in bread wheat. Theor. Appl. Genet. 129, 1843–1860. doi: 10.1007/s00122-016-2743-x
Rawat, N., Pumphrey, M. O., Liu, S., Zhang, X., Tiwari, V. K., Ando, K., et al. (2016). Wheat Fhb1 encodes a chimeric lectin with agglutinin domains and a pore-forming toxin-like domain conferring resistance to Fusarium head blight. Nat. Genet. 48, 1576–1580. doi: 10.1038/ng.3706
Ren, J., Wang, Z., Du, Z., Che, M., Zhang, Y., Quan, W., et al. (2019). Detection and validation of a novel major QTL for resistance to Fusarium head blight from Triticum aestivum in the terminal region of chromosome 7DL. Theor. Appl. Genet. 132, 241–255. doi: 10.1007/s00122-018-3213-4
Sari, E., Berraies, S., Knox, R. E., Singh, A. K., Ruan, Y., Cuthbert, R. D., et al. (2018). High density genetic mapping of Fusarium head blight resistance QTL in tetraploid wheat. PLoS One 13:e0204362. doi: 10.1371/journal.pone.0204362
Semagn, K., Skinnes, H., Bjornstad, A., Maroy, A. G., and Tarkegne, Y. (2007). Quantitative trait loci controlling Fusarium head blight resistance and low deoxynivalenol content in hexaploid wheat population from ‘Arina’ and NK93604. Crop Sci. 47, 294–303. doi: 10.2135/cropsci2006.02.0095
Shaw, L. M., Turner, A. S., Herry, L., Griffiths, S., and Laurie, D. A. (2013). Mutant alleles of Photoperiod-1 in wheat (Triticum aestivum L.) that confer a late flowering phenotype in long days. PLoS One 8:e79459. doi: 10.1371/journal.pone.0079459
Skinnes, H., Semagn, K., Tarkegne, Y., Maroy, A. G., and Bjornstad, A. (2010). The inheritance of anther extrusion in hexaploid wheat and its relationship to Fusarium head blight resistance and deoxynivalenol content. Plant Breed. 129, 149–155. doi: 10.1111/j.1439-0523.2009.01731.x
Somers, D. J., Fedak, G., and Savard, M. (2003). Molecular mapping of novel genes controlling Fusarium head blight resistance and deoxynivalenol accumulation in spring wheat. Genome 46, 555–564. doi: 10.1139/g03-033
Srinivasachary, G. N., Steed, A., Faure, S., Bayles, R., Jennings, P., and Nicholson, P. (2008). Mapping of QTL associated with Fusarium head blight in spring wheat RL4137. Czech. J. Genet. Plant 44, 147–159. doi: 10.17221/70/2008-cjgpb
Srinivasachary, G. N., Steed, A., Hollins, T., Bayles, R., Jennings, P., and Nicholson, P. (2009). Semi-dwarfing Rht-B1 and Rht-D1 loci of wheat differ significantly in their influence on resistance to Fusarium head blight. Theor. Appl. Genet. 118, 695–702. doi: 10.1007/s00122-008-0930-0
Steiner, B., Buerstmayr, M., Michel, S., Schweiger, W., Lemmens, M., and Buerstmayr, H. (2017). Breeding strategies and advances in line selection for Fusarium head blight resistance in wheat. Trop. Plant Pathol. 42, 165–174. doi: 10.1007/s40858-017-0127-7
Steiner, B., Buerstmayr, M., Wagner, C., Danler, A., Eshonkulov, B., Ehn, M., et al. (2019). Fine-mapping of the Fusarium head blight resistance QTL Qfhs.ifa-5A identifies two resistance QTL associated with anther extrusion. Theor. Appl. Genet. 132, 2039–2053. doi: 10.1007/s00122-019-03336-x
Steiner, B., Lemmens, M., Griesser, M., Scholz, U., Schondelmaier, J., and Buerstmayr, H. (2004). Molecular mapping of resistance to Fusarium head blight in the spring wheat cultivar Frontana. Theor. Appl. Genet. 109, 215–224. doi: 10.1007/s00122-004-1620-1
Su, Z., Bernardo, A., Tian, B., Chen, H., Wang, S., Ma, H., et al. (2019). A deletion mutation in TaHRC confers Fhb1 resistance to Fusarium head blight in wheat. Nat. Genet. 51, 1099–1105. doi: 10.1038/s41588-019-0425-8
Szabo-Hever, A., Lehoczki-Krsjak, S., Varga, M., Purnhauser, L., Pauk, J., Lantos, C., et al. (2014). Differential influence of QTL linked to Fusarium head blight, Fusarium-damaged kernel, deoxynivalenol contents and associated morphological traits in a Frontana-derived wheat population. Euphytica 200, 9–26. doi: 10.1007/s10681-014-1124-2
Venske, E., Dos Santos, R. S., Farias, D. D. R., Rother, V., Maia, L. C. D., Pegoraro, C., et al. (2019). Meta-analysis of the QTLome of Fusarium head blight resistance in bread wheat: refining the current puzzle. Front. Plant Sci. 10:727. doi: 10.3389/fpls.2019.00727
Wang, S., Wong, D., Forrest, K., Allen, A., Chao, S., Huang, B. E., et al. (2014). Characterization of polyploid wheat Genome diversity using a high-density 90,000 single nucleotide polymorphism array. Plant Biotechnol. J. 12, 787–796. doi: 10.1111/pbi.12183
Wang, W., Basten, C. J., and Zeng, Z. B. (2007). Windows QTL Cartographer 2.5. Raleigh, NC: Department of Statistics, NCSU.
Wolde, G. M., Trautewig, C., Mascher, M., and Schnurbusch, T. (2019). Genetic insights into morphometric inflorescence traits of wheat. Theor. Appl. Genet. 132, 1661–1676. doi: 10.1007/s00122-019-03305-4
Wu, L., Zhang, Y., He, Y., Jiang, P., Zhang, X., and Ma, H. X. (2019). Genome-wide association mapping of resistance to fusarium head blight spread and deoxynivalenol accumulation in Chinese elite wheat germplasm. Phytopathology 109, 1208–1216. doi: 10.1094/phyto-12-18-0484-r
Xu, K., He, X., Dreisigacker, S., He, Z., and Singh, P. (2020). Anther extrusion and its association with Fusarium head blight in CIMMYT wheat germplasm. Agronomy 10:47. doi: 10.3390/agronomy10010047
Xue, A. G., Butler, G., Voldeng, H. D., Fedak, G., and Savard, M. E. (2006). Comparison of the influence of inoculum sources on the development of fusarium head blight and the deoxynivalenol content in spring wheat in a disease nursery. Can. J. Plant Pathol. 28, 152–159. doi: 10.1080/07060660609507282
Xue, S., Xu, F., Tang, M., Zhou, Y., Li, G., An, X., et al. (2011). Precise mapping Fhb5, a major QTL conditioning resistance to Fusarium infection in bread wheat (Triticum aestivum L.). Theor. Appl. Genet. 123, 1055–1063. doi: 10.1007/s00122-011-1647-z
Yan, W., Li, H. B., Cai, S. B., Ma, H. X., Rebetzke, G. J., and Liu, C. J. (2011). Effects of plant height on type I and type II resistance to fusarium head blight in wheat. Plant Pathol. 60, 506–512. doi: 10.1111/j.1365-3059.2011.02426.x
Ye, Z. (2015). Analysis of Resistance to Fusarium Head Blight (FHB) in Winter Wheat and Evaluation of Genetics and Cultural Practices for FHB Mitigation. Ph.D. thesis, University of Manitoba, Manitoba.
Zhang, Q., Axtman, J. E., Faris, J. D., Chao, S., Zhang, Z., Friesen, T. L., et al. (2014). Identification and molecular mapping of quantitative trait loci for Fusarium head blight resistance in emmer and durum wheat using a single nucleotide polymorphism-based linkage map. Mol. Breed. 34, 1677–1687. doi: 10.1007/s11032-014-0180-6
Zhang, W. T., Francis, T., Gao, P., Boyle, K., Jiang, F. Y., Eudes, F., et al. (2018). Genetic characterization of type II Fusarium head blight resistance derived from transgressive segregation in a cross between eastern and western Canadian spring wheat. Mol. Breed. 38:13.
Keywords: Fusarium head blight, resistance components, QTL, genetic architecture, DON, SNP, flowering
Citation: Zhang W, Boyle K, Brûlé-Babel AL, Fedak G, Gao P, Robleh Djama Z, Polley B, Cuthbert RD, Randhawa HS, Jiang F, Eudes F and Fobert PR (2020) Genetic Characterization of Multiple Components Contributing to Fusarium Head Blight Resistance of FL62R1, a Canadian Bread Wheat Developed Using Systemic Breeding. Front. Plant Sci. 11:580833. doi: 10.3389/fpls.2020.580833
Received: 07 July 2020; Accepted: 16 September 2020;
Published: 26 October 2020.
Edited by:
Jose Luis Gonzalez Hernandez, South Dakota State University, United StatesReviewed by:
Shaobin Zhong, North Dakota State University, United StatesZhi Zheng, Plant Industry (CSIRO), Australia
Copyright © 2020 Zhang, Boyle, Brûlé-Babel, Fedak, Gao, Robleh Djama, Polley, Cuthbert, Randhawa, Jiang, Eudes and Fobert. 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: Pierre R. Fobert, UGllcnJlLkZvYmVydEBucmMtY25yYy5nYy5jYQ==; Wentao Zhang, V2VudGFvLlpoYW5nQG5yYy1jbnJjLmdjLmNh