Discovery of genes affecting resistance of barley to adapted and non-adapted powdery mildew fungi
Genome Biology volume 15, Article number: 518 (2014)
Non-host resistance, NHR, to non-adapted pathogens and quantitative host resistance, QR, confer durable protection to plants and are important for securing yield in a longer perspective. However, a more targeted exploitation of the trait usually possessing a complex mode of inheritance by many quantitative trait loci, QTLs, will require a better understanding of the most important genes and alleles.
Here we present results from a transient-induced gene silencing, TIGS, approach of candidate genes for NHR and QR in barley against the powdery mildew fungus Blumeria graminis. Genes were selected based on transcript regulation, multigene-family membership or genetic map position. Out of 1,144 tested RNAi-target genes, 96 significantly affected resistance to the non-adapted wheat- or the compatible barley powdery mildew fungus, with an overlap of four genes. TIGS results for QR were combined with transcript regulation data, allele-trait associations, QTL co-localization and copy number variation resulting in a meta-dataset of 51 strong candidate genes with convergent evidence for a role in QR.
This study represents an initial, functional inventory of approximately 3% of the barley transcriptome for a role in NHR or QR against the powdery mildew pathogen. The discovered candidate genes support the idea that QR in this Triticeae host is primarily based on pathogen-associated molecular pattern-triggered immunity, which is compromised by effector molecules produced by the compatible pathogen. The overlap of four genes with significant TIGS effects both in the NHR and QR screens also indicates shared components for both forms of durable pathogen resistance.
Plant-pathogen co-evolution has shaped a multifaceted innate immunity system triggered by the recognition of non-self-molecules via pathogen recognition receptors (PRRs) belonging to the family of receptor-like kinases (RLKs) . These non-self-molecules known as pathogen-associated molecular patterns (PAMPs) or, more generally, microbe-associated molecular patterns (MAMPs) include conserved domains of proteins such as bacterial flagellin (flg22) or chitin fragments from fungal cell walls . PAMP-triggered immunity (PTI) has been recognized as the most ancient type of plant defense sharing also components with the innate immunity system of vertebrate and invertebrate animals. Downstream of PRRs its molecular components include MAP kinases, WRKY transcription factors as well as an arsenal of downstream-responsive, (WRKY-regulated) genes encoding proteins that generate reactive oxygen species, reinforce, and break down plant and pathogen cell-walls, respectively, or catalyze the synthesis of pathogen-toxic compounds such as phytoalexins. On top of PTI plants can activate an effector-triggered immunity (ETI) response that is based on the direct or indirect recognition of avirulenve (Avr) effector molecules of some pathogen races by major R-genes encoding nucleotide-binding leucine-rich repeat (NB-LRR) proteins, and on the initiation of a very strong local defense response often culminating in host-cell death. One of the preferred targets of effectors are PRRs, which have been found to be guarded by several NB-LRR type or PRR-like proteins therefore also being involved in ETI . Durable and broad-range non-host resistance (NHR) to virtually all races of non-adapted pathogens appears to be an important manifestation of PTI in many cases , although there is also experimental evidence that NHR can - at least in grass species - be mediated by as little as one major R gene recognizing an indispensable Avr effector . Race-specificity of NHR QTL to non- or only partially adapted fungal pathogens has also been described, similar to QTL for host quantitative resistance (QR) that is another manifestation of PTI  and QR is also referred to as race-non-specific or horizontal resistance -. However, in contrast to the very robust NHR response, QR is often not very efficient suffering from effector-triggered susceptibility (ETS) brought about by small secreted proteins or peptides from adapted pathogens that are active in the plant apoplast or inside host cells . The introgression of single major R genes usually confers strong protection against specific adapted pathogen races carrying the matching avirulence (Avr) effector genes, but the trait is often overcome by rapidly evolving new pathogen races with mutated Avr effectors acting in concert with other functionally redundant effectors. In principle, QR could also be mediated by partially functional (defeated) major R-genes weakly recognizing ubiquitous Avr effectors such as ECP1 or ECP2 , but molecular evidence for this type of interactions is scarce ,.
Barley (Hordeum vulgare ssp. vulgare) is an important crop plant and exhibits genetic variability determining to which extent it is successfully colonized by a given pathogen. This opens up the possibility to improve QR as a quantitative trait by introgressing and/or combining resistance-related alleles. Often, however, QR was found to be inherited by many QTLs making the trait difficult to handle in breeding practice due to complex crossing schemes, phenotype scoring ambiguities and linkage drag problems . Knowing the genes that encode important QR components in crop plants would render targeted QR improvement by allele mining and gene marker-assisted as well as pathway-oriented introgression more efficient. One of the major diseases of barley is powdery mildew caused by the obligate biotrophic ascomycete fungus Blumeria graminis f.sp. hordei (Bgh)  that also fulfills several criteria of a model plant-pathogen interaction due to a large body of physiological, cellular, biochemical, and molecular information on changes in the host during compatible or resistant interactions -. Transient expression and gene-silencing assays such as transient-induced gene silencing (TIGS) in bombarded epidermal cells have been developed over the years and proven to be valuable tools for a better understanding of barley/powdery mildew interactions -. NHR of barley against non-adapted formae speciales or species of powdery mildew such as the wheat pathogen B. graminis f.sp. tritici (Bgt) efficiently blocks fungal penetration attempts at the epidermal cell wall and shares a number of genes with the QR pathway triggered by the adapted barley powdery mildew fungus -. Therefore, the discovery of important genes required for NHR may also reveal their importance in QR, besides opening up the fascinating option of replacing effector-suppressed defense factors by their non-host orthologues by wide crosses or in transgenic plants.
Here we describe a large-scale functional approach in barley for the identification of genes that are relevant for NHR and QR to Bgt and Bgh, respectively . By using TIGS we screened three groups of host genes that were: (1) previously found to be upregulated in powdery mildew-attacked barley epidermal cells ,,,; (2) belonging to selected multigene-families; or (3) localizing within the confidence interval of a meta-QTL for resistance to Bgh on chromosome 5H ,. The data of QR modulation by RNAi constructs targeting the corresponding transcripts were combined with meta-data of transcript regulation, SNP or gene haplotype associations with QR to Bgh, co-localization of the RNAi target genes with QTL for resistance to Bgh, and copy number variation. As main result we present a first inventory of barley genes that are likely to play a role in broad-spectrum, durable resistance against powdery mildew fungi.
Selection of gene groups
In total 1,274 TIGS constructs were bombarded into barley epidermal cells, which corresponded to 1,144 candidate target genes due to a certain number of redundant constructs targeting the same genes. In the case of redundant constructs those with the statistically most significant effects were selected for further analysis. After the bombardment, leaf segments were inoculated either with Bgt (NHR screening) or Bgh (QR screening). The sets of TIGS constructs used in both screens overlapped partially and targeted primarily transcripts found to be upregulated during host or non-host interactions of barley epidermis with Bgh or Bgt. (Table 1). Because the selection was based on preliminary transcript-profiling data (two out of four biological replicates) a number of putatively upregulated genes finally turned out to be not regulated. This resulted - together with TIGS constructs containing unintended target-gene sequences due to PCR artifacts or other errors - in a group of genes we referred to as ‘randomly selected’ in retrospect. A third group of entry genes, which were tested only in the QR screen, consisted of members of eight multigene families that have been described in barley or other plant species to contain important PTI or ETS components (see Table 2 for details of the QR screen). The fourth group consisted of 111 target genes that could be mapped to the confidence interval of a meta-QTL for QR to Bgh on the short arm of chromosome 5H. We were previously directed towards this meta-QTL by an association-genetic approach of candidate genes that revealed co-localization of several QR-associated genes within this region at genetic distances extending beyond local linkage disequilibrium ,-. Thus, the functional architecture of the resistance QTL on 5HS might be complex with more than one causative gene acting either independently or in concerto in different genotypes of the association-mapping panel.
NHR screening by TIGS
A total of 692 TIGS constructs targeting 631 genes and corresponding to 468 (68%) pathogen-regulated transcripts were bombarded into barley epidermal cells followed by challenge inoculation with Bgt (Table 1). TIGS of 44 candidate genes caused enhanced Bgt haustorium formation, and these were repeated in a total of at least five independent experiments. As shown in Table 3 the NHR screening resulted in the identification of 10 RNAi constructs that significantly enhanced non-host susceptibility, therefore presumably targeting genes required for resistance to the non-adapted Bgt fungus (Rnr1-10, for required for non-host resistance 1-10). By applying multiple-testing correction (FDR 0.1) to the significance threshold α four constructs targeting Rnr1, Rnr3, Rnr5, and Rnr9 remained with a significant effect. However, two of the RNAi constructs without FDR-corrected significant P value targeting Rnr6 and Rnr8 gave rise to stable transgenic plants with clearly enhanced susceptibility to Bgt suggesting that the applied error 1 correction was too stringent (Douchkov et al., to be published elsewhere). The TIGS construct with the most significant effect targeting Rnr3, which encodes the syntaxin protein Hv-SNAP34, had been reported before and has been used as positive control since then. Barley cv. Maythorpe was chosen for the TIGS screening because it has a high degree of NHR to Bgt (0.1% susceptible epidermal cells) as compared to the derived, widely-used stiff-straw mutant cv. Golden Promise, which is a universally susceptible model genotype exhibiting a certain degree of non-host susceptibility to initial penetration attempts by Bgt. This might explain why the breakdown of NHR in Maythorpe caused by any TIGS construct was rather weak and did not exceed approximately 5% susceptible cells (Table 3). The largest effect was caused by silencing Hv-SNAP34. Besides Hv-SNAP34 that was also found to be required in NHR of A. thaliana to Bgh , none of the Rnr genes had been identified in NHR screens of other plant species.
QR screening by TIGS
As outline above, QR is a quantitative trait depending on favorable allelic combinations at relevant QTLs. Therefore, by selecting an appropriate combination of barley genotype and Bgh isolate resulting in a moderate level of QR it should be possible to shift the interaction in both directions, that is, enhanced resistance or susceptibility, depending on the host genes or alleles introduced or silenced. In the TIGS experimental setup, the combination of cv Maythorpe and Bgh isolate CH4.8 was well suited for this purpose because the average level of initial haustorium formation, which served as readout of QR, was 0.1 to 0.2 per penetration attempt, thereby allowing to observe shifts it in both directions. A total of 1,014 TIGS constructs targeting 914 host genes were tested (Table 1). Those increasing the percentage of susceptible (haustoria bearing) bombarded cells by a factor of at least 1.5 compared to the empty-vector control or decreasing it by a factor of 2 or more were selected for four additional, independent bombardments (Additional file 1). To address potential off-target effects of these plus the 44 RNAi constructs selected for repeated bombardments in the NHR screen, we performed off-target prediction in the set of predicted high-confidence genes of barley  by using the si-Fi software (Lück et al., in preparation; ). This resulted in the identification of 69% of the TIGS constructs with zero to only one predicted off-target (Additional file 2). Where off-targets were predicted, we performed a BlastX search among the intended target and the two most significant off-targets by emphasizing on those RNAi construct that gave rise to 2% to 10% off-target-matching siRNAs compared to the main target This focus on rather weak off-targets should have maximized the chance of finding non-paralogous genes that still might be hit by a sufficiently high number of siRNAs for silencing, As a result we identified 78% of paralogous off-targets from the same multigene-family whereas 9% off-target genes appeared to be non-paralogs. Therefore, the specificity of TIGS probably allowed at least discovery of relevant gene families in about 90% of the cases.
For statistical analysis the relative susceptibility index (SI, compared to the empty-vector control that was set to 0) were log2-transformed and compared to the log2-transformed median relative SI of the entire set of 1,014 bombarded constructs, which we expected - by its complexity - to be balanced with respect to positive or negative effects on Bgh infection, both types having been described before (see  and references cited therein). Unexpectedly the median value found was -0.36, possibly reflecting a sequence-non-specific stress effect of triggering the RNAi machinery in the bombarded cells. Out of 288 RNAi target genes silenced in repeated experiments, 44 affected the relative SI to Bgh in a statistically significant manner after multiple-testing correction (P <0.05 and FDR 0.1). Without this correction, the number of genes with significant TIGS effects increased to 90. The additional 46 genes represent weaker candidates, which was taken into account by assigning them a lower TIGS score in the final meta-dataset (see below). The average relative SI values of the repeatedly bombarded TIGS constructs were plotted against the median value of all tested 1,014 constructs. This revealed more constructs reducing susceptibility in a statistically significant manner (Figure 1B, red bars) than enhancing it significantly (green bars) suggesting that a majority of the silenced host genes might be involved in negative control of stress responses rather than defense. It remains open if some of these susceptibility-related genes are also co-opted by Bgh to facilitate fungal accommodation. To address the alternative possibility that certain constructs caused cell damage or death, which would also have prevented haustoria formation, by silencing genes with essential housekeeping function we tested all resistance-enhancing constructs in a cell-death assay, as reported -. This assay is based on reduced anthocyanin accumulation after induction of the pathway by transiently expressed C1 and b-Peru transcription factors with or without co-bombarded TIGS constructs. Cellular stress was further increased by inoculating leaf segments with Bgh 3 days after co-bombardment of anthocyanin-inducing plus TIGS constructs (Additional file 3). However, no correlation between the strength of the resistance-enhancing and cell-death inducing effect was found leading to the conclusion that many of the strongly resistance-enhancing constructs may indeed target susceptibility - or accommodation - factors of the host. The only construct clearly reducing anthocyanin accumulation in repeated experiments targeted the housekeeping TCA-cycle enzyme ATP citrate lyase (arrow), besides two lethal positive-control constructs targeting polyubiquitin genes .
We compared the percentage of TIGS constructs with a significant effect on haustoria formation between the different groups of genes and found the lowest value of 2.3% to be associated with randomly selected genes whereas it was approximately 2 to 9 times higher in the defense-related groups selected based on background knowledge (Table 2). The corresponding increase in the TIGS hit rate was statistically significant for targets encoding cellulose synthase-like proteins (Csls), RLKs, and SWEET-like sugar transporters.
Target transcript regulation
Out of the 288 target genes for QR that we tested in repeated TIGS experiments 277 were represented as oligonucleotide probes on a custom 44 K Agilent transcript-profiling array , and 90 were found to be significantly regulated by B. graminis attack in peeled epidermal samples (Figure 2). Most of these transcripts were upregulated possibly because this type of regulation had been chosen as one of the input criteria for the TIGS screenings (Table 1). A larger fraction of transcripts from genes with a significant TIGS effect were found to be rapidly induced within 6 h after inoculation compared to those not affecting relative SI upon TIGS (Figure 2, clade 1 in panels A and B). On the other hand, neither of the two groups contained members with clear differential expression in susceptible host- versus immune non-host interactions. Remarkably, almost all RNAi-target genes associated with host susceptibtility, that is, enhancing resistance upon silencing, were found to be associated with upregulated transcripts supporting our speculation that these candidate genes might indeed encode bona fide host susceptibility factors that have become co-opted by the invading fungus. The transcriptional behavior of the putative susceptibility factors is also in line with the results showing that TIGS did not often trigger cell death per se (Additional file 3), which would have offered an alternative explanation to the SI-reducing effects. It is interesting to note that the rapidly upregulated clade 1 contained significantly more resistance-associated genes causing enhanced susceptibility upon TIGS than the more slowly regulated clade 2 having its regulation peak at 24 h after inoculation. when the first haustorium was established in susceptible cells (P = 0.0237; two-tailed Fisher’s exact test).
Meta-data analysis for QR to Bgh
The significant effects obtained by TIGS revealed a number of potentially important candidate genes affecting susceptibility or QR to Bgh. Including phenotypic data from transient over-expression of 11 candidate genes, for which no TIGS data were obtained , (Additional files 1 and 4), yielded four additional candidates (U35_1790, U35_2091, U35_5202, and U35_15506) for the meta-dataset. For meta-data analysis the TIGS/over-expression data were combined with results from transcript profiling reported here or derived from ,, QTL co-localization , SNP- or gene haplotype-trait association, as well as copy number variation (CNV) from own research or from the public domain. The association-genetic data were obtained in three collections of spring and in one panel of winter barley genotypes either by re-sequencing of candidate genes  or, as shown in Additional file 5, by a genome-wide association scan using the iSelect 9 K SNP chip of llumina Co. (Sharma, 2012, PhD thesis, Martin-Luther University Halle-Wittenberg). The combination of all five criteria provided a meta-dataset allowing to assign an additive score for converging evidence (CE) to each of the TIGS target genes that was bombarded in repeated experiments (Additional file 1). Genes were assigned one scoring point each for TIGS or over-expression effects with one-tailed P value <0.05, for significant transcript regulation in powdery mildew-attacked leaves, for a map position inside the confidence intervals of QTLs for resistance to Bgh, for significant association of SNPs and/or gene haplotypes with resistance or susceptibility to Bgh, and for significant CNV. One additional score was assigned if a gene exhibited a more significant TIGS or over-expression effect (two-tailed P value <0.05), or if its transcript was significantly regulated by powdery mildew attack in leaf epidermal peels, which is the tissue directly attacked by powdery-mildew fungi. The maximum CE score of seven was not reached by any candidate gene while four genes obtained the very high score of six (Table 4). The high relevance in the barley-Bgh interaction of two of these genes encoding a chorismate synthase and a germin-like protein of subfamily 4 was confirmed by independent studies in barley and wheat ,,.
It has become increasingly clear in recent years that CNV resulting in gene sub-functionalization or enhanced transcript levels represents an important mechanism for rapid evolutionary adaptation of organisms to environmental changes. Outstanding in this respect are threats imposed by biotic stressors such as fungal pathogens because these engage (plant) hosts in co-evolutionary arm races. It is therefore not astonishing that stress-related genes, especially PRRs and NB-LRR-type resistance genes are over-represented among the ones with significant CNV profiles ,. We therefore included data from a chip-based quantitative analysis of CNV in 14 genotypes of wild and cultivated barley using the cultivar Morex as reference in the meta-dataset . Fifty-seven out of 292 genes (19.5%) tested in repeated TIGS or transient over-expression experiments exhibited significant CNV between Morex and one or several of the compared genotypes (Table 5). We tested if candidate genes for QR that are associated with high CE scores of 4 to 6 (after excluding the CE score for CNV itself) exhibit CNV more frequently than those with low scores of zero to 1. Indeed, we found a significantly higher fraction (33%) of CNV in genes with a high CE score compared to the low-scoring group (17%). Only about one-third of this difference could be explained by biased CNV occurrence in gene groups of barley belonging to different GO terms (see Additional file 6 and ) thus suggesting that the observed higher frequency of CNV among genes with high CE score might indeed be causally related to their proposed role in biotic stress responses. Next we tested a selection of genes with or without significant CNV for powdery-mildew-related differences in transcript abundance. This was done in a quantitative transcript analysis in the spring barley panel used for candidate-gene re-sequencing  and split into phenotypic bulks associated with susceptibility (bulk 1), penetration resistance (bulk 2), or late colony arrest that was frequently accompanied by darkly-pigmented spots visible to the naked eye (bulk 3). As shown in Additional file 7, 18 of 26 tested genes exhibited expression differences (p(t-test) <0.1) in Bgh-inoculated leaves between the susceptible and one or both of the resistant bulks. In some cases the absolute differences in transcript levels were small, which might reflect a rather small number of genotypes per bulk contributing to the difference. Approximately 88% and 61% of the selected genes with and without significant CNV, respectively, also exhibited differences in transcript levels between the phenotypic bulks. These results, although derived from a small sample of genes, indicate that CNV is a strong predictor of transcriptional differences upon pathogen stress, whereas absence of CNV does not exclude gene-regulatory effects that might be related to promoter polymorphisms (cis-effects) or genotype-dependent differences in upstream signaling (trans-effects). One example of causally linked CNV and differential transcript accumulation may be the dense, tandemly duplicated cluster of genes on chromosome 4HL encoding the secreted germin-like protein Hv-Ger4 with superoxide-dismutase activity. For this defense-related gene cluster we have hypothesized high transcript levels as evolutionary driving force shaping the locus ,.
We found a different distribution of high versus low CE-scoring candidate genes among functional categories (Figure 3 and Additional file 1): Genes belonging to ‘signaling’ and ‘gene expression’ were strongly eQRiched in the set of high-scoring genes whereas a higher fraction of low scoring genes belonged to the categories ‘protein’ and ‘transport’. Overall, the difference in distribution among functional categories was significant (Chi-square, two-tailed P = 0.0003). Many highly CE-scoring candidate genes in the categories ‘signaling’ or ‘gene expression’ encode for RLKs or WRKY transcription factors. Thus, host factors involved in PAMP perception, transduction of corresponding signals, and execution of transcriptional programs appear indeed to be relevant for QR of barley to powdery mildew attack, in line with an eQRichment of RLK- and WRKY factor-encoding barley genes, respectively, inside meta-QTLs for resistance against powdery mildew .
Quantitative resistance of barley to Bgh is mediated by many QTLs with small to moderate effect ,,,. This implies at least as many (non) host genes to be relevant for the trait. Despite steep progress in physical and genetic mapping as well as sequencing of the barley genome , a map-based cloning approach to all these trait-determining genes would still be very laborious and time-consuming. The same complex mode of inheritance might be true for NHR to non-adapted powdery mildew, similar to what has been described in non-adapted rust interactions , but experimental data are still missing to substantiate this speculative scenario. As an alternative approach, high-throughput reverse-genetic screenings may yield larger numbers of candidate genes provided the chosen strategy keeps the risk of false discoveries reasonably low. Reverse-genetic approaches can span a range of focusing levels from hypothesis-driven testing on one or few candidate genes up to screening the whole gene space of an organism. Here we describe an intermediate strategy by pre-selecting larger gene groups based on transcript-profiling data, gene-family membership, or QTL mapping and by entering these into a TIGS screening pipeline for NHR and QR.
As a result from the two primarily phenotype-driven TIGS screens in barley for larger groups of candidate genes we identified 10 genes designated Rnr1-10 for NHR to the wheat powdery mildew fungus Bgt and up to 90 genes for QR or host susceptibility to the barley powdery mildew fungus Bgh, with an overlap of three to four genes (depending on stringency of statistical analysis) affecting both types of interactions. The overlapping genes encode: (1) for the syntaxin Hv-SNAP34 (Rnr3) that is likely to interact with the Ror2 (for required for mlo resistance 2) syntaxin engaged in vesicle-to-target membrane fusions during pathogen-induced transcytosis; (2) for an armadillo-repeat (ARM) protein representing a partial copy of a U-box/ARM E3 ubiquitin ligase (Rnr5); (3) for the cellulose-synthase like protein Hv-CslD2 that appears to be involved in cell-wall based defense (Rnr6); and (4) for a stomatin-like protein 2 (SLP-2)-related protein (Rnr10) involved in stress-induced mitochondrial hyperfusion, touch sensation, and T-cell activation in human and animals -. All four candidates are also among the 52 high-scoring genes for QR lending further support for their importance in durable pathogen resistance (Table 4). While Rnr3 and Rnr6 represent transporting and cargo components, respectively, of vesicle-mediated targeted secretion, Rnr5 and Rnr10 are likely to be involved in intracellular regulation of protein turnover and cell death. The four overlapping genes affecting both NHR and QR point at PTI as common defense mechanism, but also at co-option in the susceptible host interaction because TIGS of Rnr5 and Rnr10 resulted in enhanced resistance. By contrast, silencing of the remaining six Rnr candidate genes did not significantly affect host QR. Because all six including Rnr8 (an LRR domain-containing RLK) are conserved between barley and wheat, we favor the idea that their apparent non-host-specific function might be related to inefficient direct or indirect neutralization by effectors of the non-adapted Bgt. Because not all genes with significant TIGS effect in the QR screening were also tested in the (earlier) NHR screening, it may well be that we currently underestimate the fraction of commonly utilized barley genes for controlling attacks by adapted- as well as non-adapted powdery mildews. Overall, we found one previously described and three novel genes of barley functioning in host- as well as non-host interactions with Bgh and Bgt, besides six genes for NHR only that might have become largely neutralized by Bgh effectors.
In the TIGS screening for candidates of QR against Bgh an unexpectedly high fraction of 64% (35/55) resulted in significantly enhanced resistance suggesting that during the susceptible host interaction many barley genes support fungal accommodation (Figure 1B). A good proportion of these susceptibility-related gene candidates were associated with upregulated corresponding transcripts in attacked leaves, which might indicate their co-option by secreted Bgh effectors at the level of promoter activity or transcript stability. Three functional categories of genes were outstanding with respect to the bias for resistance-enhancing TIGS: (1) genes involved in cellular homeostasis and cell-death control; (2) transcriptional regulators including many WRKY factors; and (3) signaling genes including mostly RLKs (Figure 1C). The P values by Fisher’s exact test for deviation of resistance- versus susceptibility-enhancing effects from the null-hypothesis of 1:1 were found to be 0.12, 0.12, and 0.035, respectively, and thus only indicative for the first two categories. Nevertheless, this result proposes the genes within the three functional categories as potential targets to the identified CSEPs of Bgh . It is interesting to note in this context that a positive mid-parent heterotic effect in a wheat F1 hybrid population for susceptibility to the wheat powdery mildew fungus was recently observed suggesting a disease-supporting effect of many genes with a dominant effect in the heterozygous state . This behavior contrasted to negative mid-parent heterotic effects of the same hybrid population for susceptibility to leaf rust and Septoria tritici blotch. Generally, TIGS was found to be a reliable tool of gene discovery because we and others could reproduce TIGS-triggered changes in Bgh interaction phenotypes in stable transgenic RNAi plants or mutants of barley or Arabidopsis ,,,-. The TIGS results from the host screening with Bgh were combined with transcript profiling, association genetic, QTL co-localization as well as CNV results leading to a meta-dataset to which we assigned CE scores ranging from zero to 6. Out of 292 repeatedly bombarded genes entered into this analysis we identified 52 candidates with a CE score of at least 4 thus representing an initial inventory of genes with a proposed role in QR (Table 4). We tentatively mapped these onto cellular processes and compartments in a powdery mildew-attacked barley epidermal cell in order to get an impression of important defense- or susceptibility-related pathways (Figure 4). The 52 high-scoring candidate genes were also searched for available literature information with respect to plant-pathogen interactions (Additional file 8). This revealed 37 genes that have either been directly described in the interaction of barley or other plant species with microbial pathogens, or represent homologs and gene-family members of such functionally-characterized genes. A good proportion of these (78%) include known components of PTI or ETS pathways such as pattern recognition receptors (PRRs) encoded by RLK genes, WRKY transcription factors, the cell-death control factors Hv-Mlo, Hv-Lsd1a, and Bax-inhibitor 1, a member of the SWEET-family of sugar transporters, enzymes of the shikimate and phenylpropanoid plus oxidative pathways leading, for example, to cell-wall lignification ,,-. It therefore appears that in the Triticeae crop plant barley QR indeed reflects the difference between PTI and ETS, as suggested in model plant systems  (Figure 4 and Additional file 8). Interestingly, genes involved in early steps of PTI (Figure 4, left side of the model) appear to be often involved in mediating susceptibility (TIGS results, highlighted by green gene numbers) whereas genes positioned further downstream of the pathway (right side of the model) are more likely to encode defense-related proteins (TIGS data, highlighted by red gene numbers). It has to be considered, though, that a large proportion of genes entered the TIGS QR screen because of their transcript-up-regulation upon pathogen attack or gene-family membership, which might have biased the distribution of functional categories among low and high CE-scoring genes. However, several multigene families including ABC transporters or class III peroxidases with clearly identified, important roles in QR or NHR did not show an increased frequency of TIGS effects among their members (Table 2). This finding argues against a strong family bias introduced by the gene-selection procedure and rather suggests that functional diversification among members of eukaryotic multigene families tends to be too complex for family-wide predictions about their involvement in a particular biological process. The non-predictability of physiological gene function is probably further pronounced in di- or tritrophic interactions between host plants and attacking parasitic organisms because these always reflect the current status of a highly dynamic and interaction-specific co-evolutionary arms race.
Among the high-scoring genes 15 have no record in the literature for being involved in plant defense or pathogen accommodation, thereby offering extensions to existing models of plant-pathogen interactions, at least as far as biotrophic fungal pathogens are concerned (Additional file 8). Some of those appear especially interesting: First, an aintegumenta-like AP2/EREBP transcription factor (QR. 5, U35_11820) was found to be upregulated in Bgh-attacked epidermis and to induce strong resistance when silenced. This type of transcription factor is associated with dividing meristematic tissue and might be a first lead into effector-triggered endo-reduplication in barley epidermis as recently suggested in a susceptible A. thaliana-powdery mildew interaction ,. Second, a spliceosome component encoded by the PRP38-like gene U35_2705 might point at an important role of pre-mRNA processing to support fungal growth . Third, an ARM-repeat protein encoded by U35_3071 (Rnr5) was identified as partial gene duplicate of an Os-PUB15-related E3-ubiquitin ligase. Deeper investigation suggested that this partial protein might represent a decoy for the E3-ligase targeted by a Bgh-encoded effector (Rajaraman and Schweizer, unpublished). Fourth, silencing of the transcriptionally upregulated Rnr6 gene encoding Hv-CslD2, a member of the cellulose synthase-like (Csl) protein family, caused attenuation of QR and NHR (Tables 3 and 4). We therefore tested all available Csl unigenes of barley for a potential role in QR. This revealed Hv-CslA11 as additional candidate with a high CE score. Several members of the Csl family are known to synthesize non-cellulosic cell-wall carbohydrates such as mixed-linked β-1-3:1-4-glucans during the build-up or modification of secondary cell walls while the enzymatic activity of others is still unclear . Members of the CslA- and CslD-clades of the Csl protein family were shown or proposed to act as mannan synthases putting forward mannose-containing polysaccharides as potentially important for penetration resistance -. However, defense-related functions of mannan(s) in plants are currently not known . A few more high-scoring candidate genes are discussed in Additional file 6.
This study represents an initial, functional inventory of approximately 3% of the barley transcriptome for a role in NHR and QR against the powdery mildew pathogen. The discovered candidate genes support the idea that broad-spectrum, quantitative, and durable disease resistance in barley reflects to a large extent the difference of PAMP-triggered immunity minus effector-mediated host susceptibility. By extending the approach we expect more genes to be discovered in the future, which will require strong priorization for their labor-intensive validation in barley or related Triticeae crop plants. Meanwhile the ongoing in-depth analysis of the function of prioritized candidates will provide proof of concept for the approach of convergent evidence for the discovery of genes that are relevant to more durable forms of polygenically inherited pathogen resistance in barley.
Materials and methods
Plant and fungal material
For the TIGS screening 7-day-old seedlings of a spring barley cv. Maythorpe were used . This genotype, from which the universal susceptible cv. Golden Promise was derived by γ-ray mutagenesis, proven to be well-suited for TIGS screenings in host as well as non-host interactions with powdery mildew fungi because it was fully resistant to the wheat powdery mildew Bgt while exhibiting a moderate level of QR to Bgh, thereby allowing to detect resistance- as well as susceptibility-enhancing TIGS effects in the host interaction. The mutagenesis leading to Golden Promise appeared to have affected multiple traits besides the initially targeted stiff-straw growth habit including salt tolerance, enhanced susceptibility to powdery mildew and high efficiency of Agrobacterium-mediated transformation . Seedlings were grown in a plant incubator (Sanyo/Panasonic, address) at 20°C constant temperature, 50% rel. humidity and 16 h illumination (intensity level 5) by fluorescent tubes (OSRAM L36W/840).
For genome-wide association mapping of SNP markers with Bgh resistance, single-seed-derived 224 Genobar spring barley collection, plus 282 spring and 112 winter barley genotypes were used as described in  and . For the detached leaf assay screening of resistance to Bgh, plants were grown in trays using disease free standard greenhouse conditions at 17°C to 20°C under long day conditions (16 h).
For the RT-qPCR-based transcript analysis of phenotypic bulks of barley differing in response to Bgh, single seed-derived lines of the following accessions were used: BCC1404, BCC1412, BCC1420, BCC1430, BCC1431, BCC1450, BCC1452, BCC1468, BCC1488, BCC1498, BCC423, BCC745, BCC888, BCC893, BCC903, HOR2800, HOR3941, HOR4060, BCC852, BCC1376, (susceptible bulk); HOR261, HOR728, HOR804, HOR842, HOR1036, HOR1457, HOR1506, HOR2543, HOR2591, HOR2932, HOR3270, HOR3271, HOR3537, HOR3726, HOR3988, HOR4021, HOR4408 (penetration resistant bulk); HOR214, HOR262, HOR303, HOR683, HOR736, HOR795, HOR800, HOR1159, HOR1379, HOR1468, HOR1647, HOR1873, HOR2573, HOR3041, HOR3075, HOR3983, HOR3984, HOR4400, (late resistant bulk). These plants were grown in the greenhouse with additional light (16 h) provided by sodium halogen lamps.
Plants of cv Vada used for the transcript profiling experiments were grown in a climate chamber at 19°C with 65% relative humidity (RH) during the night and 23°C with 45% RH during the 16 h photoperiod.
Inoculation experiments for TIGS or transcript profiling were performed using Swiss field isolate CH4.8 of Bgh or Swiss field isolate FAL 92315 of Bgt. For the detached leaf assay of the genome-wide association scan, the polyvirulent German Bgh isolates D12-12 and 78P were used as described . Field data of Bgh resistance were derived from natural infection at IPK in 2009 and 2010.
Target genes for TIGS were selected and analyzed based on sequence-contig information of the HarvEST database, barley 1.83 assembly #35 . Putative off-target effects of TIGS constructs were predicted by using the si-Fi software (labtools.ipk-gatersleben.de) for finding all sequence matching putative siRNAfollowed by the application of an algorithm for siRNA guide strand selection as described . TIGS Constructs were generated and transferred into barley leaf epidermal cells by particle bombardment followed by inoculation with Bgh and Bgt 3 and 4 days after bombardment, respectively, as described . In the NHR screening, the number of susceptible cells carrying at least one Bgt haustorium was counted 48 h after inoculation. Constructs inducing at least four susceptible cells per bombardment were used for repeated bombardments resulting in a total of five independent biological replicates. Final results per TIGS construct were compared to the empty-vector control pIPKTA30 by using the Mann-Whitney non-parametric test. In the QR screening, transformed GUS-stained epidermal cells as well as haustoria-containing transformed (susceptible) cells were counted 48 h after inoculation. The susceptibility index (SI) was calculated by dividing the number of susceptible cells by the total number of transformed cells, followed by normalization to SI of the empty-vector control pIPKTA30N (rel. SI). Values of rel. SI were log(2)-transformed in order to normalize their distribution for statistical analysis by a one-sample t-test. This test was performed against the hypothetical relative susceptibility-index value ‘-0.355’ corresponding to the observed median of more than 1,000 RNAi constructs. The deviation from the control value of the empty pIPKTA30N vector (set to ‘0’) may reflect a weak and insert-non-specific side effect of triggering the cellular RNAi machinery.
For the examination of cell death-inducing TIGS effects we performed particle co-bombardment of RNAi constructs pIPKTA30N_targetX (7 μg DNA/shot) with the B-Peru/C1-expression plasmid pBC17 (7 μg DNA/shot) triggering the anthocyanin biosynthetic pathway  and pUbiGUS (7 μg DNA/shot), as described by Dong et al. . In this assay, cell death is reflected by a reduction of the GUS-normalized number of anthocyanin-accumulating cells.
Seven-day-old barley plants of cv. Vada were inoculated with Bgh or Bgt, and the abaxial epidermis of inoculated primary leaves or from non-inoculated control leaves was peeled at 6 to 74 h after inoculation, as described . Total, quality-controlled RNA was hybridized to a 44 K Agilent oligonucleotide array as described . Single-channel array processing was utilized followed by data normalization with default parameters, and significant transcript-regulation events were determined by using GeneSpring GX (v11.5.1) software (Agilent technologies Inc). Transcripts were assumed to be significantly regulated if P values corrected for false-positive rate (FDR, Benjamini-Hochberg method) were less than 0.05 and if regulation factors between inoculated and corresponding control samples harvested in parallel exceeded 2.0. All quantile-normalized signal intensities of the analyzed candidate genes are shown in Additional file 9, and the raw data from the corresponding array slides were deposited at ArrayExpress (Accession E-MTAB-2916).
The RT-qPCR-based transcript analysis of phenotypic bulks of barley accessions differing in response to Bgh was performed by using Power SYBR® Green PCR Master mix kit (Applied Biosystems, Foster City, CA, USA) on an ABI 7900HT Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA). Data were analyzed using the gene-testing standard-curve approach. Seven-day-old, greenhouse-grown (with additional light from sodium halogen lamps) seedlings were inoculated with Bgh at a density of approximately 10 to 30 conidia mm-2. Twelve hours after inoculation, total RNA was isolated . RNA was treated with DNA-free™ Kit (Ambion, Austin, TX, USA) before cDNA synthesis. One microgram total RNA was used for cDNA synthesis using the iScript cDNA Synthesis Kit (Bio-Rad Co., Munich, Germany).
Genotypic and phenotypic data of candidate-gene based association mapping by allele re-sequencing were derived from a previous study . As significance threshold for SNP-trait or haplotype-trait association we selected the (-log10)P value of 3.0 in a general linear model including marker + trait + population structure + row number (covariant). Associations were calculated using the TASSEL software package.
Phenotypic data for genome-wide association mapping were derived from detached leaf assays of seedlings as described by Altpeter et al.  and from field data of spontaneous Bgh infection at IPK Gatersleben in 2009 and 2010. For field evaluation the percentage of leaf infection on a plot basis was scored and the Restriction Estimate of Maximum Likelihood (REML) implemented in Genstat software 14th edition (VSN International, Hemel Hempstead, UK) was used to generate means over the years. For the detached leaf assay screen, two polyvirulent Bgh isolates ‘D12-12’ and ‘78P’ were used, as described in . Barley genotypes were only scored as resistant if the infection did not exceed rating class 1  with either Bgh isolate. Genotypic data for genome-wide association mapping were derived from the Illumina 9 K SNP chip of barley . SNP marker-trait associations were calculated in TASSEL 2.1  and a mixed-linear model using kinship from random markers was used to control population structure. Marker-trait associations were considered as significant if the (-log10)P value was larger than 3.0 per single SNP in the candidate gene or if it exceeded 2.0 in the candidate gene plus in at least two immediately adjacent genes (sliding window approach). Per gene, the most significant association derived from field- or detached leaf-assay data was used for assigning ‘0’ or ‘1’ AM scores (Additional file 5).
Gene functional categories
Each gene was manually assigned one out of 10 broad functional categories because the use of existing binning systems such as Gene Ontology or MapMan resulted in a large proportion of non-assigned genes that - by hand-curated BlastX analysis - could often be assigned one of the widely defined functional categories.
Primary transcript-regulation data of this study, or obtained from the public domain, were used , (PlexDB, ). Transcripts were assumed to be significantly regulated if normalized signal intensities (inoculated versus control samples) exceeded 2.0 in at least one of the analyzed time points after inoculation, and if the Benjamini-Hochberg-corrected P value for the null hypothesis was lower than 0.05. Regulation events in epidermal peels and in leaf samples were assigned 2 and 1 CE scores, respectively. If a gene was significantly regulated both in epidermal and leaf samples a score of 2 was assigned.
Candidate genes were tested for co-localization with meta-QTL for powdery-mildew resistance as described . Briefly, co-localization was assumed if the gene was positioned between the outmost flanking markers of meta-QTL for resistance to Bgh consisting of ≥3 overlapping QTL in the consensus linkage map ‘Marcel et al., integrated, 2009’ (deposited in GrainGenes 2.0 database).
Copy-number variation of candidate genes was tested in a panel of 14 barley genotypes and by using a custom Comparative Genomic Hybridization array designed by Roche NimbleGen (Roche NimbleGen, Inc., Madison, WI, USA) that used 2.2 M contigs from a whole genome shotgun (WGS) assembly of barley cv. Morex . For CNV assessment the expectation maximization algorithm was used to estimate the mixing proportion, mean, and variance associated with two predicted signal sub-distributions found within the tested genotype vs. Morex fragments. When the log2 signal ratio was positive, the variant was defined as ‘UpCNV’, while it was classified as ‘DownCNV/PAV’ when the ratio was negative ().
Jones JDG, Dangl JL: The plant immune system. Nature. 2006, 444: 323-329. 10.1038/nature05286.
Boller T, Felix G: A renaissance of elicitors: perception of microbe-associated molecular patterns and danger signals by pattern-recognition receptors. Annu Rev Plant Biol. 2009, 60: 379-406. 10.1146/annurev.arplant.57.032905.105346.
Boehm H, Albert I, Fan L, Reinhard A, Nuernberger T: Immune receptor complexes at the plant cell surface. Curr Opin Plant Biol. 2014, 20: 47-54. 10.1016/j.pbi.2014.04.007.
Fan J, Doerner P: Genetic and molecular basis of nonhost disease resistance: complex, yes; silver bullet, no. Curr Opin Plant Biol. 2012, 15: 400-406. 10.1016/j.pbi.2012.03.001.
Schulze-Lefert P, Panstruga R: A molecular evolutionary concept connecting nonhost resistance, pathogen host range, and pathogen speciation. Trends Plant Sci. 2011, 16: 117-125. 10.1016/j.tplants.2011.01.001.
Tosa Y: A model for the evolution of formae speciales and races. Phytopathology. 1992, 82: 728-730.
Niks RE, Marcel TC: Nonhost and basal resistance: how to explain specificity?. New Phytol. 2009, 182: 817-828. 10.1111/j.1469-8137.2009.02849.x.
Kou YJ, Wang SP: Broad-spectrum and durability: understanding of quantitative disease resistance. Curr Opin Plant Biol. 2010, 13: 181-185. 10.1016/j.pbi.2009.12.010.
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 Interact. 2011, 24: 1492-1501. 10.1094/MPMI-05-11-0107.
St Clair DA: Quantitative disease resistance and quantitative resistance loci in breeding. Annu Rev Phytopathol. 2010, 48: 247-268. 10.1146/annurev-phyto-080508-081904.
Koeck M, Hardham AR, Dodds PN: The role of effectors of biotrophic and hemibiotrophic fungi in infection. Cell Microbiol. 2011, 13: 1849-1857. 10.1111/j.1462-5822.2011.01665.x.
Lauge R, Joosten MH, Haanstra JP, Goodwin PH, Lindhout P, De Wit PJ: Successful search for a resistance gene in tomato targeted against a virulence factor of a fungal pathogen. Proc Natl Acad Sci U S A. 1998, 95: 9014-9018. 10.1073/pnas.95.15.9014.
Poland JA, Balint-Kurti PJ, Wisser RJ, Pratt RC, Nelson RJ: Shades of gray: the world of quantitative disease resistance. Trends Plant Sci. 2009, 14: 21-29. 10.1016/j.tplants.2008.10.006.
Dowkiw A, Bastien C: Presence of defeated qualitative resistance genes frequently has major impact on quantitative resistance to Melampsora larici-populina leaf rust in P.xinteramericana hybrid poplars. Tree Genetics Genomes. 2007, 3: 261-274. 10.1007/s11295-006-0062-0.
Glawe DA: The powdery mildews: a review of the world’s most familiar (yet poorly known) plant pathogens. Annu Rev Phytopathol. 2008, 46: 27-51. 10.1146/annurev.phyto.46.081407.104740.
Collins NC, Sadanandom A, Schulze-Lefert P: Genes and molecular mechanisms controlling powdery mildew resistance in barley. The Powdery Mildews. Edited by: Bélanger RR, Bushnell WR, Dik AJ, Carver TLW. 2002, APS Press, St. Paul, MN, 134-145.
Huckelhoven R: Cell wall - associated mechanisms of disease resistance and susceptibility. Annu Rev Phytopathol. 2007, 45: 101-127. 10.1146/annurev.phyto.45.062806.094325.
Wise RP, Lauter N, Szabo LJ, Schweizer P: Genomics of biotic interactions in the triticeae. Genetics and Genomics of the Triticeae. Edited by: Muehlbauer GJ, Feuillet C. 2009, Springer, New York, 559-589. 10.1007/978-0-387-77489-3_19. 7
Altpeter F, Varshney A, Abderhalden O, Douchkov D, Sautter C, Kumlehn J, Dudler R, Schweizer P: Stable expression of a defense-related gene in wheat epidermis under transcriptional control of a novel promoter confers pathogen resistance. Plant Mol Biol. 2005, 57: 271-283. 10.1007/s11103-004-7564-7.
Ihlow A, Schweizer P, Seiffert U: A high-throughput screening system for barley/powdery mildew interactions based on automated analysis of light micrographs. BMC Plant Biol. 2008, 8: 6-10.1186/1471-2229-8-6.
Douchkov D, Lück S, Baum T, Seiffert U, Schweizer P: Microphenomics for interactions of barley with fungal pathogens. Advances in Genomics of Plant Genetic Resources. Edited by: Tuberosa R, Graner A, Frison E. 2014, Springer, Heidelberg, 123-148. 10.1007/978-94-007-7575-6_5. 2
Huckelhoven R, Dechert C, Kogel KH: Overexpression of barley BAX inhibitor 1 induces breakdown of mlo-mediated penetration resistance to Blumeria graminis. Proc Natl Acad Sci U S A. 2003, 100: 5555-5560. 10.1073/pnas.0931464100.
Douchkov D, Nowara D, Zierold U, Schweizer P: A high-throughput gene-silencing system for the functional assessment of defense-related genes in barley epidermal cells. Mol Plant-Microbe Interact. 2005, 18: 755-761. 10.1094/MPMI-18-0755.
Elliott C, Zhou FS, 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 Interact. 2002, 15: 1069-1077. 10.1094/MPMI.2002.15.10.1069.
Zellerhoff N, Himmelbach A, Dong WB, Bieri S, Schaffrath U, Schweizer P: Nonhost resistance of barley to different fungal pathogens is associated with largely distinct, quantitative transcriptional responses. Plant Physiol. 2010, 152: 2053-2066. 10.1104/pp.109.151829.
Zierold U, Scholz U, Schweizer P: Transcriptome analysis of mlo-mediated resistance in the epidermis of barley. Mol Plant Pathol. 2005, 6: 139-151. 10.1111/j.1364-3703.2005.00271.x.
Spies A, Korzun L, Bayles R, Rajaraman J, Himmelbach A, Hedley PE, Schweizer P: Allele mining in barley genetic resources reveals genes of race-nonspecific powdery mildew resistance. Front Plant Sci Plant-Microbe Interact. 2012, 2: 113-
Aghnoum R, Marcel TC, Johrde A, Pecchioni N, Schweizer P, Niks RE: Basal host resistance of barley to powdery mildew: connecting quantitative trait loci and candidate genes. Mol Plant-Microbe Interact. 2010, 23: 91-102. 10.1094/MPMI-23-1-0091.
Douchkov D, Johrde A, Nowara D, Himmelbach A, Lueck S, Niks R, Schweizer P: Convergent evidence for a role of WIR1 proteins during the interaction of barley with the powdery mildew fungus Blumeria graminis. J Plant Physiol. 2011, 168: 20-29. 10.1016/j.jplph.2010.07.004.
Collins NC, Thordal-Christensen H, Lipka V, Bau S, Kombrink E, Qiu JL, Huckelhoven R, Stein M, Freialdenhoven A, Somerville SC, Schulze-Lefert P: SNARE-protein-mediated disease resistance at the plant cell wall. Nature. 2003, 425: 973-977. 10.1038/nature02076.
Mayer K, Consortium IBS: A physical, genetic and functional sequence assembly of the barley genome. Nature. 2012, 491: 711-716.
si - Fi is a Software for RNAi (RNA interference) off-target prediction. [http://labtools.ipk-gatersleben.de/]
Dong WB, Nowara D, Schweizer P: Protein polyubiquitination plays a role in basal host resistance of barley. Plant Cell. 2006, 18: 3321-3331. 10.1105/tpc.106.046326.
Eichmann R, Bischof M, Weis C, Shaw J, Lacomme C, Schweizer P, Duchkov D, Hensel G, Kumlehn J, Huckelhoven R: BAX INHIBITOR-1 is required for full susceptibility of barley to powdery mildew. Mol Plant-Microbe Interact. 2010, 23: 1217-1227. 10.1094/MPMI-23-9-1217.
Pliego C, Nowara D, Bonciani G, Gheorghe DM, Xu R, Surana P, Whigham E, Nettleton D, Bogdanove AJ, Wise RP, Schweizer P, Bindschedler LV, Spanu PD: Host-induced gene silencing in barley powdery mildew reveals a class of ribonuclease-like effectors. Mol Plant-Microbe Interact. 2013, 26: 633-642. 10.1094/MPMI-01-13-0005-R.
Chen XW, Hedley PE, Morris J, Liu H, Niks RE, Waugh R: Combining genetical genomics and bulked segregant analysis-based differential expression: an approach to gene localization. Theor Appl Genet. 2011, 122: 1375-1383. 10.1007/s00122-011-1538-3.
Zimmermann G, Baumlein H, Mock HP, Himmelbach A, Schweizer P: The multigene family encoding germin-like proteins of barley. Regulation and function in basal host resistance. Plant Physiol. 2006, 142: 181-192. 10.1104/pp.106.083824.
Johrde A, Schweizer P: A class III peroxidase specifically expressed in pathogen-attacked barley epidermis contributes to basal resistance. Mol Plant Pathol. 2008, 9: 687-696. 10.1111/j.1364-3703.2008.00494.x.
Moscou MJ, Lauter N, Caldo RA, Nettleton D, Wise RP: Quantitative and temporal definition of the Mla transcriptional regulon during barley-powdery mildew interactions. Mol Plant-Microbe Interact. 2011, 24: 694-705. 10.1094/MPMI-09-10-0211.
Hu PS, Meng Y, Wise RP: Functional contribution of chorismate synthase, anthranilate synthase, and chorismate mutase to penetration resistance in barley-powdery mildew interactions. Mol Plant-Microbe Interact. 2009, 22: 311-320. 10.1094/MPMI-22-3-0311.
Christensen AB, Thordal-Christensen H, Zimmermann G, Gjetting T, Lyngkjaer MF, Dudler R, Schweizer P: The germinlike protein GLP4 exhibits superoxide dismutase activity and is an important component of quantitative resistance in wheat and barley. Mol Plant Microbe Interact. 2004, 17: 109-117. 10.1094/MPMI.2004.17.1.109.
Zmienko A, Samelak A, Kozlowski P, Figlerowicz M: Copy number polymorphism in plant genomes. Theor Appl Genet. 2014, 127: 1-18. 10.1007/s00122-013-2177-7.
Muñoz‐Amatriaín M, Eichten SR, Wicker T, Richmond TA, Mascher M, Steuernagel B, Scholz U, Ariyadasa R, Spannagl M, Nussbaumer T, Mayer KFX, Taudien S, Platzer M, Jeddeloh JA, Springer NM, Muehlbauer GJ, Stein N: Distribution, functional impact, and origin mechanisms of copy number variation in the barley genome. Genome Biol. 2013, 14: R58-10.1186/gb-2013-14-6-r58.
Druka A, Kudrna D, Kannangara CG, Von Wettstein D, Kleinhofs A: Physical and genetic mapping of barley (Hordeum vulgare) germin-like cDNAs. Proc Natl Acad Sci U S A. 2002, 99: 850-855. 10.1073/pnas.022627999.
Himmelbach A, Liu L, Zierold U, Altschmied L, Maucher H, Beier F, Muller D, Hensel G, Heise A, Schutzendubel A, Kumlehn J, Schweizer P: Promoters of the barley germin-like GER4 gene cluster enable strong transgene expression in response to pathogen attack. Plant Cell. 2010, 22: 937-952. 10.1105/tpc.109.067934.
Miedaner T, Korzun V: Marker-assisted selection for disease resistance in wheat and barley breeding. Phytopathology. 2012, 102: 560-566. 10.1094/PHYTO-05-11-0157.
Aghnoum R, Niks RE: Transgressive segregation for very low and high levels of basal resistance to powdery mildew in barley. J Plant Physiol. 2011, 168: 45-50. 10.1016/j.jplph.2010.09.005.
Jafary H, Albertazzi G, Marcel TC, Niks RE: High diversity of genes for nonhost resistance of barley to heterologous rust fungi. Genetics. 2008, 178: 2327-2339. 10.1534/genetics.107.077552.
Huang MX, Gu GQ, Ferguson EL, Chalfie M: A stomatin-like protein necessary for mechanosensation in c-elegans. Nature. 1995, 378: 292-295. 10.1038/378292a0.
Tondera D, Grandemange S, Jourdain A, Karbowski M, Mattenberger Y, Herzig S, Da Cruz S, Clerc P, Raschke I, Merkwirth C, Ehses S, Krause F, Chan DC, Alexander C, Bauer C, Youle R, Langer T, Martinou JC: SLP-2 is required for stress-induced mitochondrial hyperfusion. EMBO J. 2009, 28: 1589-1600. 10.1038/emboj.2009.89.
Christie DA, Mitsopoulos P, Blagih J, Dunn SD, St-Pierre J, Jones RG, Hatch GM, Madrenas J: Stomatin-like protein 2 deficiency in T cells is associated with altered mitochondrial respiration and defective CD4+ T cell responses. J Immunol. 2012, 189: 4349-4360. 10.4049/jimmunol.1103829.
Pedersen C, Ver Loren van Themaat E, McGuffin LJ, Abbott JC, Burgis TA, Barton G, Bindschedler L, Lu X, Maekawa T, Weßling R, Cramer R, Thordal-Christensen H, Panstruga R, Spanu PD: Structure and evolution of barley powdery mildew effector candidates. BMC Genomics. 2012, 13: 649-10.1186/1471-2164-13-694.
Longin CFH, Gowda M, Muhleisen J, Ebmeyer E, Kazman E, Schachschneider R, Schacht J, Kirchhoff M, Zhao YS, Reif JC: Hybrid wheat: quantitative genetic parameters and consequences for the design of breeding programs. Theor Appl Genet. 2013, 126: 2791-2801. 10.1007/s00122-013-2172-z.
Schweizer P, Pokorny J, Schulze-Lefert P, Dudler R: Double-stranded RNA interferes with gene function at the single-cell level in cereals. Plant J. 2000, 24: 895-903. 10.1046/j.1365-313x.2000.00941.x.
Schultheiss H, Dechert C, Kogel KH, Huckelhoven R: A small GTP-binding host protein is required for entry of powdery mildew fungus into epidermal cells of barley. Plant Physiol. 2002, 128: 1447-1454. 10.1104/pp.010805.
Hoefle C, Huesmann C, Schultheiss H, Bornke F, Hensel G, Kumlehn J, Huckelhoven R: A barley ROP GTPase ACTIVATING PROTEIN associates with microtubules and regulates entry of the barley powdery mildew fungus into leaf epidermal cells. Plant Cell. 2011, 23: 2422-2439. 10.1105/tpc.110.082131.
Bueschges R, Hollricher K, Panstruga R, Simons G, Wolter M, Frijters A, Van DR, Van DLT, Diergarde P, Groenendijk J, Topsch 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.
Dietrich RA, Richberg MH, Schmidt R, Dean C, Dangl JL: A novel zinc finger protein is encoded by the Arabidopsis Lsd1 gene and functions as a negative regulator of plant cell death. Cell. 1997, 88: 685-694. 10.1016/S0092-8674(00)81911-X.
Meng Y, Wise RP: HvWRKY10, HvWRKY19, and HvWRKY28 regulate Mla-triggered immunity and basal defense to barley powdery mildew. Mol Plant-Microbe Interact. 2012, 25: 1492-1505. 10.1094/MPMI-04-12-0082-R.
Hahlbrock K, Bednarek P, Ciolkowski I, Hamberger B, Heise A, Liedgens H, Logemann E, Nurnberger T, Schmelzer E, Somssich IE, Tan JW: Non-self recognition, transcriptional reprogramming, and secondary metabolite accumulation during plant/pathogen interactions. Proc Natl Acad Sci U S A. 2003, 100: 14569-14576. 10.1073/pnas.0831246100.
Vellosillo T, Vicente J, Kulasekaran S, Hamberg M, Castresana C: Emerging complexity in reactive oxygen species production and signaling during the response of plants to pathogens. Plant Physiol. 2010, 154: 444-448. 10.1104/pp.110.161273.
Chen LQ, Hou BH, Lalonde S, Takanaga H, Hartung ML, Qu XQ, Guo WJ, Kim JG, Underwood W, Chaudhuri B, Chermak D, Antony G, White FF, Somerville SC, Mudgett MB, Frommer WB: Sugar transporters for intercellular exchange and nutrition of pathogens. Nature. 2010, 468: 527-532. 10.1038/nature09606.
Truman W, de Zabala MT, Grant M: Type III effectors orchestrate a complex interplay between transcriptional networks to modify basal defence responses during pathogenesis and resistance. Plant J. 2006, 46: 14-33. 10.1111/j.1365-313X.2006.02672.x.
Nole-Wilson S, Tranby TL, Krizek BA: AINTEGUMENTA-like (AIL) genes are expressed in young tissues and may specify meristematic or division-competent states. Plant Mol Biol. 2005, 57: 613-628. 10.1007/s11103-005-0955-6.
Chandran D, Inada N, Hather G, Kleindt CK, Wildermuth MC: Laser microdissection of Arabidopsis cells at the powdery mildew infection site reveals site-specific processes and regulators. Proc Natl Acad Sci U S A. 2010, 107: 460-465. 10.1073/pnas.0912492107.
Blanton S, Srinivasan A, Rymond BC: PRP38 encodes a yeast protein required for pre-messenger-Rna splicing and maintenance of stable U6 small nuclear-Rna levels. Mol Cell Biol. 1992, 12: 3939-3947.
Wang LQ, Guo K, Li Y, Tu YY, Hu HZ, Wang BR, Cui XC, Peng LC: Expression profiling and integrative analysis of the CESA/CSL superfamily in rice. BMC Plant Biol. 2010, 10: 282-10.1186/1471-2229-10-282.
Dhugga KS, Barreiro R, Whitten B, Stecca K, Hazebroek J, Randhawa GS, Dolan M, Kinney AJ, Tomes D, Nichols S, Anderson P: Guar seed beta-mannan synthase is a member of the cellulose synthase super gene family. Science. 2004, 303: 363-366. 10.1126/science.1090908.
Liepman AH, Wilkerson CG, Keegstra K: Expression of cellulose synthase-like (Csl) genes in insect cells reveals that CslA family members encode mannan synthases. Proc Natl Acad Sci U S A. 2005, 102: 2221-2226. 10.1073/pnas.0409179102.
Doblin MS, De Melis L, Newbigin E, Bacic A, Read SM: Pollen tubes of Nicotiana alata express two genes from different beta-glucan synthase families. Plant Physiol. 2001, 125: 2040-2052. 10.1104/pp.125.4.2040.
Hrmova M, Burton RA, Biely P, Lahnstein J, Fincher GB: Hydrolysis of (1,4)-beta-D-mannans in barley (Hordeum vulgare L.) is mediated by the concerted action of (1,4)-beta-D-mannan endohydrolase and beta-D-mannosidase. Biochem J. 2006, 399: 77-90. 10.1042/BJ20060170.
Forster BP: Mutation genetics of salt tolerance in barley: an assessment of Golden Promise and other semi-dwarf mutants. Euphytica. 2001, 120: 317-328. 10.1023/A:1017592618298.
Pasam RK, Sharma R, Malosetti M, van Eeuwijk FA, Haseneyer G, Kilian B, Graner A: Genome-wide association studies for agronomical traits in a world wide spring barley collection. BMC Plant Biol. 2012, 12: 16-10.1186/1471-2229-12-16.
Tondelli A, Xu X, Moragues M, Sharma R, Schnaithmann F, Ingvardsen C, Manninen O, Comadran J, Russell J, Waugh R, Schulman AH, Pillen K, Rasmussen SK, Kilian B, Cattivelli L, Thomas WTB, Flavell AJ: Structural and temporal variation in genetic diversity of european spring two-row barley cultivars and association mapping of quantitative traits. Plant Genome. 2013, 6: 2-10.3835/plantgenome2013.03.0007.
HarvEST, an EST database for crop plant species. [http://www.harvest-web.org]
Khvorova A, Reynolds A, Jayasena SD: Functional siRNAs and miRNAs exhibit strand bias. Cell. 2007, 131: 41-49.
Chomczynski P, Sacchi N: Single-step method of Rna isolation by acid guanidinium thiocyanate phenol chloroform extraction. Anal Biochem. 1987, 162: 156-159. 10.1016/0003-2697(87)90021-2.
Schweizer P, Vallelianbindschedler L, Mosinger E: Heat-induced resistance in barley to the powdery mildew fungus erysiphe-graminis F-Sp hordei. Physiol Mol Plant Pathol. 1995, 47: 51-66. 10.1006/pmpp.1995.1042.
Comadran J, Kilian B, Russell J, Ramsay L, Stein N, Ganal M, Shaw P, Bayer M, Thomas W, Marshall D, Hedley P, Tondelli A, Pecchioni N, Francia E, Korzun V, Walther A, Waugh R: Natural variation in a homolog of Antirrhinum CENTRORADIALIS contributed to spring growth habit and environmental adaptation in cultivated barley. Nat Genet. 2012, 44: 1388-1392. 10.1038/ng.2447.
TASSEL: Software for association mapping of complex traits in diverse samples. [http://www.maizegenetics.net]
MIAME/Plant Compliant Gene Expression Resources for Plants and Plant Pathogens. [http://www.plexdb.org]
We would like to thank Gabi Brantin and Manuela Knauft for excellent technical assistance, and Pete Hedley (James Hutton Institute) for performing Agilent array hybridazations. This work was supported by IPK (to PS), by BASF Plant Science GmbH (to PS), by DFG (project Nr. SCHW848/2-1 to PS), by the German Ministry for Education and Research (BMBF, projects GABI-non-host and GABI-phenome to PS), by ERA-net PG project EXBARDIV (to AG), and by EU-FP6 project Bioexploit (to PS).
The authors declare that they have no competing interests.
PS designed the research, analyzed meta-data, and drafted the manuscript. BK and NS designed the research. DD, SL, and DN performed TIGS experiments and analyzed data. AJ, AH, and JR performed transcript profiling experiments and analyzed the data. RS performed the association mapping experiments and analyzed the data. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: List of TIGS target genes selected for repeated bombardments with their corresponding meta-dataset leading to the CE score. (XLSX 97 KB)
Additional file 4: Effect of transient over-expression of candidate genes exhibiting expression differences between susceptible and mlo -resistant near-isogenic barley lines. (DOCX 14 KB)
Additional file 6: More detailed discussion of several candidate genes with high CE scores according to Table 4, plus additional references. (DOCX 55 KB)
Additional file 7: Summary of results from transcript quantification in phenotypic bulks of barley genotypes with different responses to Bgh attack. (XLSX 15 KB)
Additional file 8: Literature survey of 42 high-scoring candidate genes with respect to plant-pathogen interactions. (DOCX 25 KB)
Additional file 9: Primary data of 91 significantly regulated transcripts corresponding to TIGS target genes tested in repeated bombardments of the QR screening. (XLSX 41 KB)
About this article
Cite this article
Douchkov, D., Lück, S., Johrde, A. et al. Discovery of genes affecting resistance of barley to adapted and non-adapted powdery mildew fungi. Genome Biol 15, 518 (2014). https://doi.org/10.1186/s13059-014-0518-8