Functional and computational analysis of amino acid patterns predictive of type III secretion system substrates in Pseudomonas syringae

Features and algorithms: facilitating investigation of secreted effectors in Gram-negative bacteria

Gram-negative bacteria deliver effector proteins through type III, IV, or VI secretion systems (T3SSs, T4SSs, and T6SSs) into host cells, causing infections and diseases. In general, effector proteins for each of these distinct secretion systems lack homology and are difficult to identify. Sequence analysis has disclosed many common features, helping us to understand the evolution, function, and secretion mechanisms of the effectors. In combination with various algorithms, the known common features have facilitated accurate prediction of new effectors. Ensemblers or integrated pipelines achieve a better prediction of performance, which combines multiple computational models or modules with multidimensional features. Natural language processing (NLP) models also show the merits, which could enable discovery of novel features and, in turn, facilitate more precise effector prediction, extending our knowledge about each secretion mechanism.

Effectors of Puccinia striiformis f. sp. tritici Suppressing the Pathogenic-Associated Molecular Pattern-Triggered Immune Response Were Screened by Transient Expression of Wheat Protoplasts

Puccinia striiformis f. sp. tritici (Pst) is an important pathogen of wheat (Triticum aestivum L.) stripe rust, and the effector protein secreted by haustoria is a very important component involved in the pathogenic process. Although the candidate effector proteins secreted by Pst haustoria have been predicted to be abundant, few have been functionally validated. Our study confirmed that chitin and flg22 could be used as elicitors of the pathogenic-associated molecular pattern-triggered immune (PTI) reaction in wheat leaves and that TaPr-1-14 could be used as a marker gene to detect the PTI reaction. In addition, the experimental results were consistent in wheat protoplasts. A rapid and efficient method for screening and identifying the effector proteins of Pst was established by using the wheat protoplast transient expression system. Thirty-nine Pst haustorial effector genes were successfully cloned and screened for expression in the protoplast. We identified three haustorial effector proteins, PSEC2, PSEC17, and PSEC45, that may inhibit the response of wheat to PTI. These proteins are localized in the somatic cytoplasm and nucleus of wheat protoplasts and are highly expressed during the infection and parasitism of wheat.

Computational prediction of secreted proteins in gram-negative bacteria

Gram-negative bacteria harness multiple protein secretion systems and secrete a large proportion of the proteome. Proteins can be exported to periplasmic space, integrated into membrane, transported into extracellular milieu, or translocated into cytoplasm of contacting cells. It is important for accurate, genome-wide annotation of the secreted proteins and their secretion pathways. In this review, we systematically classified the secreted proteins according to the types of secretion systems in Gram-negative bacteria, summarized the known features of these proteins, and reviewed the algorithms and tools for their prediction.

Deciphering Molecular Virulence Mechanism of Mycobacterium tuberculosis Dop isopeptidase Based on Its Sequence-Structure-Function Linkage

The pupylation pathway marks proteins for prokaryotic ubiquitin-like protein (Pup)-proteasomal degradation and survival strategy of mycobacteria inside of the host macrophages. Deamidase of Pup (Dop) plays a central role in the pupylation pathway. It is still a matter of investigation to know the function of Dop in virulence of mycobacterial lineage. Hence, the present study was intended to describe the sequence-structure-function-virulence link of Dop for understanding the molecular virulence mechanism of Mycobacterium tuberculosis H37Rv (Mtb). Phylogenetic analysis of this study indicated that Dop has extensively diverged across the proteasome-harboring bacteria. The functional part of Dop was converged across the pathogenic mycobacterial lineage. The genome-wide analysis pointed out that the pupylation gene locus was identical to each other, but its genome neighborhood differed from species to species. Molecular modeling and dynamic studies proved that the predicted structure of Mtb Dop was energetically stable and low conformational freedom. Moreover, evolutionary constraints in Mtb Dop were intensively analyzed for inferring its sequence-structure-function relationships for the full virulence of Mtb. It indicated that evolutionary optimization was extensively required to stabilize its local structural environment at the side chains of mutable residues. The sequence-structure-function-virulence link of Dop might have retained in Mtb by reordering hydrophobic and hydrogen bonding patterns in the local structural environment. Thus, the results of our study provide a quest to understand the molecular virulence and pathogenesis mechanisms of Mtb during the infection process.

Phylogenetic profiling, an untapped resource for the prediction of secreted proteins and its complementation with sequence-based classifiers in bacterial type III, IV and VI secretion systems

In the establishment and maintenance of the interaction between pathogenic or symbiotic bacteria with a eukaryotic organism, protein substrates of specialized bacterial secretion systems called effectors play a critical role once translocated into the host cell. Proteins are also secreted to the extracellular medium by free-living bacteria or directly injected into other competing organisms to hinder or kill. In this work, we explore an approach based on the evolutionary dependence that most of the effectors maintain with their specific secretion system that analyzes the co-occurrence of any orthologous protein group and their corresponding secretion system across multiple genomes. We compared and complemented our methodology with sequence-based machine learning prediction tools for the type III, IV and VI secretion systems. Finally, we provide the predictive results for the three secretion systems in 1606 complete genomes at http://www.iib.unsam.edu.ar/orgsissec/.

Effector Translocation: Cya Reporter Assay

An accurate and complete roster of the Type III effector (T3E) proteins translocated by the P. syringae Type III secretion system (T3SS) into host cells is critical to understanding the pathogen's interactions with plants. The adenylate cyclase (Cya) reporter offers a highly sensitive and robust assay for monitoring the translocation of T3Es. T3Es are fused to the calmodulin-dependent adenylate-cyclase domain of CyaA. The T3E targets Cya for translocation through the T3SS into the host cell at which point it is activated by calmodulin and converts adenosine triphosphate into cyclic adenosine monophosphate (cAMP). The T3SS translocation-dependent increase in cAMP concentration in plant cells is then measured with an enzyme-linked immunosorbent assay kit. The Cya reporter can be used to determine whether a candidate protein is translocated by T3SS or to measure relative levels of T3SS translocation in a semiquantitative manner.

Behind the lines-actions of bacterial type III effector proteins in plant cells

Pathogenicity of most Gram-negative plant-pathogenic bacteria depends on the type III secretion (T3S) system, which translocates bacterial effector proteins into plant cells. Type III effectors modulate plant cellular pathways to the benefit of the pathogen and promote bacterial multiplication. One major virulence function of type III effectors is the suppression of plant innate immunity, which is triggered upon recognition of pathogen-derived molecular patterns by plant receptor proteins. Type III effectors also interfere with additional plant cellular processes including proteasome-dependent protein degradation, phytohormone signaling, the formation of the cytoskeleton, vesicle transport and gene expression. This review summarizes our current knowledge on the molecular functions of type III effector proteins with known plant target molecules. Furthermore, plant defense strategies for the detection of effector protein activities or effector-triggered alterations in plant targets are discussed.

The Type III Secretion Chaperone HpaB Controls the Translocation of Effector and Noneffector Proteins From Xanthomonas campestris pv. vesicatoria

Pathogenicity of the gram-negative bacterium Xanthomonas campestris pv. vesicatoria depends on a type III secretion (T3S) system, which translocates effector proteins into plant cells. Effector proteins contain N-terminal T3S and translocation signals and interact with the T3S chaperone HpaB, which presumably escorts effectors to the secretion apparatus. The molecular mechanisms underlying the recognition of effectors by the T3S system are not yet understood. In the present study, we analyzed T3S and translocation signals in the type III effectors XopE2 and XopJ from X. campestris pv. vesicatoria. Both effectors contain minimal translocation signals, which are only recognized in the absence of HpaB. Additional N-terminal signals promote translocation of XopE2 and XopJ in the wild-type strain. The results of translocation and interaction studies revealed that the interaction of XopE2 and XopJ with HpaB and a predicted cytoplasmic substrate docking site of the T3S system is not sufficient for translocation. In agreement with this finding, we show that the presence of an artificial HpaB-binding site does not promote translocation of the noneffector XopA in the wild-type strain. Our data, therefore, suggest that the T3S chaperone HpaB not only acts as an escort protein but also controls the recognition of translocation signals.

The TAL Effector AvrBs3 from Xanthomonas campestris pv. vesicatoria Contains Multiple Export Signals and Can Enter Plant Cells in the Absence of the Type III Secretion Translocon

Pathogenicity of the Gram-negative plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria depends on a type III secretion (T3S) system which translocates effector proteins into plant cells. Effector protein delivery is controlled by the T3S chaperone HpaB, which presumably escorts effector proteins to the secretion apparatus. One intensively studied effector is the transcription activator-like (TAL) effector AvrBs3, which binds to promoter sequences of plant target genes and activates plant gene expression. It was previously reported that type III-dependent delivery of AvrBs3 depends on the N-terminal protein region. The signals that control T3S and translocation of AvrBs3, however, have not yet been characterized. In the present study, we show that T3S and translocation of AvrBs3 depend on the N-terminal 10 and 50 amino acids, respectively. Furthermore, we provide experimental evidence that additional signals in the N-terminal 30 amino acids and the region between amino acids 64 and 152 promote translocation of AvrBs3 in the absence of HpaB. Unexpectedly, in vivo translocation assays revealed that AvrBs3 is delivered into plant cells even in the absence of HrpF, which is the predicted channel-forming component of the T3S translocon in the plant plasma membrane. The presence of HpaB- and HrpF-independent transport routes suggests that the delivery of AvrBs3 is initiated during early stages of the infection process, presumably before the activation of HpaB or the insertion of the translocon into the plant plasma membrane.

Computational approach to predict species-specific type III secretion system (T3SS) effectors using single and multiple genomes

BACKGROUND

Many gram-negative bacteria use type III secretion systems (T3SSs) to translocate effector proteins into host cells. T3SS effectors can give some bacteria a competitive edge over others within the same environment and can help bacteria to invade the host cells and allow them to multiply rapidly within the host. Therefore, developing efficient methods to identify effectors scattered in bacterial genomes can lead to a better understanding of host-pathogen interactions and ultimately to important medical and biotechnological applications.

RESULTS

We used 21 genomic and proteomic attributes to create a precise and reliable T3SS effector prediction method called Genome Search for Effectors Tool (GenSET). Five machine learning algorithms were trained on effectors selected from different organisms and a trained (voting) algorithm was then applied to identify other effectors present in the genome testing sets from the same (GenSET Phase 1) or different (GenSET Phase 2) organism. Although a select group of attributes that included the codon adaptation index, probability of expression in inclusion bodies, N-terminal disorder, and G + C content (filtered) were better at discriminating between positive and negative sets, algorithm performance was better when all 21 attributes (unfiltered) were used. Performance scores (sensitivity, specificity and area under the curve) from GenSET Phase 1 were better than those reported for six published methods. More importantly, GenSET Phase 1 ranked more known effectors (70.3%) in the top 40 ranked proteins and predicted 10-80% more effectors than three available programs in three of the four organisms tested. GenSET Phase 2 predicted 43.8% effectors in the top 40 ranked proteins when tested on four related or unrelated organisms. The lower prediction rates from GenSET Phase 2 may be due to the presence of different translocation signals in effectors from different T3SS families.

CONCLUSIONS

The species-specific GenSET Phase 1 method offers an alternative approach to T3SS effector prediction that can be used with other published programs to improve effector predictions. Additionally, our approach can be applied to predict effectors of other secretion systems as long as these effectors have translocation signals embedded in their sequences.

Type III-Dependent Translocation of HrpB2 by a Nonpathogenic hpaABC Mutant of the Plant-Pathogenic Bacterium Xanthomonas campestris pv. vesicatoria

The plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria employs a type III secretion (T3S) system to translocate effector proteins into plant cells. The T3S apparatus spans both bacterial membranes and is associated with an extracellular pilus and a channel-like translocon in the host plasma membrane. T3S is controlled by the switch protein HpaC, which suppresses secretion and translocation of the predicted inner rod protein HrpB2 and promotes secretion of translocon and effector proteins. We previously reported that HrpB2 interacts with HpaC and the cytoplasmic domain of the inner membrane protein HrcU (C. Lorenz, S. Schulz, T. Wolsch, O. Rossier, U. Bonas, and D. Büttner, PLoS Pathog 4:e1000094, 2008, http://dx.doi.org/10.1371/journal.ppat.1000094). However, the molecular mechanisms underlying the control of HrpB2 secretion are not yet understood. Here, we located a T3S and translocation signal in the N-terminal 40 amino acids of HrpB2. The results of complementation experiments with HrpB2 deletion derivatives revealed that the T3S signal of HrpB2 is essential for protein function. Furthermore, interaction studies showed that the N-terminal region of HrpB2 interacts with the cytoplasmic domain of HrcU, suggesting that the T3S signal of HrpB2 contributes to substrate docking. Translocation of HrpB2 is suppressed not only by HpaC but also by the T3S chaperone HpaB and its secreted regulator, HpaA. Deletion of hpaA, hpaB, and hpaC leads to a loss of pathogenicity but allows the translocation of fusion proteins between the HrpB2 T3S signal and effector proteins into leaves of host and non-host plants.

IMPORTANCE

The T3S system of the plant-pathogenic bacterium Xanthomonas campestris pv. vesicatoria is essential for pathogenicity and delivers effector proteins into plant cells. T3S depends on HrpB2, which is a component of the predicted periplasmic inner rod structure of the secretion apparatus. HrpB2 is secreted during the early stages of the secretion process and interacts with the cytoplasmic domain of the inner membrane protein HrcU. Here, we localized the secretion and translocation signal of HrpB2 in the N-terminal 40 amino acids and show that this region is sufficient for the interaction with the cytoplasmic domain of HrcU. Our results suggest that the T3S signal of HrpB2 is required for the docking of HrpB2 to the secretion apparatus. Furthermore, we provide experimental evidence that the N-terminal region of HrpB2 is sufficient to target effector proteins for translocation in a nonpathogenic X. campestris pv. vesicatoria strain.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.

Global analysis of the HrpL regulon in the plant pathogen Pseudomonas syringae pv. tomato DC3000 reveals new regulon members with diverse functions

The type III secretion system (T3SS) is required for virulence in the gram-negative plant pathogen Pseudomonas syringae pv. tomato DC3000. The alternative sigma factor HrpL directly regulates expression of T3SS genes via a promoter sequence, often designated as the "hrp promoter." Although the HrpL regulon has been extensively investigated in DC3000, it is not known whether additional regulon members remain to be found. To systematically search for HrpL-regulated genes, we used chromatin immunoprecipitation coupled with high-throughput sequencing (ChIP-Seq) and bulk mRNA sequencing (RNA-Seq) to identify HrpL-binding sites and likely hrp promoters. The analysis recovered 73 sites of interest, including 20 sites that represent new hrp promoters. The new promoters lie upstream of a diverse set of genes encoding potential regulators, enzymes and hypothetical proteins. PSPTO_5633 is the only new HrpL regulon member that is potentially an effector and is now designated HopBM1. Deletions in several other new regulon members, including PSPTO_5633, PSPTO_0371, PSPTO_2130, PSPTO_2691, PSPTO_2696, PSPTO_3331, and PSPTO_5240, in either DC3000 or ΔhopQ1-1 backgrounds, do not affect the hypersensitive response or in planta growth of the resulting strains. Many new HrpL regulon members appear to be unrelated to the T3SS, and orthologs for some of these can be identified in numerous non-pathogenic bacteria. With the identification of 20 new hrp promoters, the list of HrpL regulon members is approaching saturation and most likely includes all DC3000 effectors.

Consequences of flagellin export through the type III secretion system of Pseudomonas syringae reveal a major difference in the innate immune systems of mammals and the model plant Nicotiana benthamiana

Bacterial flagellin is perceived as a microbe (or pathogen)-associated molecular pattern (MAMP or PAMP) by the extracellular pattern recognition receptors, FLS2 and TLR5, of plants and mammals respectively. Flagellin accidently translocated into mammalian cells by pathogen type III secretion systems (T3SSs) is recognized by nucleotide-binding leucine-rich repeat receptor NLRC4 as a pattern of pathogenesis and induces a death-associated immune response. The non-pathogen Pseudomonas fluorescens Pf0-1, expressing a Pseudomonas syringae T3SS, and the plant pathogen P. syringae pv. tomato DC3000 were used to seek evidence of an analogous cytoplasmic recognition system for flagellin in the model plant Nicotiana benthamiana. Flagellin (FliC) was secreted in culture and translocated into plant cells by the T3SS expressed in Pf0-1 and DC3000 and in their ΔflgGHI flagellar pathway mutants. ΔfliC and ΔflgGHI mutants of Pf0-1 and DC3000 were strongly reduced in elicitation of reactive oxygen species production and in immunity induction as indicated by the ability of challenge bacteria inoculated 6 h later to translocate a type III effector-reporter and to elicit effector-triggered cell death. Agrobacterium-mediated transient expression in N. benthamiana of FliC with or without a eukaryotic export signal peptide, coupled with virus-induced gene silencing of FLS2, revealed no immune response that was not FLS2 dependent. Transiently expressed FliC from DC3000 and Pectobacterium carotovorum did notinduce cell death in N. benthamiana, tobacco or tomato leaves. Flagellin is the major Pseudomonas MAMP perceived by N. benthamiana, and although flagellin secretion through the plant cell wall by the T3SS may partially contribute to FLS2-dependent immunity, flagellin in the cytosol does not elicit immune-associated cell death. We postulate that a death response to translocated MAMPs would produce vulnerability to the many necrotrophic pathogens of plants, such as P. carotovorum, which differ from P. syringae and other (hemi)biotrophic pathogens in benefitting from death-associated immune responses.

Using weakly conserved motifs hidden in secretion signals to identify type-III effectors from bacterial pathogen genomes

Background

As one of the most important virulence factor types in gram-negative pathogenic bacteria, type-III effectors (TTEs) play a crucial role in pathogen-host interactions by directly influencing immune signaling pathways within host cells. Based on the hypothesis that type-III secretion signals may be comprised of some weakly conserved sequence motifs, here we used profile-based amino acid pair information to develop an accurate TTE predictor.

Results

For a TTE or non-TTE, we first used a hidden Markov model-based sequence searching method (i.e., HHblits) to detect its weakly homologous sequences and extracted the profile-based k-spaced amino acid pair composition (HH-CKSAAP) from the N-terminal sequences. In the next step, the feature vector HH-CKSAAP was used to train a linear support vector machine model, which we designate as BEAN (Bacterial Effector ANalyzer). We compared our method with four existing TTE predictors through an independent test set, and our method revealed improved performance. Furthermore, we listed the most predictive amino acid pairs according to their weights in the established classification model. Evolutionary analysis shows that predictive amino acid pairs tend to be more conserved. Some predictive amino acid pairs also show significantly different position distributions between TTEs and non-TTEs. These analyses confirmed that some weakly conserved sequence motifs may play important roles in type-III secretion signals. Finally, we also used BEAN to scan one plant pathogen genome and showed that BEAN can be used for genome-wide TTE identification. The webserver and stand-alone version of BEAN are available at http://protein.cau.edu.cn:8080/bean/.