- Research article
- Open Access
Characterization of QTL and eQTL controlling early Fusarium graminearum infection and deoxynivalenol levels in a Wuhan 1 x Nyubai doubled haploid wheat population
BMC Plant Biology volume 19, Article number: 536 (2019)
Fusarium head blight (FHB) is a major disease of cereal crops, caused by the fungal pathogen Fusarium graminearum and related species. Breeding wheat for FHB resistance contributes to increase yields and grain quality and to reduce the use of fungicides. The identification of genes and markers for FHB resistance in different wheat genotypes has nevertheless proven challenging.
In this study, early infection by F. graminearum was analyzed in a doubled haploid population derived from the cross of the moderately resistant wheat genotypes Wuhan 1 and Nyubai. Three quantitative trait loci (QTL) were identified: 1AL was associated with lower deoxynivalenol content, and 4BS and 5A were associated with reduced F. graminearum infection at 2 days post inoculation. Early resistance alleles were inherited from Wuhan 1 for QTL 1AL and 4BS and inherited from Nyubai for the 5A QTL. Cis and trans expression QTL (eQTL) were identified using RNA-seq data from infected head samples. Hotspots for trans eQTL were identified in the vicinity of the 1AL and 4BS QTL peaks. Among differentially expressed genes with cis eQTL within the QTL support intervals, nine genes had higher expression associated with FHB early resistance, and four genes had higher expression associated with FHB early susceptibility.
Our analysis of genotype and gene expression data of wheat infected by F. graminearum identified three QTL associated with FHB early resistance, and linked genes with eQTL and differential expression patterns to those QTL. These findings may have applications in breeding wheat for early resistance to FHB.
Wheat is the most important crop in the world in terms of area and Canada is one of the largest producers with 9 million hectares cultivated and 30 million tons of grain produced in 2017 . Allohexaploid bread wheat (Triticum aestivum L.) accounts for over 95% of the global wheat production . Fusarium head blight (FHB) is a major disease of wheat in Canada and other temperate areas of the world, caused predominantly by the fungus Fusarium graminearum Schwabe. Fusarium head blight results in yield losses and contamination of kernels by trichothecene mycotoxins including deoxynivalenol (DON) and derivatives . Strategies to prevent the disease include the cultivation of wheat varieties resistant to FHB, crop rotations and fungicide applications [4, 5]. Achieving improved resistance to FHB is a key goal of current breeding programs in major wheat-producing countries. Several quantitative trait loci (QTL) for FHB resistance and DON reduction have been identified, including Fhb1 on chromosome 3BS, Fhb2 on chromosome 6BS and Fhb5 on chromosome 5AS [6, 7]. Markers for those QTL are used to develop wheat cultivars resistant to FHB in Canada .
A previous study of a doubled haploid (DH) population derived from the cross Wuhan 1 x Nyubai in greenhouse and field trials identified QTL controlling FHB symptoms on chromosomes 2DL, 3BS, and 4B and QTL controlling the accumulation of DON on chromosomes 2DS and 5AS . We have revisited the phenotyping and genotyping of that population by performing experiments at an earlier stage of infection in a controlled environment, analyzing gene expression profiles from infected head samples, and using a genetic map combining single nucleotide polymorphism (SNP) and single sequence repeat (SSR) markers. Our analysis identified a novel QTL on 1AL associated with DON accumulation and confirmed the 4B and 5AS QTL previously identified in the same population . In addition, differentially expressed genes (DEGs) and expression QTL (eQTL) hotspots linked with reduced FHB levels were identified using RNA-seq data.
Phenotypic measures including transcript levels for the fungal genes glyceraldehyde 3-phosphate dehydrogenase (Fg-GAPDH) and β-tubulin (Fg-βTUB) and DON concentration revealed a broad range of infection levels among DH lines derived from Wuhan 1 and Nyubai (Additional file 1). There was a high correlation between the three phenotypic measures (Spearman’s rho coefficients of 0.85, 0.92 and 0.93 for Fg-βTUB vs. DON, Fg-GAPDH vs. DON and Fg-GAPDH vs. Fg-βTUB respectively). The percentages of RNA-seq reads mapping to the F. graminearum genome (%Fg reads) were also highly correlated with the above three phenotypic measures (Spearman’s rho coefficients of 0.88, 0.91 and 0.93 for %Fg reads vs. Fg-βTUB, DON and Fg-GAPDH respectively).
QTL linked with reduced F. graminearum and DON levels
QTL mapping was performed using the level of F. graminearum infection at 2 dpi, as estimated by the four measures described above (%Fg reads, Fg-GAPDH, Fg-βTUB and DON) (Fig. 1). Three QTL were identified (p ≤ 0.05, 1000 permutations): a region on 1AL (peak at 151.66 cM, LOD support interval between 147.57 and 159.79 cM) was associated with DON levels; regions on chromosome 4BS (peak at 30.71 cM, LOD support interval between 15.11 and 34.72 cM) and 5A (peak at 40.50 cM, LOD support interval between 29.76 and 59.90 cM) were associated with F. graminearum levels (%Fg reads).
The analysis of genotypes call frequencies within QTL 1AL, 4BS, and 5A in DH lines ordered according to phenotypes (%Fg reads and DON) (Fig. 2) revealed that DH lines with lower infection levels (dark blue) had genotype AA on QTL 1AL and 4BS, and genotype BB on QTL 5A, and that DH lines with higher infection levels (yellow) had the opposite pattern, which indicates that alleles associated with reduced infection levels were inherited from Wuhan 1 for QTL 1AL and 4BS and inherited from Nyubai for QTL 5A.
In order to map the physical position of genes to QTL intervals, polymorphic marker sequences were aligned against the wheat genome sequence  and a total of 6991 markers with consistent genetic distance and physical locations were identified (Additional file 2). Using these data, the QTL intervals were mapped to 577.89–589.41 Mbp, 7.03–22.50 Mbp and 21.67–461.45 Mbp for QTL 1AL, 4BS and 5A respectively (Fig. 3). Gene and marker densities decrease around the centromere, and regions around the centromere have lower recombination rates , which explains the warping between genetic distances and physical locations. The LOD support interval of 5A overlapped with the centromere, thus covering a large region of the chromosome. In total, the three QTL intervals contained 178 (1AL), 199 (4BS) and 2205 (5A) genes, respectively (Additional file 3).
Identification of eQTL hotspots
Expression QTL mapping was performed using RNA-seq data, and a total of 43,745 eQTL (8573 cis and 35,172 trans) with LOD score above significance thresholds (p ≤ 0.05, 1000 permutations) were retained, corresponding to 35,106 unique wheat genes. Among those, 240 cis eQTL and 17,401 trans eQTL were within the three QTL intervals described above (1AL, 4BS and 5A) (Additional file 4). For cis eQTL peaks, 54 were within the LOD support interval for QTL 1AL, 53 for QTL 4BS and 133 for QTL 5A. For trans eQTL peaks, 10,397 were within the LOD support interval for QTL 1AL, 6295 for QTL 4BS and 709 for QTL 5A. Trans eQTL hotspots were found within the 1AL (peak at 153.6 cM, interval between 144.79 and 161.42) and 4BS (peak at 30.5 cM, interval between 22.73 and 38. 27) QTL intervals (Fig. 4).
Intersecting sets of DEGs within QTL intervals
To identify DEGs between the two parent genotypes that were also linked with observed phenotypic differences between DH lines, three analyses were performed. First, the group of DH lines with > 90% markers corresponding to Wuhan 1 (AA genotype calls) was compared with the group of DH lines with > 90% markers corresponding to Nyubai (BB genotype calls) for each QTL interval separately. A total of 4150 unique DEGs (absolute log2FC > =1 and adjusted p-value<=0.01) were identified in this analysis (3553 for QTL 1AL, 730 for QTL 4BS and 327 for QTL 5A) (Additional file 5). For the second analysis, we compared RNA-seq data (infected heads at 2 dpi) for Wuhan 1, Nyubai and one of the DH lines (HC374) between this experiment and a separate experiment . These data revealed consistent differential gene expression (log2FC) between genotypes in the two experiments (Spearman’s rho coefficients of 0.90, 0.82 and 0.89 for Wuhan 1 vs. Nuybai, Wuhan 1 vs. HC374, and Nyubai vs. HC374 respectively) (Additional file 6). A total of 6452 DEGs were differentially expressed (absolute log2FC > =1 and adjusted p-value<=0.01) between the two parents (Additional file 7). For the third analysis, we compared expression data between two groups (10% of DH lines in each group) with the most extreme phenotypes ranked as described above using %Fg reads and DON levels: lines with the lowest infection levels (HC469, HC385, HC472, HC493, HC454, HC501, HC415, HC498) versus lines with the highest infection levels (HC382, HC401, HC373, HC491, HC391, HC520, HC400, HC405). A total of 28,254 DEGs between the two groups of DH lines were identified (absolute log2FC > =1 and adjusted p-value<=0.01) (Additional file 8). Quantitative reverse transcription-polymerase chain reaction (RT-qPCR) was used to confirm the expression profiles of five DEGs across all samples, showing correlation coefficients between 0.78 and 0.91 for RNA-seq vs. RT-qPCR (Additional file 9).
Within the three QTL intervals, a total of 13 genes were differentially expressed between genotypes, between the two parents, between the two groups of DH lines, and associated with eQTL in cis (eight genes in QTL 1AL, three genes in QTL 4BS and two genes in QTL 5A) (Table 1). Nine genes with higher expression in lines with lower F. graminearum and DON levels were considered to be associated with FHB early resistance. Among those, six genes within the 1AL QTL interval (TraesCS1A01G426000, TraesCS1A01G426500, TraesCS1A01G432900, TraesCS1A01G433000, TraesCS1A01G439000, TraesCS1A01G439100) and two genes within the 4BS QTL interval (TraesCS4B01G022400, TraesCS4B01G024600) had higher expression in genotype AA; while one gene within the 5A QTL interval (TraesCS5A01G114700) had higher expression in genotype BB. On the other hand, the four genes with higher expression in lines with higher F. graminearum and DON levels were considered to be associated with FHB early susceptibility. These included two genes within the 1AL QTL interval (TraesCS1A01G430100, TraesCS1A01G430200) and one gene within the 4BS QTL interval (TraesCS4B01G016900) that had higher expression in genotype BB, and one gene within the 5A QTL interval (TraesCS5A01G196700) that had higher expression in genotype AA.
Numerous QTL for FHB resistance have been identified in wheat and are distributed over all chromosomes [6, 13, 14]. The best characterized QTL for type II resistance (fungal spread within spikes) is Fhb1, located on the short arm of chromosome 3B . Two candidate determinants of the Fhb1 locus have been proposed: a pore-forming toxin-like gene  and a variant of a nuclear histidine-rich calcium-binding protein [17, 18]. In this study, we analyzed RNA-seq and genotype data in two genotypes (Wuhan 1 and Nyubai) with moderate resistance to F. graminearum and 81 DH lines derived from the two parents. Samples were collected at 2 dpi, which has been described as a cleavage in the host-pathogen interaction where the fungus switches from biotrophic to necrotrophic and starts producing higher amounts of DON, while the plant responds by massive transcriptional reprogramming . We have identified three QTL of interest, namely 1AL at 578–589 Mbp (DON QTL and eQTL hotspot), 4BS at 7–23 Mbp (%Fg QTL and eQTL hotspot) and 5A at 22–461 Mbp (%Fg QTL). QTL located in the same regions of chromosomes 4B and 5A have been observed in a previous study with the same population at a later stage of infection . QTL for FHB resistance were also reported in other populations for 1AL (wPt-5577-Xbarc213) , 4BS (Xhbg226-Xgwm149)  and 5A (Xgwm304-Xgwm415) ; recent work separated the latter into a major QTL mapping across the centromere and a minor QTL located on the short arm of chromosome 5A [22, 23]. QTL on 4BS (Fhb4) and 5AS (Fhb5) have both been associated with type I resistance (initial penetration) [24, 25]. As for the vast majority of FHB resistance QTL, genes underlying those QTL remain to be identified and characterized.
The strongest candidate for QTL 1AL is TraesCS1A01G426000, coding for a nucleotide-binding site leucine-rich repeat (NBS-LRR) protein; this gene had a strong differential expression between the two genotypes, higher expression in the most resistant lines, and has a regulatory function that can explain the presence of the trans eQTL peak detected at the same location. Most NBS-LRR proteins are cytoplasmic receptors that recognize specific pathogen molecules, triggering signaling cascades that lead to plant defense responses . Several wheat NBS-LRR genes contributing to disease resistance have been mapped and are being used for improvement in breeding programs . Upregulation of NBS-LRR containing genes during F. graminearum infection has also been reported in wheat genotypes [12, 28]. Other DEGs in the 1AL QTL interval with a cis eQTL close to the 1AL trans eQTL hotspot included an agenet and bromo-adjacent homology (BAH) domain-containing protein (TraesCS1A01G426500), a Na-translocating NADH-quinone reductase subunit (TraesCS1A01G432900) and a RNA-binding protein (TraesCS1A01G433000). Agenet and BAH domains have both been associated with recognition of epigenetic marks on histones, chromatin remodeling and regulation of gene expression . Very few agenet/BAH domain containing proteins have been characterized in plants; however Arabidopsis EML1, an agenet-containing protein, is required for downy mildew race-specific immunity and basal defense . The enzyme NADH-quinone reductase is involved in the metabolism of reactive oxygen species (ROS) and may play a role in plant defense knowing that an oxidative burst at the site of penetration is a common response in plants infected by fungal pathogens, including wheat . The RNA-binding protein contains a hyaluronan-binding domain (HABP4) which indicates interactions with components of the cellular matrix, known to elicit plant defense responses . Two ATP-dependent zinc metalloprotease FtsH located within the 1AL interval had higher expression in susceptible lines and were also differentially expressed between the two parents. These genes play a role in the turnover of photosystem II protein D1 and have also been associated with ROS and plant defense .
A single gene within the 4BS QTL had high differential expression across all comparisons and a cis eQTL peak located close to the 4BS trans eQTL hotspot. This gene (TraesCS4B01G024600) codes for an F-box/LRR-domain containing protein, a class associated with the degradation of proteins and hormone signalling ; this gene is a good candidate for the 4BS QTL and the data presented here points to a regulatory role for this gene in the early wheat response to FHB. Another DEG between genotypes and between most susceptible and resistant lines (TraesCS4B01G022400) codes for a DUF21 domain-containing protein, a domain found in metal transporter proteins; however the cis eQTL peak corresponding to this gene was relatively distant from the 4BS trans eQTL hotspot. Only one gene (TraesCS4B01G016900) in the 4BS QTL had higher expression associated with early susceptibility, coding for a retrovirus-related Pol polyprotein; in this case also, the cis eQTL was distant from the trans eQTL hotspot thus less likely to play a major regulatory function in the wheat response to FHB.
Our analysis did not identify a trans eQTL hotspot associated with the 5A QTL. Two genes in this region were differentially expressed between genotypes, between most resistant and susceptible DH lines, and associated with a cis eQTL. The first one, a NAD-dependent protein deacylase (TraesCS5A01G114700), had higher expression in the most resistant lines; this gene is homologous to Arabidopsis thaliana sirtuin 2 (SIRT2), which has been associated with mitochondrial energy metabolism and negative regulation of plant basal defense responses [35, 36]. The second gene (TraesCS5A01G196700), associated with susceptibility, codes for an ubiquitin. Ubiquitination controls different cell processes in plants including the regulation of plant defense responses . A recent study by Steiner et al.  associated the 5A QTL with type I resistance and anther extrusion. Although included in a list of genes covered by a 5A QTL controlling anther extrusion in wheat , it remains unclear whether DEGs presented here play a role in this phenomenon.
Altogether, results show that combining QTL, eQTL and differential gene expression analysis enhanced the identification of candidate genes controlling early resistance to FHB in wheat. Further work is needed to confirm candidates, and may include increasing the density of markers in regions of interest, additional crosses to obtain lines carrying single QTL to facilitate downstream analyses, and genetic manipulations (e.g. genome editing ) to evaluate the role of each individual gene.
We have identified three significant QTL: 1AL was associated with reduced DON levels, and 4BS and 5A were associated with reduced F. graminearum levels at 2 dpi. The resistance alleles were inherited from Wuhan 1 for the 1AL and 4BS QTL and inherited from Nyubai for the 5A QTL. Trans eQTL hotspots were identified at approximately the same location as QTL on the 1AL and 4BS chromosome arms. Candidate genes corresponding to QTL and eQTL were identified, including a NBS-LRR disease resistance protein in QTL 1AL and an F-box/LRR protein in QTL 4BS.
Eighty-one doubled haploid lines derived from the cross Wuhan 1 x Nyubai [9, 40] were used in this study. Wuhan 1 seeds were obtained from the International Maize and Wheat Improvement Center (CIMMYT), Mexico (accession BW11778) and Nyubai seeds were obtained from the National Small Grain Collection, USA (accession PI 382154). Plants were grown in controlled-environment cabinets with 16 h light at 20 °C and 8 h dark at 16 °C until mid-anthesis. Heads at mid-anthesis were point-inoculated with 10 μL of a F. graminearum (strain DAOM233423, Collection of Fungal Cultures, Agriculture and Agri-Food Canada, Ottawa, Canada) macroconidial spore suspension at 1 × 105 macroconidia/mL using a micropipette between the lemma and palea of two basal florets of each fully developed spikelet, on each treated head. Following inoculation, plants were transferred into a growth chamber at 25 °C where they were misted overhead; pots were disposed in a random order. Misting was for 2 days, 30 s every 1 h, during the light period. Inoculated heads were collected at 2 dpi and 5 to 6 whole heads were pooled into one sample per DH line. An aliquot of ground tissue from each sample was freeze-dried and weighed prior to quantification of the mycotoxin DON. DON analysis was performed using a DON-specific antibody and ELISA analysis . The reported DON concentrations correspond to the average of two technical replicates per sample.
RNA extraction, sequencing and reverse transcription quantitative PCR
Total RNA was extracted and processed for deep paired-end RNA sequencing using Illumina HiSeq 2500 as described in . Raw data were deposited in NCBI Gene Expression Omnibus under accession GSE123548. Synthesis of cDNA and reverse transcription quantitative PCR (RT-qPCR) were performed as described in . Four wheat genes were used to normalize data: glyceraldehyde-3-phosphate dehydrogenase (GAPDH, TraesCS7A01G313100), indole-3 acetaldehyde oxidase (IAAOx, TraesCS2A01G246300), amine oxidase (AOx, TraesCS2A01G327600) and heterogeneous nuclear ribonucleoprotein Q (hn-RNP-Q, TraesCS2A01G390200). For fungal biomass estimation, the expression of two F. graminearum genes was measured: β-tubulin (FGSG_09530) and GAPDH (FGSG_06257). Previous work showed that the expression of those two genes was highly correlated with the amount of F. graminearum DNA in infected plant samples and was a good estimator of fungal biomass . All Primers (Additional file 10), including genome-specific primers for the genes on chromosomes 1A and 5A, were designed and synthesized as described in .
RNA-seq data analysis
The International Wheat Genome Sequencing Consortium (IWGSC) RefSeq v1.0 wheat genome  and Ensembl Fungi (release 35) F. graminearum strain PH-1 genome  were used for RNA-seq data analysis. Wheat and F. graminearum genomes and gene annotations were combined prior to read alignment. RNA-seq reads were preprocessed as described in . Differential gene expression analysis was performed using the DESeq2 R package , using the default negative binomial GLM fitting and Wald statistics. Differentially expressed genes were selected using absolute log2FC > =1 and Benjamini-Hochberg adjusted p-value<=0.01. The same procedure and parameters were used for the analysis of a different dataset  to identify DEGs between the two parents.
Genotyping of 81 DH lines was performed using the Illumina wheat 90 K Infinium iSelect SNP array  and combined with previous genotyping using SSR markers [9, 48] into a revised genetic map (Additional file 11). A total of 12,325 SNPs and SSR markers were retained for the construction of the genetic map. Bins of co-segregating markers were identified using MSTmap . Linkage groups were created using a minimum LOD score of 4 and maximum recombination fraction (RF) of 0.25, and recombination fractions were converted into centimorgan (cM) map distances using the Kosambi mapping function.
Genotype data analysis
For QTL analysis, genotype probabilities and a genome scan by Haley-Knott regression were performed using R/qtl2 , with four phenotypic measures (%Fg reads, DON concentration, Fg-GAPDH and Fg-βTUB). Missing genotype data were imputed using conditional probabilities calculated using a hidden Markov model . LOD support intervals were determined using a drop value of one . LOD scores for eQTL were calculated using 12,325 markers and gene expression data for the 110,790 wheat genes. For this analysis, eQTL were considered cis if corresponding genes were on the same chromosome and within the LOD support interval, and trans otherwise. The distribution of eQTL peaks was evaluated using kernel densities [52, 53] using a chosen bandwidth of 10 cM, and hotspots intervals were defined using peak width at 75% height.
Sequence alignment of markers
SNP and SSR primer sequences were respectively aligned to wheat pseudomolecules using BWA-MEM  and Bowtie2 . A total of 9715 markers with only one best alignment with > = 95% identity and > =95% query length were retained. Of those, a total of 6991 markers with consistent genetic distance and physical locations were identified. Marker genetic positions (cM) were converted to physical positions (Mbp) using a cubic smoothing spline fit.
Availability of data and materials
RNA-seq data generated in this study are available in the NCBI Gene Expression Omnibus under accession GSE123548. Analyzed data are available as Additional files to this article.
- %Fg reads:
Percentage of F. graminearum reads
Differentially expressed gene(s)
days post inoculation
expression quantitative trait loci
Fusarium head blight
Glyceraldehyde 3-phosphate dehydrogenase
Logarithm of the odds
Mega base pair(s)
Nucleotide-binding site leucine-rich repeat
Quantitative trait loci
High-throughput RNA sequencing
Reactive oxygen species
Quantitative reverse transcription-polymerase chain reaction
Single nucleotide polymorphism
Single sequence repeat
FAOSTAT [www.fao.org/faostat]. Accessed 18 Jan 2019.
International Wheat Genome Sequencing Consortium (IWGSC). A chromosome-based draft sequence of the hexaploid bread wheat (Triticum aestivum) genome. Science. 2014;345(6194):1251788.
Foroud NA, Chatterton S, Reid LM, Turkington TK, Tittlemier SA, Gräfenhan T. Fusarium diseases of Canadian grain crops: impact and disease management strategies. In: Future Challenges in Crop Protection Against Fungal Pathogens. New York, NY: Springer; 2014. p. 267–316.
Bai G, Shaner G. Management and resistance in wheat and barley to Fusarium head blight. Annu Rev Phytopathol. 2004;42:135–61.
Wegulo SN, Baenziger PS, Nopsa JH, Bockus WW, Hallen-Adams H. Management of Fusarium head blight of wheat and barley. Crop Prot. 2015;73:100–7.
Buerstmayr H, Ban T, Anderson JA. QTL mapping and marker-assisted selection for Fusarium head blight resistance in wheat: a review. Plant Breed. 2009;128(1):1–26.
Xue S, Xu F, Tang M, Zhou Y, Li G, An X, Lin F, Xu H, Jia H, Zhang L. Precise mapping Fhb5, a major QTL conditioning resistance to Fusarium infection in bread wheat (Triticum aestivum L.). Theor Appl Genet. 2011;123(6):1055–63.
Randhawa HS, Asif M, Pozniak C, Clarke JM, Graf RJ, Fox SL, Humphreys DG, Knox RE, DePauw RM, Singh AK. Application of molecular markers to wheat breeding in Canada. Plant Breed. 2013;132(5):458–71.
Somers DJ, Fedak G, Savard M. Molecular mapping of novel genes controlling Fusarium head blight resistance and deoxynivalenol accumulation in spring wheat. Genome. 2003;46(4):555–64.
Appels R, Eversole K, Feuillet C, Keller B, Rogers J, Stein N, Pozniak CJ, Choulet F, Distelfeld A, Poland J. Shifting the limits in wheat research and breeding using a fully annotated reference genome. Science. 2018;361(6403):eaar7191.
Gill KS. Gene distribution in cereal genomes. In: Cereal Genomics; 2004. p. 361–85.
Pan Y, Liu Z, Rocheleau H, Fauteux F, Wang Y, McCartney C, Ouellet T. Transcriptome dynamics associated with resistance and susceptibility against fusarium head blight in four wheat genotypes. BMC Genomics. 2018;19(1):642.
Venske E, Dos Santos RS, Farias DDR, Rother V, da Maia LC, Pegoraro C, Costa de Oliveira A. Meta-analysis of the QTLome of Fusarium head blight resistance in bread wheat: refining the current puzzle. Front Plant Sci. 2019;10:727.
Steiner B, Buerstmayr M, Michel S, Schweiger W, Lemmens M, Buerstmayr H. Breeding strategies and advances in line selection for Fusarium head blight resistance in wheat. Trop Plant Pathol. 2017;42(3):165–74.
Cuthbert PA, Somers DJ, Thomas J, Cloutier S, Brule-Babel A. Fine mapping Fhb1, a major gene controlling fusarium head blight resistance in bread wheat (Triticum aestivum L.). Theor Appl Genet. 2006;112(8):1465–72.
Rawat N, Pumphrey MO, Liu S, Zhang X, Tiwari VK, Ando K, Trick HN, Bockus WW, Akhunov E, Anderson JA. Wheat Fhb1 encodes a chimeric lectin with agglutinin domains and a pore-forming toxin-like domain conferring resistance to Fusarium head blight. Nat Genet. 2016;48(12):1576–80.
Su Z, Bernardo A, Tian B, Chen H, Wang S, Ma H, Cai S, Liu D, Zhang D, Li T, et al. A deletion mutation in TaHRC confers Fhb1 resistance to Fusarium head blight in wheat. Nat Genet. 2019;51(7):1099–105.
Li G, Zhou J, Jia H, Gao Z, Fan M, Luo Y, Zhao P, Xue S, Li N, Yuan Y, et al. Mutation of a histidine-rich calcium-binding-protein gene in wheat confers resistance to Fusarium head blight. Nat Genet. 2019;51(7):1106–12.
Samad-Zamini M, Schweiger W, Nussbaumer T, Mayer KFX, Buerstmayr H. Time-course expression QTL-atlas of the global transcriptional response of wheat to Fusarium graminearum. Plant Biotechnol J. 2017;15(11):1453–64.
Semagn K, Skinnes H, Bjørnstad Å, Marøy AG, Tarkegne Y. Quantitative trait loci controlling Fusarium head blight resistance and low deoxynivalenol content in hexaploid wheat population from ‘Arina’and NK93604. Crop Sci. 2007;47(1):294–303.
Xue S, Li G, Jia H, Xu F, Lin F, Tang M, Wang Y, An X, Xu H, Zhang L, et al. Fine mapping Fhb4, a major QTL conditioning resistance to Fusarium infection in bread wheat (Triticum aestivum L.). Theor Appl Genet. 2010;121(1):147–56.
Buerstmayr M, Steiner B, Wagner C, Schwarz P, Brugger K, Barabaschi D, Volante A, Vale G, Cattivelli L, Buerstmayr H. High-resolution mapping of the pericentromeric region on wheat chromosome arm 5AS harbouring the Fusarium head blight resistance QTL Qfhs.Ifa-5A. Plant Biotechnol J. 2018;16(5):1046–56.
Steiner B, Buerstmayr M, Wagner C, Danler A, Eshonkulov B, Ehn M, Buerstmayr H. Fine-mapping of the Fusarium head blight resistance QTL Qfhs.Ifa-5A identifies two resistance QTL associated with anther extrusion. Theor Appl Genet. 2019;132(7):2039–53.
Buerstmayr H, Steiner B, Hartl L, Griesser M, Angerer N, Lengauer D, Miedaner T, Schneider B, Lemmens M. Molecular mapping of QTLs for Fusarium head blight resistance in spring wheat. II. Resistance to fungal penetration and spread. Theor Appl Genet. 2003;107(3):503–8.
Lin F, Xue SL, Zhang ZZ, Zhang CQ, Kong ZX, Yao GQ, Tian DG, Zhu HL, Li CJ, Cao Y, et al. Mapping QTL associated with resistance to Fusarium head blight in the Nanda2419 x Wangshuibai population. II: type I resistance. Theor Appl Genet. 2006;112(3):528–35.
Eitas TK, Dangl JL. NB-LRR proteins: pairs, pieces, perception, partners, and pathways. Curr Opin Plant Biol. 2010;13(4):472–7.
Marone D, Russo M, Laidò G, De Leonardis A, Mastrangelo A. Plant nucleotide binding site–leucine-rich repeat (NBS-LRR) genes: active guardians in host defense responses. Int J Mol Sci. 2013;14(4):7302–26.
Kugler KG, Siegwart G, Nussbaumer T, Ametz C, Spannagl M, Steiner B, Lemmens M, Mayer KF, Buerstmayr H, Schweiger W: Quantitative trait loci-dependent analysis of a gene co-expression network associated with Fusarium head blight resistance in bread wheat (Triticum aestivum L.). BMC genomics 2013;14(1):728.
Maurer-Stroh S, Dickens NJ, Hughes-Davies L, Kouzarides T, Eisenhaber F, Ponting CP. The Tudor domain ‘Royal Family’: tudor, plant agenet, chromo, PWWP and MBT domains. Trends Biochem Sci. 2003;28(2):69–74.
Tsuchiya T, Eulgem T. EMSY-like genes are required for full RPP7-mediated race-specific immunity and basal defense in Arabidopsis. Mol Plant-Microbe Interact. 2011;24(12):1573–81.
Greenshields DL, Liu G, Selvaraj G, Wei Y. Differential regulation of wheat quinone reductases in response to powdery mildew infection. Planta. 2005;222(5):867–75.
Galletti R, De Lorenzo G, Ferrari S. Host-derived signals activate plant innate immunity. Plant Signal Behav. 2009;4(1):33–4.
Gupta SK, Sharma M, Deeba F, Pandey V. Role of reactive oxygen species in photophosphorylation and damage to D1 protein: past and present. In: Reactive Oxygen Species Plants; 2018. p. 165.
Hellmann H, Estelle M. Plant development: regulation by protein degradation. Science. 2002;297(5582):793–7.
Wang C, Gao F, Wu J, Dai J, Wei C, Li Y. Arabidopsis putative deacetylase AtSRT2 regulates basal defense by suppressing PAD4, EDS5 and SID2 expression. Plant Cell Physiol. 2010;51(8):1291–9.
Konig AC, Hartl M, Pham PA, Laxa M, Boersema PJ, Orwat A, Kalitventseva I, Plochinger M, Braun HP, Leister D, et al. The Arabidopsis class II sirtuin is a lysine deacetylase and interacts with mitochondrial energy metabolism. Plant Physiol. 2014;164(3):1401–14.
Devoto A, Muskett PR, Shirasu K. Role of ubiquitination in the regulation of plant defence against pathogens. Curr Opin Plant Biol. 2003;6(4):307–11.
Muqaddasi QH, Reif JC, Röder MS, Basnet BR, Dreisigacker S. Genetic mapping reveals large-effect QTL for anther extrusion in CIMMYT spring wheat. Agronomy. 2019;9(7):407.
Kumar R, Kaur A, Pandey A, Mamrutha HM, Singh GP. CRISPR-based genome editing in wheat: a comprehensive review and future prospects. Mol Biol Rep. 2019;46(3):3557–69.
Somers DJ, Thomas J, Depauw R, Fox S, Humphreys G, Fedak G. Assembling complex genotypes to resist Fusarium in wheat (Triticum aestivum L.). Theor Appl Genet. 2005;111(8):1623–31.
Savard ME, Sinha RC, Lloyd Seaman W, Fedak G. Sequential distribution of the mycotoxin deoxynivalenol in wheat spikes after inoculation with Fusarium graminearum. Can J Plant Pathol. 2000;22(3):280–5.
Hu X, Rocheleau H, McCartney C, Biselli C, Bagnaresi P, Balcerzak M, Fedak G, Yan Z, Vale G, Khanizadeh S, et al. Identification and mapping of expressed genes associated with the 2DL QTL for fusarium head blight resistance in the wheat line Wuhan 1. BMC Genet. 2019;20(1):47.
Harris LJ, Balcerzak M, Johnston A, Schneiderman D, Ouellet T. Host-preferential Fusarium graminearum gene expression during infection of wheat, barley, and maize. Fungal Biol. 2016;120(1):111–23.
RefSeq v1.0 wheat genome assembly and annotation [https://urgi.versailles.inra.fr/download/iwgsc]. Accessed 30 May 2017.
Fusarium graminearum str. PH-1 genome assembly and annotation [ftp://ftp.ensemblgenomes.org/pub/fungi/release-35]. Accessed 30 May 2017.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15(12):550.
Wang S, Wong D, Forrest K, Allen A, Chao S, Huang BE, Maccaferri M, Salvi S, Milner SG, Cattivelli L. Characterization of polyploid wheat genomic diversity using a high-density 90 000 single nucleotide polymorphism array. Plant Biotechnol J. 2014;12(6):787–96.
Somers DJ, Isaac P, Edwards K. A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theor Appl Genet. 2004;109(6):1105–14.
Wu Y, Bhat PR, Close TJ, Lonardi S. Efficient and accurate construction of genetic linkage maps from the minimum spanning tree of a graph. PLoS Genet. 2008;4(10):e1000212.
Broman KW, Wu H, Sen Ś, Churchill GA. R/QTL: QTL mapping in experimental crosses. Bioinformatics. 2003;19(7):889–90.
Manichaikul A, Dupuis J, Sen S, Broman KW. Poor performance of bootstrap confidence intervals for the location of a quantitative trait locus. Genetics. 2006;174(1):481–9.
Scott DW. Multivariate density estimation: theory, practice, and visualization. Hoboken, NJ: Wiley; 2015.
Silverman BW. Density estimation for statistics and data analysis. London: Chapman and Hall; 1986.
Li H: Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. arXiv preprint arXiv:13033997 2013.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9(4):357–9.
We would like to thank the DNA Technologies Unit at National Research Council Canada, Saskatoon for performing the RNA sequencing, Dr. Barbara Blackwell and Sally Buffam for DON analyses, and the IWGSC for pre-publication access to the RefSeq v1.0 wheat genome assembly and annotation.
This project was supported by Canadian Growing Forward 2 Program (project J-000412) and Genomics Research and Development Initiative (projects J-000008 and J-001580) awarded to TO, National Wheat Improvement Program (project: Fusarium screening, genetics and marker development) awarded to CM by Western Grains Research Foundation and Agriculture and Agri-Food Canada, and National Research Council Canada’s Canadian Wheat Improvement project A1–011652 awarded to YP. The funding bodies did not play a role in the design of the study, in collection, analysis, or interpretation of data or in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Estimation of the fungal biomass in two parents and in 81 DH lines at 2 dpi. A) Percentage of F. graminearum reads in RNA-seq data. B) The mycotoxin DON measured by ELISA. C) F. graminearum GAPDH and D) β-tubulin RNA levels measured using RT-qPCR.
Genetic distance versus physical location of markers best hits on wheat chromosomes. Markers with consistent behavior between genetic distance and physical location (6991 out of 9715 hits) are colored in black. The blue line corresponds to a cubic smoothing spline fit between the genetic distance and the physical position of the 6991 markers.
Wheat genes within QTL intervals. gene_id, IWGSC RefSeq v1.0 gene ID; chr, chromosome; start, gene start position; end, gene end position; strand, gene strand; description, IWGSC RefSeq v1.0 gene annotation.
Cis and trans eQTL with LOD score above significance thresholds. Type, cis or trans; chr, peak chromosome; pos, peak position in cM; lod, LOD score; ci_lo, one-LOD interval lower limit; ci_hi, one-LOD interval higher limit; gene_id, IWGSC RefSeq v1.0 gene ID; description, IWGSC RefSeq v1.0 gene annotation.
Differentially expressed genes between genotypes within the limits of the three QTL intervals. gene_id, IWGSC RefSeq v1.0 gene ID; chr, chromosome; AA, mean of normalized counts for DH lines with genotype AA in the given QTL region; BB, mean of normalized counts for lines with genotype BB in the given QTL region; log2FC, log2 fold change; padj, adjusted p-value; description, IWGSC RefSeq v1.0 gene annotation.
Comparison between differential expression in Wuhan 1, Nyubai and HC374 at 2 days post−inoculation in two experiments. Experiment 1, Pan et al. 2018; Experiment 2, this paper. A) Wuhan 1 vs. Nyubai; B) Wuhan 1 vs. HC374; C) Nyubai vs. HC374. The red line corresponds to a linear regression between the log2 fold change in the two experiments.
Differentially expressed genes between Wuhan 1 and Nyubai at 2 dpi. gene_id, IWGSC RefSeq v1.0 gene ID; chr, chromosome; Wuhan 1, mean of normalized counts for Wuhan 1; Nyubai, mean of normalized counts for Nyubai; log2FC, log2 fold change; padj, adjusted p-value; description, IWGSC RefSeq v1.0 gene annotation.
Differentially expressed genes between DH lines with extreme phenotypes. gene_id, IWGSC RefSeq v1.0 gene ID; chr, chromosome; LF, mean of normalized counts for DH lines with lower levels of F. graminearum and DON; HF, mean of normalized counts for DH lines with higher levels of F. graminearum and DON; log2FC, log2 fold change; padj, adjusted p-value; description, IWGSC RefSeq v1.0 gene annotation.
Relative expression levels as quantified by RT-qPCR for five selected wheat genes, and correlation with RNA-seq data.
Primers used for RT-qPCR analyses.
Genetic map (cM) and genotype calls (A: Wuhan 1; B: Nyubai) for the DH population using 12,325 polymorphic markers (SSR markers and SNP from wheat 90 K Infinium beadchip array). chr, chromosome; cM, centimorgan.
About this article
Cite this article
Fauteux, F., Wang, Y., Rocheleau, H. et al. Characterization of QTL and eQTL controlling early Fusarium graminearum infection and deoxynivalenol levels in a Wuhan 1 x Nyubai doubled haploid wheat population. BMC Plant Biol 19, 536 (2019). https://doi.org/10.1186/s12870-019-2149-4
- Fusarium head blight (FHB)
- Quantitative trait loci (QTL)
- Expression quantitative trait loci (eQTL)
- Doubled haploid (DH)
- Fusarium graminearum
- Triticum aestivum L.