Cell-autonomous expression of barley Mla1 confers race-specific resistance to the powdery mildew fungus via a Rar1-independent signaling pathway

Molecular Basis of Plant-Pathogen Interactions in the Agricultural Context

Biotic stressors pose significant threats to crop yield, jeopardizing food security and resulting in losses of over USD 220 billion per year by the agriculture industry. Plants activate innate defense mechanisms upon pathogen perception and invasion. The plant immune response comprises numerous concerted steps, including the recognition of invading pathogens, signal transduction, and activation of defensive pathways. However, pathogens have evolved various structures to evade plant immunity. Given these facts, genetic improvements to plants are required for sustainable disease management to ensure global food security. Advanced genetic technologies have offered new opportunities to revolutionize and boost plant disease resistance against devastating pathogens. Furthermore, targeting susceptibility (S) genes, such as OsERF922 and BnWRKY70, through CRISPR methodologies offers novel avenues for disrupting the molecular compatibility of pathogens and for introducing durable resistance against them in plants. Here, we provide a critical overview of advances in understanding disease resistance mechanisms. The review also critically examines management strategies under challenging environmental conditions and R-gene-based plant genome-engineering systems intending to enhance plant responses against emerging pathogens. This work underscores the transformative potential of modern genetic engineering practices in revolutionizing plant health and crop disease management while emphasizing the importance of responsible application to ensure sustainable and resilient agricultural systems.

Genome-wide association study and haplotype analysis reveal novel candidate genes for resistance to powdery mildew in soybean

Powdery mildew disease (PMD) is caused by the obligate biotrophic fungus Microsphaera diffusa Cooke & Peck (M. diffusa) and results in significant yield losses in soybean (Glycine max (L.) Merr.) crops. By identifying disease-resistant genes and breeding soybean accessions with enhanced resistance, we can effectively mitigate the detrimental impact of PMD on soybeans. We analyzed PMD resistance in a diversity panel of 315 soybean accessions in two locations over 3 years, and candidate genes associated with PMD resistance were identified through genome-wide association studies (GWAS), haplotype analysis, qRT-PCR, and EMS mutant analysis. Based on the GWAS approach, we identified a region on chromosome 16 (Chr16) in which 21 genes form a gene cluster that is highly correlated with PMD resistance. In order to validate and refine these findings, we conducted haplotype analysis of 21 candidate genes and indicated there are single nucleotide polymorphisms (SNPs) and insertion-deletions (InDels) variations of Glyma.16G214000, Glyma.16G214200, Glyma.16G215100 and Glyma.16G215300 within the coding and promoter regions that exhibit a strong association with resistance against PMD. Subsequent structural analysis of candidate genes within this cluster revealed that in 315 accessions, the majority of accessions exhibited resistance to PMD when Glyma.16G214300, Glyma.16G214800 and Glyma.16G215000 were complete; however, they demonstrated susceptibility to PMD when these genes were incomplete. Quantitative real-time PCR assays (qRT-PCR) of possible candidate genes showed that 14 candidate genes (Glyma.16G213700, Glyma.16G213800, Glyma.16G213900, Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214500, Glyma.16G214585, Glyma.16G214669, Glyma.16G214700, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300) were involved in PMD resistance. Finally, we evaluated the PMD resistance of mutant lines from the Williams 82 EMS mutations library, which revealed that mutants of Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300, exhibited sensitivity to PMD. Combined with the analysis results of GWAS, haplotypes, qRT-PCR and mutants, the genes Glyma.16G214000, Glyma.16G214200, Glyma.16G214300, Glyma.16G214800, Glyma.16G215000, Glyma.16G215100 and Glyma.16G215300 were identified as highly correlated with PMD resistance. The candidate genes identified above are all NLR family genes, and these discoveries deepen our understanding of the molecular basis of PMD resistance in soybeans and will be useful for guiding breeding strategies.

Genome-wide association mapping for seedling and adult resistance to powdery mildew in barley

KEY MESSAGE

Two new major QTL were identified for powdery mildew resistance. We confirmed that the QTL on 7HS contributed mainly to the adult-plant resistance, while another one on chromosome arm 1HS made a significant contribution to the seedling resistance. Powdery mildew (PM), caused by Blumeria hordei, can occur at all post emergent stages of barley and constantly threatens crop production. To identify more genes for effective resistance to powdery mildew for use in breeding programs, 696 barley accessions collected from different regions of the world were evaluated for PM resistance at seedling and adult growth stages in three different states of Australia. These barley accessions were genotyped using DArTSeq with over 18,000 markers for a genome-wide association study (GWAS). Using the FarmCPU model, 54 markers showed significant associations with PM resistance scored at the seedling and adult-plant stages in different states of Australia. Another 40 markers showed tentative associations (LOD > 4.0) with resistance. These markers are distributed across all seven barley chromosomes. Most of them were grouped into eleven QTL regions, coinciding with the locations of most of the reported resistance genes. Two major MTAs were identified on chromosome arms 3HS and 5HL, with one on 3HS contributing to adult plant resistance and the one on 5HL to both seedling and adult plant resistance. An MTA on 7HS contributed mainly to the adult-plant resistance, while another one on chromosome arm 1HS made a significant contribution to the seedling resistance.

Barley MLA3 recognizes the host-specificity effector Pwl2 from Magnaporthe oryzae

Plant nucleotide-binding leucine-rich repeat (NLRs) immune receptors directly or indirectly recognize pathogen-secreted effector molecules to initiate plant defense. Recognition of multiple pathogens by a single NLR is rare and usually occurs via monitoring for changes to host proteins; few characterized NLRs have been shown to recognize multiple effectors. The barley (Hordeum vulgare) NLR gene Mildew locus a (Mla) has undergone functional diversification, and the proteins encoded by different Mla alleles recognize host-adapted isolates of barley powdery mildew (Blumeria graminis f. sp. hordei [Bgh]). Here, we show that Mla3 also confers resistance to the rice blast fungus Magnaporthe oryzae in a dosage-dependent manner. Using a forward genetic screen, we discovered that the recognized effector from M. oryzae is Pathogenicity toward Weeping Lovegrass 2 (Pwl2), a host range determinant factor that prevents M. oryzae from infecting weeping lovegrass (Eragrostis curvula). Mla3 has therefore convergently evolved the capacity to recognize effectors from diverse pathogens.

Identification and genomic characterization of major effect bacterial blight resistance locus (BB-13) in Upland cotton (Gossypium hirsutum L.)

Identification and genomic characterization of major resistance locus against cotton bacterial blight (CBB) using GWAS and linkage mapping to enable genomics-based development of durable CBB resistance and gene discovery in cotton. Cotton bacterial leaf blight (CBB), caused by Xanthomonas citri subsp. malvacearum (Xcm), has periodically been a damaging disease in the USA. Identification and deployment of genetic resistance in cotton cultivars is the most economical and efficient means of reducing crop losses due to CBB. In the current study, genome-wide association study (GWAS) of CBB resistance using an elite diversity panel of 380 accessions, genotyped with the cotton single nucleotide polymorphism (SNP) 63 K array, and phenotyped with race-18 of CBB, localized the CBB resistance to a 2.01-Mb region in the long arm of chromosome D02. Molecular genetic mapping using an F₆ recombinant inbred line (RIL) population showed the CBB resistance in cultivar Arkot 8102 was controlled by a single locus (BB-13). The BB-13 locus was mapped within the 0.95-cM interval near the telomeric region in the long arm of chromosome D02. Flanking SNP markers, i04890Gh and i04907Gh of the BB-13 locus, identified from the combined linkage analysis and GWAS, targeted it to a 371-Kb genomic region. Candidate gene analysis identified thirty putative gene sequences in the targeted genomic region. Nine of these putative genes and two NBS-LRR genes adjacent to the targeted region were putatively involved in plant disease resistance and are possible candidate genes for BB-13 locus. Genetic mapping and genomic targeting of the BB13 locus in the current study will help in cloning the CBB-resistant gene and establishing the molecular genetic architecture of the BB-13 locus towards developing durable resistance to CBB in cotton.

Phenotypic variation and quantitative trait loci for resistance to southern anthracnose and clover rot in red clover

High variability for and candidate loci associated with resistance to southern anthracnose and clover rot in a worldwide collection of red clover provide a first basis for genomics-assisted breeding. Red clover (Trifolium pratense L.) is an important forage legume of temperate regions, particularly valued for its high yield potential and its high forage quality. Despite substantial breeding progress during the last decades, continuous improvement of cultivars is crucial to ensure yield stability in view of newly emerging diseases or changing climatic conditions. The high amount of genetic diversity present in red clover ecotypes, landraces, and cultivars provides an invaluable, but often unexploited resource for the improvement of key traits such as yield, quality, and resistance to biotic and abiotic stresses. A collection of 397 red clover accessions was genotyped using a pooled genotyping-by-sequencing approach with 200 plants per accession. Resistance to the two most pertinent diseases in red clover production, southern anthracnose caused by Colletotrichum trifolii, and clover rot caused by Sclerotinia trifoliorum, was assessed using spray inoculation. The mean survival rate for southern anthracnose was 22.9% and the mean resistance index for clover rot was 34.0%. Genome-wide association analysis revealed several loci significantly associated with resistance to southern anthracnose and clover rot. Most of these loci are in coding regions. One quantitative trait locus (QTL) on chromosome 1 explained 16.8% of the variation in resistance to southern anthracnose. For clover rot resistance we found eight QTL, explaining together 80.2% of the total phenotypic variation. The SNPs associated with these QTL provide a promising resource for marker-assisted selection in existing breeding programs, facilitating the development of novel cultivars with increased resistance against two devastating fungal diseases of red clover.

An interolog-based barley interactome as an integration framework for immune signaling

The barley MLA nucleotide-binding leucine-rich-repeat (NLR) receptor and its orthologs confer recognition specificity to many fungal diseases, including powdery mildew, stem-, and stripe rust. We used interolog inference to construct a barley protein interactome (Hordeum vulgare predicted interactome, HvInt) comprising 66,133 edges and 7,181 nodes, as a foundation to explore signaling networks associated with MLA. HvInt was compared with the experimentally validated Arabidopsis interactome of 11,253 proteins and 73,960 interactions, verifying that the 2 networks share scale-free properties, including a power-law distribution and small-world network. Then, by successive layering of defense-specific "omics" datasets, HvInt was customized to model cellular response to powdery mildew infection. Integration of HvInt with expression quantitative trait loci (eQTL) enabled us to infer disease modules and responses associated with fungal penetration and haustorial development. Next, using HvInt and infection-time-course RNA sequencing of immune signaling mutants, we assembled resistant and susceptible subnetworks. The resulting differentially coexpressed (resistant - susceptible) interactome is essential to barley immunity, facilitates the flow of signaling pathways and is linked to mildew resistance locus a (Mla) through trans eQTL associations. Lastly, we anchored HvInt with new and previously identified interactors of the MLA coiled coli + nucleotide-binding domains and extended these to additional MLA alleles, orthologs, and NLR outgroups to predict receptor localization and conservation of signaling response. These results link genomic, transcriptomic, and physical interactions during MLA-specified immunity.

The barley immune receptor Mla recognizes multiple pathogens and contributes to host range dynamics

Crop losses caused by plant pathogens are a primary threat to stable food production. Stripe rust (Puccinia striiformis) is a fungal pathogen of cereal crops that causes significant, persistent yield loss. Stripe rust exhibits host species specificity, with lineages that have adapted to infect wheat and barley. While wheat stripe rust and barley stripe rust are commonly restricted to their corresponding hosts, the genes underlying this host specificity remain unknown. Here, we show that three resistance genes, Rps6, Rps7, and Rps8, contribute to immunity in barley to wheat stripe rust. Rps7 cosegregates with barley powdery mildew resistance at the Mla locus. Using transgenic complementation of different Mla alleles, we confirm allele-specific recognition of wheat stripe rust by Mla. Our results show that major resistance genes contribute to the host species specificity of wheat stripe rust on barley and that a shared genetic architecture underlies resistance to the adapted pathogen barley powdery mildew and non-adapted pathogen wheat stripe rust.

RACE1, a Japanese Blumeria graminis f. sp. hordei isolate, is capable of overcoming partially mlo-mediated penetration resistance in barley in an allele-specific manner

Loss-of-function mutation of the MILDEW RESISTANCE LOCUS O (Mlo) gene confers durable and broad-spectrum resistance to powdery mildew fungi in various plants, including barley. In combination with the intracellular nucleotide-binding domain and leucine-rich repeat receptor (NLR) genes, which confer the race-specific resistance, the mlo alleles have long been used in barley breeding as genetic resources that confer robust non-race-specific resistance. However, a Japanese Blumeria graminis f. sp. hordei isolate, RACE1, has been reported to have the potential to overcome partially the mlo-mediated penetration resistance, although this is yet uncertain because the putative effects of NLR genes in the tested accessions have not been ruled out. In this study, we examined the reproducibility of the earlier report and found that the infectious ability of RACE1, which partially overcomes the mlo-mediated resistance, is only exerted in the absence of NLR genes recognizing RACE1. Furthermore, using the transient-induced gene silencing technique, we demonstrated that RACE1 can partially overcome the resistance in the host cells with suppressed MLO expression but not in plants possessing the null mutant allele mlo-5.

Diploid chromosome-scale assembly of the Muscadinia rotundifolia genome supports chromosome fusion and disease resistance gene expansion during Vitis and Muscadinia divergence

Muscadinia rotundifolia, the muscadine grape, has been cultivated for centuries in the southeastern United States. M. rotundifolia is resistant to many of the pathogens that detrimentally affect Vitis vinifera, the grape species commonly used for winemaking. For this reason, M. rotundifolia is a valuable genetic resource for breeding. Single-molecule real-time reads were combined with optical maps to reconstruct the two haplotypes of each of the 20 M. rotundifolia cv. Trayshed chromosomes. The completeness and accuracy of the assembly were confirmed using a high-density linkage map. Protein-coding genes were annotated using an integrated and comprehensive approach. This included using full-length cDNA sequencing (Iso-Seq) to improve gene structure and hypothetical spliced variant predictions. Our data strongly support that Muscadinia chromosomes 7 and 20 are fused in Vitis and pinpoint the location of the fusion in Cabernet Sauvignon and PN40024 chromosome 7. Disease-related gene numbers in Trayshed and Cabernet Sauvignon were similar, but their clustering locations were different. A dramatic expansion of the Toll/Interleukin-1 Receptor-like Nucleotide-Binding Site Leucine-Rich Repeat (TIR-NBS-LRR) class was detected on Trayshed chromosome 12 at the Resistance to Uncinula necator 1 (RUN1)/Resistance to Plasmopara viticola 1 (RPV1) locus, which confers strong dominant resistance to powdery and downy mildews. A genome browser, annotation, and Blast tool for Trayshed are available at www.grapegenomics.com.

Identification and mapping of a novel resistance gene to the rice pathogen, Cercospora janseana

The genetic architecture of resistance to Cercospora janseana was examined, and a single resistance locus was identified. A SNP marker was identified and validated for utilization in U.S. breeding germplasm Cercospora janseana (Racib.) is a fungal pathogen that causes narrow brown leaf spot (NBLS) in rice. Although NBLS is a major disease in the southern United States and variation in resistance among U.S. rice germplasm exists, little is known about the genetic architecture underlying the trait. In this study, a recombinant inbred line population was evaluated for NBLS resistance under natural disease infestation in the field across three years. A single, large-effect QTL, CRSP-2.1, was identified that explained 81.4% of the phenotypic variation. The QTL was defined to a 532 kb physical interval and 13 single nucleotide polymorphisms (SNPs) were identified across the region to characterize the haplotype diversity present in U.S. rice germplasm. A panel of 387 U.S. rice germplasm was genotyped with the 13 haplotype SNPs and phenotyped over two years for NBLS resistance. Fourteen haplotypes were identified, with six haplotypes accounting for 94% of the panel. The susceptible haplotype from the RIL population was the only susceptible haplotype observed in the U.S. germplasm. A single SNP was identified that distinguished the susceptible haplotype from all resistant haplotypes, explaining 52.7% of the phenotypic variation for NBLS resistance. Pedigree analysis and haplotype characterization of historical germplasm demonstrated that the susceptible haplotype was introduced into Southern U.S. germplasm through the California line L-202 into the Louisiana variety Cypress. Cypress was extensively used as a parent over the last 25 years, resulting in the susceptible CRSP-2.1 allele increasing in frequency from zero to 44% in the modern U.S. germplasm panel.

Recognition and defence of plant-infecting fungal pathogens

Attempted infections of plants with fungi result in diverse outcomes ranging from symptom-less resistance to severe disease and even death of infected plants. The deleterious effect on crop yield have led to intense focus on the cellular and molecular mechanisms that explain the difference between resistance and susceptibility. This research has uncovered plant resistance or susceptibility genes that explain either dominant or recessive inheritance of plant resistance with many of them coding for receptors that recognize pathogen invasion. Approaches based on cell biology and phytochemistry have contributed to identifying factors that halt an invading fungal pathogen from further invasion into or between plant cells. Plant chemical defence compounds, antifungal proteins and structural reinforcement of cell walls appear to slow down fungal growth or even prevent fungal penetration in resistant plants. Additionally, the hypersensitive response, in which a few cells undergo a strong local immune reaction, including programmed cell death at the site of infection, stops in particular biotrophic fungi from spreading into surrounding tissue. In this review, we give a general overview of plant recognition and defence of fungal parasites tracing back to the early 20th century with a special focus on Triticeae and on the progress that was made in the last 30 years.

You Had Me at "MAGIC"!: Four Barley MAGIC Populations Reveal Novel Resistance QTL for Powdery Mildew

Blumeria graminis f. sp. hordei (Bgh), the causal agent of barley powdery mildew (PM), is one of the most important barley leaf diseases and is prevalent in most barley growing regions. Infection decreases grain quality and yields on average by 30%. Multi-parent advanced generation inter-cross (MAGIC) populations combine the advantages of bi-parental and association panels and offer the opportunity to incorporate exotic alleles into adapted material. Here, four barley MAGIC populations consisting of six to eight founders were tested for PM resistance in field trials in Denmark. Principle component and STRUCTURE analysis showed the populations were unstructured and genome-wide linkage disequilibrium (LD) decay varied between 14 and 38 Mbp. Genome-wide association studies (GWAS) identified 11 regions associated with PM resistance located on chromosomes 1H, 2H, 3H, 4H, 5H and 7H, of which three regions are putatively novel resistance quantitative trait locus/loci (QTL). For all regions high-confidence candidate genes were identified that are predicted to be involved in pathogen defense. Haplotype analysis of the significant SNPs revealed new allele combinations not present in the founders and associated with high resistance levels.

Apoplastic Fluid Preparation from Arabidopsis thaliana Leaves Upon Interaction with a Nonadapted Powdery Mildew Pathogen

Proteins in the extracellular space (apoplast) play a crucial role at the interface between plant cells and their proximal environment. Consequently, it is not surprising that plants actively control the apoplastic proteomic profile in response to biotic and abiotic cues. Comparative quantitative proteomics of plant apoplastic fluids is therefore of general interest in plant physiology. We here describe an efficient method to isolate apoplastic fluids from Arabidopsis thaliana leaves inoculated with a nonadapted powdery mildew pathogen.

Transcriptome Analysis Identifies Candidate Genes and Functional Pathways Controlling the Response of Two Contrasting Barley Varieties to Powdery Mildew Infection

Powdery mildew caused by Blumeria graminis f. sp. hordei (Bgh) is one of the most serious diseases in barley. The numerous barley varieties across China provide valuable genetic resources to screen the resistant germplasm and to discover the primary genes of resistance to powdery mildew. In this study, Chinese barley variety Feng 7 was identified as a highly resistant genotype which limited Bgh colonization by cell apoptosis using leaf staining assay, while another variety Hua 30 showed high susceptibility. The performance of high resistance to Bgh in F₁ plants from the two varieties suggested dominant gene(s) controlled the resistance to powdery mildew in Feng 7. To understand the host transcriptional response to Bgh infection, these two barley varieties Feng 7 and Hua 30 were inoculated with Bgh, and their transcriptional profiling using RNA sequencing (RNA-seq) at four time points (12 h post-inoculation (hpi), 24 hpi, 48 hpi, and 72 hpi) were compared. 4318 differentially expressed genes (DEGs), including 2244 upregulated and 2074 downregulated genes, were detected in Feng 7, compared with Hua 30 at 12 hpi. 4907 DEGs (2488 upregulated and 2419 downregulated) were detected at 24 hpi. 4758 DEGs (2295 upregulated and 2463 downregulated) were detected at 48 hpi. 3817 DEGs (2036 upregulated and 1781 downregulated) were detected at 72 hpi. The results showed the number of DEGs between two varieties peaked at 24 hpi (for the upregulated) or 48 hpi (for the downregulated), which is matched with the processing of Bgh infection. In addition, the number of upregulated DEGs involved in the functional pathways of plant defense (mitogen-activated protein kinase (MAPK) pathway and plant hormone signal transduction) is elevated remarkably at 24 hpi. Six candidate genes (PR13, glutaredoxin, alcohol dehydrogenase, and cytochrome P450) were identified in Feng 7. All of them present continuous expression at higher levels upon Bgh infection, compared with the performance in Hua 30, which revealed the potential contribution to Feng 7 mediate resistance to Bgh. In conclusion, the candidate genes and relevant pathways provided key information towards understanding the defense of barley to Bgh attack and the molecular mechanisms of different genetic resistance to powdery mildew.

Small RNA discovery in the interaction between barley and the powdery mildew pathogen

Background

Plants encounter pathogenic and non-pathogenic microorganisms on a nearly constant basis. Small RNAs such as siRNAs and miRNAs/milRNAs influence pathogen virulence and host defense responses. We exploited the biotrophic interaction between the powdery mildew fungus, Blumeria graminis f. sp. hordei (Bgh), and its diploid host plant, barley (Hordeum vulgare) to explore fungal and plant sRNAs expressed during Bgh infection of barley leaf epidermal cells.

Results

RNA was isolated from four fast-neutron immune-signaling mutants and their progenitor over a time course representing key stages of Bgh infection, including appressorium formation, penetration of epidermal cells, and development of haustorial feeding structures. The Cereal Introduction (CI) 16151 progenitor carries the resistance allele Mla6, while Bgh isolate 5874 harbors the AVR_a6 avirulence effector, resulting in an incompatible interaction. Parallel Analysis of RNA Ends (PARE) was used to verify sRNAs with likely transcript targets in both barley and Bgh. Bgh sRNAs are predicted to regulate effectors, metabolic genes, and translation-related genes. Barley sRNAs are predicted to influence the accumulation of transcripts that encode auxin response factors, NAC transcription factors, homeodomain transcription factors, and several splicing factors. We also identified phasing small interfering RNAs (phasiRNAs) in barley that overlap transcripts that encode receptor-like kinases (RLKs) and nucleotide-binding, leucine-rich domain proteins (NLRs).

Conclusions

These data suggest that Bgh sRNAs regulate gene expression in metabolism, translation-related, and pathogen effectors. PARE-validated targets of predicted Bgh milRNAs include both EKA (effectors homologous to AVR_k1 and AVR_a10) and CSEP (candidate secreted effector protein) families. We also identified barley phasiRNAs and miRNAs in response to Bgh infection. These include phasiRNA loci that overlap with a significant proportion of receptor-like kinases, suggesting an additional sRNA control mechanism may be active in barley leaves as opposed to predominant R-gene phasiRNA overlap in many eudicots. In addition, we identified conserved miRNAs, novel miRNA candidates, and barley genome mapped sRNAs that have PARE validated transcript targets in barley. The miRNA target transcripts are enriched in transcription factors, signaling-related proteins, and photosynthesis-related proteins. Together these results suggest both barley and Bgh control metabolism and infection-related responses via the specific accumulation and targeting of genes via sRNAs.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5947-z) contains supplementary material, which is available to authorized users.

Discovering coherency of specific gene expression and optical reflectance properties of barley genotypes differing for resistance reactions against powdery mildew

Hyperspectral imaging has proved its potential for evaluating complex plant-pathogen interactions. However, a closer link of the spectral signatures and genotypic characteristics remains elusive. Here, we show relation between gene expression profiles and specific wavebands from reflectance during three barley-powdery mildew interactions. Significant synergistic effects between the hyperspectral signal and the corresponding gene activities has been shown using the linear discriminant analysis (LDA). Combining the data sets of hyperspectral signatures and gene expression profiles allowed a more precise differentiation of the three investigated barley-Bgh interactions independent from the time after inoculation. This shows significant synergistic effects between the hyperspectral signal and the corresponding gene activities. To analyze this coherency between spectral reflectance and seven different gene expression profiles, relevant wavelength bands and reflectance intensities for each gene were computed using the Relief algorithm. Instancing, xylanase activity was indicated by relevant wavelengths around 710 nm, which are characterized by leaf and cell structures. HvRuBisCO activity underlines relevant wavebands in the green and red range, elucidating the coherency of RuBisCO to the photosynthesis apparatus and in the NIR range due to the influence of RuBisCO on barley leaf cell development. These findings provide the first insights to links between gene expression and spectral reflectance that can be used for an efficient non-invasive phenotyping of plant resistance and enables new insights into plant-pathogen interactions.

Multiple pairs of allelic MLA immune receptor-powdery mildew AVRA effectors argue for a direct recognition mechanism

Nucleotide-binding domain and leucine-rich repeat (NLR)-containing proteins in plants and animals mediate intracellular pathogen sensing. Plant NLRs typically detect strain-specific pathogen effectors and trigger immune responses often linked to localized host cell death. The barley Mla disease resistance locus has undergone extensive functional diversification in the host population and encodes numerous allelic NLRs each detecting a matching isolate-specific avirulence effector (AVRA) of the fungal pathogen Blumeria graminis f. sp. hordei (Bgh). We report here the isolation of Bgh AVRa7, AVRa9, AVRa10, and AVRa22, which encode small secreted proteins recognized by allelic MLA7, MLA9, MLA10, and MLA22 receptors, respectively. These effectors are sequence-unrelated, except for allelic AVRa10 and AVRa22 that are co-maintained in pathogen populations in the form of a balanced polymorphism. Contrary to numerous examples of indirect recognition of bacterial effectors by plant NLRs, co-expression experiments with matching Mla-AVRa pairs indicate direct detection of the sequence-unrelated fungal effectors by MLA receptors.

A Golgi-Released Subpopulation of the Trans-Golgi Network Mediates Protein Secretion in Arabidopsis

Spatiotemporal coordination of protein trafficking among organelles is essential for eukaryotic cells. The post-Golgi interface, including the trans-Golgi network (TGN), is a pivotal hub for multiple trafficking pathways. The Golgi-released independent TGN (GI-TGN) is a compartment described only in plant cells, and its cellular and physiological roles remain elusive. In Arabidopsis (Arabidopsis thaliana), the SYNTAXIN OF PLANTS (SYP) 4 group Qa-SNARE (soluble N-ethylmaleimide) membrane fusion proteins are shared components of TGN and GI-TGN and regulate secretory and vacuolar transport. Here we reveal that GI-TGNs mediate the transport of the R-SNARE VESICLE-ASSOCIATED MEMBRANE PROTEIN (VAMP) 721 to the plasma membrane. In interactions with a nonadapted powdery mildew pathogen, the SYP4 group of SNAREs is required for the dynamic relocation of VAMP721 to plant-fungus contact sites via GI-TGNs, thereby facilitating complex formation with its cognate SNARE partner PENETRATION1 to restrict pathogen entry. Furthermore, quantitative proteomic analysis of leaf apoplastic fluid revealed constitutive and pathogen-inducible secretion of cell wall-modification enzymes in a SYP4- and VAMP721-dependent manner. Hence, the GI-TGN acts as a transit compartment between the Golgi apparatus and the plasma membrane. We propose a model in which the GA-TGN matures into the GI-TGN and then into secretory vesicles by increasing the abundance of VAMP721-dependent secretory pathway components.

Evolution of Disease Defense Genes and Their Regulators in Plants

Biotic stresses do damage to the growth and development of plants, and yield losses for some crops. Confronted with microbial infections, plants have evolved multiple defense mechanisms, which play important roles in the never-ending molecular arms race of plant–pathogen interactions. The complicated defense systems include pathogen-associated molecular patterns (PAMP) triggered immunity (PTI), effector triggered immunity (ETI), and the exosome-mediated cross-kingdom RNA interference (CKRI) system. Furthermore, plants have evolved a classical regulation system mediated by miRNAs to regulate these defense genes. Most of the genes/small RNAs or their regulators that involve in the defense pathways can have very rapid evolutionary rates in the longitudinal and horizontal co-evolution with pathogens. According to these internal defense mechanisms, some strategies such as molecular switch for the disease resistance genes, host-induced gene silencing (HIGS), and the new generation of RNA-based fungicides, have been developed to control multiple plant diseases. These broadly applicable new strategies by transgene or spraying ds/sRNA may lead to reduced application of pesticides and improved crop yield.

Subfamily-Specific Specialization of RGH1/MLA Immune Receptors in Wild Barley

The barley disease resistance (R) gene locus mildew locus A (Mla) provides isolate-specific resistance against the powdery mildew fungus Blumeria graminis hordei and has been introgressed into modern cultivars from diverse germplasms, including the wild relative Hordeum spontaneum. Known Mla disease resistance specificities to B. graminis hordei appear to encode allelic variants of the R gene homolog 1 (RGH1) family of nucleotide-binding domain and leucine-rich repeat (NLR) proteins. Here, we sequenced and assembled the transcriptomes of 50 H. spontaneum accessions representing nine populations distributed throughout the Fertile Crescent. The assembled Mla transcripts exhibited rich sequence diversity, linked neither to geographic origin nor population structure, and could be grouped into two similar-sized subfamilies based on two major N-terminal coiled-coil (CC) signaling domains that are both capable of eliciting cell death. The presence of positively selected sites located mainly in the C-terminal leucine-rich repeats of both MLA subfamilies, together with the fact that both CC signaling domains mediate cell death, implies that the two subfamilies are actively maintained in the population. Unexpectedly, known MLA receptor variants that confer B. graminis hordei resistance belong exclusively to one subfamily. Thus, signaling domain divergence, potentially as adaptation to distinct pathogen populations, is an evolutionary signature of functional diversification of an immune receptor. Copyright © 2018 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Identification of QTLs for powdery mildew (Podosphaera aphanis; syn. Sphaerotheca macularis f. sp. fragariae) susceptibility in cultivated strawberry (Fragaria ×ananassa)

Strawberry powdery mildew (Podosphaera aphanis Wallr.) is a pathogen which infects the leaves, fruit, stolon and flowers of the cultivated strawberry (Fragaria ×ananassa), causing major yield losses, primarily through unmarketable fruit. The primary commercial control of the disease is the application of fungicidal sprays. However, as the use of key active ingredients of commercial fungicides is becoming increasingly restricted, interest in developing novel strawberry cultivars exhibiting resistance to the pathogen is growing rapidly. In this study, a mapping population derived from a cross between two commercial strawberry cultivars ('Sonata' and 'Babette') was genotyped with single nucleotide polymorphism (SNP) markers from the Axiom iStraw90k genotyping array and phenotyped for powdery mildew susceptibility in both glasshouse and field environments. Three distinct, significant QTLs for powdery mildew resistance were identified across the two experiments. Through comparison with previous studies and scrutiny of the F. vesca genome sequence, candidate genes underlying the genetic control of this trait were identified.

Expression of a maize NBS gene ZmNBS42 enhances disease resistance in Arabidopsis

Expression of the ZmNBS42 in Arabidopsis plants conferred resistance to bacterial pathogens, providing potential resistance enhancement of maize in further genetic breeding. Nucleotide-binding site (NBS) domain proteins play critical roles in disease resistance. In this study, we isolate a novel NBS gene ZmNBS42 from maize and systematically investigate its function on disease resistance. We find that the expression levels of ZmNBS42 in maize leaf were strikingly increased in response to Bipolaris maydis inoculation and SA treatment. The spatial expression pattern analysis reveals that, during development, ZmNBS42 is ubiquitously highly expressed in maize root, leaf, stem, internode and seed, but lowly expressed in pericarp and embryo. To better understand the roles of ZmNBS42, we overexpressed ZmNBS42 in heterologous systems. Transient overexpression of ZmNBS42 in the leaves of Nicotiana benthamiana induces a hypersensitive response. ZmNBS42 overexpression (ZmNBS42-OE) Arabidopsis plants produced more SA content than Col-0 plants, and increased the expression levels of some defense-responsive genes compared to Col-0 plants. Moreover, the ZmNBS42-OE Arabidopsis plants displayed enhanced resistance against Pseudomonas syringae pathovar tomato DC3000 (Pst DC3000). These results together suggest that ZmNBS42 can serve as an important regulator in disease resistance, thus better understanding of ZmNBS42 would benefit the resistance enhancement in maize breeding programs.

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.

Genome-Wide Association Study Identifies NBS-LRR-Encoding Genes Related with Anthracnose and Common Bacterial Blight in the Common Bean

Nucleotide-binding site and leucine-rich repeat (NBS-LRR) genes represent the largest and most important disease resistance genes in plants. The genome sequence of the common bean (Phaseolus vulgaris L.) provides valuable data for determining the genomic organization of NBS-LRR genes. However, data on the NBS-LRR genes in the common bean are limited. In total, 178 NBS-LRR-type genes and 145 partial genes (with or without a NBS) located on 11 common bean chromosomes were identified from genome sequences database. Furthermore, 30 NBS-LRR genes were classified into Toll/interleukin-1 receptor (TIR)-NBS-LRR (TNL) types, and 148 NBS-LRR genes were classified into coiled-coil (CC)-NBS-LRR (CNL) types. Moreover, the phylogenetic tree supported the division of these PvNBS genes into two obvious groups, TNL types and CNL types. We also built expression profiles of NBS genes in response to anthracnose and common bacterial blight using qRT-PCR. Finally, we detected nine disease resistance loci for anthracnose (ANT) and seven for common bacterial blight (CBB) using the developed NBS-SSR markers. Among these loci, NSSR24, NSSR73, and NSSR265 may be located at new regions for ANT resistance, while NSSR65 and NSSR260 may be located at new regions for CBB resistance. Furthermore, we validated NSSR24, NSSR65, NSSR73, NSSR260, and NSSR265 using a new natural population. Our results provide useful information regarding the function of the NBS-LRR proteins and will accelerate the functional genomics and evolutionary studies of NBS-LRR genes in food legumes. NBS-SSR markers represent a wide-reaching resource for molecular breeding in the common bean and other food legumes. Collectively, our results should be of broad interest to bean scientists and breeders.

Homologues of the RPW8 Resistance Protein Are Localized to the Extrahaustorial Membrane that Is Likely Synthesized De Novo

Upon penetration of the host cell wall, the powdery mildew fungus develops a feeding structure named the haustorium in the invaded host cell. Concomitant with haustorial biogenesis, the extrahaustorial membrane (EHM) is formed to separate the haustorium from the host cell cytoplasm. The Arabidopsis resistance protein RPW8.2 is specifically targeted to the EHM where it activates haustorium-targeted resistance against powdery mildew. RPW8.2 belongs to a small family with six members in Arabidopsis (Arabidopsis thaliana). Whether Homologs of RPW8 (HR) 1 to HR4 are also localized to the EHM and contribute to resistance has not been determined. Here, we report that overexpression of HR1, HR2, or HR3 led to enhanced resistance to powdery mildew, while genetic depletion of HR2 or HR3 resulted in enhanced susceptibility, indicating that these RPW8 homologs contribute to basal resistance. Interestingly, we found that N-terminally YFP-tagged HR1 to HR3 are also EHM-localized. This suggests that EHM-targeting is an ancestral feature of the RPW8 family. Indeed, two RPW8 homologs from Brassica oleracea tested also exhibit EHM-localization. Domain swapping analysis between HR3 and RPW8.2 suggests that sequence diversification in the N-terminal 146 amino acids of RPW8.2 probably functionally distinguishes it from other family members. Moreover, we found that N-terminally YFP-tagged HR3 is also localized to the plasma membrane and the fungal penetration site (the papilla) in addition to the EHM. Using this unique feature of YFP-HR3, we obtained preliminary evidence to suggest that the EHM is unlikely derived from invagination of the plasma membrane, rather it may be mainly synthesized de novo.

Divergent evolution of multiple virus-resistance genes from a progenitor in Capsicum spp

Plants have evolved hundreds of nucleotide-binding and leucine-rich domain proteins (NLRs) as potential intracellular immune receptors, but the evolutionary mechanism leading to the ability to recognize specific pathogen effectors is elusive. Here, we cloned Pvr4 (a Potyvirus resistance gene in Capsicum annuum) and Tsw (a Tomato spotted wilt virus resistance gene in Capsicum chinense) via a genome-based approach using independent segregating populations. The genes both encode typical NLRs and are located at the same locus on pepper chromosome 10. Despite the fact that these two genes recognize completely different viral effectors, the genomic structures and coding sequences of the two genes are strikingly similar. Phylogenetic studies revealed that these two immune receptors diverged from a progenitor gene of a common ancestor. Our results suggest that sequence variations caused by gene duplication and neofunctionalization may underlie the evolution of the ability to specifically recognize different effectors. These findings thereby provide insight into the divergent evolution of plant immune receptors.

Molecular Mapping of PMR1, a Novel Locus Conferring Resistance to Powdery Mildew in Pepper (Capsicum annuum)

Powdery mildew, caused by Leveillula taurica, is a major fungal disease affecting greenhouse-grown pepper (Capsicum annuum). Powdery mildew resistance has a complex mode of inheritance. In the present study, we investigated a novel powdery mildew resistance locus, PMR1, using two mapping populations: 102 'VK515' F_2:3 families (derived from a cross between resistant parental line 'VK515R' and susceptible parental line 'VK515S') and 80 'PM Singang' F₂ plants (derived from the F₁ 'PM Singang' commercial hybrid). Genetic analysis of the F_2:3 'VK515' and F₂ 'PM Singang' populations revealed a single dominant locus for inheritance of the powdery mildew resistance trait. Genetic mapping showed that the PMR1 locus is located on syntenic regions of pepper chromosome 4 in a 4-Mb region between markers CZ2_11628 and HRM4.1.6 in 'VK515R'. Six molecular markers including one SCAR marker and five SNP markers were localized to a region 0 cM from the PMR1 locus. Two putative nucleotide-binding site leucine-rich repeat (NBS-LRR)-type disease resistance genes were identified in this PMR1 region. Genotyping-by-sequencing (GBS) and genetic mapping analysis revealed suppressed recombination in the PMR1 region, perhaps due to alien introgression. In addition, a comparison of species-specific InDel markers as well as GBS-derived SNP markers indicated that C. baccatum represents a possible source of such alien introgression of powdery mildew resistance into 'VK515R'. The molecular markers developed in this study will be especially helpful for marker-assisted selection in pepper breeding programs for powdery mildew resistance.

Marker development using SLAF-seq and whole-genome shotgun strategy to fine-map the semi-dwarf gene ari-e in barley

Background

Barley semi-dwarf genes have been extensively explored and widely used in barley breeding programs. The semi-dwarf gene ari-e from Golden Promise is an important gene associated with some agronomic traits and salt tolerance. While ari-e has been mapped on barley chromosome 5H using traditional markers and next-generation sequencing technologies, it has not yet been finely located on this chromosome.

Results

We integrated two methods to develop molecular markers for fine-mapping the semi-dwarf gene ari-e: (1) specific-length amplified fragment sequencing (SLAF-seq) with bulked segregant analysis (BSA) to develop SNP markers, and (2) the whole-genome shotgun sequence to develop InDels. Both SNP and InDel markers were developed in the target region and used for fine-mapping the ari-e gene. Linkage analysis showed that ari-e co-segregated with marker InDel-17 and was delimited by two markers (InDel-16 and DGSNP21) spanning 6.8 cM in the doubled haploid (DH) Dash × VB9104 population. The genetic position of ari-e was further confirmed in the Hindmarsh × W1 DH population which was located between InDel-7 and InDel-17. As a result, the overlapping region of the two mapping populations flanked by InDel-16 and InDel-17 was defined as the candidate region spanning 0.58 Mb on the POPSEQ physical map.

Conclusions

The current study demonstrated the SLAF-seq for SNP discovery and whole-genome shotgun sequencing for InDel development as an efficient approach to map complex genomic region for isolation of functional gene. The ari-e gene was fine mapped from 10 Mb to 0.58 Mb interval.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3247-4) contains supplementary material, which is available to authorized users.

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.

Modulation of Plant RAB GTPase-Mediated Membrane Trafficking Pathway at the Interface Between Plants and Obligate Biotrophic Pathogens

RAB5 is a small GTPase that acts in endosomal trafficking. In addition to canonical RAB5 members that are homologous to animal RAB5, land plants harbor a plant-specific RAB5, the ARA6 group, which regulates trafficking events distinct from canonical RAB5 GTPases. Here, we report that plant RAB5, both canonical and plant-specific members, accumulate at the interface between host plants and biotrophic fungal and oomycete pathogens. Biotrophic fungi and oomycetes colonize living plant tissues by establishing specialized infection hyphae, the haustorium, within host plant cells. We found that Arabidopsis thaliana ARA6/RABF1, a plant-specific RAB5, is localized to the specialized membrane that surrounds the haustorium, the extrahaustorial membrane (EHM), formed by the A. thaliana-adapted powdery mildew fungus Golovinomyces orontii Whereas the conventional RAB5 ARA7/RABF2b was also localized to the EHM, endosomal SNARE (soluble N-ethylmaleimide-sensitive factor attachment protein receptor) and RAB5-activating proteins were not, which suggests that the EHM has modified endosomal characteristic. The recruitment of host RAB5 to the EHM was a property shared by the barley-adapted powdery mildew fungus Blumeria graminis f.sp. hordei and the oomycete Hyaloperonospora arabidopsidis, but the extrahyphal membrane surrounding the hypha of the hemibiotrophic fungus Colletotrichum higginsianum at the biotrophic stage was devoid of RAB5. The localization of RAB5 to the EHM appears to correlate with the functionality of the haustorium. Our discovery sheds light on a novel relationship between plant RAB5 and obligate biotrophic pathogens.

A Cluster of Nucleotide-Binding Site-Leucine-Rich Repeat Genes Resides in a Barley Powdery Mildew Resistance Quantitative Trait Loci on 7HL

Powdery mildew causes severe yield losses in barley production worldwide. Although many resistance genes have been described, only a few have already been cloned. A strong QTL (quantitative trait locus) conferring resistance to a wide array of powdery mildew isolates was identified in a Spanish barley landrace on the long arm of chromosome 7H. Previous studies narrowed down the QTL position, but were unable to identify candidate genes or physically locate the resistance. In this study, the exome of three recombinant lines from a high-resolution mapping population was sequenced and analyzed, narrowing the position of the resistance down to a single physical contig. Closer inspection of the region revealed a cluster of closely related NBS-LRR (nucleotide-binding site-leucine-rich repeat containing protein) genes. Large differences were found between the resistant lines and the reference genome of cultivar Morex, in the form of PAV (presence-absence variation) in the composition of the NBS-LRR cluster. Finally, a template-guided assembly was performed and subsequent expression analysis revealed that one of the new assembled candidate genes is transcribed. In summary, the results suggest that NBS-LRR genes, absent from the reference and the susceptible genotypes, could be functional and responsible for the powdery mildew resistance. The procedure followed is an example of the use of NGS (next-generation sequencing) tools to tackle the challenges of gene cloning when the target gene is absent from the reference genome.

Large-scale bioinformatic analysis of the regulation of the disease resistance NBS gene family by microRNAs in Poaceae

In the present study, we have screened 71, 713, 525, 119 and 241 mature miRNA variants from Hordeum vulgare, Oryza sativa, Brachypodium distachyon, Triticum aestivum, and Sorghum bicolor, respectively, and classified them with respect to their conservation status and expression levels. These Poaceae non-redundant miRNA species (1,669) were distributed over a total of 625 MIR families, among which only 54 were conserved across two or more plant species, confirming the relatively recent evolutionary differentiation of miRNAs in grasses. On the other hand, we have used 257 H. vulgare, 286T. aestivum, 119 B. distachyon, 269 O. sativa, and 139 S. bicolor NBS domains, which were either mined directly from the annotated proteomes, or predicted from whole genome sequence assemblies. The hybridization potential between miRNAs and their putative NBS genes targets was analyzed, revealing that at least 454 NBS genes from all five Poaceae were potentially regulated by 265 distinct miRNA species, most of them expressed in leaves and predominantly co-expressed in additional tissues. Based on gene ontology, we could assign these probable miRNA target genes to 16 functional groups, among which three conferring resistance to bacteria (Rpm1, Xa1 and Rps2), and 13 groups of resistance to fungi (Rpp8,13, Rp3, Tsn1, Lr10, Rps1-k-1, Pm3, Rpg5, and MLA1,6,10,12,13). The results of the present analysis provide a large-scale platform for a better understanding of biological control strategies of disease resistance genes in Poaceae, and will serve as an important starting point for enhancing crop disease resistance improvement by means of transgenic lines with artificial miRNAs.

Coexpression network analysis of the genes regulated by two types of resistance responses to powdery mildew in wheat

Powdery mildew disease caused by Blumeria graminis f. sp. tritici (Bgt) inflicts severe economic losses in wheat crops. A systematic understanding of the molecular mechanisms involved in wheat resistance to Bgt is essential for effectively controlling the disease. Here, using the diploid wheat Triticum urartu as a host, the genes regulated by immune (IM) and hypersensitive reaction (HR) resistance responses to Bgt were investigated through transcriptome sequencing. Four gene coexpression networks (GCNs) were developed using transcriptomic data generated for 20 T. urartu accessions showing IM, HR or susceptible responses. The powdery mildew resistance regulated (PMRR) genes whose expression was significantly correlated with Bgt resistance were identified, and they tended to be hubs and enriched in six major modules. A wide occurrence of negative regulation of PMRR genes was observed. Three new candidate immune receptor genes (TRIUR3_13045, TRIUR3_01037 and TRIUR3_06195) positively associated with Bgt resistance were discovered. Finally, the involvement of TRIUR3_01037 in Bgt resistance was tentatively verified through cosegregation analysis in a F2 population and functional expression assay in Bgt susceptible leaf cells. This research provides insights into the global network properties of PMRR genes. Potential molecular differences between IM and HR resistance responses to Bgt are discussed.

Disease Resistance Gene Analogs (RGAs) in Plants

Plants have developed effective mechanisms to recognize and respond to infections caused by pathogens. Plant resistance gene analogs (RGAs), as resistance (R) gene candidates, have conserved domains and motifs that play specific roles in pathogens' resistance. Well-known RGAs are nucleotide binding site leucine rich repeats, receptor like kinases, and receptor like proteins. Others include pentatricopeptide repeats and apoplastic peroxidases. RGAs can be detected using bioinformatics tools based on their conserved structural features. Thousands of RGAs have been identified from sequenced plant genomes. High-density genome-wide RGA genetic maps are useful for designing diagnostic markers and identifying quantitative trait loci (QTL) or markers associated with plant disease resistance. This review focuses on recent advances in structures and mechanisms of RGAs, and their identification from sequenced genomes using bioinformatics tools. Applications in enhancing fine mapping and cloning of plant disease resistance genes are also discussed.

Dominant negative RPW8.2 fusion proteins reveal the importance of haustorium-oriented protein trafficking for resistance against powdery mildew in Arabidopsis

Powdery mildew fungi form feeding structures called haustoria inside epidermal cells of host plants to extract photosynthates for their epiphytic growth and reproduction. The haustorium is encased by an interfacial membrane termed the extrahaustorial membrane (EHM). The atypical resistance protein RPW8.2 from Arabidopsis is specifically targeted to the EHM where RPW8.2 activates haustorium-targeted (thus broad-spectrum) resistance against powdery mildew fungi. EHM-specific localization of RPW8.2 suggests the existence of an EHM-oriented protein/membrane trafficking pathway during EHM biogenesis. However, the importance of this specific trafficking pathway for host defense has not been evaluated via a genetic approach without affecting other trafficking pathways. Here, we report that expression of EHM-oriented, nonfunctional RPW8.2 chimeric proteins exerts dominant negative effect over functional RPW8.2 and potentially over other EHM-localized defense proteins, thereby compromising both RPW8.2-mediated and basal resistance to powdery mildew. Thus, our results highlight the importance of the EHM-oriented protein/membrane trafficking pathway for host resistance against haustorium-forming pathogens such as powdery mildew fungi.

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.

Dynamic evolution of NBS-LRR genes in bread wheat and its progenitors

Extensive studies have focused on the largest class of disease resistance genes (nucleotide binding site-leucine-rich repeat, NBS-LRR) in various plants. However, no research on the dynamic evolution of these genes in domesticated species and their progenitors has been reported. Recently published genome sequences of bread wheat and its two ancestors provide a good opportunity for comparing NBS-encoding genes between ancestors and their progeny. Over 2000 NBS-encoding genes have been identified in bread wheat, which is the largest number having been reported so far. Compared with other grass species, its two progenitors also contained more NBS-encoding genes, indicating that there was an expansion of these genes in their common ancestor. Interestingly, the inherited relationships of NBS-LRR genes among the bread wheat and its two progenitors were ambiguous and only 3 % single-copy orthologues retained gene order in three-way genome comparisons of the three genomes. Lots of NBS-encoding genes present in the either ancestor could not be found in the bread wheat. These results indicated that NBS-LRR genes in bread wheat might have evolved rapidly through a rapid loss of ancestor genes.

Dynamic evolution of NBS-LRR genes in bread wheat and its progenitors

Extensive studies have focused on the largest class of disease resistance genes (nucleotide binding site-leucine-rich repeat, NBS-LRR) in various plants. However, no research on the dynamic evolution of these genes in domesticated species and their progenitors has been reported. Recently published genome sequences of bread wheat and its two ancestors provide a good opportunity for comparing NBS-encoding genes between ancestors and their progeny. Over 2000 NBS-encoding genes have been identified in bread wheat, which is the largest number having been reported so far. Compared with other grass species, its two progenitors also contained more NBS-encoding genes, indicating that there was an expansion of these genes in their common ancestor. Interestingly, the inherited relationships of NBS-LRR genes among the bread wheat and its two progenitors were ambiguous and only 3 % single-copy orthologues retained gene order in three-way genome comparisons of the three genomes. Lots of NBS-encoding genes present in the either ancestor could not be found in the bread wheat. These results indicated that NBS-LRR genes in bread wheat might have evolved rapidly through a rapid loss of ancestor genes.

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.

Genetic dissection of a TIR-NB-LRR locus from the wild North American grapevine species Muscadinia rotundifolia identifies paralogous genes conferring resistance to major fungal and oomycete pathogens in cultivated grapevine

The most economically important diseases of grapevine cultivation worldwide are caused by the fungal pathogen powdery mildew (Erysiphe necator syn. Uncinula necator) and the oomycete pathogen downy mildew (Plasmopara viticola). Currently, grapegrowers rely heavily on the use of agrochemicals to minimize the potentially devastating impact of these pathogens on grape yield and quality. The wild North American grapevine species Muscadinia rotundifolia was recognized as early as 1889 to be resistant to both powdery and downy mildew. We have now mapped resistance to these two mildew pathogens in M. rotundifolia to a single locus on chromosome 12 that contains a family of seven TIR-NB-LRR genes. We further demonstrate that two highly homologous (86% amino acid identity) members of this gene family confer strong resistance to these unrelated pathogens following genetic transformation into susceptible Vitis vinifera winegrape cultivars. These two genes, designated resistance to Uncinula necator (MrRUN1) and resistance to Plasmopara viticola (MrRPV1) are the first resistance genes to be cloned from a grapevine species. Both MrRUN1 and MrRPV1 were found to confer resistance to multiple powdery and downy mildew isolates from France, North America and Australia; however, a single powdery mildew isolate collected from the south-eastern region of North America, to which M. rotundifolia is native, was capable of breaking MrRUN1-mediated resistance. Comparisons of gene organization and coding sequences between M. rotundifolia and the cultivated grapevine V. vinifera at the MrRUN1/MrRPV1 locus revealed a high level of synteny, suggesting that the TIR-NB-LRR genes at this locus share a common ancestor.

HvWRKY10, HvWRKY19, and HvWRKY28 regulate Mla-triggered immunity and basal defense to barley powdery mildew

WRKY proteins represent a large family of transcription factors (TF), involved in plant development and defense. In all, 60 unique barley TF have been annotated that contain the WRKY domain; 26 of these are represented on the Barley1 GeneChip. Time-course expression profiles of these 26 HvWRKY TF were analyzed to investigate their role in mildew locus a (Mla)-mediated immunity to Blumeria graminis f. sp. hordei, causal agent of powdery mildew disease. Inoculation-responsive, Mla-specified interactions with B. graminis f. sp. hordei revealed that 12 HvWRKY were differentially expressed: 10 highly upregulated and two significantly downregulated. Barley stripe mosaic virus-induced gene silencing of HvWRKY10, HvWRKY19, and HvWRKY28 compromised resistance-gene-mediated defense to powdery mildew in genotypes harboring both Rar1-dependent and Rar1-independent Mla alleles, indicating that these WRKY TF play key roles in effector-triggered immunity. Comprehensive yeast two-hybrid analyses, however, did not reveal a direct interaction between these three nuclear-localized WRKY TF and MLA. Transient overexpression of all three WRKY TF in single cells expressing Mlo, which encodes a negative regulator of penetration resistance, significantly decreased susceptibility. Taken together, these loss- and gain-of-function studies demonstrate that HvWRKY10, HvWRKY19, and HvWRKY28 positively regulate the barley transcriptome in response to invasion by B. graminis f. sp. hordei.

PAPP2C interacts with the atypical disease resistance protein RPW8.2 and negatively regulates salicylic acid-dependent defense responses in Arabidopsis

Many fungal and oomycete pathogens differentiate a feeding structure named the haustorium to extract nutrition from the plant epidermal cell. The atypical resistance (R) protein RPW8.2 activates salicylic acid (SA)-dependent, haustorium-targeted defenses against Golovinomyces spp., the causal agents of powdery mildew diseases on multiple plant species. How RPW8.2 activates defense remains uncharacterized. Here, we report that RPW8.2 interacts with the phytochrome-associated protein phosphatase type 2C (PAPP2C) in yeast and in planta as evidenced by co-immunoprecipitation and bimolecular fluorescence complementation assays. Down-regulation of PAPP2C by RNA interference (RNAi) in Col-0 plants lacking RPW8.2 leads to leaf spontaneous cell death and enhanced disease resistance to powdery mildew via the SA-dependent signaling pathway. Moreover, down-regulation of PAPP2C by RNAi in the RPW8.2 background results in strong HR-like cell death, which correlates with elevated RPW8.2 expression. We further demonstrate that hemagglutinin (HA)-tagged PAPP2C prepared from tobacco leaf cells transiently transformed with HA-PAPP2C possesses phosphatase activity. In addition, silencing a rice gene (Os04g0452000) homologous to PAPP2C also results in spontaneous cell death in rice. Combined, our results suggest that RPW8.2 is functionally connected with PAPP2C and that PAPP2C negatively regulates SA-dependent basal defense against powdery mildew in Arabidopsis.

Comparative analysis of zinc finger proteins involved in plant disease resistance

A meta-analysis was performed to understand the role of zinc finger domains in proteins of resistance (R) genes cloned from different crops. We analyzed protein sequences of seventy R genes of various crops in which twenty six proteins were found to have zinc finger domains along with nucleotide binding sites - leucine rice repeats (NBS-LRR) domains. We identified thirty four zinc finger domains in the R proteins of nine crops and were grouped into 19 types of zinc fingers. The size of individual zinc finger domain within the R genes varied from 11 to 84 amino acids, whereas the size of proteins containing these domains varied from 263 to 1305 amino acids. The biophysical analysis revealed that molecular weight of Pi54 zinc finger was lowest whereas the highest one was found in rice Pib zinc finger named as Transposes Transcription Factor (TTF). The instability (R(2) =0.95) and the aliphatic (R(2) =0.94) indices profile of zinc finger domains follows the polynomial distribution pattern. The pairwise identity analysis showed that the Lin11, Isl-1 & Mec-3 (LIM) zinc finger domain of rice blast resistance protein pi21 have 12.3% similarity with the nuclear transcription factor, X-box binding-like 1 (NFX) type zinc finger domain of Pi54 protein. For the first time, we reported that Pi54 (Pi-k(h)-Tetep), a rice blast resistance (R) protein have a small zinc finger domain of NFX type located on the C-terminal in between NBS and LRR domains of the R-protein. Compositional analysis depicted by the helical wheel diagram revealed the presence of a hydrophobic region within this domain which might help in exposing the LRR region for a possible R-Avr interaction. This domain is unique among all other cloned plant disease resistance genes and might play an important role in broad-spectrum nature of rice blast resistance gene Pi54.

High-density genetic and physical bin mapping of wheat chromosome 1D reveals that the powdery mildew resistance gene Pm24 is located in a highly recombinogenic region

Genetic maps of wheat chromosome 1D consisting of 57 microsatellite marker loci were constructed using Chinese Spring (CS) × Chiyacao F(2) and the International Triticeae Mapping Initiative (ITMI) recombinant inbred lines (RILs) mapping populations. Marker order was consistent, but genetic distances of neighboring markers were different in two populations. Physical bin map of 57 microsatellite marker loci was generated by means of 10 CS 1D deletion lines. The physical bin mapping indicated that microsatellite marker loci were not randomly distributed on chromosome 1D. Nineteen of the 24 (79.2%) microsatellite markers were mapped in the distal 30% genomic region of 1DS, whereas 25 of the 33 (75.8%) markers were assigned to the distal 59% region of 1DL. The powdery mildew resistance gene Pm24, originating from the Chinese wheat landrace Chiyacao, was previously mapped in the vicinity of the centromere on the short arm of chromosome 1D. A high density genetic map of chromosome 1D was constructed, consisting of 36 markers and Pm24, with a total map length of 292.7 cM. Twelve marker loci were found to be closely linked to Pm24. Pm24 was flanked by Xgwm789 (Xgwm603) and Xbarc229 with genetic distances of 2.4 and 3.6 cM, respectively, whereas a microsatellite marker Xgwm1291 co-segregated with Pm24. The microsatellite marker Xgwm1291 was assigned to the bin 1DS5-0.70-1.00 of the chromosome arm 1DS. It could be concluded that Pm24 is located in the '1S0.8 gene-rich region', a highly recombinogenic region of wheat. The results presented here would provide a start point for the map-based cloning of Pm24.

A genome-wide genetic map of NB-LRR disease resistance loci in potato

Like all plants, potato has evolved a surveillance system consisting of a large array of genes encoding for immune receptors that confer resistance to pathogens and pests. The majority of these so-called resistance or R proteins belong to the super-family that harbour a nucleotide binding and a leucine-rich-repeat domain (NB-LRR). Here, sequence information of the conserved NB domain was used to investigate the genome-wide genetic distribution of the NB-LRR resistance gene loci in potato. We analysed the sequences of 288 unique BAC clones selected using filter hybridisation screening of a BAC library of the diploid potato clone RH89-039-16 (S. tuberosum ssp. tuberosum) and a physical map of this BAC library. This resulted in the identification of 738 partial and full-length NB-LRR sequences. Based on homology of these sequences with known resistance genes, 280 and 448 sequences were classified as TIR-NB-LRR (TNL) and CC-NB-LRR (CNL) sequences, respectively. Genetic mapping revealed the presence of 15 TNL and 32 CNL loci. Thirty-six are novel, while three TNL loci and eight CNL loci are syntenic with previously identified functional resistance genes. The genetic map was complemented with 68 universal CAPS markers and 82 disease resistance trait loci described in literature, providing an excellent template for genetic studies and applied research in potato.

Over-expression in the nucleotide-binding site-leucine rich repeat gene DEPG1 increases susceptibility to bacterial leaf streak disease in transgenic rice plants

Bacterial leaf streak of rice (BLS) caused by Xanthomonas oryzae pv. oryzicola (Xoc) is a widely-spread disease in the main rice-producing areas of the world. Investigating the genes that play roles in rice-Xoc interactions helps us to understand the defense signaling pathway in rice. Here we report a differentially expressed protein gene (DEPG1), which regulates susceptibility to BLS. DEPG1 is a nucleotide-binding site (NBS)-leucine rich repeat (LRR) gene, and the deduced protein sequence of DEPG1 has approximately 64% identity with that of the disease resistance gene Pi37. Phylogenetic analysis of DEPG1 and the 18 characterized NBS-LRR genes revealed that DEPG1 is more closely related to Pi37. DEPG1 protein is located to the cytoplasm, which was confirmed by transient expression of DEPG1-GFP (green fluorescent protein) fusion construct in onion epidermal cells. Semi-quantitative PCR assays showed that DEPG1 is widely expressed in rice, and is preferentially expressed in internodes, leaf blades, leaf sheaths and flag leaves. Observation of cross sections of leaves from the transgenic plants with a DEPG1-promoter::glucuronidase (GUS) fusion gene revealed that DEPG1 is also highly expressed in mesophyll tissues where Xoc mainly colonizes. Additionally, Xoc negatively regulates expression of DEPG1 at the early stage of the pathogen infection, and so do the three defense-signal compounds including salicylic acid (SA), methyl jasmonate (MeJA) and 1-aminocyclopropane-1-carboxylic-acid (ACC). Transgenic rice plants overexpressing DEPG1 exhibit enhanced susceptibility to Xoc compared to the wild-type controls. Moreover, enhanced susceptibility to Xoc may be mediated by inhibition of the expression of some SA biosynthesis-related genes and pathogenesis-related genes that may contribute to the disease resistance. Taken together, DEPG1 plays roles in the interactions between rice and BLS pathogen Xoc.