ZRK atypical kinases: emerging signaling components of plant immunity

EDR1 associates with its homologs to synergistically regulate plant immunity in Arabidopsis

ENHANCED DISEASE RESISTANCE 1 (EDR1), a Raf-like mitogen-activated protein kinase (MAPK) kinase kinase (MAPKKK), is a negative regulator of resistance. There are three homologs, RAF3/4/5, of EDR1 in Arabidopsis. However, the roles of RAF3/4/5 in resistance and their functional link with EDR1 in plant immunity remain unclear. Here, we showed that the raf3/4/5 triple mutant displayed wild-type-like phenotypes to the powdery mildew pathogen Golovinomyces cichoracearum UCSC1 and the bacterial pathogen Pseudomonas syringae pv. tomato (Pto) DC3000. However, the edr1 raf3/4/5 quadruple mutant exhibited enhanced resistance to G. cichoracearum UCSC1 and Pto DC3000 compared to edr1. Consistently, MPK3/6 kinase activity was more highly activated in edr1 raf3/4/5 than that in edr1. Moreover, the enhanced resistance of edr1 raf3/4/5 required SALICYLIC ACID INDUCTION DEFICIENT 2 (SID2), an isochorismate synthase required for salicylic acid (SA) synthesis. Additionally, unlike EDR1, RAF3/4/5 weakly and indirectly associated with MKK4/5, and EDR1 was directly associated with RAF3/4/5. Taken together, these data indicate that EDR1 associates with RAF3/4/5, and they may function together to synergistically suppress MAPK cascades activation, which reveal the complexity and importance of Raf-like MAPKKKs in plant immunity regulation.

ZAR1: Guardian of plant kinases

A key facet of innate immunity in plants entails the recognition of pathogen "effector" virulence proteins by host Nucleotide-Binding Leucine-Rich Repeat Receptors (NLRs). Among characterized NLRs, the broadly conserved ZAR1 NLR is particularly remarkable due to its capacity to recognize at least six distinct families of effectors from at least two bacterial genera. This expanded recognition spectrum is conferred through interactions between ZAR1 and a dynamic network of two families of Receptor-Like Cytoplasmic Kinases (RLCKs): ZED1-Related Kinases (ZRKs) and PBS1-Like Kinases (PBLs). In this review, we survey the history of functional studies on ZAR1, with an emphasis on how the ZAR1-RLCK network functions to trap diverse effectors. We discuss 1) the dynamics of the ZAR1-associated RLCK network; 2) the specificity between ZRKs and PBLs; and 3) the specificity between effectors and the RLCK network. We posit that the shared protein fold of kinases and the switch-like properties of their interactions make them ideal effector sensors, enabling ZAR1 to act as a broad spectrum guardian of host kinases.

Study of natural diversity in response to a key pathogenicity regulator of Ralstonia solanacearum reveals new susceptibility genes in Arabidopsis thaliana

Ralstonia solanacearum gram-negative phytopathogenic bacterium exerts its virulence through a type III secretion system (T3SS) that translocates type III effectors (T3Es) directly into the host cells. T3E secretion is finely controlled at the posttranslational level by helper proteins, T3SS control proteins, and type III chaperones. The HpaP protein, one of the type III secretion substrate specificity switch (T3S4) proteins, was previously highlighted as a virulence factor on Arabidopsis thaliana Col-0 accession. In this study, we set up a genome-wide association analysis to explore the natural diversity of response to the hpaP mutant of two A. thaliana mapping populations: a worldwide collection and a local population. Quantitative genetic variation revealed different genetic architectures in both mapping populations, with a global delayed response to the hpaP mutant compared to the GMI1000 wild-type strain. We have identified several quantitative trait loci (QTLs) associated with the hpaP mutant inoculation. The genes underlying these QTLs are involved in different and specific biological processes, some of which were demonstrated important for R. solanacearum virulence. We focused our study on four candidate genes, RKL1, IRE3, RACK1B, and PEX3, identified using the worldwide collection, and validated three of them as susceptibility factors. Our findings demonstrate that the study of the natural diversity of plant response to a R. solanacearum mutant in a key regulator of virulence is an original and powerful strategy to identify genes directly or indirectly targeted by the pathogen.

Advances in Multi-Omics Approaches for Molecular Breeding of Black Rot Resistance in Brassica oleracea L

Brassica oleracea is one of the most important species of the Brassicaceae family encompassing several economically important vegetables produced and consumed worldwide. But its sustainability is challenged by a range of pathogens, among which black rot, caused by Xanthomonas campestris pv. campestris (Xcc), is the most serious and destructive seed borne bacterial disease, causing huge yield losses. Host-plant resistance could act as the most effective and efficient solution to curb black rot disease for sustainable production of B. oleracea. Recently, 'omics' technologies have emerged as promising tools to understand the host-pathogen interactions, thereby gaining a deeper insight into the resistance mechanisms. In this review, we have summarized the recent achievements made in the emerging omics technologies to tackle the black rot challenge in B. oleracea. With an integrated approach of the omics technologies such as genomics, proteomics, transcriptomics, and metabolomics, it would allow better understanding of the complex molecular mechanisms underlying black rot resistance. Due to the availability of sequencing data, genomics and transcriptomics have progressed as expected for black rot resistance, however, other omics approaches like proteomics and metabolomics are lagging behind, necessitating a holistic and targeted approach to address the complex questions of Xcc-Brassica interactions. Genomic studies revealed that the black rot resistance is a complex trait and is mostly controlled by quantitative trait locus (QTL) with minor effects. Transcriptomic analysis divulged the genes related to photosynthesis, glucosinolate biosynthesis and catabolism, phenylpropanoid biosynthesis pathway, ROS scavenging, calcium signalling, hormonal synthesis and signalling pathway are being differentially expressed upon Xcc infection. Comparative proteomic analysis in relation to susceptible and/or resistance interactions with Xcc identified the involvement of proteins related to photosynthesis, protein biosynthesis, processing and degradation, energy metabolism, innate immunity, redox homeostasis, and defence response and signalling pathways in Xcc-Brassica interaction. Specifically, most of the studies focused on the regulation of the photosynthesis-related proteins as a resistance response in both early and later stages of infection. Metabolomic studies suggested that glucosinolates (GSLs), especially aliphatic and indolic GSLs, its subsequent hydrolysis products, and defensive metabolites synthesized by jasmonic acid (JA)-mediated phenylpropanoid biosynthesis pathway are involved in disease resistance mechanisms against Xcc in Brassica species. Multi-omics analysis showed that JA signalling pathway is regulating resistance against hemibiotrophic pathogen like Xcc. So, the bonhomie between omics technologies and plant breeding is going to trigger major breakthroughs in the field of crop improvement by developing superior cultivars with broad-spectrum resistance. If multi-omics tools are implemented at the right scale, we may be able to achieve the maximum benefits from the minimum. In this review, we have also discussed the challenges, future prospects, and the way forward in the application of omics technologies to accelerate the breeding of B. oleracea for disease resistance. A deeper insight about the current knowledge on omics can offer promising results in the breeding of high-quality disease-resistant crops.

Evolutionarily distinct resistance proteins detect a pathogen effector through its association with different host targets

Knowledge of the evolutionary processes which govern pathogen recognition is critical to understanding durable disease resistance. We determined how Phytophthora infestans effector PiAVR2 is recognised by evolutionarily distinct resistance proteins R2 and Rpi-mcq1. We employed yeast two-hybrid, co-immunoprecipitation, virus-induced gene silencing, transient overexpression, and phosphatase activity assays to investigate the contributions of BSL phosphatases to R2- and Rpi-mcq1-mediated hypersensitive response (R2 HR and Rpi-mcq1 HR, respectively). Silencing PiAVR2 target BSL1 compromises R2 HR. Rpi-mcq1 HR is compromised only when BSL2 and BSL3 are silenced. BSL1 overexpression increases R2 HR and compromises Rpi-mcq1. However, overexpression of BSL2 or BSL3 enhances Rpi-mcq1 and compromises R2 HR. Okadaic acid, which inhibits BSL phosphatase activity, suppresses both recognition events. Moreover, expression of a BSL1 phosphatase-dead (PD) mutant suppresses R2 HR, whereas BSL2-PD and BSL3-PD mutants suppress Rpi-mcq1 HR. R2 interacts with BSL1 in the presence of PiAVR2, but not with BSL2 and BSL3, whereas no interactions were detected between Rpi-mcq1 and BSLs. Thus, BSL1 activity and association with R2 determine recognition of PiAVR2 by R2, whereas BSL2 and BSL3 mediate Rpi-mcq1 perception of PiAVR2. R2 and Rpi-mcq1 utilise distinct mechanisms to detect PiAVR2 based on association with different BSLs, highlighting central roles of these effector targets for both disease and disease resistance.

A tomato receptor-like cytoplasmic kinase, SlZRK1, acts as a negative regulator in wound-induced jasmonic acid accumulation and insect resistance

Jasmonates accumulate rapidly and act as key regulators in response to mechanical wounding, but few studies have linked receptor-like cytoplasmic kinases (RLCKs) to wound-induced jasmonic acid (JA) signaling cascades. Here, we identified a novel wounding-induced RLCK-XII-2 subfamily member (SlZRK1) in tomato (Solanum lycopersicum) that was closely related to Arabidopsis HOPZ-ETI-DEFICIENT 1 (ZED1)-related kinases 1 based on phylogenetic analysis. SlZRK1 was targeted to the plasma membrane of tobacco mesophyll protoplasts as determined by transient co-expression with the plasma membrane marker mCherry-H+-ATPase. Catalytic residue sequence analysis and an in vitro kinase assay indicated that SlZRK1 may act as a pseudokinase. To further analyse the function of SlZRK1, we developed two stable knock-out mutants by CRISPR/Cas9. Loss of SlZRK1 significantly altered the expression of genes involved in JA biosynthesis, salicylic acid biosynthesis, and ethylene response. Furthermore, after mechanical wounding treatment, slzrk1 mutants increased transcription of early wound-inducible genes involved in JA biosynthesis and signaling. In addition, JA accumulation after wounding and plant resistance to herbivorous insects also were enhanced. Our findings expand plant regulatory networks in the wound-induced JA production by adding RLCKs as a new component in the wound signal transduction pathway.

Cell Death in Plant Immunity

Pathogen recognition by the plant immune system leads to defense responses that are often accompanied by a form of regulated cell death known as the hypersensitive response (HR). HR shares some features with regulated necrosis observed in animals. Genetically, HR can be uncoupled from local defense responses at the site of infection and its role in immunity may be to activate systemic responses in distal parts of the organism. Recent advances in the field reveal conserved cell death-specific signaling modules that are assembled by immune receptors in response to pathogen-derived effectors. The structural elucidation of the plant resistosome-an inflammasome-like structure that may attach to the plasma membrane on activation-opens the possibility that HR cell death is mediated by the formation of pores at the plasma membrane. Necrotrophic pathogens that feed on dead tissue have evolved strategies to trigger the HR cell death pathway as a survival strategy. Ectopic activation of immunomodulators during autoimmune reactions can also promote HR cell death. In this perspective, we discuss the role and regulation of HR in these different contexts.

Novel organ-specific genetic factors for quantitative resistance to late blight in potato

Potato, Solanum tuberosum, is one of the major consumed food in the world, being the basis of the diet of millions of people. The main limiting and destructive disease of potato is late blight, caused by Phytophtora infestans. Here, we present a multi-environmental analysis of the response to P. infestans using an association panel of 150 accessions of S. tuberosum Group Phureja, evaluated in two localities in Colombia. Disease resistance data were merged with a genotyping matrix of 83,862 SNPs obtained by 2b-restriction site–associated DNA and Genotyping by sequencing approaches into a Genome-wide association study. We are reporting 16 organ-specific QTL conferring resistance to late blight. These QTL explain between 13.7% and 50.9% of the phenotypic variance. Six and ten QTL were detected for resistance response in leaves and stem, respectively. In silico analysis revealed 15 candidate genes for resistance to late blight. Four of them have no functional genome annotation, while eleven candidate genes code for diverse proteins, including a leucine-rich repeat kinase.

Perturbations of the ZED1 pseudokinase activate plant immunity

The Pseudomonas syringae acetyltransferase HopZ1a is delivered into host cells by the type III secretion system to promote bacterial growth. However, in the model plant host Arabidopsis thaliana, HopZ1a activity results in an effector-triggered immune response (ETI) that limits bacterial proliferation. HopZ1a-triggered immunity requires the nucleotide-binding, leucine-rich repeat domain (NLR) protein, ZAR1, and the pseudokinase, ZED1. Here we demonstrate that HopZ1a can acetylate members of a family of ‘receptor-like cytoplasmic kinases’ (RLCK family VII; also known as PBS1-like kinases, or PBLs) and promote their interaction with ZED1 and ZAR1 to form a ZAR1-ZED1-PBL ternary complex. Interactions between ZED1 and PBL kinases are determined by the pseudokinase features of ZED1, and mutants designed to restore ZED1 kinase motifs can (1) bind to PBLs, (2) recruit ZAR1, and (3) trigger ZAR1-dependent immunity in planta, all independently of HopZ1a. A ZED1 mutant that mimics acetylation by HopZ1a also triggers immunity in planta, providing evidence that effector-induced perturbations of ZED1 also activate ZAR1. Overall, our results suggest that interactions between these two RLCK families are promoted by perturbations of structural features that distinguish active from inactive kinase domain conformations. We propose that effector-induced interactions between ZED1/ZRK pseudokinases (RLCK family XII) and PBL kinases (RLCK family VII) provide a sensitive mechanism for detecting perturbations of either kinase family to activate ZAR1-mediated ETI.

All plants must ward off potentially infectious microbes, and those grown in large-scale crop operations are especially vulnerable to the rapid spread of disease by successful pathogens. Although many bacteria and fungi can supress plant immune responses by producing specialized virulence proteins called ‘effectors’, these effectors can also trigger immune responses that render plants resistant to infection. We studied the molecular mechanisms underlying one such effector-triggered immune response elicited by the bacterial effector HopZ1a in the model plant host Arabidopsis thaliana. We have shown that HopZ1a promotes binding between a ZED1, a ‘pseudokinase’ required for HopZ1a-triggered immunity, and several ‘true kinases’ (known as PBLs) that are likely targets of HopZ1a activity in planta. HopZ1a-induced ZED1-PBL interactions also recruit ZAR1, an Arabidopsis ‘resistance protein’ previously implicated in HopZ1a-triggered immunity. Importantly, ZED1 mutants that restore degenerate kinase motifs can bridge interactions between PBLs and ZAR1 (independently of HopZ1a) and trigger immunity in planta. Our results suggest that equilibria between active and inactive kinase domain conformations regulate ZED1-PBL interactions and formation of ternary complexes with ZAR1. Improved models describing molecular interactions between immunity determinants, effectors and effector targets will inform efforts to exploit natural diversity for development of crops with enhanced disease resistance.

Novel organ-specific genetic factors for quantitative resistance to late blight in potato

Potato, Solanum tuberosum, is one of the major consumed food in the world, being the basis of the diet of millions of people. The main limiting and destructive disease of potato is late blight, caused by Phytophtora infestans. Here, we present a multi-environmental analysis of the response to P. infestans using an association panel of 150 accessions of S. tuberosum Group Phureja, evaluated in two localities in Colombia. Disease resistance data were merged with a genotyping matrix of 83,862 SNPs obtained by 2b-restriction site-associated DNA and Genotyping by sequencing approaches into a Genome-wide association study. We are reporting 16 organ-specific QTL conferring resistance to late blight. These QTL explain between 13.7% and 50.9% of the phenotypic variance. Six and ten QTL were detected for resistance response in leaves and stem, respectively. In silico analysis revealed 15 candidate genes for resistance to late blight. Four of them have no functional genome annotation, while eleven candidate genes code for diverse proteins, including a leucine-rich repeat kinase.

Unraveling the molecules hidden in the gray shadows of quantitative disease resistance to pathogens

Una de las preguntas más desafiantes del fitomejoramiento y de la fitopatología molecular es ¿cuáles son las bases genéticas y moleculares de la resistencia cuantitativa a enfermedades?. El escaso conocimiento de cómo este tipo de resistencia funciona ha obstaculizado que los fitomejoradores la aprovecharlo plenamente. Para superar estos obstáculos se han desarrollado nuevas metodologías para el estudio de rasgos cuantitativos. Los enfoques como el mapeo genético, la identificación de loci de rasgos cuantitativos (QTL) y el mapeo por asociaciones, incluyendo el enfoque de genes candidatos y los estudios de asociación amplia del genoma, se han llevado a cabo históricamente para describir rasgos cuantitativos y por lo tanto para estudiar QDR. Además, se han proporcionado grandes avances en la obtención de datos fenotípicos cuantitativos para mejorar estos análisis. Recientemente, algunos genes asociados a QDR han sido clonados, lo que conduce a nuevas hipótesis sobre las bases moleculares de este tipo de resistencia. En esta revisión presentamos los avances más recientes sobre QDR y la correspondiente aplicación, que han permitido postular nuevas ideas que pueden ayudar a construir nuevos modelos. Algunas de las hipótesis presentadas aquí como posibles explicaciones para QDR están relacionadas con el nivel de expresión y el splicing alternativo de algunos genes relacionados con la defensa, la acción de "alelos débiles" de genes R, la presencia de variantes alélicas en los genes implicados en la respuesta de defensa y un papel central de quinasas o pseudoqinasas. Con la información recapitulada en esta revisión es posible concluir que la distinción conceptual entre resistencia cualitativa y cuantitativa puede ser cuestionada ya que ambos comparten importantes componentes.

Evolutionary Convergence and Divergence in NLR Function and Structure

The recognition of cellular damage caused by either pathogens or abiotic stress is essential for host defense in all forms of life in the plant and animal kingdoms. The NOD-like receptors (NLRs) represent a large family of multidomain proteins that were initially discovered for their role in host defense in plants and vertebrates. Over recent years the wide distribution of NLRs among metazoans has become apparent and their origins have begun to emerge. Moreover, intense study of NLR function has shown that they play essential roles beyond pathogen recognition - in the regulation of antigen presentation, cell death, inflammation, and even in embryonic development. We summarize here the latest insights into NLR biology and discuss examples of converging and diverging evolution of NLR function and structure.

Genome-Wide Association Studies In Plant Pathosystems: Toward an Ecological Genomics Approach

The emergence and re-emergence of plant pathogenic microorganisms are processes that imply perturbations in both host and pathogen ecological niches. Global change is largely assumed to drive the emergence of new etiological agents by altering the equilibrium of the ecological habitats which in turn places hosts more in contact with pathogen reservoirs. In this context, the number of epidemics is expected to increase dramatically in the next coming decades both in wild and crop plants. Under these considerations, the identification of the genetic variants underlying natural variation of resistance is a pre-requisite to estimate the adaptive potential of wild plant populations and to develop new breeding resistant cultivars. On the other hand, the prediction of pathogen's genetic determinants underlying disease emergence can help to identify plant resistance alleles. In the genomic era, whole genome sequencing combined with the development of statistical methods led to the emergence of Genome Wide Association (GWA) mapping, a powerful tool for detecting genomic regions associated with natural variation of disease resistance in both wild and cultivated plants. However, GWA mapping has been less employed for the detection of genetic variants associated with pathogenicity in microbes. Here, we reviewed GWA studies performed either in plants or in pathogenic microorganisms (bacteria, fungi and oomycetes). In addition, we highlighted the benefits and caveats of the emerging joint GWA mapping approach that allows for the simultaneous identification of genes interacting between genomes of both partners. Finally, based on co-evolutionary processes in wild populations, we highlighted a phenotyping-free joint GWA mapping approach as a promising tool for describing the molecular landscape underlying plant - microbe interactions.

Quantitative disease resistance to the bacterial pathogen Xanthomonas campestris involves an Arabidopsis immune receptor pair and a gene of unknown function

Although quantitative disease resistance (QDR) is a durable and broad-spectrum form of resistance in plants, the identification of the genes underlying QDR is still in its infancy. RKS1 (Resistance related KinaSe1) has been reported recently to confer QDR in Arabidopsis thaliana to most but not all races of the bacterial pathogen Xanthomonas campestris pv. campestris (Xcc). We therefore explored the genetic bases of QDR in A. thaliana to diverse races of X. campestris (Xc). A nested genome-wide association mapping approach was used to finely map the genomic regions associated with QDR to Xcc12824 (race 2) and XccCFBP6943 (race 6). To identify the gene(s) implicated in QDR, insertional mutants (T-DNA) were selected for the candidate genes and phenotyped in response to Xc. We identified two major QTLs that confer resistance specifically to Xcc12824 and XccCFBP6943. Although QDR to Xcc12824 is conferred by At5g22540 encoding for a protein of unknown function, QDR to XccCFBP6943 involves the well-known immune receptor pair RRS1/RPS4. In addition to RKS1, this study reveals that three genes are involved in resistance to Xc with strikingly different ranges of specificity, suggesting that QDR to Xc involves a complex network integrating multiple response pathways triggered by distinct pathogen molecular determinants.

The Genetics Underlying Natural Variation in the Biotic Interactions of Arabidopsis thaliana: The Challenges of Linking Evolutionary Genetics and Community Ecology

In the context of global change, predicting the responses of plant communities in an ever-changing biotic environment calls for a multipronged approach at the interface of evolutionary genetics and community ecology. However, our understanding of the genetic basis of natural variation involved in mediating biotic interactions, and associated adaptive dynamics of focal plants in their natural communities, is still in its infancy. Here, we review the genetic and molecular bases of natural variation in the response to biotic interactions (viruses, bacteria, fungi, oomycetes, herbivores, and plants) in the model plant Arabidopsis thaliana as well as the adaptive value of these bases. Among the 60 identified genes are a number that encode nucleotide-binding site leucine-rich repeat (NBS-LRR)-type proteins, consistent with early examples of plant defense genes. However, recent studies have revealed an extensive diversity in the molecular mechanisms of defense. Many types of genetic variants associate with phenotypic variation in biotic interactions, even among the genes of large effect that tend to be identified. In general, we found that (i) balancing selection rather than directional selection explains the observed patterns of genetic diversity within A. thaliana and (ii) the cost/benefit tradeoffs of adaptive alleles can be strongly dependent on both genomic and environmental contexts. Finally, because A. thaliana rarely interacts with only one biotic partner in nature, we highlight the benefit of exploring diffuse biotic interactions rather than tightly associated host-enemy pairs. This challenge would help to improve our understanding of coevolutionary quantitative genetics within the context of realistic community complexity.