Fine mapping of a large-effect QTL conferring Fusarium crown rot resistance on the long arm of chromosome 3B in hexaploid wheat
- Zhi Zheng†1, 2, 3,
- Jian Ma†1, 4,
- Jiri Stiller1,
- Qiang Zhao5,
- Qi Feng5,
- Frédéric Choulet6,
- Catherine Feuillet6,
- You-Liang Zheng4,
- Yuming Wei4,
- Bin Han5,
- Guijun Yan3,
- John M. Manners1 and
- Chunji Liu1, 2Email author
© Zheng et al. 2015
Received: 16 June 2015
Accepted: 16 October 2015
Published: 23 October 2015
Fusarium crown rot (FCR) is a major cereal disease in semi-arid areas worldwide. Of the various QTL reported, the one on chromosome arm 3BL (Qcrs.cpi-3B) has the largest effect that can be consistently detected in different genetic backgrounds. Nine sets of near isogenic lines (NILs) for this locus were made available in a previous study. To identify markers that could be reliably used in tagging the Qcrs.cpi-3B locus, a NIL-derived population consisting of 774 F10 lines were generated and exploited to assess markers selected from the existing linkage map and generated from sequences of the 3B pseudomolecule.
This is the first report on fine mapping a QTL conferring FCR resistance in wheat. By three rounds of linkage mapping using the NILs and the NIL-derived population, the Qcrs.cpi-3B locus was mapped to an interval of 0.7 cM covering a physical distance of about 1.5 Mb. Seven markers co-segregating with the locus were developed. This interval contains a total of 63 gene-coding sequences based on the 3B pseudomolecule, and six of them were known to encode disease resistance proteins. Several of the genes in this interval were among those responsive to FCR infection detected in an earlier study.
The accurate localization of the Qcrs.cpi-3B locus and the development of the markers co-segregating with it should facilitate the incorporation of this large-effect QTL conferring FCR resistance into breeding programs as well as the cloning of the gene(s) underlying the QTL.
KeywordsFusarium crown rot Fine mapping Hexaploid wheat Co-segregating SSR marker
Fusarium crown rot (FCR) is a chronic and serious disease of cereals. Field surveys showed that F. pseudograminearum is the most prevalent pathogen for FCR in Queensland and New South Wales in Australia but many different species of Fusarium can cause this disease . Due most likely to the high intensity of cereal in cropping system combined with wider adoption of minimum tillage for moisture conservation, FCR has become more prevalent in many parts of the semiarid regions in recent years . A survey in 2009 found that FCR causes an estimated annual yield loss of $80 million Australia dollars in the wheat industry alone . A study in the Pacific Northwest in USA showed that FCR could significantly reduce yield of both wheat and barley . Additionally, FCR infected plants could also produce mycotoxins in grains as well as other tissues. The presence of these compounds in food and feeds can be harmful for human and livestock .
Growing resistant varieties has long been recognised as an integral part in effectively managing FCR . Working toward the breeding of varieties resistant to this disease, sources of resistance were identified by germplasm screenings [6, 7]. Over the last decade, significant effort has also been made in identifying QTL conferring FCR resistance. Several QTL have been detected from several different sources of resistance . Of the QTL reported so far, the one on 3BL consistently gave the largest effects and was identified from several genotypes [9–11]. Limited data indicates that the 3BL locus also confers field resistance to FCR and it reduces whitehead incidence significantly under field conditions .
QTL mapping has become a routine procedure in locating genes controlling quantitative traits to specific genomic regions. However, QTL mapping has only limited resolution and markers obtained from such studies may not be reliably used in tagging the targeted loci . The main reason why QTL mapping provides only limited resolution is due to the heterogeneity in genetic backgrounds in QTL mapping populations in regard to a targeted locus. Some non-targeted traits segregating in the QTL mapping populations could interfere with the accurate phenotyping of a targeted trait thus lines containing the same allele may exhibit different levels of resistance. Recent studies showed that both plant height and heading date have significant effects on FCR severity in both wheat and barley. LOD values and magnitudes of QTL conferring FCR resistance were both reduced when the effect of height and heading date were accounted for by the covariance analysis [13–15]. Thus these traits need to be fixed in segregating populations which are suitable for developing markers tightly linked to a FCR locus.
Fixing genetic backgrounds in regarding to a targeted locus can be achieved with near-isogenic lines (NILs) [16–19] or populations generated from NILs . Taken advantage of the newly released chromosome 3B pseudomolecule , we used both NILs  and a NIL-derived population in developing tightly linked markers to the large-effect FCR locus on 3BL and results obtained are reported in this paper.
Primers for new markers successfully mapped near the Qcrs.cpi-3B locus#
Forward primer/position (bp)
Reverse primer/position (bp)
Product size (bp)
To further define the FCR locus, the 3B pseudomolecule was used again to develop new markers for the targeted interval. Sixty markers were designed in this third round of marker development with an average density of about 0.2 Mbp between the markers (Additional file 2: Table S1). Seven of the 60 markers (including CS3BLCR-06, CS3BLCR-07, CS3BLCR-08. CS3BLCR-09, CS3BLCR-10, CS3BLCR-11 and CS3BLCR-12) detected polymorphism in the mapping population and they all co-segregated with the FCR locus (Fig. 2).
Unique sequences in the targeted interval harbouring the Qcrs.cpi-3B locus
Disease resistance-responsive family protein
Disease resistance rpp13-like protein 3-like
Disease resistance protein rga2
Disease resistance rpp13-like protein 3-like
Disease resistance rpp13-like protein 2-like
NBS-LRR resistance partial
Disease resistance protein rpm1
Disease resistance protein rga4
Disease resistance protein rga4
AP3-complex subunit beta-A-like
Gibberellin 2-beta-dioxygenase 8-like
Discussion and conclusions
FCR is a serious threat to cereal production in semi-arid areas worldwide and growing resistant varieties is an essential component in effectively managing the disease. Working toward the breeding of resistant varieties, sources of resistance have been identified and several QTL controlling FCR have been reported. In the study reported here, we attempted to develop markers that can be reliably used to tag the FCR locus on 3BL. Using markers developed from the chromosome 3B pseudomolecule, the large-effect QTL was successfully mapped to an interval of 0.7 cM or 1.5 Mbp based on the nine sets of NILs and a NIL-derived population consisting of 774 lines. Within this interval, seven markers co-segregating with the targeted locus were developed and 15 unique genes or contigs related to disease resistance were identified. These results should facilitate not only the incorporation of this locus into breeding programs but also the effort of cloning gene(s) underlying the FCR locus.
The 18 markers mapped in this study covers a distance of 3.5 cM (Fig. 1). These markers covered a distance of 7.1 cM in the linkage map reported previously . It is known that recombination rates at meiosis can be affected by many factors. In addition to the relateness of the parental genotypes , recombination frequencies can be regulated by major genes  and may also differ between egg and pollen mother cells [27, 28]. It is also well known that the presence of structural rearrangements between the parental genotypes can significantly reduce the recombination rates . As the NILs and the NIL-derived population used for fine mapping in this study were all obtained from crosses between bread wheat (Triticum aestivum L.) and T. spelta genotypes [9, 22], the relatively distant genetic relatedness between the parental genotypes likely contributed to the reduced linkage distances obtained in this study. A population significantly larger than the one used in this study is being produced and will be used to further define the targeted interval and identify the gene(s) conferring FCR resistance at this locus.
The availability of the chromosome 3B pseudomolecule  significantly facilitated the study reported here and there is no reason why the same approach can not be used to fine map genes of other traits on this chromosome. However, of the 98 pairs of primers designed in this study only 13 detected polymorphism and were successfully mapped in the targeted interval. This rate of success does not seem to be very high but could be even lower if populations generated from more closely related genotypes were used. The 3B pseudomolecule makes it feasible now to exploit a high-throughput genotyping approach in identifying genes controlling any trait on this chromosome. The one based on low-coverage whole genome resequencing  is very appealing especially in considering that chromosome 3B is readily sortable by flow cytometry .
Nine sets of NILs for the 3BL FCR locus derived from three segregating populations (including 120 F4 lines from ‘Aus13832’/’CSCS6’, 108 BC1F4 lines from ‘Janz’*2/’CSCR6’ and 125 F7 lines from ‘Lang’/’CSCR6’) were available  and they were used in identifying markers flanking the locus in this study. The resistance donor ‘CSCR6’ belongs to the taxon T. spelta. It is one of the genotypes which showed the best resistance among the over 2200 genotypes screened . The three susceptible parents are all Australian varieties. Detailed mapping of the 3BL locus was conducted using a NIL-derived population consisting of 774 F10 lines. This fine mapping population was developed as part of this study from a single plant of Janz/CSCR6  based on the heterogeneous inbred family (HIF) method . An SSR marker, gwm181 locating near the peak of the FCR QTL on 3BL , was used to identify heterozygous lines in developing NILs. A single hetetozygous plant was used in developing the fine mapping population (at F10 generation) using the fast-generation technique .
Genotyping and linkage analyses
Leaf tissues from the NILs and the fine mapping population were collected and stored at −80 °C until processing. Genomic DNA was extracted using the CTAB protocol . PCR reactions for the marker analyses were carried out using α [33P] dCTP (3000 ci/mmol) following the manufacturer’s protocol (Multiplex-Ready Marker User Handbook, version 2.0). The amplified products were mixed with an equal volume of loading dye, denatured at 95 °C for 10 min, and 3.8 μl amplified samples were separated on a 5 % polyacrylamide gel containing 8 M urea at 100 W for 2 h. The gels were subsequently dried using a gel dryer for 50 min at 80 °C and exposed to Kodak X-Omat X-ray film for 4–6 days.
To construct a linkage map for the 3BL region, SSR markers from the existing linkage maps [23, 35] and new PCR-based markers developed based on the sequence of the 3B pseudomolecule  were screened for polymorphism against the nine sets of NILs . Gene-coding sequences (CDSs) and 3B pseudomolecule sequence were downloaded from https://urgi.versailles.inra.fr/gb2/gbrowse/wheat_annot_3B/. Primers for all of the new markers were designed with SSRPrimerII . Linkage analysis was carried out using the computer package JoinMap 4.0 .
Evaluation of resistance to FCR
A highly aggressive F. pseudograminearum isolate, CS3096, was used in FCR assessment. This isolate was collected in northern New South Wales, Australia and maintained in CSIRO collection . The methods used for inoculums preparation, inoculation and FCR assessment were based on that described by Li et al. . Briefly, inoculums was prepared using plates of ½ strength potato dextrose agar (PDA). The inoculated plates were incubated for 7 days at room temperature before the mycelium was scraped. The plates were then incubated for a further 5–7 days under a combination of cool white and black (UVA) fluorescent lights with a 12 h photoperiod. The spores were harvested and the concentration of spore suspension was adjusted to 1 × 106 spores per millilitre in distilled water. Tween 20 was added (0.1 % v/v) to the spore suspension prior to use for inoculation.
Seeds were germinated in Petri dishes on two layers of filter paper saturated with water. The germinated seedlings were immersed in the spore suspension for 1 min and two seedlings were planted into square punnets of a 56-well tray (Rite Grow Kwik Pots, Garden City Plastics, Australia) containing stem sterilized University of California mix C (50 % sand and 50 % peat v/v). The punnets were arranged in a randomized block design in a controlled environment facility (CEF). Settings for the CEF were: 25/16 (±1) °C day/night temperature and 65/85(±5)% day/night relative humidity, and a 14-hour photoperiod with 500 μmol m−2s−1 photon flux density at the level of the plant canopy. To promote FCR development, water-stress was applied during the FCR assessment. Inoculated seedlings were watered only when wilt symptoms appeared. Three trials were conducted using the fine mapping population and five additional trials were then carried out to assess those recombinant lines identified based on markers flanking the targeted interval containing the FCR locus. Each trial contains two replicates, each with 14 seedlings. FCR severity was assessed 35 days post inoculation, using a 0 (no obvious symptom) to 5 (whole plant severely to completely nectrotic) scale as described by Li et al. .
Identification of genes of interest in the target region
Based on the fine mapping results described above, CDSs located in the targeted genomic region were retrieved from the 3B pseudomolecule . Paired RNA reads from the first of the nine sets of NILs obtained from an earlier study  were re-analysed for identifying transcripts of interest in the targeted genomic region. RNA datasets were trimmed by SolexaQA scripts (http://solexaqa.sourceforge.net/) to a minimum quality value of 30 and a minimum length of 70. Alignments of RNA reads from resistant and sensitive NILs were performed using the Biokanga suite (a tool developed at CSIRO, not published) with 2 mismatches allowed per read. Paired reads for all reads aligned to the defined region containing the FCR locus were retrieved by in house built Perl scripts and de novo assembled using CLC Genomics Workbench 7.
The authors state that the experiments conducted in this study comply with the current laws in Australia where they were conducted.
Availability of data and materials
All the supporting data are included as supplementary files.
This publication is based upon work supported by a joint CAS/CSIRO project (Project No R-1910-1) and an project from the International Science & Cooperation Program of China (No. 2015DFA30600). ZZ is grateful to the University of Western Australia and China Scholarship Council (CSC) for his PhD scholarships, and JM thanks the Sichuan Agricultural University and the China Scholarship Council for funding his visit to CSIRO Agriculture.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Akinsanmi OA, Mitter V, Simpfendorfer S, Backhouse D, Chakraborty S. Identity and pathogenicity of Fusarium spp. isolated from wheat fields in Queensland and northern New South Wales. Crop Pasture Sci. 2004;55(1):97–107.View ArticleGoogle Scholar
- Chakraborty S, Liu CJ, Mitter V, Scott JB, Akinsanmi OA, Ali S, et al. Pathogen population structure and epidemiology are keys to wheat crown rot and Fusarium head blight management. Australas Plant Pathol. 2006;35(6):643–55.View ArticleGoogle Scholar
- Murray GM, Brennan JP. Estimating disease losses to the Australian wheat industry. Australas Plant Pathol. 2009;38(6):558–70.View ArticleGoogle Scholar
- Smiley RW, Gourlie JA, Easley SA, Patterson L-M, Whittaker RG. Crop damage estimates for crown rot of wheat and barley in the Pacific Northwest. Plant Dis. 2005;89(6):595–604.View ArticleGoogle Scholar
- Mudge AM, Dill-Macky R, Dong Y, Gardiner DM, White RG, Manners JM. A role for the mycotoxin deoxynivalenol in stem colonisation during crown rot disease of wheat caused by Fusarium graminearum and Fusarium pseudograminearum. Physiol Mol Plant Pathol. 2006;69(1):73–85.View ArticleGoogle Scholar
- Purss G. Studies of varietal resistance to crown rot of wheat caused by Fusarium graminearum Schw. Queensland J Agric Anim Sci. 1966;23:475–98.Google Scholar
- Liu C, Ma J, Li H, Liu Y, Liu G, Wen S, et al. The homoeologous regions on long arms of group 3 chromosomes in wheat and barley harbour major crown rot resistance loci. Czech J Genet Plant Breed. 2011;47(SI):S109–14.Google Scholar
- Liu C, Ogbonnaya FC. Resistance to Fusarium crown rot in wheat and barley: a review. Plant Breed. 2015;134(4):365–72.View ArticleGoogle Scholar
- Ma J, Li HB, Zhang CY, Yang XM, Liu YX, Yan GJ, et al. Identification and validation of a major QTL conferring crown rot resistance in hexaploid wheat. Theor Appl Genet. 2010;120(6):1119–28.View ArticlePubMedGoogle Scholar
- Li HB, Xie GQ, Ma J, Liu GR, Wen SM, Ban T, et al. Genetic relationships between resistances to Fusarium head blight and crown rot in bread wheat (Triticum aestivum L.). Theor Appl Genet. 2010;121(5):941–50.View ArticlePubMedGoogle Scholar
- Poole GJ, Smiley RW, Paulitz TC, Walker CA, Carter AH, See DR, et al. Identification of quantitative trait loci (QTL) for resistance to Fusarium crown rot (Fusarium pseudograminearum) in multiple assay environments in the Pacific Northwestern US. Theor Appl Genet. 2012;125(1):91–107.PubMed CentralView ArticlePubMedGoogle Scholar
- Paterson AH, Lander ES, Hewitt JD, Peterson S, Lincoln SE, Tanksley SD. Resolution of quantitative traits into Mendelian factors by using a complete linkage map of restriction fragment length polymorphisms. Nature. 1988;335(6192):721–6.View ArticlePubMedGoogle Scholar
- Liu YX, Yang XM, Ma J, Wei YM, Zheng YL, Ma HX, et al. Plant height affects Fusarium crown rot severity in wheat. Phytopathology. 2010;100(12):1276–81.View ArticlePubMedGoogle Scholar
- Zheng Z, Kilian A, Yan G, Liu C. QTL conferring fusarium crown rot resistance in the elite bread wheat variety EGA Wylie. PLoS One. 2014;9(4):e96011.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen G, Liu Y, Ma J, Zheng Z, Wei Y, McIntyre CL, et al. A Novel and major quantitative trait locus for fusarium crown rot resistance in a genotype of Wild Barley (Hordeum spontaneum L.). PLoS One. 2013;8(3):e58040.PubMed CentralView ArticlePubMedGoogle Scholar
- Young ND, Zamir D, Ganal MW, Tanksley SD. Use of isogenic lines and simultaneous probing to identify DNA markers tightly linked to the tm-2a gene in tomato. Genetics. 1988;120(2):579–85.PubMed CentralPubMedGoogle Scholar
- Paterson AH, DeVerna JW, Lanini B, Tanksley SD. Fine mapping of quantitative trait loci using selected overlapping recombinant chromosomes, in an interspecies cross of tomato. Genetics. 1990;124(3):735–42.PubMed CentralPubMedGoogle Scholar
- Pumphrey MO, Bernardo R, Anderson JA. Validating the QTL for Fusarium head blight resistance in near-isogenic wheat lines developed from breeding populations. Crop Sci. 2007;47(1):200–6.View ArticleGoogle Scholar
- Kaeppler S, Phillips R, Kim T. Use of near-isogenic lines derived by backcrossing or selfing to map qualitative traits. Theor Appl Genet. 1993;87(1–2):233–7.PubMedGoogle Scholar
- Liu S, Zhang X, Pumphrey MO, Stack RW, Gill BS, Anderson JA. Complex microcolinearity among wheat, rice, and barley revealed by fine mapping of the genomic region harboring a major QTL for resistance to Fusarium head blight in wheat. Funct Integr Genomics. 2006;6(2):83–9.View ArticlePubMedGoogle Scholar
- Choulet F, Alberti A, Theil S, Glover N, Barbe V, Daron J, et al. Structural and functional partitioning of bread wheat chromosome 3B. Science. 2014;345(6194):1249721.View ArticlePubMedGoogle Scholar
- Ma J, Yan GJ, Liu CJ. Development of near-isogenic lines for a major QTL on 3BL conferring Fusarium crown rot resistance in hexaploid wheat. Euphytica. 2012;183(2):147–52.View ArticleGoogle Scholar
- Paux E, Sourdille P, Salse J, Saintenac C, Choulet F, Leroy P, et al. A physical map of the 1-gigabase bread wheat chromosome 3B. Science. 2008;322(5898):101–4.View ArticlePubMedGoogle Scholar
- Ma J, Stiller J, Zhao Q, Feng Q, Cavanagh C, Wang P, et al. Transcriptome and allele specificity associated with a 3bl locus for fusarium crown rot resistance in bread wheat. PLoS One. 2014;9(11):e113309.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu C, Devos K, Witcombe J, Pittaway T, Gale M. The effect of genome and sex on recombination rates in Pennisetum species. Theor Appl Genet. 1996;93(5–6):902–8.View ArticlePubMedGoogle Scholar
- Temin RG, Marthas M. Factors influencing the effect of segregation distortion in natural populations of Drosophila melanogaster. Genetics. 1984;107(3):375–93.PubMed CentralPubMedGoogle Scholar
- De Vicente M, Tanksley S. Genome-wide reduction in recombination of backcross progeny derived from male versus female gametes in an interspecific cross of tomato. Theor Appl Genet. 1991;83(2):173–8.View ArticlePubMedGoogle Scholar
- Busso C, Liu C, Hash C, Witcombe J, Devos K, De Wet J, et al. Analysis of recombination rate in female and male gametogenesis in pearl millet (Pennisetum glaucum) using RFLP markers. Theor Appl Genet. 1995;90(2):242–6.View ArticlePubMedGoogle Scholar
- Devos KM, Atkinson M, Chinoy C, Francis H, Harcourt R, Koebner R, et al. Chromosomal rearrangements in the rye genome relative to that of wheat. Theor Appl Genet. 1993;85(6–7):673–80.View ArticlePubMedGoogle Scholar
- Huang X, Zhao Y, Wei X, Li C, Wang A, Zhao Q, et al. Genome-wide association study of flowering time and grain yield traits in a worldwide collection of rice germplasm. Nat Genet. 2012;44(1):32–9.View ArticleGoogle Scholar
- Doležel J, Kubaláková M, Paux E, Bartoš J, Feuillet C. Chromosome-based genomics in the cereals. Chromosome Res. 2007;15(1):51–66.View ArticlePubMedGoogle Scholar
- Tuinstra M, Ejeta G, Goldsbrough P. Heterogeneous inbred family (HIF) analysis: a method for developing near-isogenic lines that differ at quantitative trait loci. Theor Appl Genet. 1997;95(5–6):1005–11.View ArticleGoogle Scholar
- Zheng Z, Wang H, Chen G, Yan G, Liu C. A procedure allowing up to eight generations of wheat and nine generations of barley per annum. Euphytica. 2013;191(2):311–6.View ArticleGoogle Scholar
- Anderson JA, Ogihara Y, Sorrells ME, Tanksley SD. Development of a chromosomal arm map for wheat based on RFLP markers. Theor Appl Genet. 1992;83(8):1035–43.View ArticlePubMedGoogle Scholar
- 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.View ArticlePubMedGoogle Scholar
- Jewell E, Robinson A, Savage D, Erwin T, Love CG, Lim GA, et al. SSRPrimer and SSR taxonomy tree: Biome SSR discovery. Nucleic Acids Res. 2006;34 suppl 2:W656–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Van Ooijen JW. JointMap 4. Software for the calculation of genetic linkage maps in experimental populations. Wageningen: Kyazma BV; 2006.Google Scholar
- Li X, Liu C, Chakraborty S, Manners JM, Kazan K. A simple method for the assessment of crown rot disease severity in wheat seedlings inoculated with Fusarium pseudograminearum. J Phytopathol. 2008;156(11–12):751–4.View ArticleGoogle Scholar