Front. Plant Sci.Frontiers in Plant ScienceFront. Plant Sci.1664-462XFrontiers Media S.A.10.3389/fpls.2019.00727Plant ScienceSystematic ReviewMeta-Analysis of the QTLome of Fusarium Head Blight Resistance in Bread Wheat: Refining the Current PuzzleVenskeEduardo1dos SantosRailson Schreinert2FariasDaniel da Rosa3RotherVianei1da MaiaLuciano Carlos1PegoraroCamila1Costa de OliveiraAntonio1*1Crop Science Department, Plant Genomics and Breeding Center, Eliseu Maciel School of Agronomy, Federal University of Pelotas, Pelotas, Brazil2Instituto Federal de Educação, Ciência e Tecnologia Farroupilha (IFFar), Alegrete, Brazil3Instituto Federal de Educação, Ciência e Tecnologia Catarinense (IFC), Araquari, Brazil
Edited by: Luigi Cattivelli, Council for Agricultural and Economics Research, Italy
Reviewed by: Agata Gadaleta, University of Bari Aldo Moro, Italy; Elias Elias, North Dakota State University, United States; Thomas Miedaner, University of Hohenheim, Germany; Aakash Chawade, Swedish University of Agricultural Sciences, Sweden
*Correspondence: Antonio Costa de Oliveira acostol@terra.com.br
This article was submitted to Plant Breeding, a section of the journal Frontiers in Plant Science
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.
Background: Fusarium Head Blight (FHB) is a worldwide devastating disease of bread wheat (Triticum aestivum L.). Genetic resistance is the most effective way to control FHB and many QTL related to this trait have been mapped on the wheat genetic map. This information, however, must be refined to be more efficiently used in breeding programs and for the advance of the basic research. The objective of the present study was to in-depth analyze the QTLome of FHB resistance in bread wheat, further integrating genetic, genomic, and transcriptomic data, aiming to find candidate genes.
Methods: An exhaustive bibliographic review on 76 scientific papers was carried out collecting information about QTL related to FHB resistance mapped on bread wheat. A dense genetic consensus map with 572,862 loci was generated for QTL projection. Meta-analysis could be performed on 323 QTL. Candidate gene mining was carried out within the most refined loci, containing genes that were cross-validated with publicly available transcriptional expression data of wheat under Fusarium infection. Most highlighted genes were investigated for protein evidence.
Results: A total of 556 QTL were found in the literature, distributed on all sub-genomes and chromosomes of wheat. Meta-analysis generated 65 meta-QTL, and this refinement allows one to find markers more tightly linked to these regions. Candidate gene mining within the most refined meta-QTL, meta-QTL 1/chr. 3B, harvested 324 genes and transcriptional data cross-validated 10 of these genes, as responsive to FHB. One is of these genes encodes a Glycosiltransferase and the other encodes for a Cytochrome P450, and these such proteins have already been verified as being responsible for FHB resistance, but the remaining eight genes still have to be further studied, as promising loci for breeding.
Conclusions: The QTLome of FHB resistance in wheat was successfully assembled and a refinement in terms of number and length of loci was obtained. The integration of the QTLome with genomic and transcriptomic data has allowed for the discovery of promising candidate genes for use in breeding programs.
molecular markersmeta-QTLgenetic architecturegenetic mapsgenometranscriptomebioinformaticsConselho Nacional de Desenvolvimento Científico e Tecnológico10.13039/501100003593Coordenação de Aperfeiçoamento de Pessoal de Nível Superior10.13039/501100002322Fundação de Amparo á Pesquisa do Estado do Rio Grande do Sul10.13039/501100004263Introduction
Fusarium Head Blight (FHB), also known as scab, is one of the most important diseases of bread or hexaploid wheat (Triticum aestivum L.), causing expressive losses every year in many countries (Buerstmayr et al., 2014; Steiner et al., 2017). It is a floral infection disease mainly caused by the fungus Fusarium graminearum, but also by F. culmorum and F. avenaceum, being highly influenced by the environment and it is difficult to control by chemical means (Parry et al., 1995; Bai and Shaner, 2004; Steiner et al., 2017). It is a serious problem for wheat production because it causes not only quantitative losses, reducing yield, but also qualitative damages, while the pathogen produces a set of mycotoxins, which are highly toxic to human and animals (Del Ponte et al., 2004; Draeger et al., 2007; Pestka, 2010). Among the most important mycotoxins, one can list Deoxynivalenol (DON) and Nivalenol (NIL), which are also pathogenicity factors (Draeger et al., 2007).
The most effective and convenient way to control this disease is through the adoption of resistant cultivars (Steiner et al., 2017). FHB resistance is a quantitative trait and it has been classically proposed and well-accepted that wheat can present five different types of resistance to this disease: viz., (I) to the initial infection of the pathogen within spike tissues; (II) to the subsequent spread of the disease throughout the spike; (III) to the accumulation of mycotoxins; (IV) to kernel infection and (V) to yield performance (Schroeder and Christensen, 1963; Miller and Arnison, 1986; Mesterhazy, 1995; Mesterhazy et al., 1999; Buerstmayr and Lemmens, 2015). Also, in some studies, resistance types I and II are assayed jointly, as at field conditions, i.e., without control regarding the inoculum presence and abundance, it is not possible to precisely determine the contribution of initial infection and fungal spread on each diseased head. Regarding methods to assay each of all these types of resistance, literature describes different protocols with this aim, from inoculation aspects to evaluation procedures, of which the descriptions, however, are not the main objective of the present study.
It has been recognized that wheat presents expressive genetic variability for FHB resistance (Steiner et al., 2017). Hundreds of QTL contributing to this trait were already mapped on the genetic map of the cereal, distributed evenly throughout all genomes and chromosomes of the species and verified in dozens of different cultivars and breeding lines (Buerstmayr et al., 2009; Kugler et al., 2013; Steiner et al., 2017). Despite of the expressive amount of QTL related to this trait mapped so far, there is a consensus among the researchers that only very few of these loci have in fact been utilized for the development of improved cultivars, resistant to this disease (Buerstmayr et al., 2014; Steiner et al., 2017).
In order to make these already mapped QTL more useful to plant breeding, as well as, for the basic research aiming to better understand FHB resistance in wheat and other small grain cereals, a comprehensive and in depth analysis of all these loci has to be carried out. However, due to the fact that most mapping experiments differ in numerous aspects, such as type of mapping population, number of lines, number, and type of molecular markers used and design and even the general quality of the experiment, an elaborated method has to be applied for this purpose. In this regard, QTL-meta analysis has been proven as an efficient approach, and has been constantly improved (Goffinet and Gerber, 2000; Salvi and Tuberosa, 2015).
In general, this analysis aims to verify whether QTL observed in isolated studies correspond to different loci or whether they actually represent a common position on the genetic map of the studied species, as well as the opposite i.e., if a QTL hitherto considered as mapped in a common region in diverse genotypes actually represents different positions. Finally, the analysis aims to stablish the occurrence of QTL “hotspots” in a consensus map, which correspond to the more precise region where these loci represent under analysis, the so-called “meta-QTL.” These refined genomic regions harbor important and, in many cases, several genes related to the trait (Goffinet and Gerber, 2000; Salvi and Tuberosa, 2015). This approach was already applied to many different crops and traits, such as cist nematode resistance in soybean (Guo et al., 2006), yield under drought conditions in rice (Swamy et al., 2011), diverse yield and quality traits in cotton (Said et al., 2015), yield in maize (Martinez et al., 2016) and yield, baking quality, and grain protein content related traits in wheat (Quraishi et al., 2017).
Three QTL meta-analysis studies were already performed for FHB resistance in wheat (Liu et al., 2009; Löffler et al., 2009; Mao et al., 2010—in chronological sequence). Among these studies, virtually simultaneous, one of special relevance was that published by Liu et al. (2009), who brought the largest amount of information generated until the year of 2009. Despite the comprehensiveness of the bibliographic review performed and the contribution in projecting distinct QTL in a common consensus map, the research could not take advantage of the current improved methods to determine meta-QTL. For the studies of Löffler et al. (2009) and Mao et al. (2010), they considered fewer QTL in the analyses, and the later focused specially in understanding the association between FHB resistance and plant height.
Many advances have been taking place that have currently improved algorithms but are still in constant evolution, and are available for QTL-meta analysis. Among these are the ones implemented in the software BioMercator (Arcade et al., 2004; Veyrieras et al., 2007; Sosnowski et al., 2012). These algorithms test, through the maximum likelihood method, the number of possible meta-QTL generated from the projection of distinct original QTL in a linkage group, and apply statistical criteria to suggest the best model to be used for deciding the actual number of these meta-loci (Goffinet and Gerber, 2000; Sosnowski et al., 2012).
The QTL meta-analysis for FHB resistance in wheat, carried out by Liu et al. (2009), Löffler et al. (2009), and Mao et al. (2010), utilized basically RFLP (Restriction Fragment Length Polymorphism) and SSR (Simple Sequence Repeat) markers. During and in the subsequent years after these studies, a revolution has taken place in terms of genotyping technology. Distinct technologies for SNP (Single Nucleotide Polymorphism) assessment have been developed (Akbari et al., 2006; Poland J. A. et al., 2012; Poland J. et al., 2012; Semagn et al., 2014; Wang et al., 2014; Winfield et al., 2016). From few RFLP or SSR loci by study, an expressive advance has taken place for the utilization of platforms which can evaluate thousands of loci at time, such as the Axiom® 820K SNP array, which can assess near to 820 thousand SNPs at once (Winfield et al., 2016). Furthermore, various studies in QTL mapping were performed since 2009 until now, thus there is an important accumulation of information to be treated or updated.
The other step-change, which was made recently in wheat research, was the completion of the genome sequence (International Wheat Genome Sequencing Consortium, 2014) and more recently, the genome annotation (International Wheat Genome Sequencing Consortium, 2018), which has been made freely available for the scientific community. In the same way, a large amount of transcriptomic data, including related to FHB response, has been made available, through a user-friendly platform (Borrill et al., 2016; Ramírez-González et al., 2018). All these resources represent a rich opportunity to unveil, in an in-depth way, the genetic architecture of FHB resistance in wheat.
Considering the exposed, the aim of this work was to perform a comprehensive meta-analysis of QTL related to FHB resistance in hexaploid wheat, further integrating genetic, genomic, and transcriptomic data, aiming to find promising candidate genes. It is expected that this work will contribute to both the improvement in strategies for wheat breeding for FHB resistance but also to the basic research aiming to better elucidate the genetic architecture underlying the trait.
Materials and Methods
The approach applied in this study comprises several steps, thus, to facilitate the understanding, a flow diagram is presented in Figure 1. Further description of each step is given as follows.
Flow diagram of the systematic review and QTL meta-analysis applied in this study, which further incorporated genomic and transcriptomic publicly available data.
Search Strategy—Collection and Characterization of QTL for FHB Resistance in Wheat
An exhaustive bibliographic review was carried out on 76 papers published between 1999 and 2017, searching for information on QTL mapped in bread wheat (T. aestivum L.) related to FHB resistance (Table 1, Supplementary File 1). These studies include the ones utilized in the meta-analysis performed by Liu et al. (2009), Löffler et al. (2009), and Mao et al. (2010).
Publications from which the data was collected to form the database for this study.
Year of publication
Number of papers
References
1999
1
Waldron et al., 1999
2000
1
Gupta et al., 2000
2001
2
Anderson et al., 2001; Gupta et al., 2001
2002
2
Buerstmayr et al., 2002; Zhou et al., 2002
2003
6
Bourdoncle and Ohm, 2003; Buerstmayr et al., 2003; Gervais et al., 2003; Shen et al., 2003a,b; Somers et al., 2003
2004
6
Bai and Shaner, 2004; Lin et al., 2004; Paillard et al., 2004; Steiner et al., 2004; Zhang et al., 2004; Zhou et al., 2004
2005
5
Gilsinger et al., 2005; Jia et al., 2005; Mardi et al., 2005; Schmolke et al., 2005; Yang J. et al., 2005; Yang Z. et al., 2005
2006
7
Chen et al., 2006; Cuthbert et al., 2006; Lin et al., 2006; Liu et al., 2006; Ma et al., 2006a,b; Mardi et al., 2006
2007
10
Cuthbert et al., 2007; Draeger et al., 2007; Häberle et al., 2007; Jiang et al., 2007a,b; Klahr et al., 2007; Li et al., 2007; Liu et al., 2007; Semagn et al., 2007; Zhang and Mergoum, 2007
2008
7
Abate et al., 2008; Handa et al., 2008; Holzapfel et al., 2008; Malla et al., 2008a,b; Srinivasachary et al., 2008; Yu et al., 2008
2009
4
Bonin and Kolb, 2009; Häberle et al., 2009; Liu et al., 2009; Srinivasachary et al., 2009
2010
3
Malla et al., 2010; Zhang et al., 2010; Zhou et al., 2010
2011
3
Chu et al., 2011; Jayatilake et al., 2011; Li et al., 2011
2012
8
Basnet et al., 2012; Li et al., 2012; Liu et al., 2012; Lu et al., 2012; Suzuki et al., 2012; Szabó-Hevér et al., 2012; Zhang et al., 2012a,b
2013
1
Cativelli et al., 2013
2014
3
Cai and Bai, 2014; Chao et al., 2014; Szabó-Hevér et al., 2014
2015
2
Buerstmayr and Buerstmayr, 2015; Zhu et al., 2015
2016
4
He et al., 2016; Islam et al., 2016; McCartney et al., 2016; Petersen et al., 2016
2017
1
Petersen et al., 2017
For each QTL, the following information was collected: (1) type of resistance (Types I, II, I and II combined, III, IV, V), which are considered different traits; (2) genotype source of resistance, (3) type of the mapping population [(F2, backcross, recombinant inbreed lines (RILs) and double-haploids (DH)]; (4) number of assayed lines; (5) LOD (logarithm of the odds) score; (6) percentage of the phenotypic variance explained by the QTL (R2); (7) the flanking or single marker(s) (for interval mapping and single marker analysis, respectively). In the few cases in which the source-paper presented the p-value statistic or simply stated that a minimum LOD score of 3 was adopted to map the QTL, instead of informing the actual value, the LOD 3 was assumed in this study. Studies of association mapping or those using other non-common mapping populations were not included in this work.
Consensus Map Development for QTL Projection
To project the diverse QTL collected, which were mapped using different types of markers, a consensus genetic map of wheat, resulting from the union of distinct individual maps, was developed, which was the basis for this study. The maps selected for the development of the consensus map, in a total of seven, were the ones available which harbor a large amount of the markers used for mapping the QTL for FHB resistance in wheat, whose information was obtained previously.
Thus, from GrainGenes website (https://wheat.pw.usda.gov/GG3/), the following maps were obtained: “Wheat, Consensus SSR, 2004,” “Wheat, Composite, 2004” e “Wheat, Synthetic x Opata, BARC,” as a source of SSR markers and few classical genes, such as Vrn1—a gene for vernalization sensitivity, which have also been used for QTL mapping in some studies. The “Wheat consensus map version 4.0,” source of DArT markers was obtained from the website of the company owner of this technology (https://www.diversityarrays.com). For chip-based SNP markers, assessed by the platforms Illumina®9K iSelect Beadchip Assay and Illumina®iSelect 90K SNP Assay, the maps were obtained from Cavanagh et al. (2013) and Wang et al. (2014), respectively. Finally, for SNPs obtained and mapped through the GBS (genotyping by sequencing) approach, the map of Saintenac et al. (2013) was used. With the exception of “Wheat, Synthetic x Opata, BARC,” all remaining maps are already consensus maps between different populations, thus they represent high quality maps. For the construction of the consensus map, the R package LPMerge was used (Endelman and Plomion, 2014).
Projection and QTL Meta-Analysis
From the 556 QTL collected, 359 could be projected onto the consensus map generated, given that the markers used to map them were present in the consensus map. However, even QTL which could not be projected or analyzed via the meta-analysis were still collected and characterized, being present in a catalog with all information collected in this study (Supplementary File 1).
For the projection of these QTL, initially the confidence interval (CI) of 95% was calculated for each locus, through the equations below, which have been modeled for each mapping population (Darvasi and Soller, 1997; Guo et al., 2006): F2 and Backcross:
CI=(530(NumberoflinesxR2))
RILs:
CI=(163(NumberoflinesxR2))
Double-haploids:
CI=(287(NumberoflinesxR2))
Next, using the software BioMercator V4.2 (Arcade et al., 2004; Sosnowski et al., 2012), the QTL, represented by their middle points, along with their calculated confidence intervals, original LOD score and R2, were directly projected onto the consensus map previously developed, and the meta-analysis was carried out, individually by chromosome, through the Veyrieras two-step algorithm in the software (Veyrieras et al., 2007). The Akaike (AIC) statistics was used to define the best model for the definition of the number of meta-QTL or “real” QTL, which best represent the original QTL. The algorithms and statistical procedures implemented in this software are well-described in the literature (Arcade et al., 2004; Veyrieras et al., 2007; Sosnowski et al., 2012). As a requirement of the method, the meta-analysis could only be performed with linkage-groups with a minimum of 10 projected QTL. Attempts to run analysis when <10 QTL were projected returned in error. Thus, 323 QTL could, in fact, be considered for the QTL meta-analysis. All files prepared to run BioMercator V.4.2, i.e., maps and QTL files for each wheat chromosome, are made available in Supplementary File 2.
Candidate Gene Mining, Transcriptional Expression Analysis, and Functional Annotation Within Meta-QTL
Among the meta-QTL obtained, the ones considered highly reliable and refined were considered for gene mining. The criteria adopted were as follows: (1) the meta-QTL is generated through the projection of at least two overlapping QTL (otherwise the meta-QTL would simply correspond to a single QTL previously known); (2) it is shorter than 1.0 cM in genetic distance; (3) it is shorter than 20 Mb in physical distance (at Chinese Spring wheat reference genome). For that, the physical positions of markers present at the flanking regions of each highly refined meta-QTL (corresponding to the end of the confidence interval) were searched at the International Wheat Genome Sequencing Consortium reference genome sequence repository at https://wheat-urgi.versailles.inra.fr/Seq-Repository/BLAST through BLAST search of their sequences or either directly through surveying the annotation browser (International Wheat Genome Sequencing Consortium, 2014, 2018). The “IWGSC RefSeq v1.0” version was used. High confidence annotated genes (HighConfidenceGenesv1.1), within each highly refined meta-QTL, were then listed and thereafter called candidate genes.
For the transcriptional expression analysis, the collected candidate genes were analyzed through the expVIP (expression Visualization and Integration Platform)/Wheat Expression Browser online resource (http://www.wheat-expression.com) (Borrill et al., 2016; Ramírez-González et al., 2018). The available data for “Fusarium head blight infected spikelets” (Kugler et al., 2013) was used. This study assayed the transcriptome profile of a set of wheat lines after F. graminearum point inoculation and control conditions (mock) at 30 and 50 h after treatment. The tissue analyzed was spikelets. From this data set, only transcripts from the varieties, CM-82036 (resistant to FHB) and a susceptible near isogenic line, were selected (a cross between CM-82036 and the susceptible variety Remus). Differentially expressed genes (considering disease × mock conditions, only) were collected, based on the standard error. Following the proposed criteria by Wagner et al. (2013), only genes presenting at least 2 transcripts per million (TPM) were considered in this search. Other transcriptomic studies approaching other sources of FHB resistance could not be used in this work as they were not made available in manageable platforms like the one we have used here. Next, differentially expressed genes were further investigated for functional annotation (i.e., protein evidence), also at the International Wheat Genome Sequencing Consortium reference genome sequence repository (https://wheat-urgi.versailles.inra.fr/Seq-Repository/BLAST).
ResultsCharacterization of QTL for FHB Resistance in Wheat
A thorough search in the literature resulted in the collection of information regarding 556 QTL described as responsive for FHB resistance, already mapped on a bread wheat genetic map, distributed on all wheat sub-genomes and chromosomes (Figure 2). The wheat B sub-genome had the largest number of QTL mapped, with 238, followed by A, with 192 and finally D, with 121 QTL. The chromosomes 3B (81 QTL), 5A (58 QTL), and 2D (57 QTL) were the ones which presented the largest amount of these mapped regions. The most common type of resistance mapped was the type II (spread of the disease throughout the spike), with 41.5%, and the less mapped trait was the type I (initial infection), with 11.5% of the collected QTL. No QTL related to type V resistance was found, the resistance type corresponding to the ability of the genotype to tolerate the pathogen without losing yield potential. Five QTL were not assigned in the original study (Szabó-Hevér et al., 2014) to any wheat sub-genome or chromosome. The QTL collected were mapped from 77 different genotypes (cultivars or breeding lines), developed in many different countries (Supplementary File 1). Also, some of these loci were verified as belonging to the susceptible or moderately susceptible parent of the cross, which generated the mapping population, e.g., as in Cai and Bai (2014).
QTL related to FHB resistance in wheat published from 1999 to 2017, classified according the type of resistance they confer and the sub-genome and chromosome of occurrence.
Generation of the Consensus Map
Seven wheat genetic maps were used for the development of the consensus map, the base for the QTL projection and meta-analysis performed in this study. The generated map was composed by SSR, DArT, SNPs (from different technologies) and few genes, presenting 572,862 loci in total. On average, each chromosome presented 27,279 markers and was 217.49 cM long (Table 2). Attempts to incorporate more maps, containing thousands of other loci, aiming to form an even larger consensus map, were performed, however numerous ordering conflicts appeared. Thus, the maps integrated in this study presented the best harmonious combination. The consensus map, separated by chromosome (or linkage group), is available through Supplementary File 2.
Consensus genetic map of wheat, generated from different maps and composed by SSR, DArT, SNP molecular markers and few genes.
Linkage group (chromosome)
Consensus map
Length (cM)
Number of markers
1A
255.6
24,278
1B
288.57
39,996
1D
208.05
16,029
2A
251.3
27,131
2B
316.13
47,062
2D
166.17
29,300
3A
207.27
23,619
3B
204.28
44,087
3D
188.33
21,470
4A
272.1
29,584
4B
135.76
18,828
4D
170.43
8,461
5A
274.91
28,558
5B
226.9
34,023
5D
223.97
15,264
6A
183.44
24,656
6B
134.46
31,794
6D
203
16,482
7A
226.6
33,174
7B
188.64
38,038
7D
241.29
21,028
Average
217.49
27,279
Total
4567.2
572,862
QTL Projection and Meta-Analysis
From the 556 QTL collected from literature, a total of 359 could be projected on the consensus genetic map generated (Figures 3–9). It was decided not to include in the meta-analysis QTL mapped with the aid of AFLP and RFLP markers, as they usually originate from less dense maps, leading to large QTL. It could impair the attempt to refine the QTL regions. The remaining QTL which could not be included in the meta-analysis were due to the absence of their flanking markers in the consensus map developed.
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 1A, 1B, and 1D. In the projection, the colors of the QTL means: black, type I resistance; red, type II; blue, type I and II combined; green, type III; orange, type IV. In the meta-analysis, the colors both in the QTL and inside the linkage groups only mean the different meta-QTL generated and how each QTL contributed to the formation of it. Numbering of the meta-QTL is located to the left side of the linkage group and to the right side is indicated the genetic distance which comprise these loci. The molecular markers and their positions are omitted in the figures for better visualization. The line in the middle of each QTL represents the LOD score of it in the original work.
Given that the method applied in this work requires the minimum of 10 QTL projected on a linkage group (chromosome) to launch a meta-analysis, it was possible to run the analysis in 15 out of the 21 linkage groups of hexaploid wheat. Thus, from the 359 QTL projected, 323 definitely composed the database for meta-QTL discovery. As a result of the meta-analysis, these 323 QTL gave origin to 65 meta-QTL, from which 33 were localized on the A sub-genome, 25 on the B sub-genome and 7 on the D sub-genome. The genetic length of these loci ranged from 0.82 to 55.7 cM with an average of 6.08 cM. For 1D, 3D, 5D, 6D, 7B, and 7D, the number of QTL projected (<10) was insufficient to perform the meta-analysis. In addition to the graphical form, which is presented here, all meta-QTL generated are available in Supplementary File 3, projected onto the genetic map, along with all markers comprising their entire length.
In chromosome 1A, 19 QTL were projected, and gave origin to 5 meta-QTL. In 1B, on the other hand, 19 QTL were also projected, however, only 2 meta-QTL were discovered, which was the smallest number of these refined loci obtained per chromosome, among all chromosomes that could undergo this analysis (Figure 3). For 1D, although QTL meta-analysis could not be performed, it is possible to verify that 5 out of the 8 projected loci co-localized in the linkage group. The homeologous group 2 (chr 2A, 2B, and 2D) was the one where the largest number of meta-QTL were verified (Figure 4). In chr 2A, from 22 QTL, six meta-QTL were calculated, and for 2B and 2D, four meta-QTL were found from each of the analyses of 12 and 26 projected QTL, respectively. In chr. 3A, 4 meta-QTL were detected from 12 projected QTL, but in 3D only the projection of 5 QTL could be performed (Figure 5). It was in chr. 3B where the largest amount of QTL was projected in this study, in a total of 62, which gave origin to 5 meta-QTL.
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 2A, 2B, and 2D. The color system and numbering follows the same as Figure 3.
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 3A, 3B, and 3D. The color system follows the same as Figure 3. The meta-QTL 1 in 3B was selected for gene mining.
Similarly, to which was verified for the homeologous group 2, it was possible to perform the meta-analysis on all 3 chr. of the group 4, being 5, 3, and 3 meta-QTL verified in 4A, 4B, and 4D, respectively (Figure 6). Regarding the group 5, chr. 5A was the second in this work in terms of number of QTL projected, in a total of 44, which allowed the detection of 5 meta-QTL (Figure 7). For chr. 5B, seven meta-QTL were obtained, the largest number among all wheat chromosomes, however for 5D, it was only possible to do the projection of 4 QTL. Chromosomes 6A and 6B presented 4 meta-QTL each, while for 6D, the projection of only three QTL was possible, which was the smallest number of QTL projected in a linkage group in this study (Figure 8). The homologous group 7 was the one in which the smallest number of QTL were projected and the smallest number of meta-QTL were obtained, as the meta-analysis could only be carried out for 7A and 4 meta-QTL were found (Figure 9). As a special feature in this linkage group, the projection of one QTL (QTL 12) with an imprecise confidence interval led to the formation of one meta-QTL (meta-QTL 4/7A) even larger than the original QTL, being the largest meta-QTL verified in this study (55.7 cM).
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 4A, 4B, and 4D. The color system follows the same as Figure 3.
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 5A, 5B, and 5D. The color system follows the same as Figure 3.
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 6A, 6B, and 6D (only projection). The color system follows the same as Figure 3.
Projection (left) and, with green borders meta-QTL of FHB resistance in wheat for chr. 7A, and only projection for 7B and 7D. The color system follows the same as Figure 3.
Candidate Gene Mining, Transcriptional Expression Within Meta-QTL, and Functional Annotation
Among the 65 meta-QTL obtained, only one fulfilled the previously established criteria for candidate gene mining, i.e., simultaneously resulting from the analysis of at least two overlapping QTL, being shorter than 1.0 cM in genetic length and also shorter than 20 Mb in physical length (based on the wheat genome reference sequence of Chinese Spring), which was the meta-QTL 1 of chromosome 3B. This meta-QTL presented the narrowest length in genetic distance among all 65 meta-QTL verified, i.e., 0.82 cM, which corresponded to 18.2 Mb in physical distance.
Candidate gene mining within meta-QTL 1/3B, using the annotation of wheat reference genome sequence of Chinese Spring, resulted in a list of 324 genes (Supplementary File 4). These genes were investigated regarding expression under control and FHB occurrence, using publicly available data (Kugler et al., 2013). A total of 10 genes were found as differentially expressed in the resistant variety CM-82036 (Figure 10). None of these genes were differentially expressed in the susceptible line assayed. From these genes, five were differentially expressed only at 30 h after inoculation, four only at 50 h after inoculation, and one at both periods. Also, five genes were up-regulated when challenged by the pathogen, four down-regulated, and one up-regulated at 30 h and down-regulated at 50 h after inoculation.
Transcriptional expression of candidate genes mined within the meta-QTL 1/3B, using publically available transcriptional data of wheat spike tissues (Kugler et al., 2013), after Fusarium graminearum inoculation and control (mock) conditions in the resistant variety CM-82036. TPM, transcripts per million.
The genes were further searched for functional annotation (protein evidence) analysis (Table 3). Five genes were only found as either uncharacterized, putative, or predicted proteins. From the five remaining genes, the first has been annotated as a Sarcoplasmic reticulum histidine-rich calcium-binding protein, in Zea mays, the second as a Glycosyltransferase, again in Z. mays, the third a Cytochrome P450 in Aegilops tauschii, Oryza sativa subsp. japonica, and Triticum Urartu and a Transposon protein in Z. mays, the fourth a Nuclease PA3, in Z. mays, and finally, the last a Metal tolerance protein, in O. sativa subsp. japonica, Triticum Urartu, and Z. mays.
Protein evidence for the candidate genes mined within Meta-QTL 1 of chromosome 3B and differentially expressed after Fusarium graminearum inoculation in the variety CM-82036 (Kugler et al., 2013).
Gene ID
Annotation at IWGSC—Protein evidence from other species
TraesCS3B02G019900
Sarcoplasmic reticulum histidine-rich calcium-binding protein (Zea mays).
TraesCS3B02G021100
Annotated only as uncharacterized, predicted or putative protein.
Metal tolerance protein (Oryza sativa subsp. japonica, Triticum urartu, Zea mays).
TraesCS3B02G041800
Annotated only as uncharacterized protein.
Discussion
During the last two decades, starting with the first studies published by Bai et al. (1999) and Waldron et al. (1999) until today, several papers have been published on QTL mapping for FHB resistance in bread wheat (Steiner et al., 2017). It has motivated the execution of the present work, in which we attempted to re-analyze all these loci jointly, aiming to refine this information in order to be better used by breeders and also for the advance of the basic research. To the best of our knowledge, this is the largest collection and meta-analysis of QTL related to FHB resistance in bread wheat performed so far, and one of the largest of this kind already carried out for this cereal. The totality of QTL here analyzed can receive the term, recently proposed, of “QTLome” (Salvi and Tuberosa, 2015; Martinez et al., 2016) of FHB resistance in wheat, as it assembles in a global analysis most of the more important loci mapped in the crop until now. We have here updated the QTL-meta analysis previously performed by Liu et al. (2009), Löffler et al. (2009), and Mao et al. (2010), but more than that, we have applied the current improved algorithms and software, not available at the time these studies were carried out. Moreover, we have integrated genomic and transcriptomic data, which were also a recent resource for wheat research.
A total of 556 QTL were collected from the literature as being responsible for FHB resistance in wheat. This is a considerably large number; however, some level of redundancy was expected within this group of loci, which is why the meta-analysis was necessary. A total of 323 QTL could be included in the meta-analysis and generated 65 meta-QTL, i.e., regions statistically validated as unique. Although the redundancy was decreased by roughly five times, the resulting number of refined loci can still be considered as large. Just to compare, another recent study of this kind performed on hexaploid wheat has found a total of 32 meta-QTL, but by combining yield, baking quality and grain protein traits (Quraishi et al., 2017). The large number of meta-QTL obtained in this study emphasizes the wealth of genetic variability present in wheat gene pool for FHB resistance, as already recognized by other researchers (Steiner et al., 2017). Interestingly, it does not agree with the fact that hexaploid wheat has been considered as a species with narrow genetic variability, an issue that has supposedly been aggravated by modern breeding practices (Fu and Somers, 2009; Charmet, 2011).
The refinement obtained in this study was variable among loci. In some situations, such as for meta-QTL 1/3B, a total of 36 QTL (whole QTL or parts) were clustered in one single locus, which is a high reduction of redundancy. On the other hand, meta-QTL 5/3B was the result of the analysis of only one original QTL. Furthermore, in six linkage groups, the meta-analysis could not be performed, due to the insufficient number of QTL projected. For these under-represented loci, more QTL have to be added to the analysis to validate these regions as responsible for the trait. Even after several years of research in the search of QTL related to FHB resistance in wheat, the reports continue in an elevated rate (Steiner et al., 2017). Just to exemplify, during the final writing of this paper, Silva et al. (2019) have mapped another four QTL for FHB resistance and DON content. However, a decisive improvement in QTL meta-analysis would be the integration of data also from genome wide association studies (GWAS), however the available methods and algorithms still do not deal with such. Through GWAS, several genotypes are assayed, so several loci can be mapped in a single study. Also, the method allows for more precision in QTL mapping, which, if joined in a meta-analysis with QTL obtained from biparental mapping population, could harvest more refined meta-QTL. In fact, important GWAS have been carried out for FHB resistance in wheat (Kollers et al., 2013; Arruda et al., 2016).
When collecting the data to carry out this work, information about QTL conferring the five different types of FHB resistance were collected. All these QTL, which represent, in fact, distinct traits, were analyzed simultaneously, as responsible by the single “meta-trait” FHB resistance. It can be verified that the majority of the meta-QTL generated throughout this study were composed by the projection of QTL related to different types of resistance to this disease. Although this classification of types of FHB resistance is broadly accepted among researchers (Steiner et al., 2017), it is also acceptable that both similar mechanisms could play a role in each type of resistance as well as that within each of these mapped QTL, therefore genes for different traits can exist. In fact, the most important locus for FHB resistance already studied, Fhb1, has been confirmed as conferring resistance to the spread of the disease (type II resistance) as well as to the mycotoxin accumulation (type III resistance) (Lemmens et al., 2005; Gunupuru et al., 2018).
Meta-QTL are possibly genomic regions highly rich in genes. For maize, it was already deeply elucidated, in which the QTLome for yield potential presented a high correlation with gene density in the genome of the species (Martinez et al., 2016). In rice, within each meta-QTL detected for drought tolerance, candidate genes were mined, through comparative genomic approaches, as most of these genes were already suggested as being related to the trait, according to previous studies (Swamy et al., 2011). For wheat, within 32 meta-QTL related to yield, baking quality and grain protein content, a total of 15,772 genes were listed, but only 37 were considered as major candidates (Quraishi et al., 2017). In the present study, gene mining was only carried out in 1 meta-QTL (meta-QTL 1/3B), in which 324 candidate genes were listed, with only 10 cross validated as differentially expressed in CM-82036 under F. graminearum infection.
The most important FHB resistance locus studied so far, namely the QTL Fhb1, firstly mapped in the Chinese cultivar Sumai 3, is localized in the short arm of the 3B chromosome (Cuthbert et al., 2006; Bernardo et al., 2012; Gunnaiah et al., 2012; Rawat et al., 2016). In this study, the genetic position of this locus corresponds to meta-QTL 1 of the chromosome 3B. Although still in dispute, a recent study has suggested that Fhb1 from Sumai 3, harbors as the main candidate for the trait expression, a gene which encodes a chimeric lectin with agglutinin domains and a pore-forming toxin-like domain (Rawat et al., 2016). In this study, it was not expected to find this gene within meta-QTL 1/3B, since we have used the sequence and annotation of Chinese Spring as a reference genome, which has been proven not to harbor an active copy of the gene (Rawat et al., 2016; Schweiger et al., 2016). Thus, since the gene evidenced by Rawat et al. (2016) is not present in meta-QTL 1/3B, all other genes found as differentially expressed within this meta-QTL could also be suggested as promising candidate genes, even for minor effect. It is important to mention that in order to examine the transcriptional expression of the listed candidate genes, we have used the available data for CM-82036, a variety derived from Sumai-3 and which harbors Fhb1. Both Sumai-3 and CM-82036 were part of the donors of QTL, which composed the analysis carried out in this study.
Five out of the 10 candidate genes which showed differential expression when CM-82036 was challenged by F. graminearum are annotated only as uncharacterized, putative or predicted proteins. These genes deserve further studies, in order to unveil their possible role in FHB resistance and future use in breeding programs. Regarding the differentially expressed genes annotated as Sarcoplasmic reticulum histidine-rich calcium-binding protein and metal tolerance protein, these genes have also been verified in other transcriptome studies, which also assayed Sumai-3 derived lines (Hofstad et al., 2016; Schweiger et al., 2016). To the best of our knowledge, nuclease PA3 and transposon proteins have not been attributed to FHB resistance mechanisms in wheat so far.
Glycosyltransferases play a vast number of roles in living organisms and have already been attributed to FHB resistance mechanisms in wheat, especially the uridinediphosphate glycosyltransferases (UDP-Glycosiltransferases) (Li et al., 2015; Xin et al., 2015; He et al., 2018). The main role of this molecule would be the Deoxynivalenol (DON) inactivation, and DON is a F. graminearum pathogenicity factor. In the transcriptome dataset we have explored, this gene was up-regulated when CM-82036 was challenged by the pathogen (Kugler et al., 2013). Regarding Cytochrome P450, previous studies have also shown the role of these proteins on wheat defense against Fusarium infection and DON accumulation. In Sumai 3 and CM-82036, the higher expression of this gene has been considered as a resistance mechanism, potentially involved in the biosynthesis of the wheat defense compound, DIMBOA (2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one), which is known to be toxic to F. graminearum (Li et al., 2010; revised by Gunupuru et al., 2018; Kazan and Gardiner (2018). However, for the available data we have analyzed, the present gene had its expression down-regulated under Fusarium infection (at 50 h after inoculation), which implies that more studies are needed to understand the specific mode-of-action of this gene. In fact, the advance in the identification of the genes behind QTL and the understanding of their mechanisms for FHB resistance in wheat are still in the initial steps, however, important projects are underway (Steiner et al., 2017).
In this study we have applied an entirely in silico approach aiming to unveil the genetic architecture of FHB resistance in bread wheat, joining the available QTL data and genomic resources of the cereal. Successfully, we could assembly an important part of the promised “whole puzzle” of FHB resistance, from QTL to differentially expressed genes, using only publicly available data. It is important to emphasize that after any in silico investigation of this kind, the confirmation of the involvement of any of these loci, especially the candidate genes listed, in FHB resistance, can only be performed through further methods, such as gene cloning, gene silencing, and even other more recent methods, such as the metabolo-transcriptomics/proteomics that have already been performed for Fhb1 and Fhb2 loci (Gunnaiah et al., 2012; Dhokane et al., 2016; Rawat et al., 2016). The results obtained through this study can support the improvement in breeding strategies for FHB resistance, through biotechnological techniques such as gene-editing of target genes (Bortesi and Fischer, 2015). Our findings can also help to drive a more precise marker assisted selection for these loci. Finally, these results might contribute for the better understanding of the mechanisms related to FHB resistance in wheat.
Author Contributions
AC, EV, RS, and DF contributed for the conception and design of the study. EV organized the database and performed all data analysis and wrote the first draft of the manuscript. VR contributed for the database organization. EV, RS, LM, CP, and AC wrote sections of the manuscript. AC performed the final revision of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.
Conflict of Interest Statement
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.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00727/full#supplementary-material
Database (QTL data) collected for the study.
Maps and QTL files prepared to run the meta-analysis performed through this study.
Meta-QTL generated positioned onto the consensus map.
Candidate genes mined within meta-QTL 1/3B.
ReferencesAbateZ. A.LiuS.McKendryA. L. (2008). Quantitative trait loci associated with deoxynivalenol content and kernel quality in the soft red winter wheat ‘Ernie’. 48, 1408–1418. 10.2135/cropsci2007.07.0411AkbariM.WenzlP.CaigV.CarlingJ.XiaL.YangS.. (2006). Diversity arrays technology (DArT) for high-throughput profiling of the hexaploid wheat genome. 113, 1409–1420. 10.1007/s00122-006-0365-417033786AndersonJ. A.StackR. W.LiuS.WaldronB. L.FjeldA. D.CoyneC.. (2001). DNA markers for Fusarium head blight resistance QTL in two wheat populations. 102, 1164–1168. 10.1007/s001220000509ArcadeA.LabourdetteA.FalqueM.ManginB.ChardonF.CharcossetA.. (2004). BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. 20, 2324–2326. 10.1093/bioinformatics/bth23015059820ArrudaM. P.BrownP.Brown-GuediraG.KrillA. M.ThurberC.MerrillK. R.. (2016). Genome-wide association mapping of Fusarium head blight resistance in wheat using genotyping-by-sequencing. 9, 1–14. 10.3835/plantgenome2015.04.002827898754BaiG.ShanerG. (2004). Management and resistance in wheat and barley to Fusarium head blight. 42, 135–161. 10.1146/annurev.phyto.42.040803.14034015283663BaiG. H.KolbF. L.ShanerG.DomierL. L. (1999). Amplified fragment length polymorphism markers linked to a major quantitative trait locus controlling scab resistance in wheat. 89, 343–348. 10.1094/PHYTO.1999.89.4.34318944781BasnetB. R.GloverK. D.IbrahimA. M.YenY.ChaoS. (2012). A QTL on chromosome 2DS of ‘Sumai 3’ increases susceptibility to Fusarium head blight in wheat. 186, 91–101. 10.1007/s10681-011-0495-xBernardoA. N.MaH.ZhangD.BaiG. (2012). Single nucleotide polymorphism in wheat chromosome region harboring Fhb1 for Fusarium head blight resistance. 29, 477–488. 10.1007/s11032-011-9565-yBoninC. M.KolbF. L. (2009). Resistance to Fusarium head blight and kernel damage in a winter wheat recombinant inbred line population. 49, 1304–1312. 10.2135/cropsci2008.08.0459BorrillP.Ramirez-GonzalezR.UauyC. (2016). expVIP: a customisable RNA-seq data analysis and visualisation platform. . 170, 2172–86. 10.1104/pp.15.01667BortesiL.FischerR. (2015). The CRISPR/Cas9 system for plant genome editing and beyond. 33, 41–52. 10.1016/j.biotechadv.2014.12.00625536441BourdoncleW.OhmH. W. (2003). Quantitative trait loci for resistance to Fusarium head blight in recombinant inbred wheat lines from the cross Huapei 57-2/Patterson. 131, 131–136. 10.1023/A:1023056207513BuerstmayrH.BanT.AndersonJ. A. (2009). QTL mapping and marker-assisted selection for Fusarium head blight resistance in wheat: a review. 128, 1–26. 10.1111/j.1439-0523.2008.01550.xBuerstmayrH.BuerstmayrM.SchweigerW.SteinerB. (2014). Breeding for resistance to head blight caused by Fusarium spp. in wheat. 9, 1–13. 10.1079/PAVSNNR20149007BuerstmayrH.LemmensM. (2015). Breeding healthy cereals: genetic improvement of Fusarium resistance and consequences for mycotoxins. 8, 591–602. 10.3920/WMJ2015.1889BuerstmayrH.LemmensM.HartlL.DoldiL.SteinerB.StierschneiderM.. (2002). Molecular mapping of QTL for Fusarium head blight resistance in spring wheat. I. Resistance to fungal spread (Type II resistance). 104, 84–91. 10.1007/s001220200009BuerstmayrH.SteinerB.HartlL.GriesserM.AngererN.LengauerD.. (2003). Molecular mapping of QTL for Fusarium head blight resistance in spring wheat. II. Resistance to fungal penetration and spread. 107, 503–508. 10.1007/s00122-003-1272-6BuerstmayrM.BuerstmayrH. (2015). Comparative mapping of quantitative trait loci for Fusarium head blight resistance and anther retention in the winter wheat population Capo × Arina. 128, 1519–1530. 10.1007/s00122-015-2527-825982129CaiJ.BaiG. (2014). Quantitative trait loci for Fusarium head blight resistance in Huangcandou × ‘Jagger’ wheat population. 54, 2520–2528. 10.2135/cropsci2013.12.0835CativelliM.LewisS.AppendinoM. L. (2013). A Fusarium head blight resistance quantitative trait locus on chromosome 7D of the spring wheat cultivar Catbird. 53, 1464–1471. 10.2135/cropsci2012.07.0435CavanaghC. R.ChaoS.WangS.HuangB. E.StephenS.KianiS.. (2013). Genome-wide comparative diversity uncovers multiple targets of selection for improvement in hexaploid wheat landraces and cultivars. 110, 8057–8062. 10.1073/pnas.121713311023630259ChaoL.SongY.GaoL.YaoQ.ZhouR.XuR.. (2014). Integration of QTL detection and marker assisted selection for improving resistance to Fusarium head blight and important agronomic traits in wheat. 2, 70–78. 10.1016/j.cj.2013.10.004CharmetG. (2011). Wheat domestication: lessons for the future. 334, 212–220. 10.1016/j.crvi.2010.12.01321377616ChenJ.GriffeyC. A.MaroofS.StrombergE. L.BiyashevR. M.ZhaoW.. (2006). Validation of two major quantitative trait loci for Fusarium head blight resistance in Chinese wheat line W14. 125, 99–101. 10.1111/j.1439-0523.2006.01182.xChuC.NiuZ.ZhongS.ChaoS.FriesenT. L.HalleyS.. (2011). Identification and molecular mapping of two QTL with major effects for resistance to Fusarium head blight in wheat. 123:1107. 10.1007/s00122-011-1652-2CuthbertP. A.SomersD. J.Brulé-BabelA. (2007). Mapping of Fhb2 on chromosome 6BS: a gene controlling Fusarium head blight field resistance in bread wheat (Triticum aestivum L.). 114, 429–437. 10.1007/s00122-006-0439-317091262CuthbertP. A.SomersD. J.ThomasJ.CloutierS.Brulé-BabelA. (2006). Fine mapping Fhb 1, a major gene controlling Fusarium head blight resistance in bread wheat (Triticum aestivum L.). 112:1465. 10.1007/s00122-006-0249-7DarvasiA.SollerM. (1997). A simple method to calculate resolving power and confidence interval of QTL map location. 27, 125–132. 10.1023/A:10256853248309145551Del PonteE. M.FernandesJ. M. C.PierobomC. R.BergstromG. C. (2004). Giberela do Trigo – Aspectos Epidemiológicos e Modelos de Previsão. 29, 587–605. 10.1590/S0100-41582004000600001DhokaneD.KarreS.KushalappaA. C.McCartneyC. (2016). Integrated metabolo-transcriptomics reveals Fusarium head blight candidate resistance genes in wheat QTL-Fhb2. 11:0155851. 10.1371/journal.pone.015585127232496DraegerR.GosmanN.SteedA.ChandlerE.ThomsettM.SchondelmaierJ.. (2007). Identification of QTL for resistance to Fusarium head blight, DON accumulation and associated traits in the winter wheat variety Arina. 115, 617–625. 10.1007/s00122-007-0592-3EndelmanJ. B.PlomionC. (2014). LPmerge: an R package for merging genetic maps by linear programming. 30, 1623–1624. 10.1093/bioinformatics/btu09124532720FuY.-B.SomersD. J. (2009). Genome-wide reduction of genetic diversity in wheat breeding. 49, 161–168. 10.2135/cropsci2008.03.0125GervaisL.DedryverF.MorlaisJ. Y.BodusseauV.NegreS.BilousM.. (2003). Mapping of quantitative trait loci for field resistance to Fusarium head blight in a European winter wheat. 106, 961–970. 10.1007/s00122-002-1160-512671743GilsingerJ.KongL.ShenX.OhmH. (2005). DNA markers associated with low Fusarium head blight incidence and narrow flower opening in wheat. 110, 1218–1225. 10.1007/s00122-005-1953-415750825GoffinetB.GerberS. (2000). Quantitative trait loci: a meta-analysis. 155, 463–473.10790417GunnaiahR.KushalappaA. C.DuggavathiR.FoxS.SomersD. J. (2012). Integrated metabolo-proteomic approach to decipher the mechanisms by which wheat QTL (Fhb1) contributes to resistance against Fusarium graminearum. 7:40695. 10.1371/journal.pone.004069522866179GunupuruL. R.ArunachalamC.MallaK. B.KahlaA.PerochonA.JiaJ.. (2018). A wheat cytochrome P450 enhances both resistance to deoxynivalenol and grain yield. 13:0204992. 10.1371/journal.pone.020499230312356GuoB.SleperD. A.LuP.ShannonJ. G.NguyenH. T.ArelliP. R. (2006). QTL associated with resistance to soybean cyst nematode in soybean: meta-analysis of QTL locations. 46, 595–602. 10.2135/cropsci2005.04-0036-2GuptaA.LippsP. E.CampbellK. G. (2000). Finding quantitative trait loci associated with Fusarium head blight of wheat using simple sequence repeat markers, in , eds CantyS. M.LewisJ.SilerL.WardR. W. (Erlanger, KY; East Lansing, MI: Michigan State University, 28–32.GuptaA.LippsP. E.CampbellK. G.SnellerC. H. (2001). Identification of QTL associated with resistance to FHB in Ning7840 and freedom, in , eds CantyS. M.LewisJ.SilerL.WardR. W. (Erlanger, KY; East Lansing, MI: Michigan State University), 180.HäberleJ.SchmolkeM.SchweizerG.KorzunV.EbmeyerE.ZimmermannG.. (2007). Effects of two major Fusarium head blight resistance QTL verified in a winter wheat backcross population. 47, 1823–1831. 10.2135/cropsci2006.12.0800HäberleJ.SchweizerG.SchondelmaierJ.ZimmermannG.HartlL. (2009). Mapping of QTL for resistance against Fusarium head blight in the winter wheat population Pelikan//Bussard/Ning8026. 128, 27–35. 10.1111/j.1439-0523.2008.01540.xHandaH.NamikiN.XuD.BanT. (2008). Dissecting of the FHB resistance QTL on the short arm of wheat chromosome 2D using a comparative genomic approach: from QTL to candidate gene. 22, 71–84. 10.1007/s11032-008-9157-7HeX.LillemoM.ShiJ.WuJ.BjornstadA.BelovaT.. (2016). QTL characterization of Fusarium head blight resistance in CIMMYT bread wheat line Soru# 1. 11:0158052. 10.1371/journal.pone.0158052HeY.AhmadD.ZhangX.ZhangY.WuL.JiangP.. (2018). Genome-wide analysis of family-1 UDP glycosyltransferases (UGT) and identification of UGT genes for FHB resistance in wheat (Triticum aestivum L.). 18:67. 10.1186/s12870-018-1286-529673318HofstadA. N.NussbaumerT.AkhunovE.ShinS.KuglerK. G.KistlerH. C.. (2016). Examining the transcriptional response in wheat Fhb1 near-isogenic lines to fusarium graminearum infection and deoxynivalenol treatment. 9. 10.3835/plantgenome2015.05.0032HolzapfelJ.VossH. H.MiedanerT.KorzunV.HäberleJ.SchweizerG.. (2008). Inheritance of resistance to Fusarium head blight in three European winter wheat populations. 117, 1119–1128. 10.1007/s00122-008-0850-z18670751International Wheat Genome Sequencing Consortium (2014). A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. 345:1251788. 10.1126/science.1251788International Wheat Genome Sequencing Consortium (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome. 17, 361. 10.1126/science.aar7191IslamM. S.Brown-GuediraG.Van SanfordD.OhmH.DongY.McKendryA. L. (2016). Novel QTL associated with the Fusarium head blight resistance in Truman soft red winter wheat. 207, 571–592. 10.1007/s10681-015-1550-9JayatilakeD. V.BaiG. H.DongY. H. (2011). A novel quantitative trait locus for Fusarium head blight resistance in chromosome 7A of wheat. 122, 1189–1198. 10.1007/s00122-010-1523-221221526JiaG.ChenP.QinG.BaiG.WangX.WangS.. (2005). QTL for Fusarium head blight response in a wheat DH population of Wangshuibai/Alondra‘s’. 146, 183–191. 10.1007/s10681-005-9001-7JiangG. L.DongY.ShiJ.WardR. W. (2007a). QTL analysis of resistance to Fusarium head blight in the novel wheat germplasm CJ 9306. II. Resistance to deoxynivalenol accumulation and grain yield loss. 115, 1043–1052. 10.1007/s00122-007-0630-117726598JiangG. L.ShiJ.WardR. W. (2007b). QTL analysis of resistance to Fusarium head blight in the novel wheat germplasm CJ 9306. I. Resistance to fungal spread. 116, 3–13. 10.1007/s00122-007-0641-y17898987KazanK.GardinerD. M. (2018). Transcriptomics of cereal–Fusarium graminearum interactions: what we have learned so far. 19, 764–778. 10.1111/mpp.1256128411402KlahrA.ZimmermannG.WenzelG.MohlerV. (2007). Effects of environment, disease progress, plant height and heading date on the detection of QTL for resistance to Fusarium head blight in an European winter wheat cross. Euphytica154, 17–28. 10.1007/s10681-006-9264-7KollersS.RodemannB.LingJ.KorzunV.EbmeyerE.ArgillierO.. (2013). Whole genome association mapping of Fusarium head blight resistance in European winter wheat (Triticum aestivum L.). 8:57500. 10.1371/journal.pone.005750023451238KuglerK. G.SiegwartG.NussbaumerT.AmetzC.SpannaglM.SteinerB.. (2013). Quantitative trait loci-dependent analysis of a gene co-expression network associated with Fusarium head blight resistance in bread wheat (Triticum aestivum L.). 14:728. 10.1186/1471-2164-14-72824152241LemmensM.ScholzU.BerthillerF.Dall'astaC.KoutnikA.SchuhmacherR.. (2005). The ability to detoxify the mycotoxin deoxynivalenol colocalizes with a major quantitative trait locus for Fusarium head blight resistance in wheat. 18, 1318–1324. 10.1094/MPMI-18-131816478051LiC.ZhuH.ZhangC.LinF.XueS.CaoY.. (2007). Mapping QTL associated with Fusarium-damaged kernels in the Nanda 2419 × Wangshuibai population. 163, 185–191. 10.1007/s10681-007-9626-9LiT.BaiG.WuS.GuS. (2011). Quantitative trait loci for resistance to Fusarium head blight in a Chinese wheat landrace Haiyanzhong. 122, 1497–1502. 10.1007/s00122-011-1549-021344182LiT.BaiG.WuS.GuS. (2012). Quantitative trait loci for resistance to Fusarium head blight in the Chinese wheat landrace Huangfangzhu. 185, 93–102. 10.1007/s10681-012-0631-2LiX.ShinS.HeinenS.Dill-MackyR.BerthillerF.NersesianN.. (2015). Transgenic wheat expressing a barley UDP-glucosyltransferase detoxifies deoxynivalenol and provides high levels of resistance to Fusarium graminearum. 28, 1237–1246. 10.1094/MPMI-03-15-0062-R26214711LiX.ZhangJ. B.SongB.LiH. P.XuH. Q.QuB.. (2010). Resistance to Fusarium head blight and seedling blight in wheat is associated with activation of a cytochrome P450 gene. 100, 183–191. 10.1094/PHYTO-100-2-018320055652LinF.KongZ. X.ZhuH. L.XueS. L.WuJ. Z.TianD. G.. (2004). Mapping QTL associated with resistance to Fusarium head blight in the Nanda 2419 × Wangshuibai population. I. Type II resistance. 109, 1504–1511. 10.1007/s00122-004-1772-zLinF.XueS. L.ZhangZ. Z.ZhangC. Q.KongZ. X.YaoG. Q.. (2006). Mapping QTL associated with resistance to Fusarium head blight in the Nanda 2419 × Wangshuibai population. II: Type I resistance. 112, 528–535. 10.1007/s00122-005-0156-3LiuS.AbateZ. A.LuH.MusketT.DavisG. L.MckendryA. L. (2007). QTL associated with Fusarium head blight resistance in the soft red winter wheat Ernie. 115, 417–427. 10.1007/s00122-007-0577-217569028LiuS.ChristopherM. D.GriffeyC. A.HallM. D.GundrumP. G.BrooksW. S. (2012). Molecular characterization of resistance to Fusarium head blight in US soft red winter wheat breeding line VA00W-38. 52, 2283–2292. 10.2135/cropsci2012.03.0144LiuS.HallM. D.GriffeyC. A.MckendryA. L. (2009). Meta-analysis of QTL associated with Fusarium head blight resistance in wheat. 49, 1955–1968. 10.2135/cropsci2009.03.0115LiuS.ZhangX.PumphreyM. O.StackR. W.GillB. S.AndersonJ. 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. 6, 83–89. 10.1007/s10142-005-0007-y16270217LöfflerM.SchönC. C.MiedanerT. (2009). Revealing the genetic architecture of FHB resistance in hexaploid wheat (Triticum aestivum L.) by QTL meta-analysis. 23, 473–488. 10.1007/s11032-008-9250-yLuQ.LillemoM.SkinnesH.HeX.ShiJ.JiF.. (2012). Anther extrusion and plant height are associated with Type I resistance to Fusarium head blight in bread wheat line ‘Shanghai-3/Catbird’. 126, 317–334. 10.1007/s00122-012-1981-923052019MaH. X.BaiG. H.ZhangX.LuW. Z. (2006a). Main effects, epistasis, and environmental interactions of quantitative trait loci for Fusarium head blight resistance in a recombinant inbred population. 96, 534–541. 10.1094/PHYTO-96-0534MaH. X.ZhangK. M.GaoL.BaiG. H.ChenH. G.CaiZ. X.. (2006b). Quantitative trait loci for resistance to Fusarium head blight and deoxynivalenol accumulation in Wangshuibai wheat under field conditions. 55, 739–745. 10.1111/j.1365-3059.2006.01447.xMallaS.IbrahimA. M. H.BerzonskyW.YenY. (2008a). Inheritance of FHB resistance in spring versus winter wheat growth habit backgrounds, in , eds CantyS. M.ClarkA.WaltonE.EllisD.MundellJ.Van SanfordD. (Lexington, KY: University of Kentucky), 182.MallaS.IbrahimA. M. H.BerzonskyW.YenY. (2008b). Mapping QTLs for Fusarium head blight from novel source–Tokai 66, in , eds CantyS. M.ClarkA.WaltonE.EllisD.MundellJ.Van SanfordD. (Lexington, KY: University of Kentucky), 183.MallaS.IbrahimA. M. H.YenY.BerzonskyW.GloverK. D.SteinJ. (2010). Quantitative Trait Loci analysis of novel Fusarium head blight resistance in Tokai 66. 5, 62–69. 10.3844/ajabssp.2010.62.69MaoS. L.WeiY. M.CaoW.LanX. J.YuM.ChenZ. M.. (2010). Confirmation of the relationship between plant height and Fusarium head blight resistance in wheat (Triticum aestivum L.) by QTL meta-analysis. 174, 343–356. 10.1007/s10681-010-0128-9MardiM.BuerstmayrH.GhareyazieB.LemmensM.MohammadiS. A.NolzR.. (2005). QTL analysis of resistance to Fusarium head blight in wheat using a ‘Wangshuibai’-derived population. 124, 329–333. 10.1111/j.1439-0523.2005.01103.xMardiM.PazoukiL.DelavarH.KazemiM. B.GhareyazieB.SteinerB.. (2006). QTL analysis of resistance to Fusarium head blight in wheat using a ‘Frontana’-derived population. 125, 313–317. 10.1111/j.1439-0523.2006.01228.xMartinezA. K.SorianoJ. M.TuberosaR.KoumproglouR.JahrmannT.SalviS. (2016). Yield QTLome distribution correlates with gene density in maize. 242, 300–309. 10.1016/j.plantsci.2015.09.02226566847McCartneyC. A.Brûlé-BabelA. L.FedakG.MartinR. A.MccallumB. D.GilbertJ.. (2016). Fusarium head blight resistance QTL in the spring wheat cross Kenyon/86ISMN 2137. 7, 1–11. 10.3389/fmicb.2016.0154227790188MesterhazyA. (1995). Types and components of resistance to Fusarium head blight of wheat. 114, 377–386. 10.1111/j.1439-0523.1995.tb00816.xMesterhazyA.BartokT.MirochaC. G.KomorozyR. (1999). Nature of wheat resistance to Fusarium head blight and role of Deoxynivalenol for breeding. 118, 7–110. 10.1046/j.1439-0523.1999.118002097.xMillerJ.ArnisonP. (1986). Degradation of de-oxynivalenol by suspension culture of the Fusarium head blight resistant wheat cultivar Frontana. 8, 147–150. 10.1080/07060668609501818PaillardS.SchnurbuschT.TiwariR.MessmerM.WinzelerM.KellerB.. (2004). QTL analysis of resistance to Fusarium head blight in Swiss winter wheat (Triticum aestivum L.). 109, 323–332. 10.1007/s00122-004-1628-615014875ParryD. W.JenkinsonP.McleodL. (1995). Fusarium ear blight (scab) in small grain cereals—a review. 44, 207–238. 10.1111/j.1365-3059.1995.tb02773.xPestkaJ. (2010). Toxicological mechanisms and potential health effects of deoxynivalenol and nivalenol. 3, 323–347. 10.3920/WMJ2010.1247PetersenS.LyerlyJ. H.MaloneyP. V.Brown-GuediraG.CowgerC.CostaJ. M.. (2016). Mapping of head blight resistance quantitative trait loci in winter wheat cultivar NC-neuse. 56, 1473–1483. 10.2135/cropsci2015.05.0312PetersenS.LyerlyJ. H.MckendryA. L.IslamM. S.Brown-GuediraG.CowgerC.. (2017). Validation of Fusarium head blight resistance QTL in US winter wheat. 57, 1–12. 10.2135/cropsci2015.07.0415PolandJ.EndelmanJ.DawsonJ.RutkoskiJ.WuS.ManesY.. (2012). Genomic selection in wheat breeding using genotyping-by-sequencing. 5, 103–113. 10.3835/plantgenome2012.06.0006PolandJ. A.BrownP. J.SorrellsM. E.JanninkJ. L. (2012). Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach. 7:32253. 10.1371/journal.pone.003225322389690QuraishiU. M.PontC.AinQ. U.FloresR.BurlotL.AlauxM.. (2017). Combined genomic and genetic data integration of major agronomical traits in bread wheat (Triticum aestivum L.). 8:1843. 10.3389/fpls.2017.0184329184557Ramírez-GonzálezR. H.BorrillP.LangD.HarringtonS. A.BrintonJ.VenturiniL.. (2018). The transcriptional landscape of hexaploid wheat across tissues and cultivars. 361:eaar6089. 10.1126/science.aar6089RawatN.PumphreyM. O.LiuS.ZhangX.TiwariV. K.AndoK.. (2016). Wheat Fhb1 encodes a chimeric lectin with agglutinin domains and a pore-forming toxin-like domain conferring resistance to Fusarium head blight. 48, 1576–1580. 10.1038/ng.370627776114SaidJ. I.SongM.WangH.LinZ.ZhangX.FangD. D.. (2015). A comparative meta-analysis of QTL between intraspecific Gossypium hirsutum and interspecific G. hirsutum × G. barbadense populations. 290, 1003–1025. 10.1007/s00438-014-0963-925501533SaintenacC.JiangD.WangS.AkhunovE. (2013). Sequence-based mapping of the polyploid wheat genome. 3, 1105–1114. 10.1534/g3.113.00581923665877SalviS.TuberosaR. (2015). The crop QTLome comes of age. 32, 179–185. 10.1016/j.copbio.2015.01.00125614069SchmolkeM.ZimmermannG.BuerstmayrH.SchweizerG.MiedanerT.KorzunV.. (2005). Molecular mapping of Fusarium head blight resistance in the winter wheat population Dream/Lynx. 111, 747–756. 10.1007/s00122-005-2060-215947905SchroederH. W.ChristensenJ. J. (1963). Factors affecting resistance of wheat to scab caused by Gibberella zeae. 53, 831–838.SchweigerW.SteinerB.VautrinS.NussbaumerT.SiegwartG.ZaminiM.. (2016). Suppressed recombination and unique candidate genes in the divergent haplotype encoding Fhb1, a major Fusarium head blight resistance locus in wheat. 129, 1607–1623. 10.1007/s00122-016-2727-x27174222SemagnK.BabuR.HearneS.OlsenM. (2014). Single nucleotide polymorphism genotyping using Kompetitive Allele Specific PCR (KASP): overview of the technology and its application in crop improvement. 33, 1–4. 10.1007/s11032-013-9917-xSemagnK.SkinnesH.BjørnstadÅ.MarøyA. G.TarkegneY. (2007). Quantitative trait loci controlling Fusarium head blight resistance and low deoxynivalenol content in hexaploid wheat population from ‘Arina’ and NK93604. 47, 294–303. 10.2135/cropsci2006.02.0095ShenX.IttuM.OhmH. W. (2003a). Quantitative trait loci conditioning resistance to Fusarium head blight in wheat line F201R. 43, 850–857. 10.2135/cropsci2003.0850ShenX.ZhouM.LuW.OhmH. (2003b). Detection of Fusarium head blight resistance QTL in a wheat population using bulked segregant analysis. 106, 1041–1047. 10.1007/s00122-002-1133-812671752SilvaC. L.FritzA.ClinesmithM.PolandJ.DowellF.PeirisK. (2019). QTL mapping of Fusarium head blight resistance and deoxynivalenol accumulation in the Kansas wheat variety ‘Everest’. 39, 1–21. 10.1007/s11032-019-0937-zSomersD. J.FedakG.SavardM. (2003). Molecular mapping of novel genes controlling Fusarium head blight resistance and deoxynivalenol accumulation in spring wheat. 46, 555–564. 10.1139/g03-03312897863SosnowskiO.CharcossetA.JoetsJ. (2012). BioMercator V3: an upgrade of genetic map compilation and quantitative trait loci meta-analysis algorithms. 28, 2082–2083. 10.1093/bioinformatics/bts31322661647Srinivasachary GosmanN.SteedA.HollinsT. W.BaylesR.JenningsP. (2009). Semi-dwarfing Rht-B1 and Rht-D1 loci of wheat differ significantly in their influence on resistance to Fusarium head blight. 118:695. 10.1007/s00122-008-0930-019034409SrinivasacharyGosman, N.SteedA.SimmondsJ.Leverington-WaiteM.WangY.. (2008). Susceptibility to Fusarium head blight is associated with the Rht-D1b semi-dwarfing allele in wheat. 116, 1145–1153. 10.1007/s00122-008-0742-218347773SteinerB.BuerstmayrM.MichelS.SchweigerW.LemmensM.BuerstmayrH. (2017). Breeding strategies and advances in line selection for Fusarium head blight resistance in wheat. 42, 165–174. 10.1007/s40858-017-0127-7SteinerB.LemmensM.GriesserM.ScholzU.SchondelmaierJ.BuerstmayrH. (2004). Molecular mapping of resistance to Fusarium head blight in the spring wheat cultivar Frontana. 109, 215–224. 10.1007/s00122-004-1620-114997302SuzukiT.SatoM.TakeuchiT. (2012). Evaluation of the effects of five QTL regions on Fusarium head blight resistance and agronomic traits in spring wheat (Triticum aestivum L.). 62, 11–17. 10.1270/jsbbs.62.1123136509SwamyB. M.VikramP.DixitS.AhmedH. U.KumarA. (2011). Meta-analysis of grain yield QTL identified during agricultural drought in grasses showed consensus. 12:319. 10.1186/1471-2164-12-31921679437Szabó-HevérÁ.Lehoczki-KrsjakS.TóthB.PurnhauserL.BuerstmayrH.SteinerB.. (2012). Identification and validation of Fusarium head blight and Fusarium-damaged kernel QTL in a Frontana/Remus DH mapping population. 34, 224–238. 10.1080/07060661.2012.676571Szabó-HevérÁ.SzabolcsL. K.MónikaV.Lászl,óP.JánosP.CsabaL.. (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. 200, 9–26. 10.1007/s10681-014-1229-7VeyrierasJ.-B.GoffinetB.CharcossetA. (2007). MetaQTL: a package of new computational methods for the meta-analysis of QTL mapping experiments. 8:49. 10.1186/1471-2105-8-4917288608WagnerG. P.KinK.LynchV. J. (2013). A model based criterion for gene expression calls using RNA-seq data. 132, 159–164. 10.1007/s12064-013-0178-323615947WaldronB. L.Moreno-SevillaB.AndersonJ. A.StackR. W.FrohbergR. C. (1999). RFLP mapping of QTL for Fusarium head blight resistance in wheat. 39, 805–811. 10.2135/cropsci1999.0011183X003900030032xWangS.WongD.ForrestK.AllenA.ChaoS.HuangB. E.. (2014). Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. 12, 787–796. 10.1111/pbi.1218324646323WinfieldM. O.AllenA. M.BurridgeA. J.BarkerG. L. A.BenbowH. R.WilkinsonP. A.. (2016). High-density SNP genotyping array for hexaploid wheat and its secondary and tertiary gene pool. 14, 1195–1206. 10.1111/pbi.1248526466852XinM. A.DuX.-Y.LiuG.-J.YangZ.-D.HouW.-Q.WangH.-W.. (2015). Cloning and characterization of a novel UDP-glycosyltransferase gene induced by DON from wheat. 14, 830–838. 10.1016/S2095-3119(14)60857-1YangJ.BaiG.ShanerG. E. (2005). Novel quantitative trait loci (QTL) for Fusarium head blight resistance in wheat cultivar Chokwang. 111, 1571–1579. 10.1007/s00122-005-0087-z16172894YangZ.GilbertJ.FedakG.SomersD. J. (2005). Genetic characterization of QTL associated with resistance to Fusarium head blight in a doubled-haploid spring wheat population. 48, 187–196. 10.1139/g04-10415838540YuJ. B.BaiG. H.ZhouW. C.DongY. H.KolbF. L. (2008). Quantitative trait loci for Fusarium head blight resistance in a recombinant inbred population of Wangshuibai/Wheaton. 98, 87–94. 10.1094/PHYTO-98-1-008718943242ZhangG.MergoumM. (2007). Molecular mapping of kernel shattering and its association with Fusarium head blight resistance in a Sumai3 derived population. 115, 757–766. 10.1007/s00122-007-0606-117646963ZhangM.ZhangR.YangJ.LuoP. (2010). Identification of a new QTL for Fusarium head blight resistance in the wheat genotype “Wang shui-bai”. 37, 1031–1035. 10.1007/s11033-009-9809-7ZhangX.BaiG.BockusW.JiX.PanH. (2012a). Quantitative trait loci for Fusarium head blight resistance in US hard winter wheat cultivar Heyne. 52, 1187–1194. 10.2135/cropsci2011.08.0418ZhangX.PanH.BaiG. (2012b). Quantitative trait loci responsible for Fusarium head blight resistance in Chinese landrace Baishanyuehuang. 125, 495–502. 10.1007/s00122-012-1848-022454145ZhangX.ZhouM.RenL.BaiG.MaH.ScholtenO. E.. (2004). Molecular characterization of Fusarium head blight resistance from wheat variety Wangshuibai. 139, 59–64. 10.1007/s10681-004-2298-9ZhouM. P.HaydenM. J.ZhangZ. Y.LuW. Z.MaH. X. (2010). Saturation and mapping of a major Fusarium head blight resistance QTL on chromosome 3BS of Sumai 3 wheat. 51, 19–25. 10.1007/BF0319570620145296ZhouW.KolbF. L.BaiG.ShanerG.DomierL. L. (2002). Genetic analysis of scab resistance QTL in wheat with microsatellite and AFLP markers. 45, 719–727. 10.1139/g02-03412175075ZhouW.KolbF. L.YuJ.BaiG.BozeL. K.DomierL. L. (2004). Molecular characterization of Fusarium head blight resistance in Wangshuibai with simple sequence repeat and amplified fragment length polymorphism markers. 47, 1137–1143. 10.1139/g04-06915644972ZhuZ.BonnettD.EllisM.HeX.HeslotN.DreisigackerS.. (2015). Characterization of Fusarium head blight resistance in a CIMMYT synthetic-derived bread wheat line. 208, 367–375. 10.1007/s10681-015-1612-z
Funding. Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq—Brazil) provided the scholarship—Pós-Doutorado Júnior (Post-Doc Junior) for the first autor. CNPq, CAPES (Coordenação de Aperfeiçoamento de Pessoal de Nível Superior—Brazil) and FAPERGS (Fundação de Amparo á Pesquisa do Estado do Rio Grande do Sul) also supported the other researchers.