Diversity at the Mla powdery mildew resistance locus from cultivated barley reveals sites of positive selection

Comparative genome analyses reveal sequence features reflecting distinct modes of host-adaptation between dicot and monocot powdery mildew

BACKGROUND

Powdery mildew (PM) is one of the most important and widespread plant diseases caused by biotrophic fungi. Notably, while monocot (grass) PM fungi exhibit high-level of host-specialization, many dicot PM fungi display a broad host range. To understand such distinct modes of host-adaptation, we sequenced the genomes of four dicot PM biotypes belonging to Golovinomyces cichoracearum or Oidium neolycopersici.

RESULTS

We compared genomes of the four dicot PM together with those of Blumeria graminis f.sp. hordei (both DH14 and RACE1 isolates), B. graminis f.sp. tritici, and Erysiphe necator infectious on barley, wheat and grapevine, respectively. We found that despite having a similar gene number (6620-6961), the PM genomes vary from 120 to 222 Mb in size. This high-level of genome size variation is indicative of highly differential transposon activities in the PM genomes. While the total number of genes in any given PM genome is only about half of that in the genomes of closely related ascomycete fungi, most (~ 93%) of the ascomycete core genes (ACGs) can be found in the PM genomes. Yet, 186 ACGs were found absent in at least two of the eight PM genomes, of which 35 are missing in some dicot PM biotypes, but present in the three monocot PM genomes, indicating remarkable, independent and perhaps ongoing gene loss in different PM lineages. Consistent with this, we found that only 4192 (3819 singleton) genes are shared by all the eight PM genomes, the remaining genes are lineage- or biotype-specific. Strikingly, whereas the three monocot PM genomes possess up to 661 genes encoding candidate secreted effector proteins (CSEPs) with families containing up to 38 members, all the five dicot PM fungi have only 116-175 genes encoding CSEPs with limited gene amplification.

CONCLUSIONS

Compared to monocot (grass) PM fungi, dicot PM fungi have a much smaller effectorome. This is consistent with their contrasting modes of host-adaption: while the monocot PM fungi show a high-level of host specialization, which may reflect an advanced host-pathogen arms race, the dicot PM fungi tend to practice polyphagy, which might have lessened selective pressure for escalating an with a particular host.

BED-domain-containing immune receptors confer diverse resistance spectra to yellow rust

Crop diseases reduce wheat yields by ~25% globally and thus pose a major threat to global food security¹. Genetic resistance can reduce crop losses in the field and can be selected through the use of molecular markers. However, genetic resistance often breaks down following changes in pathogen virulence, as experienced with the wheat yellow (stripe) rust fungus Puccinia striiformis f. sp. tritici (Pst)². This highlights the need to (1) identify genes that, alone or in combination, provide broad-spectrum resistance, and (2) increase our understanding of the underlying molecular modes of action. Here we report the isolation and characterization of three major yellow rust resistance genes (Yr7, Yr5 and YrSP) from hexaploid wheat (Triticum aestivum), each having a distinct recognition specificity. We show that Yr5, which remains effective to a broad range of Pst isolates worldwide, is closely related yet distinct from Yr7, whereas YrSP is a truncated version of Yr5 with 99.8% sequence identity. All three Yr genes belong to a complex resistance gene cluster on chromosome 2B encoding nucleotide-binding and leucine-rich repeat proteins (NLRs) with a non-canonical N-terminal zinc-finger BED domain³ that is distinct from those found in non-NLR wheat proteins. We developed diagnostic markers to accelerate haplotype analysis and for marker-assisted selection to expedite the stacking of the non-allelic Yr genes. Our results provide evidence that the BED-NLR gene architecture can provide effective field-based resistance to important fungal diseases such as wheat yellow rust.

Advances in Wheat and Pathogen Genomics: Implications for Disease Control

The gene pool of wheat and its wild and domesticated relatives contains a plethora of resistance genes that can be exploited to make wheat more resilient to pathogens. Only a few of these genes have been isolated and studied at the molecular level. In recent years, we have seen a shift from classical breeding to genomics-assisted breeding, which makes use of the enormous advancements in DNA sequencing and high-throughput molecular marker technologies for wheat improvement. These genomic advancements have the potential to transform wheat breeding in the near future and to significantly increase the speed and precision at which new cultivars can be bred. This review highlights the genomic improvements that have been made in wheat and its pathogens over the past years and discusses their implications for disease-resistance breeding.

Polymorphic residues in rice NLRs expand binding and response to effectors of the blast pathogen

Accelerated adaptive evolution is a hallmark of plant-pathogen interactions. Plant intracellular immune receptors (NLRs) often occur as allelic series with differential pathogen specificities. The determinants of this specificity remain largely unknown. Here, we unravelled the biophysical and structural basis of expanded specificity in the allelic rice NLR Pik, which responds to the effector AVR-Pik from the rice blast pathogen Magnaporthe oryzae. Rice plants expressing the Pikm allele resist infection by blast strains expressing any of three AVR-Pik effector variants, whereas those expressing Pikp only respond to one. Unlike Pikp, the integrated heavy metal-associated (HMA) domain of Pikm binds with high affinity to each of the three recognized effector variants, and variation at binding interfaces between effectors and Pikp-HMA or Pikm-HMA domains encodes specificity. By understanding how co-evolution has shaped the response profile of an allelic NLR, we highlight how natural selection drove the emergence of new receptor specificities. This work has implications for the engineering of NLRs with improved utility in agriculture.

Regulation and Evolution of NLR Genes: A Close Interconnection for Plant Immunity

NLR (NOD-like receptor) genes belong to one of the largest gene families in plants. Their role in plants' resistance to pathogens has been clearly described for many members of this gene family, and dysregulation or overexpression of some of these genes has been shown to induce an autoimmunity state that strongly affects plant growth and yield. For this reason, these genes have to be tightly regulated in their expression and activity, and several regulatory mechanisms are described here that tune their gene expression and protein levels. This gene family is subjected to rapid evolution, and to maintain diversity at NLRs, a plethora of genetic mechanisms have been identified as sources of variation. Interestingly, regulation of gene expression and evolution of this gene family are two strictly interconnected aspects. Indeed, some examples have been reported in which mechanisms of gene expression regulation have roles in promotion of the evolution of this gene family. Moreover, co-evolution of the NLR gene family and other gene families devoted to their control has been recently demonstrated, as in the case of miRNAs.

Isolate Specificity and Polygenic Inheritance of Resistance in Barley to Diverse Heterologous Puccinia striiformis Isolates

Barley is a host to Puccinia striiformis f. sp. hordei, and is an intermediate or near nonhost to the formae speciales adapted to wheat (P. striiformis f. sp. tritici) and to barley grass (P. striiformis f. sp. pseudo-hordei). The genetic basis of resistance to these forms of P. striiformis is not well understood. Accordingly, a recombinant inbred line (RIL) population was developed using a P. striiformis-susceptible accession (Biosaline-19) and the immune cultivar Pompadour. We investigated the genetic basis of resistance to four diverse P. striiformis isolates (P. striiformis f. sp. pseudo-hordei, and P. striiformis f. sp. tritici pathotypes 104 E137 A-, 134 E16 A+, and 64 E0 A-). and determined that the immunity in Pompadour at the seedling stage to the different P. striiformis isolates was due to quantitative trait loci (QTL) on chromosomes 1H, 3H, 5H, and 7H with both overlapping and distinct specificities. Further histological analysis confirmed the presence of isolate specificity. The RILs were also assessed in the field for resistance to P. striiformis f. sp. pseudo-hordei, P. striiformis f. sp. hordei, and the leaf rust pathogen (P. hordei) to identify pleiotropic QTL loci effective at the adult plant stage and determine whether the leaf rust resistance in Pompadour (Rph20) was also effective to P. striiformis. RILs that were seedling susceptible to P. striiformis f. sp. pseudo-hordei were resistant in the field, implicating the involvement of adult plant resistance (APR). Additional QTLs were identified on chromosome 7H at the same genetic position as Rph23 (APR to leaf rust), suggesting either pleiotropic resistance or the presence of a stripe rust resistance gene closely linked to or allelic with Rph23. Unlike many pleiotropic APR genes identified and isolated in wheat, our data suggest that the Rph20 locus does not confer resistance to the P. striiformis isolates used in this study (P. striiformis f. sp. hordei [χ² (independence) = 2.47 P > 0.12] and P. striiformis f. sp. pseudo-hordei [χ² (independence) = 0.42 P > 0.60]).

Distinct domains of the AVRPM3A2/F2 avirulence protein from wheat powdery mildew are involved in immune receptor recognition and putative effector function

Recognition of the AVRPM3A2/F2 avirulence protein from powdery mildew by the wheat PM3A/F immune receptor induces a hypersensitive response after co-expression in Nicotiana benthamiana. The molecular determinants of this interaction and how they shape natural AvrPm3a2/f2 allelic diversity are unknown. We sequenced the AvrPm3a2/f2 gene in a worldwide collection of 272 mildew isolates. Using the natural polymorphisms of AvrPm3a2/f2 as well as sequence information from related gene family members, we tested 85 single-residue-altered AVRPM3A2/F2 variants with PM3A, PM3F and PM3FL456P/Y458H (modified for improved signaling) in Nicotiana benthamiana for effects on recognition. An intact AvrPm3a2/f2 gene was found in all analyzed isolates and the protein variant recognized by PM3A/F occurred globally at high frequencies. Single-residue alterations in AVRPM3A2/F2 mostly disrupted, but occasionally enhanced, the recognition response by PM3A, PM3F and PM3FL456P/Y458H . Residues enhancing hypersensitive responses constituted a protein domain separate from both naturally occurring polymorphisms and positively selected residues of the gene family. These results demonstrate the utility of using gene family sequence diversity to screen residues for their role in recognition. This approach identified a putative interaction surface in AVRPM3A2/F2 not polymorphic in natural alleles. We conclude that molecular mechanisms besides recognition drive AvrPm3a2/f2 diversification.

Pathogen-inducible Ta-Lr34res expression in heterologous barley confers disease resistance without negative pleiotropic effects

Plant diseases are a serious threat to crop production. The informed use of naturally occurring disease resistance in plant breeding can greatly contribute to sustainably reduce yield losses caused by plant pathogens. The Ta-Lr34res gene encodes an ABC transporter protein and confers partial, durable, and broad spectrum resistance against several fungal pathogens in wheat. Transgenic barley lines expressing Ta-Lr34res showed enhanced resistance against powdery mildew and leaf rust of barley. While Ta-Lr34res is only active at adult stage in wheat, Ta-Lr34res was found to be highly expressed already at the seedling stage in transgenic barley resulting in severe negative effects on growth. Here, we expressed Ta-Lr34res under the control of the pathogen-inducible Hv-Ger4c promoter in barley. Sixteen independent barley transformants showed strong resistance against leaf rust and powdery mildew. Infection assays and growth parameter measurements were performed under standard glasshouse and near-field conditions using a convertible glasshouse. Two Hv-Ger4c::Ta-Lr34res transgenic events were analysed in detail. Plants of one transformation event had similar grain production compared to wild-type under glasshouse and near-field conditions. Our results showed that negative effects caused by constitutive high expression of Ta-Lr34res driven by the endogenous wheat promoter in barley can be eliminated by inducible expression without compromising disease resistance. These data demonstrate that Ta-Lr34res is agronomically useful in barley. We conclude that the generation of a large number of transformants in different barley cultivars followed by early field testing will allow identifying barley lines suitable for breeding.

Spectral Patterns Reveal Early Resistance Reactions of Barley Against Blumeria graminis f. sp. hordei

Differences in early plant-pathogen interactions are mainly characterized by using destructive methods. Optical sensors are advanced techniques for phenotyping host-pathogen interactions on different scales and for detecting subtle plant resistance responses against pathogens. A microscope with a hyperspectral camera was used to study interactions between Blumeria graminis f. sp. hordei and barley (Hordeum vulgare) genotypes with high susceptibility or resistance due to hypersensitive response (HR) and papilla formation. Qualitative and quantitative assessment of pathogen development was used to explain changes in hyperspectral signatures. Within 48 h after inoculation, genotype-specific changes in the green and red range (500 to 690 nm) and a blue shift of the red-edge inflection point were observed. Manual analysis indicated resistance-specific reflectance patterns from 1 to 3 days after inoculation. These changes could be linked to host plant modifications depending on individual host-pathogen interactions. Retrospective analysis of hyperspectral images revealed spectral characteristics of HR against B. graminis f. sp. hordei. For early HR detection, an advanced data mining approach localized HR spots before they became visible on the RGB images derived from hyperspectral imaging. The link among processes during pathogenesis and host resistance to changes in hyperspectral signatures provide evidence that sensor-based phenotyping is suitable to advance time-consuming and cost-expensive visual rating of plant disease resistances.

Interchromosomal Transfer of Immune Regulation During Infection of Barley with the Powdery Mildew Pathogen

Powdery mildew pathogens colonize over 9500 plant species, causing critical yield loss. The Ascomycete fungus, Blumeria graminis f. sp. hordei (Bgh), causes powdery mildew disease in barley (Hordeum vulgare L.). Successful infection begins with penetration of host epidermal cells, culminating in haustorial feeding structures, facilitating delivery of fungal effectors to the plant and exchange of nutrients from host to pathogen. We used expression Quantitative Trait Locus (eQTL) analysis to dissect the temporal control of immunity-associated gene expression in a doubled haploid barley population challenged with Bgh. Two highly significant regions possessing trans eQTL were identified near the telomeric ends of chromosomes (Chr) 2HL and 1HS. Within these regions reside diverse resistance loci derived from barley landrace H. laevigatum (MlLa) and H. vulgare cv. Algerian (Mla1), which associate with the altered expression of 961 and 3296 genes during fungal penetration of the host and haustorial development, respectively. Regulatory control of transcript levels for 299 of the 961 genes is reprioritized from MlLa on 2HL to Mla1 on 1HS as infection progresses, with 292 of the 299 alternating the allele responsible for higher expression, including Adaptin Protein-2 subunit μ AP2M and Vesicle Associated Membrane Protein VAMP72 subfamily members VAMP721/722. AP2M mediates effector-triggered immunity (ETI) via endocytosis of plasma membrane receptor components. VAMP721/722 and SNAP33 form a Soluble N-ethylmaleimide-sensitive factor Attachment Protein REceptor (SNARE) complex with SYP121 (PEN1), which is engaged in pathogen associated molecular pattern (PAMP)-triggered immunity via exocytosis. We postulate that genes regulated by alternate chromosomal positions are repurposed as part of a conserved immune complex to respond to different pathogen attack scenarios.

Ecological genomics of tropical trees: how local population size and allelic diversity of resistance genes relate to immune responses, cosusceptibility to pathogens, and negative density dependence

In tropical forests, rarer species show increased sensitivity to species-specific soil pathogens and more negative effects of conspecific density on seedling survival (NDD). These patterns suggest a connection between ecology and immunity, perhaps because small population size disproportionately reduces genetic diversity of hyperdiverse loci such as immunity genes. In an experiment examining seedling roots from six species in one tropical tree community, we found that smaller populations have reduced amino acid diversity in pathogen resistance (R) genes but not the transcriptome in general. Normalized R gene amino acid diversity varied with local abundance and prior measures of differences in sensitivity to conspecific soil and NDD. After exposure to live soil, species with lower R gene diversity had reduced defence gene induction, more cosusceptibility of maternal cohorts to colonization by potentially pathogenic fungi, reduced root growth arrest (an R gene-mediated response) and their root-associated fungi showed lower induction of self-defence (antioxidants). Local abundance was not related to the ability to induce immune responses when pathogen recognition was bypassed by application of salicylic acid, a phytohormone that activates defence responses downstream of R gene signalling. These initial results support the hypothesis that smaller local tree populations have reduced R gene diversity and recognition-dependent immune responses, along with greater cosusceptibility to species-specific pathogens that may facilitate disease transmission and NDD. Locally rare species may be less able to increase their equilibrium abundance without genetic boosts to defence via immigration of novel R gene alleles from a larger and more diverse regional population.

AvrPm2 encodes an RNase-like avirulence effector which is conserved in the two different specialized forms of wheat and rye powdery mildew fungus

There is a large diversity of genetically defined resistance genes in bread wheat against the powdery mildew pathogen Blumeria graminis (B. g.) f. sp. tritici. Many confer race-specific resistance to this pathogen, but until now only the mildew avirulence gene AvrPm3a2/f2 that is recognized by Pm3a/f was known molecularly. We performed map-based cloning and genome-wide association studies to isolate a candidate for the mildew avirulence gene AvrPm2. We then used transient expression assays in Nicotiana benthamiana to demonstrate specific and strong recognition of AvrPm2 by Pm2. The virulent AvrPm2 allele arose from a conserved 12 kb deletion, while there is no protein sequence diversity in the gene pool of avirulent B. g. tritici isolates. We found one polymorphic AvrPm2 allele in B. g. triticale and one orthologue in B. g. secalis and both are recognized by Pm2. AvrPm2 belongs to a small gene family encoding structurally conserved RNase-like effectors, including Avra13 from B. g. hordei, the cognate Avr of the barley resistance gene Mla13. These results demonstrate the conservation of functional avirulence genes in two cereal powdery mildews specialized on different hosts, thus providing a possible explanation for successful introgression of resistance genes from rye or other grass relatives to wheat.

The CC domain structure from the wheat stem rust resistance protein Sr33 challenges paradigms for dimerization in plant NLR proteins

Plants use intracellular immunity receptors, known as nucleotide-binding oligomerization domain-like receptors (NLRs), to recognize specific pathogen effector proteins and induce immune responses. These proteins provide resistance to many of the world's most destructive plant pathogens, yet we have a limited understanding of the molecular mechanisms that lead to defense signaling. We examined the wheat NLR protein, Sr33, which is responsible for strain-specific resistance to the wheat stem rust pathogen, Puccinia graminis f. sp. tritici We present the solution structure of a coiled-coil (CC) fragment from Sr33, which adopts a four-helix bundle conformation. Unexpectedly, this structure differs from the published dimeric crystal structure of the equivalent region from the orthologous barley powdery mildew resistance protein, MLA10, but is similar to the structure of the distantly related potato NLR protein, Rx. We demonstrate that these regions are, in fact, largely monomeric and adopt similar folds in solution in all three proteins, suggesting that the CC domains from plant NLRs adopt a conserved fold. However, larger C-terminal fragments of Sr33 and MLA10 can self-associate both in vitro and in planta, and this self-association correlates with their cell death signaling activity. The minimal region of the CC domain required for both cell death signaling and self-association extends to amino acid 142, thus including 22 residues absent from previous biochemical and structural protein studies. These data suggest that self-association of the minimal CC domain is necessary for signaling but is likely to involve a different structural basis than previously suggested by the MLA10 crystallographic dimer.

Allelic barley MLA immune receptors recognize sequence-unrelated avirulence effectors of the powdery mildew pathogen

Disease-resistance genes encoding intracellular nucleotide-binding domain and leucine-rich repeat proteins (NLRs) are key components of the plant innate immune system and typically detect the presence of isolate-specific avirulence (AVR) effectors from pathogens. NLR genes define the fastest-evolving gene family of flowering plants and are often arranged in gene clusters containing multiple paralogs, contributing to copy number and allele-specific NLR variation within a host species. Barley mildew resistance locus a (Mla) has been subject to extensive functional diversification, resulting in allelic resistance specificities each recognizing a cognate, but largely unidentified, AVR_a gene of the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). We applied a transcriptome-wide association study among 17 Bgh isolates containing different AVR_a genes and identified AVR_a1 and AVR_a13, encoding candidate-secreted effectors recognized by Mla1 and Mla13 alleles, respectively. Transient expression of the effector genes in barley leaves or protoplasts was sufficient to trigger Mla1 or Mla13 allele-specific cell death, a hallmark of NLR receptor-mediated immunity. AVR_a1 and AVR_a13 are phylogenetically unrelated, demonstrating that certain allelic MLA receptors evolved to recognize sequence-unrelated effectors. They are ancient effectors because corresponding loci are present in wheat powdery mildew. AVR_A1 recognition by barley MLA1 is retained in transgenic Arabidopsis, indicating that AVR_A1 directly binds MLA1 or that its recognition involves an evolutionarily conserved host target of AVR_A1 Furthermore, analysis of transcriptome-wide sequence variation among the Bgh isolates provides evidence for Bgh population structure that is partially linked to geographic isolation.

Molecular genetics and evolution of disease resistance in cereals

Contents 320 I. 320 II. 321 III. 321 IV. 322 V. 324 VI. 328 VII. 329 330 References 330 SUMMARY: Cereal crops produce a large part of the globally consumed food and feed. Because of the constant presence of devastating pathogens, the molecular characterization of disease resistance is a major research area and highly relevant for breeding. There has been recent and accelerating progress in the understanding of three distinct resistance mechanisms in cereals: resistance conferred by plasma membrane-localized receptor proteins; race-specific resistance conferred by intracellular immune receptors; and quantitative disease resistance. Intracellular immune receptors provide a particularly rich source for evolutionary studies, and have, for example, resulted in the recent discovery of a novel detection mechanism based on integrated decoy domains. Evolutionary studies have also revealed the origins of active resistance genes in both wild progenitors of today's cereals as well as in cultivated forms. In addition, independent evolution of orthologous genes in related cereals has resulted in resistance to different pathogen species. Quantitative resistance genes have been best characterized in wheat. The quantitative resistance genes identified so far in wheat encode transporter proteins or unusual kinase proteins. The recent discoveries in these three different resistance mechanisms have contributed to the basic molecular understanding of cereal immunity against pathogens and have suggested novel applications for resistance breeding.

Cytosolic activation of cell death and stem rust resistance by cereal MLA-family CC-NLR proteins

Plants possess intracellular immune receptors designated "nucleotide-binding domain and leucine-rich repeat" (NLR) proteins that translate pathogen-specific recognition into disease-resistance signaling. The wheat immune receptors Sr33 and Sr50 belong to the class of coiled-coil (CC) NLRs. They confer resistance against a broad spectrum of field isolates of Puccinia graminis f. sp. tritici, including the Ug99 lineage, and are homologs of the barley powdery mildew-resistance protein MLA10. Here, we show that, similarly to MLA10, the Sr33 and Sr50 CC domains are sufficient to induce cell death in Nicotiana benthamiana Autoactive CC domains and full-length Sr33 and Sr50 proteins self-associate in planta In contrast, truncated CC domains equivalent in size to an MLA10 fragment for which a crystal structure was previously determined fail to induce cell death and do not self-associate. Mutations in the truncated region also abolish self-association and cell-death signaling. Analysis of Sr33 and Sr50 CC domains fused to YFP and either nuclear localization or nuclear export signals in N benthamiana showed that cell-death induction occurs in the cytosol. In stable transgenic wheat plants, full-length Sr33 proteins targeted to the cytosol provided rust resistance, whereas nuclear-targeted Sr33 was not functional. These data are consistent with CC-mediated induction of both cell-death signaling and stem rust resistance in the cytosolic compartment, whereas previous research had suggested that MLA10-mediated cell-death and disease resistance signaling occur independently, in the cytosol and nucleus, respectively.

Geographically Distinct and Domain-Specific Sequence Variations in the Alleles of Rice Blast Resistance Gene Pib

Rice blast is caused by Magnaporthe oryzae, which is the most destructive fungal pathogen affecting rice growing regions worldwide. The rice blast resistance gene Pib confers broad-spectrum resistance against Southeast Asian M. oryzae races. We investigated the allelic diversity of Pib in rice germplasm originating from 12 major rice growing countries. Twenty-five new Pib alleles were identified that have unique single nucleotide polymorphisms (SNPs), insertions and/or deletions, in addition to the polymorphic nucleotides that are shared between the different alleles. These partially or completely shared polymorphic nucleotides indicate frequent sequence exchange events between the Pib alleles. In some of the new Pib alleles, nucleotide diversity is high in the LRR domain, whereas, in others it is distributed among the NB-ARC and LRR domains. Most of the polymorphic amino acids in LRR and NB-ARC2 domains are predicted as solvent-exposed. Several of the alleles and the unique SNPs are country specific, suggesting a diversifying selection of alleles in various geographical locations in response to the locally prevalent M. oryzae population. Together, the new Pib alleles are an important genetic resource for rice blast resistance breeding programs and provide new information on rice-M. oryzae interactions at the molecular level.

Diversity and Evolution of Disease Resistance Genes in Barley (Hordeum vulgare L.)

Plant disease resistance genes (R-genes) play a critical role in the defense response to pathogens. Barley is one of the most important cereal crops, having a genome recently made available, for which the diversity and evolution of R-genes are not well understood. The main objectives of this research were to conduct a genome-wide identification of barley Coiled-coil, Nucleotide-binding site, Leucine-rich repeat (CNL) genes and elucidate their evolutionary history. We employed a Hidden Markov Model using 52 Arabidopsis thaliana CNL reference sequences and analyzed for phylogenetic relationships, structural variation, and gene clustering. We identified 175 barley CNL genes nested into three clades, showing (a) evidence of an expansion of the CNL-C clade, primarily due to tandem duplications; (b) very few members of clade CNL-A and CNL-B; and (c) a complete absence of clade CNL-D. Our results also showed that several of the previously identified mildew locus A (MLA) genes may be allelic variants of two barley CNL genes, MLOC_66581 and MLOC_10425, which respond to powdery mildew. Approximately 23% of the barley CNL genes formed 15 gene clusters located in the extra-pericentromeric regions on six of the seven chromosomes; more than half of the clustered genes were located on chromosomes 1H and 7H. Higher average numbers of exons and multiple splice variants in barley relative to those in Arabidopsis and rice may have contributed to a diversification of the CNL-C members. These results will help us understand the evolution of R-genes with potential implications for developing durable resistance in barley cultivars.

Arabidopsis thaliana DM2h (R8) within the Landsberg RPP1-like Resistance Locus Underlies Three Different Cases of EDS1-Conditioned Autoimmunity

Plants have a large panel of nucleotide-binding/leucine rich repeat (NLR) immune receptors which monitor host interference by diverse pathogen molecules (effectors) and trigger disease resistance pathways. NLR receptor systems are necessarily under tight control to mitigate the trade-off between induced defenses and growth. Hence, mis-regulated NLRs often cause autoimmunity associated with stunting and, in severe cases, necrosis. Nucleocytoplasmic ENHANCED DISEASE SUSCEPTIBILITY1 (EDS1) is indispensable for effector-triggered and autoimmune responses governed by a family of Toll-Interleukin1-Receptor-related NLR receptors (TNLs). EDS1 operates coincidently or immediately downstream of TNL activation to transcriptionally reprogram cells for defense. We show here that low levels of nuclear-enforced EDS1 are sufficient for pathogen resistance in Arabidopsis thaliana, without causing negative effects. Plants expressing higher nuclear EDS1 amounts have the genetic, phenotypic and transcriptional hallmarks of TNL autoimmunity. In a screen for genetic suppressors of nuclear EDS1 autoimmunity, we map multiple, independent mutations to one gene, DM2h, lying within the polymorphic DANGEROUS MIX2 cluster of TNL RPP1-like genes from A. thaliana accession Landsberg erecta (Ler). The DM2 locus is a known hotspot for deleterious epistatic interactions leading to immune-related incompatibilities between A. thaliana natural accessions. We find that DM2h^Ler underlies two further genetic incompatibilities involving the RPP1-like^Ler locus and EDS1. We conclude that the DM2h^Ler TNL protein and nuclear EDS1 cooperate, directly or indirectly, to drive cells into an immune response at the expense of growth. A further conclusion is that regulating the available EDS1 nuclear pool is fundamental for maintaining homeostatic control of TNL immune pathways.

Plants tune their cellular and developmental programs to different environmental stimuli. Central players in the plant biotic stress response network are intracellular NLR receptors which intercept specific disease-inducing molecules (effectors) produced by pathogenic microbes. Variation in NLR gene repertoires between plant genetic lines is driven by pathogen selection pressure. One evolutionary question is how new, functional NLRs are assembled within a plant genome without mis-activating defense pathways, which can have strong negative effects on growth and fitness. This study focuses on a large, polymorphic sub-class of NLR receptors called TNLs present in dicotyledenous plant lineages. TNL receptors confer immunity to a broad range of pathogens. They also frequently underlie autoimmunity caused by their mis-regulation or deleterious allelic interactions with other genes in crosses between different genetic lines (hybrid incompatibility, HI). TNL pathogen-triggered and autoimmune responses require the conserved nucleocytoplasmic protein EDS1 to transcriptionally reprogram cells for defense. We discover in Arabidopsis thaliana that high levels of nuclear-enriched EDS1 induce transcriptional activation of defenses and growth inhibition without a pathogen effector stimulus. In a mutational screen, we identify one rapidly evolving TNL gene, DM2h^Ler, as a driver of nuclear EDS1 autoimmunity. DM2h^Ler also contributes to two separate cases of EDS1-dependent autoimmunity. Genetic cooperation between DM2h^Ler and EDS1 suggests a functional relationship in the transcriptional feed-forward regulation of defense pathways.

Isolation and fine mapping of Rps6: an intermediate host resistance gene in barley to wheat stripe rust

We uncouple host and nonhost resistance in barley to Puccinia striiformis ff. spp. hordei and tritici . We isolate, fine map, and physically anchor Rps6 to chromosome 7H in barley. A plant may be considered a nonhost of a pathogen if all known genotypes of a plant species are resistant to all known isolates of a pathogen species. However, if a small number of genotypes are susceptible to some known isolates of a pathogen species this plant may be considered an intermediate host. Barley (Hordeum vulgare) is an intermediate host for Puccinia striiformis f. sp. tritici (Pst), the causal agent of wheat stripe rust. We wanted to understand the genetic architecture underlying resistance to Pst and to determine whether any overlap exists with resistance to the host pathogen, Puccinia striiformis f. sp. hordei (Psh). We mapped Pst resistance to chromosome 7H and show that host and intermediate host resistance is genetically uncoupled. Therefore, we designate this resistance locus Rps6. We used phenotypic and genotypic selection on F2:3 families to isolate Rps6 and fine mapped the locus to a 0.1 cM region. Anchoring of the Rps6 locus to the barley physical map placed the region on a single fingerprinted contig spanning a physical region of 267 kb. Efforts are now underway to sequence the minimal tiling path and to delimit the physical region harboring Rps6. This will facilitate additional marker development and permit identification of candidate genes in the region.

Downregulation of a barley (Hordeum vulgare) leucine-rich repeat, non-arginine-aspartate receptor-like protein kinase reduces expression of numerous genes involved in plant pathogen defense

Pattern recognition receptors represent a first line of plant defense against pathogens. Comparing the flag leaf transcriptomes of barley (Hordeum vulgare L.) near-isogenic lines varying in the allelic state of a locus controlling senescence, we have previously identified a leucine-rich repeat receptor-like protein kinase gene (LRR-RLK; GenBank accession: AK249842), which was strongly upregulated in leaves of early-as compared to late-senescing germplasm. Bioinformatic analysis indicated that this gene codes for a subfamily XII, non-arginine-aspartate (non-RD) LRR-RLK. Virus-induced gene silencing resulted in a two-fold reduction of transcript levels as compared to controls. Transcriptomic comparison of leaves from untreated plants, from plants treated with virus only without any plant sequences (referred to as 'empty virus' control), and from plants in which AK249842 expression was knocked down identified numerous genes involved in pathogen defense. These genes were strongly induced in 'empty virus' as compared to untreated controls, but their expression was significantly reduced (again compared to 'empty virus' controls) when AK249842 was knocked down, indicating that their expression partially depends on the LRR-RLK investigated here. Expression analysis, using datasets from BarleyBase/PLEXdb, demonstrated that AK249842 transcript levels are heavily influenced by the allelic state of the well-characterized mildew resistance a (Mla) locus, and that the gene is induced after powdery mildew and stem rust infection. Together, our data suggest that AK249842 is a barley pattern recognition receptor with a tentative role in defense against fungal pathogens, setting the stage for its full functional characterization.

Avirulence Genes in Cereal Powdery Mildews: The Gene-for-Gene Hypothesis 2.0

The gene-for-gene hypothesis states that for each gene controlling resistance in the host, there is a corresponding, specific gene controlling avirulence in the pathogen. Allelic series of the cereal mildew resistance genes Pm3 and Mla provide an excellent system for genetic and molecular analysis of resistance specificity. Despite this opportunity for molecular research, avirulence genes in mildews remain underexplored. Earlier work in barley powdery mildew (B.g. hordei) has shown that the reaction to some Mla resistance alleles is controlled by multiple genes. Similarly, several genes are involved in the specific interaction of wheat mildew (B.g. tritici) with the Pm3 allelic series. We found that two mildew genes control avirulence on Pm3f: one gene is involved in recognition by the resistance protein as demonstrated by functional studies in wheat and the heterologous host Nicotiana benthamiana. A second gene is a suppressor, and resistance is only observed in mildew genotypes combining the inactive suppressor and the recognized Avr. We propose that such suppressor/avirulence gene combinations provide the basis of specificity in mildews. Depending on the particular gene combinations in a mildew race, different genes will be genetically identified as the "avirulence" gene. Additionally, the observation of two LINE retrotransposon-encoded avirulence genes in B.g. hordei further suggests that the control of avirulence in mildew is more complex than a canonical gene-for-gene interaction. To fully understand the mildew-cereal interactions, more knowledge on avirulence determinants is needed and we propose ways how this can be achieved based on recent advances in the field.

The wheat Sr50 gene reveals rich diversity at a cereal disease resistance locus

We identify the wheat stem rust resistance gene Sr50 (using physical mapping, mutation and complementation) as homologous to barley Mla, encoding a coiled-coil nucleotide-binding leucine-rich repeat (CC-NB-LRR) protein. We show that Sr50 confers a unique resistance specificity different from Sr31 and other genes on rye chromosome 1RS, and is effective against the broadly virulent Ug99 race lineage. Extensive haplotype diversity at the rye Sr50 locus holds promise for mining effective resistance genes.

Identification of novel alleles of the rice blast resistance gene Pi54

Rice blast is one of the most devastating rice diseases and continuous resistance breeding is required to control the disease. The rice blast resistance gene Pi54 initially identified in an Indian cultivar confers broad-spectrum resistance in India. We explored the allelic diversity of the Pi54 gene among 885 Indian rice genotypes that were found resistant in our screening against field mixture of naturally existing M. oryzae strains as well as against five unique strains. These genotypes are also annotated as rice blast resistant in the International Rice Genebank database. Sequence-based allele mining was used to amplify and clone the Pi54 allelic variants. Nine new alleles of Pi54 were identified based on the nucleotide sequence comparison to the Pi54 reference sequence as well as to already known Pi54 alleles. DNA sequence analysis of the newly identified Pi54 alleles revealed several single polymorphic sites, three double deletions and an eight base pair deletion. A SNP-rich region was found between a tyrosine kinase phosphorylation site and the nucleotide binding site (NBS) domain. Together, the newly identified Pi54 alleles expand the allelic series and are candidates for rice blast resistance breeding programs.

Emergence and Spread of New Races of Wheat Stem Rust Fungus: Continued Threat to Food Security and Prospects of Genetic Control

Race Ug99 (TTKSK) of Puccinia graminis f. sp. tritici, detected in Uganda in 1998, has been recognized as a serious threat to food security because it possesses combined virulence to a large number of resistance genes found in current widely grown wheat (Triticum aestivum) varieties and germplasm, leading to its potential for rapid spread and evolution. Since its initial detection, variants of the Ug99 lineage of stem rust have been discovered in Eastern and Southern African countries, Yemen, Iran, and Egypt. To date, eight races belonging to the Ug99 lineage are known. Increased pathogen monitoring activities have led to the identification of other races in Africa and Asia with additional virulence to commercially important resistance genes. This has led to localized but severe stem rust epidemics becoming common once again in East Africa due to the breakdown of race-specific resistance gene SrTmp, which was deployed recently in the 'Digalu' and 'Robin' varieties in Ethiopia and Kenya, respectively. Enhanced research in the last decade under the umbrella of the Borlaug Global Rust Initiative has identified various race-specific resistance genes that can be utilized, preferably in combinations, to develop resistant varieties. Research and development of improved wheat germplasm with complex adult plant resistance (APR) based on multiple slow-rusting genes has also progressed. Once only the Sr2 gene was known to confer slow rusting APR; now, four more genes-Sr55, Sr56, Sr57, and Sr58-have been characterized and additional quantitative trait loci identified. Cloning of some rust resistance genes opens new perspectives on rust control in the future through the development of multiple resistance gene cassettes. However, at present, disease-surveillance-based chemical control, large-scale deployment of new varieties with multiple race-specific genes or adequate levels of APR, and reducing the cultivation of susceptible varieties in rust hot-spot areas remains the best stem rust management strategy.

Natural loss-of-function mutation of EDR1 conferring resistance to tomato powdery mildew in Arabidopsis thaliana accession C24

To screen for potentially novel types of resistance to tomato powdery mildew Oidium neolycopersici, a disease assay was performed on 123 Arabidopsis thaliana accessions. Forty accessions were fully resistant, and one, C24, was analysed in detail. By quantitative trait locus (QTL) analysis of an F2 population derived from C24 × Sha (susceptible accession), two QTLs associated with resistance were identified in C24. Fine mapping of QTL-1 on chromosome 1 delimited the region to an interval of 58 kb encompassing 15 candidate genes. One of these was Enhanced Disease Resistance 1 (EDR1). Evaluation of the previously obtained edr1 mutant of Arabidopsis accession Col-0, which was identified because of its resistance to powdery mildew Golovinomyces cichoracearum, showed that it also displayed resistance to O. neolycopersici. Sequencing of EDR1 in our C24 germplasm (referred to as C24-W) revealed two missing nucleotides in the second exon of EDR1 resulting in a premature stop codon. Remarkably, C24 obtained from other laboratories does not contain the EDR1 mutation. To verify the identity of C24-W, a DNA region containing a single nucleotide polymorphism (SNP) unique to C24 was sequenced showing that C24-W contains the C24-specific nucleotide. C24-W showed enhanced resistance to O. neolycopersici compared with C24 not containing the edr1 mutation. Furthermore, C24-W displayed a dwarf phenotype, which was not associated with the mutation in EDR1 and was not caused by the differential accumulation of pathogenesis-related genes. In conclusion, we identified a natural edr1 mutant in the background of C24.

The ecological genomic basis of salinity adaptation in Tunisian Medicago truncatula

Background

As our world becomes warmer, agriculture is increasingly impacted by rising soil salinity and understanding plant adaptation to salt stress can help enable effective crop breeding. Salt tolerance is a complex plant phenotype and we know little about the pathways utilized by naturally tolerant plants. Legumes are important species in agricultural and natural ecosystems, since they engage in symbiotic nitrogen-fixation, but are especially vulnerable to salinity stress.

Results

Our studies of the model legume Medicago truncatula in field and greenhouse settings demonstrate that Tunisian populations are locally adapted to saline soils at the metapopulation level and that saline origin genotypes are less impacted by salt than non-saline origin genotypes; these populations thus likely contain adaptively diverged alleles. Whole genome resequencing of 39 wild accessions reveals ongoing migration and candidate genomic regions that assort non-randomly with soil salinity. Consistent with natural selection acting at these sites, saline alleles are typically rare in the range-wide species' gene pool and are also typically derived relative to the sister species M. littoralis. Candidate regions for adaptation contain genes that regulate physiological acclimation to salt stress, such as abscisic acid and jasmonic acid signaling, including a novel salt-tolerance candidate orthologous to the uncharacterized gene AtCIPK21. Unexpectedly, these regions also contain biotic stress genes and flowering time pathway genes. We show that flowering time is differentiated between saline and non-saline populations and may allow salt stress escape.

Conclusions

This work nominates multiple potential pathways of adaptation to naturally stressful environments in a model legume. These candidates point to the importance of both tolerance and avoidance in natural legume populations. We have uncovered several promising targets that could be used to breed for enhanced salt tolerance in crop legumes to enhance food security in an era of increasing soil salinization.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1160) contains supplementary material, which is available to authorized users.

The miR9863 family regulates distinct Mla alleles in barley to attenuate NLR receptor-triggered disease resistance and cell-death signaling

Barley (Hordeum vulgare L.) Mla alleles encode coiled-coil (CC), nucleotide binding, leucine-rich repeat (NB-LRR) receptors that trigger isolate-specific immune responses against the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh). How Mla or NB-LRR genes in grass species are regulated at post-transcriptional level is not clear. The microRNA family, miR9863, comprises four members that differentially regulate distinct Mla alleles in barley. We show that miR9863 members guide the cleavage of Mla1 transcripts in barley, and block or reduce the accumulation of MLA1 protein in the heterologous Nicotiana benthamiana expression system. Regulation specificity is determined by variation in a unique single-nucleotide-polymorphism (SNP) in mature miR9863 family members and two SNPs in the Mla miR9863-binding site that separates these alleles into three groups. Further, we demonstrate that 22-nt miR9863s trigger the biogenesis of 21-nt phased siRNAs (phasiRNAs) and together these sRNAs form a feed-forward regulation network for repressing the expression of group I Mla alleles. Overexpression of miR9863 members specifically attenuates MLA1, but not MLA10-triggered disease resistance and cell-death signaling. We propose a key role of the miR9863 family in dampening immune response signaling triggered by a group of MLA immune receptors in barley.

Plants rely on cell-surface and intracellular immune receptors to sense pathogen invasion and to mediate defense responses. However, uncontrolled activation of immune responses is harmful to plant growth and development. Small RNAs have recently been shown to fine-tune the expression of intracellular immune receptors and contribute to the regulation of defense signaling in dicot plants, while similar processes have not been well documented in monocot grain crops, such as barley and wheat. Here, we show that, in barley, some members of the miR9863 family target a subset of Mla alleles that confer race-specific disease resistance to the powdery mildew fungus. These miRNAs act on Mla transcripts by cleavage and translational repression. Production of a type of trans-acting small RNAs, designated as phasiRNAs, enhances the effects of miRNA regulation on Mla targets. We propose that Mla-mediated immune signaling is fine-tuned by the miRNAs at later stage of MLA activation to avoid overloading of immune responses in barley cells.

The Vat locus encodes for a CC-NBS-LRR protein that confers resistance to Aphis gossypii infestation and A. gossypii-mediated virus resistance

Aphis gossypii is a polyphagous sucking aphid and a vector for many viruses. In Cucumis melo, a dominant locus, Vat, confers a high level of resistance to Aphis gossypii infestation and to viruses transmitted by this vector. To investigate the mechanism underlying this double resistance, we first genetically dissected the Vat locus. We delimited the double resistance to a single gene that encodes for a coiled-coil-nucleotide-binding-site-leucine-rich repeat (CC-NBS-LRR) protein type. To validate the genetic data, transgenic lines expressing the Vat gene were generated and assessed for the double resistance. In this analysis, Vat-transgenic plants were resistant to A. gossypii infestation as well as A. gossypii-mediated virus transmission. When the plants were infected mechanically, virus infection occurred on both transgenic and non-transgenic control plants. These results confirmed that the cloned CC-NBS-LRR gene mediates both resistance to aphid infestation and virus infection using A. gossypii as a vector. This resistance also invokes a separate recognition and response phases in which the recognition phase involves the interaction of an elicitor molecule from the aphid and Vat from the plant. The response phase is not specific and blocks both aphid infestation and virus infection. Sequence analysis of Vat alleles suggests a major role of an unusual conserved LRR repeat in the recognition of A. gossypii.

Three-dimensional modeling and diversity analysis reveals distinct AVR recognition sites and evolutionary pathways in wild and domesticated wheat Pm3 R genes

The Pm3 gene confers resistance against wheat powdery mildew. Studies of Pm3 diversity have shown that Pm3 alleles isolated from southern populations of wild emmer wheat located in Lebanon, Jordan, Israel, and Syria are more diverse and more distant from bread wheat alleles than alleles from the northern wild wheat populations located in Turkey, Iran, and Iraq. Therefore, southern populations from Israel were studied extensively to reveal novel Pm3 alleles that are absent from the cultivated gene pool. Candidate Pm3 genes were isolated via a polymerase chain reaction cloning approach. Known and newly identified Pm3 genes were subjected to variation analysis and polymorphic amino acid residues were superimposed on a three-dimensional (3D) model of PM3. The region of highest interspecies diversity between Triticum aestivum and T. dicoccoides lies in leucine-rich repeats (LRR) 19 to 24, whereas most intraspecies diversity in T. aestivum is located in LRR 25 to 28. Interestingly, these two regions are separated by one large LRR whose propensity for flexibility facilitates the conformation of the PM3 LRR domain into two differently structured models. The combination of evolutionary and protein 3D structure analysis revealed that Pm3 genes in wild and domesticated wheat show different evolutionary histories which might have been triggered through different interactions with the powdery mildew pathogen.

Mla- and Rom1-mediated control of microRNA398 and chloroplast copper/zinc superoxide dismutase regulates cell death in response to the barley powdery mildew fungus

• Barley (Hordeum vulgare L.) Mildew resistance locus a (Mla) confers allele-specific interactions with natural variants of the ascomycete fungus Blumeria graminis f. sp. hordei (Bgh), the causal agent of powdery mildew disease. Significant reprogramming of Mla-mediated gene expression occurs upon infection by this obligate biotrophic pathogen. • We utilized a proteomics-based approach, combined with barley mla, required for Mla12 resistance1 (rar1), and restoration of Mla resistance1 (rom1) mutants, to identify components of Mla-directed signaling. • Loss-of-function mutations in Mla and Rar1 both resulted in the reduced accumulation of chloroplast copper/zinc superoxide dismutase 1 (HvSOD1), whereas loss of function in Rom1 re-established HvSOD1 levels. In addition, both Mla and Rom1 negatively regulated hvu-microRNA398 (hvu-miR398), and up-regulation of miR398 was coupled to reduced HvSOD1 expression. Barley stripe mosaic virus (BSMV)-mediated over-expression of both barley and Arabidopsis miR398 repressed accumulation of HvSOD1, and BSMV-induced gene silencing of HvSod1 impeded Mla-triggered H₂O₂ and hypersensitive reaction (HR) at barley-Bgh interaction sites. • These data indicate that Mla- and Rom1-regulated hvu-miR398 represses HvSOD1 accumulation, influencing effector-induced HR in response to the powdery mildew fungus.

Rye Pm8 and wheat Pm3 are orthologous genes and show evolutionary conservation of resistance function against powdery mildew

The improvement of wheat through breeding has relied strongly on the use of genetic material from related wild and domesticated grass species. The 1RS chromosome arm from rye was introgressed into wheat and crossed into many wheat lines, as it improves yield and fungal disease resistance. Pm8 is a powdery mildew resistance gene on 1RS which, after widespread agricultural cultivation, is now widely overcome by adapted mildew races. Here we show by homology-based cloning and subsequent physical and genetic mapping that Pm8 is the rye orthologue of the Pm3 allelic series of mildew resistance genes in wheat. The cloned gene was functionally validated as Pm8 by transient, single-cell expression analysis and stable transformation. Sequence analysis revealed a complex mosaic of ancient haplotypes among Pm3- and Pm8-like genes from different members of the Triticeae. These results show that the two genes have evolved independently after the divergence of the species 7.5 million years ago and kept their function in mildew resistance. During this long time span the co-evolving pathogens have not overcome these genes, which is in strong contrast to the breakdown of Pm8 resistance since its introduction into commercial wheat 70 years ago. Sequence comparison revealed that evolutionary pressure acted on the same subdomains and sequence features of the two orthologous genes. This suggests that they recognize directly or indirectly the same pathogen effectors that have been conserved in the powdery mildews of wheat and rye.

Partitioning, repressing and derepressing: dynamic regulations in MLA immune receptor triggered defense signaling

Plants and animals have evolved intracellular nucleotide-binding domain and leucine-rich repeat-containing immune receptors (NLRs) to perceive non-self and trigger immune responses. Plant NLRs detect strain-specific pathogen effectors and activate immune signaling leading to extensive transcriptional reprogramming and termination of pathogen infection. Here we review the recent findings in barley MLA immune receptor mediated immune responses against the barley powdery mildew fungus. We focus on nucleocytoplasmic partitioning of immune receptor, bifurcation of immune signaling, transcriptional repression and derepression connecting receptor activation to immune responses. We also discuss similar findings from other plant NLRs where appropriate.

Evolution and Conservation of Plant NLR Functions

In plants and animals, nucleotide-binding domain and leucine-rich repeats (NLR)-containing proteins play pivotal roles in innate immunity. Despite their similar biological functions and protein architecture, comparative genome-wide analyses of NLRs and genes encoding NLR-like proteins suggest that plant and animal NLRs have independently arisen in evolution. Furthermore, the demonstration of interfamily transfer of plant NLR functions from their original species to phylogenetically distant species implies evolutionary conservation of the underlying immune principle across plant taxonomy. In this review we discuss plant NLR evolution and summarize recent insights into plant NLR-signaling mechanisms, which might constitute evolutionarily conserved NLR-mediated immune mechanisms.

The gene Sr33, an ortholog of barley Mla genes, encodes resistance to wheat stem rust race Ug99

Wheat stem rust, caused by the fungus Puccinia graminis f. sp. tritici, afflicts bread wheat (Triticum aestivum). New virulent races collectively referred to as "Ug99" have emerged, which threaten global wheat production. The wheat gene Sr33, introgressed from the wild relative Aegilops tauschii into bread wheat, confers resistance to diverse stem rust races, including the Ug99 race group. We cloned Sr33, which encodes a coiled-coil, nucleotide-binding, leucine-rich repeat protein. Sr33 is orthologous to the barley (Hordeum vulgare) Mla mildew resistance genes that confer resistance to Blumeria graminis f. sp. hordei. The wheat Sr33 gene functions independently of RAR1, SGT1, and HSP90 chaperones. Haplotype analysis from diverse collections of Ae. tauschii placed the origin of Sr33 resistance near the southern coast of the Caspian Sea.

Mosaic genome structure of the barley powdery mildew pathogen and conservation of transcriptional programs in divergent hosts

Barley powdery mildew, Blumeria graminis f. sp. hordei (Bgh), is an obligate biotrophic ascomycete fungal pathogen that can grow and reproduce only on living cells of wild or domesticated barley (Hordeum sp.). Domestication and deployment of resistant barley cultivars by humans selected for amplification of Bgh isolates with different virulence combinations. We sequenced the genomes of two European Bgh isolates, A6 and K1, for comparative analysis with the reference genome of isolate DH14. This revealed a mosaic genome structure consisting of large isolate-specific DNA blocks with either high or low SNP densities. Some of the highly polymorphic blocks likely accumulated SNPs for over 10,000 years, well before the domestication of barley. These isolate-specific blocks of alternating monomorphic and polymorphic regions imply an exceptionally large standing genetic variation in the Bgh population and might be generated and maintained by rare outbreeding and frequent clonal reproduction. RNA-sequencing experiments with isolates A6 and K1 during four early stages of compatible and incompatible interactions on leaves of partially immunocompromised Arabidopsis mutants revealed a conserved Bgh transcriptional program during pathogenesis compared with the natural host barley despite ~200 million years of reproductive isolation of these hosts. Transcripts encoding candidate-secreted effector proteins are massively induced in successive waves. A specific decrease in candidate-secreted effector protein transcript abundance in the incompatible interaction follows extensive transcriptional reprogramming of the host transcriptome and coincides with the onset of localized host cell death, suggesting a host-inducible defense mechanism that targets fungal effector secretion or production.

Haplotype variability and identification of new functional alleles at the Rdg2a leaf stripe resistance gene locus

The barley Rdg2a locus confers resistance to the leaf stripe pathogen Pyrenophora graminea and, in the barley genotype Thibaut, it is composed of a gene family with three highly similar paralogs. Only one member of the gene family (called as Rdg2a) encoding for a CC-NB-LRR protein is able to confer resistance to the leaf stripe isolate Dg2. To study the genome evolution and diversity at the Rdg2a locus, sequences spanning the Rdg2a gene were compared in two barley cultivars, Thibaut and Morex, respectively, resistant and susceptible to leaf stripe. An overall high level of sequence conservation interrupted by several rearrangements that included three main deletions was observed in the Morex contig. The main deletion of 13,692 bp was most likely derived from unequal crossing over between Rdg2a paralogs leading to the generation of a chimeric Morex rdg2a gene which was not associated to detectable level of resistance toward leaf stripe. PCR-based analyses of genic and intergenic regions at the Rdg2a locus in 29 H. vulgare lines and one H. vulgare ssp. spontaneum accession indicated large haplotype variability in the cultivated barley gene pool suggesting rapid and recent divergence at this locus. Barley genotypes showing the same haplotype as Thibaut at the Rdg2a locus were selected for a Rdg2a allele mining through allele re-sequencing and two lines with polymorphic nucleotides leading to amino acid changes in the CC-NB and LRR encoding domains, respectively, were identified. Analysis of nucleotide diversity of the Rdg2a alleles revealed that the polymorphic sites were subjected to positive selection. Moreover, strong positively selected sites were located in the LRR encoding domain suggesting that both positive selection and divergence at homologous loci are possibly representing the molecular mechanism for the generation of high diversity at the Rdg2a locus in the barley gene pool.

Barley MLA immune receptors directly interfere with antagonistically acting transcription factors to initiate disease resistance signaling

The nucleotide binding domain and Leucine-rich repeat (NLR)-containing proteins in plants and animals mediate pathogen sensing inside host cells and mount innate immune responses against microbial pathogens. The barley (Hordeum vulgare) mildew A (MLA) locus encodes coiled-coil (CC)-type NLRs mediating disease resistance against the powdery mildew pathogen Blumeria graminis. Here, we report direct interactions between MLA and two antagonistically acting transcription factors, MYB6 and WRKY1. The N-terminal CC signaling domain of MLA interacts with MYB6 to stimulate its DNA binding activity. MYB6 functions as a positive regulator of basal and MLA-mediated immunity responses to B. graminis. MYB6 DNA binding is antagonized by direct association with WRKY1 repressor, which in turn also interacts with the MLA CC domain. The activated form of full-length MLA10 receptor is needed to release MYB6 activator from WRKY1 repression and to stimulate MYB6-dependent gene expression. This implies that, while sequestered by the WRKY1 repressor in the presence of the resting immune receptor, MYB6 acts as an immediate and positive postactivation signaling component of the active state of MLA during transcriptional reprogramming for innate immune responses.

The differentially regulated genes TvQR1 and TvPirin of the parasitic plant Triphysaria exhibit distinctive natural allelic diversity

BACKGROUND

Plant parasitism represents an extraordinary interaction among flowering plants: parasitic plants use a specialized organ, the haustorium, to invade the host vascular system to deprive host plants of water and nutrients. Various compounds present in exudates of host plants trigger haustorium development. The two most effective haustorium inducing factors (HIFs) known for the parasitic plant Triphysaria versicolor (T. versicolor) are peonidin, an antioxidant flavonoid, and 2,6-dimethoxybenzoquinone (DMBQ), an oxidative stress agent. To date, two genes involved in haustorium initiation in T. versicolor have been identified: TvQR1, a quinone oxidoreductase that generates the active HIF from DMBQ, and TvPirin, a transcription co-factor that regulates several other DMBQ- responsive and -non-responsive genes. While the expression of these genes in response to DMBQ is well characterized, their expression in response to peonidin is not. In addition, the pattern of polymorphisms in these genes is unknown, even though nucleotide changes in TvQR1 and TvPirin may have contributed to the ability of T. versicolor to develop haustoria. To gain insights into these aspects, we investigated their transcriptional responses to HIFs and non-HIF and their natural nucleotide diversity.

RESULTS

Here we show that TvQR1 and TvPirin are transcriptionally upregulated by both DMBQ and peonidin in T. versicolor roots. Yet, while TvQR1 also responded to juglone, a non-HIF quinone with toxicity comparable to that of DMBQ, TvPirin did not. We further demonstrate that TvPirin encodes a protein shorter than the one previously reported. In the T. versicolor natural population of Northern California, TvQR1 exhibited remarkably higher molecular diversity and more recombination events than TvPirin, with the highest non-synonymous substitution rate in the substrate recognition and catalytic domain of the TvQR1 protein.

CONCLUSION

Our results suggest that TvQR1 and TvPirin have most likely evolved highly distinct roles for haustorium formation. Unlike TvPirin, TvQR1 might have been under diversifying selection to maintain a diverse collection of polymorphisms, which might be related to the recognition of an assortment of HIF and non-HIF quinones as substrates for successful haustorial establishment in a wide range of host plants.

Conservation of NLR-triggered immunity across plant lineages

The nucleotide-binding domain and leucine-rich repeat (NLR) family of plant receptors detects pathogen-derived molecules, designated effectors, inside host cells and mediates innate immune responses to pathogenic invaders. Genetic evidence revealed species-specific coevolution of many NLRs with effectors from host-adapted pathogens, suggesting that the specificity of these NLRs is restricted to the host or closely related plant species. However, we report that an NLR immune receptor (MLA1) from monocotyledonous barley is fully functional in partially immunocompromised dicotyledonous Arabidopsis thaliana against the barley powdery mildew fungus, Blumeria graminis f. sp. hordei. This implies ~200 million years of evolutionary conservation of the underlying immune mechanism. A time-course RNA-seq analysis in transgenic Arabidopsis lines detected sustained expression of a large MLA1-dependent gene cluster. This cluster is greatly enriched in genes known to respond to the fungal cell wall-derived microbe-associated molecular pattern chitin. The MLA1-dependent sustained transcript accumulation could define a conserved function of the nuclear pool of MLA1 detected in barley and Arabidopsis. We also found that MLA1-triggered immunity was fully retained in mutant plants that are simultaneously depleted of ethylene, jasmonic acid, and salicylic acid signaling. This points to the existence of an evolutionarily conserved and phytohormone-independent MLA1-mediated resistance mechanism. This also suggests a conserved mechanism for internalization of B. graminis f. sp. hordei effectors into host cells of flowering plants. Furthermore, the deduced connectivity of the NLR to multiple branches of immune signaling pathways likely confers increased robustness against pathogen effector-mediated interception of host immune signaling and could have contributed to the evolutionary preservation of the immune mechanism.

A physical, genetic and functional sequence assembly of the barley genome

Barley (Hordeum vulgare L.) is among the world's earliest domesticated and most important crop plants. It is diploid with a large haploid genome of 5.1 gigabases (Gb). Here we present an integrated and ordered physical, genetic and functional sequence resource that describes the barley gene-space in a structured whole-genome context. We developed a physical map of 4.98 Gb, with more than 3.90 Gb anchored to a high-resolution genetic map. Projecting a deep whole-genome shotgun assembly, complementary DNA and deep RNA sequence data onto this framework supports 79,379 transcript clusters, including 26,159 'high-confidence' genes with homology support from other plant genomes. Abundant alternative splicing, premature termination codons and novel transcriptionally active regions suggest that post-transcriptional processing forms an important regulatory layer. Survey sequences from diverse accessions reveal a landscape of extensive single-nucleotide variation. Our data provide a platform for both genome-assisted research and enabling contemporary crop improvement.

How to build a pathogen detector: structural basis of NB-LRR function

Many plant disease resistance (R) proteins belong to the family of nucleotide-binding-leucine rich repeat (NB-LRR) proteins. NB-LRRs mediate recognition of pathogen-derived effector molecules and subsequently activate host defence. Their multi-domain structure allows these pathogen detectors to simultaneously act as sensor, switch and response factor. Structure-function analyses and the recent elucidation of the 3D structures of subdomains have provided new insight in how these different functions are combined and what the contribution is of the individual subdomains. Besides interdomain contacts, interactions with chaperones, the proteasome and effector baits are required to keep NB-LRRs in a signalling-competent, yet auto-inhibited state. In this review we explore operational models of NB-LRR functioning based on recent advances in understanding their structure.

Molecular characterization and comparative sequence analysis of defense-related gene, Oryza rufipogon receptor-like protein kinase 1

Many of the plant leucine rich repeat receptor-like kinases (LRR-RLKs) have been found to regulate signaling during plant defense processes. In this study, we selected and sequenced an LRR-RLK gene, designated as Oryza rufipogon receptor-like protein kinase 1 (OrufRPK1), located within yield QTL yld1.1 from the wild rice Oryza rufipogon (accession IRGC105491). A 2055 bp coding region and two exons were identified. Southern blotting determined OrufRPK1 to be a single copy gene. Sequence comparison with cultivated rice orthologs (OsI219RPK1, OsI9311RPK1 and OsJNipponRPK1, respectively derived from O. sativa ssp. indica cv. MR219, O. sativa ssp. indica cv. 9311 and O. sativa ssp. japonica cv. Nipponbare) revealed the presence of 12 single nucleotide polymorphisms (SNPs) with five non-synonymous substitutions, and 23 insertion/deletion sites. The biological role of the OrufRPK1 as a defense related LRR-RLK is proposed on the basis of cDNA sequence characterization, domain subfamily classification, structural prediction of extra cellular domains, cluster analysis and comparative gene expression.

Structure-function analysis of barley NLR immune receptor MLA10 reveals its cell compartment specific activity in cell death and disease resistance

Plant intracellular immune receptors comprise a large number of multi-domain proteins resembling animal NOD-like receptors (NLRs). Plant NLRs typically recognize isolate-specific pathogen-derived effectors, encoded by avirulence (AVR) genes, and trigger defense responses often associated with localized host cell death. The barley MLA gene is polymorphic in nature and encodes NLRs of the coiled-coil (CC)-NB-LRR type that each detects a cognate isolate-specific effector of the barley powdery mildew fungus. We report the systematic analyses of MLA10 activity in disease resistance and cell death signaling in barley and Nicotiana benthamiana. MLA10 CC domain-triggered cell death is regulated by highly conserved motifs in the CC and the NB-ARC domains and by the C-terminal LRR of the receptor. Enforced MLA10 subcellular localization, by tagging with a nuclear localization sequence (NLS) or a nuclear export sequence (NES), shows that MLA10 activity in cell death signaling is suppressed in the nucleus but enhanced in the cytoplasm. By contrast, nuclear localized MLA10 is sufficient to mediate disease resistance against powdery mildew fungus. MLA10 retention in the cytoplasm was achieved through attachment of a glucocorticoid receptor hormone-binding domain (GR), by which we reinforced the role of cytoplasmic MLA10 in cell death signaling. Together with our data showing an essential and sufficient nuclear MLA10 activity in disease resistance, this suggests a bifurcation of MLA10-triggered cell death and disease resistance signaling in a compartment-dependent manner.

Differential accumulation of host mRNAs on polyribosomes during obligate pathogen-plant interactions

Plant pathogens elicit dramatic changes in the expression of host genes during both compatible and incompatible interactions. Gene expression profiling studies of plant-pathogen interactions have only considered messenger RNAs (mRNAs) present in total RNA, which contains subpopulations of actively translated mRNAs associated with polyribosomes (polysomes) and non-translated mRNAs that are not associated with polysomes. The goal of this study was to enhance previous gene expression analyses by identifying host mRNAs that become differentially associated with polysomes following pathogen inoculation. Total and polysomal RNA were extracted from barley (Hordeum vulgare) plants at 32 h after inoculation with Blumeria graminis f. sp. hordei, and Arabidopsis thaliana plants at 10 days after inoculation with Turnip mosaic virus. Gene expression profiles were obtained for each pathosystem, which represent diverse plant host-obligate pathogen interactions. Using this approach, host mRNAs were identified that were differentially associated with polysomes in response to pathogen treatment. Approximately 18% and 26% of mRNAs represented by probe sets on the Affymetrix Barley1 and Arabidopsis ATH1 GeneChips, respectively, differentially accumulated in the two populations in one or more combinations of treatment and genotype. Gene ontology analysis of mRNAs sharing the same pattern of accumulation in total and polysomal RNA identified gene sets that contained a significant number of functionally related annotations, suggesting both transcript accumulation and recruitment to polyribosomes are coordinately regulated in these systems.

Regulation of basal resistance by a powdery mildew-induced cysteine-rich receptor-like protein kinase in barley

The receptor-like protein kinases (RLKs) constitute a large and diverse group of proteins controlling numerous plant physiological processes, including development, hormone perception and stress responses. The cysteine-rich RLKs (CRKs) represent a prominent subfamily of transmembrane-anchored RLKs. We have identified a putative barley (Hordeum vulgare) CRK gene family member, designated HvCRK1. The mature putative protein comprises 645 amino acids, and includes a putative receptor domain containing two characteristic 'domain 26 of unknown function' (duf26) domains in the N-terminal region, followed by a rather short 17-amino-acid transmembrane domain, which includes an AAA motif, two features characteristic of endoplasmic reticulum (ER)-targeted proteins and, finally, a characteristic putative protein kinase domain in the C-terminus. The HvCRK1 transcript was isolated from leaves inoculated with the biotrophic fungal pathogen Blumeria graminis f.sp. hordei (Bgh). HvCRK1 transcripts were observed to accumulate transiently following Bgh inoculation of susceptible barley. Transient silencing of HvCRK1 expression in bombarded epidermal cells led to enhanced resistance to Bgh, but did not affect R-gene-mediated resistance. Silencing of HvCRK1 phenocopied the effective penetration resistance found in mlo-resistant barley plants, and the possible link between HvCRK1 and MLO was substantiated by the fact that HvCRK1 induction on Bgh inoculation was dependent on Mlo. Finally, using both experimental and in silico approaches, we demonstrated that HvCRK1 localizes to the ER of barley cells. The negative effect on basal resistance against Bgh and the functional aspects of MLO- and ER-localized HvCRK1 signalling on Bgh inoculation are discussed.

Selective sweep at the Rpv3 locus during grapevine breeding for downy mildew resistance

The Rpv3 locus is a major determinant of downy mildew resistance in grapevine (Vitis spp.). A selective sweep at this locus was revealed by the DNA genotyping of 580 grapevines, which include a highly diverse set of 265 European varieties that predated the spread of North American mildews, 82 accessions of wild species, and 233 registered breeding lines with North American ancestry produced in the past 150 years. Artificial hybridisation and subsequent phenotypic selection favoured a few Rpv3 haplotypes that were introgressed from wild vines and retained in released varieties. Seven conserved haplotypes in five descent groups of resistant varieties were traced back to their founders: (1) 'Munson', a cross between two of Hermann Jaeger's selections of V. rupestris and V. lincecumii made in the early 1880s in Missouri, (2) V. rupestris 'Ganzin', first utilised for breeding in 1879 by Victor Ganzin in France, (3) 'Noah', selected in 1869 from intermingled accessions of V. riparia and V. labrusca by Otto Wasserzieher in Illinois, (4) 'Bayard', a V. rupestris × V. labrusca offspring generated in 1882 by George Couderc in France, and (5) a wild form closely related to V. rupestris accessions in the Midwestern United States and introgressed into 'Seibel 4614' in the 1880s by Albert Seibel in France. Persistence of these Rpv3 haplotypes across many of the varieties generated by human intervention indicates that a handful of vines with prominent resistance have laid the foundation for modern grape breeding. A rampant hot spot of NB-LRR genes at the Rpv3 locus has provided a distinctive advantage for the adaptation of native North American grapevines to withstand downy mildew. The coexistence of multiple resistance alleles or paralogues in the same chromosomal region but in different haplotypes counteracts efforts to pyramidise them in a diploid individual via conventional breeding.

Large-scale data integration reveals colocalization of gene functional groups with meta-QTL for multiple disease resistance in barley

Race-nonspecific and durable resistance of plant genotypes to major pathogens is highly relevant for yield stability and sustainable crop production but difficult to handle in practice due to its polygenic inheritance by quantitative trait loci (QTL). As far as the underlying genes are concerned, very little is currently known in the most important crop plants such as the cereals. Here, we integrated publicly available data for barley (Hordeum vulgare subsp. vulgare) in order to detect the most important genomic regions for QTL-mediated resistance to a number of fungal pathogens and localize specific functional groups of genes within these regions. This identified 20 meta-QTL, including eight hot spots for resistance to multiple diseases that were distributed over all chromosomes. At least one meta-QTL region for resistance to the powdery mildew fungus Blumeria graminis was found to be co-linear between barley and wheat, suggesting partial evolutionary conservation. Large-scale genetic mapping revealed that functional groups of barley genes involved in secretory processes and cell-wall reinforcement were significantly over-represented within QTL for resistance to powdery mildew. Overall, the results demonstrate added value resulting from large-scale genetic and genomic data integration and may inform genomic-selection procedures for race-nonspecific and durable disease resistance in barley.

NLR functions in plant and animal immune systems: so far and yet so close

In plants and animals, the NLR family of receptors perceives non-self and modified-self molecules inside host cells and mediates innate immune responses to microbial pathogens. Despite their similar biological functions and protein architecture, animal NLRs are normally activated by conserved microbe- or damage-associated molecular patterns, whereas plant NLRs typically detect strain-specific pathogen effectors. Plant NLRs recognize either the effector structure or effector-mediated modifications of host proteins. The latter indirect mechanism for the perception of non-self, as well as the within-species diversification of plant NLRs, maximize the capacity to recognize non-self through the use of a finite number of innate immunoreceptors. We discuss recent insights into NLR activation, signal initiation through the homotypic association of N-terminal domains and subcellular receptor dynamics in plants and compare those with NLR functions in animals.