Terminal reassortment drives the quantum evolution of type III effectors in bacterial pathogens

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.

Effector Identification in Plant Pathogens

Effectors play a central role in determining the outcome of plant-pathogen interactions. As key virulence proteins, effectors are collectively indispensable for disease development. By understanding the virulence mechanisms of effectors, fundamental knowledge of microbial pathogenesis and disease resistance have been revealed. Effectors are also considered double-edged swords because some of them activate immunity in disease resistant plants after being recognized by specific immune receptors, which evolved to monitor pathogen presence or activity. Characterization of effector recognition by their cognate immune receptors and the downstream immune signaling pathways is instrumental in implementing resistance. Over the past decades, substantial research effort has focused on effector biology, especially concerning their interactions with virulence targets or immune receptors in plant cells. A foundation of this research is robust identification of the effector repertoire from a given pathogen, which depends heavily on bioinformatic prediction. In this review, we summarize methodologies that have been used for effector mining in various microbial pathogens which use different effector delivery mechanisms. We also discuss current limitations and provide perspectives on how recently developed analytic tools and technologies may facilitate effector identification and hence generation of a more complete vision of host-pathogen interactions. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Pangenomic analysis reveals plant NAD+ manipulation as an important virulence activity of bacterial pathogen effectors

Nicotinamide adenine dinucleotide (NAD⁺) has emerged as a key component in prokaryotic and eukaryotic immune systems. The recent discovery that Toll/interleukin-1 receptor (TIR) proteins function as NAD⁺ hydrolases (NADase) links NAD⁺-derived small molecules with immune signaling. We investigated pathogen manipulation of host NAD⁺ metabolism as a virulence strategy. Using the pangenome of the model bacterial pathogen Pseudomonas syringae, we conducted a structure-based similarity search from 35,000 orthogroups for type III effectors (T3Es) with potential NADase activity. Thirteen T3Es, including five newly identified candidates, were identified that possess domain(s) characteristic of seven NAD⁺-hydrolyzing enzyme families. Most Pseudomonas syringae strains that depend on the type III secretion system to cause disease, encode at least one NAD⁺-manipulating T3E, and many have several. We experimentally confirmed the type III-dependent secretion of a novel T3E, named HopBY, which shows structural similarity to both TIR and adenosine diphosphate ribose (ADPR) cyclase. Homologs of HopBY were predicted to be type VI effectors in diverse bacterial species, indicating potential recruitment of this activity by microbial proteins secreted during various interspecies interactions. HopBY efficiently hydrolyzes NAD⁺ and specifically produces 2'cADPR, which can also be produced by TIR immune receptors of plants and by other bacteria. Intriguingly, this effector promoted bacterial virulence, indicating that 2'cADPR may not be the signaling molecule that directly initiates immunity. This study highlights a host-pathogen battleground centered around NAD⁺ metabolism and provides insight into the NAD⁺-derived molecules involved in plant immunity.

Cooperative virulence via the collective action of secreted pathogen effectors

Although virulence is typically attributed to single pathogenic strains, here we investigated whether effectors secreted by a population of non-virulent strains could function as public goods to enable the emergence of collective virulence. We disaggregated the 36 type III effectors of the phytopathogenic bacterium Pseudomonas syringae strain PtoDC3000 into a 'metaclone' of 36 coisogenic strains, each carrying a single effector in an effectorless background. Each coisogenic strain was individually unfit, but the metaclone was collectively as virulent as the wild-type strain on Arabidopsis thaliana, suggesting that effectors can drive the emergence of cooperation-based virulence through their public action. We show that independently evolved effector suits can equally drive this cooperative behaviour by transferring the effector alleles native to the strain PmaES4326 into the conspecific but divergent strain PtoDC3000. Finally, we transferred the disaggregated PtoDC3000 effector arsenal into Pseudomonas fluorescens and show that their cooperative action was sufficient to convert this rhizosphere-inhabiting beneficial bacterium into a phyllosphere pathogen. These results emphasize the importance of microbial community interactions and expand the ecological scale at which disease may be attributed.

Diversity, Evolution, and Function of Pseudomonas syringae Effectoromes

Pseudomonas syringae is an evolutionarily diverse bacterial species complex and a preeminent model for the study of plant-pathogen interactions due in part to its remarkably broad host range. A critical feature of P. syringae virulence is the employment of suites of type III secreted effector (T3SE) proteins, which vary widely in composition and function. These effectors act on a variety of plant intracellular targets to promote pathogenesis but can also be avirulence factors when detected by host immune complexes. In this review, we survey the phylogenetic diversity (PD) of the P. syringae effectorome, comprising 70 distinct T3SE families identified to date, and highlight how avoidance of host immune detection has shaped effectorome diversity through functional redundancy, diversification, and horizontal transfer. We present emerging avenues for research and novel insights that can be gained via future investigations of plant-pathogen interactions through the fusion of large-scale interaction screens and phylogenomic approaches.

The NEL Family of Bacterial E3 Ubiquitin Ligases

Some pathogenic or symbiotic Gram-negative bacteria can manipulate the ubiquitination system of the eukaryotic host cell using a variety of strategies. Members of the genera Salmonella, Shigella, Sinorhizobium, and Ralstonia, among others, express E3 ubiquitin ligases that belong to the NEL family. These bacteria use type III secretion systems to translocate these proteins into host cells, where they will find their targets. In this review, we first introduce type III secretion systems and the ubiquitination process and consider the various ways bacteria use to alter the ubiquitin ligation machinery. We then focus on the members of the NEL family, their expression, translocation, and subcellular localization in the host cell, and we review what is known about the structure of these proteins, their function in virulence or symbiosis, and their specific targets.

The evolution and role of eukaryotic-like domains in environmental intracellular bacteria: the battle with a eukaryotic cell

Intracellular pathogens that are able to thrive in different environments, such as Legionella spp. which preferentially live in protozoa in aquatic environments or environmental Chlamydiae which replicate either within protozoa or a range of animals, possess a plethora of cellular biology tools to influence their eukaryotic host. The host manipulation tools that evolved in the interaction with protozoa, confer these bacteria the capacity to also infect phylogenetically distinct eukaryotic cells, such as macrophages and thus they can also be human pathogens. To manipulate the host cell, bacteria use protein secretion systems and molecular effectors. Although these molecular effectors are encoded in bacteria, they are expressed and function in a eukaryotic context often mimicking or inhibiting eukaryotic proteins. Indeed, many of these effectors have eukaryotic-like domains. In this review we propose that the main pathways environmental intracellular bacteria need to subvert in order to establish the host eukaryotic cell as a replication niche are chromatin remodelling, ubiquitination signalling, and modulation of protein-protein interactions via tandem repeat domains. We then provide mechanistic insight into how these proteins might have evolved as molecular weapons. Finally, we highlight that in environmental intracellular bacteria the number of eukaryotic-like domains and proteins is considerably higher than in intracellular bacteria specialised to an isolated niche, such as obligate intracellular human pathogens. As mimics of eukaryotic proteins are critical components of host pathogen interactions, this distribution of eukaryotic-like domains suggests that the environment has selected them.

The ETS-ETI cycle: evolutionary processes and metapopulation dynamics driving the diversification of pathogen effectors and host immune factors

The natural diversity of pathogen effectors and host immune components represents a snapshot of the underlying evolutionary processes driving the host-pathogen arms race. In plants, this arms race is manifested by an ongoing cycle of disease and resistance driven by pathogenic effectors that promote disease (effector-triggered susceptibility; ETS) and plant resistance proteins that recognize effector activity to trigger immunity (effector-triggered immunity; ETI). Here we discuss how this ongoing ETS-ETI cycle has shaped the natural diversity of both plant resistance proteins and pathogen effectors. We focus on the evolutionary forces that drive the diversification of the molecules that determine the outcome of plant-pathogen interactions and introduce the concept of metapopulation dynamics (i.e., the introduction of genetic variation from conspecific organisms in different populations) as an alternative mechanism that can introduce and maintain diversity in both host and pathogen populations.

How to win a tug-of-war: the adaptive evolution of Phytophthora effectors

The 'zigzag' model formulates some of the fundamental principles underpinning the dynamic interactions between pathogen effectors and plant immunity. As key virulence factors, effectors often exhibit a pattern of rapid evolution, presumably as a result of the host-pathogen arms race. Here, we summarize the current knowledge of mechanisms that may accelerate effector evolution in the highly successful Phytophthora pathogens. Recent findings on epigenetic regulation of effector genes that allows evasion of host recognition and maintenance of cost/benefit balance, and a conserved structural unit in effector proteins that may promote the evolution of virulence activities are highlighted.

The Conserved Arginine Required for AvrRps4 Processing Is Also Required for Recognition of Its N-Terminal Fragment in Lettuce

Pathogens utilize a repertoire of effectors to facilitate pathogenesis, but when the host recognizes one of them, it causes effector-triggered immunity. The Pseudomonas type III effector AvrRps4 is a bipartite effector that is processed in planta into a functional 133-amino acid N-terminus (AvrRps4-N) and 88-amino acid C-terminus (AvrRps4-C). Previous studies found AvrRps4-C to be sufficient to trigger the hypersensitive response (HR) in turnip. In contrast, our recent work found that AvrRps4-N but not AvrRps4-C triggered HR in lettuce, whereas both were required for resistance induction in Arabidopsis. Here, we initially compared AvrRps4 recognition by turnip and lettuce using transient expression. By serial truncation, we identified the central conserved region consisting of 37 amino acids as essential for AvrRps4-N recognition, whereas the putative type III secretion signal peptide or the C-terminal 13 amino acids were dispensable. Surprisingly, the conserved arginine at position 112 (R112) that is required for full-length AvrRps4 processing is also required for the recognition of AvrRps4-N by lettuce. Mutating R112 to hydrophobic leucine or negatively charged glutamate abolished the HR-inducing capacity of AvrRps4-N, while a positively charged lysine at this position resulted in a slow and weak HR. Together, our results suggest an AvrRps4-N recognition-specific role of R112 in lettuce.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Keeping in Touch with Type-III Secretion System Effectors: Mass Spectrometry-Based Proteomics to Study Effector-Host Protein-Protein Interactions

Manipulation of host cellular processes by translocated bacterial effectors is key to the success of bacterial pathogens and some symbionts. Therefore, a comprehensive understanding of effectors is of critical importance to understand infection biology. It has become increasingly clear that the identification of host protein targets contributes invaluable knowledge to the characterization of effector function during pathogenesis. Recent advances in mapping protein-protein interaction networks by means of mass spectrometry-based interactomics have enabled the identification of host targets at large-scale. In this review, we highlight mass spectrometry-driven proteomics strategies and recent advances to elucidate type-III secretion system effector-host protein-protein interactions. Furthermore, we highlight approaches for defining spatial and temporal effector-host interactions, and discuss possible avenues for studying natively delivered effectors in the context of infection. Overall, the knowledge gained when unravelling effector complexation with host factors will provide novel opportunities to control infectious disease outcomes.

T3SEpp: an Integrated Prediction Pipeline for Bacterial Type III Secreted Effectors

Many Gram-negative bacteria infect hosts and cause diseases by translocating a variety of type III secreted effectors (T3SEs) into the host cell cytoplasm. However, despite a dramatic increase in the number of available whole-genome sequences, it remains challenging for accurate prediction of T3SEs. Traditional prediction models have focused on atypical sequence features buried in the N-terminal peptides of T3SEs, but unfortunately, these models have had high false-positive rates. In this research, we integrated promoter information along with characteristic protein features for signal regions, chaperone-binding domains, and effector domains for T3SE prediction. Machine learning algorithms, including deep learning, were adopted to predict the atypical features mainly buried in signal sequences of T3SEs, followed by development of a voting-based ensemble model integrating the individual prediction results. We assembled this into a unified T3SE prediction pipeline, T3SEpp, which integrated the results of individual modules, resulting in high accuracy (i.e., ∼0.94) and >1-fold reduction in the false-positive rate compared to that of state-of-the-art software tools. The T3SEpp pipeline and sequence features observed here will facilitate the accurate identification of new T3SEs, with numerous benefits for future studies on host-pathogen interactions.IMPORTANCE Type III secreted effector (T3SE) prediction remains a big computational challenge. In practical applications, current software tools often suffer problems of high false-positive rates. One of the causal factors could be the relatively unitary type of biological features used for the design and training of the models. In this research, we made a comprehensive survey on the sequence-based features of T3SEs, including signal sequences, chaperone-binding domains, effector domains, and transcription factor binding promoter sites, and assembled a unified prediction pipeline integrating multi-aspect biological features within homology-based and multiple machine learning models. To our knowledge, we have compiled the most comprehensive biological sequence feature analysis for T3SEs in this research. The T3SEpp pipeline integrating the variety of features and assembling different models showed high accuracy, which should facilitate more accurate identification of T3SEs in new and existing bacterial whole-genome sequences.

The pan-genome effector-triggered immunity landscape of a host-pathogen interaction

Effector-triggered immunity (ETI), induced by host immune receptors in response to microbial effectors, protects plants against virulent pathogens. However, a systematic study of ETI prevalence against species-wide pathogen diversity is lacking. We constructed the Pseudomonas syringae Type III Effector Compendium (PsyTEC) to reduce the pan-genome complexity of 5127 unique effector proteins, distributed among 70 families from 494 strains, to 529 representative alleles. We screened PsyTEC on the model plant Arabidopsis thaliana and identified 59 ETI-eliciting alleles (11.2%) from 19 families (27.1%), with orthologs distributed among 96.8% of P. syringae strains. We also identified two previously undescribed host immune receptors, including CAR1, which recognizes the conserved effectors AvrE and HopAA1, and found that 94.7% of strains harbor alleles predicted to be recognized by either CAR1 or ZAR1.

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/.

Molecular Evolution of Pseudomonas syringae Type III Secreted Effector Proteins

Diverse Gram-negative pathogens like Pseudomonas syringae employ type III secreted effector (T3SE) proteins as primary virulence factors that combat host immunity and promote disease. T3SEs can also be recognized by plant hosts and activate an effector triggered immune (ETI) response that shifts the interaction back toward plant immunity. Consequently, T3SEs are pivotal in determining the virulence potential of individual P. syringae strains, and ultimately help to restrict P. syringae pathogens to a subset of potential hosts that are unable to recognize their repertoires of T3SEs. While a number of effector families are known to be present in the P. syringae species complex, one of the most persistent challenges has been documenting the complex variation in T3SE contents across a diverse collection of strains. Using the entire pan-genome of 494 P. syringae strains isolated from more than 100 hosts, we conducted a global analysis of all known and putative T3SEs. We identified a total of 14,613 putative T3SEs, 4,636 of which were unique at the amino acid level, and show that T3SE repertoires of different P. syringae strains vary dramatically, even among strains isolated from the same hosts. We also find substantial diversification within many T3SE families, and in many cases find strong signatures of positive selection. Furthermore, we identify multiple gene gain and loss events for several families, demonstrating an important role of horizontal gene transfer (HGT) in the evolution of P. syringae T3SEs. These analyses provide insight into the evolutionary history of P. syringae T3SEs as they co-evolve with the host immune system, and dramatically expand the database of P. syringae T3SEs alleles.

Caspase-3 cleavage of Salmonella type III secreted effector protein SifA is required for localization of functional domains and bacterial dissemination

SifA is a bi-functional Type III Secretion System (T3SS) effector protein that plays an important role in Salmonella virulence. The N-terminal domain of SifA binds SifA-Kinesin-Interacting-Protein (SKIP), and via an interaction with kinesin, forms tubular membrane extensions called Sif filaments (Sifs) that emanate from the Salmonella Containing Vacuole (SCV). The C-terminal domain of SifA harbors a WxxxE motif that functions to mimic active host cell GTPases. Taken together, SifA functions in inducing endosomal tubulation in order to maintain the integrity of the SCV and promote bacterial dissemination. Since SifA performs multiple, unrelated functions, the objective of this study was to determine how each functional domain of SifA becomes processed. Our work demonstrates that a linker region containing a caspase-3 cleavage motif separates the two functional domains of SifA. To test the hypothesis that processing of SifA by caspase-3 at this particular site is required for function and proper localization of the effector protein domains, we developed two tracking methods to analyze the intracellular localization of SifA. We first adapted a fluorescent tag called phiLOV that allowed for type-III secretion system (T3SS) mediated delivery of SifA and observation of its intracellular colocalization with caspase-3. Additionally, we created a dual-tagging strategy that permitted tracking of each of the SifA functional domains following caspase-3 cleavage to different subcellular locations. The results of this study reveal that caspase-3 cleavage of SifA is required for the proper localization of functional domains and bacterial dissemination. Considering the importance of these events in Salmonella pathogenesis, we conclude that caspase-3 cleavage of effector proteins is a more broadly applicable effector processing mechanism utilized by Salmonella to invade and persist during infection.

Recombination of ecologically and evolutionarily significant loci maintains genetic cohesion in the Pseudomonas syringae species complex

Background

Pseudomonas syringae is a highly diverse bacterial species complex capable of causing a wide range of serious diseases on numerous agronomically important crops. We examine the evolutionary relationships of 391 agricultural and environmental strains using whole-genome sequencing and evolutionary genomic analyses.

Results

We describe the phylogenetic distribution of all 77,728 orthologous gene families in the pan-genome, reconstruct the core genome phylogeny using the 2410 core genes, hierarchically cluster the accessory genome, identify the diversity and distribution of type III secretion systems and their effectors, predict ecologically and evolutionary relevant loci, and establish the molecular evolutionary processes operating on gene families. Phylogenetic and recombination analyses reveals that the species complex is subdivided into primary and secondary phylogroups, with the former primarily comprised of agricultural isolates, including all of the well-studied P. syringae strains. In contrast, the secondary phylogroups include numerous environmental isolates. These phylogroups also have levels of genetic diversity typically found among distinct species. An analysis of rates of recombination within and between phylogroups revealed a higher rate of recombination within primary phylogroups than between primary and secondary phylogroups. We also find that “ecologically significant” virulence-associated loci and “evolutionarily significant” loci under positive selection are over-represented among loci that undergo inter-phylogroup genetic exchange.

Conclusions

While inter-phylogroup recombination occurs relatively rarely, it is an important force maintaining the genetic cohesion of the species complex, particularly among primary phylogroup strains. This level of genetic cohesion, and the shared plant-associated niche, argues for considering the primary phylogroups as a single biological species.

Electronic supplementary material

The online version of this article (10.1186/s13059-018-1606-y) contains supplementary material, which is available to authorized users.

The species-spanning family of LPX-motif harbouring effector proteins

The delivery of effector proteins into infected eukaryotic cells represents a key virulence feature of many microbial pathogens in order to derail essential cellular processes and effectively counter the host defence system. Although bacterial effectors are truly numerous and exhibit a wide range of biochemical activities, commonalities in terms of protein structure and function shared by many bacterial pathogens exist. Recent progress has shed light on a species-spanning family of bacterial effectors containing an LPX repeat motif as a subtype of the leucine-rich repeat superfamily, partially combined with a novel E3 ubiquitin ligase domain. This review highlights the immunomodulatory effects of LPX effector proteins, with particular emphasis on the exploitation of the host ubiquitin system.

Functional Analysis of PsAvr3c Effector Family From Phytophthora Provides Probes to Dissect SKRP Mediated Plant Susceptibility

PsAvr3c is an effector identified from oomycete plant pathogen Phytophthora sojae that causes soybean root and stem rot disease. Earlier studies have demonstrated that PsAvr3c binds to a novel soybean spliceosomal complex protein, GmSKRP, to reprogram the splicing of hundreds of pre-mRNAs and consequently subvert host immunity. PsAvr3c family genes are present in some other Phytophthora species, but their function remains unknown. Here, we characterized the functions of PsAvh27b (PsAvr3c paralog from P. sojae), ProbiAvh89 and PparvAvh214 (orthologs from P. cinnamomi var. robiniae and Phytophthora parvispora, respectively). The study reveals that both PsAvh27b and ProbiAvh89 interact with GmSKRPs in vitro, and stabilize GmSKRP1 in vivo. However, PparvAvh214 cannot interact with GmSKRPs proteins. The qRT-PCR result illustrates that the alternative splicing of pre-mRNAs of several soybean defense-related genes are altered in PsAvh27b and ProbiAvh89 when over-expressed on soybean hairy roots. Moreover, PsAvr3c family members display differences in promoting Phytophthora infection in a SKRP-dependent manner. Overall, this study highlights that the effector-mediated host pre-mRNA alternative splicing occurs in other pathosystems, thus providing new probes to further dissect SKRP-mediated plant susceptibility.

The bacterial type III-secreted protein AvrRps4 is a bipartite effector

Bacterial effector proteins secreted into host plant cells manipulate those cells to the benefit of the pathogen, but effector-triggered immunity (ETI) occurs when effectors are recognized by host resistance proteins. The RPS4/RRS1 pair recognizes the Pseudomonas syringae pv. pisi effector AvrRps4. AvrRps4 is processed in planta into AvrRps4^N (133 amino acids), homologous to the N-termini of other effectors including the native P. syringae pv. tomato strain DC3000 effector HopK1, and AvrRps4^C (88 amino acids). Previous data suggested that AvrRps4^C alone is necessary and sufficient for resistance when overexpressed in heterologous systems. We show that delivering AvrRps4^C from DC3000, but not from a DC3000 hopK1^- strain, triggers resistance in the Arabidopsis accession Col-0. Delivering AvrRps4^C in tandem with AvrRps4^N, or as a chimera with HopK1^N, fully complements AvrRps4-triggered immunity. AvrRps4^N in the absence of AvrRps4^C enhances virulence in Col-0. In addition, AvrRps4^N triggers a hypersensitive response in lettuce that is attenuated by coexpression of AvrRps4^C, further supporting the role of AvrRps4^N as a bona fide effector domain. Based on these results we propose that evolutionarily, fusion of AvrRps4^C to AvrRps4^N may have counteracted recognition of AvrRps4^N, and that the plant RPS4/RRS1 resistance gene pair was selected as a countermeasure. We conclude that AvrRps4 represents an unusual chimeric effector, with recognition in Arabidopsis by RPS4/RRS1 requiring the presence of both processed effector moieties.

An important component of the plant immune system relies on the detection of pathogen-derived effectors by immune receptors called resistance proteins. Bacterial pathogens inject effectors into the host cell via a dedicated secretion system to suppress defenses and to manipulate the physiology of the host cell to the pathogen's advantage. Usually, a single resistance protein recognizes a single effector, but an increasing number of exceptions and elaborations on this one-to-one relationship are known. The plant Arabidopsis uses a pair of resistance proteins, RRS1 and RPS4, to detect the effector AvrRps4. After injection into the cell, AvrRps4 is cleaved into two protein parts, and it had been assumed that only the C-terminal part needs to be present to trigger RPS4/RRS1. We show here that both AvrRps4 parts are required for triggering resistance in Arabidopsis, and that the N-terminal part, which previously had been assumed to only function in effector secretion into the host cell, in fact on its own has some functions of an effector. This conclusion is supported by the observation that the N-terminal part of AvrRps4 is sufficient to trigger resistance in lettuce. The fusion of the two AvrRps4 parts may have arisen to counteract plant defenses.

A bacterial toxin-antitoxin module is the origin of inter-bacterial and inter-kingdom effectors of Bartonella

Host-targeting type IV secretion systems (T4SS) evolved from conjugative T4SS machineries that mediate interbacterial plasmid transfer. However, the origins of effectors secreted by these virulence devices have remained largely elusive. Previous work showed that some effectors exhibit homology to toxins of bacterial toxin-antitoxin modules, but the evolutionary trajectories underlying these ties had not been resolved. We previously reported that FicT toxins of FicTA toxin-antitoxin modules disrupt cellular DNA topology via their enzymatic FIC (filamentation induced by cAMP) domain. Intriguingly, the FIC domain of the FicT toxin VbhT of Bartonella schoenbuchensis is fused to a type IV secretion signal–the BID (Bep intracellular delivery) domain—similar to the Bartonella effector proteins (Beps) that are secreted into eukaryotic host cells via the host-targeting VirB T4SS. In this study, we show that the VbhT toxin is an interbacterial effector protein secreted via the conjugative Vbh T4SS that is closely related to the VirB T4SS and encoded by plasmid pVbh of B. schoenbuchensis. We therefore propose that the Vbh T4SS together with its effector VbhT represent an evolutionary missing link on a path that leads from a regular conjugation system and FicTA toxin-antitoxin modules to the VirB T4SS and the Beps. Intriguingly, phylogenetic analyses revealed that the fusion of FIC and BID domains has probably occurred independently in VbhT and the common ancestor of the Beps, suggesting parallel evolutionary paths. Moreover, several other examples of TA module toxins that are bona fide substrates of conjugative T4SS indicate that their recruitment as interbacterial effectors is prevalent and serves yet unknown biological functions in the context of bacterial conjugation. We propose that the adaptation for interbacterial transfer favors the exaptation of FicT and other TA module toxins as inter-kingdom effectors and may thus constitute an important stepping stone in the evolution of host-targeted effector proteins.

Many bacterial pathogens use secretion systems to translocate effector proteins into host cells where they manipulate cell functions in favor of the pathogen. It is well-known that these secretion systems evolved from ancestors with functions in genuine bacterial contexts, but the origins of their secreted effectors have largely remained elusive. In this article we studied the evolutionary history of a host-targeting effector secretion system of the mammalian pathogen Bartonella that belongs to a group of machineries descended from secretion systems originally mediating DNA transfer between bacterial cells. Intriguingly, we found that such a DNA transfer machinery closely related to the host-targeting secretion system of Bartonella has recruited a bacterial protein involved in modulating DNA topology as an interbacterial effector protein that is translocated together with the DNA into recipient cells. The overall setup of this interbacterial effector is remarkably similar to the host-targeted effectors of Bartonella, and we propose that it represents an evolutionary missing link on the path from a genuine bacterial protein to effectors that manipulates host cell functioning. Further analyses showed that interbacterial effectors in DNA transfer may be a more common phenomenon and represent an important reservoir for the evolution of new host-targeted effectors.

The HopF family of Pseudomonas syringae type III secreted effectors

Pseudomonas syringae is a bacterial phytopathogen that utilizes the type III secretion system to inject effector proteins into plant host cells. Pseudomonas syringae can infect a wide range of plant hosts, including agronomically important crops such as tomatoes and beans. The ability of P. syringae to infect such numerous hosts is caused, in part, by the diversity of effectors employed by this phytopathogen. Over 60 different effector families exist in P. syringae; one such family is HopF, which contains over 100 distinct alleles. Despite this diversity, research has focused on only two members of this family: HopF1 from P. syringae pathovar phaseolicola 1449B and HopF2 from P. syringae pathovar tomato DC3000. In this study, we review the research on HopF family members, including their host targets and molecular mechanisms of immunity suppression, and their enzymatic function. We also provide a phylogenetic analysis of this expanding effector family which provides a basis for a proposed nomenclature to guide future research. The extensive genetic diversity that exists within the HopF family presents a great opportunity to study how functional diversification on an effector family contributes to host specialization.

3'-NADP and 3'-NAADP, Two Metabolites Formed by the Bacterial Type III Effector AvrRxo1

An arsenal of effector proteins is injected by bacterial pathogens into the host cell or its vicinity to increase virulence. The commonly used top-down approaches inferring the toxic mechanism of individual effector proteins from the host's phenotype are often impeded by multiple targets of different effectors as well as by their pleiotropic effects. Here we describe our bottom-up approach, showing that the bacterial type III effector AvrRxo1 of plant pathogens is an authentic phosphotransferase that produces two novel metabolites by phosphorylating nicotinamide/nicotinic acid adenine dinucleotide at the adenosine 3′-hydroxyl group. Both products of AvrRxo1, 3′-NADP and 3′-nicotinic acid adenine dinucleotide phosphate (3′-NAADP), are substantially different from the ubiquitous co-enzyme 2′-NADP and the calcium mobilizer 2′-NAADP. Interestingly, 3′-NADP and 3′-NAADP have previously been used as inhibitors or signaling molecules but were regarded as “artificial” compounds so far. Our findings now necessitate a shift in thinking about the biological importance of 3′-phosphorylated NAD derivatives.

Using Ecology, Physiology, and Genomics to Understand Host Specificity in Xanthomonas

How pathogens coevolve with and adapt to their hosts are critical to understanding how host jumps and/or acquisition of novel traits can lead to new disease emergences. The Xanthomonas genus includes Gram-negative plant-pathogenic bacteria that collectively infect a broad range of crops and wild plant species. However, individual Xanthomonas strains usually cause disease on only a few plant species and are highly adapted to their hosts, making them pertinent models to study host specificity. This review summarizes our current understanding of the molecular basis of host specificity in the Xanthomonas genus, with a particular focus on the ecology, physiology, and pathogenicity of the bacterium. Despite our limited understanding of the basis of host specificity, type III effectors, microbe-associated molecular patterns, lipopolysaccharides, transcriptional regulators, and chemotactic sensors emerge as key determinants for shaping host specificity.

A De-Novo Genome Analysis Pipeline (DeNoGAP) for large-scale comparative prokaryotic genomics studies

Background

Comparative analysis of whole genome sequence data from closely related prokaryotic species or strains is becoming an increasingly important and accessible approach for addressing both fundamental and applied biological questions. While there are number of excellent tools developed for performing this task, most scale poorly when faced with hundreds of genome sequences, and many require extensive manual curation.

Results

We have developed a de-novo genome analysis pipeline (DeNoGAP) for the automated, iterative and high-throughput analysis of data from comparative genomics projects involving hundreds of whole genome sequences. The pipeline is designed to perform reference-assisted and de novo gene prediction, homolog protein family assignment, ortholog prediction, functional annotation, and pan-genome analysis using a range of proven tools and databases. While most existing methods scale quadratically with the number of genomes since they rely on pairwise comparisons among predicted protein sequences, DeNoGAP scales linearly since the homology assignment is based on iteratively refined hidden Markov models. This iterative clustering strategy enables DeNoGAP to handle a very large number of genomes using minimal computational resources. Moreover, the modular structure of the pipeline permits easy updates as new analysis programs become available.

Conclusion

DeNoGAP integrates bioinformatics tools and databases for comparative analysis of a large number of genomes. The pipeline offers tools and algorithms for annotation and analysis of completed and draft genome sequences. The pipeline is developed using Perl, BioPerl and SQLite on Ubuntu Linux version 12.04 LTS. Currently, the software package accompanies script for automated installation of necessary external programs on Ubuntu Linux; however, the pipeline should be also compatible with other Linux and Unix systems after necessary external programs are installed. DeNoGAP is freely available at https://sourceforge.net/projects/denogap/.

Electronic supplementary material

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

AnnoTALE: bioinformatics tools for identification, annotation, and nomenclature of TALEs from Xanthomonas genomic sequences

Transcription activator-like effectors (TALEs) are virulence factors, produced by the bacterial plant-pathogen Xanthomonas, that function as gene activators inside plant cells. Although the contribution of individual TALEs to infectivity has been shown, the specific roles of most TALEs, and the overall TALE diversity in Xanthomonas spp. is not known. TALEs possess a highly repetitive DNA-binding domain, which is notoriously difficult to sequence. Here, we describe an improved method for characterizing TALE genes by the use of PacBio sequencing. We present ‘AnnoTALE’, a suite of applications for the analysis and annotation of TALE genes from Xanthomonas genomes, and for grouping similar TALEs into classes. Based on these classes, we propose a unified nomenclature for Xanthomonas TALEs that reveals similarities pointing to related functionalities. This new classification enables us to compare related TALEs and to identify base substitutions responsible for the evolution of TALE specificities.

Crystal Structure of Xanthomonas AvrRxo1-ORF1, a Type III Effector with a Polynucleotide Kinase Domain, and Its Interactor AvrRxo1-ORF2

Xanthomonas oryzae pv. oryzicola (Xoc) causes bacterial leaf streak (BLS) disease on rice plants. Xoc delivers a type III effector AvrRxo1-ORF1 into rice plant cells that can be recognized by disease resistance (R) protein Rxo1, and triggers resistance to BLS disease. However, the mechanism and virulence role of AvrRxo1 is not known. In the genome of Xoc, AvrRxo1-ORF1 is adjacent to another gene AvrRxo1-ORF2, which was predicted to encode a molecular chaperone of AvrRxo1-ORF1. We report the co-purification and crystallization of the AvrRxo1-ORF1:AvrRxo1-ORF2 tetramer complex at 1.64 Å resolution. AvrRxo1-ORF1 has a T4 polynucleotide kinase domain, and expression of AvrRxo1-ORF1 suppresses bacterial growth in a manner dependent on the kinase motif. Although AvrRxo1-ORF2 binds AvrRxo1-ORF1, it is structurally different from typical effector-binding chaperones, in that it has a distinct fold containing a novel kinase-binding domain. AvrRxo1-ORF2 functions to suppress the bacteriostatic activity of AvrRxo1-ORF1 in bacterial cells.

BEAN 2.0: an integrated web resource for the identification and functional analysis of type III secreted effectors

Gram-negative pathogenic bacteria inject type III secreted effectors (T3SEs) into host cells to sabotage their immune signaling networks. Because T3SEs constitute a meeting-point of pathogen virulence and host defense, they are of keen interest to host-pathogen interaction research community. To accelerate the identification and functional understanding of T3SEs, we present BEAN 2.0 as an integrated web resource to predict, analyse and store T3SEs. BEAN 2.0 includes three major components. First, it provides an accurate T3SE predictor based on a hybrid approach. Using independent testing data, we show that BEAN 2.0 achieves a sensitivity of 86.05% and a specificity of 100%. Second, it integrates a set of online sequence analysis tools. Users can further perform functional analysis of putative T3SEs in a seamless way, such as subcellular location prediction, functional domain scan and disorder region annotation. Third, it compiles a database covering 1215 experimentally verified T3SEs and constructs two T3SE-related networks that can be used to explore the relationships among T3SEs. Taken together, by presenting a one-stop T3SE bioinformatics resource, we hope BEAN 2.0 can promote comprehensive understanding of the function and evolution of T3SEs.

Host-microbe and microbe-microbe interactions in the evolution of obligate plant parasitism

Research on obligate biotrophic plant parasites, which reproduce only on living hosts, has revealed a broad diversity of filamentous microbes that have independently acquired complex morphological structures, such as haustoria. Genome studies have also demonstrated a concerted loss of genes for metabolism and lytic enzymes, and gain of diversity of genes coding for effectors involved in host defense suppression. So far, these traits converge in all known obligate biotrophic parasites, but unexpected genome plasticity remains. This plasticity is manifested as transposable element (TE)-driven increases in genome size, observed to be associated with the diversification of virulence genes under selection pressure. Genome expansion could result from the governing of the pathogen response to ecological selection pressures, such as host or nutrient availability, or to microbial interactions, such as competition, hyperparasitism and beneficial cooperations. Expansion is balanced by alternating sexual and asexual cycles, as well as selfing and outcrossing, which operate to control transposon activity in populations. In turn, the prevalence of these balancing mechanisms seems to be correlated with external biotic factors, suggesting a complex, interconnected evolutionary network in host-pathogen-microbe interactions. Therefore, the next phase of obligate biotrophic pathogen research will need to uncover how this network, including multitrophic interactions, shapes the evolution and diversity of pathogens.

A Ralstonia solanacearum type III effector directs the production of the plant signal metabolite trehalose-6-phosphate

The plant pathogen Ralstonia solanacearum possesses two genes encoding a trehalose-6-phosphate synthase (TPS), an enzyme of the trehalose biosynthetic pathway. One of these genes, named ripTPS, was found to encode a protein with an additional N-terminal domain which directs its translocation into host plant cells through the type 3 secretion system. RipTPS is a conserved effector in the R. solanacearum species complex, and homologues were also detected in other bacterial plant pathogens. Functional analysis of RipTPS demonstrated that this type 3 effector synthesizes trehalose-6-phosphate and identified residues essential for this enzymatic activity. Although trehalose-6-phosphate is a key signal molecule in plants that regulates sugar status and carbon assimilation, the disruption of ripTPS did not alter the virulence of R. solanacearum on plants. However, heterologous expression assays showed that this effector specifically elicits a hypersensitive-like response on tobacco that is independent of its enzymatic activity and is triggered by the C-terminal half of the protein. Recognition of this effector by the plant immune system is suggestive of a role during the infectious process.

Ralstonia solanacearum, the causal agent of bacterial wilt disease, infects more than two hundred plant species, including economically important crops. The type III secretion system plays a major role in the pathogenicity of this bacterium, and approximately 70 effector proteins have been shown to be translocated into host plant cells. This study provides the first description of a type III effector endowed with a trehalose-6-phosphate synthase enzymatic activity and illustrates a new mechanism by which the bacteria may manipulate the plant metabolism upon infection. In recent years, trehalose-6-phosphate has emerged as an essential signal molecule in plants, connecting plant metabolism and development. The finding that a bacterial pathogen could induce the production of trehalose-6-phosphate in plant cells further highlights the importance of this metabolite in multiple aspects of the molecular physiology of plants.

New insights into the role of Bartonella effector proteins in pathogenesis

The facultative intracellular bacteria Bartonella spp. share a common infection strategy to invade and colonize mammals in a host-specific manner. Following transmission by blood-sucking arthropods, Bartonella are inoculated in the derma and then spread, via two sequential enigmatic niches, to the blood stream where they cause a long-lasting intra-erythrocytic bacteraemia. The VirB/VirD4 type IV secretion system (VirB/D4 T4SS) is essential for the pathogenicity of most Bartonella species by injecting an arsenal of effector proteins into host cells. These bacterial effector proteins share a modular architecture, comprising domains and/or motifs that confer an array of functions. Here, we review recent advances in understanding the function and evolutionary origin of this fascinating repertoire of host-targeted bacterial effectors.

Bacterial secreted effectors and caspase-3 interactions

Apoptosis is a critical process that intrinsically links organism survival to its ability to induce controlled death. Thus, functional apoptosis allows organisms to remove perceived threats to their survival by targeting those cells that it determines pose a direct risk. Central to this process are apoptotic caspases, enzymes that form a signalling cascade, converting danger signals via initiator caspases into activation of the executioner caspase, caspase-3. This enzyme begins disassembly of the cell by activating DNA degrading enzymes and degrading the cellular architecture. Interaction of pathogenic bacteria with caspases, and in particular, caspase-3, can therefore impact both host cell and bacterial survival. With roles outside cell death such as cell differentiation, control of signalling pathways and immunomodulation also being described for caspase-3, bacterial interactions with caspase-3 may be of far more significance in infection than previously recognized. In this review, we highlight the ways in which bacterial pathogens have evolved to subvert caspase-3 both through effector proteins that directly interact with the enzyme or by modulating pathways that influence its activation and activity.

Mobile effector proteins on phage genomes

Bacteriophage genomes found in a range of bacterial pathogens encode a diverse array of virulence factors ranging from superantigens or pore forming lysins to numerous exotoxins. Recent studies have uncovered an entirely new class of bacterial virulence factors, called effector proteins or effector toxins, which are encoded within phage genomes that reside among several pathovars of Escherichia coli and Salmonella enterica. These effector proteins have multiple domains resulting in proteins that can be multifunctional. The effector proteins encoded within phage genomes are translocated directly from the bacterial cytosol into their eukaryotic target cells by specialized bacterial type three secretion systems (T3SSs). In this review, we will give an overview of the different types of effector proteins encoded within phage genomes and examine their roles in bacterial pathogenesis.

Genomic plasticity enables phenotypic variation of Pseudomonas syringae pv. tomato DC3000

Whole genome sequencing revealed the presence of a genomic anomaly in the region of 4.7 to 4.9 Mb of the Pseudomonas syringae pv. tomato (Pst) DC3000 genome. The average read depth coverage of Pst DC3000 whole genome sequencing results suggested that a 165 kb segment of the chromosome had doubled in copy number. Further analysis confirmed the 165 kb duplication and that the two copies were arranged as a direct tandem repeat. Examination of the corresponding locus in Pst NCPPB1106, the parent strain of Pst DC3000, suggested that the 165 kb duplication most likely formed after the two strains diverged via transposition of an ISPsy5 insertion sequence (IS) followed by unequal crossing over between ISPsy5 elements at each end of the duplicated region. Deletion of one copy of the 165 kb region demonstrated that the duplication facilitated enhanced growth in some culture conditions, but did not affect pathogenic growth in host tomato plants. These types of chromosomal structures are predicted to be unstable and we have observed resolution of the 165 kb duplication to single copy and its subsequent re-duplication. These data demonstrate the role of IS elements in recombination events that facilitate genomic reorganization in P. syringae.

Variable suites of non-effector genes are co-regulated in the type III secretion virulence regulon across the Pseudomonas syringae phylogeny

Pseudomonas syringae is a phylogenetically diverse species of Gram-negative bacterial plant pathogens responsible for crop diseases around the world. The HrpL sigma factor drives expression of the major P. syringae virulence regulon. HrpL controls expression of the genes encoding the structural and functional components of the type III secretion system (T3SS) and the type three secreted effector proteins (T3E) that are collectively essential for virulence. HrpL also regulates expression of an under-explored suite of non-type III effector genes (non-T3E), including toxin production systems and operons not previously associated with virulence. We implemented and refined genome-wide transcriptional analysis methods using cDNA-derived high-throughput sequencing (RNA-seq) data to characterize the HrpL regulon from six isolates of P. syringae spanning the diversity of the species. Our transcriptomes, mapped onto both complete and draft genomes, significantly extend earlier studies. We confirmed HrpL-regulation for a majority of previously defined T3E genes in these six strains. We identified two new T3E families from P. syringae pv. oryzae 1_6, a strain within the relatively underexplored phylogenetic Multi-Locus Sequence Typing (MLST) group IV. The HrpL regulons varied among strains in gene number and content across both their T3E and non-T3E gene suites. Strains within MLST group II consistently express the lowest number of HrpL-regulated genes. We identified events leading to recruitment into, and loss from, the HrpL regulon. These included gene gain and loss, and loss of HrpL regulation caused by group-specific cis element mutations in otherwise conserved genes. Novel non-T3E HrpL-regulated genes include an operon that we show is required for full virulence of P. syringae pv. phaseolicola 1448A on French bean. We highlight the power of integrating genomic, transcriptomic, and phylogenetic information to drive concise functional experimentation and to derive better insight into the evolution of virulence across an evolutionarily diverse pathogen species.

Host immune responses accelerate pathogen evolution

Pathogens face a hostile and often novel environment when infecting a new host, and adaptation is likely to be an important determinant of the success in colonization and establishment. We hypothesized that resistant hosts will impose stronger selection on pathogens than susceptible hosts, which should accelerate pathogen evolution through selection biased toward effector genes. To test this hypothesis, we conducted an experimental evolution study on Xanthomonas citri subsp. citri (Xcc) in a susceptible plant species and a resistant plant species. We performed 55 rounds of repeated reinoculation of Xcc through susceptible host grapefruit (isolates G1, G2, G3) and resistant host kumquat (isolates K1, K2, K3). Consequently, only K1 and K3 isolates lost their ability to elicit a hypersensitive response (HR) in kumquat. Illumina sequencing of the parental and descendant strains P, G1, G2, G3, K1, K2 and K3 revealed that fixed mutations were biased toward type three secretion system effectors in isolates K1 and K3. Parallel evolution was observed in the K1 and K3 strains, suggesting that the mutations result from selection rather than by random drift. Our results support our hypothesis and suggest that repeated infection of resistant hosts by pathogens should be prevented to avoid selecting for adaptive pathogens.

Evolution of RXLR-class effectors in the oomycete plant pathogen Phytophthora ramorum

Phytophthora plant pathogens contain many hundreds of effectors potentially involved in infection of host plants. Comparative genomic analyses have shown that these effectors evolve rapidly and have been subject to recent expansions. We examined the recent sequence evolution of RXLR-class effector gene families in the sudden oak death pathogen, P. ramorum. We found that P. ramorum RXLR effectors have taken multiple evolutionary paths, including loss or gain of repeated domains, recombination or gene conversion among paralogs, and selection on point mutations. Sequencing of homologs from two subfamilies in P. ramorum’s closest known relatives revealed repeated gene duplication and divergence since speciation with P. lateralis. One family showed strong signatures of recombination while the other family has evolved primarily by point mutation. Comparison of a small number of the hundreds of RXLR-class effectors across three clonal lineages of P. ramorum shows striking divergence in alleles among lineages, suggesting the potential for functional differences between lineages. Our results suggest future avenues for examination of rapidly evolving effectors in P. ramorum, including investigation of the functional and coevolutionary significance of the patterns of sequence evolution that we observed.

Genome sequence of Xanthomonas fuscans subsp. fuscans strain 4834-R reveals that flagellar motility is not a general feature of xanthomonads

BACKGROUND

Xanthomonads are plant-associated bacteria responsible for diseases on economically important crops. Xanthomonas fuscans subsp. fuscans (Xff) is one of the causal agents of common bacterial blight of bean. In this study, the complete genome sequence of strain Xff 4834-R was determined and compared to other Xanthomonas genome sequences.

RESULTS

Comparative genomics analyses revealed core characteristics shared between Xff 4834-R and other xanthomonads including chemotaxis elements, two-component systems, TonB-dependent transporters, secretion systems (from T1SS to T6SS) and multiple effectors. For instance a repertoire of 29 Type 3 Effectors (T3Es) with two Transcription Activator-Like Effectors was predicted. Mobile elements were associated with major modifications in the genome structure and gene content in comparison to other Xanthomonas genomes. Notably, a deletion of 33 kbp affects flagellum biosynthesis in Xff 4834-R. The presence of a complete flagellar cluster was assessed in a collection of more than 300 strains representing different species and pathovars of Xanthomonas. Five percent of the tested strains presented a deletion in the flagellar cluster and were non-motile. Moreover, half of the Xff strains isolated from the same epidemic than 4834-R was non-motile and this ratio was conserved in the strains colonizing the next bean seed generations.

CONCLUSIONS

This work describes the first genome of a Xanthomonas strain pathogenic on bean and reports the existence of non-motile xanthomonads belonging to different species and pathovars. Isolation of such Xff variants from a natural epidemic may suggest that flagellar motility is not a key function for in planta fitness.

The role of protein effectors in plant-aphid interactions

Aphid salivary proteins, which are injected into the phloem sieve elements during feeding, play a central role in plant-aphid interactions. Among the dozens of known salivary proteins, many have no homology to proteins from other organisms. These aphid-specific proteins likely have evolved as effectors that inhibit plant defenses, prevent phloem sieve-element occlusion, and otherwise promote the unique phloem feeding style. However, aphid salivary proteins also are recognized by plants to mount defense responses and are likely a major factor in limiting the host range of particular aphid species and biotypes. Newly developed research tools provide excellent opportunities for analyzing the mostly unknown functions of aphid salivary proteins and elucidating their contribution to the complex interactions between aphids and their host plants.

Gene duplication and fragment recombination drive functional diversification of a superfamily of cytoplasmic effectors in Phytophthora sojae

Phytophthora and other oomycetes secrete a large number of putative host cytoplasmic effectors with conserved FLAK motifs following signal peptides, termed crinkling and necrosis inducing proteins (CRN), or Crinkler. Here, we first investigated the evolutionary patterns and mechanisms of CRN effectors in Phytophthora sojae and compared them to two other Phytophthora species. The genes encoding CRN effectors could be divided into 45 orthologous gene groups (OGG), and most OGGs unequally distributed in the three species, in which each underwent large number of gene gains or losses, indicating that the CRN genes expanded after species evolution in Phytophthora and evolved through pathoadaptation. The 134 expanded genes in P. sojae encoded family proteins including 82 functional genes and expressed at higher levels while the other 68 genes encoding orphan proteins were less expressed and contained 50 pseudogenes. Furthermore, we demonstrated that most expanded genes underwent gene duplication or/and fragment recombination. Three different mechanisms that drove gene duplication or recombination were identified. Finally, the expanded CRN effectors exhibited varying pathogenic functions, including induction of programmed cell death (PCD) and suppression of PCD through PAMP-triggered immunity or/and effector-triggered immunity. Overall, these results suggest that gene duplication and fragment recombination may be two mechanisms that drive the expansion and neofunctionalization of the CRN family in P. sojae, which aids in understanding the roles of CRN effectors within each oomycete pathogen.

Effective identification of bacterial type III secretion signals using joint element features

Type III secretion system (T3SS) plays important roles in bacteria and host cell interactions by specifically translocating type III effectors into the cytoplasm of the host cells. The N-terminal amino acid sequences of the bacterial type III effectors determine their specific secretion via type III secretion conduits. It is still unclear as to how the N-terminal sequences guide this specificity. In this work, the amino acid composition, secondary structure, and solvent accessibility in the N-termini of type III and non-type III secreted proteins were compared and contrasted. A high-efficacy mathematical model based on these joint features was developed to distinguish the type III proteins from the non-type III ones. The results indicate that secondary structure and solvent accessibility may make important contribution to the specific recognition of type III secretion signals. Analysis also showed that the joint feature of the N-terminal 6^th–10^th amino acids are especially important for guiding specific type III secretion. Furthermore, a genome-wide screening was performed to predict Salmonella type III secreted proteins, and 8 new candidates were experimentally validated. Interestingly, type III secretion signals were also predicted in gram-positive bacteria and yeasts. Experimental validation showed that two candidates from yeast can indeed be secreted through Salmonella type III secretion conduit. This research provides the first line of direct evidence that secondary structure and solvent accessibility contain important features for guiding specific type III secretion. The new software based on these joint features ensures a high accuracy (general cross-validation sensitivity of ∼96% at a specificity of ∼98%) in silico identification of new type III secreted proteins, which may facilitate our understanding about the specificity of type III secretion and the evolution of type III secreted proteins.

The apoptogenic toxin AIP56 is a metalloprotease A-B toxin that cleaves NF-κb P65

AIP56 (apoptosis-inducing protein of 56 kDa) is a major virulence factor of Photobacterium damselae piscicida (Phdp), a Gram-negative pathogen that causes septicemic infections, which are among the most threatening diseases in mariculture. The toxin triggers apoptosis of host macrophages and neutrophils through a process that, in vivo, culminates with secondary necrosis of the apoptotic cells contributing to the necrotic lesions observed in the diseased animals. Here, we show that AIP56 is a NF-κB p65-cleaving zinc-metalloprotease whose catalytic activity is required for the apoptogenic effect. Most of the bacterial effectors known to target NF-κB are type III secreted effectors. In contrast, we demonstrate that AIP56 is an A-B toxin capable of acting at distance, without requiring contact of the bacteria with the target cell. We also show that the N-terminal domain cleaves NF-κB at the Cys³⁹-Glu⁴⁰ peptide bond and that the C-terminal domain is involved in binding and internalization into the cytosol.

The apoptosis inducing protein of 56 kDa (AIP56) is a key virulence factor secreted by Photobacterium damselae piscicida (Phdp), a Gram-negative bacterium that causes septicaemic infections in economically important marine fish species. It is known that AIP56 induces massive destruction of the phagocytic cells of the infected host, allowing the extracellular multiplication of the bacteria and contributing to the genesis of the pathology. Here we show that AIP56 acts by cleaving NF-κB p65. The NF-κB family of transcription factors is evolutionarily conserved and plays a central role in the host responses to microbial pathogen invasion, regulating the expression of inflammatory and anti-apoptotic genes. Pathogenic bacteria have evolved complex strategies to interfere with NF-κB signalling, usually by injecting protein effectors directly into the cell's cytosol through bacterial secretion machineries that require contact with host cells. In contrast, AIP56 acts at distance and has an intrinsic ability to reach the cytosol due to the presence of a C-terminal domain that functions as “delivery module.”

The mangotoxin biosynthetic operon (mbo) is specifically distributed within Pseudomonas syringae genomospecies 1 and was acquired only once during evolution

Mangotoxin production was first described in Pseudomonas syringae pv. syringae strains. A phenotypic characterization of 94 P. syringae strains was carried out to determine the genetic evolution of the mangotoxin biosynthetic operon (mbo). We designed a PCR primer pair specific for the mbo operon to examine its distribution within the P. syringae complex. These primers amplified a 692-bp DNA fragment from 52 mangotoxin-producing strains and from 7 non-mangotoxin-producing strains that harbor the mbo operon, whereas 35 non-mangotoxin-producing strains did not yield any amplification. This, together with the analysis of draft genomes, allowed the identification of the mbo operon in five pathovars (pathovars aptata, avellanae, japonica, pisi, and syringae), all of which belong to genomospecies 1, suggesting a limited distribution of the mbo genes in the P. syringae complex. Phylogenetic analyses using partial sequences from housekeeping genes differentiated three groups within genomospecies 1. All of the strains containing the mbo operon clustered in groups I and II, whereas those lacking the operon clustered in group III; however, the relative branching order of these three groups is dependent on the genes used to construct the phylogeny. The mbo operon maintains synteny and is inserted in the same genomic location, with high sequence conservation around the insertion point, for all the strains in groups I and II. These data support the idea that the mbo operon was acquired horizontally and only once by the ancestor of groups I and II from genomospecies 1 within the P. syringae complex.

Promoter propagation in prokaryotes

Transcriptional activation or ‘rewiring’ of silent genes is an important, yet poorly understood, phenomenon in prokaryotic genomes. Anecdotal evidence coming from experimental evolution studies in bacterial systems has shown the promptness of adaptation upon appropriate selective pressure. In many cases, a partial or complete promoter is mobilized to silent genes from elsewhere in the genome. We term hereafter such recruited regulatory sequences as Putative Mobile Promoters (PMPs) and we hypothesize they have a large impact on rapid adaptation of novel or cryptic functions. Querying all publicly available prokaryotic genomes (1362) uncovered >4000 families of highly conserved PMPs (50 to 100 long with ≥80% nt identity) in 1043 genomes from 424 different genera. The genomes with the largest number of PMP families are Anabaena variabilis (28 families), Geobacter uraniireducens (27 families) and Cyanothece PCC7424 (25 families). Family size varied from 2 to 93 homologous promoters (in Desulfurivibrio alkaliphilus). Some PMPs are present in particular species, but some are conserved across distant genera. The identified PMPs represent a conservative dataset of very recent or conserved events of mobilization of non-coding DNA and thus they constitute evidence of an extensive reservoir of recyclable regulatory sequences for rapid transcriptional rewiring.

Computational analysis of interactomes: current and future perspectives for bioinformatics approaches to model the host-pathogen interaction space

Bacterial and viral pathogens affect their eukaryotic host partly by interacting with proteins of the host cell. Hence, to investigate infection from a systems' perspective we need to construct complete and accurate host-pathogen protein-protein interaction networks. Because of the paucity of available data and the cost associated with experimental approaches, any construction and analysis of such a network in the near future has to rely on computational predictions. Specifically, this challenge consists of a number of sub-problems: First, prediction of possible pathogen interactors (e.g. effector proteins) is necessary for bacteria and protozoa. Second, the prospective host binding partners have to be determined and finally, the impact on the host cell analyzed. This review gives an overview of current bioinformatics approaches to obtain and understand host-pathogen interactions. As an application example of the methods covered, we predict host-pathogen interactions of Salmonella and discuss the value of these predictions as a prospective for further research.

Host specificity determinants as a genetic continuum

Host specificity is an important concept that underlies the interaction of all clinically and agriculturally relevant microbes with their hosts. Changes in the host specificity of animal pathogens, in particular, are often of greatest concern due to their immediate and unexpected impact on human health. Host switching or host jumps can often be traced to modification of key microbial pathogenicity factors that facilitate the formation of particular host associations. An increase in the number of genome-level studies has begun revealing that almost any type of change, from the simplest to the most complex, can potentially impact host specificity. This review highlights examples of host specificity determinants of viruses, bacteria and fungi, and presents them from within a genetic continuum that spans from the single residue through to entire genomic islands.

E622, a miniature, virulence-associated mobile element

Miniature inverted terminal repeat elements (MITEs) are nonautonomous mobile elements that have a significant impact on bacterial evolution. Here we characterize E622, a 611-bp virulence-associated MITE from Pseudomonas syringae, which contains no coding region but has almost perfect 168-bp inverted repeats. Using an antibiotic coupling assay, we show that E622 is transposable and can mobilize an antibiotic resistance gene contained between its borders. Its predicted parent element, designated TnE622, has a typical transposon structure with a three-gene operon, consisting of resolvase, integrase, and exeA-like genes, which is bounded by the same terminal inverted repeats as E622. A broader genome level survey of the E622/TnE622 inverted repeats identified homologs in Pseudomonas, Salmonella, Shewanella, Erwinia, Pantoea, and the cyanobacteria Nostoc and Cyanothece, many of which appear to encompass known virulence genes, including genes encoding toxins, enzymes, and type III secreted effectors. Its association with niche-specific genetic determinants, along with its persistence and evolutionary diversification, indicates that this mobile element family has played a prominent role in the evolution of many agriculturally and clinically relevant pathogenic bacteria.

Involvement of Type IV Pili in Pathogenicity of Plant Pathogenic Bacteria

Type IV pili (T4P) are hair-like appendages found on the surface of a wide range of bacteria belonging to the β-, γ-, and δ-Proteobacteria, Cyanobacteria and Firmicutes. They constitute an efficient device for a particular type of bacterial surface motility, named twitching, and are involved in several other bacterial activities and functions, including surface adherence, colonization, biofilm formation, genetic material uptake and virulence. Tens of genes are involved in T4P synthesis and regulation, with the majority of them being generally named pil/fim genes. Despite the multiple functionality of T4P and their well-established role in pathogenicity of animal pathogenic bacteria, relatively little attention has been given to the role of T4P in plant pathogenic bacteria. Only in recent years studies have begun to examine with more attention the relevance of these surface appendages for virulence of plant bacterial pathogens. The aim of this review is to summarize the current knowledge about T4P genetic machinery and its role in the interactions between phytopathogenic bacteria and their plant hosts.

Sequence and role in virulence of the three plasmid complement of the model tumor-inducing bacterium Pseudomonas savastanoi pv. savastanoi NCPPB 3335

Pseudomonas savastanoi pv. savastanoi NCPPB 3335 is a model for the study of the molecular basis of disease production and tumor formation in woody hosts, and its draft genome sequence has been recently obtained. Here we closed the sequence of the plasmid complement of this strain, composed of three circular molecules of 78,357 nt (pPsv48A), 45,220 nt (pPsv48B), and 42,103 nt (pPsv48C), all belonging to the pPT23A-like family of plasmids widely distributed in the P. syringae complex. A total of 152 coding sequences were predicted in the plasmid complement, of which 38 are hypothetical proteins and seven correspond to putative virulence genes. Plasmid pPsv48A contains an incomplete Type IVB secretion system, the type III secretion system (T3SS) effector gene hopAF1, gene ptz, involved in cytokinin biosynthesis, and three copies of a gene highly conserved in plant-associated proteobacteria, which is preceded by a hrp box motif. A complete Type IVA secretion system, a well conserved origin of transfer (oriT), and a homolog of the T3SS effector gene hopAO1 are present in pPsv48B, while pPsv48C contains a gene with significant homology to isopentenyl-diphosphate delta-isomerase, type 1. Several potential mobile elements were found on the three plasmids, including three types of MITE, a derivative of IS801, and a new transposon effector, ISPsy30. Although the replication regions of these three plasmids are phylogenetically closely related, their structure is diverse, suggesting that the plasmid architecture results from an active exchange of sequences. Artificial inoculations of olive plants with mutants cured of plasmids pPsv48A and pPsv48B showed that pPsv48A is necessary for full virulence and for the development of mature xylem vessels within the knots; we were unable to obtain mutants cured of pPsv48C, which contains five putative toxin-antitoxin genes.