Genetic basis of qualitative and quantitative resistance to powdery mildew in wheat: from consensus regions to candidate genes
BMC Genomics volume 14, Article number: 562 (2013)
Powdery mildew (Blumeria graminis f. sp. tritici) is one of the most damaging diseases of wheat. The objective of this study was to identify the wheat genomic regions that are involved in the control of powdery mildew resistance through a quantitative trait loci (QTL) meta-analysis approach. This meta-analysis allows the use of collected QTL data from different published studies to obtain consensus QTL across different genetic backgrounds, thus providing a better definition of the regions responsible for the trait, and the possibility to obtain molecular markers that will be suitable for marker-assisted selection.
Five QTL for resistance to powdery mildew were identified under field conditions in the durum-wheat segregating population Creso × Pedroso. An integrated map was developed for the projection of resistance genes/ alleles and the QTL from the present study and the literature, and to investigate their distribution in the wheat genome. Molecular markers that correspond to candidate genes for plant responses to pathogens were also projected onto the map, particularly considering NBS-LRR and receptor-like protein kinases. More than 80 independent QTL and 51 resistance genes from 62 different mapping populations were projected onto the consensus map using the Biomercator statistical software. Twenty-four MQTL that comprised 2–6 initial QTL that had widely varying confidence intervals were found on 15 chromosomes. The co-location of the resistance QTL and genes was investigated. Moreover, from analysis of the sequences of DArT markers, 28 DArT clones mapped on wheat chromosomes have been shown to be associated with the NBS-LRR genes and positioned in the same regions as the MQTL for powdery mildew resistance.
The results from the present study provide a detailed analysis of the genetic basis of resistance to powdery mildew in wheat. The study of the Creso × Pedroso durum-wheat population has revealed some QTL that had not been previously identified. Furthermore, the analysis of the co-localization of resistance loci and functional markers provides a large list of candidate genes and opens up a new perspective for the fine mapping and isolation of resistance genes, and for the marker-assisted improvement of resistance in wheat.
Bread wheat (Triticum aestivum L.) and durum wheat [Triticum turgidum (L.) subsp. turgidum (L.) convar. durum (Desf.)] are among the most important cultivated crops worldwide. Powdery mildew, which is caused by the fungus Blumeria graminis f.sp tritici, results in important yield losses and affects stable wheat production in areas with cool or maritime climates. Breeding of resistant cultivars is the most economical and environmentally sound method to decrease the use of fungicides and to reduce crop losses due to this disease.
Both qualitatively and quantitatively inherited resistances have been reported in wheat against powdery mildew. Qualitative resistance to powdery mildew is controlled by major race-specific genes that are effective only against some isolates of powdery mildew. Unfortunately, qualitative resistance is usually of short durability, due to frequent changes in the pathogen population . Consequently, new resistance genes are continuously needed to replace the defeated ones.
To date, more than 60 powdery mildew resistance genes/ alleles have been reported in common and durum wheat . Some of these were transferred from domesticated as well as wild relatives, such as Triticum turgidum var. dicoccoides (Körn.) and var. dicoccum (Schrank), T. timopheevii (Zhuk.), T. monococcum (L.), T. tauschii (Schmalh), and Aegilops speltoides (Tausch), or from more distant species, like Secale cereale (L.) . Molecular markers have largely been used for mapping a number of these genes to specific chromosomes or chromosome regions , and some of these genes have been cloned. In particular, Pm3b from hexaploid wheat is a member of the coiled-coil nucleotide binding site leucine-rich repeat (NBS-LRR) class of disease resistance genes, and it determines complete resistance . Very often partial resistance is a polygenic trait, but this is not always the case. An example of monogenic partial resistance is the gene Mlo. Homologs of the barley [Hordeum vulgare (L.)] gene Mlo were found in syntenic positions in all three genomes of hexaploid wheat [5–8]. The Mlo gene was isolated by positional cloning, and its deduced amino-acid sequence revealed no homologies to other characterized plant R proteins . This sequence defines an integral membrane protein with seven transmembrane helices and two casein kinase II motifs. A putative serine/ threonine protein kinase gene (Stpk-V) in the Pm21 locus was also characterized as conferring durable resistance; it is located on chromosome 6 V of Haynaldia villosa (L.) and was transferred to wheat as a 6VS · 6AL translocation .
The polygenic nature of partial resistance to powdery mildew makes it more complicated to handle in any breeding programs, compared to race-specific resistance. Many reports on high-density linkage maps and quantitative trait loci (QTL) that govern this trait are available in the literature, e.g., [11–13]. These QTL have been mapped against a single genetic background, and they have been evaluated in a limited number of environments. Moreover, it is difficult to predict the usefulness of QTL for marker-assisted selection based only on the QTL performance in an individual genetic background in any particular study. Goffinet and Gerber  proposed a method to combine the results from independent studies, called QTL meta-analysis, which allows the numbers and consensus positions of ‘real’ QTL involved in the control of a certain trait to be obtained. These QTL regions from the QTL meta-analysis (‘meta-QTL’, or MQTL) have refined, or consensus, positions, and they can facilitate the identification of positional candidate genes. In wheat, some examples of MQTL have been reported. Quraishi et al.  provided an overall view of MQTL for nitrogen use efficiency and grain dietary fiber content in bread wheat, with an unraveling of the candidate genes for these traits. Stable QTL that have provided MQTL have been defined also for yield, yield components, and crop height, which are very important traits in wheat breeding [3, 16]. Ear emergence and resistance against Fusarium head blight have also been studied in wheat, which resulted in increased precision in the QTL position estimations and the identification of molecular markers linked to them [17–20].
In the present study, we report on the MQTL for resistance to powdery mildew in wheat. The objectives were: 1) to identify new sources for resistance to powdery mildew in durum wheat; 2) to develop a high-density integrated map for the projection of resistance genes/ alleles and QTL, and to investigate their distribution in the wheat genome; 3) to identify MQTL across studies to validate results from independent studies and reliable markers for marker-assisted selection; and 4) to identify the homologous regions of MQTL in other crops using a comparative genomics approach, to provide candidate genes that drive this trait.
Results and discussion
QTL mapping for resistance to powdery mildew in the durum-wheat Creso × Pedroso RIL population
The durum-wheat cultivars Creso and Pedroso have been evaluated for different traits and have shown some interesting phenotypes, in particular for disease resistance. A major QTL for resistance to leaf rust was previously mapped on the long arm of chromosome 7B of Creso . In the present study, we performed a combined ANOVA in the recombinant inbred line (RIL) population for powdery mildew disease severity over two environments. The mean squares due to RILs, environments, and RIL × environment interactions were highly significant. The means of the parental lines, the means and ranges of the RIL population, and the variance components estimated for each trait over the two environments are reported in Additional file 1. The parents of the segregating population showed a significant, although small, difference, with higher disease severity observed for Pedroso. This difference was consistent for both of the environments, although the infection pressure was higher in Italy than in Spain. The continuous phenotypic frequency distributions across the two environments of the 123 RILs for disease severity score showed a normal distribution, which suggests a quantitative genetic basis for this trait (data not shown).
A total of 5 QTL were detected in the present study, on chromosomes 6A, 2B, 3B and 6B, with logarithm of odds (LOD) values between 3.0 and 5.0 (Additional file 2). The observed variability for the phenotype ranged from 10.6% to 18.5%. Interestingly, the two QTL named CP 3 and located at 16 cM on chromosome 6B could correspond to the same genomic region, which would be common to both environments. This idea is supported because the two QTL fall within the same region on the integrated map (Additional file 3: Figure S1). All of the other QTL identified in the present study were expressed only in one environment, probably due to differences in the natural populations of powdery mildew present in the two environments. The alleles of resistance to powdery mildew were contributed by both parents. In particular, based on the sign of the additive effects, the allele of Creso was effective in decreasing disease severity for QTL CP 2; for the other QTL, the resistant phenotype was contributed by the allele of Pedroso. The limited difference in resistance between the two parents justifies the possibility that some resistant alleles are also contributed by the cv. Pedroso. Furthermore, the low disease pressure registered in the two environments can allow to detect a qualitative resistance more than a quantitative one. Another possible explanation of our results is that a major QTL or other QTL contributed by the cv. Creso could be present in genomic regions that are not covered by molecular markers in the Creso x Pedroso genetic map. In this case they have not been detected in the present study, while only minor QTL are evident. The map positions of the QTL are reported in Additional file 3: Figure S1, except for QTL CP 4, for which the projection was not possible, as only one marker (Xbarc134 - 121.6 cM) was in common between the Creso × Pedroso map and the integrated map developed in the present study. As will be described in detail in the following paragraphs, the QTL CP 2, 3, 4, and 5 did not fall into MQTL regions, even if they were at a distance of <10 cM from other individual QTL. For this reason we cannot be sure that they correspond to regions not previously identified. Interestingly, QTL CP 1 on chromosome 6A identified a genomic region that is characterized by the absence of previously described QTL. Therefore, the results from the analysis of the Creso × Pedroso population have revealed a new source of quantitatively inherited resistance to powdery mildew.
Development of a wheat integrated map
A dense consensus map was obtained by merging the recent durum-wheat high-density consensus map developed by Marone et al.  with two well-saturated bread-wheat maps [23, 24]. The resulting integrated map represents a useful tool for MQTL, by starting from genetic loci that were identified in both tetraploid and hexaploid wheat. After that, a number of partial genetic maps related to small chromosomal regions were projected onto the integrated map; these regions contained genes/ QTL of resistance to powdery mildew, as revealed by the present study and previous studies. Furthermore, as a large proportion of diversity arrays technology (DArT) markers corresponding to expressed sequences was revealed to code for resistance proteins, such as NBS–LRR proteins or LRR kinases , the chromosomal regions from the genetic maps available in the literature and containing these DArT markers were also projected onto the map.
The final consensus map was composed of 3,618 markers (DArT, simple sequence repeat [SSR], expressed sequence tag [EST]-derived, sequence-tagged sites, and restriction fragment length polymorphism) and spanned a total map length of 3,723.9 cM (Table 1; Additional file 3: Figure S1; Additional file 4). A total of 1,260 markers were positioned on genome A, 1,876 on genome B, and only 482 on genome D, even if the total lengths were not very different (Table 1). The average distance between markers on the whole genome was 1 cM per marker, thus enabling the more precise location of the QTL.
Out of 3,618 markers mapped, 295 corresponded to putatively expressed sequences, and for most of these (234), a putative function could be proposed based on similarity searches (Additional file 4). In many cases, this function was related to plant responses to diseases. Also in this case, the A and B genomes were much more represented than the D genome. In particular, one PCR-based molecular marker was found to correspond to a NBS-LRR gene (Xcdo244 on 2B - 129.4 cM), and five other markers to three protein kinases (TC91851 on 5A - 85.4 cM -, Xbcd1088 on 5A - 95 cM -, and Xmag1759 on 7A - 206.6 cM) and two pathogenesis-related proteins (TC77302 on 3A - 151.8 cM - and TC92445 on 7A - 226.9 cM). Finally, out of 105 mapped DArT markers, 72 corresponded to NBS-LRR proteins (Additional file 3: Figure S1, purple) and 33 to protein kinases, with most of these containing a LRR domain (Additional file 3: Figure S1, brown). These markers were spread all over the genome, except for chromosomes 3D, 4B, 4D, and 5D. The number per chromosome was from 1 (3A and 5A) to 13 (2B and 7A). Eleven of the NBS-LRR genes showed similarity with proteins belonging to the subgroup of coiled-coil NBS-LRR genes (wPt-1862 on chromosome 1A - 45.1 cM -, wPt-6059 - 72.6 cM - and wPt-7421 - 81.9 cM - on chromosome 1D, wPt-6064 on chromosome 2A - 80.3 cM -, wPt-4368 on chromosome 2B - 118.7 cM -, wPt-4077 on chromosome 3A - 82.5 cM -, wPt-1516 - 11.7 cM - and wPt-8845 - 114.6 cM - on chromosome 3B, wPt-3729 on chromosome 4A - 177.6 cM -, wPt-4936 on chromosome 5B - 120.8 cM -, and wPt-0833 on chromosome 7D - 16.3 cM). Furthermore, four DArT markers were similar to NBS-LRR and also contained an additional BED zinc finger domain (wPt-2120 on 2B - 73.3 cM -, wPt-0950 on 2B - 93 cM -, wPt-0189 on 2B - 113.3 cM -, and wPt-4660 on 4A - 73.7 cM).
The sequences of the markers corresponding to R genes were searched against the database of the wheat separate chromosome arms promoted by the International Wheat Genome Sequencing Consortium (http://urgi.versailles.inra.fr/srs83/displayTool.do?toolName=BlastN), to study other genes eventually present in the same genomic clone. The results of the searches are reported in Additional file 5. In many cases, other genes were present, and in particular, genes corresponding to transposons and retrotransposons, exonucleases, acyl transferase, RNaseH, and other unknown proteins. In particular, transposons are recognized as one of the factors that drive gene duplication and amplification events that are at the basis of the evolution of this class of genes . However, in other cases, more copies of resistance genes were observed. The marker wPt-4077 (82.5 cM) corresponded to a genomic region on the long arm of chromosome 3A that contained two copies in tandem (at a distance of about 2000 bp) and in the same orientation as a NBS-LRR gene similar to the Rp1-like gene of Brachypodium distachion. There were also two copies of a NBS-LRR gene in the region of chromosome 4A, where the marker wPt-6303 (55.6 cM) was positioned. In this case, the sequences were similar to a gene of Zea mays, and they showed opposite orientation. The distance was again ca. 2000 bp. Similar cases were observed for the markers wPt-0992a (176.9 cM) and wPt-3729 (177.6 cM) on chromosome 4A, and wPt-2782 (83.4 cM) on chromosome 6D. Two genes coding for receptor-like kinases were found within the genomic clone corresponding to the marker wPt-8238 (70.6 cM) on chromosome 3B.
In some cases, three copies of a gene were found in the clone. This was the case for the markers wPt-1314 (chromosome 6D - 90.4 cM), which corresponds to a NBS-LRR gene, and wPt-5011 (chromosome 1B - 59.1 cM), which corresponds to a LRR receptor-like protein. For this marker, the genomic clone on chromosome 1B contained three genes, the first in opposite orientation with respect to the second and third ones. There was also a genomic clone in a homoelogous position on the long arm of chromosome 1A, and also in this case there were three copies of the gene, with the same orientation as the corresponding genes in the homoeologous locus.
In other cases, clusters of resistance genes of different types were observed. The marker wPt-2147 (44.8 cM) on the short arm of chromosome 2A fell into a region containing a NBS-LRR gene and a gene coding for a probable receptor-like kinase. On the same chromosome, the marker wPt-1601a (48 cM) corresponded to a genomic clone that contained a protein kinase and a NBS-LRR gene. R genes are organized in large families that are subjected to continuous evolution, because of their constant interactions with the evolving pathogens. They often occur in clusters at specific loci following gene duplication and amplification events, and these clusters can combine similar or different genes, as has been shown for a number of species . The results reported in the present study confirm this feature also in wheat.
An abundance of R genes in plant genomes has been described in many species. As an example, around 600 R gene loci were identified by Shang et al.  in a genome-wide comparison of the major class of R genes between the indica and japonica rice varieties. Knowledge of the map positions of these genes is of great importance for genetic studies that are aimed at the identification and fine mapping of resistance determinants. Furthermore, positioning these genes in genetic maps rich in molecular markers offers the chance to select many molecular markers closely linked to resistance determinants, to be used in marker-assisted selection programs . Forty DArT sequences that correspond directly or indirectly to NBS-LRR genes and protein kinases were mapped in durum wheat in a previous study, and many of these were positioned in regions where determinants of resistance to a number of plant pathogens have been identified . The construction of a dense integrated map allowed us to extend the number of markers corresponding to R genes mapped on the wheat genome: a total of 117 genes were mapped (104 markers positioned on the map, plus 13 additional genes found in the same genomic clone of some markers), which comprised 79 NBS-LRR genes and 38 protein kinases, most of which contained LRR domains.
Meta-analysis provides the estimates of QTL consensus positions and the confidence intervals (CIs) where the causative polymorphisms of the QTL are most likely to be positioned . MQTL have been used for different traits in wheat and other crops, such as for wheat-grain morphology , for digestibility and cell-wall traits in maize , for grain yield in rice  and for multiple disease resistance in barley .
In total, this analysis used 20 publications that described 96 QTL for powdery mildew resistance, plus the 5 QTL described in the present study (101 in total). The characteristics of the phenotyping experiments and mapping populations are reported in Additional file 6. The chromosomal locations of 52 resistance genes/ alleles collected from 48 different mapping populations are also reported in this study (Additional file 3: Figure S1; Additional file 7). Three additional genes for resistance to powdery mildew have been published over the last few months, and therefore they are not included in the present analysis, which was already finalized at that stage. The gene Ml5323 has been mapped to chromosome 2B  in a region that corresponds to the gene MlIW170 that was mapped in the present study. Moreover, another gene (Pm46) has been reported in the literature, in the same region as Pm45, under the individual QTL QPm.inra-5D. The last gene, Pm24b, was reported as tightly associated to Pm24 by Xue et al. , which was mapped in the present study on chromosome 1D.
The analysis resulted in 24 MQTL, each of which resulted from integration of at least two initial QTL. Thirty-seven out of the initial 101 QTL remained as singletons. The MQTL, along with their Akaike Information Criterion values, CIs, flanking markers, and number of initial QTL involved are reported in Table 2. The number of clustered initial QTL ranged from 2 to 6, whereas the 95% CI of the MQTL varied from 0.2 to 48.5 cM, with an average of 7.3 cM. MQTL were found on all chromosomes, except for 1D, 3B, 3D, 4D, 6B, and 6D, and their number per chromosome ranged from 1 to 3 (Additional file 3: Figure S1, Table 2).
In the present study, resistance ‘hot spots’ were highlighted, where the QTL detected in different studies are located within the same genomic region, together with the Pm genes, such as on chromosomes of groups 1, 2, 4, 5, and 7. The chromosomes that are richer in loci for resistance to powdery mildew are: chromosome 2B, with three MQTL, one individual QTL and 6 Pm genes, with two of these positioned under MQTL8; and chromosome 7A, with two MQTL, two individual QTL, and the highest number of Pm genes (13). No genes were mapped on chromosomes of group 4 and on 1B, 3A and 3D.
In more detail, the region comprised between MQTL1 and MQTL2 on chromosome 1A is characterized by the presence of four genes (Pm3g, Pm3a, Pm3e, and Mlar), one of which (Pm3g) is under MQTL1, and the others are a few cM from MQTL1 and MQTL2. On chromosome 2A, MQTL4 co-maps with PmHNK54, while Pm4b is under MQTL5. Four genes (PmJM22, PmPS5B, MlAB10, and MlZec1) are located in the region of MQTL8 on chromosome 2B, but only the first one co-maps. On chromosome 2D, the gene Pm43 is localized near (1.4 cM from) MQTL9. Chromosome 5D has MQTL18, which includes PmY212 and two other genes (Pm35 and PmY201) that are mapped very close to MQTL19 (less than 1 cM away). The homoeologous group 7 is characterized by gene/ MQTL co-localization on all of the three chromosomes. MQTL22 on chromosome 7A includes NCA6Pm. Two genes (Pm5d and PmTm4) are under MQTL23 on chromosome 7B, and another two (mlxbd and Pm5e) are at very short distance from MQTL23. Finally, the gene Lr34/Yr18/Pm is included in MQTL24 on chromosome 7D.
Genes mapped at very small distances can represent allelic series at the same locus, or clusters of R genes. The three Pm3 genes on chromosome 1A might represent an allelic series at the Pm3 locus. Similar cases are Pm4 (chromosome 2A) with three alleles, and Pm5 (chromosome 7B) with two alleles. On chromosome 7A, two genes, Mlm80 and Mlm2033, are mapped at a distance of 0.7 cM from each other on the consensus map, and Yao et al.  described these as two alleles of the same locus; however, in other cases, allelism tests might be necessary to indicate either different alleles at the same locus or tightly linked resistance genes, even with different resistance specificities. Schweizer and Stein  recently reported that loci for resistance to different diseases can associate in the same locus in barley. We found similar results in wheat, as we can identify the regions on chromosomes 4B and 4D that correspond to TaMlo loci, homologs of barley Mlo, based on common markers between our map and that developed by the International Triticeae Mapping Initiative. In particular, these correspond to the MQTL14 region on chromosome 4B and the single QTL QPm.caas on chromosome 4D. For chromosome 4A, two DArT markers that correspond to the barley NBS2-RDG2a gene for resistance to leaf stripe (wPt-0992 - 176.9 cM ; ) and to a rust resistance Rp1-like protein of T. aestivum (wPt-3729 - 177.6 cM ; ) are positioned in the region corresponding to MQTL13 on chromosome 4A. The Pm1 gene is associated with the leaf and stem rust resistance genes Lr20 and Sr15 on chromosome 7A . Many genes/ alleles are also located under or very close to a single QTL, as for chromosomes 1A, 1D, 2A, 3B, 5A, 5B, 5D, 6A, 6B, 6D, 7A, and 7B (Additional file 3: Figure S1). Resistance determinants are generally located in the telomeric and sub-telomeric regions, according to Schmolke et al. .
The chromosomal regions involved in the control of powdery mildew resistance were also investigated to find candidate genes for this trait, by searching in public databases for possible functions for these sequenced markers that putatively correspond to expressed genes. Many functional categories were identified, as annotated in Additional file 4. Most of the candidate genes encode NBS-LRR proteins, LRR-protein kinases, receptor-like protein kinases, WRKY or NAC transcription factors, pathogenesis-related proteins, cell transporters, or glutathione transferases (involved in the metabolism of reactive oxygen species), or proteins involved in lipid metabolism, amino-acid biosynthesis or cell wall modifications upon powdery mildew attack. As examples, on chromosome 1A there are genes that correspond to putative ketol-acid reductoisomerase, an enzyme induced by plant–pathogen interactions  near MQTL2 (Xcdo1160 - 53.6 cM) and under the individual QTL QPm_Lan (Xmwg632 - 80.3 cM; ). Three other genes of interest are a thioredoxin, 1,4-benzoquinone reductase, and a protein kinase (markers Swes578 - 51.6 cM, CA651264 - 55.9 cM -, and wPt-5011 - 59.1 cM, respectively); these are positioned under the individual QTL QPm.osu-1B. There are numerous markers that correspond to NBS-LRR genes positioned under individual QTL or MQTL (Additional file 8). These have particular importance as candidate genes. Indeed, some of these are similar to cloned R genes: the marker wPt-1862 that is mapped under MQTL1 on chromosome 1A (45.1 cM) is very near to the Pm3g allele and has 44% identity at the amino-acid level with the Pm3 gene. Three DArT markers, wPt-1912, wPt-4107, and wPt-1560, are mapped within the region of MQTL3 on chromosome 1B (15 cM, 16.1 cM, and 16.5 cM, respectively); the first two of these putatively correspond to NBS-LRR proteins, and the third to a protein kinase (Table 3). The marker wPt-1912 in particular shows 69% identity at the amino-acid level with the R gene Lr21 of T. aestivum. There are four other DArT markers that putatively code for NBS-LRR proteins within the region of MQTL13 on chromosome 4A (Table 3). One of these, wPt-3729 (177.6 cM), shows 83% identity with the rust resistance Rp1-like gene of T. aestivum. Many genes of interest are associated with individual QTL on chromosome 7A. The markers wPt-4487b (36.2 cM) and wPt-7491b (43.1 cM) (LRR protein kinases), wPt-6966 (36.6 cM), wPt-3648 (41 cM), and wPt-3434 (43.6 cM) (NBS-LRR proteins), and wPt-1441 (37.9 cM) (acyl-protein synthetase) are positioned under QTL-7A on the short arm. The marker wPt-6966, in particular, shows 49% identity at the amino-acid level to the B. distachyon gene Lr21, for resistance to leaf rust. Examples like these make the NBS-LRR genes good candidates for the R genes and for the MQTL, and suitable markers to follow resistant phenotypes in a breeding program.
As well as the information on putative functions, candidate genes were also searched for on the basis of the expression data available in the PLEXdb database. The available nucleotide sequences of markers positioned on the integrated map developed in the present study were searched for the corresponding probe-set on the wheat 61 K chip (http://www.plexdb.org/modules/tools/plexdb_blast.php), and the expression of each probe-set was then evaluated in two experiments in which wheat plants were subjected to infection with powdery mildew (TA34, TA39; Additional file 4). Some probe-sets were up-regulated in response to powdery mildew infection, and in particular those that correspond to the markers: Xmag974 (155.9 cM), close to MQTL12 on chromosome 4; wPt-5766 (24 cM), under QPm.inra-5D; TC85303a (73.1 cM); and CA677684, under QPm.caas.6B; TC77994 (100.1 cM), very close to MQTL21 on chromosome 7A; and wPt-4902 (160.4 cM), under QPm_RF on chromosome 7B.
The expression of some other probe-sets was down-regulated by powdery mildew infection, as for markers: Xbcd1355 (87.9 cM), under MQTL15 on chromosome 5A; and TC70722 (158.9 cM), under QPm_RF on chromosome 7B. The probe-sets that correspond to markers TC85037 (76.4 cM) and TC65966 (78.2 cM) under QPm.caas.6B were up-regulated and down-regulated in the two experiments taken into consideration. The probe-sets that correspond to five markers also showed differential expression between the resistant and susceptible genotypes studied in the framework of experiment TA39 (TC85303a - 73.1 cM on 6B, TC85037 - 76.4 cM on 6B, CA677684b - 98.2 cM on 6B, TC77994 -100.1 cM on 7A, and TC70722 - 158.9 cM on 7B). More details are reported in Additional file 4.
Some chromosomes were characterized by the presence of resistance QTL and MQTL, but not of Pm genes, and therefore these appeared to bring only loci for quantitative resistance to powdery mildew, even if single genes responsible for qualitative resistance to other diseases have been reported previously, such as Lr28 and Lr30, and Sr7, and Stb7, for resistance to leaf rust, stem rust and septoria, respectively, on chromosome 4A .
The availability of markers with a known sequence offers an opportunity to expand the collinearity analysis across grass genomes, moving from the wheat integrated map. The sequences of markers that correspond to expressed genes were used as queries in BLASTX searches, to identify segments of collinearity with the Brachypodium and rice genomes (http://www.phytozome.net/). We report the rice sequences with the ID reported by Phytozome, and we also add in round brackets the new Os ID, when available, according to RAP-DB (http://rapdb.dna.affrc.go.jp/) to make it easy to find the right gene position. Additional file 4 includes the results of the collinearity analysis for all of the DArT markers and the EST-SSRs positioned under individual QTL or MQTL, and for which there was a significant match on the Brachypodium and rice genomes. Syntenic regions were identified on durum-wheat chromosomes 2A, 4A, 5A, 6A, 7A, 1B, 2B, 3B, 5B, 6B, and 7B. Overall, many of the wheat genomic regions identified showed correspondence only with the Brachypodium genome; in other cases, correspondence was maintained also with the rice genome.
On chromosome 3A, there was correspondence across nine markers mapped on the wheat integrated map (Xbcd22 - 27.2 cM, BJ213673 - 45 cM, wPt-2938 - 48.4 cM, Xbcd828 - 50.1 cM, TC74823 - 64.8 cM, wPt-2202 - 72.7 cM, wPt-6422 - 76.7 cM, wPt-4077 - 82.5 cM, and Xbcd372 - 100.6 cM) and regions on chromosomes 2 and 1 of Brachypodium and rice, respectively (Additional file 4). This region on the map comprised a DArT marker that corresponds to a NBS-LRR gene (wPt-4077). The region of chromosome 1 of rice comprised between the loci Os01g47540 (Os01g0665300) and Os01g67800 contains 1,430 genes with descriptions. Expression data were searched for the genes included in this region of rice chromosome 1, to identify genes that might be involved in the rice reaction to powdery mildew. In particular, the genes where there was a 4-fold change in expression during non-host interactions with Blumeria graminis f. sp. hordei were selected from the plant expression database (experiment OS92 - PLEXdb; http://www.plexdb.org/modules/PD_browse/experiment_browser.php?experiment=OS92). A total of 15 genes with a 4-fold change were identified (Table 4), which included a putative LRK1 protein, a MDR-like ABC transporter, a NPK1-related protein kinase, a putative RING-H2 zinc finger protein, a putative peroxidase, and an S-adenosylmethionine-dependent methyltransferase. Most of the selected sequences belong to gene families known to have key roles in plant responses to pathogen attack.
The same approach was used for a region on chromosome 2B comprised between the markers wPt-2120 (74.4 cM) and Xcdo244 (129.9 cM). There are corresponding regions on chromosome 5 (Bd5g03654-Bd5g27500) of Brachypodium and on chromosome 4 of rice [Os04g22100-Os04g59494(Os04g0691500)]. The rice genomic region contains 3,413 predicted genes, and based on expression data, 10 genes were identified with 4-fold changes, which included a pentatricopeptide repeat protein, a heavy-metal transport/ detoxification protein, a receptor-like protein kinase, and two retrotransposon proteins.
Three markers under MQTL14 on chromosome 4B, Xbcd110 (107.4 cM), wPt-7062 (109.6 cM), and wPt-9223 (112.1 cM), identified collinear regions on chromosome 1 of Brachypodium (Bradi1g68460-Bradi1g75960) and chromosome 3 of rice [Os03g04410 (Os03g0136900)-Os03g16310]. Also in this case, seven genes that were in this region were differentially expressed in response to powdery mildew in rice, including a cyclin-dependent protein kinase inhibitor, an Fe oxygenase, a CCAAT-binding transcription factor, a sugar transporter, a glyoxal oxidase, and a pentatricopeptide repeat protein.
When possible, information on the physical mapping of genes for resistance to powdery mildew was used for the collinearity search. Blanco et al.  assigned the Pm36 gene to chromosome bin 5BL6-0.29-0.76. A list of wheat ESTs was selected that were previously mapped to the same chromosomal region [44, 45]; these were used to search the corresponding loci in rice (Additional file 9), and strong correspondence with chromosome 9 was identified. All of the predicted genes contained within the interval Os09g01690 (Os09g0104300)-Os09g38790 (Os09g0560900) were inspected for expression data, and there were nine genes with 4-fold changes, including a protein belonging to the Ulp1 protease family, an ethylene-responsive protein containing an AP2/ERF domain, an F-box-like protein, blight-associated protein p12, and a ribosomal protein. Furthermore, in a very recent study, six ESTs that derive from a wheat ‘totipotent’ cDNA library (AJ609811, AJ610871, AJ611689, AJ614358, AJ716441, and FM208374) were differentially expressed in two durum-wheat near-isogenic lines that differ in their resistance to powdery mildew, and some of them (AJ610871, AJ611689, AJ716441, and FM208374) were physically mapped to the centromeric bin region of chromosome 5BL, where the Pm36 gene was localized previously . Based on similarity searches, these sequences appear to correspond to the EARLY RESPONSIVE TO DEHYDRATION 15 (ERD15) transcription factor (AJ611689), an ATPase (AJ716441), and the heat-shock protein HSP90C (FM208374).
Following the approaches described above, a list of candidate genes was obtained, with each gene showing three fundamental features: (i) a map position coincident with a gene/ QTL/ MQTL of resistance to powdery mildew; (ii) a putative function related to disease response, and (iii) modulation of the expression levels by the pathogen. In particular, there were proteins belonging to the receptor-like kinase (RLK) family, genes involved in cell-wall modifications, NBS-LRR proteins, genes belonging to the ‘secretion and transport’ category, WRKY transcription factors, ABC transporters, peroxidases, and proteins involved in lipid and amino-acid metabolism.
Homoelogous relationship between MQTL
Many markers that revealed multiple loci were positioned in the integrated map elaborated in the present study. The detailed analysis of homoelogous and nonhomoeologous positions for each marker enables us to propose some hypotheses relating to the evolution of loci for resistance to powdery mildew in wheat.
The markers Xcfd15 (55.2 cM), Xwmc818a (57.4 cM), XksuG9 (68.3 cM), and Xgwm666 (72.8 cM), are positioned in the region of MQTL2 on chromosome 1A, and they show correspondence with homoelogous loci on chromosome 1B, under MQTL3. Furthermore, the markers XksuG9 (45.1 cM), Xgwm608 (55.2 cM), and Xbarc119 (48.6 cM) reveal a third locus on chromosome 1D, under the individual QTL QPm.inra.1D.1. Altogether, these data suggest a homoelogous relationship across MQTL2, MQTL3, and the QTL QPm.inra.1D.1 in group 1. The absence of homoelogous loci between the regions of MQTL1 and MQTL3 supports the results of the meta-analysis and that MQTL1 is separated from MQTL2. A group of markers (Xgwm165a - 13.9 cM; b - 105.2 cM, Xwmc89a - 15.5 cM; b - 106.6 cM, Xwmc491a − 15.8 cM; b - 107.2 cM, Xwmc617b - 28.5 cM; a - 103.8 cM, and Xgwm192a - 38.5 cM; b − 107 cM - between 4A and 4B; and Xwmc617a; c - 0.0 cM, Xgwm165b; c - 38 cM, Xgwm608d - 105.5 cM; e - 31 cM, Xwmc206b - 106.6 cM; e - 54 cM, Xcfd39a - 103.1 cM; b - 57 cM, and Xwmc89b; c - 27 cM - between 4B and 4D) identified loci in homoelogous positions on chromosomes 4A, 4B and 4D, which suggests that there might be a homoelogous relationship across MQTL11, MQTL14, and the individual QTL QPm.caas.4D on chromosomes 4A, 4B, and 4D, respectively. The same relationship was found for MQTL20 and the individual QTL PMm on chromosomes 6A and 6B, respectively, based on markers Xabc175 (b - 81.9 cM; a - 98.1 cM), Xwmc179 (b - 64.7 cM; c - 95.5 cM), wPt-3191 (a - 114.6 cM; b - 98.2 cM), Xwm417 (a - 121.7 cM; b - 111.9 cM), Xcdo1091 (a - 125.4 cM; b - 104.3 cM), and Xpsr546 (a - 138.2 cM; b - 108 cM), and between MQTL22 and MQTL23 on chromosomes 7A and 7B, respectively, based on markers Xgwm746 (a - 172.7 cM; b - 110.1 cM), wPt-6869 (a - 200.7 cM; b - 148.9 cM), wPt-3439 (a - 202.5 cM; b - 151.2 cM), Xcfa2040 (b - 202.7 cM; a - 151.5 cM), Xgwm344 (b - 208.8 cM and c 153 cM; a - 234.6 cM and d - 182.7 cM), and Xgwm1061 (a - 210.7 cM; b - 162.1 cM). In all of these cases, resistance determinants might have originated before diversification of homoelogous genomes.
The case of MQTL6 was different, which was mapped on chromosome 2B. No QTL or Pm genes were identified in the homoelogous region on chromosome 2A based on the markers Xgwm630 (a - 62.1 cM; b - 81.3 cM), XksuF37 (a - 61.3 cM; b - 82.5 cM), Xwmc474,(a - 56.1 cM; b - 82.7 cM), and wPt-2120 (a - 57.4 cM; b - 73.3 cM). In this case, there might be a homoelogous QTL, although it has not been mapped yet, or a resistance source might have evolved after diversification of the homoelogous genomes, or even before, but followed by the loss of the resistance determinant on chromosome 2A.
Groups of multi-locus markers that reveal loci mapped on nonhomoeologous positions were considered to be involved in putative translocations, as described by Marone et al. . In particular, the translocations that involve chromosome arms 4AL, 5AL, and 7BS have been firmly established [22, 47–56]. At the diploid level, an exchange of the terminal segments of chromosomes 4AL and 5AL took place, followed in tetraploid wheat by the exchange of the distal portion of chromosome 5A segment on chromosome 4AL with a terminal segment from chromosome 7BS. In the present study, a group of seven markers reported in Figure 1 identified the loci on the short arm of chromosome 7AS and on the long arm of chromosome 4A, instead of chromosome 7B; this allowed us to identify the translocation event that took place between homoelogous groups 4 and 7 on the integrated map developed in the present study. As previously shown , the order of the markers was inverted in the two regions. More interestingly, individual QTL and MQTL were positioned in the translocated region. In particular, a correspondence can be seen between the region of MQTL12 on chromosome 4A and MQTL21 on chromosome 7A. Furthermore, the position of the two loci revealed by the marker wPt-4487 (a - 173.7 cM; b - 36.2 cM) suggests a correspondence between MQTL13 on 4AL and the individual QTL-7A on 7AS. These results strongly suggest that the loci for resistance to powdery mildew evolved during very ancient times, before this translocation event took place in wheat. Local rearrangements might then have driven the evolution of new resistance specificities in the described loci.
Another translocation event is underlined by a group of markers that reveal loci on chromosomes 5B and 7B (Xbarc4a - 100.5 cM - and b - 80 cM, wPt-9814a - 102.9 cM - and b - 81.8 cM, wPt-5737a - 103.1 cM - and b - 81.7 cM, Xgwm68a - 107.4 cM - and b - 72.6 cM, wPt-1457a - 108.1 cM - and b - 78.4 cM, and Xbarc176a - 109.5 cM - and b - 89.3 cM). All of these markers are within the region of MQTL17 on chromosome 5B, and correspond to the individual QTL QPmV.inra.7B on chromosome 7B.
Analysis of sources of resistance to powdery mildew in wheat
The genetic background of resistant varieties, as the parents of the populations used to map genes and QTL responsible for powdery mildew resistance, was analyzed by searching for information in the Wheat Pedigree and Identified Alleles of Genes database (http://genbank.vurv.cz/wheat/pedigree/). A detailed investigation of the pedigree data of resistant varieties that is reported in Additional file 10 allowed us to hypothesize the origin of resistance and its occurrence in other cultivars characterized by common ancestors. Twenty resistance genes mapped in the present study are already known from the literature to be transferred from diploid and tetraploid relatives of wheat, including T. turgidum var. dicoccoides (MlIW170, Pm42, MlAB10, MlZec1, Pm41, Pm16, Ml3D232, PmAS846, Pm36, PmG3M, PmG16, and MlIW72), T. timopheevii (Pm27, Pm37, PmNCAG11, and MlAG12), T. carthlicum (Nevski in Kom.) (Pm6 and Pm33), T. monococcum (Pm4d, Pm2026, NCA6Pm, Mlm2033, and PmNCA4), T. boeoticum (Boiss) (PmTb7A.1, PmTb7A.2, and Mlm80), T. urartu (Tum.) (PmU) and several Aegilops spp. (PmY212, Pm35, PmY201, Pm34, and Pm12) or from related genera, such as Elytrigia (Pm40), Secale (PmHNK54 and PmTm4), and Thinopyrum intermedium (Host) (Pm43) (Additional file 10). In addition, there were MQTL that grouped single QTL that were contributed by the same wild ancestor. For example, MQTL2 combined the two QTL mapped in a population that derives from a cross between T. militinae (Zhuk.) and a wheat cultivar [57, 58], and one QTL where the resistance resulted from a variety that had T. timopheevii in its genetic background. As T. militinae derives from T. timopheevii, the resistance source might be the same. MQTL3, 13, 16, and 17 were also characterized by T. timopheevii-derived resistance. In contrast, MQTL that derive specifically from T. aestivum were also identified, such as MQTL1 (Additional file 10). This demonstrates that it is still possible to exploit the cultivated wheat to search for new resistance genes.
The coincidence of single QTL that co-map with resistance genes that derive from the same species was verified, as for chromosome 1D on which Pm24 was located under the QTL QPm.inra , where both derive from T. aestivum. An analogous case was on chromosome 6B, where Pm27 and the QTL PmM  co-map, and both derive from T. timopheevii.
The convergence in the same region of QTL and genes that derive from different species was the most represented case along the chromosomes of the wheat genome, which leads to the hypothesis that different sources of resistance might have contributed to these regions. Chromosome 1A was characterized by the presence of introgression from T. militinae, which started from the marker Xpsp2999 (47.2 cM), as proposed by Jakobson et al. , where MQTL1 and MQTL2 were located, also including Pm3g, Pm3a, Mlar, and Pm3e, where their resistance derives instead from T. aestivum. A gene-rich region was identified on chromosome 2A, corresponding to MQTL4. This region appears to have been subjected to multiple introgression events from wild species, as genes and QTL that derive from different sources were mapped (T. aestivum, T. dicoccum, T. monococcum, T. timopheevii, and T. carthilicum). There are other intriguing regions in which different sources of resistance contributed to the same chromosomal region on chromosomes 2B (T. aestivum, T. carthlicum, and T. dicoccoides), 5A (T. aestivum, T. monococcum, and T. militinae), 5B (T. aestivum and T. dicoccoides), and 6B (T. aestivum, Ae. Speltoides, T. timopheevi, and T. dicoccoides). In particular, the long arms of chromosomes 7A and 7B are characterized by a large number of genes/ QTL with different resistant sources, which are localized at small distances. Finally, the analysis of pedigree performed in the present study reveals that common wheat varieties recur in the genetic background of most of the resistant varieties used to map powdery mildew resistance genes/ QTL. In particular, the American cultivar ‘Mediterranean’, the Japanese landrace ‘Akakomugi’, the French cultivar ‘Chiddam d’automne’ and the Japanese cultivars ‘Norin10’ and ‘29’ recur in most of the pedigrees analyzed (Additional file 10).
This QTL meta-analysis has allowed us to reduce the CIs of the resulting MQTL relative to those of the initial QTL by a coefficient of reduction of 3.2 (mean initial CI/ mean MQTL CI), thus facilitating the search for candidate genes and providing markers that are more closely associated with the corresponding MQTL.
Meta-analysis of QTL for powdery mildew resistance is an effective approach to identify precise consensus QTL, which allows the confusion that exists due to redundancy in the number of QTL in overlapping genomic regions to be overcome. The development of an integrated map that is very dense and rich in markers that correspond to expressed genes that have putative roles in plant responses to pathogens represents a useful resource for the analysis of multiple components of resistance to diseases. The 24 MQTL identified that show small genetic intervals represent an important tool that can be used for marker-assisted selection/ pyramiding in wheat-breeding programs, or for map-based cloning. The presence of clusters of genes in the regions involved in the control of this trait indicates that a well-coordinated response of many genes is fundamental to achieve pathogen resistance. The comparative genomics approach to identify the consistency of QTL for resistance against B. graminis across grass genomes reveals the conservation and evolutionary significance of some of these loci. The information generated in this study will be of great use for future studies aimed to improving the powdery mildew resistance in wheat.
QTL mapping for powdery mildew resistance
The Creso × Pedroso mapping population was composed of 123 RILs and was evaluated during the growing season of 2005–2006 in Córdoba, southern Spain, and in Foggia, southern Italy, two of the major durum-wheat-growing areas in the Mediterranean basin. These 123 RILs and the parental cultivars Creso and Pedroso were sown in single 1-m-long rows, with 25 cm between rows, and organized in a randomized block design with three replications. No artificial inoculations were performed, as heavy infections of powdery mildew usually occur in these areas. All of the recommended agronomic practices were followed according to local standards. The disease severity was estimated visually when the disease was maximally spread, as the percentage of leaf area covered by powdery mildew in adult plants, at the level of the whole plant canopy. The parents Creso and Pedroso were characterized by a nearly identical heading time (data not shown), and the small differences observed across the segregating populations did not produce relevant effects of growth stage on the evaluation of disease severity. Phenotypic data were used for QTL mapping, together with the Creso × Pedroso genetic map that was previously developed [22, 25]. Genome-wide QTL searches were conducted using the MapQTL software package, version 5.0 , using both the simple interval mapping and the multiple QTL mapping functions. The LOD profiles from simple interval mapping were examined, and the marker closest to each LOD peak was selected as the cofactor to perform the multiple QTL mapping analysis. The LOD significance threshold levels of the respective traits and the corresponding map were calculated with the permutation test option provided in MapQTL, using 10,000 permutations. The LOD threshold was 3.0.
Bibliographic collection of QTL/ gene mapping
Twenty-three previously published studies were identified that reported on QTL for resistance against powdery mildew in bread wheat and durum wheat. Out of these, 20 were based on 19 different segregating populations, and these provided sufficient information on mapping and QTL characteristics for map projection and the MQTL (Additional file 6). Moreover, the information on the mapping data of resistance genes was collected from 58 published reports. Out of these, the projection of 52 genes was performed according to the common markers between the integrated map and the single genetic maps. For each gene, we defined the mapping population and the source of resistance, when known (Additional file 7).
Development of a consensus map
An integrated genetic map of the A, B, and D wheat genomes was constructed using a bread-wheat consensus map developed by Crossa et al.  as the reference map, on which two other consensus maps were projected: the Somers consensus map  and the most recent durum-wheat consensus map . All of the calculations for the creation of the integrated map were performed with the Biomercator software, version 2.1 . Maps of chromosomes with fewer than two common markers to the reference map were excluded before the creation of the consensus map. The marker order was finally verified according to the single genetic maps, and possible inversions were filtered out by discarding inconsistent loci. Genetic maps containing genes and QTL for resistance to powdery mildew were also integrated, to determine their positions on the consensus map. Moreover, the available wheat chromosomal regions containing DArT markers for which a function associated to plant pathogenesis was previously suggested  were collected [62–71]; these were integrated into the consensus map to facilitate the search for candidate genes.
QTL projection and meta-analysis
Quantitative trait loci were projected onto the consensus map using the Biomercator software. The 95% CIs of the initial QTL on their original maps were calculated through the empirical formula proposed by Darvasi and Soller , and extended by Guo et al.  according to the population type. For backcross, F2, and double-haploid populations, the QTL CI was 530/ N × R2, where N is the population size, and R2 is the proportion of the phenotypic variance explained by the QTL. The formula CI = 163/ N × R2 was used when a RIL population was considered. This approach is important, to assess the CIs using the same method for each QTL and study, and to estimate these when they were not published. In some cases, the calculated CI for a single QTL did not correspond to that reported in the literature, probably due to the small population size used to map the QTL, the very low percentage of explained variability, the method used for the QTL detection, or a sparse genetic map, which thus resulted in a broader interval on the integrated map [74, 75]. Co-localizing QTL that derive from the same experimental population used in different studies were considered if the fungal race or plant growth stage analyzed were not identical. This reduced the pre-selected QTL (148) to 96 (Additional file 6).
Meta-analysis was carried out on the QTL cluster on each chromosome separately, and the MQTL were obtained from the midpoint positions of the overlapping QTL. For n individual QTL, the Biomercator software tests the most likely assumption between 1, 2, 3, 4, and n QTL. The Akaike Information Criterion (AIC) was considered to select the best QTL model that indicated the number of MQTL. The model with the lowest AIC value was considered the best fit . When the n-model was the most likely model, the meta-analysis was performed again, choosing a subset of the QTL.
Analysis of collinearity
Collinearity was investigated for all of the markers positioned on the integrated map for which a gene sequence was available. The sequences were used as queries in a BLASTX search against the Brachypodium distachion and rice protein sequences available in the Phytozome database (http://www.phytozome.net). Information regarding the expression of the rice genes comprised in synthenic regions was retrieved from the Plant Expression Database (PLEXdb; http://www.plexdb.org/modules/PD_browse/experiment_browser.php?experiment=OS92).
Hsam SLK, Zeller FJ: Breeding for powdery mildew resistance in common wheat Triticum aestivum L.). The powdery mildews: A comprehensive treatise, Am. Phytopath. Soc. Edited by: Berlanger RR, Bushnell WR, Dik AJ, Carver DL. 2002, MN USA: St. Paul, 219-238.
Alam A, Xue F, Wang C, Ji C: Powdery mildew resistance genes in wheat: identification and genetic analysis. J Mol Biol Res. 2011, 1: 20-39.
Zhang H, Guan H, Li J, Zhu J, Xie C, Zhou Y, Duan X, Yang T, Sun Q, Liu Z: Genetic and comparative genomics mapping reveals that a powdery mildew resistance gene Ml3D232 originating from wild emmer co-segregates with an NBS-LRR analog in common wheat (Triticum aestivum L.). Theor Appl Genet. 2010, 121: 1613-1621. 10.1007/s00122-010-1414-6.
Yahiaoui N, Srichumpa P, Dudler R, Keller B: Genome analysis at different ploidy levels allows cloning of the powdery mildew resistance gene Pm3b from hexaploid wheat. Plant J. 2004, 37: 528-538. 10.1046/j.1365-313X.2003.01977.x.
Elliott C, Zhou F, Spielmeyer W, Panstruga R, Schulze-Lefert P: Functional conservation of wheat and rice mlo orthologs in defense modulation to the powdery mildew fungus. Mol Plant-Microbe Interac. 2002, 15: 1069-1077. 10.1094/MPMI.2002.15.10.1069.
Salmeron JM, Weislo LJ, Strawn LJ, Kramer CM, Wang HX, Vernooij BT, Levin JZ, Heifetz PB, Patton DA, Que Q: Mlo-genes controlling diseases. 2000, U.S: Patent WO, 0078799-A 3-
Niu J, He D: Molecular basis of powdery mildew resistance in wheat (Triticum aestivum L.). Afr J Biotechnol. 2009, 8: 4708-4716.
Konishi S, Sasakuma T, Sasanuma T: Identification of novel Mlo family members in wheat and their genetic characterization. Genes Genet Syst. 2010, 85: 167-175. 10.1266/ggs.85.167.
Büschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, van Daelen R, van der Lee T, Diergaarde P, Groenendijk J, Töpsch S, Vos P, Salamini F, Schulze-Lefert P: The barley Mlo gene: a novel control element of plant pathogen resistance. Cell. 1997, 88: 695-705. 10.1016/S0092-8674(00)81912-1.
Cao A, Xing L, Wang X, Yang X, Wang W, Sun Y, Qian C, Ni J, Chen Y, Liu D, Wang X, Chen P: Serine/threonine kinase gene Stpk-V, a key member of powdery mildew resistance gene Pm21, confers powdery mildew resistance in wheat. Proc Natl Acad Sci U S A. 2011, 108: 7727-7732. 10.1073/pnas.1016981108.
Zhang KP, Zhao L, Hai Y, Chen GF, Tian JC: QTL mapping for adult-plant resistance to powdery mildew, lodging resistance and internode length below spike in wheat. Acta Agron Sin. 2008, 34: 1350-1357.
Lan C, Ni X, Yan J, Zhang Y, Xia X, Chen X, He Z: Quantitative trait loci mapping of adult-plant resistance to powdery mildew in Chinese wheat cultivar Lumai 21. Mol Breed. 2010, 25: 615-622. 10.1007/s11032-009-9358-8.
Muranty H, Pavoine MT, Doussinault G, Barloy D: Origin of powdery mildew resistance factors in RE714, a wheat breeding line obtained from two interspecific crosses. Plant Breed. 2010, 129: 465-471.
Goffinet B, Gerber S: Quantitative trait loci: a meta-analysis. Genetics. 2000, 155: 463-473.
Quraishi UM, Abrouk M, Murat F, Pont C, Foucrier S, Desmaizieres G, Confolent C, Rivière N, Charmet G, Paux E, Murigneux A, Guerreiro L, Lafarge S, Le Gouis J, Feuillet C, Salse J: Cross-genome map based dissection of a nitrogen use efficiency ortho-metaQTL in bread wheat unravels concerted cereal genome evolution. Plant J. 2011, 65: 745-756. 10.1111/j.1365-313X.2010.04461.x.
Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Snape J: Meta-QTL analysis of the genetic control of crop height in elite European winter wheat germplasm. Mol Breeding. 2010, 29: 159-171.
Hanocq E, Laperche A, Jaminon O, Laine AL, Le Gouis J: Most significant genome regions involved in the control of earliness traits in bread wheat, as revealed by QTL meta-analysis. Theor Appl Genet. 2007, 114: 569-584. 10.1007/s00122-006-0459-z.
Griffiths S, Simmonds J, Leverington M, Wang Y, Fish L, Sayers L, Alibert L, Orford S, Wingen L, Herry L, Faure S, Laurie D, Bilham L, Snape J: Meta-QTL analysis of the genetic control of ear emergence in elite European winter wheat germplasm. Theor Appl Genet. 2009, 119: 383-395. 10.1007/s00122-009-1046-x.
Häberle J, Holzapfel J, Schweizer G, Hartl L: A major QTL for resistance against Fusarium head blight in European winter wheat. Theor Appl Genet. 2009, 119: 325-332. 10.1007/s00122-009-1041-2.
Löffler M, Schön CC, Miedaner T: Revealing the genetic architecture of FHB resistance in hexaploid wheat (Triticum aestivum L.) by QTL meta-analysis. Mol Breed. 2009, 23: 473-488. 10.1007/s11032-008-9250-y.
Marone D, Del Olmo AI, Laidò G, Sillero JC, Emeran AA, Russo MA, Ferragonio P, Giovanniello V, Mazzucotelli E, De Leonardis AM, De Vita P, Blanco A, Cattivelli L, Rubiales D, Mastrangelo AM: Genetic analysis of durable resistance against leaf rust in durum wheat. Mol Breed. 2009, 24: 25-39. 10.1007/s11032-009-9268-9.
Marone D, Laido’ G, Gadaleta A, Colasuonno P, Ficco DBM, Giancaspro A, Giove S, Panio G, Russo MA, De Vita P, Cattivelli L, Papa R, Blanco A, Mastrangelo AM: A high-density consensus map of A and B wheat genomes. Theor Appl Genet. 2012, 125: 1619-1638. 10.1007/s00122-012-1939-y.
Crossa J, Burgueño J, Dreisigacker S, Vargas M, Herrera-Foessel SA, Lillemo M, Singh RP, Trethowan R, Warburton M, Franco J, Reynolds M, Crouch JH, Ortiz R: Association analysis of historical bread wheat germplasm using additive genetic covariance of relatives and population structure. Genetics. 2007, 177: 1889-1913. 10.1534/genetics.107.078659.
Somers JD, Isaac P, Edwards K: A high-density microsatellite consensus map for bread wheat (Triticum aestivum L.). Theor Appl Genet. 2004, 109: 1105-1114. 10.1007/s00122-004-1740-7.
Marone D, Panio G, Ficco DBM, Russo MA, De Vita P, Papa R, Rubiales D, Cattivelli L, Mastrangelo AM: Characterization of wheat DArT markers: genetic and functional features. Mol Genet Genomics. 2012, 287: 741-753. 10.1007/s00438-012-0714-8.
Leister D: Tandem and segmental gene duplication and recombination in the evolution of plant disease resistance genes. Trends Genet. 2004, 20: 116-122. 10.1016/j.tig.2004.01.007.
Marone D, Russo MA, Laidò G, De Leonardis AM, Mastrangelo AM: Plant nucleotide binding site–leucine-rich repeat (NBS-LRR) genes: active guardians in host defense responses. Int J Mol Sci. 2013, 14: 7302-7326. 10.3390/ijms14047302.
Shang J, Tao Y, Chen X, Zou Y, Lei C, Wang J, Li X, Zhao X, Zhang M, Lu Z, Xu J, Cheng Z, Wan J, Zhu L: Identification of a new rice blast resistance gene, pid3, by genomewide comparison of paired nucleotide-binding site–leucine-rich repeat genes and their pseudogene alleles between the two sequenced rice genomes. Genetics. 2009, 182: 1303-1311. 10.1534/genetics.109.102871.
Truntzler M, Barrière Y, Sawkins MC, Lespinasse D, Betran J, Charcosset A, Moreau L: Meta-analysis of QTL involved in silage quality of maize and comparison with the position of candidate genes. Theor Appl Genet. 2010, 121: 1465-1482. 10.1007/s00122-010-1402-x.
Gegas VC, Nazari A, Griffiths S, Simmonds J, Fish L, Orford S, Sayers L, Doonan JH, Snape JW: A genetic framework for grain size and shape variation in wheat. Plant Cell. 2010, 22: 1046-1056. 10.1105/tpc.110.074153.
Swamy BPM, Vikram P, Dixit S, Ahmed HU, Kumar A: Meta-analysis of grain yield QTL identified during agricultural drought in grasses showed consensus. BMC Genomics. 2011, 12: 319-10.1186/1471-2164-12-319.
Schweizer P, Stein N: Large-scale data integration reveals colocalization of gene functional groups with Meta-QTL for multiple disease resistance in barley. Mol Plant-Microbe Interac. 2011, 24: 1492-1501. 10.1094/MPMI-05-11-0107.
Piarulli L, Gadaleta A, Mangini G, Signorile MA, Pasquini M, Blanco A, Simeone R: Molecular identification of a new powdery mildew resistance gene on chromosome 2BS from Triticum turgidum ssp. dicoccum. Plant Sci. 2012, 196: 101-106.
Gao H, Zhu F, Jiang Y, Wu J, Yan W, Zhang Q, Jacobi A, Cai S: Genetic analysis and molecular mapping of a new powdery mildew resistant gene Pm46 in common wheat. Theor Appl Genet. 2012, 125: 967-973. 10.1007/s00122-012-1886-7.
Xue F, Wang C, Li C, Duan X, Zhou Y, Zhao N, Wang Y, Ji W: Molecular mapping of a powdery mildew resistance gene in common wheat landrace Baihulu and its allelism with Pm24. Theor Appl Genet. 2012, 125: 1425-1432. 10.1007/s00122-012-1923-6.
Yao G, Zhang J, Yang L, Xu H, Jiang Y, Xiong L, Zhang C, Zhang Z, Ma Z, Sorrells ME: Genetic mapping of two powdery mildew resistance genes in einkorn (Triticum monococcum L.) accessions. Theor Appl Genet. 2007, 114: 351-358.
Neu C, Stein N, Keller B: Genetic mapping of the Lr20–Pm1 resistance locus reveals suppressed recombination on chromosome arm 7AL in hexaploid wheat. Genome. 2002, 45: 737-744. 10.1139/g02-040.
Schmolke M, Mohler V, Hartl L, Zeller FJ, Hsam SLK: A new powdery mildew resistance allele at the Pm4 wheat locus transferred from einkorn (Triticum monococcum). Mol Breed. 2012, 29: 449-456. 10.1007/s11032-011-9561-2.
Freitas-Astúa J, Bastianel M, Locali-Fabris EC, Novelli VM, Silva-Pinhati AC, Basílio-Palmieri AC, Targon MLPN, Machado MA: Differentially expressed stress-related genes in the compatible citrus-Citrus leprosis virus interaction. Genet Mol Biol. 2007, 30: 980-990. 10.1590/S1415-47572007000500026.
Lan CX, Liang SS, Wang ZL, Yan J, Zhang Y, Xia XC, He ZH: Quantitative trait loci mapping for adult-plant resistance to powdery mildew in Chinese wheat cultivar Bainong 64. Phytopathology. 2009, 99: 1121-1126. 10.1094/PHYTO-99-10-1121.
Chen Y, Hunger RM, Carver BF, Zhang H, Yan L: Genetic characterization of powdery mildew resistance in U.S. hard winter wheat. Mol Breed. 2009, 24: 141-152. 10.1007/s11032-009-9279-6.
Beavis WD: QTL analysis: Power, precision, and accuracy. Molecular dissection of complex traits. Edited by: Paterson AH. 1998, Boca Raton, FL: CRC Press, 145-161.
Blanco A, Gadaleta A, Cenci A, Carluccio AV, Abdelbacki AMM, Simeone R: Molecular mapping of the novel powdery mildew resistance gene Pm36 introgressed from Triticum turgidum var. dicoccoides in durum wheat. Theor Appl Genet. 2008, 117: 135-142. 10.1007/s00122-008-0760-0.
Sorrels ME, La Rota M, Bermudes-Kandianis C, Greene RA, Kantety R, Munkvold JD, Miftahudin NI, Mahmoud A, Ma X, Gustafson PJ, Qi LL, Echalier B, Gill BS, Matthews DE, Lazo GR, Chao S, Anderson OD, Edwards H, Linkiewicz AM, Dubcovsky J, Akhunov ED, Dvorak J, Zhang D, Nguyen HT, Peng J, Lapitan NL, Gonzalez-Hernandez JL, Anderson JA, Hossain K, Kalavacharla V, Kianian SF, Choi DW, Close TJ, Dilbirligi M, Gill KS, Steber C, Walker-Simmons MK, McGuire PE, Qualset CO: Comparative DNA sequence analysis of wheat and rice genomes. Genome Res. 2003, 13: 1818-1827.
La Rota M, Kantety RV, Yu JK, Sorrells ME: Nonrandom distribution and frequencies of genomic and EST-derived microsatellite markers in rice, wheat, and barley. BMC Genomics. 2005, 6: 23-10.1186/1471-2164-6-23.
Cifarelli RA, D’Onofrio O, Grillo R, Mango T, Cellini F, Piarulli L, Simeone R, Giancaspro A, Colasuonno P, Blanco A, Gadaleta A: Development of a new wheat microarray from a durum wheat totipotent cDNA library used for a powdery mildew resistance study. Cell Mol Biol Lett. 2013, 18: 231-248. 10.2478/s11658-013-0086-z.
Chao S, Sharp PJ, Worland AJ, Warham EJ, Koebner RMD, Gale MD: RFLP-based genetic maps of wheat homoeologous group 7 chromosomes. Theor Appl Genet. 1989, 78: 495-504. 10.1007/BF00290833.
Naranjo T: Chromosome structure of durum wheat. Theor Appl Genet. 1990, 79: 397-400. 10.1007/BF01186085.
Liu CJ, Devos KM, Chinoy CN, Atkinson MD, Gale MD: Non-homoeologous translocations between group 4, 5 and 7 chromosomes in wheat and rye. Theor Appl Genet. 1992, 83: 305-312.
Chen J, Gustafson JP: Physical mapping of genetically mapped molecular markers in homoeologous group 7 chromosomes of wheat by in-situ hybridization. J Jiangsu Agric Coll. 1994, 15: 1-9.
Chen J, Gustafson JP: Chromosomal rearrangement of wheat (T. aestivum) chromosome 4A by in situ hybridization. Chin J Genet. 1997, 24: 39-47.
Blanco A, Bellomo MP, Cenci A, De Giovanni C, D’Ovidio R, Iacono E, Laddomada B, Pagnotta MA, Porceddu E, Sciancalepore A, Simeone R, Tanzarella OA: A genetic linkage map of durum wheat. Theor Appl Genet. 1998, 97: 721-728. 10.1007/s001220050948.
Devos KM, Dubcovsky J, Dvorak J, Chinoy CN, Gale MD: Structural evolution of wheat chromosomes 4A, 5A and 7B and its impact on recombination. Theor Appl Genet. 1995, 91: 282-288. 10.1007/BF00220890.
Mickelson-Young L, Endo TR, Gill BS: A cytogenetic ladder-map of the wheat homoeologous group 4 chromosomes. Theor Appl Genet. 1995, 90: 1007-1011.
Nelson JC, Sorrels ME, Van Deynze AE, Lu YH, Atkinson M, Bernard M, Leroy P, Faris JD, Anderson JA: Molecular mapping of wheat: major genes and rearrangements in homoeologous groups 4, 5 and 7. Genetics. 1995, 141: 721-731.
Quarrie SA, Steed A, Calestani C, Semikhodskii A, Lebreton C, Chinoy C, Steele N, Pljevljakusic D, Waterman E, Weyen J, Schondelmaier J, Habash DZ, Farmer P, Saker L, Clarkson DT, Abugalieva A, Yessimbekova M, Turuspekov Y, Abugalieva S, Tuberosa R, Sanguineti MC, Hollington PA, Araguès R, Royo A, Dodig D: A high-density genetic map of hexaploid wheat (Triticum aestivum L.) from the cross Chinese Spring × SQ1 and its use to compare QTLs for grain yield across a range of environments. Theor Appl Genet. 2005, 110: 865-880. 10.1007/s00122-004-1902-7.
Jakobson I, Timofejeva HPL, Järve K: Adult plant and seedling resistance to powdery mildew in a Triticum aestivum × Triticum militinae hybrid line. Theor Appl Genet. 2006, 112: 760-769. 10.1007/s00122-005-0181-2.
Jakobson I, Reis D, Tiidema A, Peusha H, Timofejeva L, Valárik M, Kladivová M, Šimková H, Doležel J, Järve K: Fine mapping, phenotypic characterization and validation of non-race-specific resistance to powdery mildew in a wheat–Triticum militinae introgression line. Theor Appl Genet. 2012, 125: 609-623. 10.1007/s00122-012-1856-0.
Bougot Y, Lemoine J, Pavoine MT, Guyomar’ch H, Gautier V, Muranty H, Barloy D: A major QTL effect controlling resistance to powdery mildew in winter wheat at the adult plant stage. Plant Breed. 2006, 125: 550-556. 10.1111/j.1439-0523.2006.01308.x.
Lu Q, Bjørnstad Å, Ren Y, Asad MA, Xia X, Chen X, Ji F, Shi J, Lillemo M: Partial resistance to powdery mildew in German spring wheat ‘Naxos’ is based on multiple genes with stable effects in diverse environments. Theor Appl Genet. 2012, 125: 297-309. 10.1007/s00122-012-1834-6.
Van Ooijen JW: MapQTL5, Software for the mapping of quantitative trait loci in experimental populations. 2004, Kyazma B.V: Wageningen
Francki MG, Walker E, Crawford AC, Broughton S, Ohm HW, Barclay I, Wilson RE, McLean R: Comparison of genetic and cytogenetic maps of hexaploid wheat (Triticum aestivum L.) using SSR and DArT markers. Mol Genet Genomics. 2009, 281: 181-191. 10.1007/s00438-008-0403-9.
Peleg Z, Saranga Y, Suprunova T, Ronin Y, Röder MS, Kilian A, Korol AB, Fahima T: High-density genetic map of durum wheat x wild emmer wheat based on SSR and DArT markers. Theor Appl Genet. 2008, 117: 103-115. 10.1007/s00122-008-0756-9.
Båga M, Fowler DB, Chibbar RN: Identification of genomic regions determining the phenological development leading to floral transition in wheat (Triticum aestivum L.). J Exp Bot. 2009, 60: 3575-3585. 10.1093/jxb/erp199.
Jing HC, Bayon C, Kanyuka K, Berry S, Wenzl P, Huttner E, Kilian A, Hammond-Kosack KE: DArT markers: diversity analyses, genomes comparison, mapping and integration with SSR markers in Triticum monococcum. BMC Genomics. 2009, 10: 458-10.1186/1471-2164-10-458.
Munkvold JD, Tanaka J, Benscher D, Sorrells ME: Mapping quantitative trait loci for preharvest sprouting resistance in white wheat. Theor Appl Genet. 2009, 119: 1223-1235. 10.1007/s00122-009-1123-1.
Bovill WD, Horne M, Herde D, Davis M, Wildermuth GB, Sutherland MW: Pyramiding QTL increases seedling resistance to crown rot (Fusarium pseudograminearum) of wheat (Triticum aestivum). Theor Appl Genet. 2010, 121: 127-136. 10.1007/s00122-010-1296-7.
Genc Y, Oldach K, Verbyla AP, Lott G, Hassan M, Tester M, Wallwork H, McDonald GK: Sodium exclusion QTL associated with improved seedling growth in bread wheat under salinity stress. Theor Appl Genet. 2010, 121: 877-894. 10.1007/s00122-010-1357-y.
Tsilo TJ, Hareland GA, Simsek S, Chao S, Anderson JA: Genome mapping of kernel characteristics in hard red spring wheat breeding lines. Theor Appl Genet. 2010, 121: 717-730. 10.1007/s00122-010-1343-4.
Wenzl P, Suchánková P, Carling J, Šimková H, Huttner E, Kubaláková M, Sourdille P, Paul E, Feuillet C, Kilian A, Doležel J: Isolated chromosomes as a new and efficient source of DArT markers for the saturation of genetic maps. Theor Appl Genet. 2010, 121: 465-474. 10.1007/s00122-010-1323-8.
Lowe I, Jankuloski L, Chao S, Chen X, See D, Dubcovsky J: Mapping and validation of QTL which confer partial resistance to broadly virulent post-2000 North American races of stripe rust in hexaploid wheat. Theor Appl Genet. 2011, 123: 143-157. 10.1007/s00122-011-1573-0.
Darvasi A, Soller M: A simple method to calculate resolving power and confidence interval of QTL map location. Behav Genet. 1997, 27: 125-132. 10.1023/A:1025685324830.
Guo B, Sleper DA, Lu P, Shannon JG, Nguyen HT, Arelli PR: QTLs associated with resistance to soybean cyst nematode in soybean: meta-analysis of QTL locations. Crop Sci. 2006, 46: 595-602. 10.2135/cropsci2005.04-0036-2.
Schön C, Utz H, Groh S, Truberg B, Openshaw S, Melchinger A: Quantitative trait locus mapping based on resampling in a vast maize testcross experiment and its relevance to quantitative genetics for complex traits. Genetics. 2004, 167: 485-498. 10.1534/genetics.167.1.485.
Arcade A, Labourdette A, Falque M, Mangin B, Chardon F, Charcosset A, Joets J: BioMercator: integrating genetic maps and QTL towards discovery of candidate genes. Bioinformatics. 2004, 20: 2324-2326. 10.1093/bioinformatics/bth230.
This study was supported by the Italian Ministry of Agriculture (MiPAAF), with the special grants MAPPA 5A, CANADAIR and ESPLORA, and by the Italian Ministry of Education, Universities and Research (MIUR), with the special grant ISCOCEM. We are grateful to Dr. Christopher Berrie for scientific English language editorial assistance. We would also like to thank the IWGSC for providing early, pre-publication BLAST access to the survey sequences developed within the IWGSC Survey Sequence Initiative.
The authors declare no competing financial interests.
DM and MAR developed the integrated map and carried out the QTL meta-analysis. GL, AG and AMM carried out the collinearity study. DR developed the Creso × Pedroso segregating population. PDV and DR carried out the phenotyping of the Creso × Pedroso segregating population. RP and AB provided general guidance for the study. DM and AMM drafted the manuscript. AMM conceived and coordinated the study. All of the authors have read and approved this version of the manuscript.
Electronic supplementary material
Additional file 1:Phenotypic variation among the parental lines and RILs from the durum-wheat Creso × Pedroso population for disease severity caused by powdery mildew.(DOCX 22 KB)
Additional file 2:Results of the QTL analysis for resistance to powdery mildew in the Creso × Pedroso population.(DOCX 24 KB)
Additional file 3: Figure S1: Integrated map in wheat of the QTL and the MQTL identified by meta-analysis, for powdery mildew resistance. Vertical lines on the right of chromosomes indicate the confidence intervals, horizontal lines indicate the peak marker positions, where the length represents the percentage of variability explained by the QTL. The MQTL are in bold, with the single QTL in gray. The names of the QTL grouped in the same MQTL are in the same color. Pm genes are blue, markers that correspond to NBS-LRR proteins are purple, markers that correspond to kinases are brown. Markers where there is correspondence with probe-sets differentially expressed following powdery mildew infection are circled in black. (PDF 4 MB)
Additional file 4:The integrated map (in Excel format) with the markers for which a sequence is available and a putative function is annotated in public databases. The information given is for the functional markers positioned within the individual QTL and MQTL regarding the collinearity analysis and the correspondence with probe-sets of the wheat 61 K array. (XLS 1 MB)
Additional file 5:Results of the similarity search against the database of wheat separate chromosome arms promoted by the International Wheat Genome Sequencing Consortium.(XLSX 75 KB)
Additional file 10:Summary of the pedigree and sources of resistance for the QTL and genes mapped onto the integrated map.(XLS 58 KB)
About this article
Cite this article
Marone, D., Russo, M.A., Laidò, G. et al. Genetic basis of qualitative and quantitative resistance to powdery mildew in wheat: from consensus regions to candidate genes. BMC Genomics 14, 562 (2013). https://doi.org/10.1186/1471-2164-14-562