Solution NMR structures of Pyrenophora tritici-repentis ToxB and its inactive homolog reveal potential determinants of toxin activity

Evolution of the Toxb Gene in Pyrenophora tritici-repentis and Related Species

Pyrenophora tritici-repentis (tan spot) is a destructive foliar pathogen of wheat with global impact. This ascomycete fungus possesses a highly plastic open pangenome shaped by the gain and loss of effector genes. This study investigated the allelic variations in the chlorosis-encoding gene ToxB across 422 isolates representing all identified pathotypes and worldwide origins. To gain better insights into ToxB evolution, we examined its presence and variability in other Pyrenophora spp. A ToxB haplotype network was constructed, revealing the evolutionary relationships of this gene (20 haplotypes) across four Pyrenophora species. Notably, toxb, the homolog of ToxB, was detected for the first time in the barley pathogen Pyrenophora teres. The ToxB/toxb genes display evidence of selection that is characterized by loss of function, duplication, and diverse mutations. Within the ToxB/toxb open reading frame, 72 mutations were identified, including 14 synonymous, 55 nonsynonymous, and 3 indel mutations. Remarkably, a, ∼5.6-kb Copia-like retrotransposon, named Copia-1_Ptr, was found inserted in the toxb gene of a race 3 isolate. This insert disrupted the ToxB gene's function, a first case of effector gene disruption by a transposable element in P. tritici-repentis. Additionally, a microsatellite with 25 nucleotide repeats (0 to 10) in the upstream region of ToxB suggested a potential mechanism influencing ToxB expression and regulation. Exploring ToxB-like protein distribution in other ascomycetes revealed the presence of ToxB-like proteins in 19 additional species, including the Leotiomycetes class for the first time. The presence/absence pattern of ToxB-like proteins defied species relatedness compared with a phylogenetic tree, suggesting a past horizontal gene transfer event during the evolution of the ToxB gene. [Formula: see text] Copyright © 2024 His Majesty the King in Right of Canada, as represented by the Minister of Agriculture and Agri-Food. This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

The Molecular Docking of MAX Fungal Effectors with Plant HMA Domain-Binding Proteins

Fungal effector proteins are important in mediating disease infections in agriculturally important crops. These secreted small proteins are known to interact with their respective host receptor binding partners in the host, either inside the cells or in the apoplastic space, depending on the localisation of the effector proteins. Consequently, it is important to understand the interactions between fungal effector proteins and their target host receptor binding partners, particularly since this can be used for the selection of potential plant resistance or susceptibility-related proteins that can be applied to the breeding of new cultivars with disease resistance. In this study, molecular docking simulations were used to characterise protein-protein interactions between effector and plant receptors. Benchmarking was undertaken using available experimental structures of effector-host receptor complexes to optimise simulation parameters, which were then used to predict the structures and mediating interactions of effector proteins with host receptor binding partners that have not yet been characterised experimentally. Rigid docking was applied for both the so-called bound and unbound docking of MAX effectors with plant HMA domain protein partners. All bound complexes used for benchmarking were correctly predicted, with 84% being ranked as the top docking pose using the ZDOCK scoring function. In the case of unbound complexes, a minimum of 95% of known residues were predicted to be part of the interacting interface on the host receptor binding partner, and at least 87% of known residues were predicted to be part of the interacting interface on the effector protein. Hydrophobic interactions were found to dominate the formation of effector-plant protein complexes. An optimised set of docking parameters based on the use of ZDOCK and ZRANK scoring functions were established to enable the prediction of near-native docking poses involving different binding interfaces on plant HMA domain proteins. Whilst this study was limited by the availability of the experimentally determined complexed structures of effectors and host receptor binding partners, we demonstrated the potential of molecular docking simulations to predict the likely interactions between effectors and their respective host receptor binding partners. This computational approach may accelerate the process of the discovery of putative interacting plant partners of effector proteins and contribute to effector-assisted marker discovery, thereby supporting the breeding of disease-resistant crops.

Does a Similar 3D Structure Mean a Similar Folding Pathway? The Presence of a C-Terminal α-Helical Extension in the 3D Structure of MAX60 Drastically Changes the Folding Pathway Described for Other MAX-Effectors from Magnaporthe oryzae

Does a similar 3D structure mean a similar folding pathway? This question is particularly meaningful when it concerns proteins sharing a similar 3D structure, but low sequence identity or homology. MAX effectors secreted by the phytopathogenic fungus Magnaporthe oryzae present such characteristics. They share a common 3D structure, a ß-sandwich with the same topology for all the family members, but an extremely low sequence identity/homology. In a previous study, we have investigated the folding of two MAX effectors, AVR-Pia and AVR-Pib, using High-Hydrostatic-Pressure NMR and found that they display a similar folding pathway, with a common folding intermediate. In the present work, we used a similar strategy to investigate the folding conformational landscape of another MAX effector, MAX60, and found a very different folding intermediate. Our analysis strongly supports that the presence of a C-terminal α-helical extension in the 3D structure of MAX60 could be responsible for its different folding pathway.

A pathogen effector FOLD diversified in symbiotic fungi

Pathogenic fungi use secreted effector proteins to suppress immunity and support their infection, but effectors have also been reported from fungi that engage in nutritional symbioses with plants. Sequence-based effector comparisons between pathogens and symbiotic arbuscular mycorrhizal (AM) fungi are hampered by the huge diversity of effector sequences even within closely related microbes. To find sequence-divergent but structurally similar effectors shared between symbiotic and pathogenic fungi, we compared secreted protein structure models of the AM fungus Rhizophagus irregularis to known pathogen effectors. We identified proteins with structural similarity to known Fusarium oxysporum f. sp. lycopersici dual domain (FOLD) effectors, which occur in low numbers in several fungal pathogens. Contrastingly, FOLD genes from AM fungi (MycFOLDs) are found in enlarged and diversified gene families with higher levels of positive selection in their C-terminal domains. Our structure model comparison suggests that MycFOLDs are similar to carbohydrate-binding motifs. Different MycFOLD genes are expressed during colonisation of different hosts and MycFOLD-17 transcripts accumulate in plant intracellular arbuscules. The exclusive presence of MycFOLDs across unrelated plant-colonising fungi, their inducible expression, lineage-specific sequence diversification and transcripts in arbuscules suggest that FOLD proteins act as effectors during plant colonisation of symbiotic and pathogenic fungi.

Ab Initio Modelling of the Structure of ToxA-like and MAX Fungal Effector Proteins

Pathogenic fungal diseases in crops are mediated by the release of effector proteins that facilitate infection. Characterising the structure of these fungal effectors is vital to understanding their virulence mechanisms and interactions with their hosts, which is crucial in the breeding of plant cultivars for disease resistance. Several effectors have been identified and validated experimentally; however, their lack of sequence conservation often impedes the identification and prediction of their structure using sequence similarity approaches. Structural similarity has, nonetheless, been observed within fungal effector protein families, creating interest in validating the use of computational methods to predict their tertiary structure from their sequence. We used Rosetta ab initio modelling to predict the structures of members of the ToxA-like and MAX effector families for which experimental structures are known to validate this method. An optimised approach was then used to predict the structures of phenotypically validated effectors lacking known structures. Rosetta was found to successfully predict the structure of fungal effectors in the ToxA-like and MAX families, as well as phenotypically validated but structurally unconfirmed effector sequences. Interestingly, potential new effector structural families were identified on the basis of comparisons with structural homologues and the identification of associated protein domains.

The Venturia inaequalis effector repertoire is dominated by expanded families with predicted structural similarity, but unrelated sequence, to avirulence proteins from other plant-pathogenic fungi

Background

Scab, caused by the biotrophic fungus Venturia inaequalis, is the most economically important disease of apples worldwide. During infection, V. inaequalis occupies the subcuticular environment, where it secretes virulence factors, termed effectors, to promote host colonization. Consistent with other plant-pathogenic fungi, many of these effectors are expected to be non-enzymatic proteins, some of which can be recognized by corresponding host resistance proteins to activate plant defences, thus acting as avirulence determinants. To develop durable control strategies against scab, a better understanding of the roles that these effector proteins play in promoting subcuticular growth by V. inaequalis, as well as in activating, suppressing, or circumventing resistance protein-mediated defences in apple, is required.

Results

We generated the first comprehensive RNA-seq transcriptome of V. inaequalis during colonization of apple. Analysis of this transcriptome revealed five temporal waves of gene expression that peaked during early, mid, or mid-late infection. While the number of genes encoding secreted, non-enzymatic proteinaceous effector candidates (ECs) varied in each wave, most belonged to waves that peaked in expression during mid-late infection. Spectral clustering based on sequence similarity determined that the majority of ECs belonged to expanded protein families. To gain insights into function, the tertiary structures of ECs were predicted using AlphaFold2. Strikingly, despite an absence of sequence similarity, many ECs were predicted to have structural similarity to avirulence proteins from other plant-pathogenic fungi, including members of the MAX, LARS, ToxA and FOLD effector families. In addition, several other ECs, including an EC family with sequence similarity to the AvrLm6 avirulence effector from Leptosphaeria maculans, were predicted to adopt a KP6-like fold. Thus, proteins with a KP6-like fold represent another structural family of effectors shared among plant-pathogenic fungi.

Conclusions

Our study reveals the transcriptomic profile underpinning subcuticular growth by V. inaequalis and provides an enriched list of ECs that can be investigated for roles in virulence and avirulence. Furthermore, our study supports the idea that numerous sequence-unrelated effectors across plant-pathogenic fungi share common structural folds. In doing so, our study gives weight to the hypothesis that many fungal effectors evolved from ancestral genes through duplication, followed by sequence diversification, to produce sequence-unrelated but structurally similar proteins.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01442-9.

1H, 13C, 15 N backbone and side-chain NMR assignments for three MAX effectors from Magnaporthe oryzae

Effectors are small and very diverse proteins secreted by fungi and translocated in plant cells during infection. Among them, MAX effectors (for Magnaporthe Avrs and ToxB) were identified as a family of effectors that share an identical fold topology despite having highly divergent sequences. They are mostly secreted by ascomycetes from the Magnaporthe genus, a fungus that causes the rice blast, a plant disease leading to huge crop losses. As rice is the first source of calories in many countries, especially in Asia and Africa, this constitutes a threat for world food security. Hence, a better understanding of these effectors, including structural and functional characterization, constitutes a strategic milestone in the fight against phytopathogen fungi and may give clues for the development of resistant varieties of rice. We report here the near complete ¹H, ¹⁵ N and ¹³C NMR resonance assignment of three new putative MAX effectors (MAX47, MAX60 and MAX67). Secondary structure determination using TALOS-N and CSI.3 demonstrates a high content of β-strands in all the three proteins, in agreement with the canonic ß-sandwich structure of MAX effectors. This preliminary study provides foundations for further structural characterization, that will help in turn to improve sequence predictions of other MAX effectors through data mining.

Seeing is believing: Exploiting advances in structural biology to understand and engineer plant immunity

Filamentous plant pathogens cause disease in numerous economically important crops. These pathogens secrete virulence proteins, termed effectors, that modulate host cellular processes and promote infection. Plants have evolved immunity receptors that detect effectors and activate defence pathways, resulting in resistance to the invading pathogen. This leads to an evolutionary arms race between pathogen and host that is characterised by highly diverse effector repertoires in plant pathogens. Here, we review the recent advances in understanding host-pathogen co-evolution provided by the structural determination of effectors alone, and in complex with immunity receptors. We highlight the use of recent advances in structural prediction within this field and its role for future development of designer resistance proteins.

A Conserved Hypothetical Gene Is Required but Not Sufficient for Ptr ToxC Production in Pyrenophora tritici-repentis

The fungus Pyrenophora tritici-repentis causes tan spot, an important foliar disease of wheat worldwide. The fungal pathogen produces three necrotrophic effectors, namely Ptr ToxA, Ptr ToxB, and Ptr ToxC to induce necrosis or chlorosis in wheat. Both Ptr ToxA and Ptr ToxB are proteins, and their encoding genes have been cloned. Ptr ToxC was characterized as a low-molecular weight molecule 20 years ago but the one or more genes controlling its production in P. tritici-repentis are unknown. Here, we report the genetic mapping, molecular cloning, and functional analysis of a fungal gene that is required for Ptr ToxC production. The genetic locus controlling the production of Ptr ToxC, termed ToxC, was mapped to a subtelomeric region using segregating biparental populations, genome sequencing, and association analysis. Additional marker analysis further delimited ToxC to a 173-kb region. The predicted genes in the region were examined for presence/absence polymorphism in different races and isolates leading to the identification of a single candidate gene. Functional validation showed that this gene was required but not sufficient for Ptr ToxC production, thus it is designated as ToxC1. ToxC1 encoded a conserved hypothetical protein likely located on the vacuole membrane. The gene was highly expressed during infection, and only one haplotype was identified among 120 isolates sequenced. Our work suggests that Ptr ToxC is not a protein and is likely produced through a cascade of biosynthetic pathway. The identification of ToxC1 is a major step toward revealing the Ptr ToxC biosynthetic pathway and studying its molecular interactions with host factors.[Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Ca2+ sensor-mediated ROS scavenging suppresses rice immunity and is exploited by a fungal effector

Plant immunity is activated upon pathogen perception and often affects growth and yield when it is constitutively active. How plants fine-tune immune homeostasis in their natural habitats remains elusive. Here, we discover a conserved immune suppression network in cereals that orchestrates immune homeostasis, centering on a Ca²⁺-sensor, RESISTANCE OF RICE TO DISEASES1 (ROD1). ROD1 promotes reactive oxygen species (ROS) scavenging by stimulating catalase activity, and its protein stability is regulated by ubiquitination. ROD1 disruption confers resistance to multiple pathogens, whereas a natural ROD1 allele prevalent in indica rice with agroecology-specific distribution enhances resistance without yield penalty. The fungal effector AvrPiz-t structurally mimics ROD1 and activates the same ROS-scavenging cascade to suppress host immunity and promote virulence. We thus reveal a molecular framework adopted by both host and pathogen that integrates Ca²⁺ sensing and ROS homeostasis to suppress plant immunity, suggesting a principle for breeding disease-resistant, high-yield crops.

Remote homology clustering identifies lowly conserved families of effector proteins in plant-pathogenic fungi

Plant diseases caused by fungal pathogens are typically initiated by molecular interactions between 'effector' molecules released by a pathogen and receptor molecules on or within the plant host cell. In many cases these effector-receptor interactions directly determine host resistance or susceptibility. The search for fungal effector proteins is a developing area in fungal-plant pathology, with more than 165 distinct confirmed fungal effector proteins in the public domain. For a small number of these, novel effectors can be rapidly discovered across multiple fungal species through the identification of known effector homologues. However, many have no detectable homology by standard sequence-based search methods. This study employs a novel comparison method (RemEff) that is capable of identifying protein families with greater sensitivity than traditional homology-inference methods, leveraging a growing pool of confirmed fungal effector data to enable the prediction of novel fungal effector candidates by protein family association. Resources relating to the RemEff method and data used in this study are available from https://figshare.com/projects/Effector_protein_remote_homology/87965.

The crystal structure of SnTox3 from the necrotrophic fungus Parastagonospora nodorum reveals a unique effector fold and provides insight into Snn3 recognition and pro-domain protease processing of fungal effectors

Plant pathogens cause disease through secreted effector proteins, which act to promote infection. Typically, the sequences of effectors provide little functional information and further targeted experimentation is required. Here, we utilized a structure/function approach to study SnTox3, an effector from the necrotrophic fungal pathogen Parastagonospora nodorum, which causes cell death in wheat-lines carrying the sensitivity gene Snn3. We developed a workflow for the production of SnTox3 in a heterologous host that enabled crystal structure determination and functional studies. We show this approach can be successfully applied to study effectors from other pathogenic fungi. The β-barrel fold of SnTox3 is a novel fold among fungal effectors. Structure-guided mutagenesis enabled the identification of residues required for Snn3 recognition. SnTox3 is a pre-pro-protein, and the pro-domain of SnTox3 can be cleaved in vitro by the protease Kex2. Complementing this, an in silico study uncovered the prevalence of a conserved motif (LxxR) in an expanded set of putative pro-domain-containing fungal effectors, some of which can be cleaved by Kex2 in vitro. Our in vitro and in silico study suggests that Kex2-processed pro-domain (designated here as K2PP) effectors are common in fungi and this may have broad implications for the approaches used to study their functions.

The crystal structure of SnTox3 from the necrotrophic fungus Parastagonospora nodorum reveals a unique effector fold and provides insight into Snn3 recognition and pro-domain protease processing of fungal effectors

Plant pathogens cause disease through secreted effector proteins, which act to promote infection. Typically, the sequences of effectors provide little functional information and further targeted experimentation is required. Here, we utilized a structure/function approach to study SnTox3, an effector from the necrotrophic fungal pathogen Parastagonospora nodorum, which causes cell death in wheat-lines carrying the sensitivity gene Snn3. We developed a workflow for the production of SnTox3 in a heterologous host that enabled crystal structure determination and functional studies. We show this approach can be successfully applied to study effectors from other pathogenic fungi. The β-barrel fold of SnTox3 is a novel fold among fungal effectors. Structure-guided mutagenesis enabled the identification of residues required for Snn3 recognition. SnTox3 is a pre-pro-protein, and the pro-domain of SnTox3 can be cleaved in vitro by the protease Kex2. Complementing this, an in silico study uncovered the prevalence of a conserved motif (LxxR) in an expanded set of putative pro-domain-containing fungal effectors, some of which can be cleaved by Kex2 in vitro. Our in vitro and in silico study suggests that Kex2-processed pro-domain (designated here as K2PP) effectors are common in fungi and this may have broad implications for the approaches used to study their functions.

Characterization of Effector-Target Interactions in Necrotrophic Pathosystems Reveals Trends and Variation in Host Manipulation

Great strides have been made in defining the details of the plant defense response involving biotrophic fungal and bacterial pathogens. The groundwork for the current model was laid by H.H. Flor and others who defined the gene-for-gene hypothesis, which is now known to involve effector-triggered immunity (ETI). PAMP-triggered immunity (PTI) is also a highly effective response to most pathogens because of the recognition of common pathogen molecules by pattern recognition receptors. In this article, we consider the three pathogens that make up the foliar disease complex of wheat, Zymoseptoria tritici, Pyrenophora tritici-repentis, and Parastagonospora nodorum, to review the means by which necrotrophic pathogens circumvent, or outright hijack, the ETI and PTI pathways to cause disease.

Tactics of host manipulation by intracellular effectors from plant pathogenic fungi

Fungal pathogens can secrete hundreds of effectors, some of which are known to promote host susceptibility. This biological complexity, together with the lack of genetic tools in some fungi, presents a substantial challenge to develop a broad picture of the mechanisms these pathogens use for host manipulation. Nevertheless, recent advances in understanding individual effector functions are beginning to flesh out our view of fungal pathogenesis. This review discusses some of the latest findings that illustrate how effectors from diverse species use similar strategies to modulate plant physiology to their advantage. We also summarize recent breakthroughs in the identification of effectors from challenging systems, like obligate biotrophs, and emerging concepts such as the 'iceberg model' to explain how the activation of plant immunity can be turned off by effectors with suppressive activity.

Deciphering the Infectious Process of Colletotrichum lupini in Lupin through Transcriptomic and Proteomic Analysis

The fungal phytopathogen Colletotrichum lupini is responsible for lupin anthracnose, resulting in significant yield losses worldwide. The molecular mechanisms underlying this infectious process are yet to be elucidated. This study proposes to evaluate C. lupini gene expression and protein synthesis during lupin infection, using, respectively, an RNAseq-based transcriptomic approach and a mass spectrometry-based proteomic approach. Patterns of differentially-expressed genes in planta were evaluated from 24 to 84 hours post-inoculation, and compared to in vitro cultures. A total of 897 differentially-expressed genes were identified from C. lupini during interaction with white lupin, of which 520 genes were predicted to have a putative function, including carbohydrate active enzyme, effector, protease or transporter-encoding genes, commonly described as pathogenicity factors for other Colletotrichum species during plant infection, and 377 hypothetical proteins. Simultaneously, a total of 304 proteins produced during the interaction were identified and quantified by mass spectrometry. Taken together, the results highlight that the dynamics of symptoms, gene expression and protein synthesis shared similarities to those of hemibiotrophic pathogens. In addition, a few genes with unknown or poorly-described functions were found to be specifically associated with the early or late stages of infection, suggesting that they may be of importance for pathogenicity. This study, conducted for the first time on a species belonging to the Colletotrichum acutatum species complex, presents an opportunity to deepen functional analyses of the genes involved in the pathogenicity of Colletotrichum spp. during the onset of plant infection.

Host-specificity factors in plant pathogenic fungi

Fortunately, no fungus can cause disease on all plant species, and although some plant-pathogenic fungi have quite a broad host range, most are highly limited in the range of plant species or even cultivars that they cause disease in. The mechanisms of host specificity have been extensively studied in many plant-pathogenic fungi, especially in fungal pathogens causing disease on economically important crops. Specifically, genes involved in host specificity have been identified during the last few decades. In this overview, we describe and discuss these host-specificity genes. These genes encode avirulence (Avr) proteins, proteinaceous host-specific toxins or secondary metabolites. We discuss the genomic context of these genes, their expression, polymorphism, horizontal transfer and involvement in pathogenesis.

Necrotrophic Exploitation and Subversion of Plant Defense: A Lifestyle or Just a Phase, and Implications in Breeding Resistance

Breeding disease-resistant plants is a critical, environmentally friendly component of any strategy to sustainably feed and clothe the 9.8 billion people expected to live on Earth by 2050. Here, I review current literature detailing plant defense responses as they relate to diverse biological outcomes; disease resistance, susceptibility, and establishment of mutualistic plant-microbial relationships. Of particular interest is the degree to which these outcomes are a function of plant-associated microorganisms' lifestyles; biotrophic, hemibiotrophic, necrotrophic, or mutualistic. For the sake of brevity, necrotrophic pathogens and the necrotrophic phase of pathogenicity are emphasized in this review, with special attention given to the host-specific pathogens that exploit defense. Defense responses related to generalist necrotrophs and mutualists are discussed in the context of excellent reviews by others. In addition, host evolutionary trade-offs of disease resistance with other desirable traits are considered in the context of breeding for durable disease resistance.

Heterologous Expression of the Pyrenophora tritici-repentis Effector Proteins ToxA and ToxB, and the Prevalence of Effector Sensitivity in Australian Cereal Crops

Here, we evaluate the expression of the proteinaceous effectors ToxA and ToxB, produced by the necrotrophic fungal pathogen Pyrenophora tritici-repentis, which confer tan spot disease susceptibility on wheat. These necrotrophic effectors were expressed in two heterologous systems: Escherichia coli and Pichia pastoris. The E. coli SHuffle system was demonstrated to be superior to P. pastoris in generating high-levels of recombinant proteins that were soluble and stable. In addition, protein extracts from P. pastoris induced non-specific chlorosis on wheat, postulated to be caused by co-purified glucanases secreted by the host. Up to 79.6 μg/ml of ToxB was obtained using the SHuffle system in the absence of the native signal peptide, whilst the ToxA yield was considerably lower at 3.2 μg/ml. Results indicated that a histidine tag at the ToxA C-terminus interfered with effector functionality. Heterologously expressed ToxA and ToxB were tested on a panel of Australian cereals, including 122 varieties of bread wheat, 16 durum, 20 triticale and 5 barley varieties, as well as common plant model species including tobacco and Arabidopsis thaliana. A varying degree of effector sensitivities was observed, with a higher ToxB sensitivity and prevalence in the durum and triticale varieties. ToxB-induced chlorosis was also detected on barley. The heterologous expression of effectors that are easily scalable, will facilitate effector-assisted selection of varieties in wheat breeding programs as well as the investigation of P. tritici-repentis effectors in host and non-host interactions.

A positive-charged patch and stabilized hydrophobic core are essential for avirulence function of AvrPib in the rice blast fungus

Fungal avirulence effectors, a key weapon utilized by pathogens to promote their infection, are recognized by immune receptors to boost host R gene-mediated resistance. Many avirulence effectors share sparse sequence homology to proteins with known functions, and their molecular and biochemical functions together with the evolutionary relationship among different members remain largely unknown. Here, the crystal structure of AvrPib, an avirulence effector from Magnaporthe oryzae, was determined and showed a high degree of similarity to the M. oryzae Avrs and ToxB (MAX) effectors. Compared with other MAX effectors, AvrPib has a distinct positive-charge patch formed by five positive-charged residues (K29, K30, R50, K52 and K70) on the surface. These five key residues were essential to avirulence function of AvrPib and affected its nuclear localization into host cells. Moreover, residues V39 and V58, which locate in the hydrophobic core of the structure, cause loss of function of AvrPib by single-point mutation in natural isolates. In comparison with the wild-type AvrPib, the V39A or V58A mutations resulted in a partial or entire loss of secondary structure elements. Taken together, our results suggest that differences in the surface charge distribution of avirulence proteins could be one of the major bases for the variation in effector-receptor specificity, and that destabilization of the hydrophobic core is one of the major mechanisms employed by AvrPib for the fungus to evade recognition by resistance factors in the host cell.

Lessons in Effector and NLR Biology of Plant-Microbe Systems

A diversity of plant-associated organisms secrete effectors-proteins and metabolites that modulate plant physiology to favor host infection and colonization. However, effectors can also activate plant immune receptors, notably nucleotide-binding domain and leucine-rich repeat region (NLR)-containing proteins, enabling plants to fight off invading organisms. This interplay between effectors, their host targets, and the matching immune receptors is shaped by intricate molecular mechanisms and exceptionally dynamic coevolution. In this article, we focus on three effectors, AVR-Pik, AVR-Pia, and AVR-Pii, from the rice blast fungus Magnaporthe oryzae (syn. Pyricularia oryzae), and their corresponding rice NLR immune receptors, Pik, Pia, and Pii, to highlight general concepts of plant-microbe interactions. We draw 12 lessons in effector and NLR biology that have emerged from studying these three little effectors and are broadly applicable to other plant-microbe systems.

Lessons in Effector and NLR Biology of Plant-Microbe Systems

A diversity of plant-associated organisms secrete effectors-proteins and metabolites that modulate plant physiology to favor host infection and colonization. However, effectors can also activate plant immune receptors, notably nucleotide-binding domain and leucine-rich repeat region (NLR)-containing proteins, enabling plants to fight off invading organisms. This interplay between effectors, their host targets, and the matching immune receptors is shaped by intricate molecular mechanisms and exceptionally dynamic coevolution. In this article, we focus on three effectors, AVR-Pik, AVR-Pia, and AVR-Pii, from the rice blast fungus Magnaporthe oryzae (syn. Pyricularia oryzae), and their corresponding rice NLR immune receptors, Pik, Pia, and Pii, to highlight general concepts of plant-microbe interactions. We draw 12 lessons in effector and NLR biology that have emerged from studying these three little effectors and are broadly applicable to other plant-microbe systems.

Isolation and Characterization of Avirulence Genes in Magnaporthe oryzae

Magnaporthe oryzae is a fungal pathogen contributing to rice blast diseases globally via their Avr (avirulence) gene. Although the occurrence of M. oryzae has been reported in Sarawak since several decades ago, however, none has focused specifically on Avr genes, which confer resistance against pathogen associated molecular pattern-triggered immunity (PTI) in host. The objective of this study is to isolate Avr genes from M. oryzae 7’ (a Sarawak isolate) that may contribute to susceptibility of rice towards diseases. In this study, AvrPiz-t, AVR-Pik, Avr-Pi54, and AVR-Pita1 genes were isolated via PCR and cloning approaches. The genes were then compared with set of similar genes from related isolates derived from NCBI. Results revealed that all eight Avr genes (including four other global isolates) shared similar N-myristoylation site and a novel motif. 3D modeling revealed similar β-sandwich structure in AvrPiz-t and AVR-Pik despite sequence dissimilarities. In conclusion, it is confirmed of the presence of these genes in the Sarawak (M. oryzae) isolate. This study implies that Sarawak isolate may confer similar avirulence properties as their counterparts worldwide. Further R/Avr gene-for-gene relationship studies may aid in strategic control of rice blast diseases in future.

Inverse gene-for-gene interactions contribute additively to tan spot susceptibility in wheat

Tan spot susceptibility is conferred by multiple interactions of necrotrophic effector and host sensitivity genes. Tan spot of wheat, caused by Pyrenophora tritici-repentis, is an important disease in almost all wheat-growing areas of the world. The disease system is known to involve at least three fungal-produced necrotrophic effectors (NEs) that interact with the corresponding host sensitivity (S) genes in an inverse gene-for-gene manner to induce disease. However, it is unknown if the effects of these NE-S gene interactions contribute additively to the development of tan spot. In this work, we conducted disease evaluations using different races and quantitative trait loci (QTL) analysis in a wheat recombinant inbred line (RIL) population derived from a cross between two susceptible genotypes, LMPG-6 and PI 626573. The two parental lines each harbored a single known NE sensitivity gene with LMPG-6 having the Ptr ToxC sensitivity gene Tsc1 and PI 626573 having the Ptr ToxA sensitivity gene Tsn1. Transgressive segregation was observed in the population for all races. QTL mapping revealed that both loci (Tsn1 and Tsc1) were significantly associated with susceptibility to race 1 isolates, which produce both Ptr ToxA and Ptr ToxC, and the two genes contributed additively to tan spot susceptibility. For isolates of races 2 and 3, which produce only Ptr ToxA and Ptr ToxC, only Tsn1 and Tsc1 were associated with tan spot susceptibility, respectively. This work clearly demonstrates that tan spot susceptibility in this population is due primarily to two NE-S interactions. Breeders should remove both sensitivity genes from wheat lines to obtain high levels of tan spot resistance.

Effectors of Filamentous Plant Pathogens: Commonalities amid Diversity

Fungi and oomycetes are filamentous microorganisms that include a diversity of highly developed pathogens of plants. These are sophisticated modulators of plant processes that secrete an arsenal of effector proteins to target multiple host cell compartments and enable parasitic infection. Genome sequencing revealed complex catalogues of effectors of filamentous pathogens, with some species harboring hundreds of effector genes. Although a large fraction of these effector genes encode secreted proteins with weak or no sequence similarity to known proteins, structural studies have revealed unexpected similarities amid the diversity. This article reviews progress in our understanding of effector structure and function in light of these new insights. We conclude that there is emerging evidence for multiple pathways of evolution of effectors of filamentous plant pathogens but that some families have probably expanded from a common ancestor by duplication and diversification. Conserved folds, such as the oomycete WY and the fungal MAX domains, are not predictive of the precise function of the effectors but serve as a chassis to support protein structural integrity while providing enough plasticity for the effectors to bind different host proteins and evolve unrelated activities inside host cells. Further effector evolution and diversification arise via short linear motifs, domain integration and duplications, and oligomerization.

Recognition of the Magnaporthe oryzae Effector AVR-Pia by the Decoy Domain of the Rice NLR Immune Receptor RGA5

Nucleotide binding domain and leucine-rich repeat proteins (NLRs) are important receptors in plant immunity that allow recognition of pathogen effectors. The rice (Oryza sativa) NLR RGA5 recognizes the Magnaporthe oryzae effector AVR-Pia through direct interaction. Here, we gained detailed insights into the molecular and structural bases of AVR-Pia-RGA5 interaction and the role of the RATX1 decoy domain of RGA5. NMR titration combined with in vitro and in vivo protein-protein interaction analyses identified the AVR-Pia interaction surface that binds to the RATX1 domain. Structure-informed AVR-Pia mutants showed that, although AVR-Pia associates with additional sites in RGA5, binding to the RATX1 domain is necessary for pathogen recognition but can be of moderate affinity. Therefore, RGA5-mediated resistance is highly resilient to mutations in the effector. We propose a model that explains such robust effector recognition as a consequence, and an advantage, of the combination of integrated decoy domains with additional independent effector-NLR interactions.

Efficient production of (2)H, (13)C, (15)N-enriched industrial enzyme Rhizopus chinensis lipase with native disulfide bonds

BACKGROUND

In order to use most modern methods of NMR spectroscopy to study protein structure and dynamics, isotope-enriched protein samples are essential. Especially for larger proteins (>20 kDa), perdeuterated and Ile (δ1), Leu, and Val methyl-protonated protein samples are required for suppressing nuclear relaxation to provide improved spectral quality, allowing key backbone and side chain resonance assignments needed for protein structure and dynamics studies. Escherichia coli and Pichia pastoris are two of the most popular expression systems for producing isotope-enriched, recombinant protein samples for NMR investigations. The P. pastoris system can be used to produce (13)C, (15)N-enriched and even (2)H,(13)C, (15)N-enriched protein samples, but efficient methods for producing perdeuterated proteins with Ile (δ1), Leu and Val methyl-protonated groups in P. pastoris are still unavailable. Glycosylation heterogeneity also provides challenges to NMR studies. E. coli expression systems are efficient for overexpressing perdeuterated and Ile (δ1), Leu, Val methyl-protonated protein samples, but are generally not successful for producing secreted eukaryotic proteins with native disulfide bonds.

RESULTS

The 33 kDa protein-Rhizopus chinensis lipase (RCL), an important industrial enzyme, was produced using both P. pastoris and E. coli BL21 trxB (DE3) systems. Samples produced from both systems exhibit identical native disulfide bond formation and similar 2D NMR spectra, indicating similar native protein folding. The yield of (13)C, (15)N-enriched r27RCL produced using P. pastoris was 1.7 times higher that obtained using E. coli, while the isotope-labeling efficiency was ~15 % lower. Protein samples produced in P. pastoris exhibit O-glycosylation, while the protein samples produced in E. coli were not glycosylated. The specific activity of r27RCL from P. pastoris was ~1.4 times higher than that produced in E. coli.

CONCLUSIONS

These data demonstrate efficient production of (2)H, (13)C, (15)N-enriched, Ile (δ1), Leu, Val methyl-protonated eukaryotic protein r27RCL with native disulfides using the E. coli BL21 trxB (DE3) system. For certain NMR studies, particularly efforts for resonance assignments, structural studies, and dynamic studies, E. coli provides a cost-effective system for producing isotope-enriched RCL. It should also be potential for producing other (2)H, (13)C, (15)N-enriched, Ile (δ1), Leu, Val methyl-protonated eukaryotic proteins with native disulfide bonds.

Structure Analysis Uncovers a Highly Diverse but Structurally Conserved Effector Family in Phytopathogenic Fungi

Phytopathogenic ascomycete fungi possess huge effector repertoires that are dominated by hundreds of sequence-unrelated small secreted proteins. The molecular function of these effectors and the evolutionary mechanisms that generate this tremendous number of singleton genes are largely unknown. To get a deeper understanding of fungal effectors, we determined by NMR spectroscopy the 3-dimensional structures of the Magnaporthe oryzae effectors AVR1-CO39 and AVR-Pia. Despite a lack of sequence similarity, both proteins have very similar 6 β-sandwich structures that are stabilized in both cases by a disulfide bridge between 2 conserved cysteins located in similar positions of the proteins. Structural similarity searches revealed that AvrPiz-t, another effector from M. oryzae, and ToxB, an effector of the wheat tan spot pathogen Pyrenophora tritici-repentis have the same structures suggesting the existence of a family of sequence-unrelated but structurally conserved fungal effectors that we named MAX-effectors (MagnaportheAvrs and ToxB like). Structure-informed pattern searches strengthened this hypothesis by identifying MAX-effector candidates in a broad range of ascomycete phytopathogens. Strong expansion of the MAX-effector family was detected in M. oryzae and M. grisea where they seem to be particularly important since they account for 5–10% of the effector repertoire and 50% of the cloned avirulence effectors. Expression analysis indicated that the majority of M. oryzae MAX-effectors are expressed specifically during early infection suggesting important functions during biotrophic host colonization. We hypothesize that the scenario observed for MAX-effectors can serve as a paradigm for ascomycete effector diversity and that the enormous number of sequence-unrelated ascomycete effectors may in fact belong to a restricted set of structurally conserved effector families.

Fungal plant pathogens are of outstanding economic and ecological importance and cause destructive diseases on many cultivated and wild plants. Effector proteins that are secreted during infection to manipulate the host and to promote disease are a key element in fungal virulence. Phytopathogenic fungi possess huge effector repertoires that are dominated by hundreds of sequence-unrelated small secreted proteins. The molecular functions of this most important class of fungal effectors and the evolutionary mechanisms that generate this tremendous numbers of apparently unrelated proteins are largely unknown. By investigating the 3-dimensional structures of effectors from the rice blast fungus M. oryzae, we discovered an effector family comprising structurally conserved but sequence-unrelated effectors from M. oryzae and the phylogenetically distant wheat pathogen Pyrenophora tritici-repentis that we named MAX-effectors (M. oryzaeAvrs and ToxB). Structure-informed searches of whole genome sequence databases suggest that MAX-effectors are present at low frequencies and with a patchy phylogenetic distribution in many ascomycete phytopathogens. They underwent strong lineage-specific expansion in fungi of the Pyriculariae family that contains M. oryzae where they seem particularly important during biotrophic plant colonization and account for 50% of the cloned Avr effectors and 5–10% of the effector repertoire. Based on our results on the MAX-effectors and the widely accepted concept that fungal effectors evolve according to a birth-and-death model we propose the hypothesis that the majority of the immense numbers of different ascomycete effectors could in fact belong to a limited set of structurally defined families whose members are phylogenetically related.

Structural basis of pathogen recognition by an integrated HMA domain in a plant NLR immune receptor

Plants have evolved intracellular immune receptors to detect pathogen proteins known as effectors. How these immune receptors detect effectors remains poorly understood. Here we describe the structural basis for direct recognition of AVR-Pik, an effector from the rice blast pathogen, by the rice intracellular NLR immune receptor Pik. AVR-PikD binds a dimer of the Pikp-1 HMA integrated domain with nanomolar affinity. The crystal structure of the Pikp-HMA/AVR-PikD complex enabled design of mutations to alter protein interaction in yeast and in vitro, and perturb effector-mediated response both in a rice cultivar containing Pikp and upon expression of AVR-PikD and Pikp in the model plant Nicotiana benthamiana. These data reveal the molecular details of a recognition event, mediated by a novel integrated domain in an NLR, which initiates a plant immune response and resistance to rice blast disease. Such studies underpin novel opportunities for engineering disease resistance to plant pathogens in staple food crops.

Residue-level global and local ensemble-ensemble comparisons of protein domains

Many methods of protein structure generation such as NMR-based solution structure determination and template-based modeling do not produce a single model, but an ensemble of models consistent with the available information. Current strategies for comparing ensembles lose information because they use only a single representative structure. Here, we describe the ENSEMBLATOR and its novel strategy to directly compare two ensembles containing the same atoms to identify significant global and local backbone differences between them on per-atom and per-residue levels, respectively. The ENSEMBLATOR has four components: eePREP (ee for ensemble-ensemble), which selects atoms common to all models; eeCORE, which identifies atoms belonging to a cutoff-distance dependent common core; eeGLOBAL, which globally superimposes all models using the defined core atoms and calculates for each atom the two intraensemble variations, the interensemble variation, and the closest approach of members of the two ensembles; and eeLOCAL, which performs a local overlay of each dipeptide and, using a novel measure of local backbone similarity, reports the same four variations as eeGLOBAL. The combination of eeGLOBAL and eeLOCAL analyses identifies the most significant differences between ensembles. We illustrate the ENSEMBLATOR's capabilities by showing how using it to analyze NMR ensembles and to compare NMR ensembles with crystal structures provides novel insights compared to published studies. One of these studies leads us to suggest that a "consistency check" of NMR-derived ensembles may be a useful analysis step for NMR-based structure determinations in general. The ENSEMBLATOR 1.0 is available as a first generation tool to carry out ensemble-ensemble comparisons.