Metaeffector interactions modulate the type III effector-triggered immunity load of Pseudomonas syringae.

The bacterial plant pathogen Pseudomonas syringae requires type III secreted effectors (T3SEs) for pathogenesis. However, a major facet of plant immunity entails the recognition of a subset of P. syringae's T3SEs by intracellular host receptors in a process called Effector-Triggered Immunity (ETI). Prior work has shown that ETI-eliciting T3SEs are pervasive in the P. syringae species complex raising the question of how P. syringae mitigates its ETI load to become a successful pathogen. While pathogens can evade ETI by T3SE mutation, recombination, or loss, there is increasing evidence that effector-effector (a.k.a., metaeffector) interactions can suppress ETI. To study the ETI-suppression potential of P. syringae T3SE repertoires, we compared the ETI-elicitation profiles of two genetically divergent strains: P. syringae pv. tomato DC3000 (PtoDC3000) and P. syringae pv. maculicola ES4326 (PmaES4326), which are both virulent on Arabidopsis thaliana but harbour largely distinct effector repertoires. Of the 529 T3SE alleles screened on A. thaliana Col-0 from the P. syringae T3SE compendium (PsyTEC), 69 alleles from 21 T3SE families elicited ETI in at least one of the two strain backgrounds, while 50 elicited ETI in both backgrounds, resulting in 19 differential ETI responses including two novel ETI-eliciting families: AvrPto1 and HopT1. Although most of these differences were quantitative, three ETI responses were completely absent in one of the pathogenic backgrounds. We performed ETI suppression screens to test if metaeffector interactions contributed to these ETI differences, and found that HopQ1a suppressed AvrPto1m-mediated ETI, while HopG1c and HopF1g suppressed HopT1b-mediated ETI. Overall, these results show that P. syringae strains leverage metaeffector interactions and ETI suppression to overcome the ETI load associated with their native T3SE repertoires.

Predictive modeling of Pseudomonas syringae virulence on bean using gradient boosted decision trees.

Pseudomonas syringae is a genetically diverse bacterial species complex responsible for numerous agronomically important crop diseases. Individual P. syringae isolates are assigned pathovar designations based on their host of isolation and the associated disease symptoms, and these pathovar designations are often assumed to reflect host specificity although this assumption has rarely been rigorously tested. Here we developed a rapid seed infection assay to measure the virulence of 121 diverse P. syringae isolates on common bean (Phaseolus vulgaris). This collection includes P. syringae phylogroup 2 (PG2) bean isolates (pathovar syringae) that cause bacterial spot disease and P. syringae phylogroup 3 (PG3) bean isolates (pathovar phaseolicola) that cause the more serious halo blight disease. We found that bean isolates in general were significantly more virulent on bean than non-bean isolates and observed no significant virulence difference between the PG2 and PG3 bean isolates. However, when we compared virulence within PGs we found that PG3 bean isolates were significantly more virulent than PG3 non-bean isolates, while there was no significant difference in virulence between PG2 bean and non-bean isolates. These results indicate that PG3 strains have a higher level of host specificity than PG2 strains. We then used gradient boosting machine learning to predict each strain's virulence on bean based on whole genome k-mers, type III secreted effector k-mers, and the presence/absence of type III effectors and phytotoxins. Our model performed best using whole genome data and was able to predict virulence with high accuracy (mean absolute error = 0.05). Finally, we functionally validated the model by predicting virulence for 16 strains and found that 15 (94%) had virulence levels within the bounds of estimated predictions. This study strengthens the hypothesis that P. syringae PG2 strains have evolved a different lifestyle than other P. syringae strains as reflected in their lower level of host specificity. It also acts as a proof-of-principle to demonstrate the power of machine learning for predicting host specific adaptation.

Pseudomonas syringae coffee blight is associated with the horizontal transfer of plasmid-encoded type III effectors.

The emergence of new pathogens is an ongoing threat to human health and agriculture. While zoonotic spillovers received considerable attention, the emergence of crop diseases is less well studied. Here, we identify genomic factors associated with the emergence of Pseudomonas syringae bacterial blight of coffee. Fifty-three P. syringae strains from diseased Brazilian coffee plants were sequenced. Comparative and evolutionary analyses were used to identify loci associated with coffee blight. Growth and symptomology assays were performed to validate the findings. Coffee isolates clustered in three lineages, including primary phylogroups PG3 and PG4, and secondary phylogroup PG11. Genome-wide association study of the primary PG strains identified 37 loci, including five effectors, most of which were encoded on a plasmid unique to the PG3 and PG4 coffee strains. Evolutionary analyses support the emergence of coffee blight in PG4 when the coffee-associated plasmid and associated effectors derived from a divergent plasmid carried by strains associated with other hosts. This plasmid was only recently transferred into PG3. Natural diversity and CRISPR-Cas9 plasmid curing were used to show that strains with the coffee-associated plasmid grow to higher densities and cause more severe disease symptoms in coffee. This work identifies possible evolutionary mechanisms underlying the emergence of a new lineage of coffee pathogens.

An extracellular network of Arabidopsis leucine-rich repeat receptor kinases.

The cells of multicellular organisms receive extracellular signals using surface receptors. The extracellular domains (ECDs) of cell surface receptors function as interaction platforms, and as regulatory modules of receptor activation. Understanding how interactions between ECDs produce signal-competent receptor complexes is challenging because of their low biochemical tractability. In plants, the discovery of ECD interactions is complicated by the massive expansion of receptor families, which creates tremendous potential for changeover in receptor interactions. The largest of these families in Arabidopsis thaliana consists of 225 evolutionarily related leucine-rich repeat receptor kinases (LRR-RKs), which function in the sensing of microorganisms, cell expansion, stomata development and stem-cell maintenance. Although the principles that govern LRR-RK signalling activation are emerging, the systems-level organization of this family of proteins is unknown. Here, to address this, we investigated 40,000 potential ECD interactions using a sensitized high-throughput interaction assay, and produced an LRR-based cell surface interaction network (CSILRR) that consists of 567 interactions. To demonstrate the power of CSILRR for detecting biologically relevant interactions, we predicted and validated the functions of uncharacterized LRR-RKs in plant growth and immunity. In addition, we show that CSILRR operates as a unified regulatory network in which the LRR-RKs most crucial for its overall structure are required to prevent the aberrant signalling of receptors that are several network-steps away. Thus, plants have evolved LRR-RK networks to process extracellular signals into carefully balanced responses.

Publisher Correction: An extracellular network of Arabidopsis leucine-rich repeat receptor kinases.

In this Letter, an incorrect version of the Supplementary Information file was inadvertently used, which contained several errors. The details of references 59-65 were missing from the end of the Supplementary Discussion section on page 4. In addition, the section 'Text 3. Y2H on ICD interactions' incorrectly referred to 'Extended Data Fig. 4d' instead of 'Extended Data Fig. 3d' on page 3. Finally, the section 'Text 4. Interaction network analysis' incorrectly referred to 'Fig. 1b and Extended Data Fig. 6' instead of 'Fig. 2b and Extended Data Fig. 7' on page 3. These errors have all been corrected in the Supplementary Information.

The small molecule Zaractin activates ZAR1-mediated immunity in Arabidopsis.

Pathogenic effector proteins use a variety of enzymatic activities to manipulate host cellular proteins and favor the infection process. However, these perturbations can be sensed by nucleotide-binding leucine-rich-repeat (NLR) proteins to activate effector-triggered immunity (ETI). Here we have identified a small molecule (Zaractin) that mimics the immune eliciting activity of the Pseudomonas syringae type III secreted effector (T3SE) HopF1r and show that both HopF1r and Zaractin activate the same NLR-mediated immune pathway in Arabidopsis Our results demonstrate that the ETI-inducing action of pathogenic effectors can be harnessed to identify synthetic activators of the eukaryotic immune system.

RecPD: A Recombination-aware measure of phylogenetic diversity.

A critical step in studying biological features (e.g., genetic variants, gene families, metabolic capabilities, or taxa) is assessing their diversity and distribution among a sample of individuals. Accurate assessments of these patterns are essential for linking features to traits or outcomes of interest and understanding their functional impact. Consequently, it is of crucial importance that the measures employed for quantifying feature diversity can perform robustly under any evolutionary scenario. However, the standard measures used for quantifying and comparing the distribution of features, such as prevalence, phylogenetic diversity, and related approaches, either do not take into consideration evolutionary history, or assume strictly vertical patterns of inheritance. Consequently, these approaches cannot accurately assess diversity for features that have undergone recombination or horizontal transfer. To address this issue, we have devised RecPD, a novel recombination-aware phylogenetic-diversity statistic for measuring the distribution and diversity of features under all evolutionary scenarios. RecPD utilizes ancestral-state reconstruction to map the presence / absence of features onto ancestral nodes in a species tree, and then identifies potential recombination events in the evolutionary history of the feature. We also derive several related measures from RecPD that can be used to assess and quantify evolutionary dynamics and correlation of feature evolutionary histories. We used simulation studies to show that RecPD reliably reconstructs feature evolutionary histories under diverse recombination and loss scenarios. We then applied RecPD in two diverse real-world scenarios including a preliminary study type III effector protein families secreted by the plant pathogenic bacterium Pseudomonas syringae and growth phenotypes of the Pseudomonas genus and demonstrate that prevalence is an inadequate measure that obscures the potential impact of recombination. We believe RecPD will have broad utility for revealing and quantifying complex evolutionary processes for features at any biological level.

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.

Diversity and Evolution of Type III Secreted Effectors: A Case Study of Three Families.

A broad range of Gram-negative bacteria employ a type III secretion system (T3SS) to deliver virulence proteins termed type III secreted effectors directly into the cytoplasm of eukaryotic host cells. While effectors can contribute to the colonization of eukaryotic hosts by bacterial symbionts and pathogens, they can also elicit host immune responses that restrict bacterial growth. These opposing selective pressures have shaped the evolution of effector families and may be responsible for their incredible diversity in biochemical function, mechanism of action, and taxonomic distribution. In this chapter, we focus on three distinct effector families whose members are distributed among both plant and animal pathogens. We first discuss the LRR-NEL and YopJ families of effectors. These two effector families possess ubiquitin ligase and acetyltransferase activity, respectively, which in both cases can be directed against host innate immune signal transduction pathways to promote infection. Finally, we discuss the TALE family of transcription activator-like effectors that serve to reprogram host immunity transcriptional responses. This chapter aims to highlight the diversity within these three effector families that results from the strong and dynamic evolutionary forces shaping the interface between host and bacterium.

A host-pathogen interactome uncovers phytopathogenic strategies to manipulate plant ABA responses.

The phytopathogen Pseudomonas syringae delivers into host cells type III secreted effectors (T3SEs) that promote virulence. One virulence mechanism employed by T3SEs is to target hormone signaling pathways to perturb hormone homeostasis. The phytohormone abscisic acid (ABA) influences interactions between various phytopathogens and their plant hosts, and has been shown to be a target of P. syringae T3SEs. In order to provide insight into how T3SEs manipulate ABA responses, we generated an ABA-T3SE interactome network (ATIN) between P. syringae T3SEs and Arabidopsis proteins encoded by ABA-regulated genes. ATIN consists of 476 yeast-two-hybrid interactions between 97 Arabidopsis ABA-regulated proteins and 56 T3SEs from four pathovars of P. syringae. We demonstrate that T3SE interacting proteins are significantly enriched for proteins associated with transcription. In particular, the ETHYLENE RESPONSIVE FACTOR (ERF) family of transcription factors is highly represented. We show that ERF105 and ERF8 displayed a role in defense against P. syringae, supporting our overall observation that T3SEs of ATIN converge on proteins that influence plant immunity. In addition, we demonstrate that T3SEs that interact with a large number of ABA-regulated proteins can influence ABA responses. One of these T3SEs, HopF3Pph6 , inhibits the function of ERF8, which influences both ABA-responses and plant immunity. These results provide a potential mechanism for how HopF3Pph6 manipulates ABA-responses to promote P. syringae virulence, and also demonstrate the utility of ATIN as a resource to study the ABA-T3SE interface.

A High-Throughput, Seedling Screen for Plant Immunity.

An understanding of how biological diversity affects plant-microbe interactions is becoming increasingly important, particularly with respect to components of the pathogen effector arsenal and the plant immune system. Although technological improvements have greatly advanced our ability to examine molecular sequences and interactions, relatively few advances have been made that facilitate high-throughput, in vivo pathology screens. Here, we present a high-throughput, microplate-based, nondestructive seedling pathology assay, and apply it to identify Arabidopsis thaliana effector-triggered immunity (ETI) responses against Pseudomonas syringae type III secreted effectors. The assay was carried out in a 48-well microplate format with spray inoculation, and disease symptoms were quantitatively recorded in a semiautomated manner, thereby greatly reducing both time and costs. The assay requires only slight modifications of common labware and uses no proprietary software. We validated the assay by recapitulating known ETI responses induced by P. syringae in Arabidopsis. We also demonstrated that we can quantitatively differentiate responses from a diversity of plant genotypes grown in the same microplate. Finally, we showed that the results obtained from our assay can be used to perform genome-wide association studies to identify host immunity genes, recapitulating results that have been independently obtained with mature plants.

Oh, the places they'll go! A survey of phytopathogen effectors and their host targets.

Phytopathogens translocate effector proteins into plant cells where they sabotage the host cellular machinery to promote infection. An individual pathogen can translocate numerous distinct effectors during the infection process to target an array of host macromolecules (proteins, metabolites, DNA, etc.) and manipulate them using a variety of enzymatic activities. In this review, we have surveyed the literature for effector targets and curated them to convey the range of functions carried out by phytopathogenic proteins inside host cells. In particular, we have curated the locations of effector targets, as well as their biological and molecular functions and compared these properties across diverse phytopathogens. This analysis validates previous observations about effector functions (e.g. immunosuppression), and also highlights some interesting features regarding effector specificity as well as functional diversification of phytopathogen virulence strategies.

A High-Sensitivity, Microtiter-Based Plate Assay for Plant Pattern-Triggered Immunity.

The first step in the plant immune response to pathogen challenge involves the perception of conserved epitopes, called microbe-associated molecular patterns (MAMPs), by cell-surface pattern recognition receptors (PRRs). Given the key roles that MAMPs and PRRs play in plant innate immunity, great effort has been expended to identify these molecules. Current methods for assaying these immune responses are often limited in their resolution and throughput, and consequently, there is a need for medium- to high-throughput methodologies. Here, we describe the development of a 96-well microtiter plate-based assay for plant pattern-triggered immunity that measures the activity of plant peroxidase (POX) enzymes produced in response to treatment with bacterial MAMPs. The system has been optimized to minimize both the amount of plant tissue and MAMPs required and displays up to three orders of magnitude greater sensitivity than the traditional luminol-based reactive oxygen species assay when measuring the plant response to treatment with the bacterial MAMP flg22, reaching detection limits in the picomolar range. This high sensitivity opens the possibility of evaluating the immune-eliciting effects of weaker elicitors. The throughput and material requirements of the assay make it ideal for screens involving quantitative measurement of the plant innate immune response to MAMPs.

Arabidopsis ETHYLENE RESPONSE FACTOR 8 (ERF8) has dual functions in ABA signaling and immunity.

BACKGROUND: ETHYLENE RESPONSE FACTOR (ERF) 8 is a member of one of the largest transcription factor families in plants, the APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF) superfamily. Members of this superfamily have been implicated in a wide variety of processes such as development and environmental stress responses. RESULTS: In this study we demonstrated that ERF8 is involved in both ABA and immune signaling. ERF8 overexpression induced programmed cell death (PCD) in Arabidopsis and Nicotiana benthamiana. This PCD was salicylic acid (SA)-independent, suggesting that ERF8 acts downstream or independent of SA. ERF8-induced PCD was abolished by mutations within the ERF-associated amphiphilic repression (EAR) motif, indicating ERF8 induces cell death through its transcriptional repression activity. Two immunity-related mitogen-activated protein kinases, MITOGEN-ACTIVATED PROTEIN KINASE 4 (MPK4) and MPK11, were identified as ERF8-interacting proteins and directly phosphorylated ERF8 in vitro. Four putative MPK phosphorylation sites were identified in ERF8, one of which (Ser103) was determined to be the predominantly phosphorylated residue in vitro, while mutation of all four putative phosphorylation sites partially suppressed ERF8-induced cell death in N. benthamiana. Genome-wide transcriptomic analysis and pathogen growth assays confirmed a positive role of ERF8 in mediating immunity, as ERF8 knockdown or overexpression lines conferred compromised or enhanced resistance against the hemibiotrophic bacterial pathogen Pseudomonas syringae, respectively. CONCLUSIONS: Together these data reveal that the ABA-inducible transcriptional repressor ERF8 has dual roles in ABA signaling and pathogen defense, and further highlight the complex influence of ABA on plant-microbe interactions.

New BAR tools for mining expression data and exploring Cis-elements in Arabidopsis thaliana.

Identifying sets of genes that are specifically expressed in certain tissues or in response to an environmental stimulus is useful for designing reporter constructs, generating gene expression markers, or for understanding gene regulatory networks. We have developed an easy-to-use online tool for defining a desired expression profile (a modification of our Expression Angler program), which can then be used to identify genes exhibiting patterns of expression that match this profile as closely as possible. Further, we have developed another online tool, Cistome, for predicting or exploring cis-elements in the promoters of sets of co-expressed genes identified by such a method, or by other methods. We present two use cases for these tools, which are freely available on the Bio-Analytic Resource at http://BAR.utoronto.ca.

The ERF11 Transcription Factor Promotes Internode Elongation by Activating Gibberellin Biosynthesis and Signaling.

The phytohormone gibberellin (GA) plays a key role in promoting stem elongation in plants. Previous studies show that GA activates its signaling pathway by inducing rapid degradation of DELLA proteins, GA signaling repressors. Using an activation-tagging screen in a reduced-GA mutant ga1-6 background, we identified AtERF11 to be a novel positive regulator of both GA biosynthesis and GA signaling for internode elongation. Overexpression of AtERF11 partially rescued the dwarf phenotype of ga1-6 AtERF11 is a member of the ERF (ETHYLENE RESPONSE FACTOR) subfamily VIII-B-1a of ERF/AP2 transcription factors in Arabidopsis (Arabidopsis thaliana). Overexpression of AtERF11 resulted in elevated bioactive GA levels by up-regulating expression of GA3ox1 and GA20ox genes. Hypocotyl elongation assays further showed that overexpression of AtERF11 conferred elevated GA response, whereas loss-of-function erf11 and erf11 erf4 mutants displayed reduced GA response. In addition, yeast two-hybrid, coimmunoprecipitation, and transient expression assays showed that AtERF11 enhances GA signaling by antagonizing the function of DELLA proteins via direct protein-protein interaction. Interestingly, AtERF11 overexpression also caused a reduction in the levels of another phytohormone ethylene in the growing stem, consistent with recent finding showing that AtERF11 represses transcription of ethylene biosynthesis ACS genes. The effect of AtERF11 on promoting GA biosynthesis gene expression is likely via its repressive function on ethylene biosynthesis. These results suggest that AtERF11 plays a dual role in promoting internode elongation by inhibiting ethylene biosynthesis and activating GA biosynthesis and signaling pathways.

Image-Based Quantification of Plant Immunity and Disease.

Measuring the extent and severity of disease is a critical component of plant pathology research and crop breeding. Unfortunately, existing visual scoring systems are qualitative, subjective, and the results are difficult to transfer between research groups, while existing quantitative methods can be quite laborious. Here, we present plant immunity and disease image-based quantification (PIDIQ), a quantitative, semi-automated system to rapidly and objectively measure disease symptoms in a biologically relevant context. PIDIQ applies an ImageJ-based macro to plant photos in order to distinguish healthy tissue from tissue that has yellowed due to disease. It can process a directory of images in an automated manner and report the relative ratios of healthy to diseased leaf area, thereby providing a quantitative measure of plant health that can be statistically compared with appropriate controls. We used the Arabidopsis thaliana-Pseudomonas syringae model system to show that PIDIQ is able to identify both enhanced plant health associated with effector-triggered immunity as well as elevated disease symptoms associated with effector-triggered susceptibility. Finally, we show that the quantitative results provided by PIDIQ correspond to those obtained via traditional in planta pathogen growth assays. PIDIQ provides a simple and effective means to nondestructively quantify disease from whole plants and we believe it will be equally effective for monitoring disease on excised leaves and stems.

The Arabidopsis ZED1 pseudokinase is required for ZAR1-mediated immunity induced by the Pseudomonas syringae type III effector HopZ1a.

Plant and animal pathogenic bacteria can suppress host immunity by injecting type III secreted effector (T3SE) proteins into host cells. However, T3SEs can also elicit host immunity if the host has evolved a means to recognize the presence or activity of specific T3SEs. The diverse YopJ/HopZ/AvrRxv T3SE superfamily, which is found in both animal and plant pathogens, provides examples of T3SEs playing this dual role. The T3SE HopZ1a is an acetyltransferase carried by the phytopathogen Pseudomonas syringae that elicits effector-triggered immunity (ETI) when recognized in Arabidopsis thaliana by the nucleotide-binding leucine-rich repeat (NB-LRR) protein ZAR1. However, recognition of HopZ1a does not require any known ETI-related genes. Using a forward genetics approach, we identify a unique ETI-associated gene that is essential for ZAR1-mediated immunity. The hopZ-ETI-deficient1 (zed1) mutant is specifically impaired in the recognition of HopZ1a, but not the recognition of other unrelated T3SEs or in pattern recognition receptor (PRR)-triggered immunity. ZED1 directly interacts with both HopZ1a and ZAR1 and is acetylated on threonines 125 and 177 by HopZ1a. ZED1 is a nonfunctional kinase that forms part of small genomic cluster of kinases in Arabidopsis. We hypothesize that ZED1 acts as a decoy to lure HopZ1a to the ZAR1-resistance complex, resulting in ETI activation.

A bacterial acetyltransferase destroys plant microtubule networks and blocks secretion.

The eukaryotic cytoskeleton is essential for structural support and intracellular transport, and is therefore a common target of animal pathogens. However, no phytopathogenic effector has yet been demonstrated to specifically target the plant cytoskeleton. Here we show that the Pseudomonas syringae type III secreted effector HopZ1a interacts with tubulin and polymerized microtubules. We demonstrate that HopZ1a is an acetyltransferase activated by the eukaryotic co-factor phytic acid. Activated HopZ1a acetylates itself and tubulin. The conserved autoacetylation site of the YopJ / HopZ superfamily, K289, plays a critical role in both the avirulence and virulence function of HopZ1a. Furthermore, HopZ1a requires its acetyltransferase activity to cause a dramatic decrease in Arabidopsis thaliana microtubule networks, disrupt the plant secretory pathway and suppress cell wall-mediated defense. Together, this study supports the hypothesis that HopZ1a promotes virulence through cytoskeletal and secretory disruption.